Serial analysis of ribosomal sequence tags (SARST): a high-throughput method for profiling complex microbial communities

Two decades of culture-independent studies have confirmed that microbial communities represent the most complex and concentrated pool of phylogenetic diversity on the planet. There remains a need for innovative molecular tools that can further our knowledge of microbial diversity and its functional implications. We present the method and application of serial analysis of ribosomal sequence tags (SARST) as a novel tool for elucidating complex microbial communities, such as those found in soils and sediments. Serial analysis of ribosomal sequence tags uses a series of enzymatic reactions to amplify and ligate ribosomal sequence tags (RSTs) from bacterial small subunit rRNA gene (SSU rDNA) V1-regions into concatemers that are cloned and sequenced. This approach offers a significant increase in throughput over traditional SSU rDNA clone libraries, as up to 20 RSTs are obtained from each sequencing reaction. To test SARST and measure the bias associated with this approach, RST libraries were prepared from a defined mixture of pure cultures and from duplicate arctic soil DNA samples. The actual RST distribution reflected the theoretical composition of the original defined mixture. Data from duplicate soil libraries (1345 and 1217 RSTs, with 525 and 505 unique RSTs, respectively) indicated that replication provides a strongly correlated RST profile (r(2) = 0.80) and division-level distribution of RSTs (r(2) = 0.99). Using sequence data from abundant soil RSTs, we designed specific primers that successfully amplified a larger portion of the SSU rDNA for further phylogenetic analysis. These results suggest that SARST is a powerful approach for reproducible high-throughput profiling of microbial diversity amenable to medical, industrial or environmental microbiology applications.     

PCR扩增循环数对细菌群落多样性测序分析的影响

运用高通量测序技术分析复杂样品中微生物群落组成及变化趋势,已经成为目前微生物研究领域的热点之一。本研究以复杂土壤样品和应用范围较广的瘤胃食糜样品为对象,选取20、25和30三个扩增循环数分别对样品的16S r RNA基因的V3区进行扩增,然后进行文库构建和测序。最后通过数据分析比较不同的扩增循环数对细菌多样性测定结果的影响。结果表明,扩增循环数越多,捕获到的细菌数量和种类越多;但并非循环数越多,群落中的微生物组成比例最优。整体来看,当扩增循环数为25时,样品中物种的数量和组成是最优的。     

Shotgun microbial profiling of fossil remains

Millions to billions of DNA sequences can now be generated from ancient skeletal remains thanks to the massive throughput of next‐generation sequencing platforms. Except in cases of exceptional endogenous DNA preservation, most of the sequences isolated from fossil material do not originate from the specimen of interest, but instead reflect environmental organisms that colonized the specimen after death. Here, we characterize the microbial diversity recovered from seven c. 200‐ to 13 000‐year‐old horse bones collected from northern Siberia. We use a robust, taxonomy‐based assignment approach to identify the microorganisms present in ancient DNA extracts and quantify their relative abundance. Our results suggest that molecular preservation niches exist within ancient samples that can potentially be used to characterize the environments from which the remains are recovered. In addition, microbial community profiling of the seven specimens revealed site‐specific environmental signatures. These microbial communities appear to comprise mainly organisms that colonized the fossils recently. Our approach significantly extends the amount of useful data that can be recovered from ancient specimens using a shotgun sequencing approach. In future, it may be possible to correlate, for example, the accumulation of postmortem DNA damage with the presence and/or abundance of particular microbes.     

引物应用策略：基于2对古菌16S rRNA基因通用引物对环境样品古菌多样性分析

【目的】找到适宜的16S rRNA基因通用引物应用策略,应对复杂环境微生物多样性调查,尤其目前高速发展的高通量测序技术带来的巨大挑战。【方法】用Oligocheck软件分别将两对应试的古菌16S rRNA基因通用引物与RDP(Ribosomal database project)数据库中古菌16S rRNA基因序列进行匹配比对。用两对应试引物分别构建海洋沉积物样品的古菌16S rRNA基因文库。【结果】软件匹配结果显示引物f109/r958与目的基因的匹配程度高于引物f21/r958。该结果与古菌16S rRNA基因文库RFLP分析、古菌多样性指数分析结果相吻合。数据还表明,2对引物的综合文库能更好满足该沉积物样品的古菌多样性分析。【结论】选用与数据库中目的基因匹配性高的通用引物和多个引物的联合使用,可以有效提高环境样品微生物多样性调查的分辨率。     

