RiboFR-Seq: a novel approach to linking 16S rRNA amplicon profiles to metagenomes

16S rRNA amplicon analysis and shotgun metagenome sequencing are two main culture-independent strategies to explore the genetic landscape of various microbial communities. Recently, numerous studies have employed these two approaches together, but downstream data analyses were performed separately, which always generated incongruent or conflict signals on both taxonomic and functional classifications. Here we propose a novel approach, RiboFR-Seq (Ribosomal RNA gene flanking region sequencing), for capturing both ribosomal RNA variable regions and their flanking protein-coding genes simultaneously. Through extensive testing on clonal bacterial strain, salivary microbiome and bacterial epibionts of marine kelp, we demonstrated that RiboFR-Seq could detect the vast majority of bacteria not only in well-studied microbiomes but also in novel communities with limited reference genomes. Combined with classical amplicon sequencing and shotgun metagenome sequencing, RiboFR-Seq can link the annotations of 16S rRNA and metagenomic contigs to make a consensus classification. By recognizing almost all 16S rRNA copies, the RiboFR-seq approach can effectively reduce the taxonomic abundance bias resulted from 16S rRNA copy number variation. We believe that RiboFR-Seq, which provides an integrated view of 16S rRNA profiles and metagenomes, will help us better understand diverse microbial communities.     

Versatile methanotrophs form an active methane biofilter in the oxycline of a seasonally stratified coastal basin

The potential and drivers of microbial methane removal in the water column of seasonally stratified coastal ecosystems and the importance of the methanotrophic community composition for ecosystem functioning are not well explored. Here, we combined depth profiles of oxygen and methane with 16S rRNA gene amplicon sequencing, metagenomics and methane oxidation rates at discrete depths in a stratified coastal marine system (Lake Grevelingen, The Netherlands). Three amplicon sequence variants (ASVs) belonging to different genera of aerobic Methylomonadaceae and the corresponding three methanotrophic metagenome-assembled genomes (MOB-MAGs) were retrieved by 16S rRNA sequencing and metagenomic analysis, respectively. The abundances of the different methanotrophic ASVs and MOB-MAGs peaked at different depths along the methane oxygen counter-gradient and the MOB-MAGs show a quite diverse genomic potential regarding oxygen metabolism, partial denitrification and sulphur metabolism. Moreover, potential aerobic methane oxidation rates indicated high methanotrophic activity throughout the methane oxygen counter-gradient, even at depths with low in situ methane or oxygen concentration. This suggests that niche-partitioning with high genomic versatility of the present Methylomonadaceae might contribute to the functional resilience of the methanotrophic community and ultimately the efficiency of methane removal in the stratified water column of a marine basin.     

VITCOMIC2: visualization tool for the phylogenetic composition of microbial communities based on 16S rRNA gene amplicons and metagenomic shotgun sequencing

Background

The 16S rRNA gene-based amplicon sequencing analysis is widely used to determine the taxonomic composition of microbial communities. Once the taxonomic composition of each community is obtained, evolutionary relationships among taxa are inferred by a phylogenetic tree. Thus, the combined representation of taxonomic composition and phylogenetic relationships among taxa is a powerful method for understanding microbial community structure; however, applying phylogenetic tree-based representation with information on the abundance of thousands or more taxa in each community is a difficult task. For this purpose, we previously developed the tool VITCOMIC (VIsualization tool for Taxonomic COmpositions of MIcrobial Community), which is based on the genome-sequenced microbes’ phylogenetic information. Here, we introduce VITCOMIC2, which incorporates substantive improvements over VITCOMIC that were necessary to address several issues associated with 16S rRNA gene-based analysis of microbial communities.

Results

