Microbiome profiling by illumina sequencing of combinatorial sequence-tagged PCR products

We developed a low-cost, high-throughput microbiome profiling method that uses combinatorial sequence tags attached to PCR primers that amplify the rRNA V6 region. Amplified PCR products are sequenced using an Illumina paired-end protocol to generate millions of overlapping reads. Combinatorial sequence tagging can be used to examine hundreds of samples with far fewer primers than is required when sequence tags are incorporated at only a single end. The number of reads generated permitted saturating or near-saturating analysis of samples of the vaginal microbiome. The large number of reads allowed an in-depth analysis of errors, and we found that PCR-induced errors composed the vast majority of non-organism derived species variants, an observation that has significant implications for sequence clustering of similar high-throughput data. We show that the short reads are sufficient to assign organisms to the genus or species level in most cases. We suggest that this method will be useful for the deep sequencing of any short nucleotide region that is taxonomically informative; these include the V3, V5 regions of the bacterial 16S rRNA genes and the eukaryotic V9 region that is gaining popularity for sampling protist diversity.     

Development of a 9600-clone procedure for oligonucleotide fingerprinting of rRNA genes: utilization to identify soil bacterial rRNA genes that correlate in abundance with the development of avocado root rot

Oligonucleotide fingerprinting of rRNA genes (OFRG) is an array-based method that generates microbial community profiles through analysis of rRNA gene clone libraries. The original OFRG method allowed 1536 clones to be analyzed per experiment. This report describes a procedure for analyzing 9600 clones per experiment, including a new probe set for bacterial analysis, and improved data processing and statistical analysis tools. The software tools are available at the OFRG website (). Use of the 9600-clone procedure was demonstrated by examining the bacterial rRNA gene compositions of soils subjected to various temperature treatments. These treatments produced a series of soils with a range of abilities to suppress avocado root rot, enabling the identification of bacterial rRNA genes that correlate in abundance with root rot suppressiveness. OFRG analysis of these soils produced 8876 bacterial rRNA gene fingerprints grouped into 5123 clusters, or operational taxonomic units (OTUs). Eleven OTUs exhibited a positive correlation between the number of clones and the percentage of healthy roots. An in silico analysis was performed to examine the relationship between the number of rRNA genes analyzed and the number of correlates (rRNA gene-avocado root rot symptoms) identified. As the number of clones decreased, fewer correlates were identified. To further increase the throughput of the OFRG method, use of a glass slide-fluorescent probe microarray format was also explored.     

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.     

Horizontal transfer of segments of the 16S rRNA genes between species of the Streptococcus anginosus group

The nature in variation of the 16S rRNA gene of members of the Streptococcus anginosus group was investigated by hybridization and DNA sequencing. A collection of 708 strains was analyzed by reverse line blot hybridization. This revealed the presence of distinct reaction patterns representing 11 different hybridization groups. The 16S rRNA genes of two strains of each hybridization group were sequenced to near-completion, and the sequence data confirmed the reverse line blot hybridization results. Closer inspection of the sequences revealed mosaic-like structures, strongly suggesting horizontal transfer of segments of the 16S rRNA gene between different species belonging to the Streptococcus anginosus group. Southern blot hybridization further showed that within a single strain all copies of the 16S rRNA gene had the same composition, indicating that the apparent mosaic structures were not PCR-induced artifacts. These findings indicate that the highly conserved rRNA genes are also subject to recombination and that these events may be fixed in the population. Such recombination may lead to the construction of incorrect phylogenetic trees based on the 16S rRNA genes.     

Analysis of bacterial community composition by oligonucleotide fingerprinting of rRNA genes

