mTAGs: taxonomic profiling using degenerate consensus reference sequences of ribosomal RNA genes

Database size positively correlates with the loss of species-level taxonomic resolution for the 16S rRNA and other prokaryotic marker genes

For decades, the 16S rRNA gene has been used to taxonomically classify prokaryotic species and to taxonomically profile microbial communities. However, the 16S rRNA gene has been criticized for being too conserved to differentiate between distinct species. We argue that the inability to differentiate between species is not a unique feature of the 16S rRNA gene. Rather, we observe the gradual loss of species-level resolution for other nearly-universal prokaryotic marker genes as the number of gene sequences increases in reference databases. This trend was strongly correlated with how represented a taxonomic group was in the database and indicates that, at the gene-level, the boundaries between many species might be fuzzy. Through our study, we argue that any approach that relies on a single marker to distinguish bacterial taxa is fraught even if some markers appear to be discriminative in current databases.

Gut microbiota carcinogen metabolism causes distal tissue tumours

Exposure to environmental pollutants and human microbiome composition are important predisposition factors for tumour development1,2. Similar to drug molecules, pollutants are typically metabolized in the body, which can change their carcinogenic potential and affect tissue distribution through altered toxicokinetics3. Although recent studies demonstrated that human-associated microorganisms can chemically convert a wide range of xenobiotics and influence the profile and tissue exposure of resulting metabolites4,5, the effect of microbial biotransformation on chemical-induced tumour development remains unclear. Here we show that the depletion of the gut microbiota affects the toxicokinetics of nitrosamines, which markedly reduces the development and severity of nitrosamine-induced urinary bladder cancer in mice6,7. We causally linked this carcinogen biotransformation to specific gut bacterial isolates in vitro and in vivo using individualized bacterial culture collections and gnotobiotic mouse models, respectively. We tested gut communities from different human donors to demonstrate that microbial carcinogen metabolism varies between individuals and we showed that this metabolic activity applies to structurally related nitrosamine carcinogens. Altogether, these results indicate that gut microbiota carcinogen metabolism may be a contributing factor for chemical-induced carcinogenesis, which could open avenues to target the microbiome for improved predisposition risk assessment and prevention of cancer.

Metagenomic 18S rDNA reads revealed zonation of eukaryotic communities in the Yongle blue hole

The Yongle blue hole (YBH), situated in the South China Sea, represents a compelling subject of study in marine microbiology due to its unique redox-layered microbial ecosystems. However, the diversity and ecology of microbial eukaryotes within the YBH remains underexplored. This study endeavors to bridge this gap through the application of the in situ microbial filtration and fixation (ISMIFF) device to collect 0.22-30 μm microbial samples from 21 water layers of YBH. Subsequent extraction of 18S rRNA metagenomic reads of 21 metagenomes and 10 metatranscriptomes facilitated a comprehensive analysis of community structures. Findings revealed a pronounced superiority in the diversity and richness of eukaryotic microorganisms in the oxic zone compared to its suboxic and anoxic counterparts. Notably, Dinophyceae and Maxillopoda emerged as the predominant taxa based on the analysis of the 18S rRNA reads for the V4 and V9 regions, which showed stratification In their relative abundance and suggested their potential role in the thermo-halocline boundaries and oxic-anoxic interface. Specifically, In these eukaryotic microbial communities, Dinophyceae exhibited significant abundance at 20 m (20.01%) and 105 m (26.13%) water depths, while Maxillopoda was prevalent at 40 m (22.84%), 80 m (23.19%), and 100 m (15.42%) depths. A part of these organisms, identified as larvae and protists, were likely attracted by swarming chemosynthetic bacterial prey prevailing at the thermo-halocline boundaries and oxic-anoxic interface. Furthermore, the phylogenetic relationships of the major 18S operational taxonomic units (OTUs) showed a close adjacency to known species, except for three Dinophyceae OTUs. In conclusion, this study provides critical insights into the vertical distribution and transcriptional activity of <30-μm eukaryotic microbes, shedding light on the taxonomic novelty of eukaryotic microorganisms within the semi-enclosed blue holes.

Marine picoplankton metagenomes and MAGs from eleven vertical profiles obtained by the Malaspina Expedition

The Ocean microbiome has a crucial role in Earth's biogeochemical cycles. During the last decade, global cruises such as Tara Oceans and the Malaspina Expedition have expanded our understanding of the diversity and genetic repertoire of marine microbes. Nevertheless, there are still knowledge gaps regarding their diversity patterns throughout depth gradients ranging from the surface to the deep ocean. Here we present a dataset of 76 microbial metagenomes (MProfile) of the picoplankton size fraction (0.2-3.0 µm) collected in 11 vertical profiles covering contrasting ocean regions sampled during the Malaspina Expedition circumnavigation (7 depths, from surface to 4,000 m deep). The MProfile dataset produced 1.66 Tbp of raw DNA sequences from which we derived: 17.4 million genes clustered at 95% sequence similarity (M-GeneDB-VP), 2,672 metagenome-assembled genomes (MAGs) of Archaea and Bacteria (Malaspina-VP-MAGs), and over 100,000 viral genomic sequences. This dataset will be a valuable resource for exploring the functional and taxonomic connectivity between the photic and bathypelagic tropical and sub-tropical ocean, while increasing our general knowledge of the Ocean microbiome.