Molecular analysis of carbon monoxide-oxidizing bacteria associated with recent Hawaiian volcanic deposits

Genomic DNA extracts from four sites at Kilauea Volcano were used as templates for PCR amplification of the large subunit (coxL) of aerobic carbon monoxide dehydrogenase. The sites included a 42-year-old tephra deposit, a 108-year-old lava flow, a 212-year-old partially vegetated ash-and-tephra deposit, and an approximately 300-year-old forest. PCR primers amplified coxL sequences from the OMP clade of CO oxidizers, which includes isolates such as Oligotropha carboxidovorans, Mycobacterium tuberculosis, and Pseudomonas thermocarboxydovorans. PCR products were used to create clone libraries that provide the first insights into the diversity and phylogenetic affiliations of CO oxidizers in situ. On the basis of phylogenetic and statistical analyses, clone libraries for each site were distinct. Although some clone sequences were similar to coxL sequences from known organisms, many sequences appeared to represent phylogenetic lineages not previously known to harbor CO oxidizers. On the basis of average nucleotide diversity and average pairwise difference, a forested site supported the most diverse CO-oxidizing populations, while an 1894 lava flow supported the least diverse populations. Neither parameter correlated with previous estimates of atmospheric CO uptake rates, but both parameters correlated positively with estimates of microbial biomass and respiration. Collectively, the results indicate that the CO oxidizer functional group associated with recent volcanic deposits of the remote Hawaiian Islands contains substantial and previously unsuspected diversity.     

16S rRNA基因在微生物生态学中的应用

16S rRNA(Small subunit ribosomal RNA)基因是对原核微生物进行系统进化分类研究时最常用的分子标志物(Biomarker),广泛应用于微生物生态学研究中。近些年来随着高通量测序技术及数据分析方法等的不断进步,大量基于16S rRNA基因的研究使得微生物生态学得到了快速发展,然而使用16S rRNA基因作为分子标志物时也存在诸多问题,比如水平基因转移、多拷贝的异质性、基因扩增效率的差异、数据分析方法的选择等,这些问题影响了微生物群落组成和多样性分析时的准确性。对当前使用16S rRNA基因分析微生物群落组成和多样性的进展情况做一总结,重点讨论当前存在的主要问题以及各种分析方法的发展,尤其是与高通量测序技术有关的实验和数据处理问题。     

免培养法研究野生川金丝猴肠道内生细菌多样性

【目的】了解野生川金丝猴(Rhinopithecus roxellana)肠道内生细菌的组成及其多样性。【方法】提取川金丝猴肠道内生细菌总DNA,选用细菌通用引物799F和1492R对总DNA进行16S rRNA基因特异性扩增,构建川金丝猴肠道内生细菌16S rRNA基因克隆文库,对阳性克隆进行限制性内切酶片段长度多态性(PCR-RFLP)分析,并对HaeⅢ酶切带谱菌株进行测序,构建系统发育树。【结果】根据酶切带谱分析和测序结果,将随机挑取的157个阳性克隆归为27个不同的可操作分类单元(OTUs)。系统发育分析表明这些克隆序列有62.10%属于厚壁菌门(Firmicutes),其中包括梭菌属(Clostridium)、Cellulosilyticum属、Robinsoniella属、Anaerofustis属、Blautia属和Anaerovorax属,有37.90%属于未培养细菌。【结论】川金丝猴肠道内生细菌多样性丰富,并且可能存在新的分类单元。     

Bacterial communities of Bartonella-positive fleas: diversity and community assembly patterns

