Search

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., (© 2024. The Author(s).)

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., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2024 The Author(s). Published by Elsevier Inc. All rights reserved.)

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., (© 2023. The Author(s), under exclusive licence to Springer Nature America, Inc.)

Coastal marine macrophytes exhibit some of the highest rates of primary productivity in the world. They have been found to host a diverse set of microbes, many of which may impact the biology of their hosts through metabolisms that are unique to microbial taxa. Here, we characterized the metabolic functions of macrophyte-associated microbial communities using metagenomes collected from 2 species of kelp ( Laminaria setchellii and Nereocystis luetkeana ) and 3 marine angiosperms ( Phyllospadix scouleri , P. serrulatus , and Zostera marina ), including the rhizomes of two surfgrass species ( Phyllospadix spp.), the seagrass Zostera marina , and the sediments surrounding P. scouleri and Z. marina . Using metagenomic sequencing, we describe 63 metagenome-assembled genomes (MAGs) that potentially benefit from being associated with macrophytes and may contribute to macrophyte fitness through their metabolic activity. Host-associated metagenomes contained genes for the use of dissolved organic matter from hosts and vitamin (B 1 , B 2 , B 7 , B 12 ) biosynthesis in addition to a range of nitrogen and sulfur metabolisms that recycle dissolved inorganic nutrients into forms more available to the host. The rhizosphere of surfgrass and seagrass contained genes for anaerobic microbial metabolisms, including nifH genes associated with nitrogen fixation, despite residing in a well-mixed and oxygenated environment. The range of oxygen environments engineered by macrophytes likely explains the diversity of both oxidizing and reducing microbial metabolisms and contributes to the functional capabilities of microbes and their influences on carbon and nitrogen cycling in nearshore ecosystems. IMPORTANCE Kelps, seagrasses, and surfgrasses are ecosystem engineers on rocky shorelines, where they show remarkably high levels of primary production. Through analysis of their associated microbial communities, we found a variety of microbial metabolisms that may benefit the host, including nitrogen metabolisms, sulfur oxidation, and the production of B vitamins. In turn, these microbes have the genetic capabilities to assimilate the dissolved organic compounds released by their macrophyte hosts. We describe a range of oxygen environments associated with surfgrass, including low-oxygen microhabitats in their rhizomes that host genes for nitrogen fixation. The tremendous productivity of coastal seaweeds and seagrasses is likely due in part to the activities of associated microbes, and an increased understanding of these associations is needed.

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., (© 2022 The Authors. Molecular Ecology Resources published by John Wiley & Sons Ltd.)

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., (© 2021. The Author(s).)

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.

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.

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., (© 2020 Chen et al.; Published by Cold Spring Harbor Laboratory Press.)

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.

Little is known about evolutionary drivers of microbial populations in the warm subseafloor of deep-sea hydrothermal vents. Here we reconstruct 73 metagenome-assembled genomes (MAGs) from two geochemically distinct vent fields in the Mid-Cayman Rise to investigate patterns of genomic variation within subseafloor populations. Low-abundance populations with high intra-population diversity coexist alongside high-abundance populations with low genomic diversity, with taxonomic differences in patterns of genomic variation between the mafic Piccard and ultramafic Von Damm vent fields. Populations from Piccard are significantly enriched in nonsynonymous mutations, suggesting stronger purifying selection in Von Damm relative to Piccard. Comparison of nine Sulfurovum MAGs reveals two high-coverage, low-diversity MAGs from Piccard enriched in unique genes related to the cellular membrane, suggesting these populations were subject to distinct evolutionary pressures that may correlate with genes related to nutrient uptake, biofilm formation, or viral invasion. These results are consistent with distinct evolutionary histories between geochemically different vent fields, with implications for understanding evolutionary processes in subseafloor microbial populations.

We introduce DESMAN for De novo Extraction of Strains from Metagenomes. Large multi-sample metagenomes are being generated but strain variation results in fragmentary co-assemblies. Current algorithms can bin contigs into metagenome-assembled genomes but are unable to resolve strain-level variation. DESMAN identifies variants in core genes and uses co-occurrence across samples to link variants into haplotypes and abundance profiles. These are then searched for against non-core genes to determine the accessory genome of each strain. We validated DESMAN on a complex 50-species 210-genome 96-sample synthetic mock data set and then applied it to the Tara Oceans microbiome.

Background: The epidemiology of bacterial vaginosis (BV) suggests it is sexually transmissible, yet no transmissible agent has been identified. It is probable that BV-associated bacterial communities are transferred from male to female partners during intercourse; however, the microbiota of sexual partners has not been well-studied., Results: Pyrosequencing analysis of PCR-amplified 16S rDNA was used to examine BV-associated bacteria in monogamous couples with and without BV using vaginal, male urethral, and penile skin specimens. The penile skin and urethral microbiota of male partners of women with BV was significantly more similar to the vaginal microbiota of their female partner compared to the vaginal microbiota of non-partner women with BV. This was not the case for male partners of women with normal vaginal microbiota. Specific BV-associated species were concordant in women with BV and their male partners., Conclusions: In monogamous heterosexual couples in which the woman has BV, the significantly higher similarity between the vaginal microbiota and the penile skin and urethral microbiota of the male partner, supports the hypothesis that sexual exchange of BV-associated bacterial taxa is common.

Objectives: Exploring associations between the gut microbiota and colonic inflammation and assessing sequential changes during exclusive enteral nutrition (EEN) may offer clues into the microbial origins of Crohn's disease (CD)., Methods: Fecal samples (n=117) were collected from 23 CD and 21 healthy children. From CD children fecal samples were collected before, during EEN, and when patients returned to their habitual diets. Microbiota composition and functional capacity were characterized using sequencing of the 16S rRNA gene and shotgun metagenomics., Results: Microbial diversity was lower in CD than controls before EEN (P=0.006); differences were observed in 36 genera, 141 operational taxonomic units (OTUs), and 44 oligotypes. During EEN, the microbial diversity of CD children further decreased, and the community structure became even more dissimilar than that of controls. Every 10 days on EEN, 0.6 genus diversity equivalents were lost; 34 genera decreased and one increased during EEN. Fecal calprotectin correlated with 35 OTUs, 14 of which accounted for 78% of its variation. OTUs that correlated positively or negatively with calprotectin decreased during EEN. The microbiota of CD patients had a broader functional capacity than healthy controls, but diversity decreased with EEN. Genes involved in membrane transport, sulfur reduction, and nutrient biosynthesis differed between patients and controls. The abundance of genes involved in biotin (P=0.005) and thiamine biosynthesis decreased (P=0.017), whereas those involved in spermidine/putrescine biosynthesis (P=0.031), or the shikimate pathway (P=0.058), increased during EEN., Conclusions: Disease improvement following treatment with EEN is associated with extensive modulation of the gut microbiome.