Toward automatic reconstruction of a highly resolved tree of life

Novel misos shape distinct microbial ecologies: opportunities for flavourful sustainable food innovation

Fermentation is resurgent around the world as people seek healthier, more sustainable, and tasty food options. This study explores the microbial ecology of miso, a traditional Japanese fermented paste, made with novel regional substrates to develop new plant-based foods. Eight novel miso varieties were developed using different protein-rich substrates: yellow peas, Gotland lentils, and fava beans (each with two treatments: standard and nixtamalisation), as well as rye bread and soybeans. The misos were produced at Noma, a restaurant in Copenhagen, Denmark. Samples were analysed with biological and technical triplicates at the beginning and end of fermentation. We also incorporated in this study six samples of novel misos produced following the same recipe at Inua, a former affiliate restaurant of Noma in Tokyo, Japan. To analyse microbial community structure and diversity, metabarcoding (16S and ITS) and shotgun metagenomic analyses were performed. The misos contain a greater range of microbes than is currently described for miso in the literature. The composition of the novel yellow pea misos was notably similar to the traditional soybean ones, suggesting they are a good alternative, which supports our culinary collaborators' sensory conclusions. For bacteria, we found that overall substrate had the strongest effect, followed by time, treatment (nixtamalisation), and geography. For fungi, there was a slightly stronger effect of geography and a mild effect of substrate, and no significant effects for treatment or time. Based on an analysis of metagenome-assembled genomes (MAGs), strains of Staphylococccus epidermidis differentiated according to substrate. Carotenoid biosynthesis genes in these MAGs appeared in strains from Japan but not from Denmark, suggesting a possible gene-level geographical effect. The benign and possibly functional presence of S. epidermidis in these misos, a species typically associated with the human skin microbiome, suggests possible adaptation to the miso niche, and the flow of microbes between bodies and foods in certain fermentation as more common than is currently recognised. This study improves our understanding of miso ecology, highlights the potential for developing novel misos using diverse local ingredients, and suggests how fermentation innovation can contribute to studies of microbial ecology and evolution.

The utilization of an unconventional approach to introduce basic bacteriology in a medical school bridge program

Bacteria form an intense portion of reading and learning for students enrolled in microbiology education. As a part of the foundational course outline of bacteriology, bacterial classification is a significant topic of discussion. The purpose of our study was to analyze whether bacterial classification can be taught with a phylogenetic tree approach that might be more engaging and beneficial to student learners of microbiology. This methodology is unique compared to the conventional approach applied in introductory lectures of bacteriology that relies on morphology and Gram-staining to classify bacteria. The participants of this study were students enrolled in a two-semester medical school bridge program that offers a Master's degree in Pre-clinical Sciences. We presented bacterial origin and classification in the light of evolution and used a phylogenetic tree to signify clinically relevant groups of bacteria. Students were also taught the traditional bacterial classification using Gram stains and morphology. Both methods of classification were delivered in a didactic classroom session considering equal time spent and utilizing the same format. An online survey was distributed to the students after the session to collect their feedback. The results from the survey showed that 74% of participants would prefer learning bacterial classification using a combined approach that includes both Gram-staining and morphology as well as the phylogenetic tree. When asked if the study of bacterial classification through an evolutionary tree diagram is a clear and concise way of understanding bacteria, 79% of the students either agreed or strongly agreed with this statement. Interestingly, the alternative phylogenetic tree approach was considered more engaging and regarded as a means to expand the clinical knowledge of bacteria by 78% and 71% of the students, respectively. Overall, our study strongly supports the use of tree-based classification as an additional method to improve the learning of medically important groups of bacteria at varying levels of education.

Streamlining Protein Fractional Synthesis Rates Using SP3 Beads and Stable Isotope Mass Spectrometry: A Case Study on the Plant Ribosome

Ribosomes are an archetypal ribonucleoprotein assembly. Due to ribosomal evolution and function, r-proteins share specific physicochemical similarities, making the riboproteome particularly suited for tailored proteome profiling methods. Moreover, the structural proteome of ribonucleoprotein assemblies reflects context-dependent functional features. Thus, characterizing the state of riboproteomes provides insights to uncover the context-dependent functionality of r-protein rearrangements, as they relate to what has been termed the ribosomal code, a concept that parallels that of the histone code, in which chromatin rearrangements influence gene expression. Compared to high-resolution ribosomal structures, omics methods lag when it comes to offering customized solutions to close the knowledge gap between structure and function that currently exists in riboproteomes. Purifying the riboproteome and subsequent shot-gun proteomics typically involves protein denaturation and digestion with proteases. The results are relative abundances of r-proteins at the ribosome population level. We have previously shown that, to gain insight into the stoichiometry of individual proteins, it is necessary to measure by proteomics bound r-proteins and normalize their intensities by the sum of r-protein abundances per ribosomal complex, i.e., 40S or 60S subunits. These calculations ensure that individual r-protein stoichiometries represent the fraction of each family/paralog relative to the complex, effectively revealing which r-proteins become substoichiometric in specific physiological scenarios. Here, we present an optimized method to profile the riboproteome of any organism as well as the synthesis rates of r-proteins determined by stable isotope-assisted mass spectrometry. Our method purifies the r-proteins in a reversibly denatured state, which offers the possibility for combined top-down and bottom-up proteomics. Our method offers a milder native denaturation of the r-proteome via a chaotropic GuHCl solution as compared with previous studies that use irreversible denaturation under highly acidic conditions to dissociate rRNA and r-proteins. As such, our method is better suited to conserve post-translational modifications (PTMs). Subsequently, our method carefully considers the amino acid composition of r-proteins to select an appropriate protease for digestion. We avoid non-specific protease cleavage by increasing the pH of our standardized r-proteome dilutions that enter the digestion pipeline and by using a digestion buffer that ensures an optimal pH for a reliable protease digestion process. Finally, we provide the R package ProtSynthesis to study the fractional synthesis rates of r-proteins. The package uses physiological parameters as input to determine peptide or protein fractional synthesis rates. Once the physiological parameters are measured, our equations allow a fair comparison between treatments that alter the biological equilibrium state of the system under study. Our equations correct peptide enrichment using enrichments in soluble amino acids, growth rates, and total protein accumulation. As a means of validation, our pipeline fails to find "false" enrichments in non-labeled samples while also filtering out proteins with multiple unique peptides that have different enrichment values, which are rare in our datasets. These two aspects reflect the accuracy of our tool. Our method offers the possibility of elucidating individual r-protein family/paralog abundances, PTM status, fractional synthesis rates, and dynamic assembly into ribosomal complexes if top-down and bottom-up proteomic approaches are used concomitantly, taking one step further into mapping the native and dynamic status of the r-proteome onto high-resolution ribosome structures. In addition, our method can be used to study the proteomes of all macromolecular assemblies that can be purified, although purification is the limiting step, and the efficacy and accuracy of the proteases may be limited depending on the digestion requirements. Key features • Efficient purification of the ribosomal proteome: streamlined procedure for the specific purification of the ribosomal proteome or complex Ome. • Accurate calculation of fractional synthesis rates: robust method for calculating fractional protein synthesis rates in macromolecular complexes under different physiological steady states. • Holistic ribosome methodology focused on plants: comprehensive approach that provides insights into the ribosomes and translational control of plants, demonstrated using cold acclimation [1]. • Tailored strategies for stable isotope labeling in plants: methodology focusing on materials and labeling considerations specific to free and proteinogenic amino acid analysis [2].

Architecture, Function, Regulation, and Evolution of α-Glucans Metabolic Enzymes in Prokaryotes

Bacteria have acquired sophisticated mechanisms for assembling and disassembling polysaccharides of different chemistry. α-d-Glucose homopolysaccharides, so-called α-glucans, are the most widespread polymers in nature being key components of microorganisms. Glycogen functions as an intracellular energy storage while some bacteria also produce extracellular assorted α-glucans. The classical bacterial glycogen metabolic pathway comprises the action of ADP-glucose pyrophosphorylase and glycogen synthase, whereas extracellular α-glucans are mostly related to peripheral enzymes dependent on sucrose. An alternative pathway of glycogen biosynthesis, operating via a maltose 1-phosphate polymerizing enzyme, displays an essential wiring with the trehalose metabolism to interconvert disaccharides into polysaccharides. Furthermore, some bacteria show a connection of intracellular glycogen metabolism with the genesis of extracellular capsular α-glucans, revealing a relationship between the storage and structural function of these compounds. Altogether, the current picture shows that bacteria have evolved an intricate α-glucan metabolism that ultimately relies on the evolution of a specific enzymatic machinery. The structural landscape of these enzymes exposes a limited number of core catalytic folds handling many different chemical reactions. In this Review, we present a rationale to explain how the chemical diversity of α-glucans emerged from these systems, highlighting the underlying structural evolution of the enzymes driving α-glucan bacterial metabolism.

Hot springs viruses at Yellowstone National Park have ancient origins and are adapted to thermophilic hosts

Geothermal springs house unicellular red algae in the class Cyanidiophyceae that dominate the microbial biomass at these sites. Little is known about host-virus interactions in these environments. We analyzed the virus community associated with red algal mats in three neighboring habitats (creek, endolithic, soil) at Lemonade Creek, Yellowstone National Park (YNP), USA. We find that despite proximity, each habitat houses a unique collection of viruses, with the giant viruses, Megaviricetes, dominant in all three. The early branching phylogenetic position of genes encoded on metagenome assembled virus genomes (vMAGs) suggests that the YNP lineages are of ancient origin and not due to multiple invasions from mesophilic habitats. The existence of genomic footprints of adaptation to thermophily in the vMAGs is consistent with this idea. The Cyanidiophyceae at geothermal sites originated ca. 1.5 Bya and are therefore relevant to understanding biotic interactions on the early Earth.