We investigated the bacterial communities of nine Bartonella-positive fleas (n = 6 Oropsylla hirsuta fleas and n = 3 Oropsylla montana fleas), using universal primers, clone libraries, and DNA sequencing. DNA sequences were used to classify bacteria detected in a phylogenetic context, to explore community assembly patterns within individual fleas, and to survey diversity patterns in dominant lineages.     

Latent Dirichlet Allocation reveals spatial and taxonomic structure in a DNA‐based census of soil biodiversity from a tropical forest

High‐throughput sequencing of amplicons from environmental DNA samples permits rapid, standardized and comprehensive biodiversity assessments. However, retrieving and interpreting the structure of such data sets requires efficient methods for dimensionality reduction. Latent Dirichlet Allocation (LDA) can be used to decompose environmental DNA samples into overlapping assemblages of co‐occurring taxa. It is a flexible model‐based method adapted to uneven sample sizes and to large and sparse data sets. Here, we compare LDA performance on abundance and occurrence data, and we quantify the robustness of the LDA decomposition by measuring its stability with respect to the algorithm's initialization. We then apply LDA to a survey of 1,131 soil DNA samples that were collected in a 12‐ha plot of primary tropical forest and amplified using standard primers for bacteria, protists, fungi and metazoans. The analysis reveals that bacteria, protists and fungi exhibit a strong spatial structure, which matches the topographical features of the plot, while metazoans do not, confirming that microbial diversity is primarily controlled by environmental variation at the studied scale. We conclude that LDA is a sensitive, robust and computationally efficient method to detect and interpret the structure of large DNA‐based biodiversity data sets. We finally discuss the possible future applications of this approach for the study of biodiversity.     

Systematic Design of 18S rRNA Gene Primers for Determining Eukaryotic Diversity in Microbial Consortia

High-throughput sequencing of ribosomal RNA gene (rDNA) amplicons has opened up the door to large-scale comparative studies of microbial community structures. The short reads currently produced by massively parallel sequencing technologies make the choice of sequencing region crucial for accurate phylogenetic assignments. While for 16S rDNA, relevant regions have been well described, no truly systematic design of 18S rDNA primers aimed at resolving eukaryotic diversity has yet been reported. Here we used 31,862 18S rDNA sequences to design a set of broad-taxonomic range degenerate PCR primers. We simulated the phylogenetic information that each candidate primer pair would retrieve using paired- or single-end reads of various lengths, representing different sequencing technologies. Primer pairs targeting the V4 region performed best, allowing discrimination with paired-end reads as short as 150 bp (with 75% accuracy at genus level). The conditions for PCR amplification were optimised for one of these primer pairs and this was used to amplify 18S rDNA sequences from isolates as well as from a range of environmental samples which were then Illumina sequenced and analysed, revealing good concordance between expected and observed results. In summary, the reported primer sets will allow minimally biased assessment of eukaryotic diversity in different microbial ecosystems.     

metaBIT,an integrative and automated metagenomic pipeline for analysing microbial profiles from high‐throughput sequencing shotgun data

Micro‐organisms account for most of the Earth's biodiversity and yet remain largely unknown. The complexity and diversity of microbial communities present in clinical and environmental samples can now be robustly investigated in record times and prices thanks to recent advances in high‐throughput DNA sequencing (HTS). Here, we develop metaBIT, an open‐source computational pipeline automatizing routine microbial profiling of shotgun HTS data. Customizable by the user at different stringency levels, it performs robust taxonomy‐based assignment and relative abundance calculation of microbial taxa, as well as cross‐sample statistical analyses of microbial diversity distributions. We demonstrate the versatility of metaBIT within a range of published HTS data sets sampled from the environment (soil and seawater) and the human body (skin and gut), but also from archaeological specimens. We present the diversity of outputs provided by the pipeline for the visualization of microbial profiles (barplots, heatmaps) and for their characterization and comparison (diversity indices, hierarchical clustering and principal coordinates analyses). We show that metaBIT allows an automatic, fast and user‐friendly profiling of the microbial DNA present in HTS shotgun data sets. The applications of metaBIT are vast, from monitoring of laboratory errors and contaminations, to the reconstruction of past and present microbiota, and the detection of candidate species, including pathogens.     

