Antibiotic-Induced Gut Microbiota Dysbiosis Modulates Host Transcriptome and m6A Epitranscriptome via Bile Acid Metabolism

Gut microbiota can influence host gene expression and physiology through metabolites. Besides, the presence or absence of gut microbiome can reprogram host transcriptome and epitranscriptome as represented by N6-methyladenosine (m6A), the most abundant mammalian mRNA modification. However, which and how gut microbiota-derived metabolites reprogram host transcriptome and m6A epitranscriptome remain poorly understood. Here, investigation is conducted into how gut microbiota-derived metabolites impact host transcriptome and m6A epitranscriptome using multiple mouse models and multi-omics approaches. Various antibiotics-induced dysbiotic mice are established, followed by fecal microbiota transplantation (FMT) into germ-free mice, and the results show that bile acid metabolism is significantly altered along with the abundance change in bile acid-producing microbiota. Unbalanced gut microbiota and bile acids drastically change the host transcriptome and the m6A epitranscriptome in multiple tissues. Mechanistically, the expression of m6A writer proteins is regulated in animals treated with antibiotics and in cultured cells treated with bile acids, indicating a direct link between bile acid metabolism and m6A biology. Collectively, these results demonstrate that antibiotic-induced gut dysbiosis regulates the landscape of host transcriptome and m6A epitranscriptome via bile acid metabolism pathway. This work provides novel insights into the interplay between microbial metabolites and host gene expression.

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.

Diverse plasmid systems and their ecology across human gut metagenomes revealed by PlasX and MobMess

Plasmids alter microbial evolution and lifestyles by mobilizing genes that often confer fitness in changing environments across clades. Yet our ecological and evolutionary understanding of naturally occurring plasmids is far from complete. Here we developed a machine-learning model, PlasX, which identified 68,350 non-redundant plasmids across human gut metagenomes and organized them into 1,169 evolutionarily cohesive 'plasmid systems' using our sequence containment-aware network-partitioning algorithm, MobMess. Individual plasmids were often country specific, yet most plasmid systems spanned across geographically distinct human populations. Cargo genes in plasmid systems included well-known determinants of fitness, such as antibiotic resistance, but also many others including enzymes involved in the biosynthesis of essential nutrients and modification of transfer RNAs, revealing a wide repertoire of likely fitness determinants in complex environments. Our study introduces computational tools to recognize and organize plasmids, and uncovers the ecological and evolutionary patterns of diverse plasmids in naturally occurring habitats through plasmid systems.

A cryptic plasmid is among the most numerous genetic elements in the human gut

Plasmids are extrachromosomal genetic elements that often encode fitness-enhancing features. However, many bacteria carry "cryptic" plasmids that do not confer clear beneficial functions. We identified one such cryptic plasmid, pBI143, which is ubiquitous across industrialized gut microbiomes and is 14 times as numerous as crAssphage, currently established as the most abundant extrachromosomal genetic element in the human gut. The majority of mutations in pBI143 accumulate in specific positions across thousands of metagenomes, indicating strong purifying selection. pBI143 is monoclonal in most individuals, likely due to the priority effect of the version first acquired, often from one's mother. pBI143 can transfer between Bacteroidales, and although it does not appear to impact bacterial host fitness in vivo, it can transiently acquire additional genetic content. We identified important practical applications of pBI143, including its use in identifying human fecal contamination and its potential as an alternative approach to track human colonic inflammatory states.

Bile acid fitness determinants of a Bacteroides fragilis isolate from a human pouchitis patient

IMPORTANCE

The Gram-negative bacterium Bacteroides fragilis is a common member of the human gut microbiota that colonizes multiple host niches and can influence human physiology through a variety of mechanisms. Identification of genes that enable B. fragilis to grow across a range of host environments has been impeded in part by the relatively limited genetic tractability of this species. We have developed a high-throughput genetic resource for a B. fragilis strain isolated from a UC pouchitis patient. Bile acids limit microbial growth and are altered in abundance in UC pouches, where B. fragilis often blooms. Using this resource, we uncovered pathways and processes that impact B. fragilis fitness in bile and that may contribute to population expansions during bouts of gut inflammation.

Dietary- and host-derived metabolites are used by diverse gut bacteria for anaerobic respiration

Respiratory reductases enable microorganisms to use molecules present in anaerobic ecosystems as energy-generating respiratory electron acceptors. Here we identify three taxonomically distinct families of human gut bacteria (Burkholderiaceae, Eggerthellaceae and Erysipelotrichaceae) that encode large arsenals of tens to hundreds of respiratory-like reductases per genome. Screening species from each family (Sutterella wadsworthensis, Eggerthella lenta and Holdemania filiformis), we discover 22 metabolites used as respiratory electron acceptors in a species-specific manner. Identified reactions transform multiple classes of dietary- and host-derived metabolites, including bioactive molecules resveratrol and itaconate. Products of identified respiratory metabolisms highlight poorly characterized compounds, such as the itaconate-derived 2-methylsuccinate. Reductase substrate profiling defines enzyme-substrate pairs and reveals a complex picture of reductase evolution, providing evidence that reductases with specificities for related cinnamate substrates independently emerged at least four times. These studies thus establish an exceptionally versatile form of anaerobic respiration that directly links microbial energy metabolism to the gut metabolome.

Microbial-enrichment method enables high-throughput metagenomic characterization from host-rich samples

Host-microbe interactions have been linked to health and disease states through the use of microbial taxonomic profiling, mostly via 16S ribosomal RNA gene sequencing. However, many mechanistic insights remain elusive, in part because studying the genomes of microbes associated with mammalian tissue is difficult due to the high ratio of host to microbial DNA in such samples. Here we describe a microbial-enrichment method (MEM), which we demonstrate on a wide range of sample types, including saliva, stool, intestinal scrapings, and intestinal mucosal biopsies. MEM enabled high-throughput characterization of microbial metagenomes from human intestinal biopsies by reducing host DNA more than 1,000-fold with minimal microbial community changes (roughly 90% of taxa had no significant differences between MEM-treated and untreated control groups). Shotgun sequencing of MEM-treated human intestinal biopsies enabled characterization of both high- and low-abundance microbial taxa, pathways and genes longitudinally along the gastrointestinal tract. We report the construction of metagenome-assembled genomes directly from human intestinal biopsies for bacteria and archaea at relative abundances as low as 1%. Analysis of metagenome-assembled genomes reveals distinct subpopulation structures between the small and large intestine for some taxa. MEM opens a path for the microbiome field to acquire deeper insights into host-microbe interactions by enabling in-depth characterization of host-tissue-associated microbial communities.

Bile acid fitness determinants of a Bacteroides fragilis isolate from a human pouchitis patient

Bacteroides fragilis comprises 1-5% of the gut microbiota in healthy humans but can expand to >50% of the population in ulcerative colitis (UC) patients experiencing inflammation. The mechanisms underlying such microbial blooms are poorly understood, but the gut of UC patients has physicochemical features that differ from healthy patients and likely impact microbial physiology. For example, levels of the secondary bile acid deoxycholate (DC) are highly reduced in the ileoanal J-pouch of UC colectomy patients. We isolated a B. fragilis strain from a UC patient with pouch inflammation (i.e. pouchitis) and developed it as a genetic model system to identify genes and pathways that are regulated by DC and that impact B. fragilis fitness in DC and crude bile. Treatment of B. fragilis with a physiologically relevant concentration of DC reduced cell growth and remodeled transcription of one-quarter of the genome. DC strongly induced expression of chaperones and select transcriptional regulators and efflux systems and downregulated protein synthesis genes. Using a barcoded collection of ≈50,000 unique insertional mutants, we further defined B. fragilis genes that contribute to fitness in media containing DC or crude bile. Genes impacting cell envelope functions including cardiolipin synthesis, cell surface glycosylation, and systems implicated in sodium-dependent bioenergetics were major bile acid fitness factors. As expected, there was limited overlap between transcriptionally regulated genes and genes that impacted fitness in bile when disrupted. Our study provides a genome-scale view of a B. fragilis bile response and genetic determinants of its fitness in DC and crude bile.