We developed VITCOMIC2 to provide (i) sequence identity searches against broad reference taxa including uncultured taxa; (ii) normalization of 16S rRNA gene copy number differences among taxa; (iii) rapid sequence identity searches by applying the graphics processing unit-based sequence identity search tool CLAST; (iv) accurate taxonomic composition inference and nearly full-length 16S rRNA gene sequence reconstructions for metagenomic shotgun sequencing; and (v) an interactive user interface for simultaneous representation of the taxonomic composition of microbial communities and phylogenetic relationships among taxa. We validated the accuracy of processes (ii) and (iv) by using metagenomic shotgun sequencing data from a mock microbial community.

Conclusions

The improvements incorporated into VITCOMIC2 enable users to acquire an intuitive understanding of microbial community composition based on the 16S rRNA gene sequence data obtained from both metagenomic shotgun and amplicon sequencing.

     

Shallow shotgun sequencing of the microbiome recapitulates 16S amplicon results and provides functional insights

Prevailing 16S rRNA gene-amplicon methods for characterizing the bacterial microbiome of wildlife are economical, but result in coarse taxonomic classifications, are subject to primer and 16S copy number biases, and do not allow for direct estimation of microbiome functional potential. While deep shotgun metagenomic sequencing can overcome many of these limitations, it is prohibitively expensive for large sample sets. Here we evaluated the ability of shallow shotgun metagenomic sequencing to characterize taxonomic and functional patterns in the faecal microbiome of a model population of feral horses (Sable Island, Canada). Since 2007, this unmanaged population has been the subject of an individual-based, long-term ecological study. Using deep shotgun metagenomic sequencing, we determined the sequencing depth required to accurately characterize the horse microbiome. In comparing conventional vs. high-throughput shotgun metagenomic library preparation techniques, we validate the use of more cost-effective laboratory methods. Finally, we characterize similarities between 16S amplicon and shallow shotgun characterization of the microbiome, and demonstrate that the latter recapitulates biological patterns first described in a published amplicon data set. Unlike for amplicon data, we further demonstrate how shallow shotgun metagenomic data provide useful insights regarding microbiome functional potential which support previously hypothesized diet effects in this study system.     

Survey of (Meta)genomic Approaches for Understanding Microbial Community Dynamics

Advancement in the next generation sequencing technologies has led to evolution of the field of genomics and metagenomics in a slim duration with nominal cost at precipitous higher rate. While metagenomics and genomics can be separately used to reveal the culture-independent and culture-based microbial evolution, respectively, (meta)genomics together can be used to demonstrate results at population level revealing in-depth complex community interactions for specific ecotypes. The field of metagenomics which started with answering “who is out there?” based on 16S rRNA gene has evolved immensely with the precise organismal reconstruction at species/strain level from the deeply covered metagenome data outweighing the need to isolate bacteria of which 99% are de facto non-cultivable. In this review we have underlined the appeal of metagenomic-derived genomes in providing insights into the evolutionary patterns, growth dynamics, genome/gene-specific sweeps, and durability of environmental pressures. We have demonstrated the use of culture-based genomics and environmental shotgun metagenome data together to elucidate environment specific genome modulations via metagenomic recruitments in terms of gene loss/gain, accessory and core-genome extent. We further illustrated the benefit of (meta)genomics in the understanding of infectious diseases by deducing the relationship between human microbiota and clinical microbiology. This review summarizes the technological advances in the (meta)genomic strategies using the genome and metagenome datasets together to increase the resolution of microbial population studies.     

Inverted microbial community stratification and spatial–temporal stability in hypersaline anaerobic sediments from the S’Avall solar salterns

The anaerobic hypersaline sediments of an ephemeral pond from the S’Avall solar salterns constituted an excellent study system because of their easy accessibility, as well as the analogy of their microbial assemblages with some known deep-sea hypersaline anaerobic brines. By means of shotgun metagenomics and 16S rRNA gene amplicon sequencing, the microbial composition of the sediment was shown to be stable in time and space. The communities were formed by prokaryote representatives with a clear inferred anaerobic metabolism, mainly related to the methane, sulfur and nitrate cycles. The most conspicuous finding was the inverted nature of the vertical stratification. Contrarily to what could be expected, a methanogenic archaeal metabolism was found to dominate in the upper layers, whereas Bacteria with fermentative and anaerobic respiration metabolisms increased with depth. We could demonstrate the methanogenic nature of the members of candidate lineages DHVE2 and MSBL1, which were present in high abundance in this system, and described, for the first time, viruses infecting these lineages. Members of the putatively active aerobic genera Salinibacter and Halorubrum were detected especially in the deepest layers for which we hypothesize that either oxygen could be sporadically available, or they could perform anaerobic metabolisms. We also report a novel repertoire of virus species thriving in these sediments, which had special relevance because of their lysogenic lifestyles.     