Successional dynamics of soil fungal diversity along a restoration chronosequence post‐coal mining

Soil disruption from open‐cut mining practices can adversely impact microbial communities and the ecosystem services that they mediate. Despite this, assessment of impacts of soil disruption, and the subsequent recovery of microbial communities is rarely studied. Monitoring of ecological restoration success on mine sites has traditionally focused on vegetation; however, most plants rely, at least in part, on associations with soil fungi for enhanced nutrient and water acquisition. Here, we used high‐throughput phylogenetic marker gene sequencing to characterize the diversity of soil fungal communities along a restoration chronosequence ranging from 3 to 23 years at a rehabilitated mine site. We used nonmined analogue sites as a baseline for comparative purposes and examined the associations of soil fungal communities with soil physicochemical and aboveground vegetation variables. Fungal richness on rehabilitated sites was significantly larger than on nonmined sites, suggesting that mixing of topsoil during stockpiling resulted in a composite microbial community. Fungal community composition was significantly influenced by edaphic variables and the length of rehabilitation, with mined sites becoming more similar to nonmined sites over time. Fungal populations associated with ectomycorrhizae were relatively more abundant than those associated with arbuscular mycorrhizae and declined in response to disturbance, but recovered over time on the woody dominated sites indicating a strong coupling of these fungi with aboveground vegetation. Our data indicate that soil fungal diversity is a useful bioindicator of soil restoration in mined sites and may complement more traditional vegetation‐based surveys.     

Use of High Throughput Sequencing and Light Microscopy Show Contrasting Results in a Study of Phytoplankton Occurrence in a Freshwater Environment

Assessing phytoplankton diversity is of primary importance for both basic and applied ecological studies. Following the advances in molecular methods, phytoplankton studies are switching from using classical microscopy to high throughput sequencing approaches. However, methodological comparisons of these approaches have rarely been reported. In this study, we compared the two methods, using a unique dataset of multiple water samples taken from a natural freshwater environment. Environmental DNA was extracted from 300 water samples collected weekly during 20 years, followed by high throughput sequencing of amplicons from the 16S and 18S rRNA hypervariable regions. For each water sample, phytoplankton diversity was also estimated using light microscopy. Our study indicates that species compositions detected by light microscopy and 454 high throughput sequencing do not always match. High throughput sequencing detected more rare species and picoplankton than light microscopy, and thus gave a better assessment of phytoplankton diversity. However, when compared to light microscopy, high throughput sequencing of 16S and 18S rRNA amplicons did not adequately identify phytoplankton at the species level. In summary, our study recommends a combined strategy using both morphological and molecular techniques.     

Salt marsh sediment diversity: a test of the variability of the rare biosphere among environmental replicates

Much of the phylogenetic diversity in microbial systems arises from rare taxa that comprise the long tail of taxon rank distribution curves. This vast diversity presents a challenge to testing hypotheses about the effects of perturbations on microbial community composition because variability of rare taxa among environmental replicates may be sufficiently large that it would require a prohibitive degree of sequencing to discern differences between samples. In this study we used pyrosequencing of 16S rRNA tags to examine the diversity and within-site variability of salt marsh sediment bacteria. Our goal was to determine whether pyrosequencing could produce similar patterns in community composition among replicate environmental samples from the same location. We hypothesized that repeated sampling from the same location would produce different snapshots of the rare community due to incomplete sequencing of the taxonomically rich rare biosphere. We demonstrate that the salt marsh sediments we sampled contain a remarkably diverse array of bacterial taxa and, in contrast to our hypothesis, repeated sampling from within the same site produces reliably similar patterns in bacterial community composition, even among rare organisms. These results demonstrate that deep sequencing of 16s tags is well suited to distinguish site-specific similarities and differences among rare taxa and is a valuable tool for hypothesis testing in microbial ecology.     

Endophyte microbiome diversity in micropropagated Atriplex canescens and Atriplex torreyi var griffithsii