Species-specific responses of marine bacteria to environmental perturbation

Environmental perturbations shape the structure and function of microbial communities. Oil spills are a major perturbation and resolving spills often requires active measures like dispersant application that can exacerbate the initial disturbance. Species-specific responses of microorganisms to oil and dispersant exposure during such perturbations remain largely unknown. We merged metatranscriptomic libraries with pangenomes to generate Core-Accessory Metatranscriptomes (CA-Metatranscriptomes) for two microbial hydrocarbon degraders that played important roles in the aftermath of the Deepwater Horizon oil spill. The Colwellia CA-Metatranscriptome illustrated pronounced dispersant-driven acceleration of core (~41%) and accessory gene (~59%) transcription, suggesting an opportunistic strategy. Marinobacter responded to oil exposure by expressing mainly accessory genes (~93%), suggesting an effective hydrocarbon-degrading lifestyle. The CA-Metatranscriptome approach offers a robust way to identify the underlying mechanisms of key microbial functions and highlights differences of specialist-vs-opportunistic responses to environmental disturbance.

Microbes with higher metabolic independence are enriched in human gut microbiomes under stress

A wide variety of human diseases are associated with loss of microbial diversity in the human gut, inspiring a great interest in the diagnostic or therapeutic potential of the microbiota. However, the ecological forces that drive diversity reduction in disease states remain unclear, rendering it difficult to ascertain the role of the microbiota in disease emergence or severity. One hypothesis to explain this phenomenon is that microbial diversity is diminished as disease states select for microbial populations that are more fit to survive environmental stress caused by inflammation or other host factors. Here, we tested this hypothesis on a large scale, by developing a software framework to quantify the enrichment of microbial metabolisms in complex metagenomes as a function of microbial diversity. We applied this framework to over 400 gut metagenomes from individuals who are healthy or diagnosed with inflammatory bowel disease (IBD). We found that high metabolic independence (HMI) is a distinguishing characteristic of microbial communities associated with individuals diagnosed with IBD. A classifier we trained using the normalized copy numbers of 33 HMI-associated metabolic modules not only distinguished states of health versus IBD, but also tracked the recovery of the gut microbiome following antibiotic treatment, suggesting that HMI is a hallmark of microbial communities in stressed gut environments.

pWCP is a widely distributed and highly conserved Wolbachia plasmid in Culex pipiens and Culex quinquefasciatus mosquitoes worldwide

Mosquitoes represent the most important pathogen vectors and are responsible for the spread of a wide variety of poorly treatable diseases. Wolbachia are obligate intracellular bacteria that are widely distributed among arthropods and collectively represents one of the most promising solutions for vector control. In particular, Wolbachia has been shown to limit the transmission of pathogens, and to dramatically affect the reproductive behavior of their host through its phage WO. While much research has focused on deciphering and exploring the biocontrol applications of these WO-related phenotypes, the extent and potential impact of the Wolbachia mobilome remain poorly appreciated. Notably, several Wolbachia plasmids, carrying WO-like genes and Insertion Sequences (IS), thus possibly interrelated to other genetic units of the endosymbiont, have been recently discovered. Here we investigated the diversity and biogeography of the first described plasmid of Wolbachia in Culex pipiens (pWCP) in several islands and continental countries around the world-including Cambodia, Guadeloupe, Martinique, Thailand, and Mexico-together with mosquito strains from colonies that evolved for 2 to 30 years in the laboratory. We used PCR and qPCR to determine the presence and copy number of pWCP in individual mosquitoes, and highly accurate Sanger sequencing to evaluate potential variations. Together with earlier observation, our results show that pWCP is omnipresent and strikingly conserved among Wolbachia populations within mosquitoes from distant geographies and environmental conditions. These data suggest a critical role for the plasmid in Wolbachia ecology and evolution, and the potential of a great tool for further genetic dissection and possible manipulation of this endosymbiont.

Metabolic independence drives gut microbial colonization and resilience in health and disease

BACKGROUND

Changes in microbial community composition as a function of human health and disease states have sparked remarkable interest in the human gut microbiome. However, establishing reproducible insights into the determinants of microbial succession in disease has been a formidable challenge.

RESULTS

Here we use fecal microbiota transplantation (FMT) as an in natura experimental model to investigate the association between metabolic independence and resilience in stressed gut environments. Our genome-resolved metagenomics survey suggests that FMT serves as an environmental filter that favors populations with higher metabolic independence, the genomes of which encode complete metabolic modules to synthesize critical metabolites, including amino acids, nucleotides, and vitamins. Interestingly, we observe higher completion of the same biosynthetic pathways in microbes enriched in IBD patients.

CONCLUSIONS

These observations suggest a general mechanism that underlies changes in diversity in perturbed gut environments and reveal taxon-independent markers of "dysbiosis" that may explain why widespread yet typically low-abundance members of healthy gut microbiomes can dominate under inflammatory conditions without any causal association with disease.

A highly conserved and globally prevalent cryptic plasmid is among the most numerous mobile genetic elements in the human gut

Plasmids are extrachromosomal genetic elements that often encode fitness enhancing features. However, many bacteria carry 'cryptic' plasmids that do not confer clear beneficial functions. We identified one such cryptic plasmid, pBI143, which is ubiquitous across industrialized gut microbiomes, and is 14 times as numerous as crAssphage, currently established as the most abundant genetic element in the human gut. The majority of mutations in pBI143 accumulate in specific positions across thousands of metagenomes, indicating strong purifying selection. pBI143 is monoclonal in most individuals, likely due to the priority effect of the version first acquired, often from one's mother. pBI143 can transfer between Bacteroidales and although it does not appear to impact bacterial host fitness in vivo, can transiently acquire additional genetic content. We identified important practical applications of pBI143, including its use in identifying human fecal contamination and its potential as an inexpensive alternative for detecting human colonic inflammatory states.

Transient Suppression of Bacterial Populations Associated with Gut Health is Critical in Success of Exclusive Enteral Nutrition for Children with Crohn's Disease

BACKGROUND AND AIMS

Exclusive enteral nutrition [EEN] is a dietary intervention to induce clinical remission in children with active luminal Crohn's disease [CD]. While changes in the gut microbial communities have been implicated in achieving this remission, a precise understanding of the role of microbial ecology in the restoration of gut homeostasis is lacking.

METHODS

Here we reconstructed genomes from the gut metagenomes of 12 paediatric subjects who were sampled before, during and after EEN. We then classified each microbial population into distinct 'phenotypes' or patterns of response based on changes in their relative abundances throughout the therapy on a per-individual basis.

RESULTS

Our data show that children achieving clinical remission during therapy were enriched with microbial populations that were either suppressed or that demonstrated a transient bloom as a function of EEN. In contrast, this ecosystem-level response was not observed in cases of EEN failure. Further analysis revealed that populations that were suppressed during EEN were significantly more prevalent in healthy children and adults across the globe compared with those that bloomed ephemerally during the therapy.