One of the first steps in characterizing an ecosystem is to describe the organisms inhabiting it. For microbial studies, experimental limitations have hindered the ability to depict diverse communities. Here we describe oligonucleotide fingerprinting of rRNA genes (OFRG), a method that permits identification of arrayed rRNA genes (rDNA) through a series of hybridization experiments using small DNA probes. To demonstrate this strategy, we examined the bacteria inhabiting two different soils. Analysis of 1,536 rDNA clones revealed 766 clusters grouped into five major taxa: Bacillus, Actinobacteria, Proteobacteria, and two undefined assemblages. Soil-specific taxa were identified and then independently confirmed through cluster-specific PCR of the original soil DNA. Near-species-level resolution was obtained by this analysis as clones with average sequence identities of 97% were grouped in the same cluster. A comparison of these OFRG results with the results obtained in a denaturing gradient gel electrophoresis analysis of the same two soils demonstrated the significance of this methodological advance. OFRG provides a cost-effective means to extensively analyze microbial communities and should have applications in medicine, biotechnology, and ecosystem studies.     

Serial analysis of V6 ribosomal sequence tags (SARST-V6): a method for efficient, high-throughput analysis of microbial community composition

Serial analysis of ribosomal sequence tags (SARST) is a novel technique for characterizing microbial community composition. The SARST method captures sequence information from concatemers of short 16S rDNA polymerase chain reaction (PCR) amplicons from complex populations of DNA. Here, we describe a similar method, serial analysis of V6 ribosomal sequence tags (SARST-V6), which targets the V6 hypervariable region of bacterial 16S rRNA genes. The SARST-V6 technique exploits internal primer sequences to generate compatible restriction digest overhangs, thereby improving upon the efficiency of SARST. Serial analysis of V6 ribosomal sequence tags of bacterial community composition in hydrothermal marine sediments from Guaymas Basin resembled results of cloning and sequencing of single, full-length PCR products from ribosomal RNA genes of the same microbial community. Both methods identified the same major bacterial groups, but only SARST-V6 recovered thermodesulfobacteria and gamma-proteobacteria sequences, while only full-length PCR product cloning recovered candidate division OP11 se-quences. There were differences in the relative frequencies of some phylotypes. The disparities reflect differences in the amplicon pool obtained during initial amplification that may result from different primer affinities or DNA degradation. These results demonstrate the utility of SARST-V6 in collecting taxonomically informative data for high-throughput analysis of microbial communities.     

Analysis of Bacterial Community Composition by Oligonucleotide Fingerprinting of rRNA Genes

One of the first steps in characterizing an ecosystem is to describe the organisms inhabiting it. For microbial studies, experimental limitations have hindered the ability to depict diverse communities. Here we describe oligonucleotide fingerprinting of rRNA genes (OFRG), a method that permits identification of arrayed rRNA genes (rDNA) through a series of hybridization experiments using small DNA probes. To demonstrate this strategy, we examined the bacteria inhabiting two different soils. Analysis of 1,536 rDNA clones revealed 766 clusters grouped into five major taxa: Bacillus, Actinobacteria, Proteobacteria, and two undefined assemblages. Soil-specific taxa were identified and then independently confirmed through cluster-specific PCR of the original soil DNA. Near-species-level resolution was obtained by this analysis as clones with average sequence identities of 97% were grouped in the same cluster. A comparison of these OFRG results with the results obtained in a denaturing gradient gel electrophoresis analysis of the same two soils demonstrated the significance of this methodological advance. OFRG provides a cost-effective means to extensively analyze microbial communities and should have applications in medicine, biotechnology, and ecosystem studies.     

Frequency of formation of chimeric molecules as a consequence of PCR coamplification of 16S rRNA genes from mixed bacterial genomes.

PCR is routinely used in amplification and cloning of rRNA genes from environmental DNA samples for studies of microbial community structure and identification of novel organisms. There have been concerns about generation of chimeric sequences as a consequence of PCR coamplification of highly conserved genes, because such sequences may lead to reports of nonexistent organisms. To quantify the frequency of chimeric molecule formation, mixed genomic DNAs from eight actinomycete species whose 16S rRNA sequences had been determined were used for PCR coamplification of 16S rRNA genes. A large number of cloned 16S ribosomal DNAs were examined by sequence analysis, and chimeric molecules were identified by multiple-sequence alignment with reference species. Here, we report that the level of occurrence of chimeric sequences after 30 cycles of PCR amplification was 32%. We also show that PCR-induced chimeras were formed between different rRNA gene copies from the same organism. Because of the wide use of PCR for direct isolation of 16S rRNA sequences from environmental DNA to assess microbial diversity, the extent of chimeric molecule formation deserves serious attention.     