Expanded expression of pro-neurogenic factor SoxB1 during larval development of gastropod Lymnaea stagnalis suggests preadaptation to prolonged neurogenesis in Mollusca

Introduction

The remarkable diversity observed in the structure and development of the molluscan nervous system raises intriguing questions regarding the molecular mechanisms underlying neurogenesis in Mollusca. The expression of SoxB family transcription factors plays a pivotal role in neuronal development, thereby offering valuable insights into the strategies of neurogenesis.

Methods

In this study, we conducted gene expression analysis focusing on SoxB-family transcription factors during early neurogenesis in the gastropod Lymnaea stagnalis. We employed a combination of hybridization chain reaction in situ hybridization (HCR-ISH), immunocytochemistry, confocal microscopy, and cell proliferation assays to investigate the spatial and temporal expression patterns of LsSoxB1 and LsSoxB2 from the gastrula stage to hatching, with particular attention to the formation of central ring ganglia.

Results

Our investigation reveals that LsSoxB1 demonstrates expanded ectodermal expression from the gastrula to the hatching stage, whereas expression of LsSoxB2 in the ectoderm ceases by the veliger stage. LsSoxB1 is expressed in the ectoderm of the head, foot, and visceral complex, as well as in forming ganglia and sensory cells. Conversely, LsSoxB2 is mostly restricted to the subepithelial layer and forming ganglia cells during metamorphosis. Proliferation assays indicate a uniform distribution of dividing cells in the ectoderm across all developmental stages, suggesting the absence of distinct neurogenic zones with increased proliferation in gastropods.

Discussion

Our findings reveal a spatially and temporally extended pattern of SoxB1 expression in a gastropod representative compared to other lophotrochozoan species. This prolonged and widespread expression of SoxB genes may be interpreted as a form of transcriptional neoteny, representing a preadaptation to prolonged neurogenesis. Consequently, it could contribute to the diversification of nervous systems in gastropods and lead to an increase in the complexity of the central nervous system in Mollusca.

Rich microbial and depolymerising diversity in Antarctic krill gut

With almost a quadrillion individuals, the Antarctic krill processes five million tons of organic carbon every day during austral summer. This high carbon flux requires a broad range of hydrolytic enzymes to decompose the diverse food-derived biopolymers. While krill itself possesses numerous such enzymes, it is unclear, to what extent the endogenous microbiota contribute to the hydrolytic potential of the gut environment. Here we applied amplicon sequencing, shotgun metagenomics, cultivation, and physiological assays to characterize the krill gut microbiota. The broad bacterial diversity (273 families, 919 genera, and 2,309 species) also included a complex potentially anaerobic sub-community. Plate-based assays with 198 isolated pure cultures revealed widespread capacities to utilize lipids (e.g., tributyrin), followed by proteins (casein) and to a lesser extent by polysaccharides (e.g., alginate and chitin). While most isolates affiliated with the genera Pseudoalteromonas and Psychrobacter, also Rubritalea spp. (Verrucomicrobia) were observed. The krill gut microbiota growing on marine broth agar plates possess 13,012 predicted hydrolyses; 15-fold more than previously predicted from a transcriptome-proteome compendium of krill. Cultivation-independent and -dependent approaches indicated members of the families Flavobacteriaceae and Pseudoalteromonadaceae to dominate the capacities for lipid/protein hydrolysis and to provide a plethora of carbohydrate-active enzymes, sulfatases, and laminarin- or porphyrin-depolymerizing hydrolases. Notably, also the potential to hydrolyze plastics such as polyethylene terephthalate and polylactatide was observed, affiliating mostly with Moraxellaceae. Overall, this study shows extensive microbial diversity in the krill gut, and suggests that the microbiota likely play a significant role in the nutrient acquisition of the krill by enriching its hydrolytic enzyme repertoire.IMPORTANCEThe Antarctic krill (Euphausia superba) is a keystone species of the Antarctic marine food web, connecting the productivity of phyto- and zooplankton with the nutrition of the higher trophic levels. Accordingly, krill significantly contributes to biomass turnover, requiring the decomposition of seasonally varying plankton-derived biopolymers. This study highlights the likely role of the krill gut microbiota in this ecosystem function by revealing the great number of diverse hydrolases that microbes contribute to the krill gut environment. The here resolved repertoire of hydrolytic enzymes could contribute to the overall nutritional resilience of krill and to the general organic matter cycling under changing environmental conditions in the Antarctic sea water. Furthermore, the krill gut microbiome could serve as a valuable resource of cold-adapted hydrolytic enzymes for diverse biotechnological applications.

Mitochondrial genome comparison and phylogenetic position of Fannia pusio among the Calyptratae flies

Fannia pusio, the chicken dung fly species, remains unexplored despite its forensic, sanitary, and veterinary importance in the Nearctic and Neotropical regions. In this study, we obtained the complete mitochondrial genome of Fannia pusio for the first time using next-generation sequencing. We compared it with previously published mitogenomes of the genus from the Palearctic region, and its phylogenetic position was studied based on the concatenated protein-coding genes (PCGs) dataset of Calyptratae flies. The circular mitochondrial genome of F. pusio is 16,176 bp in length, with a high A + T content (78.3%), whose gene synteny, codon usage analysis, and amino acid frequency are similar to previously reported Fannia mitogenomes. All PCGs underwent purifying selection except the nad2 gene. Interspecific K2P distances of PCGs of Fannia yielded an average of 12.4% (8.1%-21.1%). The Fannia genus is monophyletic and closely related to Muscidae based on molecular data. Further taxonomic sampling is required to deep into the phylogenetic relationships of the originally proposed species-groups and subgroups within the genus. These results provide a valuable dataset for studying the mitochondrial genome evolution and a resource for the taxonomy and systematics of Fannia.

Exploring Andean High-Altitude Lake Extremophiles through Advanced Proteotyping

Quickly identifying and characterizing isolates from extreme environments is currently challenging while very important to explore the Earth's biodiversity. As these isolates may, in principle, be distantly related to known species, techniques are needed to reliably identify the branch of life to which they belong. Proteotyping these environmental isolates by tandem mass spectrometry offers a rapid and cost-effective option for their identification using their peptide profiles. In this study, we document the first high-throughput proteotyping approach for environmental extremophilic and halophilic isolates. Microorganisms were isolated from samples originating from high-altitude Andean lakes (3700-4300 m a.s.l.) in the Chilean Altiplano, which represent environments on Earth that resemble conditions on other planets. A total of 66 microorganisms were cultivated and identified by proteotyping and 16S rRNA gene amplicon sequencing. Both the approaches revealed the same genus identification for all isolates except for three isolates possibly representing not yet taxonomically characterized organisms based on their peptidomes. Proteotyping was able to indicate the presence of two potentially new genera from the families of Paracoccaceae and Chromatiaceae/Alteromonadaceae, which have been overlooked by 16S rRNA amplicon sequencing approach only. The paper highlights that proteotyping has the potential to discover undescribed microorganisms from extreme environments.

Exploring the microbial diversity of novel misos with metagenomics

Interest in fermented foods, especially plant-based ones, has increased considerably in the last decade. Miso-a Japanese paste traditionally fermented with soybeans, salt, and kōji (Aspergillus oryzae grown on grains or beans)-has gained attention among chefs for its rich flavour and versatility. Some chefs have even been experimenting with making novel misos with untraditional substrates to create new flavours. Such novel fermented foods also offer new scientific opportunities. To explore the microbial diversity of these new traditional foods, we sampled six misos made by the team at a leading restaurant called Noma in Copenhagen (Denmark), using yellow peas (including a nixtamalised treatment), lupin seeds, Swedish Vreta peas, grey peas, and Gotland lentils as substrates. All misos were made with the same recipe and fermented for 3 months at 28 °C. Samples were collected at the end of fermentation for subsequent shotgun metagenomic sequencing and a genome-resolved metagenomic analysis. The taxonomic profile of the samples revealed the presence of kōji mould (A. oryzae) and Bacillus amyloliquefaciens in all misos. Various species of the genera Latilactobacillus, Lactiplantibacillus, Pediococcus and Staphylococcus were also detected. The Metagenome-Assembled Genomes (MAGs) revealed genomic sequences belonging to 12 different species and functional analyses of these MAGs were performed. Notably, we detected the presence of Exiguobacterium-the first reported instance of the genus in miso-and Average Nucleotide Identity (ANI) analyses suggest a potentially new species. We hope these results will improve the scientific literature on misos and contribute to developing novel fermented plant-based foods.

Taxonomic and functional assembly cues enrich the endophytic tobacco microbiota across epiphytic compartments

The plant microbiome plays a critical role in plant growth, development, and health, with endophytes being recognized as essential members due to their close interactions with host plants. However, knowledge gaps remain in understanding the mechanisms driving the colonization and establishment of endophytic communities. To address this issue, we investigated the microbiota of tobacco roots and leaves, including both epiphytic and endophytic microorganisms. We found that Actinobacteria and Alphaproteobacteria were significantly enriched in the root endosphere. Additionally, we identified higher abundances of functional traits involved in antibiotic synthesis, plant cell wall degradation, iron metabolism, secretion systems, and nicotine degradation enzymes in the endosphere. We further studied metagenome-assembled genomes from the rhizosphere and root endosphere, revealing a greater diversity of secondary metabolites in bacteria within the root endosphere. Together, this study provides insights into the taxonomic and functional assembly cues that may contribute to shaping the endophytic plant microbiota.IMPORTANCEThe presence of diverse microorganisms within plant tissues under natural conditions is a well-established fact. However, due to the plant immune system's barrier and the unique microhabitat of the plant interior, it remains unclear what specific characteristics bacteria require to successfully colonize and thrive in the plant endosphere. Recognizing the significance of unraveling these functional features, our study focused on investigating the enriched traits in the endophytic microbiota compared to the epiphytes. Through our research, we have successfully identified the taxonomic and functional assembly cues that drive the enrichment of the endophytic microbiota across the epiphytic compartments. These findings shed new light on the intricate mechanisms of endophyte colonization, thereby deepening our understanding of plant-microbe interactions and paving the way for further advancements in microbiome manipulation.