CONCLUSIONS

These observations taken together suggest that successful outcomes of EEN are marked by a temporary emergence of microbial populations that are rare in healthy individuals, and a concomitant reduction in microbes that are commonly associated with gut homeostasis. Our work is a first attempt to highlight individual-specific, complex environmental factors that influence microbial response in EEN. This model offers a novel, alternative viewpoint to traditional taxonomic strategies used to characterize associations with health and disease states.

Structure-informed microbial population genetics elucidate selective pressures that shape protein evolution

Comprehensive sampling of natural genetic diversity with metagenomics enables highly resolved insights into the interplay between ecology and evolution. However, resolving adaptive, neutral, or purifying processes of evolution from intrapopulation genomic variation remains a challenge, partly due to the sole reliance on gene sequences to interpret variants. Here, we describe an approach to analyze genetic variation in the context of predicted protein structures and apply it to a marine microbial population within the SAR11 subclade 1a.3.V, which dominates low-latitude surface oceans. Our analyses reveal a tight association between genetic variation and protein structure. In a central gene in nitrogen metabolism, we observe decreased occurrence of nonsynonymous variants from ligand-binding sites as a function of nitrate concentrations, revealing genetic targets of distinct evolutionary pressures maintained by nutrient availability. Our work yields insights into the governing principles of evolution and enables structure-aware investigations of microbial population genetics.

Guts of the Urban Ecosystem: Microbial Ecology of Sewer Infrastructure

Microbes have inhabited the oceans and soils for millions of years and are uniquely adapted to their habitat. In contrast, sewer infrastructure in modern cities dates back only ~150 years. Sewer pipes transport human waste and provide a view into public health, but the resident organisms that likely modulate these features are relatively unexplored. Here, we show that the bacterial assemblages sequenced from untreated wastewater in 71 U.S. cities were highly coherent at a fine sequence level, suggesting that urban infrastructure separated by great spatial distances can give rise to strikingly similar communities. Within the overall microbial community structure, temperature had a discernible impact on the distribution patterns of closely related amplicon sequence variants, resulting in warm and cold ecotypes. Two bacterial genera were dominant in most cities regardless of their size or geographic location; on average, Arcobacter accounted for 11% and Acinetobacter 10% of the entire community. Metagenomic analysis of six cities revealed these highly abundant resident organisms carry clinically important antibiotic resistant genes blaCTX-M, blaOXA, and blaTEM. In contrast, human fecal bacteria account for only ~13% of the community; therefore, antibiotic resistance gene inputs from human sources to the sewer system could be comparatively small, which will impact measurement capabilities when monitoring human populations using wastewater. With growing awareness of the metabolic potential of microbes within these vast networks of pipes and the ability to examine the health of human populations, it is timely to increase our understanding of the ecology of these systems. IMPORTANCE Sewer infrastructure is a relatively new habitat comprised of thousands of kilometers of pipes beneath cities. These wastewater conveyance systems contain large reservoirs of microbial biomass with a wide range of metabolic potential and are significant reservoirs of antibiotic resistant organisms; however, we lack an adequate understanding of the ecology or activity of these communities beyond wastewater treatment plants. The striking coherence of the sewer microbiome across the United States demonstrates that the sewer environment is highly selective for a particular microbial community composition. Therefore, results from more in-depth studies or proven engineering controls in one system could be extrapolated more broadly. Understanding the complex ecology of sewer infrastructure is critical for not only improving our ability to treat human waste and increasing the sustainability of our cities but also to create scalable and effective sewage microbial observatories, which are inevitable investments of the future to monitor health in human populations.

Functional repertoire convergence of distantly related eukaryotic plankton lineages abundant in the sunlit ocean

Marine planktonic eukaryotes play critical roles in global biogeochemical cycles and climate. However, their poor representation in culture collections limits our understanding of the evolutionary history and genomic underpinnings of planktonic ecosystems. Here, we used 280 billion Tara Oceans metagenomic reads from polar, temperate, and tropical sunlit oceans to reconstruct and manually curate more than 700 abundant and widespread eukaryotic environmental genomes ranging from 10 Mbp to 1.3 Gbp. This genomic resource covers a wide range of poorly characterized eukaryotic lineages that complement long-standing contributions from culture collections while better representing plankton in the upper layer of the oceans. We performed the first, to our knowledge, comprehensive genome-wide functional classification of abundant unicellular eukaryotic plankton, revealing four major groups connecting distantly related lineages. Neither trophic modes of plankton nor its vertical evolutionary history could completely explain the functional repertoire convergence of major eukaryotic lineages that coexisted within oceanic currents for millions of years.

Microbial metabolites in the marine carbon cycle

One-quarter of photosynthesis-derived carbon on Earth rapidly cycles through a set of short-lived seawater metabolites that are generated from the activities of marine phytoplankton, bacteria, grazers and viruses. Here we discuss the sources of microbial metabolites in the surface ocean, their roles in ecology and biogeochemistry, and approaches that can be used to analyse them from chemistry, biology, modelling and data science. Although microbial-derived metabolites account for only a minor fraction of the total reservoir of marine dissolved organic carbon, their flux and fate underpins the central role of the ocean in sustaining life on Earth.

Unifying the known and unknown microbial coding sequence space

Genes of unknown function are among the biggest challenges in molecular biology, especially in microbial systems, where 40–60% of the predicted genes are unknown. Despite previous attempts, systematic approaches to include the unknown fraction into analytical workflows are still lacking. Here, we present a conceptual framework, its translation into the computational workflow AGNOSTOS and a demonstration on how we can bridge the known-unknown gap in genomes and metagenomes. By analyzing 415,971,742 genes predicted from 1749 metagenomes and 28,941 bacterial and archaeal genomes, we quantify the extent of the unknown fraction, its diversity, and its relevance across multiple organisms and environments. The unknown sequence space is exceptionally diverse, phylogenetically more conserved than the known fraction and predominantly taxonomically restricted at the species level. From the 71 M genes identified to be of unknown function, we compiled a collection of 283,874 lineage-specific genes of unknown function for Cand. Patescibacteria (also known as Candidate Phyla Radiation, CPR), which provides a significant resource to expand our understanding of their unusual biology. Finally, by identifying a target gene of unknown function for antibiotic resistance, we demonstrate how we can enable the generation of hypotheses that can be used to augment experimental data.

It is estimated that scientists do not know what half of microbial genes actually do. When these genes are discovered in microorganisms grown in the lab or found in environmental samples, it is not possible to identify what their roles are. Many of these genes are excluded from further analyses for these reasons, meaning that the study of microbial genes tends to be limited to genes that have already been described.

These limitations hinder research into microbiology, because information from newly discovered genes cannot be integrated to better understand how these organisms work. Experiments to understand what role these genes have in the microorganisms are labor-intensive, so new analytical strategies are needed.

To do this, Vanni et al. developed a new framework to categorize genes with unknown roles, and a computational workflow to integrate them into traditional analyses. When this approach was applied to over 400 million microbial genes (both with known and unknown roles), it showed that the share of genes with unknown functions is only about 30 per cent, smaller than previously thought. The analysis also showed that these genes are very diverse, revealing a huge space for future research and potential applications. Combining their approach with experimental data, Vanni et al. were able to identify a gene with a previously unknown purpose that could be involved in antibiotic resistance.

