Direct sequencing of the human microbiome readily reveals community differences

In the past few years, the availability of improved sequencing methods, including pyrosequencing [1], has revolutionized what we know about the microbes that inhabit our bodies. Although it has been known for decades that our microbial symbionts outnumber our own cellsby about a factor of 10 [2], the differences in the repertoires ofsymbiontsharbored by different healthy individuals, different siteswithin the individual, and by individuals over time are only now coming to light. Initially, it was assumed that a 'core microbiome' existed; that is, that a substantial number of microbial species was shared in each body habitat in all or most humans, and that the genomes of these core species could be used as scaffolds to assemble fragmentary data from short-read shotgun sequencing of microbial community DNA [3].

The first three individuals whose gut microbiomes were surveyed using substantial numbers of 16S rRNA genesequences shared few of their species, however [4]. Similarly, observations that a person's left and right hands have only 17% of bacterial species in common, and that two different people's hands share only 13% [5], cast doubt on the concept of a substantial core set of microbial species shared by all or most people. This doubt has been reinforced by recent work that redefines core lineages or genes as 'core' even if shared by relatively few people [6, 7]. In fact, on the basis of 16S rRNA geneanalyses we can rule out the possibility that, even within relatively homogeneous small populations of fewer than 100 individuals, everyone's skin-surface communities or gut communities share more than a tiny fraction of species [6–8]. This unanticipated variability in shared community membership, and also in other important aspects of the human microbiome, poses substantial conceptual and computational challenges.

Of particular importance for microbiome studies is the following question: what is the effect size? That is, using standard terminology from statistics, how distinguishable are two communities or groups of communities? Obtaining an answer is essential for addressing many practical concerns with experimental design. For example, the effect size determineshow many individuals need to be recruited for a given study, and how many sequences need to be collected per sample to observe differences if they exist. These considerations are particularly importantfor the study ofsystemic disorders such as diabetes or some autoimmune disorders, which are expected to influence the microbiomein multiple body habitats. We need a sense of how much variation exists among different body habitats, how much variation is observed among healthy individuals for the same body habitat, and how much of a shift occurs due to a pathophysiologic state. It is also importantto define the most appropriate method for determining the magnitude of similarity or difference between communities, as the choice of method has a large influence on the results of community comparisons [9–12]. A general discussion of the pros and cons of different metrics of community overlap is beyond the scope of this paper (see [9–12] for reviews). Here, we summarize the types and sizes of effects found in studies that used various methods of comparing groups of samples, and look for large-scale patterns that can give information on the number of individuals and sequences that are needed to observe different types of effects (Figure 1).

A variety of interrelated features differentiate microbial communities. These features include the the relative abundance of specific taxa (the proportion of the bacteria in the sample that are Firmicutes, for example), the level of species richness or diversity observed within a community (alpha diversity), and the degree to which different communities share membership or structure (beta diversity). A major challenge in comparing studies is that there is no consistent way in which the size of community differences is reported, as the type of difference that is relevant depends on the study. For example, lean and obese mice and humans differ in their ratios of prominent bacterial phyla (Bacteroidetes (which include the common gut commensal Bacteroides), Firmicutes (Gram-positive bacteria, including Lactobacillus and Clostridium), and Actinobacteria (which include Corynebacteria and Mycobacteria) [13–15]); men's and women's hands differ in the number of species-level phylotypes (defined as organisms with 16S sequence identity >97%) observed on average [5]; and samples from the same or similar sites on the bodies of different individuals cluster together using UniFrac-based principal coordinates analysis [4, 16, 17]. UniFrac is a metric for comparing microbial communities using phylogenetic information, which has been implemented in several tools.

Because of the diverse ways in which microbial communities respond to various environmental factors, it is difficult to compare effect sizes across different studies or systems, as an analysis that highlights differences in one system may obscure them in another. Thus, in what follows, we review effect types and sizes as reported by the authors of individual studies. We focus on variation in human-associated microbial community diversity as assessed by 16S rRNA gene sequence surveys of abundant lineages, using various measures of both within- and between-sample diversity (alpha and beta diversity, respectively). We review comparisons of microbial communities in relationship toboth sampling depth (that is, number of sequences per sample) and breadth (that is, number of samples or individuals). We then perform simulations using an atlas of microbes associated with different sites in the human body to ask how many sequences per sample are needed in order to detect differences across individuals, time, and locations within the body.

Although a wide range of studies in healthy subjects have identified substantial interpersonal variation in overall microbial community composition, how do these effect sizes compare with differences correlated with disease, or in response to treatments ofvarious environmental samples? To address this question, we reviewed culture-independent, 16S rRNAgene-based surveys associated with different physiological conditions (Table 1b) and associated with experimental manipulations in non-human environments (which were surprisingly scarce; Table 1c).