Microbial diversity associated with micropropagated Atriplex species was assessed using microscopy, isolate culturing, and sequencing. Light, electron, and confocal microscopy revealed microbial cells in aseptically regenerated leaves and roots. Clone libraries and tag-encoded FLX amplicon pyrosequencing (TEFAP) analysis amplified sequences from callus homologous to diverse fungal and bacterial taxa. Culturing isolated some seed borne endophyte taxa which could be readily propagated apart from the host. Microbial cells were observed within biofilm-like residues associated with plant cell surfaces and intercellular spaces. Various universal primers amplified both plant and microbial sequences, with different primers revealing different patterns of fungal diversity. Bacterial and fungal TEFAP followed by alignment with sequences from curated databases revealed 7 bacterial and 17 ascomycete taxa in A. canescens, and 5 bacterial taxa in A. torreyi. Additional diversity was observed among isolates and clone libraries. Micropropagated Atriplex retains a complex, intimately associated microbiome which includes diverse strains well poised to interact in manners that influence host physiology. Microbiome analysis was facilitated by high throughput sequencing methods, but primer biases continue to limit recovery of diverse sequences from even moderately complex communities.     

枯草芽胞杆菌菌肥对有机冬瓜根区土壤微生态的影响

【背景】微生物肥料已广泛应用于我国有机作物的种植,其对有机种植土壤微生态的影响尚需科学评测。【目的】高通量测序技术可用于精确分析土壤微生物群落,从细菌、真菌群落结构和多样性的角度阐释枯草芽胞杆菌菌肥对有机农田根区土壤微生物群落的影响。【方法】在有机农田轮作种植条件下,施用枯草芽胞杆菌菌肥后提取冬瓜根区土壤基因组DNA,通过PCR扩增建立文库,利用IlluminaMiSeq高通量测序技术,并结合相关生物信息学方法分析土壤细菌16SrRNA基因V3-V4区和真菌ITS1区的多样性指数及群落结构;测定根区土壤化学性质及酶活性,分析有机冬瓜果实品质,并作相关分析。【结果】从6个有机冬瓜根区土壤样本中获得14199个细菌操作分类单元(OTU)和3378个真菌OTU,细菌和真菌文库测序覆盖率分别在98%、99%以上。枯草芽胞杆菌菌肥会在一定程度上提高土壤细菌种群多样性而降低真菌种群多样性,丰富了细菌群落结构,但显著降低了真菌群落丰富度(P0.05);并减少了根区土壤特有细菌和真菌物种。变形菌门、厚壁菌门和放线菌门是优势细菌,子囊菌门是优势真菌;枯草芽胞杆菌菌肥会提高绿弯菌门和子囊菌门的相对丰度,比例分别为46.23%、10.01%;降低变形菌门和担子菌门的相对丰度,比例分别为11.14%、74.72%。枯草芽胞杆菌菌肥显著降低了土壤pH,显著提高了有机冬瓜果实总氨基酸、可溶性固形物等营养成分含量(P0.05)。【结论】施用枯草芽胞杆菌菌肥改变有机冬瓜根区土壤细菌和真菌的丰富度和多样性,降低了土壤pH,提高了有机冬瓜果实品质。     

Molecular Analysis of Carbon Monoxide-Oxidizing Bacteria Associated with Recent Hawaiian Volcanic Deposits

Genomic DNA extracts from four sites at Kilauea Volcano were used as templates for PCR amplification of the large subunit (coxL) of aerobic carbon monoxide dehydrogenase. The sites included a 42-year-old tephra deposit, a 108-year-old lava flow, a 212-year-old partially vegetated ash-and-tephra deposit, and an approximately 300-year-old forest. PCR primers amplified coxL sequences from the OMP clade of CO oxidizers, which includes isolates such as Oligotropha carboxidovorans, Mycobacterium tuberculosis, and Pseudomonas thermocarboxydovorans. PCR products were used to create clone libraries that provide the first insights into the diversity and phylogenetic affiliations of CO oxidizers in situ. On the basis of phylogenetic and statistical analyses, clone libraries for each site were distinct. Although some clone sequences were similar to coxL sequences from known organisms, many sequences appeared to represent phylogenetic lineages not previously known to harbor CO oxidizers. On the basis of average nucleotide diversity and average pairwise difference, a forested site supported the most diverse CO-oxidizing populations, while an 1894 lava flow supported the least diverse populations. Neither parameter correlated with previous estimates of atmospheric CO uptake rates, but both parameters correlated positively with estimates of microbial biomass and respiration. Collectively, the results indicate that the CO oxidizer functional group associated with recent volcanic deposits of the remote Hawaiian Islands contains substantial and previously unsuspected diversity.     

