A framework for human microbiome research.

A variety of microbial communities and their genes (the microbiome) exist throughout the human body, with fundamental roles in human health and disease. The National Institutes of Health (NIH)-funded Human Microbiome Project Consortium has established a population-scale framework to develop metagenomic protocols, resulting in a broad range of quality-controlled resources and data including standardized methods for creating, processing and interpreting distinct types of high-throughput metagenomic data available to the scientific community. Here we present resources from a population of 242 healthy adults sampled at 15 or 18 body sites up to three times, which have generated 5,177 microbial taxonomic profiles from 16S ribosomal RNA genes and over 3.5 terabases of metagenomic sequence so far. In parallel, approximately 800 reference strains isolated from the human body have been sequenced. Collectively, these data represent the largest resource describing the abundance and variety of the human microbiome, while providing a framework for current and future studies.

Introduction

Advances in sequencing technologies coupled with novel bioinformatic developments have allowed the scientific community to commence assessment of the uncultivated majority; the microbes that inhabit our oceans, soils, human body and other locations [1]. Microbes associated with the human body include eukaryotes, archaea, bacteria, and viruses, with bacteria alone estimated to outnumber human cells within an individual by an order of magnitude. Our knowledge of these communities and their gene content, referred to collectively as the human microbiome, has to date been limited by a lack of population-scale data detailing their composition and function.