Microbial communities of Auka hydrothermal sediments shed light on vent biogeography and the evolutionary history of thermophily

Hydrothermal vents have been key to our understanding of the limits of life, and the metabolic and phylogenetic diversity of thermophilic organisms. Here we used environmental metagenomics combined with analysis of physicochemical data and 16S rRNA gene amplicons to characterize the sediment-hosted microorganisms at the recently discovered Auka vents in the Gulf of California. We recovered 325 metagenome assembled genomes (MAGs) representing 54 phyla, over 30% of those currently known, showing the microbial community in Auka hydrothermal sediments is highly diverse. 16S rRNA gene amplicon screening of 224 sediment samples across the vent field indicates that the MAGs retrieved from a single site are representative of the microbial community in the vent field sediments. Metabolic reconstruction of a vent-specific, deeply branching clade within the Desulfobacterota suggests these organisms metabolize sulfur using novel octaheme cytochrome-c proteins related to hydroxylamine oxidoreductase. Community-wide comparison between Auka MAGs and MAGs from Guaymas Basin revealed a remarkable 20% species-level overlap, suggestive of long-distance species transfer over 400 km and subsequent sediment colonization. Optimal growth temperature prediction on the Auka MAGs, and thousands of reference genomes, shows that thermophily is a trait that has evolved frequently. Taken together, our Auka vent field results offer new perspectives on our understanding of hydrothermal vent microbiology.     

Comparison of Methods to Identify Pathogens and Associated Virulence Functional Genes in Biosolids from Two Different Wastewater Treatment Facilities in Canada

The use of treated municipal wastewater residues (biosolids) as fertilizers is an attractive, inexpensive option for growers and farmers. Various regulatory bodies typically employ indicator organisms (fecal coliforms, E. coli and Salmonella) to assess the adequacy and efficiency of the wastewater treatment process in reducing pathogen loads in the final product. Molecular detection approaches can offer some advantages over culture-based methods as they can simultaneously detect a wider microbial species range, including non-cultivable microorganisms. However, they cannot directly assess the viability of the pathogens. Here, we used bacterial enumeration methods together with molecular methods including qPCR, 16S rRNA and cpn60 gene amplicon sequencing and shotgun metagenomic sequencing to compare pre- and post-treatment biosolids from two Canadian wastewater treatment plants (WWTPs). Our results show that an anaerobic digestion WWTP was unsuccessful at reducing the live indicator organism load (coliforms, generic E. coli and Salmonella) below acceptable regulatory criteria, while biosolids from a dewatering/pelletization WWTP met these criteria. DNA from other pathogens was detected by the molecular methods, but these species were considered less abundant. Clostridium DNA increased significantly following anaerobic digestion treatments. In addition to pathogen DNA, genes related to virulence and antibiotic resistance were identified in treated biosolids. Shotgun metagenomics revealed the widest range of pathogen DNA and, among the approaches used here, was the only approach that could access functional gene information in treated biosolids. Overall, our results highlight the potential usefulness of amplicon sequencing and shotgun metagenomics as complementary screening methods that could be used in parallel with culture-based methods, although more detailed comparisons across a wider range of sites would be needed.     

Testing the limits of 454 pyrotag sequencing: reproducibility, quantitative assessment and comparison to T-RFLP fingerprinting of aquifer microbes