New DGGE strategies for the analyses of methanotrophic microbial communities using different combinations of existing 16S rRNA-based primers

Methane-oxidising microbial communities are studied intensively because of their importance for global methane cycling. A suite of molecular microbial techniques has been applied to the study of these communities. Denaturing gradient gel electrophoresis (DGGE) is a diversity screening tool combining high sample throughput with phylogenetic information of high resolution. The existing 16S rRNA-based DGGE assays available for methane-oxidising bacteria suffer from low-specificity, low phylogentic information due to the length of the amplified fragments and/or from lack of resolving power. In the present study we developed new combinations of existing primers and applied these on methane-oxidising microbial communities in a freshwater wetland marsh. The designed strategies comprised nested as well as direct amplification of environmental DNA. Successful application of direct amplification using combinations of universal and specific primers circumvents the nested designs currently used. All developed assays resulted in identical community profiles in wetland soil cores with Methylobacter sp. and Methylocystis sp.-related sequences. Changes in the occurrence of Methylobacter-related sequences with depth in the soil profile may be related to the decrease in methane-oxidizing activity.     

污水处理活性污泥微生物群落多样性研究

为研究污水处理活性污泥微生物多样性,提取了活性污泥宏基因组DNA,并采用细菌通用引物27F和1492R扩增了上海污泥厂活性污泥细菌16S rDNA片段,构建了细菌16S rDNA克隆文库,并对该文库中的微生物群落进行了分析。共获得200条高质量序列并建立系统发育树,结果显示活性污泥主要的细菌类群为变形菌门(Proteobacteria)(91.9%)、厚壁菌门(Firmicures)(4.6%)、拟杆菌门(Bacteroidetes)(2%)、绿弯菌门(Chloroflexi)(0.5%)、硝化螺菌门(Nitrospirae)(1%)。其中,明显的优势菌群为Alcaligenes feacalis(55%)、Pseudomonas aeruginosa(12.8%)和Stenotrophomonas(12.8%),优势菌的产酶能力在活性污泥中显示生态修复功能菌的作用。     

Evaluation of QIAamp DNA Stool Mini Kit for ecological studies of gut microbiota

Cell lysis efficiency and the quality of DNA extracts from complex bacterial ecosystems are two major concerns in molecular ecological studies of gut microbiota. In this study, we use PCR-denaturing gradient gel electrophoresis (DGGE) DNA profiling, random cloning and sequence analysis of 16S rRNA genes to compare the QIAamp DNA Stool Mini Kit with the bead beating technique in the preparation of DNA extracts from gut microbiota of pigs. We also developed a washing procedure that can release more than 93% of bacterial cells attached to the gut mucosa. Both the QIAamp kit and bead beating method lysed approximately 95% of bacterial cells. PCR-DGGE DNA profiles of ileal and cecal microbiota from both digesta and mucosa that were generated from the DNA extracts using the two methods were nearly identical. Random cloning and sequence analysis also demonstrated the high quality of DNA extracts using the two methods. Two random clone sets of 16S rRNA genes generated from the DNA extracts had a similar degree of bacterial diversity. Different preparations of DNA extract from a single sample using the QIAamp kit consistently produced similar PCR-DGGE DNA profiles with similarity indexes higher than 99%. Our data suggest the appropriateness of the QIAamp DNA Stool Mini Kit for the studies of gut microbial ecology and the effectiveness of the QIAamp kit in processing multiple samples for cell lysis and DNA extraction.