Protein structures unravel the signatures and patterns of deep time evolution

The formulation and testing of hypotheses using 'big biology data' often lie at the interface of computational biology and structural biology. The Protein Data Bank (PDB), which was established about 50 years ago, catalogs three-dimensional (3D) shapes of organic macromolecules and showcases a structural view of biology. The comparative analysis of the structures of homologs, particularly of proteins, from different species has significantly improved the in-depth analyses of molecular and cell biological questions. In addition, computational tools that were developed to analyze the 'protein universe' are providing the means for efficient resolution of longstanding debates in cell and molecular evolution. In celebrating the golden jubilee of the PDB, much has been written about the transformative impact of PDB on a broad range of fields of scientific inquiry and how structural biology transformed the study of the fundamental processes of life. Yet, the transforming influence of PDB on one field of inquiry of fundamental interest-the reconstruction of the distant biological past-has gone almost unnoticed. Here, I discuss the recent advances to highlight how insights and tools of structural biology are bearing on the data required for the empirical resolution of vigorously debated and apparently contradicting hypotheses in evolutionary biology. Specifically, I show that evolutionary characters defined by protein structure are superior compared to conventional sequence characters for reliable, data-driven resolution of competing hypotheses about the origins of the major clades of life and evolutionary relationship among those clades. Since the better quality data unequivocally support two primary domains of life, it is imperative that the primary classification of life be revised accordingly.

Genetic Diversity of Salmonella enterica subsp. enterica Serovar Enteritidis from Human and Non-Human Sources in Portugal

Salmonella enterica subsp. enterica serovar Enteritidis (S. Enteritidis) is one of the leading causes of foodborne infections associated with broilers and laying hens. Portugal has had the lowest notification rates of salmonellosis in recent years, due to the vaccinations of layer and breeder flocks and strict compliance with biosecurity measures. However, data about the genetic diversity of S. Enteritidis in Portugal are scarce. In this study, 102 S. Enteritidis isolates selected from human (n = 63) and non-human sources (n = 39) were characterized by serotyping, antimicrobial susceptibility, and whole genome sequencing. The S. Enteritidis population was mainly resistant to fluoroquinolones, and a sole isolate showed resistance to extended-spectrum cephalosporins. ST11 was the most frequent sequence type, and three novel STs from human isolates (ST9236, ST4457, and ST9995) were assigned. Several Salmonella pathogenic islands (SPI) and Putative SPI were present in the genomes, namely SPI-1, 2, 3, 4, 5, 9, 10, 12, 13, and 14, C63PI, CS54_island, and 170 virulence genes were identified. The phylogenetic analysis revealed that strains from Portugal are genetically heterogeneous regarding sample type, collection date, and genetic content. This study increases the available data, essential to a better characterization of strains in a global context.

Inference of phylogenetic trees directly from raw sequencing reads using Read2Tree

Current methods for inference of phylogenetic trees require running complex pipelines at substantial computational and labor costs, with additional constraints in sequencing coverage, assembly and annotation quality, especially for large datasets. To overcome these challenges, we present Read2Tree, which directly processes raw sequencing reads into groups of corresponding genes and bypasses traditional steps in phylogeny inference, such as genome assembly, annotation and all-versus-all sequence comparisons, while retaining accuracy. In a benchmark encompassing a broad variety of datasets, Read2Tree is 10-100 times faster than assembly-based approaches and in most cases more accurate-the exception being when sequencing coverage is high and reference species very distant. Here, to illustrate the broad applicability of the tool, we reconstruct a yeast tree of life of 435 species spanning 590 million years of evolution. We also apply Read2Tree to >10,000 Coronaviridae samples, accurately classifying highly diverse animal samples and near-identical severe acute respiratory syndrome coronavirus 2 sequences on a single tree. The speed, accuracy and versatility of Read2Tree enable comparative genomics at scale.

Gene Transfer-Based Phylogenetics: Analytical Expressions and Additivity via Birth-Death Theory

The genomic era has opened up vast opportunities in molecular systematics, one of which is deciphering the evolutionary history in fine detail. Under this mass of data, analyzing the point mutations of standard markers is often too crude and slow for fine-scale phylogenetics. Nevertheless, genome dynamics (GD) events provide alternative, often richer information. The synteny index (SI) between a pair of genomes combines gene order and gene content information, allowing the comparison of genomes of unequal gene content, together with order considerations of their common genes. Recently, genome dynamics has been modeled as a continuous-time Markov process, and gene distance in the genome as a birth-death-immigration process. Nevertheless, due to complexities arising in this setting, no precise and provably consistent estimators could be derived, resulting in heuristic solutions. Here, we extend this modeling approach by using techniques from birth-death theory to derive explicit expressions of the system's probabilistic dynamics in the form of rational functions of the model parameters. This, in turn, allows us to infer analytically accurate distances between organisms based on their SI. Subsequently, we establish additivity of this estimated evolutionary distance (a desirable property yielding phylogenetic consistency). Applying the new measure in simulation studies shows that it provides accurate results in realistic settings and even under model extensions such as gene gain/loss or over a tree structure. In the real-data realm, we applied the new formulation to unique data structure that we constructed-the ordered orthology DB-based on a new version of the EggNOG database, to construct a tree with more than 4.5K taxa. To the best of our knowledge, this is the largest gene-order-based tree constructed and it overcomes shortcomings found in previous approaches. Constructing a GD-based tree allows to confirm and contrast findings based on other phylogenetic approaches, as we show.

Ecological and molecular associations between neotropical wild felids and Taenia (Cestoda: Taeniidae) in the Atlantic Forest: a new report for Taenia omissa

Ecological associations between wild felids and parasites from the Taeniidae family are related to predator-prey interactions, where felids act as definitive hosts while their prey, herbivores and/or omnivores, act as intermediate hosts. In the Atlantic Forest, six neotropical felid species coexist in sympatry, but the ecological parasite-host interactions remain poorly studied. Taenia omissa is a tapeworm that parasitizes cougars (Puma concolor) as its only definitive host and their ungulate prey as intermediate hosts. The aim of this study was to identify tapeworms present in road-killed fauna using both molecular and morphological characteristics and their predator-prey relationship. Adult tapeworms found in a cougar, a jaguarundi (Herpailurus yagouaroundi), and two ocelots (Leopardus pardalis); and metacestodes found in a red brocket deer (Mazama americana) and a wild guinea pig (Cavia aperea) were analyzed. Through morphological analysis of rostellar hooks and molecular analysis of the mitochondrial genetic marker cox1, Taenia omissa adult individuals were identified in the cougar, and metacestodes in the red brocket deer, proving the existence of a full host-parasite life cycle in the Atlantic Forest region. This new report reveals the southernmost record of T. omissa and broadens its geographic distribution. In addition, isolates of the Taenia genus divergent from those described so far in molecular databases were reported and suggested a wild cycle that involves the jaguarundi and agouti (Dasyprocta asarae) as definitive and intermediate hosts, respectively. These results highlight the complexity of the tapeworm population in the region and the need to study it with both morphological and molecular approaches.

Sequence-discrete species for prokaryotes and other microbes: A historical perspective and pending issues

Whether prokaryotes, and other microorganisms, form distinct clusters that can be recognized as species remains an issue of paramount theoretical as well as practical consequence in identifying, regulating, and communicating about these organisms. In the past decade, comparisons of thousands of genomes of isolates and hundreds of metagenomes have shown that prokaryotic diversity may be predominantly organized in such sequence-discrete clusters, albeit organisms of intermediate relatedness between the identified clusters are also frequently found. Accumulating evidence suggests, however, that the latter "intermediate" organisms show enough ecological and/or functional distinctiveness to be considered different species. Notably, the area of discontinuity between clusters often-but not always-appears to be around 85%-95% genome-average nucleotide identity, consistently among different taxa. More recent studies have revealed remarkably similar diversity patterns for viruses and microbial eukaryotes as well. This high consistency across taxa implies a specific mechanistic process that underlies the maintenance of the clusters. The underlying mechanism may be a substantial reduction in the efficiency of homologous recombination, which mediates (successful) horizontal gene transfer, around 95% nucleotide identity. Deviations from the 95% threshold (e.g., species showing lower intraspecies diversity) may be caused by ecological differentiation that imposes barriers to otherwise frequent gene transfer. While this hypothesis that clusters are driven by ecological differentiation coupled to recombination frequency (i.e., higher recombination within vs. between groups) is appealing, the supporting evidence remains anecdotal. The data needed to rigorously test the hypothesis toward advancing the species concept are also outlined.

Phylogenetics and biomineralization of a novel magnetotactic Gammaproteobacterium from a freshwater lake in Beijing, China

