Comparative metagenomic and rRNA microbial diversity characterization using archaeal and bacterial synthetic communities

Next‐generation sequencing has dramatically changed the landscape of microbial ecology, large‐scale and in‐depth diversity studies being now widely accessible. However, determining the accuracy of taxonomic and quantitative inferences and comparing results obtained with different approaches are complicated by incongruence of experimental and computational data types and also by lack of knowledge of the true ecological diversity. Here we used highly diverse bacterial and archaeal synthetic communities assembled from pure genomic DNAs to compare inferences from metagenomic and SSU rRNA amplicon sequencing. Both Illumina and 454 metagenomic data outperformed amplicon sequencing in quantifying the community composition, but the outcome was dependent on analysis parameters and platform. New approaches in processing and classifying amplicons can reconstruct the taxonomic composition of the community with high reproducibility within primer sets, but all tested primers sets lead to significant taxon‐specific biases. Controlled synthetic communities assembled to broadly mimic the phylogenetic richness in target environments can provide important validation for fine‐tuning experimental and computational parameters used to characterize natural communities.     

Reproducibility and quantitation of amplicon sequencing-based detection

To determine the reproducibility and quantitation of the amplicon sequencing-based detection approach for analyzing microbial community structure, a total of 24 microbial communities from a long-term global change experimental site were examined. Genomic DNA obtained from each community was used to amplify 16S rRNA genes with two or three barcode tags as technical replicates in the presence of a small quantity (0.1% wt/wt) of genomic DNA from Shewanella oneidensis MR-1 as the control. The technical reproducibility of the amplicon sequencing-based detection approach is quite low, with an average operational taxonomic unit (OTU) overlap of 17.2%±2.3% between two technical replicates, and 8.2%±2.3% among three technical replicates, which is most likely due to problems associated with random sampling processes. Such variations in technical replicates could have substantial effects on estimating β-diversity but less on α-diversity. A high variation was also observed in the control across different samples (for example, 66.7-fold for the forward primer), suggesting that the amplicon sequencing-based detection approach could not be quantitative. In addition, various strategies were examined to improve the comparability of amplicon sequencing data, such as increasing biological replicates, and removing singleton sequences and less-representative OTUs across biological replicates. Finally, as expected, various statistical analyses with preprocessed experimental data revealed clear differences in the composition and structure of microbial communities between warming and non-warming, or between clipping and non-clipping. Taken together, these results suggest that amplicon sequencing-based detection is useful in analyzing microbial community structure even though it is not reproducible and quantitative. However, great caution should be taken in experimental design and data interpretation when the amplicon sequencing-based detection approach is used for quantitative analysis of the β-diversity of microbial communities.     

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.     

A comparison of DNA profiling techniques for monitoring nutrient impact on microbial community composition during bioremediation of petroleum-contaminated soils

Amplicon length heterogeneity PCR (LH-PCR) and terminal restriction fragment length polymorphisms (TRFLP) were used to monitor the impact that nutrient amendments had on microbial community dynamics and structural diversity during bioremediation of petroleum-contaminated soils. Slurried soils contaminated with petroleum hydrocarbons were treated in airlift bench-scale bioreactors and were either amended with optimal inorganic nutrients or left unamended. Direct DNA extraction and PCR amplification of whole eubacterial community DNA were performed with universal primers that bracketed the first two or three hypervariable regions of the 16S rDNA gene sequences. The LH-PCR method profiled a more diverse microbial community than did the TRFLP method. The LH-PCR method also tracked differences between the communities due to nutrient amendments. An in silico database search for bacterial genera with amplicon lengths represented in the community fingerprints was performed. It was possible to qualitatively identify different groups in the microbial community based on the amplicon length variations. A similar "virtual" search was performed for the TRFLP fragments using the web-based TAP-TRFLP program. Cloning and sequencing of the PCR products confirmed the in silico database matches. The application of the LH-PCR method as a monitoring tool for bioremediation could greatly enhance and extend the current understanding of the microbial community dynamics during the biodegradation of environmental contaminants.     

Does peptide-nucleic acid (PNA) clamping of host plant DNA benefit ITS1 amplicon-based characterization of the fungal endophyte community?