This system could be useful for other scientists studying microorganisms to get a more complete view of microbial systems. In future, it may also be used to analyze the genetics of other organisms, such as plants and animals.

tRNA modification dynamics from individual organisms to metaepitranscriptomics of microbiomes

tRNA is the most extensively modified RNA in cells. On average, a bacterial tRNA contains 8 modifications per molecule and a eukaryotic tRNA contains 13 modifications per molecule. Recent studies reveal that tRNA modifications are highly dynamic and respond extensively to environmental conditions. Functions of tRNA modification dynamics include enhanced, on-demand decoding of specific codons in response genes and regulation of tRNA fragment biogenesis. This review summarizes recent advances in the studies of tRNA modification dynamics in biological processes, tRNA modification erasers, and human-associated bacteria. Furthermore, we use the term "metaepitranscriptomics" to describe the potential and approach of tRNA modification studies in natural biological communities such as microbiomes. tRNA is highly modified in cells, and tRNA modifications respond extensively to environmental conditions to enhance translation of specific genes and produce tRNA fragments on demand. We review recent advances in tRNA sequencing methods, tRNA modification dynamics in biological processes, and tRNA modification studies in natural communities such as the microbiomes.

Eco-evolutionary significance of domesticated retroelements in microbial genomes

Since the first discovery of reverse transcriptase in bacteria, and later in archaea, bacterial and archaeal retroelements have been defined by their common enzyme that coordinates diverse functions. Yet, evolutionary refinement has produced distinct retroelements across the tree of microbial life that are perhaps best described in terms of their programmed RNA—a compact sequence that preserves core information for a sophisticated mechanism. From this perspective, reverse transcriptase has been selected as the modular tool for carrying out nature’s instructions in various RNA templates. Beneficial retroelements—those that can provide a fitness advantage to their host—evolved to their extant forms in a wide array of microorganisms and their viruses, spanning nearly all habitats. Within each specialized retroelement class, several universal features seem to be shared across diverse taxa, while specific functional and mechanistic insights are based on only a few model retroelement systems from clinical isolates. Currently, little is known about the diversity of cellular functions and ecological significance of retroelements across different biomes. With increasing availability of isolate, metagenome-assembled, and single-amplified genomes, the taxonomic and functional breadth of prokaryotic retroelements is coming into clearer view. This review explores the recently characterized classes of beneficial, yet accessory retroelements of bacteria and archaea. We describe how these specialized mechanisms exploit a form of fixed mobility, whereby the retroelements do not appear to proliferate selfishly throughout the genome. Moreover, we discuss computational approaches for systematic identification of retroelements from vast sequence repositories and highlight recent discoveries in terms of their apparent distribution and ecological significance in nature. Lastly, we present a new perspective on the eco-evolutionary significance of these genetic elements in marine bacteria and demonstrate approaches that enable the characterization of their environmental diversity through metagenomics.

High molecular weight DNA extraction strategies for long-read sequencing of complex metagenomes

By offering extremely long contiguous characterization of individual DNA molecules, rapidly emerging long‐read sequencing strategies offer comprehensive insights into the organization of genetic information in genomes and metagenomes. However, successful long‐read sequencing experiments demand high concentrations of highly purified DNA of high molecular weight (HMW), which limits the utility of established DNA extraction kits designed for short‐read sequencing. The challenges associated with input DNA quality intensify further when working with complex environmental samples of low microbial biomass, which requires new protocols that are tailored to study metagenomes with long‐read sequencing. Here, we use human tongue scrapings to benchmark six HMW DNA extraction strategies that are based on commercially available kits, phenol–chloroform (PC) extraction and agarose encasement followed by agarase digestion. A typical end goal of HMW DNA extractions is to obtain the longest possible reads during sequencing, which is often achieved by PC extractions, as demonstrated in sequencing of cultured cells. Yet our analyses that consider overall read‐size distribution, assembly performance and the number of circularized elements found in sequencing results suggest that column‐based kits with enzyme supplementation, rather than PC methods, may be more appropriate for long‐read sequencing of metagenomes.

The mosquito microbiome includes habitat-specific but rare symbionts

Microbial communities are known to influence mosquito lifestyles by modifying essential metabolic and behavioral processes that affect reproduction, development, immunity, digestion, egg survival, and the ability to transmit pathogens. Many studies have used 16S rRNA gene amplicons to characterize mosquito microbiota and investigate factors that influence host-microbiota dynamics. However, a relatively low taxonomic resolution due to clustering methods based on arbitrary threshold and the overall dominance of Wolbachia or Asaia symbionts obscured the investigation of rare members of mosquito microbiota in previous studies. Here, we used high resolution Shannon entropy-based oligotyping approach to analyze the microbiota of Culex pipiens, Culex quinquefasciatus and Aedes individuals from continental Southern France and overseas Guadeloupe as well as from laboratories with or without antibiotics treatment. Our experimental design that resulted in a series of mosquito samples with a gradient of Wolbachia density and relative abundance along with high-resolution analyses of amplicon sequences enabled the recovery of a robust signal from typically less accessible bacterial taxa. Our data confirm species-specific mosquito-bacteria associations with geography as a primary factor that influences bacterial community structure. But interestingly, they also reveal co-occurring symbiotic bacterial variants within single individuals for both Elizabethkingia and Erwinia genera, distinct and specific Asaia and Chryseobacterium in continental and overseas territories, and a putative rare Wolbachia variant. Overall, our study reveals the presence of previously overlooked microdiversity and multiple closely related symbiotic strains within mosquito individuals with a remarkable habitat-specificity.

Interferon inducible pseudouridine modification in human mRNA by quantitative nanopore profiling

Pseudouridine (Ψ) is an abundant mRNA modification in mammalian transcriptome, but its functions have remained elusive due to the difficulty of transcriptome-wide mapping. We develop a nanopore native RNA sequencing method for quantitative Ψ prediction (NanoPsu) that utilizes native content training, machine learning modeling, and single-read linkage analysis. Biologically, we find interferon inducible Ψ modifications in interferon-stimulated gene transcripts which are consistent with a role of Ψ in enabling efficacy of mRNA vaccines.

Heterotrophic bacterial diazotrophs are more abundant than their cyanobacterial counterparts in metagenomes covering most of the sunlit ocean

Biological nitrogen fixation contributes significantly to marine primary productivity. The current view depicts few cyanobacterial diazotrophs as the main marine nitrogen fixers. Here, we used 891 Tara Oceans metagenomes derived from surface waters of five oceans and two seas to generate a manually curated genomic database corresponding to free-living, filamentous, colony-forming, particle-attached, and symbiotic bacterial and archaeal populations. The database provides the genomic content of eight cyanobacterial diazotrophs including a newly discovered population related to known heterocystous symbionts of diatoms, as well as 40 heterotrophic bacterial diazotrophs that considerably expand the known diversity of abundant marine nitrogen fixers. These 48 populations encapsulate 92% of metagenomic signal for known nifH genes in the sunlit ocean, suggesting that the genomic characterization of the most abundant marine diazotrophs may be nearing completion. Newly identified heterotrophic bacterial diazotrophs are widespread, express their nifH genes in situ, and also occur in large planktonic size fractions where they might form aggregates that provide the low-oxygen microenvironments required for nitrogen fixation. Critically, we found heterotrophic bacterial diazotrophs to be more abundant than cyanobacterial diazotrophs in most metagenomes from the open oceans and seas, emphasizing the importance of a wide range of heterotrophic populations in the marine nitrogen balance.

