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.     

Comparison of DNA extraction methods for microbial community profiling with an application to pediatric bronchoalveolar lavage samples

Barcoded amplicon sequencing is rapidly becoming a standard method for profiling microbial communities, including the human respiratory microbiome. While this approach has less bias than standard cultivation, several steps can introduce variation including the type of DNA extraction method used. Here we assessed five different extraction methods on pediatric bronchoalveolar lavage (BAL) samples and a mock community comprised of nine bacterial genera to determine method reproducibility and detection limits for these typically low complexity communities. Additionally, using the mock community, we were able to evaluate contamination and select a relative abundance cut-off threshold based on the geometric distribution that optimizes the trade off between detecting bona fide operational taxonomic units and filtering out spurious ones. Using this threshold, the majority of genera in the mock community were predictably detected by all extraction methods including the hard-to-lyse Gram-positive genus Staphylococcus. Differences between extraction methods were significantly greater than between technical replicates for both the mock community and BAL samples emphasizing the importance of using a standardized methodology for microbiome studies. However, regardless of method used, individual patients retained unique diagnostic profiles. Furthermore, despite being stored as raw frozen samples for over five years, community profiles from BAL samples were consistent with historical culturing results. The culture-independent profiling of these samples also identified a number of anaerobic genera that are gaining acceptance as being part of the respiratory microbiome. This study should help guide researchers to formulate sampling, extraction and analysis strategies for respiratory and other human microbiome samples.     

Effect of PCR amplicon size on assessments of clone library microbial diversity and community structure

PCR-based surveys of microbial communities commonly use regions of the small-subunit ribosomal RNA (SSU rRNA) gene to determine taxonomic membership and estimate total diversity. Here we show that the length of the target amplicon has a significant effect on assessments of microbial richness and community membership. Using operational taxonomic unit (OTU)- and taxonomy-based tools, we compared the V6 hypervariable region of the bacterial SSU rRNA gene of three amplicon libraries of c. 100, 400 and 1000 base pairs (bp) from each of two hydrothermal vent fluid samples. We found that the smallest amplicon libraries contained more unique sequences, higher diversity estimates and a different community structure than the other two libraries from each sample. We hypothesize that a combination of polymerase dissociation, cloning bias and mispriming due to secondary structure accounts for the differences. While this relationship is not linear, it is clear that the smallest amplicon libraries contained more different types of sequences, and accordingly, more diverse members of the community. Because divergent and lower abundant taxa can be more readily detected with smaller amplicons, they may provide better assessments of total community diversity and taxonomic membership than longer amplicons in molecular studies of microbial communities.     

High-Density Universal 16S rRNA Microarray Analysis Reveals Broader Diversity than Typical Clone Library When Sampling the Environment

Molecular approaches aimed at detection of a broad-range of prokaryotes in the environment routinely rely on classifying heterogeneous 16S rRNA genes amplified by polymerase chain reaction (PCR) using primers with broad specificity. The general method of sampling and categorizing DNA has been to clone then sequence the PCR products. However, the number of clones required to adequately catalog the majority of taxa in a sample is unwieldy. Alternatively, hybridizing target sequences to a universal 16S rRNA gene microarray may provide a more rapid and comprehensive view of prokaryotic community composition. This study investigated the breadth and accuracy of a microarray in detecting diverse 16S rRNA gene sequence types compared to clone-and-sequencing using three environmental samples: urban aerosol, subsurface soil, and subsurface water. PCR products generated from universal 16S rRNA gene-targeted primers were classified by using either the clone-and-sequence method or by hybridization to a novel high-density microarray of 297,851 probes complementary to 842 prokaryotic subfamilies. The three clone libraries comprised 1391 high-quality sequences. Approximately 8% of the clones could not be placed into a known subfamily and were considered novel. The microarray results confirmed the majority of clone-detected subfamilies and additionally demonstrated greater amplicon diversity extending into phyla not observed by the cloning method. Sequences matching operational taxonomic units within the phyla Nitrospira, Planctomycetes, and TM7, which were uniquely detected by the array, were verified with specific primers and subsequent amplicon sequencing. Subfamily richness detected by the array corresponded well with nonparametric richness predictions extrapolated from clone libraries except in the water community where clone-based richness predictions were greatly exceeded. It was concluded that although the microarray is unreliable in identifying novel prokaryotic taxa, it reveals greater diversity in environmental samples than sequencing a typically sized clone library. Furthermore, the microarray allowed samples to be rapidly evaluated with replication, a significant advantage in studies of microbial ecology.     

