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.     

Shining a Light on Dark Sequencing: Characterising Errors in Ion Torrent PGM Data

The Ion Torrent Personal Genome Machine (PGM) is a new sequencing platform that substantially differs from other sequencing technologies by measuring pH rather than light to detect polymerisation events. Using re-sequencing datasets, we comprehensively characterise the biases and errors introduced by the PGM at both the base and flow level, across a combination of factors, including chip density, sequencing kit, template species and machine. We found two distinct insertion/deletion (indel) error types that accounted for the majority of errors introduced by the PGM. The main error source was inaccurate flow-calls, which introduced indels at a raw rate of 2.84% (1.38% after quality clipping) using the OneTouch 200 bp kit. Inaccurate flow-calls typically resulted in over-called short-homopolymers and under-called long-homopolymers. Flow-call accuracy decreased with consecutive flow cycles, but we also found significant periodic fluctuations in the flow error-rate, corresponding to specific positions within the flow-cycle pattern. Another less common PGM error, high frequency indel (HFI) errors, are indels that occur at very high frequency in the reads relative to a given base position in the reference genome, but in the majority of instances were not replicated consistently across separate runs. HFI errors occur approximately once every thousand bases in the reference, and correspond to 0.06% of bases in reads. Currently, the PGM does not achieve the accuracy of competing light-based technologies. However, flow-call inaccuracy is systematic and the statistical models of flow-values developed here will enable PGM-specific bioinformatics approaches to be developed, which will account for these errors. HFI errors may prove more challenging to address, especially for polymorphism and amplicon applications, but may be overcome by sequencing the same DNA template across multiple chips.     

Effects of Continuous Cropping of Rye on Soil Biota and Biochemistry

Long-term studies on the ecological effects of continuous rye cultivations carried out in Poland are summarized. It was shown that in continuous cropping of rye, despite the decrease of crop yields, no significant difference was observed in annual primary production rates compared with estimates found for rye fields cultivated in diversified crop rotation patterns. In continuous cultivation of rye fauna impoverishment was observed, while at the same time the number of crop pests increased. The animal groups reduced by continuous cropping of rye belong to Protozoa, Lumbricidae, Acarina, and winged insects. The cultivation of rye in continuous cropping influenced the increase of fungi with bacteria and actinomycetes becoming less abundant. In soils under continuous rye the average concentration of phenolic acids was 400% higher than under rye in diversified rotation. Five different phenolic acids were identified. The mean content of bound amino acids in soils under diversified rotation was 417?mg×kg^?1and in soils under continuous rye cropping was 371?mg×kg^?1. The energetic cost of total soil animal community maintenance calculated per unit of biomass under continuous cropping of rye was almost twice as high as under the diversified rotation. The accumulation of toxic metabolites caused negative effects among biota under continuous cropping of rye.     

Microbial 16S rRNA Ion Tag and community metagenome sequencing using the Ion Torrent (PGM) Platform

Here we demonstrate a cost effective and scalable microbial ecology sequencing platform using the Ion Torrent Personal Genome Machine (PGM). We assessed both PCR amplified 16S rRNA and shotgun metagenomic approaches and generated 100,000+ to 1,000,000+ reads using 'post-light' based sequencing technology within different sized semi-conductor chips. Further development of Golay barcoded Ion Tags allowed multiplex analyses of microbial communities with substantially reduced costs compared with platforms such as 454/GS-FLX. Using these protocols we assessed the bacterial and archaeal dynamics within covered anaerobic digesters used to treat piggery wastes. Analysis of these sequence data showed that these novel methanogenic waste treatment systems are dominated by bacterial taxa, in particular Clostridium, Synergistia and Bacteroides that were maintained as a stable community over extended time periods. Archaeal community dynamics were more stochastic with the key methanogenic taxa more difficult to resolve, principally due to the poor congruence seen between community structures generated either by nested PCR or metagenomic approaches for archaeal analyses. Our results show that for microbial community structure and function analyses, the PGM platform provides a low cost, scalable and high throughput solution for both Tag sequencing and metagenomic analyses.     

Pyrosequencing Reveals Contrasting Soil Bacterial Diversity and Community Structure of Two Main Winter Wheat Cropping Systems in China