Early-Life Microbial Restitution Reduces Colitis Risk Promoted by Antibiotic-Induced Gut Dysbiosis in Interleukin 10-/- Mice

BACKGROUND & AIMS

Perturbations in the early-life gut microbiome are associated with increased risk for complex immune disorders like inflammatory bowel diseases. We previously showed that maternal antibiotic-induced gut dysbiosis vertically transmitted to offspring increases experimental colitis risk in interleukin (IL) 10 gene deficient (IL10-/-) mice, a finding that may result from the loss/lack of essential microbes needed for appropriate immunologic education early in life. Here, we aimed to identify key microbes required for proper development of the early-life gut microbiome that decrease colitis risk in genetically susceptible animals.

METHODS

Metagenomic sequencing followed by reconstruction of metagenome-assembled genomes was performed on fecal samples of IL10-/- mice with and without antibiotic-induced dysbiosis to identify potential missing microbial members needed for immunologic education. One high-value target strain was then engrafted early and/or late into the gut microbiomes of IL10-/- mice with antibiotic-induced dysbiosis.

RESULTS

Early-, but not late-, life engraftment of a single dominant Bacteroides strain of non-antibiotic-treated IL10-/- mice was sufficient to restore the development of the gut microbiome, promote immune tolerance, and prevent colitis in IL10-/- mice that had antibiotic-induced dysbiosis.

CONCLUSIONS

Restitution of a keystone microbial strain missing in the early-life antibiotic-induced gut dysbiosis results in recovery of the microbiome, proper development of immune tolerance, and reduced risk for colitis in genetically prone hosts.

STRONG: metagenomics strain resolution on assembly graphs

We introduce STrain Resolution ON assembly Graphs (STRONG), which identifies strains de novo, from multiple metagenome samples. STRONG performs coassembly, and binning into metagenome assembled genomes (MAGs), and stores the coassembly graph prior to variant simplification. This enables the subgraphs and their unitig per-sample coverages, for individual single-copy core genes (SCGs) in each MAG, to be extracted. A Bayesian algorithm, BayesPaths, determines the number of strains present, their haplotypes or sequences on the SCGs, and abundances. STRONG is validated using synthetic communities and for a real anaerobic digestor time series generates haplotypes that match those observed from long Nanopore reads.

Differential prevalence and host-association of antimicrobial resistance traits in disinfected and non-disinfected drinking water systems

Antimicrobial resistance (AMR) in drinking water has received less attention than its counterparts in the urban water cycle. While culture-based techniques or gene-centric PCR have been used to probe the impact of treatment approaches (e.g., disinfection) on AMR in drinking water, to our knowledge there is no systematic comparison of AMR trait distribution and prevalence between disinfected and disinfectant residual-free drinking water systems. We used metagenomics to assess the associations between disinfectant residuals and AMR prevalence and its host association in full-scale drinking water distribution systems (DWDSs) with and without disinfectant residuals. While the differences in AMR profiles between DWDSs were associated with the presence or absence of disinfectant, they were also associated with overall water chemistry and more importantly with microbial community structure. AMR genes and mechanisms differentially abundant in disinfected systems were primarily associated with nontuberculous mycobacteria (NTM). Finally, evaluation of metagenome assembled genomes (MAGs) also suggests that NTM possessing AMR genes conferring intrinsic resistance to key antibiotics were prevalent in disinfected systems, whereas such NTM genomes were not detected in disinfectant residual free DWDSs. Altogether, our findings provide insights into the drinking water resistome and its association with potential opportunistic pathogens, particularly in systems with disinfectant residual.

Functional and genetic markers of niche partitioning among enigmatic members of the human oral microbiome

INTRODUCTION

Microbial residents of the human oral cavity have long been a major focus of microbiology due to their influence on host health and intriguing patterns of site specificity amidst the lack of dispersal limitation. However, the determinants of niche partitioning in this habitat are yet to be fully understood, especially among taxa that belong to recently discovered branches of microbial life.

RESULTS

Here, we assemble metagenomes from tongue and dental plaque samples from multiple individuals and reconstruct 790 non-redundant genomes, 43 of which resolve to TM7, a member of the Candidate Phyla Radiation, forming six monophyletic clades that distinctly associate with either plaque or tongue. Both pangenomic and phylogenomic analyses group tongue-specific clades with other host-associated TM7 genomes. In contrast, plaque-specific TM7 group with environmental TM7 genomes. Besides offering deeper insights into the ecology, evolution, and mobilome of cryptic members of the oral microbiome, our study reveals an intriguing resemblance between dental plaque and non-host environments indicated by the TM7 evolution, suggesting that plaque may have served as a stepping stone for environmental microbes to adapt to host environments for some clades of microbes. Additionally, we report that prophages are widespread among oral-associated TM7, while absent from environmental TM7, suggesting that prophages may have played a role in adaptation of TM7 to the host environment.

CONCLUSIONS

Our data illuminate niche partitioning of enigmatic members of the oral cavity, including TM7, SR1, and GN02, and provide genomes for poorly characterized yet prevalent members of this biome, such as uncultivated Flavobacteriaceae.

Metapangenomics of the oral microbiome provides insights into habitat adaptation and cultivar diversity

BACKGROUND

The increasing availability of microbial genomes and environmental shotgun metagenomes provides unprecedented access to the genomic differences within related bacteria. The human oral microbiome with its diverse habitats and abundant, relatively well-characterized microbial inhabitants presents an opportunity to investigate bacterial population structures at an ecosystem scale.

RESULTS

Here, we employ a metapangenomic approach that combines public genomes with Human Microbiome Project (HMP) metagenomes to study the diversity of microbial residents of three oral habitats: tongue dorsum, buccal mucosa, and supragingival plaque. For two exemplar taxa, Haemophilus parainfluenzae and the genus Rothia, metapangenomes reveal distinct genomic groups based on shared genome content. H. parainfluenzae genomes separate into three distinct subgroups with differential abundance between oral habitats. Functional enrichment analyses identify an operon encoding oxaloacetate decarboxylase as diagnostic for the tongue-abundant subgroup. For the genus Rothia, grouping by shared genome content recapitulates species-level taxonomy and habitat preferences. However, while most R. mucilaginosa are restricted to the tongue as expected, two genomes represent a cryptic population of R. mucilaginosa in many buccal mucosa samples. For both H. parainfluenzae and the genus Rothia, we identify not only limitations in the ability of cultivated organisms to represent populations in their native environment, but also specifically which cultivar gene sequences are absent or ubiquitous.

CONCLUSIONS

Our findings provide insights into population structure and biogeography in the mouth and form specific hypotheses about habitat adaptation. These results illustrate the power of combining metagenomes and pangenomes to investigate the ecology and evolution of bacteria across analytical scales.

B cell superantigens in the human intestinal microbiota

IgA is prominently secreted at mucosal surfaces and coats a fraction of the commensal microbiota, a process that is critical for intestinal homeostasis. However, the mechanisms of IgA induction and the molecular targets of these antibodies remain poorly understood, particularly in humans. Here, we demonstrate that microbiota from a subset of human individuals encode two protein "superantigens" expressed on the surface of commensal bacteria of the family Lachnospiraceae such as Ruminococcus gnavus that bind IgA variable regions and stimulate potent IgA responses in mice. These superantigens stimulate B cells expressing human VH3 or murine VH5/6/7 variable regions and subsequently bind their antibodies, allowing these microbial organisms to become highly coated with IgA in vivo. These findings demonstrate a previously unappreciated role for commensal superantigens in host-microbiota interactions. Furthermore, as superantigen-expressing strains show an uneven distribution across human populations, they should be systematically considered in studies evaluating human B cell responses and microbiota during homeostasis and disease.