Magnetotactic bacteria (MTB) have the remarkable capability of producing intracellularly membrane-enveloped magnetic nanocrystals (i.e. magnetosomes) and swimming along geomagnetic field lines. Despite more than 50 years of research, bacterial diversity and magnetosome biomineralization within MTB are relatively less known in the Gammaproteobacteria class than other groups. This is incompatible with the status of Gammaproteobacteria as the most diverse class of gram-negative bacteria with a number of ecologically important bacteria. Here, we identify a novel MTB strain YYHR-1 affiliated with the Gammaproteobacteria class of the Pseudomonadota phylum from a freshwater lake. In YYHR-1, most magnetosome crystals are organized into a long chain aligned along the cell long axis; unusually, a few small superparamagnetic crystals are located at the side of the chain, off the main chain axis. Micromagnetic simulations indicate that magnetostatic interactions among adjacent crystals within a chain reduce the Gibbs energy to enhance chain stability. Genomic analysis suggests that duplication of magnetosome gene clusters may result in off-chain magnetosomes formation. By integrating available genomic data from Gammaproteobacteria, the phylogenetic position of MTB in this class is reassigned here. Our new findings expand knowledge about MTB diversity and magnetosome biomineralization, and deepen understanding of the phylogenetics of the Gammaproteobacteria.

Bacteriophages suppress CRISPR-Cas immunity using RNA-based anti-CRISPRs

Many bacteria use CRISPR-Cas systems to combat mobile genetic elements, such as bacteriophages and plasmids1. In turn, these invasive elements have evolved anti-CRISPR proteins to block host immunity2,3. Here we unveil a distinct type of CRISPR-Cas Inhibition strategy that is based on small non-coding RNA anti-CRISPRs (Racrs). Racrs mimic the repeats found in CRISPR arrays and are encoded in viral genomes as solitary repeat units4. We show that a prophage-encoded Racr strongly inhibits the type I-F CRISPR-Cas system by interacting specifically with Cas6f and Cas7f, resulting in the formation of an aberrant Cas subcomplex. We identified Racr candidates for almost all CRISPR-Cas types encoded by a diverse range of viruses and plasmids, often in the genetic context of other anti-CRISPR genes5. Functional testing of nine candidates spanning the two CRISPR-Cas classes confirmed their strong immune inhibitory function. Our results demonstrate that molecular mimicry of CRISPR repeats is a widespread anti-CRISPR strategy, which opens the door to potential biotechnological applications6.

INTERPIN: A repository for intrinsic transcription termination hairpins in bacteria

The large-scale detection of putative intrinsic transcription terminators is limited to only a few bacteria currently. We discovered a group of hairpins, called cluster hairpins, present within 15 nucleotides from each other. These are expected to work in tandem to cause intrinsic transcription termination (ITT), while the single hairpin can do the same alone. Therefore, exploring these ITT sites and the hairpins across bacterial genomes becomes highly desirable. INTERPIN is the largest archived collection of in silico inferred ITT hairpins in bacteria, covering 12745 bacterial genomes and encompassing ten bacterial phyla for ∼25 million hairpins. Users can obtain details on operons, individual cluster, and single ITT hairpins that were screened therein. Integrated Genome Viewer (IGV) software interactively visualizes hairpin secondary and tertiary structures in the genomic context. We also discuss statistics for the occurrence of cluster or single hairpins and other termination alternatives while showing the validation of predicted hairpins against in vivo detected hairpins. The database is freely available at http://pallab.cds.iisc.ac.in/INTERPIN/. INTERPIN (database and software) can make predictions for both AT and GC-rich genomes, which has not been achieved by any other program so far. It can also be used to improve genome annotation as well as to get predictions to improve the understanding of the ITT pathway by further analysis.

Four functional profiles for fibre and mucin metabolism in the human gut microbiome

BACKGROUND

With the emergence of metagenomic data, multiple links between the gut microbiome and the host health have been shown. Deciphering these complex interactions require evolved analysis methods focusing on the microbial ecosystem functions. Despite the fact that host or diet-derived fibres are the most abundant nutrients available in the gut, the presence of distinct functional traits regarding fibre and mucin hydrolysis, fermentation and hydrogenotrophic processes has never been investigated.

RESULTS

After manually selecting 91 KEGG orthologies and 33 glycoside hydrolases further aggregated in 101 functional descriptors representative of fibre and mucin degradation pathways in the gut microbiome, we used nonnegative matrix factorization to mine metagenomic datasets. Four distinct metabolic profiles were further identified on a training set of 1153 samples, thoroughly validated on a large database of 2571 unseen samples from 5 external metagenomic cohorts and confirmed with metatranscriptomic data. Profiles 1 and 2 are the main contributors to the fibre-degradation-related metagenome: they present contrasted involvement in fibre degradation and sugar metabolism and are differentially linked to dysbiosis, metabolic disease and inflammation. Profile 1 takes over Profile 2 in healthy samples, and unbalance of these profiles characterize dysbiotic samples. Furthermore, high fibre diet favours a healthy balance between profiles 1 and profile 2. Profile 3 takes over profile 2 during Crohn's disease, inducing functional reorientations towards unusual metabolism such as fucose and H2S degradation or propionate, acetone and butanediol production. Profile 4 gathers under-represented functions, like methanogenesis. Two taxonomic makes up of the profiles were investigated, using either the covariation of 203 prevalent genomes or metagenomic species, both providing consistent results in line with their functional characteristics. This taxonomic characterization showed that profiles 1 and 2 were respectively mainly composed of bacteria from the phyla Bacteroidetes and Firmicutes while profile 3 is representative of Proteobacteria and profile 4 of methanogens.

CONCLUSIONS

Integrating anaerobic microbiology knowledge with statistical learning can narrow down the metagenomic analysis to investigate functional profiles. Applying this approach to fibre degradation in the gut ended with 4 distinct functional profiles that can be easily monitored as markers of diet, dysbiosis, inflammation and disease. Video Abstract.

Bridging Evolutionary History and Conservation of New World Vultures

The New World Vultures (Cathartidae) include seven species of obligate scavengers that, despite their ecological relevance, present critical information gaps around their evolutionary history and conservation. Insights into their phylogenetic relationships in recent years has enabled the addressing of such information gaps through approaches based on phylogeny. We reconstructed the ancestral area in America of the current species using two regionalization schemes and methods: Biogeography with Bayesian Evolutionary Analysis (BioGeoBears) and Bayesian Binary Model-Monte Carlo Markov Chains (BBM-MCMC). Then, we identified the priority species and areas for conservation by means of the Evolutionary Distinctiveness index (ED), as a proxy of the uniqueness of species according to phylogeny, and the Global Endangerment index (GE), mapping phylogenetic diversity. We found that the ancestral area of New World Vultures in America corresponds to South America, with dispersal processes that led to a recolonization of North America by Coragyps atratus, Gymnogyps californianus and Cathartes aura. We identified the Black Vulture, G. californianus and Vultur gryphus as priority species based on ED and "Evolutionary Distinct Globally Endangered" (EDGE) indexes, and the lowlands of Amazon River basin and the Orinoco basin and some tributaries areas of the Guiana Shield were identified as the priority areas when mapping the phylogenetic diversity. This study highlights the importance of filling knowledge gaps of species of conservation concern through the integration of evolutionary and ecological information and tools and, thus, developing adequate strategies to enhance the preservation of these species in the face of the current loss of biodiversity.

Intrinsically disordered CsoS2 acts as a general molecular thread for α-carboxysome shell assembly

Carboxysomes are a paradigm of self-assembling proteinaceous organelles found in nature, offering compartmentalisation of enzymes and pathways to enhance carbon fixation. In α-carboxysomes, the disordered linker protein CsoS2 plays an essential role in carboxysome assembly and Rubisco encapsulation. Its mechanism of action, however, is not fully understood. Here we synthetically engineer α-carboxysome shells using minimal shell components and determine cryoEM structures of these to decipher the principle of shell assembly and encapsulation. The structures reveal that the intrinsically disordered CsoS2 C-terminus is well-structured and acts as a universal "molecular thread" stitching through multiple shell protein interfaces. We further uncover in CsoS2 a highly conserved repetitive key interaction motif, [IV]TG, which is critical to the shell assembly and architecture. Our study provides a general mechanism for the CsoS2-governed carboxysome shell assembly and cargo encapsulation and further advances synthetic engineering of carboxysomes for diverse biotechnological applications.

Concatenated 16S rRNA sequence analysis improves bacterial taxonomy

Background: Microscopic, biochemical, molecular, and computer-based approaches are extensively used to identify and classify bacterial populations. Advances in DNA sequencing and bioinformatics workflows have facilitated sophisticated genome-based methods for microbial taxonomy although sequencing of the 16S rRNA gene is widely employed to identify and classify bacterial communities as a cost-effective and single-gene approach. However, the 16S rRNA sequence-based species identification accuracy is limited because of the occurrence of multiple copies of the 16S rRNA gene and higher sequence identity between closely related species. The availability of the genomes of several bacterial species provided an opportunity to develop comprehensive species-specific 16S rRNA reference libraries. Methods: Sequences of the 16S rRNA genes were retrieved from the whole genomes available in the Genome databases. With defined criteria, four 16S rRNA gene copy variants were concatenated to develop a species-specific reference library. The sequence similarity search was performed with a web-based BLAST program, and MEGA software was used to construct the phylogenetic tree. Results: Using this approach, species-specific 16S rRNA gene libraries were developed for four closely related  Streptococcus species ( S. gordonii,  S. mitis,  S. oralis, and  S. pneumoniae). Sequence similarity and phylogenetic analysis using concatenated 16S rRNA copies yielded better resolution than single gene copy approaches. Conclusions: The approach is very effective in classifying genetically closely related bacterial species and may reduce misclassification of bacterial species and genome assemblies.