Microbes are key components of the soil environment, playing an important role in maintaining soil health, sustainability, and productivity. The composition and structure of soil bacterial communities were examined in winter wheat–rice (WR) and winter wheat–maize (WM) cropping systems derived from five locations in the Low-Middle Yangtze River plain and the Huang-Huai-Hai plain by pyrosequencing of the 16S ribosomal RNA gene amplicons. A total of 102,367 high quality sequences were used for multivariate statistical analysis and to test for correlation between community structure and environmental variables such as crop rotations, soil properties, and locations. The most abundant phyla across all soil samples were Proteobacteria, Acidobacteria, and Bacteroidetes. Similar patterns of bacterial diversity and community structure were observed within the same cropping systems, and a higher relative abundance of anaerobic bacteria was found in WR compared to WM cropping systems. Variance partitioning analysis revealed complex relationships between bacterial community and environmental variables. The effect of crop rotations was low but significant, and interactions among soil properties, locations, and crop rotations accounted for most of the explained variation in the structure of bacterial communities. Soil properties such as pH, available P, and available K showed higher correlations (positive or negative) with the majority of the abundant taxa. Bacterial diversity (the Shannon index) and richness (Chao1 and ACE) were higher under WR than WM cropping systems.     

Depicting more accurate pictures of protistan community complexity using pyrosequencing of hypervariable SSU rRNA gene regions

Initial environmental pyrosequencing studies suggested highly complex protistan communities with phylotype richness decisively higher than previously estimated. However, recent studies on individual bacteria or artificial bacterial communities evidenced that pyrosequencing errors may skew our view of the true complexity of microbial communities. We pyrosequenced two diversity markers (hypervariable regions V4 and V9 of the small-subunit rDNA) of an intertidal protistan model community, using the Roche GS-FLX and the most recent GS-FLX Titanium sequencing systems. After pyrosequencing 24 reference sequences we obtained up to 2039 unique tags (from 3879 V4 GS-FLX Titanium reads), 77% of which were singletons. Even binning sequences that share 97% similarity still emulated a pseudodiversity exceeding the true complexity of the model community up to three times (V9 GS-FLX). Pyrosequencing error rates were higher for V4 fragments compared with the V9 domain and for the GS-FLX Titanium compared with the GS-FLX system. Furthermore, this experiment revealed that error rates are taxon-specific. As an outcome of this study we suggest a fast and efficient strategy to discriminate pyrosequencing signals from noise in order to more realistically depict the structure of protistan communities using simple tools that are implemented in standard tag data-processing pipelines.     

Performance comparison of benchtop high-throughput sequencing platforms

Three benchtop high-throughput sequencing instruments are now available. The 454 GS Junior (Roche), MiSeq (Illumina) and Ion Torrent PGM (Life Technologies) are laser-printer sized and offer modest set-up and running costs. Each instrument can generate data required for a draft bacterial genome sequence in days, making them attractive for identifying and characterizing pathogens in the clinical setting. We compared the performance of these instruments by sequencing an isolate of Escherichia coli O104:H4, which caused an outbreak of food poisoning in Germany in 2011. The MiSeq had the highest throughput per run (1.6 Gb/run, 60 Mb/h) and lowest error rates. The 454 GS Junior generated the longest reads (up to 600 bases) and most contiguous assemblies but had the lowest throughput (70 Mb/run, 9 Mb/h). Run in 100-bp mode, the Ion Torrent PGM had the highest throughput (80–100 Mb/h). Unlike the MiSeq, the Ion Torrent PGM and 454 GS Junior both produced homopolymer-associated indel errors (1.5 and 0.38 errors per 100 bases, respectively).     

Soil eukaryotic microorganism succession as affected by continuous cropping of peanut--pathogenic and beneficial fungi were selected