The characterization of microbial community structure via 16S rRNA gene profiling has been greatly advanced in recent years by the introduction of amplicon pyrosequencing. The possibility of barcoding gives the opportunity to massively screen multiple samples from environmental or clinical sources for community details. However, an on-going debate questions the reproducibility and semi-quantitative rigour of pyrotag sequencing, similar to the early days of community fingerprinting. In this study we demonstrate the reproducibility of bacterial 454 pyrotag sequencing over biological and technical replicates of aquifer sediment bacterial communities. Moreover, we explore the potential of recovering specific template ratios via quantitatively defined template spiking to environmental DNA. We sequenced pyrotag libraries of triplicate sediment samples taken in annual sampling campaigns at a tar oil contaminated aquifer in Düsseldorf, Germany. The abundance of dominating lineages was highly reproducible with a maximal standard deviation of ~4% read abundance across biological, and ~2% across technical replicates. Our workflow also allows for the linking of read abundances within defined assembled pyrotag contigs to that of specific 'in vivo' fingerprinting signatures. Thus we demonstrate that both terminal restriction fragment length polymorphism (T-RFLP) analysis and pyrotag sequencing are capable of recovering highly comparable community structure. Overall diversity was roughly double in amplicon sequencing. Pyrotag libraries were also capable of linearly recovering increasing ratios (up to 20%) of 16S rRNA gene amendments from a pure culture of Aliivibrio fisheri spiked to sediment DNA. Our study demonstrates that 454 pyrotag sequencing is a robust and reproducible method, capable of reliably recovering template abundances and overall community structure within natural microbial communities.     

Comparing Apples and Oranges?: Next Generation Sequencing and Its Impact on Microbiome Analysis

Rapid advancements in sequencing technologies along with falling costs present widespread opportunities for microbiome studies across a vast and diverse array of environments. These impressive technological developments have been accompanied by a considerable growth in the number of methodological variables, including sampling, storage, DNA extraction, primer pairs, sequencing technology, chemistry version, read length, insert size, and analysis pipelines, amongst others. This increase in variability threatens to compromise both the reproducibility and the comparability of studies conducted. Here we perform the first reported study comparing both amplicon and shotgun sequencing for the three leading next-generation sequencing technologies. These were applied to six human stool samples using Illumina HiSeq, MiSeq and Ion PGM shotgun sequencing, as well as amplicon sequencing across two variable 16S rRNA gene regions. Notably, we found that the factor responsible for the greatest variance in microbiota composition was the chosen methodology rather than the natural inter-individual variance, which is commonly one of the most significant drivers in microbiome studies. Amplicon sequencing suffered from this to a large extent, and this issue was particularly apparent when the 16S rRNA V1-V2 region amplicons were sequenced with MiSeq. Somewhat surprisingly, the choice of taxonomic binning software for shotgun sequences proved to be of crucial importance with even greater discriminatory power than sequencing technology and choice of amplicon. Optimal N50 assembly values for the HiSeq was obtained for 10 million reads per sample, whereas the applied MiSeq and PGM sequencing depths proved less sufficient for shotgun sequencing of stool samples. The latter technologies, on the other hand, provide a better basis for functional gene categorisation, possibly due to their longer read lengths. Hence, in addition to highlighting methodological biases, this study demonstrates the risks associated with comparing data generated using different strategies. We also recommend that laboratories with particular interests in certain microbes should optimise their protocols to accurately detect these taxa using different techniques.     

Vaginotypes of the human vaginal microbiome

The human vaginal environment harbours a community of bacteria that plays an important role in maintaining vaginal health and in protecting this environment from various urogenital infections. This bacterial population, also known as vaginal microbiota, has been demonstrated to be dominated by members of the Lactobacillus genus. Several studies employing 16S rRNA gene-based amplicon sequencing have classified the vaginal microbiota into five distinct community state types (CSTs) or vaginotypes. To deepen our understanding of the vaginal microbiota we performed an in-depth meta-analysis of 1312 publicly available datasets concerning healthy vaginal microbiome information obtained by metagenomics sequencing. The analysis confirmed the predominance of taxa belonging to the Lactobacillus genus, followed by members of the genera Gardnerella, Vibrio and Atopobium. Moreover, the statistical robustness offered by this meta-analysis allowed us to disentangle the species-level composition of dominant and accessory taxa constituting each vaginotype and to revisit and refine the previously proposed CST classification. In addition, a functional characterization of the metagenomic datasets revealed particular genetic features associated with each assigned vaginotype.     