Droplet-based high-throughput cultivation for accurate screening of antibiotic resistant gut microbes

Traditional cultivation approaches in microbiology are labor-intensive, low-throughput, and yield biased sampling of environmental microbes due to ecological and evolutionary factors. New strategies are needed for ample representation of rare taxa and slow-growers that are often outcompeted by fast-growers in cultivation experiments. Here we describe a microfluidic platform that anaerobically isolates and cultivates microbial cells in millions of picoliter droplets and automatically sorts them based on colony density to enhance slow-growing organisms. We applied our strategy to a fecal microbiota transplant (FMT) donor stool using multiple growth media, and found significant increase in taxonomic richness and larger representation of rare and clinically relevant taxa among droplet-grown cells compared to conventional plates. Furthermore, screening the FMT donor stool for antibiotic resistance revealed 21 populations that evaded detection in plate-based assessment of antibiotic resistance. Our method improves cultivation-based surveys of diverse microbiomes to gain deeper insights into microbial functioning and lifestyles.

Disinfection exhibits systematic impacts on the drinking water microbiome

Limiting microbial growth during drinking water distribution is achieved either by maintaining a disinfectant residual or through nutrient limitation without using a disinfectant. The impact of these contrasting approaches on the drinking water microbiome is not systematically understood. We use genome-resolved metagenomics to compare the structure, metabolic traits, and population genomes of drinking water microbiome samples from bulk drinking water across multiple full-scale disinfected and non-disinfected drinking water systems. Microbial communities cluster at the structural- and functional potential-level based on the presence/absence of a disinfectant residual. Disinfectant residual alone explained 17 and 6.5% of the variance in structure and functional potential of the drinking water microbiome, respectively, despite including multiple drinking water systems with variable source waters and source water communities and treatment strategies. The drinking water microbiome is structurally and functionally less diverse and variable across disinfected compared to non-disinfected systems. While bacteria were the most abundant domain, archaea and eukaryota were more abundant in non-disinfected and disinfected systems, respectively. Community-level differences in functional potential were driven by enrichment of genes associated with carbon and nitrogen fixation in non-disinfected systems and γ-aminobutyrate metabolism in disinfected systems likely associated with the recycling of amino acids. Genome-level analyses for a subset of phylogenetically-related microorganisms suggests that disinfection selects for microorganisms capable of using fatty acids, presumably from microbial decay products, via the glyoxylate cycle. Overall, we find that disinfection exhibits systematic selective pressures on the drinking water microbiome and may select for microorganisms able to utilize microbial decay products originating from disinfection-inactivated microorganisms. Video abstract.

Accurate and complete genomes from metagenomes

Genomes are an integral component of the biological information about an organism; thus, the more complete the genome, the more informative it is. Historically, bacterial and archaeal genomes were reconstructed from pure (monoclonal) cultures, and the first reported sequences were manually curated to completion. However, the bottleneck imposed by the requirement for isolates precluded genomic insights for the vast majority of microbial life. Shotgun sequencing of microbial communities, referred to initially as community genomics and subsequently as genome-resolved metagenomics, can circumvent this limitation by obtaining metagenome-assembled genomes (MAGs); but gaps, local assembly errors, chimeras, and contamination by fragments from other genomes limit the value of these genomes. Here, we discuss genome curation to improve and, in some cases, achieve complete (circularized, no gaps) MAGs (CMAGs). To date, few CMAGs have been generated, although notably some are from very complex systems such as soil and sediment. Through analysis of about 7000 published complete bacterial isolate genomes, we verify the value of cumulative GC skew in combination with other metrics to establish bacterial genome sequence accuracy. The analysis of cumulative GC skew identified potential misassemblies in some reference genomes of isolated bacteria and the repeat sequences that likely gave rise to them. We discuss methods that could be implemented in bioinformatic approaches for curation to ensure that metabolic and evolutionary analyses can be based on very high-quality genomes.

Global phylogeography and ancient evolution of the widespread human gut virus crAssphage

Microbiomes are vast communities of microorganisms and viruses that populate all natural ecosystems. Viruses have been considered to be the most variable component of microbiomes, as supported by virome surveys and examples of high genomic mosaicism. However, recent evidence suggests that the human gut virome is remarkably stable compared with that of other environments. Here, we investigate the origin, evolution and epidemiology of crAssphage, a widespread human gut virus. Through a global collaboration, we obtained DNA sequences of crAssphage from more than one-third of the world's countries and showed that the phylogeography of crAssphage is locally clustered within countries, cities and individuals. We also found fully colinear crAssphage-like genomes in both Old-World and New-World primates, suggesting that the association of crAssphage with primates may be millions of years old. Finally, by exploiting a large cohort of more than 1,000 individuals, we tested whether crAssphage is associated with bacterial taxonomic groups of the gut microbiome, diverse human health parameters and a wide range of dietary factors. We identified strong correlations with different clades of bacteria that are related to Bacteroidetes and weak associations with several diet categories, but no significant association with health or disease. We conclude that crAssphage is a benign cosmopolitan virus that may have coevolved with the human lineage and is an integral part of the normal human gut virome.

Global phylogeography and ancient evolution of the widespread human gut virus crAssphage

Microbiomes are vast communities of microorganisms and viruses that populate all natural ecosystems. Viruses have been considered to be the most variable component of microbiomes, as supported by virome surveys and examples of high genomic mosaicism. However, recent evidence suggests that the human gut virome is remarkably stable compared with that of other environments. Here, we investigate the origin, evolution and epidemiology of crAssphage, a widespread human gut virus. Through a global collaboration, we obtained DNA sequences of crAssphage from more than one-third of the world's countries and showed that the phylogeography of crAssphage is locally clustered within countries, cities and individuals. We also found fully colinear crAssphage-like genomes in both Old-World and New-World primates, suggesting that the association of crAssphage with primates may be millions of years old. Finally, by exploiting a large cohort of more than 1,000 individuals, we tested whether crAssphage is associated with bacterial taxonomic groups of the gut microbiome, diverse human health parameters and a wide range of dietary factors. We identified strong correlations with different clades of bacteria that are related to Bacteroidetes and weak associations with several diet categories, but no significant association with health or disease. We conclude that crAssphage is a benign cosmopolitan virus that may have coevolved with the human lineage and is an integral part of the normal human gut virome.

Single-amino acid variants reveal evolutionary processes that shape the biogeography of a global SAR11 subclade

Members of the SAR11 order Pelagibacterales dominate the surface oceans. Their extensive diversity challenges emerging operational boundaries defined for microbial 'species' and complicates efforts of population genetics to study their evolution. Here, we employed single-amino acid variants (SAAVs) to investigate ecological and evolutionary forces that maintain the genomic heterogeneity within ubiquitous SAR11 populations we accessed through metagenomic read recruitment using a single isolate genome. Integrating amino acid and protein biochemistry with metagenomics revealed that systematic purifying selection against deleterious variants governs non-synonymous variation among very closely related populations of SAR11. SAAVs partitioned metagenomes into two main groups matching large-scale oceanic current temperatures, and six finer proteotypes that connect distant oceanic regions. These findings suggest that environmentally-mediated selection plays a critical role in the journey of cosmopolitan surface ocean microbial populations, and the idea 'everything is everywhere but the environment selects' has credence even at the finest resolutions.