Fungal endophyte community amplicon sequencing can lose a significant number of informative reads due to host-plant co-amplification. Blocking of plant-specific sequences with peptide nucleic acid (PNA) clamps has been shown to improve metrics of detected microbial diversity in studies targeting 16S and 18S regions of rRNA genes. However, PNA clamping has not been applied to the plant ITS region of rRNA gene – a widely accepted fungal marker. By applying PNA clamping technique to ITS amplicon sequencing of the endophytic fungal community of elderberry this study shows that PNA clamping significantly reduces host-plant co-amplification with the universal ITS1/ITS4 primer set. However, PNA clamping in combination with the discriminatory ITS1F/ITS2 primer set did not improve the metrics of fungal endophyte community ITS amplicon Illumina sequencing. This study shows that PNA clamping does not add practical benefit to taxonomic profiling of plant-associated fungal communities if the primers are already specific enough to exclude amplification of host DNA.     

Dissecting the dominant hot spring microbial populations based on community-wide sampling at single-cell genomic resolution

With advances in DNA sequencing and miniaturized molecular biology workflows, rapid and affordable sequencing of single-cell genomes has become a reality. Compared to 16S rRNA gene surveys and shotgun metagenomics, large-scale application of single-cell genomics to whole microbial communities provides an integrated snapshot of community composition and function, directly links mobile elements to their hosts, and enables analysis of population heterogeneity of the dominant community members. To that end, we sequenced nearly 500 single-cell genomes from a low diversity hot spring sediment sample from Dewar Creek, British Columbia, and compared this approach to 16S rRNA gene amplicon and shotgun metagenomics applied to the same sample. We found that the broad taxonomic profiles were similar across the three sequencing approaches, though several lineages were missing from the 16S rRNA gene amplicon dataset, likely the result of primer mismatches. At the functional level, we detected a large array of mobile genetic elements present in the single-cell genomes but absent from the corresponding same species metagenome-assembled genomes. Moreover, we performed a single-cell population genomic analysis of the three most abundant community members, revealing differences in population structure based on mutation and recombination profiles. While the average pairwise nucleotide identities were similar across the dominant species-level lineages, we observed differences in the extent of recombination between these dominant populations. Most intriguingly, the creek’s Hydrogenobacter sp. population appeared to be so recombinogenic that it more closely resembled a sexual species than a clonally evolving microbe. Together, this work demonstrates that a randomized single-cell approach can be useful for the exploration of previously uncultivated microbes from community composition to population structure.Subject terms: Environmental microbiology, Microbial ecology, Population dynamics     

Critical review of NGS analyses for de novo genotyping multigene families

The genotyping of highly polymorphic multigene families across many individuals used to be a particularly challenging task because of methodological limitations associated with traditional approaches. Next‐generation sequencing (NGS) can overcome most of these limitations, and it is increasingly being applied in population genetic studies of multigene families. Here, we critically review NGS bioinformatic approaches that have been used to genotype the major histocompatibility complex (MHC) immune genes, and we discuss how the significant advances made in this field are applicable to population genetic studies of gene families. Increasingly, approaches are introduced that apply thresholds of sequencing depth and sequence similarity to separate alleles from methodological artefacts. We explain why these approaches are particularly sensitive to methodological biases by violating fundamental genotyping assumptions. An alternative strategy that utilizes ultra‐deep sequencing (hundreds to thousands of sequences per amplicon) to reconstruct genotypes and applies statistical methods on the sequencing depth to separate alleles from artefacts appears to be more robust. Importantly, the ‘degree of change’ (DOC) method avoids using arbitrary cut‐off thresholds by looking for statistical boundaries between the sequencing depth for alleles and artefacts, and hence, it is entirely repeatable across studies. Although the advances made in generating NGS data are still far ahead of our ability to perform reliable processing, analysis and interpretation, the community is developing statistically rigorous protocols that will allow us to address novel questions in evolution, ecology and genetics of multigene families. Future developments in third‐generation single molecule sequencing may potentially help overcome problems that still persist in de novo multigene amplicon genotyping when using current second‐generation sequencing approaches.     

Rapid phylogenetic dissection of prokaryotic community structure in tidal flat using pyrosequencing