Ironing out the wrinkles in the rare biosphere through improved OTU clustering

Deep sequencing of PCR amplicon libraries facilitates the detection of low‐abundance populations in environmental DNA surveys of complex microbial communities. At the same time, deep sequencing can lead to overestimates of microbial diversity through the generation of low‐frequency, error‐prone reads. Even with sequencing error rates below 0.005 per nucleotide position, the common method of generating operational taxonomic units (OTUs) by multiple sequence alignment and complete‐linkage clustering significantly increases the number of predicted OTUs and inflates richness estimates. We show that a 2% single‐linkage preclustering methodology followed by an average‐linkage clustering based on pairwise alignments more accurately predicts expected OTUs in both single and pooled template preparations of known taxonomic composition. This new clustering method can reduce the OTU richness in environmental samples by as much as 30–60% but does not reduce the fraction of OTUs in long‐tailed rank abundance curves that defines the rare biosphere.     

“Pitfalls in relative abundance estimation using eDNA metabarcoding”

Metabarcoding provides unparalleled information of overall ecosystem ecology and diversity, from occurrence to relative abundances, from micro‐ to macro‐organisms, including cryptic species and larval stages that would otherwise be impossible to identify. PCR is still the most cost‐effective and routinely used method to extract information from environmental samples coupled with high‐throughput sequencing (HTS). Nichols et al. ( 2018 ) metabarcoded nine soil environmental samples from Alaska and six mixtures of synthetic oligonucleotides, using a short fragment of the chloroplast trnL gene. The authors tested the effect of PCR replicates, different polymerase mixes and amplicon GC content on diversity estimates. Results show that occurrence data was not biased for any of the parameters tested and only insufficient PCR replicates could impair the detection of rare taxa. For abundance data, polymerases were found to differentially bias estimates of the relative amplicon abundances and this depended on GC content.     

Assessment of a 16S rRNA amplicon Illumina sequencing procedure for studying the microbiome of a symbiont‐rich aphid genus

The bacterial communities inhabiting arthropods are generally dominated by a few endosymbionts that play an important role in the ecology of their hosts. Rather than comparing bacterial species richness across samples, ecological studies on arthropod endosymbionts often seek to identify the main bacterial strains associated with each specimen studied. The filtering out of contaminants from the results and the accurate taxonomic assignment of sequences are therefore crucial in arthropod microbiome studies. We aimed here to validate an Illumina 16S rRNA gene sequencing protocol and analytical pipeline for investigating endosymbiotic bacteria associated with aphids. Using replicate DNA samples from 12 species (Aphididae: Lachninae, Cinara) and several controls, we removed individual sequences not meeting a minimum threshold number of reads in each sample and carried out taxonomic assignment for the remaining sequences. With this approach, we show that (i) contaminants accounted for a negligible proportion of the bacteria identified in our samples; (ii) the taxonomic composition of our samples and the relative abundance of reads assigned to a taxon were very similar across PCR and DNA replicates for each aphid sample; in particular, bacterial DNA concentration had no impact on the results. Furthermore, by analysing the distribution of unique sequences across samples rather than aggregating them into operational taxonomic units (OTUs), we gained insight into the specificity of endosymbionts for their hosts. Our results confirm that Serratia symbiotica is often present in Cinara species, in addition to the primary symbiont, Buchnera aphidicola. Furthermore, our findings reveal new symbiotic associations with Erwinia‐ and Sodalis‐related bacteria. We conclude with suggestions for generating and analysing 16S rRNA gene sequences for arthropod‐endosymbiont studies.     