Oligonucleotide fingerprinting of rRNA genes for analysis of fungal community composition

Thorough assessments of fungal diversity are currently hindered by technological limitations. Here we describe a new method for identifying fungi, oligonucleotide fingerprinting of rRNA genes (OFRG). ORFG sorts arrayed rRNA gene (ribosomal DNA [rDNA]) clones into taxonomic clusters through a series of hybridization experiments, each using a single oligonucleotide probe. A simulated annealing algorithm was used to design an OFRG probe set for fungal rDNA. Analysis of 1,536 fungal rDNA clones derived from soil generated 455 clusters. A pairwise sequence analysis showed that clones with average sequence identities of 99.2% were grouped into the same cluster. To examine the accuracy of the taxonomic identities produced by this OFRG experiment, we determined the nucleotide sequences for 117 clones distributed throughout the tree. For all but two of these clones, the taxonomic identities generated by this OFRG experiment were consistent with those generated by a nucleotide sequence analysis. Eighty-eight percent of the clones were affiliated with Ascomycota, while 12% belonged to BASIDIOMYCOTA: A large fraction of the clones were affiliated with the genera Fusarium (404 clones) and Raciborskiomyces (176 clones). Smaller assemblages of clones had high sequence identities to the Alternaria, Ascobolus, Chaetomium, Cryptococcus, and Rhizoctonia clades.     

Assessing the Fecal Microbiota: An Optimized Ion Torrent 16S rRNA Gene-Based Analysis Protocol

Assessing the distribution of 16S rRNA gene sequences within a biological sample represents the current state-of-the-art for determination of human gut microbiota composition. Advances in dissecting the microbial biodiversity of this ecosystem have very much been dependent on the development of novel high-throughput DNA sequencing technologies, like the Ion Torrent. However, the precise representation of this bacterial community may be affected by the protocols used for DNA extraction as well as by the PCR primers employed in the amplification reaction. Here, we describe an optimized protocol for 16S rRNA gene-based profiling of the fecal microbiota.     

东南极达尔克冰川-183m深处冰层融水细菌多样性研究

【目的】目前对于南极冰层微生物研究较少,而且研究手段多为纯培养和高通量测序,对于其中的微生物群落多样性仍知之甚少。本研究拟研究东南极达尔克冰川冰层中微生物群落组成。【方法】采用纯培养法、单细胞分选和高通量测序3种方法对东南极达尔克冰川冰层微生物进行研究,探究该生境中微生物的群落组成。【结果】从达尔克冰川中分离出10门19纲94属。其中,变形菌门(Proteobacteria)为优势菌门,α-变形菌纲(Alphaproteobacteria)为优势纲,鞘氨醇单胞菌属(Sphingomonas)为优势属,结果显示冰层中存在较为丰富的微生物多样性。其中,纯培养法分离出25株细菌。单细胞分选方法分离得到24株细菌。高通量测序共得到55 183条序列,聚类出116个操作分类单元(operational taxonomic unit, OTU)。3种研究方法得出的优势种群不尽相同。通过单细胞分选和纯培养法共获得7株细菌,它们与数据库最近源16S rRNA基因序列的相似度小于98.65%,其中有2株菌株与最近源16S rRNA基因序列的相似度小于95.00%,推测可能有2个潜在新属,5个潜在新种。【结论】本研究通过纯培养法、单细胞分选以及高通量测序3种方法对东南极达尔克冰川冰层微生物多样性进行研究发现,该生境中细菌多样性复杂。对比3种方法,其各有优势和局限性。这意味着结合使用多种研究方法研究微生物多样性,可获得更加全面的微生物群落组成。研究结果可为挖掘南极微生物资源提供数据基础。     