Dissection of prokaryotic community structure is prerequisite to understand their ecological roles. Various methods are available for such a purpose which amplification and sequencing of 16S rRNA genes gained its popularity. However, conventional methods based on Sanger sequencing technique require cloning process prior to sequencing, and are expensive and labor-intensive. We investigated prokaryotic community structure in tidal flat sediments, Korea, using pyrosequencing and a subsequent automated bioinformatic pipeline for the rapid and accurate taxonomic assignment of each amplicon. The combination of pyrosequencing and bioinformatic analysis showed that bacterial and archaeal communities were more diverse than previously reported in clone library studies. Pyrosequencing analysis revealed 21 bacterial divisions and 37 candidate divisions. Proteobacteria was the most abundant division in the bacterial community, of which Gamma-and Delta-Proteobacteria were the most abundant. Similarly, 4 archaeal divisions were found in tidal flat sediments. Euryarchaeota was the most abundant division in the archaeal sequences, which were further divided into 8 classes and 11 unclassified euryarchaeota groups. The system developed here provides a simple, in-depth and automated way of dissecting a prokaryotic community structure without extensive pretreatment such as cloning.     

Resolving allele dosage in duplicated loci using genotyping‐by‐sequencing data: A path forward for population genetic analysis

Whole‐genome duplications have occurred in the recent ancestors of many plants, fish and amphibians. Signals of these whole‐genome duplications still exist in the form of paralogous loci. Recent advances have allowed reliable identification of paralogs in genotyping‐by‐sequencing (GBS) data such as that generated from restriction‐site‐associated DNA sequencing (RADSeq); however, excluding paralogs from analyses is still routine due to difficulties in genotyping. This exclusion of paralogs may filter a large fraction of loci, including loci that may be adaptively important or informative for population genetic analyses. We present a maximum‐likelihood method for inferring allele dosage in paralogs and assess its accuracy using simulated GBS, empirical RADSeq and amplicon sequencing data from Chinook salmon. We accurately infer allele dosage for some paralogs from a RADSeq data set and show how accuracy is dependent upon both read depth and allele frequency. The amplicon sequencing data set, using RADSeq‐derived markers, achieved sufficient depth to infer allele dosage for all paralogs. This study demonstrates that RADSeq locus discovery combined with amplicon sequencing of targeted loci is an effective method for incorporating paralogs into population genetic analyses.     

Metabarcoding as a tool for investigating arthropod diversity in Nepenthes pitcher plants

The biodiversity of tropical forests consists primarily of small organisms that are difficult to detect and characterize. Next‐generation sequencing (NGS) methods can facilitate analyses of these arthropod and microbial communities, leading to a better understanding of existing diversity and factors influencing community assembly. The pitchers of carnivorous pitcher plants often house surprisingly discrete communities and provide ideal systems for analysis using an NGS approach. The plants digest insects in order to access essential nutrients while growing in poor soils; however, the pitchers are also home to communities of living organisms, called inquilines. Certain arthropods appear to have coevolved with their pitcher plant hosts and are not found in other environments. We used Illumina amplicon sequencing of 18S rDNA to characterize the eukaryotes in three species of Nepenthes (Nepenthaceae) pitcher plants – N. gracilis, N. rafflesiana and N. ampullaria – in each of three different parks in Singapore. The data reveal an unexpected diversity of eukaryotes, significant differences in community diversity among host species, variation in host specificity of inquilines and the presence of gregarine parasites. Counts of whole inquiline arthropods from the first collection year were roughly correlated with scaled 18S sequence abundances, indicating that amplicon sequencing is an effective means of gauging community structure. We barcoded a subset of the dipteran larvae using COI primers, and the resulting phylogenetic tree is mostly congruent with that found using the 18S locus, with the exception of one of five morphospecies. For many 18S and COI sequences, the best BLASTn matches showed low sequence identity, illustrating the need for better databases of Southeast Asian dipterans. Finally, networks of core arthropods and their host species were used to investigate degree of host specificity across multiple hosts, and this revealed significant specialization of certain arthropod fauna.     

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.     

Marine microbial community structure assessed from combined metagenomic analysis and ribosomal amplicon deep-sequencing