High-Resolution Melt Analysis for Rapid Comparison of Bacterial Community Compositions

In the study of bacterial community composition, 16S rRNA gene amplicon sequencing is today among the preferred methods of analysis. The cost of nucleotide sequence analysis, including requisite computational and bioinformatic steps, however, takes up a large part of many research budgets. High-resolution melt (HRM) analysis is the study of the melt behavior of specific PCR products. Here we describe a novel high-throughput approach in which we used HRM analysis targeting the 16S rRNA gene to rapidly screen multiple complex samples for differences in bacterial community composition. We hypothesized that HRM analysis of amplified 16S rRNA genes from a soil ecosystem could be used as a screening tool to identify changes in bacterial community structure. This hypothesis was tested using a soil microcosm setup exposed to a total of six treatments representing different combinations of pesticide and fertilization treatments. The HRM analysis identified a shift in the bacterial community composition in two of the treatments, both including the soil fumigant Basamid GR. These results were confirmed with both denaturing gradient gel electrophoresis (DGGE) analysis and 454-based 16S rRNA gene amplicon sequencing. HRM analysis was shown to be a fast, high-throughput technique that can serve as an effective alternative to gel-based screening methods to monitor microbial community composition.     

Reconstruction of ribosomal RNA genes from metagenomic data

Direct sequencing of environmental DNA (metagenomics) has a great potential for describing the 16S rRNA gene diversity of microbial communities. However current approaches using this 16S rRNA gene information to describe community diversity suffer from low taxonomic resolution or chimera problems. Here we describe a new strategy that involves stringent assembly and data filtering to reconstruct full-length 16S rRNA genes from metagenomicpyrosequencing data. Simulations showed that reconstructed 16S rRNA genes provided a true picture of the community diversity, had minimal rates of chimera formation and gave taxonomic resolution down to genus level. The strategy was furthermore compared to PCR-based methods to determine the microbial diversity in two marine sponges. This showed that about 30% of the abundant phylotypes reconstructed from metagenomic data failed to be amplified by PCR. Our approach is readily applicable to existing metagenomic datasets and is expected to lead to the discovery of new microbial phylotypes.     

Simultaneous analysis of seven 16S rRNA hypervariable gene regions increases efficiency in marine bacterial diversity detection

Environmental DNA sequencing is the gold standard to reveal microbial community structures. In most applications, a one-fragment PCR approach is applied to amplify a taxonomic marker gene, usually a hypervariable region of the 16S rRNA gene. We used a new reverse complement (RC)-PCR-based assay that amplifies seven out of the nine hypervariable regions of the 16S rRNA gene, to interrogate bacterial communities in sediment samples collected from different coastal marine sites with an impact gradient. In parallel, we employed a traditional one-fragment analysis of the hypervariable V3–V4 region to investigate whether the RC-PCR reveals more of the ‘unseen’ diversity obtained by the one-fragment approach. As a benchmark for the full deck of diversity, we subjected the samples to PCR-free metagenomic sequencing. None of the two PCR-based approaches recorded the full taxonomic repertoire obtained from the metagenomics datasets. However, the RC-PCR approach detected 2.8 times more bacterial genera compared to the near-saturation sequenced V3–V4 samples. RC-PCR is an ideal compromise between the standard one-fragment approach and metagenomics sequencing and may guide future environmental sequencing studies, in which bacterial diversity is a central subject.     

OTU Analysis Using Metagenomic Shotgun Sequencing Data