Peanut is an important oil crop worldwide and shows considerable adaptability but growth and yield are negatively affected by continuous cropping. Soil micro-organisms are efficient bio-indicators of soil quality and plant health and are critical to the sustainability of soil-based ecosystem function and to successful plant growth. In this study, 18S rRNA gene clone library analyses were employed to study the succession progress of soil eukaryotic micro-organisms under continuous peanut cultivation. Eight libraries were constructed for peanut over three continuous cropping cycles and its representative growth stages. Cluster analyses indicated that soil micro-eukaryotic assemblages obtained from the same peanut cropping cycle were similar, regardless of growth period. Six eukaryotic groups were found and fungi predominated in all libraries. The fungal populations showed significant dynamic change and overall diversity increased over time under continuous peanut cropping. The abundance and/or diversity of clones affiliated with Eurotiales, Hypocreales, Glomerales, Orbiliales, Mucorales and Tremellales showed an increasing trend with continuous cropping but clones affiliated with Agaricales, Cantharellales, Pezizales and Pyxidiophorales decreased in abundance and/or diversity over time. The current data, along with data from previous studies, demonstrated that the soil microbial community was affected by continuous cropping, in particular, the pathogenic and beneficial fungi that were positively selected over time, which is commonplace in agro-ecosystems. The trend towards an increase in fungal pathogens and simplification of the beneficial fungal community could be important factors contributing to the decline in peanut growth and yield over many years of continuous cropping.     

Effects of Polyethylene Glycol 400 Monolaurate and Sodium Chloride Mixtures on the Antibacterial Activity of Solubilized Trichlorocarbanilide

S ummary . Polyethylene glycol 400 monolaurate (PGM), a nonionic surfactant, affected the solubility of trichlorocarbanilide (TCC), the amount solubilized being directly proportional to the concentration of PGM in the system. The bacteriostatic activity of TCC against Staphylococcus aureus was affected by high but not low concentrations of PGM; proportions of PGM to the germicide of up to 10: 1 did not interfere with the minimum inhibitory concentration (MIC) of TCC (0·08 μg/ml). The bactericidal activity against Staph. aureus of systems consisting of TCC and different proportions of PGM with and without 0·5% of NaCl was much higher than that of saturated aqueous solutions of TCC alone. However, TCC solutions containing NaCl and PGM were more bactericidal than corresponding systems without NaCl; the death rates in TCC-PGM-NaCl systems were 2–4 times that in TCC-PGM systems. The marked increase in the death rate was obtained with systems containing 2·5% of PGM.     

冬季覆盖作物秸秆还田对双季稻田根际土壤微生物群落功能多样性的影响

微生物群落功能多样性是土壤质量变化重要的指标,不同作物类型的秸秆还田措施对土壤微生物群落功能多样性具有明显的影响。以位于双季稻主产区不同冬季覆盖作物-双季稻种植模式大田定位试验田为研究对象,以冬闲-双季稻种植模式为对照(CK),应用Biolog-GN技术开展黑麦草-双季稻(Ry)、紫云英-双季稻(Mv)、油菜-双季稻(Ra)和马铃薯-双季稻(Po)种植模式条件下不同冬季覆盖作物秸秆还田后对双季稻田根际土壤微生物功能多样性影响的研究。研究结果表明,早稻和晚稻成熟期,与CK处理相比,冬季覆盖作物秸秆还田处理增加了稻田土壤碳源平均颜色变化率(AWCD),以Po处理AWCD均为最高,均显著高于Ry和CK处理。不同冬季覆盖作物秸秆还田处理土壤微生物代谢多样性指数表现出明显的差异,早稻成熟期,Po处理的Richness、Shannon和McIntosh指数均为最高,其次为Ry、Mv和Ra处理,CK处理最低;晚稻成熟期,各处理的Richness、Shannon和McIntosh指数大小顺序均表现为PoRaMvRyCK。土壤微生物碳源利用的主成分分析结果表明,各冬季覆盖作物秸秆还田处理根际土壤微生物利用的主要碳源为氨基酸类和糖类物质,不同处理间碳源利用类型有差异。冬季覆盖作物秸秆还田措施有利于提高双季稻田根际土壤微生物对碳源的利用能力、物种丰富度和均匀度。     

Groundtruthing next-gen sequencing for microbial ecology-biases and errors in community structure estimates from PCR amplicon pyrosequencing