Small secreted peptides (SSPs) in tomato and their potential roles in drought stress response

Tomato (Solanum lycopersicum) is one of the most important vegetable crops in the world and abiotic stresses often cause serious problems in tomato production. It is thus important to identify new regulators in stress response and to devise new approaches to promote stress tolerance in tomato. Previous studies have shown that small secreted peptides (SSPs) are important signal molecules regulating plant growth and stress response by mediating intercellular communication. However, little is known about tomato SSPs, especially their roles in responding to abiotic stresses. Here we report the identification of 1,050 putative SSPs in the tomato genome, 557 of which were classified into 38 known SSP families based on their conserved domains. GO and transcriptome analyses revealed that a large proportion of SlSSPs might be involved in abiotic stress response. Further analysis indicated that stress response related cis-elements were present on the SlCEP promotors and a number of SlCEPs were significantly upregulated by drought treatments. Among the drought-inducible SlCEPs, SlCEP10 and SlCEP11b were selected for further analysis via exogenous application of synthetic peptides. The results showed that treatments with both SlCEP10 and SlCEP11b peptides enhanced tomato drought stress tolerance, indicating the potential roles of SlSSPs in abiotic stress response.

Complete genome sequence of a novel bacteriophage vB_Pci_PCMW57 infecting phytobacteria pseudomonas cichorii

BACKGROUND

A novel virulent bacteriophage infecting phytobacteria Pseudomonas cichorii (P. cichorii) was isolated from leafy vegetables in Brazil. P. cichorii is a Gram-negative soil phytobacterium, the causal agent of a number of economically important plant diseases worldwide.

METHODS AND RESULTS

In this study, a new phage specific for P. cichorii was isolated from solid samples (lettuce, chicory and cabbage), designated vB_Pci_PCMW57. Electron microscopy revealed a small virion (~ 50-nm-diameter icosahedral capsid) with a short, non-contractile tail. The genome of vB_Pci_PCMW57 is 40,117 bp in size, with a GC content of 57.6% and encodes 49 open reading frames. The phage is genetically similar to P. syringae phages Pst_GM1 and Pst_GIL1, and the P. fluorescens phages WRT and KNP. According to electron microscopy and whole-genome sequence analysis, vB_Pci_PCMW57 should be classified as a Caudoviticetes, family Autographiviridae, subfamily Studiervirinae.

CONCLUSIONS

The complete phage genome was annotated, and the sequence identity of the virus with other Pseudomonas viruses was higher than 95%. To our knowledge, this is the first report of a bacteriophage infecting Pseudomonas cichorii.

New viruses of Cladosporium sp. expand considerably the taxonomic structure of Gammapartitivirus genus

Despite the fact that Cladosporium sp. are ubiquitous fungi, their viromes have been little studied. By analysing a collection of Cladosporium fungi, two new partitiviruses named Cladosporium cladosporioides partitivirus 1 (CcPV1) and Cladosporium cladosporioides partitivirus 2 (CcPV2) co-infecting a strain of Cladosporium cladosporioides were identified. Their complete genome consists of two monocistronic dsRNA segments (RNA1 and RNA2) with a high percentage of pairwise identity on 5' and 3' end. The RNA directed RNA polymerase (RdRp) of both viruses and the capsid protein (CP) of CcPV1 display the classic characteristics required for their assignment to the Gammapartitivirus genus. In contrast, CcPV2 RNA2 encodes for a 41 KDa CP that is unusually smaller when aligned to CPs of other viruses classified in this genus. The structural role of this protein is confirmed by electrophoresis on acrylamide gel of purified viral particles. Despite the low percentage of identity between the capsid proteins of CcPV1 and CcPV2, their three-dimensional structures predicted by AlphaFold2 show strong similarities and confirm functional proximity. Fifteen similar viral sequences of unknown function were annotated using the CcPV2 CP sequence. The phylogeny of the CP was highly consistent with the phylogeny of their corresponding RdRp, supporting the organization of Gammapartitiviruses into three distinct clades despite stretching the current demarcation criteria. It is proposed that a new subgenus be created within the genus Gammapartitivirus for this new group.

Genome-wide phylodynamic approach reveals the epidemic dynamics of the main Mycoplasma bovis subtype circulating in France

Mycoplasma bovis is a major aetiological agent of bovine respiratory disease worldwide. Genome-based analyses are increasingly being used to monitor the genetic diversity and global distribution of M. bovis, complementing existing subtyping schemes based on locus sequencing. However, these analyses have so far provided limited information on the spatiotemporal and population dynamics of circulating subtypes. Here we applied a genome-wide phylodynamic approach to explore the epidemic dynamics of 88 French M. bovis strains collected between 2000 and 2019 in France and belonging to the currently dominant polC subtype 2 (st2). A strong molecular clock signal detected in the genomic data enabled robust phylodynamic inferences, which estimated that the M. bovis st2 population in France is composed of two lineages that successively emerged from independent introductions of international strains. The first lineage appeared around 2000 and supplanted the previously established antimicrobial-susceptible polC subtype 1. The second lineage, which is likely more transmissible, progressively replaced the first M. bovis st2 lineage population from 2005 onward and became predominant after 2010. Analyses also showed a brief decline in this second M. bovis st2 lineage population in around 2011, possibly due to the challenge from the concurrent emergence of M. bovis polC subtype 3 in France. Finally, we identified non-synonymous mutations in genes associated with lineages, which raises prospects for identifying new surveillance molecular markers. A genome-wide phylodynamic approach provides valuable resources for monitoring the evolution and epidemic dynamics of circulating M. bovis subtypes, and may prove critical for developing more effective surveillance systems and disease control strategies.

UACG: Up-to-Date Archaeal Core Genes and Software for Phylogenomic Tree Reconstruction

In the post-genomic era, phylogenomics is a powerful and routinely-used tool to discover evolutionary relationships between microorganisms. Inferring phylogenomic trees by concatenating core gene sequences into a supermatrix is the standard method. The previously released up-to-date bacterial core gene (UBCG) tool provides a pipeline to infer phylogenomic trees using single-copy core genes for the Bacteria domain. In this study, we established up-to-date archaeal core gene (UACG), comprising 128 genes suitable for inferring archaeal phylogenomic trees. To test the gene set, we selected the Haloarcula genus and scrutinized its phylogeny. The phylogeny inferred using the UACG tool was consistent with the orthoANIu dendrogram, whereas the 16S rRNA gene phylogeny showed high intragenomic heterogeneity resulting in phylogenetic discrepancies. The software tool using the UACG set is available at https://www.ezbiocloud.net/tools/uacg .

Improving In Vitro Detection of Sensitization to Lipid Transfer Proteins: A New Molecular Multiplex IgE Assay

SCOPE

LTP-syndrome is characterized by sensitization (IgE) to multiple non-specific lipid transfer proteins (nsLTPs) with a variable clinical outcome. The treatment is primarily based on offending food avoidance. However, the determination of Pru p 3-specific IgE is currently the main diagnostic tool to assess sensitization to nsLTPs. Herein, the study evaluates improvement of LTP-syndrome diagnosis and clinical management using a new IgE multiplex-immunoblot assay with a high diversity of food nsLTPs.

METHODS AND RESULTS

An EUROLINE-LTP strip with 28 recombinant nsLTPs from 18 allergenic sources is designed. In total the study investigates 38 patients with LTP-syndrome and compares results from the nsLTPs (LTP-strip) with the respective food extracts of Prick-by-prick (PbP) testing. The agreement exceeds 70% for most nsLTPs, e.g., Pru p 3 (100%), Mal d 3 (97%), Pru av 3 (89%), Pha v 3 isoforms (87%/84%), Ara h 9 (82%), Cor a 8 (82%), and Jug r 3 (82%). The functionality and allergenic relevance of nine recombinant nsLTPs are proven by Basophil activation testing (BAT).

CONCLUSIONS

The new IgE multiplex-immunoblot nsLTP assay shows a good diagnostic performance allowing culprit food assessment. Negative results from LTP-strip may indicate potentially tolerable foods, improving diet intervention and patients' quality of life.

Fads2b Plays a Dominant Role in ∆6/∆5 Desaturation Activities Compared with Fads2a in Common Carp (Cyprinus carpio)

Highly unsaturated fatty acids (HUFAs) are essential for mammalian health, development and growth. However, most mammals, including humans, are incapable of synthesizing n-6 and n-3 HUFAs. Fish can convert C18 unsaturated fatty acids into n-6 and n-3 HUFAs via fatty acid desaturase (Fads), in which Fads2 is a key enzyme in HUFA biosynthesis. The allo-tetraploid common carp theoretically encode two duplicated fads2 genes. The expression patterns and desaturase functions of these two homologous genes are still unknown. In this study, the full length of the fads2a and fads2b were identified in common carp (Cyprinus carpio). Expression analyses indicate that both genes were mainly expressed in the liver and the expression of fads2b is higher than fads2a at different developmental stages in carp embryos. Heterogenous expression and 3D docking analyses suggested that Fads2b demonstrated stronger ∆6 and ∆5 desaturase activities than Fads2a. The core promotor regions of fads2a and fads2b were characterized and found to have different potential transcriptional binding sites. These results revealed the same desaturase functions, but different activities of two homologues of fasd2 genes in common carp. The data showed that fads2b played a more important role in HUFA synthesis through both expression and functional analyses.

RNA Encapsulation Mode and Evolutionary Insights from the Crystal Structure of Emaravirus Nucleoprotein