The microbial taxonomic composition of the three domains of life in two coastal plankton samples was assessed by random total community metagenomic sequencing and PCR-based rDNA amplicon deep-sequencing in order to compare the resulting diversity and investigate possible limitations and complementarities of each method. The various universal primer sets, used to amplify different hypervariable rDNA regions, revealed the same major high-level taxonomic groups in Bacteria and unicellular Eukaryota, and showed a scarce Archaea apparent richness. However, significant differences were found between the different primer sets (p-value < 0.05, with the Kolmogorov–Smirnov test), regarding both operational taxonomic unit (OTU) richness and relative abundance of the major high-level taxonomic groups detected. Based on the metagenomic approach, the phylum Bacteroidetes dominated the prokaryotic community, followed by Proteobacteria, while the detected eukaryotic unicellular taxa belonged to the groups of Alveolata, Fungi, Chlorophyta, Stramenopiles and Phaeophyceae. These groups were found to carry genes typically found in microbial communities, which are linked to DNA, RNA and protein metabolism and the synthesis of nucleotides, amino acids, carbohydrates and vitamins. Although our findings suggest that the total community metagenomic approach can provide a more comprehensive picture of the planktonic microbial community structure, a number of issues associated with this approach emerged. These issues include the still relatively high cost compared to amplicon sequencing, the possible low coverage of the full marine diversity, the insufficiency of databases for other gene markers than the small subunit gene, and the bias towards bacterial sequences because of their higher abundance relative to eukaryotes in marine environments.     

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.

     

Experimental evidence of microbial inheritance in plants and transmission routes from seed to phyllosphere and root

While the environment is considered the primary origin of the plant microbiome, the potential role of seeds as a source of transmitting microorganisms has not received much attention. Here we tested the hypothesis that the plant microbiome is partially inherited through vertical transmission. An experimental culturing device was constructed to grow oak seedlings in a microbe-free environment while keeping belowground and aboveground tissues separated. The microbial communities associated with the acorn's embryo and pericarp and the developing seeding's phyllosphere and root systems were analysed using amplicon sequencing of fungal ITS and bacterial 16S rDNA. Results showed that the seed microbiome is diverse and non-randomly distributed within an acorn. The microbial composition of the phyllosphere was diverse and strongly resembled the composition found in the embryo, whereas the roots and pericarp each had a less diverse and distinct microbial community. Our findings demonstrate a high level of microbial diversity and spatial partitioning of the fungal and bacterial community within both seed and seedling, indicating inheritance, niche differentiation and divergent transmission routes for the establishment of root and phyllosphere communities.     

Analysis of 16S rRNA amplicon sequencing options on the Roche/454 next-generation titanium sequencing platform

Background

16S rRNA gene pyrosequencing approach has revolutionized studies in microbial ecology. While primer selection and short read length can affect the resulting microbial community profile, little is known about the influence of pyrosequencing methods on the sequencing throughput and the outcome of microbial community analyses. The aim of this study is to compare differences in output, ease, and cost among three different amplicon pyrosequencing methods for the Roche/454 Titanium platform

Methodology/Principal Findings

The following three pyrosequencing methods for 16S rRNA genes were selected in this study: Method-1 (standard method) is the recommended method for bi-directional sequencing using the LIB-A kit; Method-2 is a new option designed in this study for unidirectional sequencing with the LIB-A kit; and Method-3 uses the LIB-L kit for unidirectional sequencing. In our comparison among these three methods using 10 different environmental samples, Method-2 and Method-3 produced 1.5–1.6 times more useable reads than the standard method (Method-1), after quality-based trimming, and did not compromise the outcome of microbial community analyses. Specifically, Method-3 is the most cost-effective unidirectional amplicon sequencing method as it provided the most reads and required the least effort in consumables management.

Conclusions

Our findings clearly demonstrated that alternative pyrosequencing methods for 16S rRNA genes could drastically affect sequencing output (e.g. number of reads before and after trimming) but have little effect on the outcomes of microbial community analysis. This finding is important for both researchers and sequencing facilities utilizing 16S rRNA gene pyrosequencing for microbial ecological studies.     

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.