Because of technological limitations, the primer and amplification biases in targeted sequencing of 16S rRNA genes have veiled the true microbial diversity underlying environmental samples. However, the protocol of metagenomic shotgun sequencing provides 16S rRNA gene fragment data with natural immunity against the biases raised during priming and thus the potential of uncovering the true structure of microbial community by giving more accurate predictions of operational taxonomic units (OTUs). Nonetheless, the lack of statistically rigorous comparison between 16S rRNA gene fragments and other data types makes it difficult to interpret previously reported results using 16S rRNA gene fragments. Therefore, in the present work, we established a standard analysis pipeline that would help confirm if the differences in the data are true or are just due to potential technical bias. This pipeline is built by using simulated data to find optimal mapping and OTU prediction methods. The comparison between simulated datasets revealed a relationship between 16S rRNA gene fragments and full-length 16S rRNA sequences that a 16S rRNA gene fragment having a length >150 bp provides the same accuracy as a full-length 16S rRNA sequence using our proposed pipeline, which could serve as a good starting point for experimental design and making the comparison between 16S rRNA gene fragment-based and targeted 16S rRNA sequencing-based surveys possible.     

Methanogen diversity evidenced by molecular characterization of methyl coenzyme M reductase A (mcrA) genes in hydrothermal sediments of the Guaymas Basin

The methanogenic community in hydrothermally active sediments of Guaymas Basin (Gulf of California, Mexico) was analyzed by PCR amplification, cloning, and sequencing of methyl coenzyme M reductase (mcrA) and 16S rRNA genes. Members of the Methanomicrobiales and Methanosarcinales dominated the mcrA and 16S rRNA clone libraries from the upper 15 cm of the sediments. Within the H2/CO2- and formate-utilizing family Methanomicrobiales, two mcrA and 16S rRNA lineages were closely affiliated with cultured species of the genera Methanoculleus and Methanocorpusculum. The most frequently recovered mcrA PCR amplicons within the Methanomicrobiales did not branch with any cultured genera. Within the nutritionally versatile family Methanosarcinales, one 16S rRNA amplicon and most of the mcrA PCR amplicons were affiliated with the obligately acetate utilizing species Methanosaeta concilii. The mcrA clone libraries also included phylotypes related to the methyl-disproportionating genus Methanococcoides. However, two mcrA and two 16S rRNA lineages within the Methanosarcinales were unrelated to any cultured genus. Overall, the clone libraries indicate a diversified methanogen community that uses H2/CO2, formate, acetate, and methylated substrates. Phylogenetic affiliations of mcrA and 16S rRNA clones with thermophilic and nonthermophilic cultured isolates indicate a mixed mesophilic and thermophilic methanogen community in the surficial Guaymas sediments.     

Performance Comparison of Illumina and Ion Torrent Next-Generation Sequencing Platforms for 16S rRNA-Based Bacterial Community Profiling

High-throughput sequencing of the taxonomically informative 16S rRNA gene provides a powerful approach for exploring microbial diversity. Here we compare the performances of two common “benchtop” sequencing platforms, Illumina MiSeq and Ion Torrent Personal Genome Machine (PGM), for bacterial community profiling by 16S rRNA (V1-V2) amplicon sequencing. We benchmarked performance by using a 20-organism mock bacterial community and a collection of primary human specimens. We observed comparatively higher error rates with the Ion Torrent platform and report a pattern of premature sequence truncation specific to semiconductor sequencing. Read truncation was dependent on both the directionality of sequencing and the target species, resulting in organism-specific biases in community profiles. We found that these sequencing artifacts could be minimized by using bidirectional amplicon sequencing and an optimized flow order on the Ion Torrent platform. Results of bacterial community profiling performed on the mock community and a collection of 18 human-derived microbiological specimens were generally in good agreement for both platforms; however, in some cases, results differed significantly. Disparities could be attributed to the failure to generate full-length reads for particular organisms on the Ion Torrent platform, organism-dependent differences in sequence error rates affecting classification of certain species, or some combination of these factors. This study demonstrates the potential for differential bias in bacterial community profiles resulting from the choice of sequencing platform alone.