Enveloped RNA viruses are rare among plant viruses. Fimoviridae is a newly founded family of plant viruses within the Bunyavirales order that inflicts diverse crop losses worldwide. The fig mosaic virus (FMV), the representative member of the Fimoviridae family, was shown to be a causative agent for the fig mosaic disease. Like all bunyaviruses, FMV has a segmented, negative-sense, single-stranded RNA (ssRNA) genome that is encapsulated by the viral nucleoprotein (N). Here, we present high-resolution crystal structures of FMV N in its RNA-free and RNA-bound forms, revealing a "paper fortune teller" structural transition between the two states. The tightly packed tetramer of FNV N is similar to the structures of other N proteins of different members of the Bunyavirales order. In its RNA-bound form, the tetramer reorganizes to adopt a more open state that allows the accommodation of the RNA. Despite the low sequence similarity to N proteins of animal-infecting bunyaviruses, there is a striking structural resemblance between FMV N and nucleoproteins from members of the Peribunyaviridae, an animal-infecting family of viruses. This structural homology implies that enveloped plant viruses and animal-infecting viruses might have a common ancestor from which they diverged. IMPORTANCE Most insect-born viruses circulate within the Animalia kingdom, whereas plant-infecting RNA viruses are cross-kingdom pathogens. Many plant-infecting viruses cause devastating crop damage that leads to food security endangerment. The evolutionary crossroads of interkingdom circulation and infection are poorly understood. Thus, we took the structural approach to understand the similarities and differences between interkingdom-infecting viruses and viruses that circulate within one kingdom of life. Using our structures of FMV N in its free form and in complex with a single-stranded RNA (ssRNA), we dissected the mechanism by which FMV N binds to the RNA and revealed the conformational changes associated with the binding. The resemblance of our structure to N proteins from members of the Peribunyaviridae family and their recently published ribonucleoprotein (RNP) pseudoatomic resolution assembly model suggests that the FMV genome is similarly encapsulated. Thus, our finding unveils yet another bridge by which plant- and animal-infecting viruses are interconnected.

Physiological insight into the conserved properties of Caenorhabditis elegans acid-sensing degenerin/epithelial sodium channels

Acid-sensing ion channels (ASICs) are members of the diverse family of degenerin/epithelial sodium channels (DEG/ENaCs). They perform a wide range of physiological roles in healthy organisms, including in gut function and synaptic transmission, but also play important roles in disease, as acidosis is a hallmark of painful inflammatory and ischaemic conditions. We performed a screen for acid sensitivity on all 30 subunits of the Caenorhabditis elegans DEG/ENaC family using two-electrode voltage clamp in Xenopus oocytes. We found two groups of acid-sensitive DEG/ENaCs characterised by being either inhibited or activated by increasing proton concentrations. Three of these acid-sensitive C. elegans DEG/ENaCs were activated by acidic pH, making them functionally similar to the vertebrate ASICs. We also identified three new members of the acid-inhibited DEG/ENaC group, giving a total of seven additional acid-sensitive channels. We observed sensitivity to the anti-hypertensive drug amiloride as well as modulation by the trace element zinc. Acid-sensitive DEG/ENaCs were found to be expressed in both neurons and non-neuronal tissue, highlighting the likely functional diversity of these channels. Our findings provide a framework to exploit the C. elegans channels as models to study the function of these acid-sensing channels in vivo, as well as to study them as potential targets for anti-helminthic drugs. KEY POINTS: Acidosis plays many roles in healthy physiology, including synaptic transmission and gut function, but is also a key feature of inflammatory pain, ischaemia and many other conditions. Cells monitor acidosis of their surroundings via pH-sensing channels, including the acid-sensing ion channels (ASICs). These are members of the degenerin/epithelial sodium channel (DEG/ENaC) family, along with, as the name suggests, vertebrate ENaCs and degenerins of the roundworm Caenorhabditis elegans. By screening all 30 C. elegans DEG/ENaCs for pH dependence, we describe, for the first time, three acid-activated members, as well as three additional acid-inhibited channels. We surveyed both groups for sensitivity to amiloride and zinc; like their mammalian counterparts, their currents can be blocked, enhanced or unaffected by these modulators. Likewise, they exhibit diverse ion selectivity. Our findings underline the diversity of acid-sensitive DEG/ENaCs across species and provide a comparative resource for better understanding the molecular basis of their function.

Molecular recognition of SARS-CoV-2 spike protein with three essential partners: exploring possible immune escape mechanisms of viral mutants

OBJECTIVE

The COVID-19 epidemic is raging around the world, with the emergence of viral mutant strains such as Delta and Omicron, posing severe challenges to people's health and quality of life. A full understanding life cycle of the virus in host cells helps to reveal inactivation mechanism of antibody and provide inspiration for the development of a new-generation vaccines.

METHODS

In this work, molecular recognitions and conformational changes of SARS-CoV-2 spike protein mutants (i.e., Delta, Mu, and Omicron) and three essential partners (i.e., membrane receptor hACE2, protease TMPRSS2, and antibody C121) both were compared and analyzed using molecular simulations.

RESULTS

Water basin and binding free energy calculations both show that the three mutants possess higher affinity for hACE2 than WT, exhibiting stronger virus transmission. The descending order of cleavage ability by TMPRSS2 is Mu, Delta, Omicron, and WT, which is related to the new S1/S2 cutting site induced by transposition effect. The inefficient utilization of TMPRSS2 by Omicron is consistent with its primary entry into cells via the endosomal pathway. In addition, RBD-directed antibody C121 showed obvious resistance to Omicron, which may have originated from high fluctuation of approaching angles, high flexibility of I472-F490 loop, and reduced binding ability.

CONCLUSIONS

According to the overall characteristics of the three mutants, high infectivity, high immune escape, and low virulence may be the future evolutionary selection of SARS-CoV-2. In a word, this work not only proposes the possible resistance mechanism of SARS-CoV-2 mutants, but also provides theoretical guidance for the subsequent drug design against COVID-19 based on S protein structure.

β-Cryptoxanthin Production in Escherichia coli by Optimization of the Cytochrome P450 CYP97H1 Activity

Cytochromes P450, forming a superfamily of monooxygenases containing heme as a cofactor, show great versatility in substrate specificity. Metabolic engineering can take advantage of this feature to unlock novel metabolic pathways. However, the cytochromes P450 often show difficulty being expressed in a heterologous chassis. As a case study in the prokaryotic host Escherichia coli, the heterologous synthesis of β-cryptoxanthin was addressed. This carotenoid intermediate is difficult to produce, as its synthesis requires a monoterminal hydroxylation of β-carotene whereas most of the classic carotene hydroxylases are dihydroxylases. This study was focused on the optimization of the in vivo activity of CYP97H1, an original P450 β-carotene monohydroxylase. Engineering the N-terminal part of CYP97H1, identifying the matching redox partners, defining the optimal cellular background and adjusting the culture and induction conditions improved the production by 400 times compared to that of the initial strain, representing 2.7 mg/L β-cryptoxanthin and 20% of the total carotenoids produced.

Eurypsychrophilic acidophiles: From (meta)genomes to low-temperature biotechnologies

Low temperature and acidic environments encompass natural milieus such as acid rock drainage in Antarctica and anthropogenic sites including drained sulfidic sediments in Scandinavia. The microorganisms inhabiting these environments include polyextremophiles that are both extreme acidophiles (defined as having an optimum growth pH < 3), and eurypsychrophiles that grow at low temperatures down to approximately 4°C but have an optimum temperature for growth above 15°C. Eurypsychrophilic acidophiles have important roles in natural biogeochemical cycling on earth and potentially on other planetary bodies and moons along with biotechnological applications in, for instance, low-temperature metal dissolution from metal sulfides. Five low-temperature acidophiles are characterized, namely, Acidithiobacillus ferriphilus, Acidithiobacillus ferrivorans, Acidithiobacillus ferrooxidans, “Ferrovum myxofaciens,” and Alicyclobacillus disulfidooxidans, and their characteristics are reviewed. Our understanding of characterized and environmental eurypsychrophilic acidophiles has been accelerated by the application of “omics” techniques that have aided in revealing adaptations to low pH and temperature that can be synergistic, while other adaptations are potentially antagonistic. The lack of known acidophiles that exclusively grow below 15°C may be due to the antagonistic nature of adaptations in this polyextremophile. In conclusion, this review summarizes the knowledge of eurypsychrophilic acidophiles and places the information in evolutionary, environmental, biotechnological, and exobiology perspectives.

Aspergillus fumigatus Can Display Persistence to the Fungicidal Drug Voriconazole

Aspergillus fumigatus is a filamentous fungus that can infect the lungs of patients with immunosuppression and/or underlying lung diseases. The mortality associated with chronic and invasive aspergillosis infections remain very high, despite availability of antifungal treatments. In the last decade, there has been a worrisome emergence and spread of resistance to the first-line antifungals, the azoles. The mortality caused by resistant isolates is even higher, and patient management is complicated as the therapeutic options are reduced. Nevertheless, treatment failure is also common in patients infected with azole-susceptible isolates, which can be due to several non-mutually exclusive reasons, such as poor drug absorption. In addition, the phenomena of tolerance or persistence, where susceptible pathogens can survive the action of an antimicrobial for extended periods, have been associated with treatment failure in bacterial infections, and their occurrence in fungal infections already proposed. Here, we demonstrate that some isolates of A. fumigatus display persistence to voriconazole. A subpopulation of the persister isolates can survive for extended periods and even grow at low rates in the presence of supra-MIC of voriconazole and seemingly other azoles. Persistence cannot be eradicated with adjuvant drugs or antifungal combinations and seemed to reduce the efficacy of treatment for certain individuals in a Galleria mellonella model of infection. Furthermore, persistence implies a distinct transcriptional profile, demonstrating that it is an active response. We propose that azole persistence might be a relevant and underestimated factor that could influence the outcome of infection in human aspergillosis. IMPORTANCE The phenomena of antibacterial tolerance and persistence, where pathogenic microbes can survive for extended periods in the presence of cidal drug concentrations, have received significant attention in the last decade. Several mechanisms of action have been elucidated, and their relevance for treatment failure in bacterial infections demonstrated. In contrast, our knowledge of antifungal tolerance and, in particular, persistence is still very limited. In this study, we have characterized the response of the prominent fungal pathogen Aspergillus fumigatus to the first-line therapy antifungal voriconazole. We comprehensively show that some isolates display persistence to this fungicidal antifungal and propose various potential mechanisms of action. In addition, using an alternative model of infection, we provide initial evidence to suggest that persistence may cause treatment failure in some individuals. Therefore, we propose that azole persistence is an important factor to consider and further investigate in A. fumigatus.