Biogeographical distribution and ecological ranges of benthic cyanobacteria in East Antarctic lakes

For the first time, the cyanobacterial diversity from microbial mats in lakes of Eastern Antarctica was investigated using microscopic and molecular approaches. The present study assessed the biogeographical distribution of cyanobacteria in Antarctica. Five samples were taken from four lakes spanning a range of different ecological environments in Larsemann Hills, Vestfold Hills and Rauer Islands to evaluate the influence of lake characteristics on the cyanobacterial diversity. Seventeen morphospecies and 28 16S rRNA gene-based operational taxonomic units belonging to the Oscillatoriales, Nostocales and Chroococcales were identified. The internal transcribed spacer was evaluated to complement the 16S rRNA gene data and showed similar but more clear-cut tendencies. The molecular approach suggested that potential Antarctic endemic species, including a previously undiscovered diversity, are more abundant than has been estimated by morphological methods. Moreover, operational taxonomic units, also found outside Antarctica, were more widespread over the continent than potential endemics. The cyanobacterial diversity of the most saline lakes was found to differ from the others, and correlations between the sampling depth and the cyanobacterial communities can also be drawn. Comparison with database sequences illustrated the ubiquity of several cyanobacterial operational taxonomic units and their remarkable range of tolerance to harsh environmental conditions.     

Analysis of rhizobacterial communities in perennial Graminaceae from polluted water meadow soil, and screening of metal-resistant, potentially plant growth-promoting bacteria

This study investigates the impact of long-term heavy metal contamination on the culturable, heterotrophic, functional and genetic diversity of rhizobacterial communities of perennial grasses in water meadow soil. The culturable heterotrophic diversity was investigated by colony appearance on solid LB medium. Genetic diversity was measured as bands in denaturing gradient gel electrophoresis (DGGE) obtained directly from rhizosphere soil and rhizoplane DNA extracts, and from the corresponding culturable communities. In the two rhizospheric fractions the DGGE profiles of the direct DNA extracts were similar and stable among replicates, whereas in the enriched cultures the profiles of the fractions differed, but among the replicates they were similar. One hundred isolates were collected into 33 different operational taxonomic units by use of amplified internal transcribed spacers and into 19 heavy metal-resistant phenotypes. The phylogenetic position of strains belonging to 18 operational taxonomic units, representing more than 80% of the isolates, was determined by 16S rRNA gene sequencing. Several heavy metal-resistant strains were isolated from rhizoplane. Finally, metal-resistant rhizobacteria were tested for plant growth-promoting characteristics; some were found to contain 1-aminocyclopropane-1-carboxylic acid deaminase and/or to produce indole acetic acid and siderophores. Two strains resistant to cadmium and zinc, Pseudomonas tolaasii RP23 and Pseudomonas fluorescens RS9, had all three plant growth-promoting characteristics. Our findings suggest that bacteria can respond to soil metal contamination, and the described methodological approach appears promising for targeting potential plant growth-promoting rhizobacteria.     

肠道微生物测序研究中不同实验来源差异的综合评价以及应用

肠道微生物测序研究具有将微生物结果转化为人类健康的巨大潜力。16S扩增子测序和宏基因组鸟枪测序(whole-metagenome shotgun, WMS)是微生物组研究中的两大主要方法,各具优势。然而,研究样品的异质性、测序仪和文库制备方法的差异如何影响肠道微生物测序结果的可重复性仍有待深入研究。该研究旨在通过比较粪便样本中微生物组成的差异,为测序技术的选择提供参考标准。3种广泛采用的测序仪的结果显示,WMS法中技术重复相关性(r=0.94)较高,而生物学重复相关性(r=0.69)较低。Bray-Curtis距离表明,生物学重复的差异大于技术重复(P<0.001)。此外,16S和WMS数据集间具有明显的分类学图谱差异。研究结果表明,同质化是样品DNA提取前的一个必要步骤;测序仪对分类学变异的贡献小于文库制备方法。我们开发了经验贝叶斯方法,即在计算中“借用信息”,利用标准化数据和(非)参数先验分布分析批次效应参数,提高了16S和WMS之间的群体可比性,为进一步应用于融合分析已发表的16S和微生物数据集提供了依据。     