Abundant and diverse fungal microbiota in the murine intestine

Enteric microbiota play a variety of roles in intestinal health and disease. While bacteria in the intestine have been broadly characterized, little is known about the abundance or diversity of enteric fungi. This study utilized a culture-independent method termed oligonucleotide fingerprinting of rRNA genes (OFRG) to describe the compositions of fungal and bacterial rRNA genes from small and large intestines (tissue and luminal contents) of restricted-flora and specific-pathogen-free mice. OFRG analysis identified rRNA genes from all four major fungal phyla: Ascomycota, Basidiomycota, Chytridiomycota, and Zygomycota. The largest assemblages of fungal rRNA sequences were related to the genera Acremonium, Monilinia, Fusarium, Cryptococcus/Filobasidium, Scleroderma, Catenomyces, Spizellomyces, Neocallimastix, Powellomyces, Entophlyctis, Mortierella, and Smittium and the order Mucorales. The majority of bacterial rRNA gene clones were affiliated with the taxa Bacteroidetes, Firmicutes, Acinetobacter, and Lactobacillus. Sequence-selective PCR analyses also detected several of these bacterial and fungal rRNA genes in the mouse chow. Fluorescence in situ hybridization analysis with a fungal small-subunit rRNA probe revealed morphologically diverse microorganisms resident in the mucus biofilm adjacent to the cecal and proximal colonic epithelium. Hybridizing organisms comprised about 2% of the DAPI (4',6-diamidino-2-phenylindole, dihydrochloride)-positive organisms in the mucus biofilm, but their abundance in fecal material may be much lower. These data indicate that diverse fungal taxa are present in the intestinal microbial community. Their abundance suggests that they may play significant roles in enteric microbial functions.     

Phylogenetic screening of a bacterial,metagenomic library using homing endonuclease restriction and marker insertion

Metagenomics provides access to the uncultured majority of the microbial world. The approaches employed in this field have, however, had limited success in linking functional genes to the taxonomic or phylogenetic origin of the organism they belong to. Here we present an efficient strategy to recover environmental DNA fragments that contain phylogenetic marker genes from metagenomic libraries. Our method involves the cleavage of 23S ribsosmal RNA (rRNA) genes within pooled library clones by the homing endonuclease I-CeuI followed by the insertion and selection of an antibiotic resistance cassette. This approach was applied to screen a library of 6500 fosmid clones derived from the microbial community associated with the sponge Cymbastela concentrica. Several fosmid clones were recovered after the screen and detailed phylogenetic and taxonomic assignment based on the rRNA gene showed that they belong to previously unknown organisms. In addition, compositional features of these fosmid clones were used to classify and taxonomically assign a dataset of environmental shotgun sequences. Our approach represents a valuable tool for the analysis of rapidly increasing, environmental DNA sequencing information.     

Relative effect of glyphosate on glyphosate-tolerant maize rhizobacterial communities is not altered by soil properties

The rhizobacterial composition varies according to the soil properties. To test if the effect of herbicides on the rhizobacterial communities of genetically modified NK603 glyphosate-tolerant maize varies according to different soil locations, a comparison was made between the effects of glyphosate (Roundup Plus), a post-emergence applied herbicide, and a pre-emergence applied herbicide (GTZ) versus untreated soil. The potential effect was monitored by direct amplification, cloning, and sequencing of the soil DNA encoding 16S rRNA, and high-throughput DNA pyrosequencing of the bacterial DNA coding for the 16S rRNA hypervariable V6 region. The results obtained using three different methods to analyze the herbicide effect on the rhizobacterial communities of genetically modified NK603 maize were comparable to those previously obtained when glyphosate-tolerant maize was grown in soil with different characteristics. Both herbicides decreased the bacterial diversity in the rhizosphere, with Actinobacteria being the taxonomic group most affected. The results suggest that both herbicides affected the structure of the maize rhizobacterial community, but glyphosate was environmentally less aggressive.     

Phylogeny of Intestinal Ciliates,Including Charonina ventriculi,and Comparison of Microscopy and 18S rRNA Gene Pyrosequencing for Rumen Ciliate Community Structure Analysis

The development of high-throughput methods, such as the construction of 18S rRNA gene clone or pyrosequencing libraries, has allowed evaluation of ciliate community composition in hundreds of samples from the rumen and other intestinal habitats. However, several genera of mammalian intestinal ciliates have been described based only on morphological features and, to date, have not been identified using molecular methods. Here, we isolated single cells of one of the smallest but widely distributed intestinal ciliates, Charonina ventriculi, and sequenced its 18S rRNA gene. We verified the sequence in a full-cycle rRNA approach using fluorescence in situ hybridization and thereby assigned an 18S rRNA gene sequence to this species previously known only by its morphology. Based on its full-length 18S rRNA gene sequence, Charonina ventriculi was positioned within the phylogeny of intestinal ciliates in the subclass Trichostomatia. The taxonomic framework derived from this phylogeny was used for taxonomic assignment of trichostome ciliate 18S rRNA gene sequence data stemming from high-throughput amplicon pyrosequencing of rumen-derived DNA samples. The 18S rRNA gene-based ciliate community structure was compared to that obtained from microscopic counts using the same samples. Both methods allowed identification of dominant members of the ciliate communities and classification of the rumen ciliate community into one of the types first described by Eadie in 1962. Notably, each method is associated with advantages and disadvantages. Microscopy is a highly accurate method for evaluation of total numbers or relative abundances of different ciliate genera in a sample, while 18S rRNA gene pyrosequencing represents a valuable alternative for comparison of ciliate community structure in a large number of samples from different animals or treatment groups.     

DNA复合条形码在太白山土壤动物多样性研究中的应用

DNA复合条形码技术(metabarcoding)将DNA条形码与高通量测序技术相结合,快速便捷地鉴定群落混合样本中的物种,成为监测群落中物种组成和丰富度的可靠方法。采用这一方法分析了秦岭太白山5种不同生境的中小型土壤动物多样性,共得到土壤动物3门9纲28目199科。群落组成分析显示生境的变化对土壤动物群落组成有一定的影响。α多样性分析显示土壤动物群落丰富度指数最高的生境为针叶林,最低的为农田;土壤动物群落多样性指数最高的生境为针叶林,最低的为落叶小叶林。群落相似性分析显示高山草甸、针叶林和农田3种生境的土壤动物群落组成相似性较高,落叶小叶林和落叶阔叶林的土壤动物群落组成与这3种生境的差异较大,落叶小叶林与落叶阔叶林的土壤动物群落组成差异也较大。     

A direct PCR approach to accelerate analyses of human-associated microbial communities

Since the composition of the human microbiome is highly variable both within and between individuals, researchers are increasingly reliant on high-throughput molecular approaches to identify linkages between the composition of these communities and human health. While new sequencing technologies have made it increasingly feasible to analyze large numbers of human-associated samples, the extraction of DNA from samples often remains a bottleneck in the process. Here we tested a direct PCR approach using the Extract-N-Amp Plant PCR Kit to accelerate the 16S rRNA gene-based analyses of human-associated bacterial communities, directly comparing this method to a more commonly-used approach whereby DNA is first extracted and purified from samples using a series of steps prior to PCR amplification. We used both approaches on replicate samples collected from each of five body habitats (tongue surface, feces, forehead skin, underarm skin, and forearm skin) from four individuals. With the exception of the tongue samples, there were few significant differences in the estimates of taxon richness or phylogenetic diversity obtained using the two approaches. Perhaps more importantly, there were no significant differences between the methods in their ability resolve body habitat differences or inter-individual differences in bacterial community composition and the estimates of the relative abundances of individual taxa were nearly identical with the two methods. Overall, the two methods gave very similar results and the direct PCR approach is clearly advantageous for many studies exploring the diversity and composition of human-associated bacterial communities given that large numbers of samples can be processed far more quickly and efficiently.     

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.