MMR Deficiency Defines Distinct Molecular Subtype of Breast Cancer with Histone Proteomic Networks

Mismatch repair (MMR) alterations are important prognostic and predictive biomarkers in a variety of cancer subtypes, including colorectal and endometrial. However, in breast cancer (BC), the distinction and clinical significance of MMR are largely unknown. This may be due in part to the fact that genetic alterations in MMR genes are rare and only seen to occur in around 3% of BCs. In the present study, we analyzed TCGA data using a multi-sample protein–protein interaction (PPI) analysis tool, Proteinarium, and showed a distinct separation between specific MMR-deficient and -intact networks in a cohort of 994 BC patients. In the PPI networks specific to MMR deficiency, highly connected clusters of histone genes were identified. We also found the distribution of MMR-deficient BC to be more prevalent in HER2-enriched and triple-negative (TN) BC subtypes compared to luminal BCs. We recommend defining MMR-deficient BC by next-generation sequencing (NGS) when any somatic mutation is detected in one of the seven MMR genes.

Molecular Mechanisms of Resistance to Ceftazidime/Avibactam in Clinical Isolates of Enterobacterales and Pseudomonas aeruginosa in Latin American Hospitals

Ceftazidime-avibactam (CZA) is the combination of a third-generation cephalosporin and a new non-β-lactam β-lactamase inhibitor capable of inactivating class A, C, and some D β-lactamases. From a collection of 2,727 clinical isolates of Enterobacterales (n = 2,235) and P. aeruginosa (n = 492) that were collected between 2016 and 2017 from five Latin American countries, we investigated the molecular resistance mechanisms to CZA of 127 (18/2,235 [0.8%] Enterobacterales and 109/492 [22.1%] P. aeruginosa). First, by qPCR for the presence of genes encoding KPC, NDM, VIM, IMP, OXA-48-like, and SPM-1 carbapenemases, and second, by whole-genome sequencing (WGS). From the CZA-resistant isolates, MBL-encoding genes were detected in all 18 Enterobacterales and 42/109 P. aeruginosa isolates, explaining their resistant phenotype. Resistant isolates that yielded a negative qPCR result for any of the MBL encoding genes were subjected to WGS. The WGS analysis of the 67 remaining P. aeruginosa isolates showed mutations in genes previously associated with reduced susceptibility to CZA, such as those involved in the MexAB-OprM efflux pump and AmpC (PDC) hyperproduction, PoxB (bla_OXA-50-like), FtsI (PBP3), DacB (PBP4), and OprD. The results presented here offer a snapshot of the molecular epidemiological landscape for CZA resistance before the introduction of this antibiotic into the Latin American market. Therefore, these results serve as a valuable comparison tool to trace the evolution of the resistance to CZA in this carbapenemase-endemic geographical region. IMPORTANCE In this manuscript, we determine the molecular mechanisms of ceftazidime-avibactam resistance in Enterobacterales and P. aeruginosa isolates from five Latin American countries. Our results reveal a low rate of resistance to ceftazidime-avibactam among Enterobacterales; in contrast, resistance in P. aeruginosa has proven to be more complex, as it might involve multiple known and possibly unknown resistance mechanisms.

Lebanese fermented goat milk products: From tradition to meta-omics

Ambriss, Serdaleh and Labneh El Darff are traditional Lebanese products made from fermented goat's milk. A questionnaire completed by 50 producers of these products showed that they are prepared by periodic percolation either by milk or by Laban in amphora or goat skins during the lactation season. Production is carried out on a small scale and in a limited number of production units, often by elderly people, resulting in a real risk of disappearance of these products and loss of the corresponding microbial resources. In this study, 34 samples from 18 producers were characterized by culture-dependent and -independent analyses. The results obtained from these two methods were radically different, the latter revealing in Ambriss and Serdaleh the co-dominance of Lactobacillus kefiranofaciens, a fastidious-growing species, and Lactococcus lactis in a viable but not culturable state. Overall, their composition is reminiscent of kefir grains. Phylogenomic and functional analyses of the genomes of the key species Lb. kefiranofaciens have revealed differences from those found in kefir, particularly in their polysaccharide genes, which may explain the absence of grains. However, Labneh El Darff displayed a dominance of Lactobacillus delbrueckii, probably due to the addition of Laban. In addition, the study identified several zoonotic pathogens, including Streptococcus parasuis, which dominated in one sample. Metagenome-Assembled Genome (MAG) analysis indicated that this pathogen acquired lactose utilization genes through horizontal gene transfer. The contamination of the herd with Mycoplasmopsis agalactiae in the Chouf region was also revealed by MAG analysis of the Serdaleh samples. Antibiotic resistance genes were detected in most of the samples, particularly in the Serdaleh ones, where the dominant L. lactis strains possessed a plasmid with a multi-resistance island. Finally, this study paves the way for further analyses to shed light on the resilience of these ecosystems established in amphora or in goatskins and to improve hygiene practices for milk production.

Plant-Pathogenic Ralstonia Phylotypes Evolved Divergent Respiratory Strategies and Behaviors To Thrive in Xylem

Bacterial pathogens in the Ralstonia solanacearum species complex (RSSC) infect the water-transporting xylem vessels of plants, causing bacterial wilt disease. Strains in RSSC phylotypes I and III can reduce nitrate to dinitrogen via complete denitrification. The four-step denitrification pathway enables bacteria to use inorganic nitrogen species as terminal electron acceptors, supporting their growth in oxygen-limited environments such as biofilms or plant xylem. Reduction of nitrate, nitrite, and nitric oxide all contribute to the virulence of a model phylotype I strain. However, little is known about the physiological role of the last denitrification step, the reduction of nitrous oxide to dinitrogen by NosZ. We found that phylotypes I and III need NosZ for full virulence. However, strains in phylotypes II and IV are highly virulent despite lacking NosZ. The ability to respire by reducing nitrate to nitrous oxide does not greatly enhance the growth of phylotype II and IV strains. These partial denitrifying strains reach high cell densities during plant infection and cause typical wilt disease. However, unlike phylotype I and III strains, partial denitrifiers cannot grow well under anaerobic conditions or form thick biofilms in culture or in tomato xylem vessels. Furthermore, aerotaxis assays show that strains from different phylotypes have different oxygen and nitrate preferences. Together, these results indicate that the RSSC contains two subgroups that occupy the same habitat but have evolved divergent energy metabolism strategies to exploit distinct metabolic niches in the xylem. IMPORTANCE Plant-pathogenic Ralstonia spp. are a heterogeneous globally distributed group of bacteria that colonize plant xylem vessels. Ralstonia cells multiply rapidly in plants and obstruct water transport, causing fatal wilting and serious economic losses of many key food security crops. The virulence of these pathogens depends on their ability to grow to high cell densities in the low-oxygen xylem environment. Plant-pathogenic Ralstonia can use denitrifying respiration to generate ATP. The last denitrification step, nitrous oxide reduction by NosZ, contributes to energy production and virulence for only one of the three phytopathogenic Ralstonia species. These complete denitrifiers form thicker biofilms in culture and in tomato xylem, suggesting they are better adapted to hypoxic niches. Strains with partial denitrification physiology form less biofilm and are more often planktonic. They are nonetheless highly virulent. Thus, these closely related bacteria have adapted their core metabolic functions to exploit distinct microniches in the same habitat.

Native structure of mosquito salivary protein uncovers domains relevant to pathogen transmission

Female mosquitoes inject saliva into vertebrate hosts during blood feeding. This process transmits mosquito-borne human pathogens that collectively cause ~1,000,000 deaths/year. Among the most abundant and conserved proteins secreted by female salivary glands is a high-molecular weight protein called salivary gland surface protein 1 (SGS1) that facilitates pathogen transmission, but its mechanism remains elusive. Here, we determine the native structure of SGS1 by the cryoID approach, showing that the 3364 amino-acid protein has a Tc toxin-like Rhs/YD shell, four receptor domains, and a set of C-terminal daisy-chained helices. These helices are partially shielded inside the Rhs/YD shell and poised to transform into predicted transmembrane helices. This transformation, and the numerous receptor domains on the surface of SGS1, are likely key in facilitating sporozoite/arbovirus invasion into the salivary glands and manipulating the host's immune response.

Cytochrome P450 Surface Domains Prevent the β-Carotene Monohydroxylase CYP97H1 of Euglena gracilis from Acting as a Dihydroxylase