Characterization of the fecal bacteria communities of forage-fed horses by pyrosequencing of 16S rRNA V4 gene amplicons

The diversity of the equine fecal bacterial community was evaluated using pyrosequencing of 16S rRNA gene amplicons. Fecal samples were obtained from horses fed cool-season grass hay. Fecal bacteria were characterized by amplifying the V4 region of bacterial 16S rRNA gene. Of 5898 mean unique sequences, a mean of 1510 operational taxonomic units were identified in the four fecal samples. Equine fecal bacterial richness was higher than that reported in humans, but lower than that reported in either cattle feces or soil. Bacterial classified sequences were assigned to 16 phyla, of which 10 were present in all samples. The largest number of reads belonged to Firmicutes (43.7% of total bacterial sequences), Verrucomicrobia (4.1%), Proteobacteria (3.8%), and Bacteroidetes (3.7%). The less abundant Actinobacteria, Cyanobacteria, and TM7 phyla presented here have not been previously described in the gut contents or feces of horses. Unclassified sequences represented 38.1% of total bacterial sequences; therefore, the equine fecal microbiome diversity is likely greater than that described. This is the first study to characterize the fecal bacterial community in horses by the use of 16S rRNA gene amplicon pyrosequencing, expanding our knowledge of the fecal microbiota of forage-fed horses.     

A fast,highly sensitive double‐nested PCR‐based method to screen fish immunobiomes

Efficient methods for constructing 16S tag amplicon libraries for pyrosequencing are needed for the rapid and thorough screening of infectious bacterial diversity from host tissue samples. Here we have developed a double‐nested PCR methodology that generates 16S tag amplicon libraries from very small amounts of bacteria/host samples. This methodology was tested for 133 kidney samples from the lake whitefish Coregonus clupeaformis (Salmonidae) sampled in five different lake populations. The double‐nested PCR efficiency was compared with two other PCR strategies: single primer pair amplification and simple nested PCR. The double‐nested PCR was the only amplification strategy to provide highly specific amplification of bacterial DNA. The resulting 16S amplicon libraries were synthesized and pyrosequenced using 454 FLX technology to analyse the variation of pathogenic bacteria abundance. The proportion of the community sequenced was very high (Good’s coverage estimator; mean = 95.4%). Furthermore, there were no significant differences of sequence coverage among samples. Finally, the occurrence of chimeric amplicons was very low. Therefore, the double‐nested PCR approach provides a rapid, informative and cost‐effective method for screening fish immunobiomes and most likely applicable to other low‐density microbiomes as well.     

Evaluating the effects of laboratory protocols on eDNA detection probability for an endangered freshwater fish

The effectiveness and accuracy of detection using environmental DNA (eDNA) is dependent on understanding the influence laboratory methods such as DNA extraction and PCR strategies have on detection probability. Ideally choice of sampling and extraction method will maximize eDNA yield and detection probability. Determining the survey effort required to reach a satisfactory detection probability (via increased PCR replicates or more sampling) could compensate for a lower eDNA yield if the sampling and extraction method has other advantages for a study, species or system. I analysed the effect of three different sampling and extraction methods on eDNA yield, detection probability and PCR replication for detecting the endangered freshwater fish Macquaria australasica from water samples. The impact of eDNA concentration, PCR strategy, target amplicon size and two marker regions: 12S (a mitochondrial gene) and 18S (a nuclear gene) was also assessed. The choice of sampling and extraction method and PCR strategy, rather than amplicon size and marker region, had the biggest effect on detection probability and PCR replication. The PCR replication effort required to achieve a detection probability of 0.95, ranged from 2 to 6 PCR replicates depending on the laboratory method used. As all methods yielded eDNA from which M. australasica was detected using the three target amplicons, differences in eDNA yield and detection probability between the three methods could be mitigated by determining the appropriate PCR replication effort. Evaluating the effect sampling and extraction methods will have on the detection probability and determining the laboratory protocols and PCR replication required to maximize detection and minimize false positives and negatives is a useful first step for eDNA occupancy studies.     