One of the best-characterized effects of health status on the gut microbiome is the association between obesity and the proportional representation of Bacteroidetes, Firmicutes and Actinobacteria [6, 13–15]. Studies in mice indicate that the microbiota contributes to the obese state by providing the host with a greater amount of energy from the diet compared with the microbiota of a lean host [15], as well as by manipulating host genes that regulate the deposition of energy in adipocytes [19]. The obesity-associated microbiomes of humans (and mice) are enriched in functional genes for certain types of carbohydrate metabolism, and this is directly attributable to the reduction in the numbers of genomes of members of the Bacteroidetes [6, 15].

However, even the size of the differences in gut bacterial community composition of obese versus lean hosts is debated, as different studies using different methodologies have returned varied results [20]. The impact of methodology is particularly evident in a study of twins concordant for obesity or leanness, in which the observed relative abundances of Bacteroidetes, Actinobacteria and Firmicutes, as judged by sequencing of different regions of 16S rRNA clones, depended on the sequencing approach - pyrosequencing of PCR products, Sanger sequencing of 16S rRNA clones, or shotgun sequencing and phylogenetic classification of reads [6]. However, the direction of the effect was consistent across methodologies, and detectable with as few as a couple of hundred sequences per sample.

Observable phenotypes such as obesity may be caused by a variety of underlying factors, and which of those factors is responsible for shifts in the host's microbiota is difficult to address in such correlative studies. Experimental manipulations of microbial communities, however, allow determination of the relative effects of specific variables on overall community composition or the abundance of particular taxa, and as such, allow researchers to draw conclusions regarding cause and effect. Examples of experimental manipulations of non-human environments that used 16S rRNA gene sequencing approaches (either clone libraries or pyrosequencing) and that were well enough replicated to allow statistical analysis are shown in Table 1c. For soil samples, three to four replicates with 70 to 100 sequences were sufficient to observe differences in microbial communities due to land use and moisture regimes [21, 22]. For piglet gut microbiota, the effects of antibiotics on overall community composition were evident with as few as 96 sequences per sample [23]. It would be fascinating to test whether similar antibiotic-induced effects in outbred populations of humans with diverse diets [24] can be found with relatively few sequences. Similarly, it would be important to consider sampling depth under human physiological conditions in cases where the effect size is known to be large, for example, in the development of the infant gut microbiota [25].

The literature reviewed in Table 1 reports how many sequences were used to reveal a variety of different effects. Could the same results have been achieved with less sequencing? To begin to address this question, we carried out a limited reanalysis of a study of multiple body habitats by Costello et al. [16], which encompasses variability explained by nested factors with different effect sizes (Box 1).

In conclusion, the results described here, and previously reported [8, 37], show that arbitrarily choosing to generate large numbers of sequences may not be the most cost-effective way to identify changes in microbial communities associated with different physiological or pathophysiological states. Instead, we call for a few standardized methods to assess differences among microbial communities, which will allow for effect size and power calculations, and therefore a considered assessment of the number of individuals and sequences required to differentiate among given communities. The following four methods have been successful in a range of studies: differences in alpha diversity (number of phylotypes observed or extrapolated); differences in abundance of specific lineages; differences in location on a principal coordinates plot obtained from UniFrac distances or other metrics; and the F _ST measure described in the previous section.

The rapid increase in sequencing capacity provides a spectacular opportunity to advance the field in ways that were unimaginable even 3 years ago. How can individual investigators, or groups of investigators, use these resources most wisely at this unique moment of democratization of the ability to perform sequence-based studies? The data summarized here suggest that study designs consisting of tens of thousands of samples sequenced at shallow coverage will be highly informative (depending on the effect size), and such studies are possible with the instruments available today. Given recent observations that inter-habitat and inter-personal variations are large effects, we believe that individual researchers can and should sieze the opportunity provided by these findings to analyze vast numbers of samples at low-coverage (for example, 100 to 1,000 sequences). At this number of samples, detailed exploration of spatial and temporal dynamics of microbial communities will be possible, as will comparisons of large patient populations. In addition, replicate samples can be acquired and analyzed without too strongly impairing the breadth of an investigation, allowing more robust experimental designs to be implemented. One can envisage that perhaps within the next few years, a group of motivated high-school students might, for a science-fair project, be able to track movements in microbes between humans and their pets and livestock across the planet. These studies, especially when combined with hypothesis-driven approches to understanding the effects of factors such as diet and antibiotic exposure, could go far beyond even the largest purely observational studies being contemplated today.

Such studies will yield an overall map of variation within the human microbial ecosystem, and relate differences to specific physiological states within and between individuals in a manner that is replicated across individuals. These studies will serve as a framework to identify and compare the shifts that take place in the microbial community that are related to specific disorders.