Molecular biodiversity results from branched metabolic pathways driven by enzymatic regioselectivities. An additional complexity occurs in metabolites with an internal structural symmetry, offering identical extremities to the enzymes. For example, in the terpene family, β-carotene presents two identical terminal closed-ring structures. Theses cycles can be hydroxylated by cytochrome P450s from the CYP97 family. Two sequential hydroxylations lead first to the formation of monohydroxylated β-cryptoxanthin and subsequently to that of dihydroxylated zeaxanthin. Among the CYP97 dihydroxylases, CYP97H1 from Euglena gracilis has been described as the only monohydroxylase. This study aims to determine which enzymatic domains are involved in this regioselectivity, conferring unique monohydroxylase activity on a substrate offering two identical sites for hydroxylation. We explored the effect of truncations, substitutions and domain swapping with other CYP97 members and found that CYP97H1 harbours a unique N-terminal globular domain. This CYP97H1 N-terminal domain harbours a hydrophobic patch at the entrance of the substrate channel, which is involved in the monohydroxylase activity of CYP97H1. This domain, at the surface of the enzyme, highlights the role of distal and non-catalytic domains in regulating enzyme specificity.

Ancestral Sequence Reconstruction of the Ribosomal Protein uS8 and Reduction of Amino Acid Usage to a Smaller Alphabet

Understanding the origin and early evolution of proteins is important for unveiling how the RNA world developed into an RNA-protein world. Because the composition of organic molecules in the Earth's primitive environment was plausibly not as diverse as today, the number of different amino acids used in early protein synthesis is likely to be substantially less than the current 20 proteinogenic residues. In this study, we have explored the thermal stability and RNA binding of ancestral variants of the ribosomal protein uS8 constructed from a reduced-alphabet of amino acids. First, we built a phylogenetic tree based on the amino acid sequences of uS8 from multiple extant organisms and used the tree to infer two plausible amino acid sequences corresponding to the last bacterial common ancestor of uS8. Both ancestral proteins were thermally stable and bound to an RNA fragment. By eliminating individual amino acid letters and monitoring thermal stability and RNA binding in the resulting proteins, we reduced the size of the amino acid set constituting one of the ancestral proteins, eventually finding that convergent sequences consisting of 15- or 14-amino acid alphabets still folded into stable structures that bound to the RNA fragment. Furthermore, a simplified variant reconstructed from a 13-amino-acid alphabet retained affinity for the RNA fragment, although it lost conformational stability. Collectively, RNA-binding activity may be achieved with a subset of the current 20 amino acids, raising the possibility of a simpler composition of RNA-binding proteins in the earliest stage of protein evolution.

proGenomes3: approaching one million accurately and consistently annotated high-quality prokaryotic genomes

The interpretation of genomic, transcriptomic and other microbial 'omics data is highly dependent on the availability of well-annotated genomes. As the number of publicly available microbial genomes continues to increase exponentially, the need for quality control and consistent annotation is becoming critical. We present proGenomes3, a database of 907 388 high-quality genomes containing 4 billion genes that passed stringent criteria and have been consistently annotated using multiple functional and taxonomic databases including mobile genetic elements and biosynthetic gene clusters. proGenomes3 encompasses 41 171 species-level clusters, defined based on universal single copy marker genes, for which pan-genomes and contextual habitat annotations are provided. The database is available at http://progenomes.embl.de/.

Incompatibility and Interchangeability in Molecular Evolution

There is remarkable variation in the rate at which genetic incompatibilities in molecular interactions accumulate. In some cases, minor changes-even single-nucleotide substitutions-create major incompatibilities when hybridization forces new variants to function in a novel genetic background from an isolated population. In other cases, genes or even entire functional pathways can be horizontally transferred between anciently divergent evolutionary lineages that span the tree of life with little evidence of incompatibilities. In this review, we explore whether there are general principles that can explain why certain genes are prone to incompatibilities while others maintain interchangeability. We summarize evidence pointing to four genetic features that may contribute to greater resistance to functional replacement: (1) function in multisubunit enzyme complexes and protein-protein interactions, (2) sensitivity to changes in gene dosage, (3) rapid rate of sequence evolution, and (4) overall importance to cell viability, which creates sensitivity to small perturbations in molecular function. We discuss the relative levels of support for these different hypotheses and lay out future directions that may help explain the striking contrasts in patterns of incompatibility and interchangeability throughout the history of molecular evolution.

Characterization of a foxtail mosaic virus vector for gene silencing and analysis of innate immune responses in Sorghum bicolor

Sorghum is vulnerable to many biotic and abiotic stresses, which cause considerable yield losses globally. Efforts to genetically characterize beneficial sorghum traits, including disease resistance, plant architecture, and tolerance to abiotic stresses, are ongoing. One challenge faced by sorghum researchers is its recalcitrance to transformation, which has slowed gene validation efforts and utilization for cultivar development. Here, we characterize the use of a foxtail mosaic virus (FoMV) vector for virus-induced gene silencing (VIGS) by targeting two previously tested marker genes: phytoene desaturase (PDS) and ubiquitin (Ub). We additionally demonstrate VIGS of a subgroup of receptor-like cytoplasmic kinases (RLCKs) and report the role of these genes as positive regulators of early defence signalling. Silencing of subgroup 8 RLCKs also resulted in higher susceptibility to the bacterial pathogens Pseudomonas syringae pv. syringae (B728a) and Xanthomonas vasicola pv. holcicola, demonstrating the role of these genes in host defence against bacterial pathogens. Together, this work highlights the utility of FoMV-induced gene silencing in the characterization of genes mediating defence responses in sorghum. Moreover, FoMV was able to systemically infect six diverse sorghum genotypes with high efficiency at optimal temperatures for sorghum growth and therefore could be extrapolated to study additional traits of economic importance.

Characterization, comparative, and functional analysis of arylacetamide deacetylase from Gnathostomata organisms

BACKGROUND

Arylacetamide deacetylase (AADAC) is a lipolytic enzyme involved in xenobiotic metabolism. The characterization in terms of activity and substrate preference has been limited to a few mammalian species. The potential role and catalytic activities of AADAC from other organisms are still poorly understood. Therefore, in this work, the physicochemical properties, proteomic analysis, and protein-protein interactions from Gnathostomata organisms were investigated.

RESULTS

The analysis were performed with 142 orthologue sequences with ~ 48-100% identity with human AADAC. The catalytic motif HGG[A/G] tetrapeptide block was conserved through all AADAC orthologues. Four variations were found in the consensus pentapeptide GXSXG sequence (GDSAG, GESAG, GDSSG, and GSSSG), and a novel motif YXLXP was found. The prediction of N-glycosylation sites projected 4, 1, 6, and 4 different patterns for amphibians, birds, mammals, and reptiles, respectively. The transmembrane regions of AADAC orthologues were not conserved among groups, and variations in the number and orientation of the active site and C-terminal carboxyl were observed among the sequences studied. The protein-protein interaction of AADAC orthologues were related to cancer, lipid, and xenobiotic metabolism genes.

CONCLUSION

The findings from this computational analysis offer new insight into one of the main enzymes involved in xenobiotic metabolism from mammals, reptiles, amphibians, and birds and its potential use in medical and veterinarian biotechnological approaches.

Read2Tree: scalable and accurate phylogenetic trees from raw reads

The inference of phylogenetic trees is foundational to biology. However, state-of-the-art phylogenomics requires running complex pipelines, at significant computational and labour costs, with additional constraints in sequencing coverage, assembly and annotation quality. To overcome these challenges, we present Read2Tree, which directly processes raw sequencing reads into groups of corresponding genes. In a benchmark encompassing a broad variety of datasets, our assembly-free approach was 10-100x faster than conventional approaches, and in most cases more accurate-the exception being when sequencing coverage was high and reference species very distant. To illustrate the broad applicability of the tool, we reconstructed a yeast tree of life of 435 species spanning 590 million years of evolution. Applied to Coronaviridae samples, Read2Tree accurately classified highly diverse animal samples and near-identical SARS-CoV-2 sequences on a single tree-thereby exhibiting remarkable breadth and depth. The speed, accuracy, and versatility of Read2Tree enables comparative genomics at scale.

Detection and characterization of ESBL-producing Escherichia coli and additional co-existence with mcr genes from river water in northern Thailand

Background

Extended-spectrum β-lactamase producing Escherichia coli (ESBL-producing E. coli) have emerged, causing human and animal infections worldwide. This study was conducted to investigate the prevalence and molecular genetic features of ESBL-producing and multidrug-resistant (MDR) E. coli in river water.

Methods

A total of 172 E. coli samples were collected from the Kok River and Kham River in Chiang Rai, Thailand, during a 10-month period (2020-2021).

Results

We detected 45.3% of E. coli to be MDR. The prevalence of ESBL-producers was 22%. Among those ESBL-producing strains, CTX-M-15 (44.7%) was predominantly found, followed by CTX-M-55 (26.3%), CTX-M-14 (18.4%), and CTX-M-27 (10.5%). The bla _TEM-1 and bla _TEM-116 genes were found to be co-harbored with the bla _CTX-M genes. Mobile elements, i.e., ISEcp1 and Tn3, were observed. Twelve plasmid replicons were found, predominantly being IncF (76.3%) and IncFIB (52.6%). Whole genome sequencing of ten selected isolates revealed the co-existence of ESBL with mcr genes in two ESBL-producing E. coli. A wide diversity of MLST classifications was observed. An mcr-1.1-pap2 gene cassette was found to disrupt the PUF2806 domain-containing gene, while an mcr-3.4 contig on another isolate contained the nimC/nimA-mcr-3.4-dgkA core segment.

Discussion

In conclusion, our data provides compelling evidence of MDR and ESBL-producing E. coli, co-existing with mcr genes in river water in northern Thailand, which may be disseminated into other environments and so cause increased risks to public health.