Communities of wood-inhabiting fungi in dead pine logs along a geographical gradient in Japan

Fungi are the main agents of coarse woody debris decomposition in forest ecosystems. We examined the associations of environmental variables with fungal community structures in dead pine logs at 12 geographically distant sites using amplicon pyrosequencing of fungal ITS rDNA. A total of 575 operational taxonomic units (OTUs) were identified based on clustering at 97% similarity. Among the known fungal ecological groups, saprotrophic fungi generally showed highest frequency of occurrence and were positively associated with mean annual temperature (MAT) and log diameter. Wood decay fungi with unknown decay type were positively associated with pine wilt disease and negatively associated with log diameter. Ordination analysis of the 42 most prevalent OTUs showed that MAT and annual precipitation significantly explained the observed fungal community structure. These results suggested that climate conditions and site history differentially effect structure fungal communities in pine logs among different ecological groups.     

Microbiological Survey of the Human Gastric Ecosystem Using Culturing and Pyrosequencing Methods

Stomach mucosa biopsies and gastric juices samples of 12 healthy persons were analysed by culturing in selective- and non-selective-rich media. Microbial DNA from four mucosal samples was also amplified by nested PCR using universal bacterial primers, and the 16S rDNA amplicons pyrosequenced. The total number of cultivable microorganisms recovered from the samples ranged from 10² to 10⁴?cfu/g or ml. The isolates were identified at the species level by PCR amplification and sequencing of the 16S rDNA. Isolates belonged mainly to four genera; Propionibacterium, Lactobacillus, Streptococcus and Staphylococcus. A total of 15,622 high-quality 16S rDNA sequence reads were obtained by pyrosequencing from the four mucosal samples. Sequence analysis grouped the reads into 59 families and 69 genera, revealing wide bacterial diversity. Considerable differences in the composition of the gastric microbiota were observed among the subjects, although in all samples the most abundant operational taxonomic units belonged to Streptococcus, Propionibacterium and Lactobacillus. Comparison of the stomach microbiota with that present in other parts of the human gastrointestinal tract revealed distinctive microbial communities. This is the first study in which a combination of culture and culture-independent techniques has been used to explore the bacterial diversity of the human stomach.     

Interaction and assembly processes of abundant and rare microbial communities during a diatom bloom process

Diatom blooms can significantly influence the dynamics of microbial communities, yet little is known about the interaction and assembly mechanisms of abundant and rare taxa during bloom process. Here, using 16S rRNA gene amplicon sequencing, we investigated the co-occurrence patterns and assembly processes of abundant and rare microbial communities during an early spring diatom bloom in Xiangshan bay. Our results showed that α-diversity indices in the rare subcommunity (RS) were significantly higher than those in the abundant and common subcommunities. β-Diversity of the RS was the highest among three subcommunities, and the variation of β-diversity in the three subcommunities was mainly induced by species turnover, which was also the highest in the RS. The assembly of microbial communities was mainly driven by the neutral processes, but the roles of neutral processes might differ in each subcommunity. Co-occurrence network analysis revealed that abundant and common operational taxonomic units were more often located in central positions within the network. Most of the modules in the network were specific to a particular bloom stage, owing to the succession of Skeletonema costatum. Overall, these findings expand current understanding of the microbial interaction and assembly mechanisms in marine environment suffering harmful algal bloom disturbance.     

Genetic diversity in marine algal virus communities as revealed by sequence analysis of DNA polymerase genes.