Genome-resolved insights into a novel Spiroplasma symbiont of the Wheat Stem Sawfly (Cephus cinctus)

Arthropods often have obligate relationships with symbiotic microbes, and recent investigations have demonstrated that such host-microbe relationships could be exploited to suppress natural populations of vector carrying mosquitos. Strategies that target the interplay between agricultural pests and their symbionts could decrease the burden caused by agricultural pests; however, the lack of comprehensive genomic insights into naturally occurring microbial symbionts presents a significant bottleneck. Here we employed amplicon surveys, genome-resolved metagenomics, and scanning electron microscopy to investigate symbionts of the wheat stem sawfly (Cephus cinctus), a major pest that causes an estimated $350 million dollars or more in wheat yield losses in the northwestern United States annually. Through 16S rRNA gene sequencing of two major haplotypes and life stages of wheat stem sawfly, we show a novel Spiroplasma species is ever-present and predominant, with phylogenomic analyses placing it as a member of the ixodetis clade of mollicutes. Using state-of-the-art metagenomic assembly and binning strategies we were able to reconstruct a 714 Kb, 72.7%-complete Spiroplasma genome, which represents just the second draft genome from the ixodetis clade of mollicutes. Functional annotation of the Spiroplasma genome indicated carbohydrate-metabolism involved PTS-mediated import of glucose and fructose followed by glycolysis to lactate, acetate, and propionoate. The bacterium also encoded biosynthetic pathways for essential vitamins B2, B3, and B9. We identified putative Spiroplasma virulence genes: cardiolipin and chitinase. These results identify a previously undescribed symbiosis between wheat stem sawfly and a novel Spiroplasma sp., availing insight into their molecular relationship, and may yield new opportunities for microbially-mediated pest control strategies.

Biochemical diversity of glycosphingolipid biosynthesis as a driver of Coccolithovirus competitive ecology

Coccolithoviruses (EhVs) are large, double-stranded DNA-containing viruses that infect the single-celled, marine coccolithophore Emiliania huxleyi. Given the cosmopolitan nature and global importance of E. huxleyi as a bloom-forming, calcifying, photoautotroph, E. huxleyi-EhV interactions play a key role in oceanic carbon biogeochemistry. Virally-encoded glycosphingolipids (vGSLs) are virulence factors that are produced by the activity of virus-encoded serine palmitoyltransferase (SPT). Here, we characterize the dynamics, diversity and catalytic production of vGSLs in an array of EhV strains in relation to their SPT sequence composition and explore the hypothesis that they are a determinant of infectivity and host demise. vGSL production and diversity was positively correlated with increased virulence, virus replication rate and lytic infection dynamics in laboratory experiments, but they do not explain the success of less-virulent EhVs in natural EhV communities. The majority of EhV-derived SPT amplicon sequences associated with infected cells in the North Atlantic derived from slower infecting, less virulent EhVs. Our lab-, field- and mathematical model-based data and simulations support ecological scenarios whereby slow-infecting, less-virulent EhVs successfully compete in North Atlantic populations of E. huxleyi, through either the preferential removal of fast-infecting, virulent EhVs during active infection or by having access to a broader host range.

The Wolbachia mobilome in Culex pipiens includes a putative plasmid

Wolbachia is a genus of obligate intracellular bacteria found in nematodes and arthropods worldwide, including insect vectors that transmit dengue, West Nile, and Zika viruses. Wolbachia's unique ability to alter host reproductive behavior through its temperate bacteriophage WO has enabled the development of new vector control strategies. However, our understanding of Wolbachia's mobilome beyond its bacteriophages is incomplete. Here, we reconstruct near-complete Wolbachia genomes from individual ovary metagenomes of four wild Culex pipiens mosquitoes captured in France. In addition to viral genes missing from the Wolbachia reference genome, we identify a putative plasmid (pWCP), consisting of a 9.23-kbp circular element with 14 genes. We validate its presence in additional Culex pipiens mosquitoes using PCR, long-read sequencing, and screening of existing metagenomes. The discovery of this previously unrecognized extrachromosomal element opens additional possibilities for genetic manipulation of Wolbachia.

Co-occurring genomic capacity for anaerobic methane and dissimilatory sulfur metabolisms discovered in the Korarchaeota

Phylogenetic and geological evidence supports the hypothesis that life on Earth originated in thermal environments and conserved energy through methanogenesis or sulfur reduction. Here we describe two populations of the deeply rooted archaeal phylum Korarchaeota, which were retrieved from the metagenome of a circumneutral, suboxic hot spring that contains high levels of sulfate, sulfide, methane, hydrogen and carbon dioxide. One population is closely related to 'Candidatus Korarchaeum cryptofilum OPF8', while the more abundant korarchaeote, 'Candidatus Methanodesulfokores washburnensis', contains genes that are necessary for anaerobic methane and dissimilatory sulfur metabolisms. Phylogenetic and ancestral reconstruction analyses suggest that methane metabolism originated in the Korarchaeota, whereas genes for dissimilatory sulfite reduction were horizontally transferred to the Korarchaeota from the Firmicutes. Interactions among enzymes involved in both metabolisms could facilitate exergonic, sulfite-dependent, anaerobic oxidation of methane to methanol; alternatively, 'Ca. M. washburnensis' could conduct methanogenesis and sulfur reduction independently. Metabolic reconstruction suggests that 'Ca. M. washburnensis' is a mixotroph, capable of amino acid uptake, assimilation of methane-derived carbon and/or CO₂ fixation by archaeal type III-b RuBisCO for scavenging ribose carbon. Our findings link anaerobic methane metabolism and dissimilatory sulfur reduction within a single deeply rooted archaeal population and have implications for the evolution of these traits throughout the Archaea.

Microbiome characterization by high-throughput transfer RNA sequencing and modification analysis

Advances in high-throughput sequencing have facilitated remarkable insights into the diversity and functioning of naturally occurring microbes; however, current sequencing strategies are insufficient to reveal physiological states of microbial communities associated with protein translation dynamics. Transfer RNAs (tRNAs) are core components of protein synthesis machinery, present in all living cells, and are phylogenetically tractable, which make them ideal targets to gain physiological insights into environmental microbes. Here we report a direct sequencing approach, tRNA-seq, and a software suite, tRNA-seq-tools, to recover sequences, abundance profiles, and post-transcriptional modifications of microbial tRNA transcripts. Our analysis of cecal samples using tRNA-seq distinguishes high-fat- and low-fat-fed mice in a comparable fashion to 16S ribosomal RNA gene amplicons, and reveals taxon- and diet-dependent variations in tRNA modifications. Our results provide taxon-specific in situ insights into the dynamics of tRNA gene expression and post-transcriptional modifications within complex environmental microbiomes.

Fecal source identification using random forest

BACKGROUND

Clostridiales and Bacteroidales are uniquely adapted to the gut environment and have co-evolved with their hosts resulting in convergent microbiome patterns within mammalian species. As a result, members of Clostridiales and Bacteroidales are particularly suitable for identifying sources of fecal contamination in environmental samples. However, a comprehensive evaluation of their predictive power and development of computational approaches is lacking. Given the global public health concern for waterborne disease, accurate identification of fecal pollution sources is essential for effective risk assessment and management. Here, we use random forest algorithm and 16S rRNA gene amplicon sequences assigned to Clostridiales and Bacteroidales to identify common fecal pollution sources. We benchmarked the accuracy, consistency, and sensitivity of our classification approach using fecal, environmental, and artificial in silico generated samples.