The US National Institutes of Health (NIH) funded Human Microbiome Project (HMP) Consortium allied a broad collection of scientific experts to explore these microbial communities and their relationships with their human hosts. As such, the HMP [2] has focused on reference genomes (viral, bacterial and eukaryotic) which provide a critical framework for subsequent metagenomic annotation and analysis, and on generating a baseline of microbial community structure and function from an adult cohort defined by a carefully delineated set of clinical inclusion and exclusion criteria which we term “healthy” in this study (http://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/GetPdf.cgi?id=phd002854.2). Investigations of the microbiome from this cohort incorporated several complementary analyses including: 16S rRNA gene sequence (16S) and taxonomic profiles, whole genome shogun (WGS) or metagenomic sequencing of whole community DNA, and alignment of the assembled sequences to the reference microbial genomes from the human body [3, 4]. Thus, the HMP complements other large-scale sequence-based human microbiome projects such as the MetaHIT project [5] which focused on examination of the gut microbiome using WGS data including samples from cohorts exhibiting a wide range of health status and physiological characteristics.

Additional projects supported by the HMP are investigating the association of specific components and dynamics of the microbiome with a variety of disease conditions, developing tools and technology including isolating and sequencing uncultured organisms, and studying the ethical, legal, and social implications of human microbiome research (http://commonfund.nih.gov/hmp/fundedresearch.aspx). A comprehensive list of current publications from HMP projects is available at http://commonfund.nih.gov/hmp/publications.aspx.

Here we detail the resources created to date by the HMP initiative including: clinical specimens (samples), reference genomes, sequencing and annotation protocols, methods and analyses. We describe the thousands of samples obtained from 15 to 18 distinct body sites from 242 donors over multiple time points that were processed at two clinical centers (Baylor College of Medicine (BCM) and Washington University School of Medicine). We also describe the laboratory and computational protocols developed for reliably generating and interpreting the human microbiome data. HMP resources include both protocols for, and the subsequent data generated from, 16S and metagenomic sequencing of human microbiome samples. During this study, these protocols were rigorously standardized and quality controlled for simultaneous use across four sequencing centers (BCM Human Genome Sequencing Center, The Broad Institute of Harvard and MIT, J. Craig Venter Institute, The Genome Institute at Washington University). In particular, we focus on the production of the first phase of metagenomic data sets (Phase I) used for subsequent in depth analyses and we summarize standards and recommendations based on our experiences generating and analyzing these data. An additional set of publications (many included in this communication’s references and in [4]) will describe in further detail the microbial ecology and microbiological implications of these data. Collectively these resources and analyses represent an important framework for human microbiome research.

Main

HMP resource organization

Supplementary Figure 1 summarizes organization of the HMP, including the data processing and analytical steps, and the scientific entities gathered to conduct the project. An overview of available HMP data sets and additional resources are provided in Supplementary Tables 1, 2 and 3. Donors were recruited and enrolled into the HMP through the two clinical centers. Over 240 adults were carefully screened and phenotyped prior to sampling one to three times at 15 (male) or 18 (female) body sites using a common sampling protocol (http://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/GetPdf.cgi?id=phd003190.2). All included subjects were between the ages of 18 and 40 years and had passed a screening for systemic health based on oral, cutaneous, and body mass exclusion criteria (http://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/GetPdf.cgi?id=phd002854.2) [6].

A Data Analysis and Coordination Center (DACC) was created to serve as the central repository for all HMP WGS, 16S, and reference genome sequence information generated by the four sequencing centers. The DACC supports access to analysis software, biological samples, clinical protocols, news, publication announcements and project statistics, and performed centralized analysis of HMP reference genome and WGS annotation in cooperation with the sequencing centers. Unless otherwise noted, all data sets and protocols described here are available to the scientific community at the DACC (http://hmpdacc.org).

Phase I 16S and WGS sequencing overview

A set of 5,298 samples were collected from 242 adults [6] (Table 1, Supplementary Table 4) from which 16S and WGS data were generated for a total of 5,177 taxonomically characterized communities (16S) and 681 WGS samples describing the microbial communities from habitats within the human airways, skin, oral cavity, gut, and vagina. For a subset of 560 samples, both data types were generated (Table 1). These efforts constitute our initial primary metagenomic data sets (Phase I) described in more detail below. Additional efforts are ongoing to sequence and analyze the remaining samples from the complete HMP collection (11,174 primary specimens in total from 300 individuals sampled up to three times over 22 months) [6].

Table 1

HMP donor samples examined by 16S and WGS

Body Region	Body Site	Total Samples	Total 16S Samples	V13 Samples	V13 Read Depth (M)*	V35 Samples	V35 Read depth (M)*	Samples V13&V35	Total WGS Samples	Total Read Depth (G)†	% Filtered Reads±	% Human Reads¥	Remaining Read Depth (G)†	Samples 16S&WGS
Gut	Stool	352	337	193	1.4	328	2.38	184	139	1720.7	15	1	1450.6	124
	Buccal mucosa	346	330	184	1.3	314	1.74	168	107	1438	9	82	136.68	91
	Hard palate	325	325	179	1.2	310	1.67	164	1	10.87	20	25	5.92	1
	Keratinized gingiva	335	329	183	1.3	319	1.74	173	6	72.3	5	47	34.42	0
	Palatine Tonsils	337	332	189	1.2	315	1.87	172	6	74.75	2	80	13.45	1
Oral Cavity	Saliva	315	310	166	0.9	292	1.45	148	5	55.69	1	91	4.24	0
	Subgingival plaque	334	328	186	1.2	314	1.84	172	7	92.06	5	79	15.29	1
	Supragingival plaque	345	331	192	1.3	316	1.88	177	115	1500.7	15	40	674.81	101
	Throat	331	325	176	1	312	1.67	163	7	78.78	4	79	13.57	1
	Tongue dorsum	348	332	193	1.3	320	2.04	181	122	1620.1	15	19	1084.3	106
Airway	Anterior nares	316	302	169	1	283	1.17	150	84	1129.9	3	96	14.31	70
	Left Antecubital fossa	269	269	158	0.7	221	0.47	110	0	na	na	na	0	na
	Left Retroauricular crease	313	312	188	1.6	295	1.46	171	9	126.34	9	73	22.07	8
Skin	Right Antecubital fossa	274	274	158	0.7	229	0.52	113	0	na	na	na	0	na
	Right Retroauricular crease	319	316	190	1.4	304	1.56	178	15	181.94	18	59	42.38	12
	Mid vagina	145	143	91	0.6	140	0.96	88	2	22.58	0	99	0.18	0
Vagina	Posterior fornix	152	142	89	0.6	136	0.98	83	53	702.13	6	90	25.24	43
	Vaginal introitus	142	140	87	0.6	131	0.85	78	3	36.48	1	98	0.58	1
total		5298	5177	2971	19	4879	26.3	2673	681	8863.3	11	49	3538.1	560

1×106 reads post-processing with the mothur pipeline (Supplementary Information)

1×109 reads post-processing with the mothur pipeline (Supplementary Information)

Fraction of reads with low quality bases that were removed (Supplementary Information)

Fraction of human reads that were removed (Supplementary Information)

16S standards development and sequencing

The goals of the HMP required that 16S sequences and profiles from data produced at the four participating sequencing centers be comparable in a variety of downstream analyses however, no suitable methodology was available at the commencement of the project. While establishing 16S protocols, we determined that many components of data production and processing can contribute errors and artifacts. We investigated methods that avert these errors and their subsequent effects on taxonomic classification and OTU-based community structure. The results are discussed in detail in Supplementary Information and [7]. Thus, multiple evaluations of 16S protocols were undertaken before adopting a single standardized protocol that ensured consistency in the high throughput production.

To maximize accuracy and consistency, protocols were evaluated primarily using a synthetic mock community (MC) of 21 known organisms (Supplementary Table 5, [7]). Additional testing of the protocol was carried out on a subset of HMP samples (Supplementary Table 1). Collectively, these efforts resulted in adoption of a protocol to amplify and sequence samples using the Roche-454 FLX Titanium platform ([7], http://www.hmpdacc.org/doc/HMP_MDG_454_16S_Protocol_V4_2_102109.pdf). The HMP created both cell mixtures and genomic DNA extracts of the MC (Supplementary Table 3). A large body of metagenomic data (both 16S and WGS) (RES:HMMC) from these and other calibration experiments are available to the community to facilitate further benchmarking of new molecular and analytical approaches (Supplementary Table 2).

The majority of the sample collection was targeted for 16S sequencing using the 454 FLX Titanium based strategy [7]. The nucleotide sequence of the 16S rRNA gene consists of regions of highly conserved sequence which alternate with nine regions or windows of variable nucleotide sequence that constitute the most informative portions of the gene sequence for use in taxonomic classification. A window covering the number three (V3) through five (V5) variable regions (V35) of the 16S rRNA gene was chosen as the target for 4,879 samples. Sequence of a V1 to V3 (V13) window was also included for a subset of 2,971 samples to provide a complementary view of taxonomic profiles (RES:HMR16S, Supplementary Figs. 2 and 3, Supplementary Information, [7]).

Following adoption of the 16S protocol including removal of multiple sources of potential artifacts or bias generated by 16S sequencing using pyrosequencing [8, 9], a variety of approaches for accurate diversity estimation were developed and compared [10]. A 16S data processing pipeline was established using the mothur software package [11] (Supplementary Information) that includes two optional low and high stringency approaches. The former provides an output favoring longer read lengths tailored towards taxonomic classification, the latter an output with more aggressive sequence error reduction tailored towards Operational Taxonomic Unit (OTU) construction (RES:HMMCP). A third complementary pipeline was also developed using the QIIME software package [12] (Supplementary Information), which processes these data using an OTU-binning strategy to which taxonomic classification is added (RES:HMQCP). All pipelines result in highly comparable views of the human microbiome.

Metagenomic assembly and gene cataloging

Approximately 749 samples representing targeted body sites were chosen for WGS sequencing using the Illumina GAIIx platform with 101 bp paired-end reads. From a high quality set of 681 samples an average depth of 13 Gb (±4.3) was achieved per sample, collectively producing a total of 8.8 Tb (RES:HMIWGS) (Table 1). Theoretically, these per sample data are sufficient to cover a 3Mb bacterial genome present at only 0.8% abundance with a probability of 90% [13]. In addition, 12 stool samples were simultaneously sequenced using the 454 FLX Titanium platform (RES:HM4WGS). Comparisons between the centers demonstrated high consistency of target sequencing depth and success rates [4]. Following development of a protocol for removing reads resulting from human DNA contamination (Supplementary Information), 49% of the reads were targeted for removal as human (for information on authorized access to these reads see Supplementary Information). Samples collected from soft tissue tended to have higher human contamination (e.g., mid vagina (96%), anterior nares (82%), throat (75%)). Preparations from saliva were also high in human DNA sequence (80%), while stool contained a relatively low abundance of human reads (up to 1%) (Supplementary Fig. 4).

After application of a quality control protocol that includes human sequence removal, quality filtering and trimming of reads, (Supplementary Information), the remaining 3.5 Tb from 681 samples were subjected to a three-tiered complementary analysis strategy (Supplementary Information), of reference genome mapping (that was able to use ~57% of the data), assembly and gene prediction (~50% of the data), and metabolic reconstruction (~36% of the data). This combined strategy facilitated the extraction of maximal organismal and functional information.

Metagenomic assemblies were generated for all available samples using an optimized SOAPdenovo protocol with parameters designed to produce substrates for downstream analyses such as gene and function prediction, resulting in a total of 41 million contigs (RES:HMASM) (Supplementary Information). Reads that remained unassembled were pooled across individual body sites and re-assembled using the same approach, resulting in an additional 4,200,672 contigs (RES:HMBSA). These body site-specific assemblies are aimed at reconstructing organisms that represent too small a fraction in any individual sample to assemble but are found among many individuals. For 12 stool samples both Illumina and 454 FLX Titanium data (RES:HM4WGS) were generated allowing a hybrid assembly approach, using Newbler (Supplementary Information) (RES:HMHASM). Overall, the assembly statistics recovered varied substantially by body site and community complexity (Supplementary Fig. 5). However, our results indicate that, for the assembly strategy we employed, metagenomic assembly quality plateaus at approximately 6 Gb of microbial sequence coverage for a sample possessing a microbial community structure similar to that of stool samples (Supplementary Fig. 6).

A WGS based perspective of community membership was obtained by aligning the reads to a set of 1,742 finished bacterial, 131 archaeal, 3,683 viral, and 326 lower eukaryotic reference genomes (RES:HMREFG) (Supplementary Information, [14]) representing a broad taxonomic range from each of these four domains. A total of 57.6% of the high quality microbial reads could be associated with a known genome (ranging from 33 to 77%, for anterior nares and posterior fornix, respectively) (RES:HMSCP). The overwhelming majority of mapped sequences originated from bacteria (99.7%), while the remaining reads mapped to microeukaryotes (0.3%) or archaea (< 0.01%) (Supplementary Information).

Two complementary approaches were used to summarize overall function and metabolism of the human microbiome producing two primary data sets of annotations (Supplementary Information, RES:HMMRC, RES:HMGI) and additional secondary analyses (Supplementary Information, RES:HMGS, HMHGI, HMGC, HMGOI) available to the community for further interrogation. The first primary data set of annotations was produced by mapping individual shotgun reads to characterized protein families (RES:HMMRC) [15]. The second, was produced from functionally annotated gene predictions generated from the metagenomic assemblies (RES:HMGI) which were subsequently grouped according to high-level biological processes and to selected additional processes specific to metabolism and regulation (RES:HMGS) [16] (Supplementary Tables 6 and 7, Supplementary Fig. 7).

HMP data generation and analysis lessons

A key manner in which the HMP resources will serve to guide future studies of the microbiome is by enabling informed decisions regarding sampling protocols and genomic DNA preparation [6], sequencing depth, [13] statistical power, [17] and metagenomic data type. As indicated in Table 1, the consortium successfully amplified 16S sequences to our target depth at all 18 body sites, with the fewest sequences recovered consistently from the antecubital fossae. The amount of host human DNA recovered and the finest level of OTU resolution varied for 16S sequences among body sites (Supplementary Figs. 3-4, [7]).

From our WGS investigations, a series of protocols (http://hmpdacc.org/tools_protocols/tools_protocols.php) have been established to process large volumes of short read WGS data and to annotate and examine these data through both a multi-tiered assembly approach and as single reads [18]. An investigator’s choice of metagenomic technologies can thus be guided not only by a 16S versus WGS dichotomy, but also by the expected fraction of host sequence and appropriate 16S region targeting the dominant taxa at each body site (Supplementary Figs. 2-6, 8).

Together, these datasets represent comprehensive and complementary views of the human microbiome as shown by comparing organismal (Fig. 1a) and gene (Fig. 1b) catalogues, and the ratio of genes contributed per OTU (Fig. 1c). The discovery rate of new gene clusters (as determined by annotation of assembled WGS data) is in general detected more slowly relative to organismal discovery (as determined by OTU data) due to the fragmentary nature of these community reads and assemblies despite high sequence depth (Fig. 1a and b, Supplementary Fig. 9) and the number of genes contributed per OTU varies by body site (Fig. 1c, Supplementary Information). However in general, these results highlight an important point for consideration of further microbiome investigations using these data sets as they suggest that the majority of the common taxa and genes present in this reference population have been detected.

Figure 1

Rates of gene and OTU discovery from HMP taxonomic and metagenomic data

Accumulation curves for a, OTU counts from 16S data (all body sites) b, clustered gene index counts from metagenomic data (all applicable body sites) and c, the ratio of average unique genes contributed versus unique OTUs encountered with increasing sample counts (Supplementary Information). Ratios given for each curve in c represent the average number of unique genes contributed per unique OTU at the final sample count. Curves for stool, buccal mucosa and anterior nares suggest that the proportion of gene-to-taxa discovery has stabilized. In contrast, the curve for supragingival plaque suggests relatively fewer new genes are being contributed per additional OTU. Error bars represent 95% confidence intervals.

We additionally compared the gut community gene catalogue sampled by the HMP with that of MetaHIT in terms of total detected gene counts. The HMP recovered more total non-redundant gene counts (5,140,472) than reported by MetaHIT (3,299,822) [5] likely reflecting a combination of the increased sequence depth obtained by the HMP (11.7 Gb HMP, 4.5Gb MetaHIT on average) and differences in data generation and processing [5].

The two non-redundant sets of gene sequences were subsequently combined and compared by matches to a database of orthologous groups [19] of functionally annotated genes. Approximately 57% of the orthologous groups recovered by this method overlapped between the data sets, while an additional 34% versus 10% were unique to the HMP and MetaHIT, respectively (Supplementary Fig. 10, Supplementary Table 8, Supplementary Information). After removal of genes that received any orthologous group assignment, the remaining novel genes were subsequently clustered [20]. Approximately 79% of the HMP-derived novel gene clusters were orthologous to one or more clusters in MetaHIT, while an additional 16% were unique to this study versus 5% for MetaHIT-derived data (Supplementary Fig. 11, Supplementary Table 8, Supplementary Information, [5]). These results suggest that for this body habitat, relatively similar gene catalogues were recovered despite differences in experimental design and protocols. However, a greater proportion of both annotated and unique novel genes were detected in the HMP data set, emphasizing the utility of sequencing depth in recovering gene function and in particular, deriving rare function. These results further underscore the importance of large-scale sequenced based studies of the microbiome to better characterize its gene content and diversity.

Human microbiome reference genomes

The current goal for the reference genome component of the HMP is to sequence at least 3,000 reference microbial genomes associated with the human body. Thus far, more than 800 genomes have been sequenced and are available from NCBI and the DACC (http://hmpdacc.org/HMRGD). From an alignment of WGS reads to reference genomes (RES:HMREFG), approximately 26% from the total read set (46% of all reads that could be aligned) were matched to a subset of 223 HMP reference genomes (Supplementary Information, Supplementary Data).

We continue to solicit community feedback for strains that will best benefit our attempts at understanding the breadth of human microbiome diversity. For example, a prioritized list of the “most wanted” HMP taxa is currently being maintained (http://hmpdacc.org/most_wanted/) with the goal of targeting these difficult to obtain organisms using both culture-based and single-cell approaches.

A catalogue of all HMP reference genomes along with custom filtering, viewing, graphing, and download options can be found at the DACC Project Catalogue (http://www.hmpdacc-resources.org/hmp_catalog/main.cgi). In addition, comparative analyses of reference genomes are provided by the data warehouse and analytical systems, IMG/HMP (http://www.hmpdacc-resources.org/cgi-bin/imgm_hmp/main.cgi). Cultures of all HMP reference strains are required to be made publicly available through, the Biodefense and Emerging Infections Research Resources Repository (BEI). Information on strain acquisition can be found at the DACC (http://hmpdacc.org/reference_genomes/reference_genomes.php) and BEI (http://www.beiresources.org/tabid/1901/stabid/1901/CollectionLinkID/4/Default.aspx).

Conclusion

An overarching goal of this multi-year, multi-center project is the generation of a community resource to advance research efforts related to the microbiome. The result is a collection of 11,174 primary biological specimens representing the human microbiome as well as corresponding blood samples from the human donors which are being reserved for sequencing at a future date and from which cell lines will be developed. A variety of new protocols were developed to enable a project of this scope; these include methods for donor recruitment, laboratory and sequence processing, and analysis of 16S and WGS sequence and profiles. These resources serve as models to guide the design of similar projects. Studies with a primary focus on disease can use this reference for comparative purposes including detecting shifts in microbial taxonomic and functional profiles, or identification of new species not present in healthy cohorts that appear under disease conditions. The catalogue described in this study is the largest reference set to date of human microbiome data associated with healthy individuals. Collectively the data represents a treasure trove that can be continually mined to identify new organisms, gene functions, and metabolic and regulatory networks as well as correlations between microbial community structure and health and disease [4]. Among other future benefits, this resource may promote the development of novel prophylactic strategies such as the application of prebiotics and probiotics to foster human health.