Algal-virus-specific PCR primers were used to amplify DNA polymerase gene (pol) fragments (683 to 689 bp) from the virus-sized fraction (0.02 to 0.2 microns) concentrated from inshore and offshore water samples collected from the Gulf of Mexico. Algal-virus-like DNA pol genes were detected in five samples collected from the surface and deep chlorophyll maximum. PCR products from an offshore station were cloned, and the genetic diversity of 33 fragments was examined by restriction fragment length polymorphism and sequence analysis. The five different genotypes or operational taxonomic units (OTUs) that were identified on the basis of restriction fragment length polymorphism banding patterns were present in different relative abundances (9 to 34%). One clone from each OTU was sequenced, and phylogenetic analysis showed that all of the OTUs fell within the family Phycodnaviridae. Four of the OTUs fell within a group of viruses (MpV) which infect the photosynthetic picoplankter Micromonas pusilla. The genetic diversity among these genotypes was as large as that previously found for MpV isolates from different oceans. The remaining genotype formed its own clade between viruses which infect M. pusilla and Chrysochromulina brevifilum. These results imply that marine virus communities contain a diverse assemblage of MpV-like viruses, as well as other unknown members of the Phycodnaviridae.     

Exploring the accuracy of amplicon‐based internal transcribed spacer markers for a fungal community

With the continual improvement in high‐throughput sequencing technology and constant updates to fungal reference databases, the use of amplicon‐based DNA markers as a tool to reveal fungal diversity and composition in various ecosystems has become feasible. However, both primer selection and the experimental procedure require meticulous verification. Here, we computationally and experimentally evaluated the accuracy and specificity of three widely used or newly designed internal transcribed spacer (ITS) primer sets (ITS1F/ITS2, gITS7/ITS4 and 5.8S‐Fun/ITS4‐Fun). In silico evaluation revealed that primer coverage varied at different taxonomic levels due to differences in degeneracy and the location of primer sets. Using even and staggered mock community standards, we identified different proportions of chimeric and mismatch reads generated by different primer sets, as well as great variation in species abundances, suggesting that primer selection would affect the results of amplicon‐based metabarcoding studies. Choosing proofreading and high‐fidelity polymerase (KAPA HiFi) could significantly reduce the percentage of chimeric and mismatch sequences, further reducing inflation of operational taxonomic units. Moreover, for two types of environmental fungal communities, plant endophytic and soil fungi, it was demonstrated that the three primer sets could not reach a consensus on fungal community composition or diversity, and that primer selection, not experimental treatment, determines observed soil fungal community diversity and composition. Future DNA marker surveys should pay greater attention to potential primer effects and improve the experimental scheme to increase credibility and accuracy.     

Comparing diversity levels in environmental samples: DNA sequence capture and metabarcoding approaches using 18S and COI genes

Environmental DNA studies targeting multiple taxa using metabarcoding provide remarkable insights into levels of species diversity in any habitat. The main drawbacks are the presence of primer bias and difficulty in identifying rare species. We tested a DNA sequence‐capture method in parallel with the metabarcoding approach to reveal possible advantages of one method over the other. Both approaches were performed using the same eDNA samples and the same 18S and COI regions, followed by high throughput sequencing. Metabarcoded eDNA libraries were PCR amplified with one primer pair from 18S and COI genes. DNA sequence‐capture libraries were enriched with 3,639 baits targeting the same gene regions. We tested amplicon sequence variants (ASVs) and operational taxonomic units (OTUs) in silico approaches for both markers and methods, using for this purpose the metabarcoding data set. ASVs methods uncovered more species for the COI gene, whereas the opposite occurred for the 18S gene, suggesting that clustering reads into OTUs could bias diversity richness especially using 18S with relaxed thresholds. Additionally, metabarcoding and DNA sequence‐capture recovered 80%–90% of the control sample species. DNA sequence‐capture was 8x more expensive, nonetheless it identified 1.5x more species for COI and 13x more genera for 18S than metabarcoding. Both approaches offer reliable results, sharing ca. 40% species and 72% families and retrieve more taxa when nuclear and mitochondrial markers are combined. eDNA metabarcoding is quite well established and low‐cost, whereas DNA‐sequence capture for biodiversity assessment is still in its infancy, is more time‐consuming but provides more taxonomic assignments.