RESULTS

Clostridiales and Bacteroidales classifiers were composed mainly of sequences that displayed differential distributions (host-preferred) among sewage, cow, deer, pig, cat, and dog sources. Each classifier correctly identified human and individual animal sources in approximately 90% of the fecal and environmental samples tested. Misclassifications resulted mostly from false-positive identification of cat and dog fecal signatures in host animals not used to build the classifiers, suggesting characterization of additional animals would improve accuracy. Random forest predictions were highly reproducible, reflecting the consistency of the bacterial signatures within each of the animal and sewage sources. Using in silico generated samples, we could detect fecal bacterial signatures when the source dataset accounted for as little as ~ 0.5% of the assemblage, with ~ 0.04% of the sequences matching the classifiers. Finally, we developed a proxy to estimate proportions among sources, which allowed us to determine which sources contribute the most to observed fecal pollution.

CONCLUSION

Random forest classification with 16S rRNA gene amplicons offers a rapid, sensitive, and accurate solution for identifying host microbial signatures to detect human and animal fecal contamination in environmental samples.

New Biological Insights Into How Deforestation in Amazonia Affects Soil Microbial Communities Using Metagenomics and Metagenome-Assembled Genomes

Deforestation in the Brazilian Amazon occurs at an alarming rate, which has broad effects on global greenhouse gas emissions, carbon storage, and biogeochemical cycles. In this study, soil metagenomes and metagenome-assembled genomes (MAGs) were analyzed for alterations to microbial community composition, functional groups, and putative physiology as it related to land-use change and tropical soil. A total of 28 MAGs were assembled encompassing 10 phyla, including both dominant and rare biosphere lineages. Amazon Acidobacteria subdivision 3, Melainabacteria, Microgenomates, and Parcubacteria were found exclusively in pasture soil samples, while Candidatus Rokubacteria was predominant in the adjacent rainforest soil. These shifts in relative abundance between land-use types were supported by the different putative physiologies and life strategies employed by the taxa. This research provides unique biological insights into candidate phyla in tropical soil and how deforestation may impact the carbon cycle and affect climate change.

Minimizing confounders and increasing data quality in murine models for studies of the gut microbiome

Murine models are widely used to explore host-microbe interactions because of the challenges and limitations inherent to human studies. However, microbiome studies in murine models are not without their nuances. Inter-individual variations in gut microbiota are frequent even in animals housed within the same room. We therefore sought to find an efficient and effective standard operating procedure (SOP) to minimize these effects to improve consistency and reproducibility in murine microbiota studies. Mice were housed in a single room under specific-pathogen free conditions. Soiled cage bedding was routinely mixed weekly and distributed among all cages from weaning (three weeks old) until the onset of the study. Females and males were separated by sex and group-housed (up to five mice/cage) at weaning. 16S rRNA gene analyses of fecal samples showed that this protocol significantly reduced pre-study variability of gut microbiota amongst animals compared to other conventional measures used to normalize microbiota when large experimental cohorts have been required. A significant and consistent effect size was observed in gut microbiota when mice were switched from regular chow to purified diet in both sexes. However, sex and aging appeared to be independent drivers of gut microbial assemblage and should be taken into account in studies of this nature. In summary, we report a practical and effective pre-study SOP for normalizing the gut microbiome of murine cohorts that minimizes inter-individual variability and resolves co-housing problems inherent to male mice. This SOP may increase quality, rigor, and reproducibility of data acquisition and analysis.

Vaginal microbiome in early pregnancy and subsequent risk of spontaneous preterm birth: a case-control study

OBJECTIVES

To explore differences in the vaginal microbiome between preterm and term deliveries.

DESIGN

Nested case-control study in 3D cohort (design, develop, discover).

SETTING

Quebec, Canada.

SAMPLE

Ninety-four women with spontaneous preterm birth as cases [17 early (<34 weeks) and 77 late (34-36 weeks) preterm birth] and 356 women as controls with term delivery (≥37 weeks).

METHODS

To assess the vaginal microbiome by sequencing the V4 region of the 16S ribosomal RNA (rRNA) gene in swabs self-collected during early pregnancy.

MAIN OUTCOME MEASURES

Comparison of relative abundance of bacterial operational taxonomic units and oligotypes and identifying vaginal community state types (CSTs) in early or late spontaneous preterm and term deliveries.

RESULTS

Lactobacillus gasseri/ Lactobacillus johnsonii (coefficient -5.36, 95% CI -8.07 to -2.65), Lactobacillus crispatus (99%)/ Lactobacillus acidophilus (99%) (-4.58, 95% CI -6.20 to -2.96), Lactobacillus iners (99%)/ Ralstonia solanacearum (99%) (-3.98, 95% CI -6.48 to -1.47) and Bifidobacterium longum/ Bifidobacterium breve (-8.84, 95% CI -12.96 to -4.73) were associated with decreased risk of early but not late preterm birth. Six vaginal CSTs were identified: four dominated by Lactobacillus; one with presence of bacterial vaginosis-associated bacteria (Gardnerella vaginalis, Atopobium vaginae and Veillonellaceae bacterium) (CST IV); and one with nondominance of Lactobacillus (CST VI). CST IV was associated with increased risk of early (4.22, 95% CI 1.24-24.85) but not late (1.63, 95% CI 0.68-5.04) preterm birth, compared with CST VI.

CONCLUSIONS

Lactobacillus gasseri/L. johnsonii, L. crispatus/L. acidophilus, L. iners/R. solanacearum and B. longum/B. breve may be associated with decreased risk of early preterm birth. A bacterial vaginosis-related vaginal CST versus a CST nondominated by Lactobacillus may be associated with increased risk of early preterm birth.

TWEETABLE ABSTRACT

Largest study of its kind finds certain species of vaginal Lactobacillus + Bifidobacterium may relate to lower risk of preterm birth.

Nitrogen-fixing populations of Planctomycetes and Proteobacteria are abundant in surface ocean metagenomes

Nitrogen fixation in the surface ocean impacts global marine nitrogen bioavailability and thus microbial primary productivity. Until now, cyanobacterial populations have been viewed as the main suppliers of bioavailable nitrogen in this habitat. Although PCR amplicon surveys targeting the nitrogenase reductase gene have revealed the existence of diverse non-cyanobacterial diazotrophic populations, subsequent quantitative PCR surveys suggest that they generally occur in low abundance. Here, we use state-of-the-art metagenomic assembly and binning strategies to recover nearly one thousand non-redundant microbial population genomes from the TARA Oceans metagenomes. Among these, we provide the first genomic evidence for non-cyanobacterial diazotrophs inhabiting surface waters of the open ocean, which correspond to lineages within the Proteobacteria and, most strikingly, the Planctomycetes. Members of the latter phylum are prevalent in aquatic systems, but have never been linked to nitrogen fixation previously. Moreover, using genome-wide quantitative read recruitment, we demonstrate that the discovered diazotrophs were not only widespread but also remarkably abundant (up to 0.3% of metagenomic reads for a single population) in both the Pacific Ocean and the Atlantic Ocean northwest. Our results extend decades of PCR-based gene surveys, and substantiate the importance of heterotrophic bacteria in the fixation of nitrogen in the surface ocean.