Analysis of microbial communities by high-throughput pyrosequencing of SSU rRNA gene PCR amplicons has transformed microbial ecology research and led to the observation that many communities contain a diverse assortment of rare taxa-a phenomenon termed the Rare Biosphere. Multiple studies have investigated the effect of pyrosequencing read quality on operational taxonomic unit (OTU) richness for contrived communities, yet there is limited information on the fidelity of community structure estimates obtained through this approach. Given that PCR biases are widely recognized, and further unknown biases may arise from the sequencing process itself, a priori assumptions about the neutrality of the data generation process are at best unvalidated. Furthermore, post-sequencing quality control algorithms have not been explicitly evaluated for the accuracy of recovered representative sequences and its impact on downstream analyses, reducing useful discussion on pyrosequencing reads to their diversity and abundances. Here we report on community structures and sequences recovered for in vitro-simulated communities consisting of twenty 16S rRNA gene clones tiered at known proportions. PCR amplicon libraries of the V3-V4 and V6 hypervariable regions from the in vitro-simulated communities were sequenced using the Roche 454 GS FLX Titanium platform. Commonly used quality control protocols resulted in the formation of OTUs with >1% abundance composed entirely of erroneous sequences, while over-aggressive clustering approaches obfuscated real, expected OTUs. The pyrosequencing process itself did not appear to impose significant biases on overall community structure estimates, although the detection limit for rare taxa may be affected by PCR amplicon size and quality control approach employed. Meanwhile, PCR biases associated with the initial amplicon generation may impose greater distortions in the observed community structure.     

Hot-Alkaline DNA Extraction Method for Deep-Subseafloor Archaeal Communities

A prerequisite for DNA-based microbial community analysis is even and effective cell disruption for DNA extraction. With a commonly used DNA extraction kit, roughly two-thirds of subseafloor sediment microbial cells remain intact on average (i.e., the cells are not disrupted), indicating that microbial community analyses may be biased at the DNA extraction step, prior to subsequent molecular analyses. To address this issue, we standardized a new DNA extraction method using alkaline treatment and heating. Upon treatment with 1 M NaOH at 98°C for 20 min, over 98% of microbial cells in subseafloor sediment samples collected at different depths were disrupted. However, DNA integrity tests showed that such strong alkaline and heat treatment also cleaved DNA molecules into short fragments that could not be amplified by PCR. Subsequently, we optimized the alkaline and temperature conditions to minimize DNA fragmentation and retain high cell disruption efficiency. The best conditions produced a cell disruption rate of 50 to 80% in subseafloor sediment samples from various depths and retained sufficient DNA integrity for amplification of the complete 16S rRNA gene (i.e., ∼1,500 bp). The optimized method also yielded higher DNA concentrations in all samples tested compared with extractions using a conventional kit-based approach. Comparative molecular analysis using real-time PCR and pyrosequencing of bacterial and archaeal 16S rRNA genes showed that the new method produced an increase in archaeal DNA and its diversity, suggesting that it provides better analytical coverage of subseafloor microbial communities than conventional methods.     

Modelling herbivore movement and colonization: pest management potential of intercropping and trap cropping

1 Using a stochastic simulation model, we explored the effects of agroecosystem diversity on herbivore densities. 2 Using parameters that included reproduction, colonization, and local movement rates, we simulated an insect herbivore population colonizing rows of plants in an agricultural setting. 3 Plant rows were comprised of either principal crop, intercrop, or trap crop. Herbivore parameters varied for different plant types. 4 Percent crop cover and movement rates were varied, and ensuing herbivore densities on crop rows were recorded. 5 In trap cropping schemes, both percent crop cover and movement rates were critical in determining herbivore densities. Intercropping schemes were governed primarily by colonization rates. 6 These results suggest that trap cropping schemes merit more attention than intercropping systems in the design and analysis of mixed cropping systems.     

Diversity and dynamics of the DNA and cDNA-derived bacterial compost communities throughout the <Emphasis Type="Italic">Agaricus bisporus</Emphasis> mushroom cropping process