Disentangling the effects of sulfate and other seawater ions on microbial communities and greenhouse gas emissions in a coastal forested wetland

Seawater intrusion into freshwater wetlands causes changes in microbial communities and biogeochemistry, but the exact mechanisms driving these changes remain unclear. Here we use a manipulative laboratory microcosm experiment, combined with DNA sequencing and biogeochemical measurements, to tease apart the effects of sulfate from other seawater ions. We examined changes in microbial taxonomy and function as well as emissions of carbon dioxide, methane, and nitrous oxide in response to changes in ion concentrations. Greenhouse gas emissions and microbial richness and composition were altered by artificial seawater regardless of whether sulfate was present, whereas sulfate alone did not alter emissions or communities. Surprisingly, addition of sulfate alone did not lead to increases in the abundance of sulfate reducing bacteria or sulfur cycling genes. Similarly, genes involved in carbon, nitrogen, and phosphorus cycling responded more strongly to artificial seawater than to sulfate. These results suggest that other ions present in seawater, not sulfate, drive ecological and biogeochemical responses to seawater intrusion and may be drivers of increased methane emissions in soils that received artificial seawater addition. A better understanding of how the different components of salt water alter microbial community composition and function is necessary to forecast the consequences of coastal wetland salinization.

Ecology of Endozoicomonadaceae in three coral genera across the Pacific Ocean

Health and resilience of the coral holobiont depend on diverse bacterial communities often dominated by key marine symbionts of the Endozoicomonadaceae family. The factors controlling their distribution and their functional diversity remain, however, poorly known. Here, we study the ecology of Endozoicomonadaceae at an ocean basin-scale by sampling specimens from three coral genera (Pocillopora, Porites, Millepora) on 99 reefs from 32 islands across the Pacific Ocean. The analysis of 2447 metabarcoding and 270 metagenomic samples reveals that each coral genus harbored a distinct new species of Endozoicomonadaceae. These species are composed of nine lineages that have distinct biogeographic patterns. The most common one, found in Pocillopora, appears to be a globally distributed symbiont with distinct metabolic capabilities, including the synthesis of amino acids and vitamins not produced by the host. The other lineages are structured partly by the host genetic lineage in Pocillopora and mainly by the geographic location in Porites. Millepora is more rarely associated to Endozoicomonadaceae. Our results show that different coral genera exhibit distinct strategies of host-Endozoicomonadaceae associations that are defined at the bacteria lineage level.

Cultivation-independent genomes greatly expand taxonomic-profiling capabilities of mOTUs across various environments

BACKGROUND

Taxonomic profiling is a fundamental task in microbiome research that aims to detect and quantify the relative abundance of microorganisms in biological samples. Available methods using shotgun metagenomic data generally depend on the deposition of sequenced and taxonomically annotated genomes, usually from cultures of isolated strains, in reference databases (reference genomes). However, the majority of microorganisms have not been cultured yet. Thus, a substantial fraction of microbial community members remains unaccounted for during taxonomic profiling, particularly in samples from underexplored environments. To address this issue, we developed the mOTU profiler, a tool that enables reference genome-independent species-level profiling of metagenomes. As such, it supports the identification and quantification of both "known" and "unknown" species based on a set of select marker genes.

RESULTS

We present mOTUs3, a command line tool that enables the profiling of metagenomes for >33,000 species-level operational taxonomic units. To achieve this, we leveraged the reconstruction of >600,000 draft genomes, most of which are metagenome-assembled genomes (MAGs), from diverse microbiomes, including soil, freshwater systems, and the gastrointestinal tract of ruminants and other animals, which we found to be underrepresented by reference genomes. Overall, two thirds of all species-level taxa lacked a reference genome. The cumulative relative abundance of these newly included taxa was low in well-studied microbiomes, such as the human body sites (6-11%). By contrast, they accounted for substantial proportions (ocean, freshwater, soil: 43-63%) or even the majority (pig, fish, cattle: 60-80%) of the relative abundance across diverse non-human-associated microbiomes. Using community-developed benchmarks and datasets, we found mOTUs3 to be more accurate than other methods and to be more congruent with 16S rRNA gene-based methods for taxonomic profiling. Furthermore, we demonstrate that mOTUs3 increases the resolution of well-known microbial groups into species-level taxa and helps identify new differentially abundant taxa in comparative metagenomic studies.

CONCLUSIONS

We developed mOTUs3 to enable accurate species-level profiling of metagenomes. Compared to other methods, it provides a more comprehensive view of prokaryotic community diversity, in particular for currently underexplored microbiomes. To facilitate comparative analyses by the research community, it is released with >11,000 precomputed profiles for publicly available metagenomes and is freely available at: https://github.com/motu-tool/mOTUs . Video Abstract.