PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of ismejLink to Publisher's site
 
ISME J. Aug 2012; 6(8): 1621–1624.
Published online Mar 8, 2012. doi:  10.1038/ismej.2012.8
PMCID: PMC3400413
Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms
J Gregory Caporaso,1 Christian L Lauber,2 William A Walters,3 Donna Berg-Lyons,2 James Huntley,4 Noah Fierer,2,5 Sarah M Owens,6 Jason Betley,7 Louise Fraser,7 Markus Bauer,7 Niall Gormley,7 Jack A Gilbert,6,8 Geoff Smith,7 and Rob Knight9,10*
1Department of Computer Science, Northern Arizona University, Flagstaff, AZ, USA
2Cooperative Institute for Research in Environmental Sciences, UCB 216, University of Colorado, Boulder, CO, USA
3Department of Molecular, Cellular and Developmental Biology, UCB 347, University of Colorado, Boulder, CO, USA
4Colorado Initiative in Molecular Biotechnology, UCB 347, University of Colorado, Boulder, CO, USA
5Department of Ecology and Evolutionary Biology, UCB 334, University of Colorado, Boulder, Colorado, USA
6Argonne National Laboratory, Argonne, IL, USA
7Illumina Cambridge Ltd., Chesterford Research Park, Saffron Walden, Essex, UK
8Department of Ecology and Evolution, University of Chicago, Chicago, IL, USA
9Department of Chemistry and Biochemistry, UCB 215, University of Colorado, Boulder, CO, USA
10Howard Hughes Medical Institute, University of Colorado at Boulder, UCB 215, Boulder, CO, USA
*Howard Hughes Medical Institute, University of Colorado at Boulder, UCB 215, Boulder, CO 80309, USA. E-mail: rob/at/spot.colorado.edu
Received September 12, 2011; Revised January 13, 2012; Accepted January 19, 2012.
DNA sequencing continues to decrease in cost with the Illumina HiSeq2000 generating up to 600 Gb of paired-end 100 base reads in a ten-day run. Here we present a protocol for community amplicon sequencing on the HiSeq2000 and MiSeq Illumina platforms, and apply that protocol to sequence 24 microbial communities from host-associated and free-living environments. A critical question as more sequencing platforms become available is whether biological conclusions derived on one platform are consistent with what would be derived on a different platform. We show that the protocol developed for these instruments successfully recaptures known biological results, and additionally that biological conclusions are consistent across sequencing platforms (the HiSeq2000 versus the MiSeq) and across the sequenced regions of amplicons.
Keywords: illumine, barcoded sequencing, QIIME
DNA sequencing cost continues to decline: a vast price per sequence decrease on Illumina HiSeq2000 and MiSeq platforms further supports democratization of sequencing (Tringe and Hugenholtz, 2008). Interest in amplicon sequencing on Illumina is growing (Bartram et al., 2011; Caporaso et al., 2011; Zhou et al., 2011), largely due to lower cost per sequence than other platforms, enabling high-throughput microbial ecology at the greatest coverage yet possible. Although some technical issues exist with community sequencing, such as PCR primer biases and differential DNA extraction efficiency from different organisms in complex communities, these techniques continue to vastly expand our understanding of the microbial world.
Here we present an amplicon sequencing protocol for the HiSeq2000 and MiSeq platforms, and apply this protocol to sequence host-associated and free-living microbial communities to verify that biological conclusions drawn from the data are consistent across platforms and sequence reads. The HiSeq and MiSeq platforms differ markedly in scale. The HiSeq2000 produces >50 Gb per day, and in the course of a 10.8 day run produces 1.6 billion 100-base paired-end reads. By contrast, the MiSeq is for single-day experiments, and generates 1.5 Gb per day from 5 million 150-base paired-end reads. Our results capture known differences between microbial communities on each platform; biological conclusions drawn are consistent across platforms and sequence reads. This protocol is therefore ready for widespread use in microbial community analysis, such as by the Earth Microbiome Project (Gilbert et al., 2010), which has adopted it for amplicon sequencing. Details on the sequencing protocol are provided as Supplementary Methods.
Twenty-four samples were sequenced on three paired-end Illumina HiSeq2000 lanes, and in one paired-end MiSeq run. The samples represented soil (source: USA; n=8) and several host-associated environment types: human feces (source: USA; n=2), mouth (source: USA; n=2) and skin (source: USA; n=6); canine feces (source: USA; n=1) mouth (source: USA; n=1) and skin (source: USA; n=4). These four paired-end lanes (three on HiSeq and one on MiSeq) resulted in eight sets of reads, corresponding to 5′ and 3′ reads from each lane. These sets of reads were treated as independent replicates to assess the reproducibility of the results.
We were primarily interested in whether known differences between microbial communities could be recaptured on these Illumina platforms to determine their suitability for large-scale surveys of microbial communities. We observed several expected results in principal coordinates plots of weighted UniFrac distances (Figure 1). First, we observed primary separation of samples based on whether they were derived from a free-living environment (soil; cyan) or host-associated environment (all other colors) (Ley et al., 2008). Next we observed separation of fecal samples (yellow; red) from all other host-associated sample types (Costello et al., 2009).
Figure 1
Figure 1
Procrustes plots comparing: (a) 5′ reads from HiSeq lane 6 to 5′ reads from HiSeq lane 8; (b) 5′ reads from HiSeq lane 6 to 3′ reads from HiSeq lane 8; (c) 5′ reads from HiSeq lane 6 to 5′ MiSeq reads; ( (more ...)
We were additionally interested in reproducibility across lanes and reads within and between each platform. To test this, we ran the 24 samples on three HiSeq paired-end lanes and 1 MiSeq paired-end lane, and analyzed each resulting set of reads independently. As our biological conclusions are frequently driven by the results of principal coordinates analyses based on weighted UniFrac distances, we compared these plots using Procrustes analysis (Gower, 1975; Figure 1; Table 1) as implemented in QIIME and found that the observations were highly reproducible across lanes, read directions and platforms. All 28 possible lane/read pair combinations produced highly significant P-values based on 10 000 Monte Carlo iterations (P<0.0001; Bonferroni-adjusted α0.01=0.0004).
Table 1
Table 1
M2 and Monte Carlo P-values for all Procrustes comparisons
Taken together, these results suggest that the protocol previously developed for high-throughput community sequencing on the Illumina GAIIx has been successfully adapted for the HiSeq2000 and MiSeq platforms, again greatly decreasing the cost per sequence of amplicon sequencing to ~15 000 single-end reads per USD$1 on the HiSeq2000. For example, based on our lowest high-quality sequence per lane count of 22 928 291 reads (Supplementary File 2, HiSeq 3′ lane 6), if using all 2167 barcodes in each of 15 lanes on the HiSeq2000, leaving one lane for a control, then it is possible to sequence 32 505 samples in a week at a depth of 10 580 sequences per sample for approximately $22 000 in sequencing costs. Longer barcodes could additionally be developed to facilitate more sequences per sample at a lower depth of sequencing. On the basis of the lowest high-quality sequence count on the MiSeq of 1 603 532 reads (Supplementary File 2, MiSeq 3′), if using all 2167 barcodes, it is possible to sequence 2167 samples in a 12 h run at a depth of 740 sequence per sample for approximately $800 in sequencing costs.
A relevant question is whether the decreased cost of sequencing should be applied to obtain deeper coverage of samples, or to increase the number of samples that are sequenced. Figure 1c compares the results of sequencing the same samples on the HiSeq 2000 at a median depth of 1207 709 sequences per sample and the MiSeq platform at a depth of 43 271 sequences per sample. The highly significant Procrustes result (P<0.0001) implies that we draw the same beta diversity conclusions from either sequencing run, despite a two order of magnitude increase in sequencing depth on the HiSeq2000. Similarly, when sampling to only 10 sequences per sample Procrustes results are still highly significant (P<0.0001; Supplementary Figure 1), although the higher M2 value indicates that the correlation is not as strong as when sampling to 100 sequences per sample. These observations, in agreement with studies that have addressed this question directly (Kuczynski et al., 2010), suggest that increasing the sequencing depth is not likely to provide additional insight into questions of beta diversity, and we therefore argue that (for questions of beta diversity in particular) the decreased cost of sequencing should be applied to study microbial systems using many more samples, for example, in dense temporal or spatial analyses, rather than with many more sequences per sample. Of course, if the objective is to identify taxa that are very rare in communities, deeper sequencing will be advantageous. Additionally we note that while as few as 10 sequences per sample may be useful for differentiating very different environment types (for example, soil and feces), as environments become more similar (for example, two soil samples of different pH) more sequences will be required to differentiate them.
As sequencing costs continue to decrease our studies of the microbial world can continue to increase in scope. The protocol presented here opens the HiSeq2000 and MiSeq Illumina platforms to community amplicon sequencing. The data generated by each is similar, but differs in scale and therefore support different applications. For large projects where time is less of an issue but cost per sequence is a major concern, the HiSeq platform allows massively parallel sequencing at the lowest cost. Here we show that comparable data can be generated on the MiSeq for smaller projects where it is important to process samples quickly, for example, in routine environmental or patient monitoring or in preliminary investigations for larger projects. We expect that this is another step toward the era of ubiquitous DNA sequencing, when sequencers become standard equipment in research and clinical laboratories. Finally, we show that technical replicates run on different sequencing platforms and from sequencing of different regions of amplicons should yield the same biological conclusions: critical information as more sequencing platforms become available.
Acknowledgments
We wish to thank the National Ecological Observatory Network (a project sponsored by the National Science Foundation and managed under cooperative agreement by NEON, Inc.) for donation of the soil samples; and Aurelie Breton and Joshua Quick for running the MiSeq instrument. This work was funded in part by Amazon Web Services, NIH, Crohn's and Colitis Foundation of America, The Bill and Melinda Gates Foundation, and the Howard Hughes Medical Institute.
Notes
Several authors on this manuscript are employees of Illumina, Inc., whose technology is tested in this study.
Footnotes
Supplementary Information accompanies the paper on The ISME Journal website (http://www.nature.com/ismej)
Supplementary Material
Supplementary Figure 1
Supplementary File 1
Supplementary File 2
Supplementary Information
Supplementary Methods
  • Bartram AK, Lynch MD, Stearns JC, Moreno-Hagelsieb G, Neufeld JD. Generation of multimillion-sequence 16S rRNA gene libraries from complex microbial communities by assembling paired-end illumina reads. Appl Environ Microbiol. 2011;77:3846–3852. [PMC free article] [PubMed]
  • Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Lozupone CA, Turnbaugh PJ, et al. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc Natl Acad Sci USA. 2011;108 (Suppl 1:4516–4522. [PubMed]
  • Costello EK, Lauber CL, Hamady M, Fierer N, Gordon JI, Knight R. Bacterial community variation in human body habitats across space and time. Science. 2009;326:1694–1697. [PubMed]
  • Gilbert JA, Meyer F, Jansson J, Gordon J, Pace N, Tiedje J, et al. The Earth Microbiome Project: meeting report of the ‘1 EMP meeting on sample selection and acquisition' at Argonne National Laboratory October 6 2010. Stand Genomic Sci. 2010;3:249–253. [PMC free article] [PubMed]
  • Gower JC. Generalized procrustes analysis. Psychometrika. 1975;40:33–51.
  • Kuczynski J, Liu Z, Lozupone C, McDonald D, Fierer N, Knight R. Microbial community resemblance methods differ in their ability to detect biologically relevant patterns. Nat Methods. 2010;7:813–819. [PMC free article] [PubMed]
  • Ley RE, Lozupone CA, Hamady M, Knight R, Gordon JI. Worlds within worlds: evolution of the vertebrate gut microbiota. Nat Rev Microbiol. 2008;6:776–788. [PMC free article] [PubMed]
  • Tringe SG, Hugenholtz P. A renaissance for the pioneering 16S rRNA gene. Curr Opin Microbiol. 2008;11:442–446. [PubMed]
  • Zhou HW, Li DF, Tam NF, Jiang XT, Zhang H, Sheng HF, et al. BIPES, a cost-effective high-throughput method for assessing microbial diversity. ISME J. 2011;5:741–749. [PMC free article] [PubMed]
Articles from The ISME Journal are provided here courtesy of
Nature Publishing Group