The cultivation of Agaricus bisporus involves the conversion of agricultural materials via fermentation into utilisable simple sugars as a nutrient source for the fungal crop during mushroom cropping. However, little is currently known about the role of the bacterial community contributing to the fermentation process. In this investigation we characterised the composition and dynamics of the DNA and cDNA-derived bacterial populations throughout a commercial mushroom cropping process using MiSeq sequencing. Both methods indicated substantial changes in the bacterial community structure after the first flush of the mushroom crop. However, differences were observed between the composition of the bacterial community determined by each of the two methods. The cDNA-derived community indicated that thermotolerant genera with known sulphur-reducing characteristics were highly active up to the first flush. Activity of the phyla Actinobacteria and Firmicutes was observed to increase as fermentation progressed, indicating that the members of these phyla played prominent roles in the conversion of wheat straw into utilisable sugars. The cDNA-derived community comprised genera with roles in the nitrification process that became highly active at post flush 1. Subsequent chemical analysis of extractable nitrate indicated that substantial nitrification occurred up until the termination of the cropping process. This study has demonstrated that a highly dynamic bacterial community is present throughout the mushroom cropping process.     

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.     

稻田轮作系统的生态学分析

针对目前稻田单一化的连作耕作制度对农田环境造成的不利因素,通过5a田间定位试验,从土壤理化性状、作物产量变化、病虫害发生发展规律、能流及养分平衡状况等角度对稻田连作耕作制度和轮作系统进行生态学分析.结果表明,与连作耕作制度相比较,在试验设计范围内,稻田轮作系统明显改善了土壤的理化性状,使得土壤随着耕种年限增加,容重下降,而孔隙度增加,固相比率下降,气相比率上升,气液比值增大,土壤通透性大大增强,有效阻止土壤次生潜育化和土壤酸化,提高土壤pH值.轮作系统不但提高了作物产量,而且总初级生产力、光能利用率、辅助能利用率分别比连作系统高17.47%、9.87%和5.0%.N、P、K的养分利用率也同样明显高于连作系统.提出了合理轮作的优化模式和复种模式.     

江苏省棉区棉田杂草群落发生分布规律的数量分析

在对江苏省四大主产棉区８２个样点７２７块田地共５７．５ｈｍ＾２棉田杂草群落及草害进行７级目测法调查采集数据后,对其进行主成分分析（ＰＣＡ）,并赋以生态学意义的解释。研究表明,江苏省棉田杂草的发生和分布与轮作种植制度和地理区域性密切相关,其中轮作制度导致的田间水分的巨大差异是决定杂草群落结构特征的最深刻的原因,导致全省水旱轮作棉田的杂草群落有趋向性。而地理区域构成的土壤、气候等生态因子的显著影响表现     

Assessment of metagenomic assembly using simulated next generation sequencing data

Due to the complexity of the protocols and a limited knowledge of the nature of microbial communities, simulating metagenomic sequences plays an important role in testing the performance of existing tools and data analysis methods with metagenomic data. We developed metagenomic read simulators with platform-specific (Sanger, pyrosequencing, Illumina) base-error models, and simulated metagenomes of differing community complexities. We first evaluated the effect of rigorous quality control on Illumina data. Although quality filtering removed a large proportion of the data, it greatly improved the accuracy and contig lengths of resulting assemblies. We then compared the quality-trimmed Illumina assemblies to those from Sanger and pyrosequencing. For the simple community (10 genomes) all sequencing technologies assembled a similar amount and accurately represented the expected functional composition. For the more complex community (100 genomes) Illumina produced the best assemblies and more correctly resembled the expected functional composition. For the most complex community (400 genomes) there was very little assembly of reads from any sequencing technology. However, due to the longer read length the Sanger reads still represented the overall functional composition reasonably well. We further examined the effect of scaffolding of contigs using paired-end Illumina reads. It dramatically increased contig lengths of the simple community and yielded minor improvements to the more complex communities. Although the increase in contig length was accompanied by increased chimericity, it resulted in more complete genes and a better characterization of the functional repertoire. The metagenomic simulators developed for this research are freely available.     

江苏省棉区棉田杂草群落发生分布规律的数量分析

在对江苏省四大主产棉区８２个样点７２７块田地共５７．５ｈｍ＾２棉田杂草群落及草害进行了７级目测法调查采集数据后,对其进行主成分分析,并赋以生态学意义的解释,研究表明,江苏省棉田杂草的发生和分布与轮作种植制度和区域性化作制度导致的田间水分的巨大差异是决定杂草群落结构特征的最深刻的原因,导致全省水旱轮作棉 杂草群落有趋同性。