Skip to main content

Experimental design and quantitative analysis of microbial community multiomics


Studies of the microbiome have become increasingly sophisticated, and multiple sequence-based, molecular methods as well as culture-based methods exist for population-scale microbiome profiles. To link the resulting host and microbial data types to human health, several experimental design considerations, data analysis challenges, and statistical epidemiological approaches must be addressed. Here, we survey current best practices for experimental design in microbiome molecular epidemiology, including technologies for generating, analyzing, and integrating microbiome multiomics data. We highlight studies that have identified molecular bioactives that influence human health, and we suggest steps for scaling translational microbiome research to high-throughput target discovery across large populations.


Population-scale studies of the human microbiome now have at their disposal a remarkable range of culture-independent and other molecular and cellular biology technologies, but the identification of elements of the microbiome that are functionally important for human health remains challenging. This is in part due to the variety of tools available and the diversity of processes that they measure: microbial community composition [1,2,3], species and strain diversity [4,5,6,7], genomic elements [8, 9], transcription, translation, and metabolism [10,11,12], along with the corresponding human molecular processes in multiple epithelial, immune, and other cell types [13,14,15]. Research challenges also arise, however, at the intersection of microbial ecology and molecular epidemiology, as population-scale microbiome study designs and methods that adequately account for human variability, environmental exposures, and technical reproducibility are also still in the early stages of development [14,17,, 1618].

Existing technologies for population-scale microbiome studies share many similarities with molecular epidemiology techniques for human gene expression and genome-wide association studies [19, 20]. Human-associated microbial communities are most often profiled in terms of their composition, for example by sequencing the 16S ribosomal RNA (rRNA) genes to yield phylogenetic or taxonomic profiles (abbreviated here as 16S amplicon profiling) [21]. 16S and other amplicon-based technologies [22] are limited in their phylogenetic ranges; for example, 16S rRNA gene studies primarily target bacteria, with some crossover, whereas 18S or internal transcribed spacer (ITS) studies typically target fungi. Although highly sensitive, these technologies also suffer from contamination, amplification, and extraction biases [23]. A subset of these issues are shared by whole-community shotgun metagenomic sequencing approaches, which can further describe the functional genetic potential of the entire community, but do not tell us what portion of this genetic potential is actively transcribed or translated in any particular environment [24, 25]. Community metatranscriptomics, metabolomics, and metaproteomics techniques are emerging to link nucleotide sequence-based profiles to their bioactive products [26, 27], as are complementary technologies such as immunoglobulin A gene sequencing (IgA-seq), immunoprofiling, and human cell screening techniques to jointly profile microbial and human host activities [13, 28, 29]. When combined with culture-based microbial characterization [30], recent advances in the resulting experimental toolkit have greatly improved our ability to identify relevant components of host–microbiome interactions.

Translational applications of the microbiome at the population scale, however, require careful experimental, computational, and statistical considerations, combining lessons learned from earlier molecular epidemiology with challenges unique to microbiome profiling. First, the identification of relevant human or microbial cellular and molecular mechanisms requires sufficiently precise technologies; if bioactivity is due to a particular microbial strain or transcript, for example, it is unlikely to be identified by amplicon sequencing. Next, the identification of signals that are sufficiently reproducible for clinical actionability requires well-powered experimental designs and, ideally, meta-analysis among studies—both challenging for current microbiome protocols. Many environmental exposures and covariates, such as diet or medications, must also be measured because the microbiome (unlike the human genome) can both modify and be modified by these factors. Finally, appropriate computational and statistical methods must be used during analysis, as many standard approaches can be prone to surprising false positive or negative rates. In this review, we thus detail the current best practices in this field with respect to these challenges, delineate methods and computational tools (or lack thereof) for addressing these challenges, and discuss potential future directions for conducting integrated multiomics studies in microbiome molecular epidemiology.

Microbial strain as the fundamental epidemiological unit for microbiome taxonomic profiles

It has become increasingly apparent that many, although not all, analyses of translational activities in the human microbiome will require the identification and characterization of microbial taxa at the strain level. Many current culture-independent tools profile microbial community membership by delineating genera or species, but microbial epidemiologists have long recognized that not all strains within a species are equally functional, particularly with respect to pathogenicity. For example, Escherichia coli may be neutral to the host, enterohemorrhagic [9], or probiotic [31], and epidemiologists have long employed methods such as serotyping, phage typing, or pulse gel electrophoresis to reveal and track the relationships between microbial strains within single species (as opposed to communities) of interest. Indeed, there is enormous genomic variation within E. coli alone; studies suggest a pangenome of well over 16,000 genes, with ~ 3000 gene families present in most strains and fewer than 2000 universal genes [32, 33]. While more comprehensively characterized for Escherichia than for other genera, this variability is not atypical of many microbial species.

Critically, such inter-strain variation has phenotypic consequences for human health, even in such well-studied organisms as E. coli. For instance, the probiotic strain E. coli Nissle was isolated during World War I due to its ability to confer resistance to Shigella upon its host [31], despite the close relationship of this strain to the uropathogenic strain CFT073 [34]. Escherichia is not unique among human commensals in having a large pangenome with a relatively small core. The Staphylococcus aureus pangenome is also approximately five times larger than its core genome [35], and this variation likewise has important consequences in differentiating commensal staphylococci from methicillin-resistant S. aureus (MRSA) [36]. Even gut commensals that are not traditionally associated with pathogenicity, such as Bacteroides vulgatus [6, 37], may show large intra-species genomic variation. Like those of better-characterized pathogens, these genomic differences within commensal microbe species may have consequences for the host; for example, not only was Prevotella copri recently correlated with new-onset rheumatoid arthritis, but specific gene differences among P. copri strains were also correlated with this phenotype [38].

Although strain differences can have profound implications for human health, culture-independent tools have only recently begun to distinguish among strains during taxonomic profiling (Fig. 1ac). For example, amplicon analyses are fundamentally limited in their ability to differentiate strains because critical functionality may arise from differences that occur outside of the otherwise-identical amplified gene regions (e.g., plasmids in Escherichia and Shigella). Both shotgun metagenomics and, when possible, 16S-based approaches can now be used to discriminate strains (Table 1), although both (especially the former) require care during such analyses. Most traditional operational taxonomic unit (OTU) clustering approaches for amplicon data, for example, differentiate only among taxa above some nucleotide identity threshold (e.g., 97% similarity). Likewise, metagenomic assembly protocols may intentionally avoid nucleotide-level variants. For 16S data, newer approaches [39,40,41] employ novel algorithms to distinguish between biological signal and sequencing error, and can discriminate small sequence differences corresponding to large phenotypic differences, such as sponge symbionts and their choice of host [39], or the specific ecological niches of human oral taxa [42]. Recent progress in developing bioinformatic tools further improves this resolution, revealing strain-level differentiation within the 16S region that can be as small as a single nucleotide [43,44,45].

Fig. 1

Strategies for detailed strain and molecular functional profiling of the microbiome in human population studies. a Culture-independent analysis methods can now identify members of the microbiome at the strain level using any of several related techniques. This is important in population studies as strains are often the functional units at which specific members of microbial communities can be causal in human health outcomes. b Among different approaches, reference-based methods can require less metagenomic sequence coverage (as little as ~ 1×), but are limited to identifying variation that is based on genes or single nucleotide variants (SNVs) related to available reference genomes. c Assembly-based methods can additionally resolve syntenic information across multiple markers at the cost of higher coverage (≥10×, Table 1). d , e Metatranscriptomic analysis, another emerging tool for characterizing microbiome function in human health, reveals over- or under-expression of microbial features with respect to their genomic content, both on d the population and e the individual level. ORF open reading frame

Table 1 Tools for metagenomic strain analysis

Algorithms for strain identification from shotgun metagenomic sequences generally rely on one or both of two techniques: calling single nucleotide variants (SNVs, within a community or between community members and reference genomes) or identifying variable regions (such as gained or lost genomic elements; Table 1). Community SNV identification, like microbial isolate or human genetic profiling, requires sufficiently deep coverage (typically 10× or more) of every microbial strain to be differentiated [5], but can delineate closely related strains very precisely. SNVs can be assessed either extrinsically, with respect to one or more reference sequences (e.g., by mapping metagenomic sequences to that of reference and calling SNVs) [5], or intrinsically, by aligning sequences directly from one or more metagenomes and identifying SNVs among them [4]. Finally, as microbial strains often differ dramatically in their carriage of different core or pangenome elements or genomic islands (unlike most populations within eukaryote species [46]), strains can also be identified by the presence or absence of one or more genes or genomic regions [6]. This requires less sequencing depth (and is thus sensitive to less abundant members of a community), but can be more susceptible to noise and unable to delineate closely related strains.

Although strain identification, characterization, and phylogenetics are well-developed for microbial isolates [47], the use of culture-independent amplicon or metagenomic sequence data to perform such tasks is still in its infancy and can suffer from a variety of drawbacks. Amplicon methods in particular require variation to exist in the targeted region, and detecting the few variants that might exist in such short sequences requires extremely careful data generation and analysis protocols to distinguish biological from technical variation [39, 40]. Metagenomic strain identification is typically only accurate for the single most dominant strain of any one organism in complex communities, requiring extreme sequencing depths (e.g., tens to hundreds of gigabases) to differentiate secondary strains except when only one or a few organisms dominate [5]. Finally, as in other areas of microbial genomics, metagenomic strain identification is sensitive to the definition of a 'strain', which can vary from clonality at all genomic loci (possibly including plasmids), clonality at all sequenced locations (possibly only within an amplified region), or allowing some non-zero degree of nucleotide-level divergence [48].

Metatranscriptomics enables characterization of context-specific, dynamic, biomolecular activity in microbial communities

Taxonomic profiling, at any level of resolution, is increasingly accompanied by functional profiling—pairing a community's organismal makeup with its gene and/or pathway catalog [9]. Metagenomic DNA sequencing, however, yields information only regarding the community's functional potential—which organisms, at what abundances, might be able to carry out which biological processes (and not necessarily which genes are being transcribed under current conditions). Metatranscriptomic RNA sequencing is arguably the first scalable, culture-independent technology to overcome this limitation, although its application to the human microbiome at an epidemiological scale still presents unique design and analysis challenges. Microbiome samples for metatranscriptomics must be collected in a manner that preserves RNA for sequencing, and they are (by definition) much more sensitive to the exact circumstances and timing of sample collection (Box 1) [17]. The associated protocols for nucleotide extraction are generally more challenging and sensitive to technical variability [49]. The resulting metatranscriptomes must generally be accompanied by paired metagenomes in order to allow interpretation of the data, otherwise changes in DNA copy number (i.e., microbial growth) cannot be differentiated from changes in transcriptional activity [24]. This is particularly true for amplicon-based rRNA metatranscriptomics, a proposed proxy for organismal growth or metabolic activity within a community [50]. In such settings, it is not yet clear how we could account for 16S rRNA gene copy number variation, differences in ribosomal transcription rates, or even the exact biological interpretation of 16S rRNA transcript abundances (as opposed to gene abundances as profiled by typical DNA amplicon sequencing).

By contrast, shotgun metatranscriptome studies provide biological information that complements metagenome studies, including detection of RNA viruses and quantification of rare but functional genes that might remain undetected in DNA-based metagenomic surveys [51] (Fig. 1d and e, and Table 2). Metatranscriptomic sequencing can also highlight the taxon- and strain-specific transcriptional activity of a community, providing a comprehensive overview of the functional ecology of the microbiome (Box 2). A typical metatranscriptomic study, such as a single-microbe RNA-seq study [52], consists of several steps, including: 1) transcript mapping and/or assembly; 2) annotation with functional and/or taxonomic information; 3) normalization; and 4) differential expression analysis. When processing reads, a metatranscriptomic analysis pipeline typically either maps reads to a reference genome or performs de novo assembly of the reads into transcript contigs. The first approach (mapping to a reference genome) is limited by the information in the reference database, whereas the second approach (de novo assembly) is limited by the difficulty of assembling long contigs of highly variable transcriptional coverage from complex metagenomes. Downstream bioinformatic analysis of metatranscriptomic expression profiles must further account for taxonomic composition variations and for technical biases associated with RNA-seq experiments. In particular, taxon-specific rescaling (RNA transcript abundance normalized to its DNA copy number) is a necessary step in order to ascertain whether apparent shifts in transcript levels are concordant with changes in taxon abundances. Finally, to conduct differential gene expression analysis post-normalization, off-the-shelf tools from single-organism RNA-seq can be used, some of which have already been adapted to microbial community settings [53].

Table 2 Tools for primary processing of metatranscriptomes
Table 3 Tools for feature-wise differential abundance analysis in microbial community taxonomic profiles

Microbiome-associated metabolomics as an emerging opportunity to characterize bioactivity

Although several other culture-independent molecular methods are now joining metatranscriptomics for human microbiome profiling, non-targeted metabolomics may represent one of the most successful to date in explaining the mechanisms of bioactivity [26, 68]. This includes a range of nuclear magnetic resonance (NMR) and mass-spectrometry technologies for profiling small molecules from stool [26, 68], skin [69], circulating metabolites [70, 71], or coupled with other human-associated microbial communities. In many of these environments, it has been estimated that over 10% of small molecules may be of microbial origin or microbially modified [72], highlighting the need to associate specific microbial strains or genetic elements with the specific small molecules that, in turn, mediate human health phenotypes. The associated study designs have as yet seen limited application at the population scale, with some success stories highlighted below, and it remains to be seen which microbiome-associated metabolites are appropriate for predicting or modulating population health outcomes. The resulting data share similar strengths and weaknesses to metatranscriptomics; protocols are often still technically challenging, and while the resulting data may be more difficult to characterize at the molecular level, when possible they represent measurements that are often more directly causal (e.g., small molecules responsible for a specific bioactivity).

Statistical questions, issues, and practice in modern epidemiological microbiome studies

In all of these approaches—amplicon-based, shotgun sequencing, or other technologies—the persistent goal of microbiome epidemiology has been to determine whether and how microbial and molecular feature abundances are associated with the certain characteristics of the samples, such as donor health, disease status or outcome, donor dietary intake, donor medication, or environment (Fig. 2ad). This translation of molecular epidemiology to the setting of the microbiome is challenging for several reasons. Among these is the technical nature of data associated with microbial communities, which typically consist of counts that have a compositional structure. That is, microbiome sample data (of most types) are frequently represented as vectors of fractional relative abundances (the total of all features in a sample sum to a value such as 1 or 100%). When typical statistical inference methods are used on compositional data, false positives result as a consequence of spurious correlation. This problem is exacerbated in population-scale microbiome studies by high data dimensionality (up to tens of thousands of samples containing potentially millions of microbial features), sparsity (made more challenging as the result of a mix of true zeros and undersampling events), and mean-variance dependency (variance of counts changes with the value of the mean) [63]. Failure to account for these specific characteristics of microbiome count data during statistical analysis can lead to strong biases in results; in particular, false positives outcomes are common, leading to irreproducible associations even (or especially) in large cohorts [73].

Fig. 2

Microbiome molecular epidemiology. a Multiomic profiling of host and microbiota enables in-depth characterization of community properties from multiple culture-independent data types (including metagenomics, metatranscriptomics, metaproteomics, and metametabolomics) to address questions concerning the microbiome’s composition and function. b As in host-targeted molecular epidemiology, metagenomic and other metaomic data types can be integrated and associated with the available metadata to provide a comprehensive mechanistic understanding of the microbiome. c A wide range of early-stage data analysis choices can strongly affect microbial community data analysis, including the quality control of the raw data, the normalization of the raw data, choice of host and microbial features to extract, and algorithms to profile them. A hypothetical example of four taxonomic features is shown derived from four samples with differing metagenomic sequencing depths (top). Features with the same relative abundances may thus appear to be different on an absolute scale because larger sequencing depth can generate larger read counts (top). Normalization also corrects for potential batch effects and helps to preserve meaningful signal between cases and controls (bottom). Note that the precise methods used for global visualizations, such as the ordination method, can dramatically affect how the data are summarized, as can important parameters in the process, such as the (dis)similarity measures used to compare features or samples. d Within an individual study, the integration of multiple metaomic data types can provide stronger collective support for a hypothesis. Here, a hypothetical disease association is shown at the DNA, RNA, and protein or metabolite levels, providing a more complete picture of the disease pathogenesis. e When they differ between datasets, the strong technical effects that the choices mentioned above have on individual studies can impede multi-study meta-analyses, making this type of population-scale analysis difficult in the microbiome. When possible, the meta-analysis of host and microbial features with respect to shared phenotypes of interest can allow more confidence in prioritizing microbial taxa, gene products, or small molecules that have statistically significant roles in disease relative to covariates. f Finally, as with genome-wide association studies, it is critical to validate putative associations of top candidate microbial features with follow-up experimentation. In the microbiome, this can include studies involving animal models (such as gnotobiotic mice), mammalian cell systems, and/or microbial cultures

Several analysis methods have been developed to specifically address these problems in tests for differential feature abundance in the microbiome (Table 3 and Box 3). Virtually all of these methods rely on some form of normalization, and they differ primarily in the choice of the data transformation, statistical model, and null distribution (or equivalent) for p value calculation. For example, metagenomeSeq [74] takes raw read counts as input and accounts for possible biases using a zero-inflated Gaussian mixture model to integrate normalization and differential abundance analysis of log-counts. MaAsLin [75] uses a variance-stabilizing arcsine square root transformation to create continuous abundance profiles that can be analyzed by regular linear models. Apart from these community-specific tools, methods developed for differential expression analysis of similar RNA-seq data—such as edgeR [76], DESeq2 [77], and limma-voom [78]—have been adopted in microbiome research. These methods are typically based on a negative binomial statistical model of the normalized counts (with the exception of limma-voom, which applies an empirical Bayes linear model to the normalized counts) [53, 79]. Apart from these parametric approaches, several non-parametric alternatives have also been developed, such as LEfSe [80], Metastats [81], and ANCOM [82]. These methods make minimal assumptions about the data and estimate the null distribution for inference from ranks or from the observed data alone.

Normalization plays a crucial role in differential abundance analysis because variation in sequencing depth can make read counts incomparable across samples. Directly comparing read counts among samples with different sequencing depths may lead to the false conclusion that features are differentially abundant even when they have the same composition. In addition to simple total sum scaling (TSS) or rarefaction, this has led to the development of a variety of normalization approaches, such as trimmed mean of M-values (TMM) [83], relative log expression (RLE) [84], and cumulative sum scaling (CSS) [74], that aim to address the heteroscedasticity of the samples by variance stabilization and robustification or filtering [53]. Rarefaction is not ideal for many purposes because of its lack of statistical power and the existence of more appropriate methods [53], but it is fast and can be reasonably accurate in approximating a reliable normalization when necessary, especially given sufficient sequencing depth.

Table 4 Tools for compositionality-aware ecological network inference

Given the prominence of multivariate metadata in modern epidemiological cohorts, the availability of multivariable analysis tools is becoming increasingly important in the microbiome research community (Boxes 3 and 4). Some methods for differential abundance testing can only detect univariate associations, whereas other methods, such as edgeR, DESeq2, metagenomeSeq, limma-voom, and MaAsLin, can perform multivariable association. Future microbiome analytical tools must further leverage the hierarchical, spatial, and temporal nature of modern study designs, which typically result from repeated measurements across subjects, body sites, and time points. Several recent studies have taken initial steps to address one or both of these issues. One avenue of research aims to capture the correlation among repeated measurements by using random effects [75, 78, 85, 86]; other studies have relied on dynamical system or probabilistic spline modeling [87] of microbiome time-series data to study the temporal dynamics and stability of microbial ecosystems. Despite these innovations, the longitudinal modeling of microbiome data is still in its infancy, particularly in combination with multiple covariates in large human populations. There is a dearth of systematic studies aimed at the evaluation of multiple-covariate, repeated-measure methods for microbiome epidemiology, with no clear consensus to date. As microbiome data continue to accumulate, there is a pressing need for a rigorous comparison of these multivariable tools to help guide experimental designers and meta-analysts.

Many current microbiome epidemiology studies also use unsupervised models or visualizations to reveal structural patterns. Ordination is a particularly common visualization technique [21] that aims to plot samples in a low-dimensional space (usually no more than three axes) that also reflects their overall community similarities. This enables intuitive but rough inspection of strong signals in microbiome data (for example, an analyst might quickly identify samples with certain common characteristics that also have similar microbial compositions). Clustering analysis, also referred to as enterotyping or identifying community state types [88,89,90], is a related unsupervised technique for separating samples that have distinct profiles into different groups ('clusters'), and is appropriate only when distinct microbial sub-classes reliably exist in the data. Both methods have been heavily explored in high-dimensional biological datasets, such as gene expression and single-cell sequencing datasets, and while they can provide powerful tools for data overview and hypothesis generation, it is also important to recognize their limitations. First, both ordination and clustering analyses rely on a sample-against-sample dissimilarity (i.e., beta-diversity) matrix as input, and are thus sensitive to the choice of dissimilarity measure [73]. Second, as unsupervised approaches, both come with a wide variety of tunable parameters that are difficult to evaluate objectively. Third, for clustering analysis, distinguishing between discrete and continuous sample distribution patterns can be challenging when sample size is limited and/or signal is weak. Under such circumstances, quantitative examination of clustering strength is important to ensure that the identified clusters actually exist [89]. Finally, both methods are best suited to identifying the strongest patterns driven by population-level characteristics, both for microbiome data and in other ’omics settings [21]. To identify microbial associations with an outcome variable, supervised analysis [91] provides the resolution needed to identify patterns that might not be captured by the single strongest axis of variation, as well as rigorous, statistically justified quantification of such associations.

To this end, several families of omnibus test assess whether overall patterns of microbial variation in a community associate with covariates by some significance model (e.g., PERMANOVA [92], MiRKAT [93], ANOSIM [94]), typically with the ability to adjust for additional covariates. These tests are complementary to the supervised per-feature epidemiological association tests described above. They also take beta-diversity matrices as input, and they adopt statistically justified procedures to evaluate significance against the null hypothesis that covariates are not associated with overall microbiome composition. This is in contrast to the use of multiple individual tests for each microbial feature (species, clade, pathway, and so on) independently with respect to covariates, as described above. Similarly to ordination, the choice of dissimilarity measure can affect results, and some methods [93, 95] have correspondingly developed extensions to incorporate multiple metrics simultaneously in order to improve robustness. Another limitation of the omnibus testing methods is that, in some cases, only statistical significance (i.e., p values) are provided as output; newer methods aimed at assigning more interpretable effect sizes are under development [96]. Finally, omnibus testing procedures by definition do not identify what variation in a microbial community might be associated with an outcome of interest. Thus, although they may require smaller sample sizes than per-feature tests to be well-powered, they provide less actionable information as a result. Nevertheless, omnibus tests are an important accompaniment to unsupervised visualization in providing a quantitative model in support of qualitative data exploration by ordination.

Integration of studies needs to address confounding effects that are unique to microbiome data

Meta-analyses of microbiome features are becoming more desirable and common, particularly when scaled to large human populations in order to achieve reliability and power for translational findings (Fig. 2e and f). Meta-analysis [91] is, in general, the quantitative integration of findings from multiple studies, and it is crucial in any molecular 'omics field for verifying true, biological associations and improving power. Meta-analyses of most types of microbiome data face major challenges because of strong, batch- and study-specific biases that arise in most stages of data generation (sample collection, DNA extraction, PCR amplification, sequencing, and bioinformatics [17, 104]). Previous multi-cohort studies have confirmed the driving effect of study-specific protocols on the clustering of sample-specific microbial profiles (i.e., on population structure discovery). In the absence of active efforts to normalize protocols among meta-analyzed studies, the effects of these batch differences may be surpassed in strength only by a few extreme microbial phenotypes (such as body site of origin) and can easily mask even strong biological factors such as antibiotics usage and disease subtype [105].

Changes in protocol can thus heavily influence both overall community configuration and the abundances of individual features [23], making analyses such as meta-analytic differential abundance tests challenging. This does not, of course, prevent sufficiently strong effects from being observed across studies (for example, in inflammatory bowel disease patients). Although such issues are generally acknowledged in the microbiome research community, efforts to address them have been limited to date. From an experimental design point of view, sharing among studies one or more 'mock communities', comprised of reference material and/or pre-determined collections of microbial strains in known proportions, can provide a reference for identifying and estimating sources of bias [106]. Likewise, the publication of negative control sequencing results in a consistent manner would allow background subtraction and contaminant identification among studies. However, such controls need to be incorporated during the early stages of a study and cannot be added in retrospect. They have the potential to make meta-analysis much easier when included. Mock communities can also be technically challenging to generate and, of course, incur additional costs during data generation, but they are likely to be of high value if included systematically in multiple studies within and across projects.

To enable true meta-analysis of microbial community surveys, quantitative protocols to adjust for batch- and study-specific effects must be developed. For population structure identification and adjustment, additional steps are necessary to correct for and reduce such effects before comparing and aggregating samples from different studies. Existing popular methods in RNA-seq whole-transcriptome profiling—such as ComBat [107] and limma [108]—may be potential candidates, though they should be modified to account for the zero-inflated and compositional (or count) nature of microbial abundances. For single-feature differential abundance analysis, study-specific effects may alternatively be addressed by adopting a unified model with identically defined effect sizes, which can then be compared and combined across studies using existing proper statistical methods (for example, mixed-effects models [86, 109]). Another promising direction is high-dimensional predictive modeling techniques (that is, using subjects’ microbial profiles as predictors for outcomes of interests), such as random forests, neural networks, and support vector machines, which are often successful in reproducibly predicting phenotype across multiple cohorts [91, 110]. The results obtained to date suggest avenues by which discriminative machine-learning models can be applied in microbial community settings to robustly associate features across multiple studies with outcomes of interest.


Like existing molecular epidemiology technologies, the translation of population studies of the human microbiome will require complex processes in order to achieve observational discovery, reproducibility across cohorts, and mechanistic validation (typically in models or in vitro). To date, a small number of studies have achieved this goal. For example, combining mouse models with a small cohort of 20 human subjects, Haiser and colleagues [111] built on decades of work linking Eggerthella lenta to inactivation of digoxin [112] to identify an operon that is expressed in a strain-specific manner in a subset of human microbiome carriers. As a further example, it has been shown that early-life exposure to distinct forms of taxon-specific lipopolysaccharide correlate with immune development and type 1 diabetes (T1D) risk, a result that was subsequently confirmed in mouse models (Box 5) [16]. Finally, in Clostridium difficile infection, models linking antibiotic exposure to bacterial species that are responsible for secondary bile acid synthesis in the gut have been successful in reducing recurrence [113]. In each of these cases, a combination of human population surveys with appropriate statistical modeling and mechanistic follow-up was able to identify specific bioactive microbes and, often, molecules. Further examples are emerging, particularly in the area of cancer immunotherapy, which can be dramatically modulated by the microbiome [114].

One of the outstanding gaps in translational population-scale microbiome studies is the lack of frameworks integrating host and microbiome functional properties at scale. For example, functional profiling of microbiome metagenomes and metatranscriptomes might be combined with cell-circuit reconstructions of immune cell subsets [115] and with electronic medical records for precision medicine. At the methodological level, few profiles of the microbiome have been carried out with scale and precision appropriate for advanced machine-learning tools such as causal inference and mediation analysis. Indeed, it is not yet clear which covariates should be collected to disambiguate cause from effect in the highly modifiable microbiome, particularly to facilitate risk-prediction models or clinical decision-making tools incorporating microbiome profiles. The microbiome has shown a remarkable combination of long-term persistence (e.g., strain retention for months or years [41, 116, 117]) with modifiability by a wide range of environmental factors (diet, pharmaceuticals, physical activity, age, and so on), making population structure and unobserved confounders a risk in large cohort studies.

Finally, human population studies provide a starting point for the follow-up characterization of microbial biochemical mechanisms, which can integrate characterization techniques such as culture-based physiology, microbial metabolism, co-culture, and interactions. Several of the most successful translational microbiome studies to date have—as in other areas of molecular epidemiology—begun with a population-level observation that was, eventually, traced back to one or more specific molecular mechanisms. In the case of the microbiome, this provides unique opportunities not only for prioritization of novel human drug targets, but also for the modulation of microbial activities by small molecules, diet or prebiotics, targeted probiotics, or engineered microbes or communities. To achieve these goals, studies of the microbiome must continue to refine the multiomic tools in the setting of population-scale epidemiology with rich study designs that can fully realize the therapeutic and diagnostic potential of the microbiome.



Ribosomal RNA


Single nucleotide variant


Type 1 diabetes


  1. 1.

    Kostic AD, Gevers D, Siljander H, Vatanen T, Hyotylainen T, Hamalainen AM, et al. The dynamics of the human infant gut microbiome in development and in progression toward type 1 diabetes. Cell Host Microbe. 2015;17:260–73.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  2. 2.


    CAS  PubMed  PubMed Central  Google Scholar 

  3. 3.


    CAS  PubMed  PubMed Central  Article  Google Scholar 

  4. 4.

    Luo C, Knight R, Siljander H, Knip M, Xavier RJ, Gevers D. ConStrains identifies microbial strains in metagenomic datasets. Nat Biotechnol. 2015;33:1045–52.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  5. 5.


    CAS  PubMed  PubMed Central  Article  Google Scholar 

  6. 6.


    CAS  PubMed  Article  Google Scholar 

  7. 7.


    Article  CAS  Google Scholar 

  8. 8.


    CAS  PubMed  PubMed Central  Article  Google Scholar 

  9. 9.


    CAS  PubMed  PubMed Central  Article  Google Scholar 

  10. 10.

    Verberkmoes NC, Russell AL, Shah M, Godzik A, Rosenquist M, Halfvarson J, et al. Shotgun metaproteomics of the human distal gut microbiota. ISME J. 2009;3:179–89.

    CAS  PubMed  Article  Google Scholar 

  11. 11.

    Zeevi D, Korem T, Zmora N, Israeli D, Rothschild D, Weinberger A, et al. Personalized nutrition by prediction of glycemic responses. Cell. 2015;163:1079–94.

    CAS  PubMed  Article  Google Scholar 

  12. 12.

    Suez J, Korem T, Zeevi D, Zilberman-Schapira G, Thaiss CA, Maza O, et al. Artificial sweeteners induce glucose intolerance by altering the gut microbiota. Nature. 2014;514:181–6.

    CAS  PubMed  Google Scholar 

  13. 13.

    Schirmer M, Smeekens SP, Vlamakis H, Jaeger M, Oosting M, Franzosa EA, et al. Linking the human gut microbiome to inflammatory cytokine production capacity. Cell. 2016;167:1125–36.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  14. 14.

    Zhernakova A, Kurilshikov A, Bonder MJ, Tigchelaar EF, Schirmer M, Vatanen T, et al. Population-based metagenomics analysis reveals markers for gut microbiome composition and diversity. Science. 2016;352:565–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  15. 15.

    Rooks MG, Garrett WS. Gut microbiota, metabolites and host immunity. Nat Rev Immunol. 2016;16:341–52.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  16. 16.

    Vatanen T, Kostic AD, d'Hennezel E, Siljander H, Franzosa EA, Yassour M, et al. Variation in microbiome LPS immunogenicity contributes to autoimmunity in humans. Cell. 2016;165:842–53.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  17. 17.

    Sinha R, Abnet CC, White O, Knight R, Huttenhower C. The microbiome quality control project: baseline study design and future directions. Genome Biol. 2015;16:276.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  18. 18.

    Falony G, Joossens M, Vieira-Silva S, Wang J, Darzi Y, Faust K, et al. Population-level analysis of gut microbiome variation. Science. 2016;352:560–4.

    CAS  PubMed  Article  Google Scholar 

  19. 19.

    Boutros PC. The path to routine use of genomic biomarkers in the cancer clinic. Genome Res. 2015;25:1508–13.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  20. 20.

    Ward LD, Kellis M. Interpreting noncoding genetic variation in complex traits and human disease. Nat Biotechnol. 2012;30:1095–106.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  21. 21.

    Hamady M, Knight R. Microbial community profiling for human microbiome projects: tools, techniques, and challenges. Genome Res. 2009;19:1141–52.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  22. 22.

    Underhill DM, Iliev ID. The mycobiota: interactions between commensal fungi and the host immune system. Nat Rev Immunol. 2014;14:405–16.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  23. 23.

    Brooks JP, Edwards DJ, Harwich Jr MD, Rivera MC, Fettweis JM, Serrano MG, et al. The truth about metagenomics: quantifying and counteracting bias in 16S rRNA studies. BMC Microbiol. 2015;15:66.

    PubMed  PubMed Central  Article  Google Scholar 

  24. 24.

    Franzosa EA, Morgan XC, Segata N, Waldron L, Reyes J, Earl AM, et al. Relating the metatranscriptome and metagenome of the human gut. Proc Natl Acad Sci U S A. 2014;111:E2329–38.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  25. 25.

    Booijink CC, Boekhorst J, Zoetendal EG, Smidt H, Kleerebezem M, de Vos WM. Metatranscriptome analysis of the human fecal microbiota reveals subject-specific expression profiles, with genes encoding proteins involved in carbohydrate metabolism being dominantly expressed. Appl Environ Microbiol. 2010;76:5533–40.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  26. 26.

    McHardy IH, Goudarzi M, Tong M, Ruegger PM, Schwager E, Weger JR, et al. Integrative analysis of the microbiome and metabolome of the human intestinal mucosal surface reveals exquisite inter-relationships. Microbiome. 2013;1:17.

    PubMed  PubMed Central  Article  Google Scholar 

  27. 27.

    Grassl N, Kulak NA, Pichler G, Geyer PE, Jung J, Schubert S, et al. Ultra-deep and quantitative saliva proteome reveals dynamics of the oral microbiome. Genome Med. 2016;8:44.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  28. 28.


    CAS  PubMed  PubMed Central  Article  Google Scholar 

  29. 29.


    CAS  PubMed  Article  Google Scholar 

  30. 30.

    Rajilic-Stojanovic M, de Vos WM. The first 1000 cultured species of the human gastrointestinal microbiota. FEMS Microbiol Rev. 2014;38:996–1047.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  31. 31.

    Scaldaferri F, Gerardi V, Mangiola F, Lopetuso LR, Pizzoferrato M, Petito V, et al. Role and mechanisms of action of Escherichia coli Nissle 1917 in the maintenance of remission in ulcerative colitis patients: an update. World J Gastroenterol. 2016;22:5505–11.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  32. 32.

    Kaas RS, Friis C, Ussery DW, Aarestrup FM. Estimating variation within the genes and inferring the phylogeny of 186 sequenced diverse Escherichia coli genomes. BMC Genomics. 2012;13:577.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  33. 33.

    Salipante SJ, Roach DJ, Kitzman JO, Snyder MW, Stackhouse B, Butler-Wu SM, et al. Large-scale genomic sequencing of extraintestinal pathogenic Escherichia coli strains. Genome Res. 2015;25:119–28.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  34. 34.

    Dobrindt U, Chowdary MG, Krumbholz G, Hacker J. Genome dynamics and its impact on evolution of Escherichia coli. Med Microbiol Immunol. 2010;199:145–54.

    CAS  PubMed  Article  Google Scholar 

  35. 35.

    Bosi E, Monk JM, Aziz RK, Fondi M, Nizet V, Palsson BO. Comparative genome-scale modelling of Staphylococcus aureus strains identifies strain-specific metabolic capabilities linked to pathogenicity. Proc Natl Acad Sci U S A. 2016;113:E3801–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  36. 36.


    CAS  PubMed  Article  Google Scholar 

  37. 37.


    CAS  PubMed  PubMed Central  Article  Google Scholar 

  38. 38.


    PubMed  PubMed Central  Article  CAS  Google Scholar 

  39. 39.

    Eren AM, Morrison HG, Lescault PJ, Reveillaud J, Vineis JH, Sogin ML. Minimum entropy decomposition: unsupervised oligotyping for sensitive partitioning of high-throughput marker gene sequences. ISME J. 2015;9:968–79.

    CAS  PubMed  Article  Google Scholar 

  40. 40.

    Tikhonov M, Leach RW, Wingreen NS. Interpreting 16S metagenomic data without clustering to achieve sub-OTU resolution. ISME J. 2015;9:68–80.

    PubMed  Article  CAS  Google Scholar 

  41. 41.

    Faith JJ, Guruge JL, Charbonneau M, Subramanian S, Seedorf H, Goodman AL, et al. The long-term stability of the human gut microbiota. Science. 2013;341:1237439.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  42. 42.

    Eren AM, Borisy GG, Huse SM, Mark Welch JL. Oligotyping analysis of the human oral microbiome. Proc Natl Acad Sci U S A. 2014;111:E2875–84.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  43. 43.


    CAS  PubMed  PubMed Central  Article  Google Scholar 

  44. 44.

    Edgar RC. UNOISE2: improved error-correction for Illumina 16S and ITS amplicon sequencing. bioRxiv. 2016; doi:

  45. 45.

    Amir A, McDonald D, Navas-Molina JA, Kopylova E, Morton JT, Zech Xu Z, et al. Deblur rapidly resolves single-nucleotide community sequence patterns. mSystems. 2017;2.

  46. 46.

    McInerney JO, McNally A, O'Connell MJ. Why prokaryotes have pangenomes. Nat Microbiol. 2017;2:17040.

    CAS  PubMed  Article  Google Scholar 

  47. 47.

    Loman NJ, Pallen MJ. Twenty years of bacterial genome sequencing. Nat Rev Microbiol. 2015;13:787–94.

    CAS  PubMed  Article  Google Scholar 

  48. 48.

    Konstantinidis KT, Ramette A, Tiedje JM. The bacterial species definition in the genomic era. Philos Trans R Soc Lond B Biol Sci. 2006;361:1929–40.

    PubMed  PubMed Central  Article  Google Scholar 

  49. 49.


    CAS  PubMed  PubMed Central  Article  Google Scholar 

  50. 50.

    Blazewicz SJ, Barnard RL, Daly RA, Firestone MK. Evaluating rRNA as an indicator of microbial activity in environmental communities: limitations and uses. ISME J. 2013;7:2061–8.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  51. 51.

    Franzosa EA, Hsu T, Sirota-Madi A, Shafquat A, Abu-Ali G, Morgan XC, Huttenhower C. Sequencing and beyond: integrating molecular 'omics' for microbial community profiling. Nat Rev Microbiol. 2015;13:360–72.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  52. 52.

    Conesa A, Madrigal P, Tarazona S, Gomez-Cabrero D, Cervera A, McPherson A, et al. A survey of best practices for RNA-seq data analysis. Genome Biol. 2016;17:13.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  53. 53.

    McMurdie PJ, Holmes S. Waste not, want not: why rarefying microbiome data is inadmissible. PLoS Comput Biol. 2014;10:e1003531.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  54. 54.


    PubMed  PubMed Central  Article  CAS  Google Scholar 

  55. 55.

    Stephen AM, Cummings JH. The microbial contribution to human faecal mass. J Med Microbiol. 1980;13:45–56.

    CAS  PubMed  Article  Google Scholar 

  56. 56.


    Article  CAS  Google Scholar 

  57. 57.

    Hang J, Desai V, Zavaljevski N, Yang Y, Lin X, Satya RV, et al. 16S rRNA gene pyrosequencing of reference and clinical samples and investigation of the temperature stability of microbiome profiles. Microbiome. 2014;2:31.

    PubMed  PubMed Central  Article  Google Scholar 

  58. 58.

    Song SJ, Amir A, Metcalf JL, Amato KR, Xu ZZ, Humphrey G, Knight R. Preservation methods differ in fecal microbiome stability, affecting suitability for field studies. mSystems 2016;1.

  59. 59.

    Faust K, Raes J. Microbial interactions: from networks to models. Nat Rev Microbiol. 2012;10:538–50.

    CAS  PubMed  Article  Google Scholar 

  60. 60.

    Mitri S, Foster KR. The genotypic view of social interactions in microbial communities. Annu Rev Genet. 2013;47:247–73.

    CAS  PubMed  Article  Google Scholar 

  61. 61.

    Tan J, Zuniga C, Zengler K. Unraveling interactions in microbial communities—from co-cultures to microbiomes. J Microbiol. 2015;53:295–305.

    PubMed  Article  Google Scholar 

  62. 62.

    Yu Z, Krause SM, Beck DA, Chistoserdova L. A synthetic ecology perspective: how well does behavior of model organisms in the laboratory predict microbial activities in natural habitats? Front Microbiol. 2016;7:946.

    PubMed  PubMed Central  Google Scholar 

  63. 63.

    Tsilimigras MC, Fodor AA. Compositional data analysis of the microbiome: fundamentals, tools, and challenges. Ann Epidemiol. 2016;26:330–5.

    PubMed  Article  Google Scholar 

  64. 64.

    Faust K, Sathirapongsasuti JF, Izard J, Segata N, Gevers D, Raes J, Huttenhower C. Microbial co-occurrence relationships in the human microbiome. PLoS Comput Biol. 2012;8:e1002606.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  65. 65.


    CAS  PubMed  PubMed Central  Article  Google Scholar 

  66. 66.

    Fang H, Huang C, Zhao H, Deng M. CCLasso: correlation inference for compositional data through Lasso. Bioinformatics. 2015;31:3172–80.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  67. 67.


    PubMed  PubMed Central  Article  CAS  Google Scholar 

  68. 68.


    CAS  PubMed  Article  Google Scholar 

  69. 69.


    CAS  PubMed  PubMed Central  Article  Google Scholar 

  70. 70.

    Org E, Blum Y, Kasela S, Mehrabian M, Kuusisto J, Kangas AJ, et al. Relationships between gut microbiota, plasma metabolites, and metabolic syndrome traits in the METSIM cohort. Genome Biol. 2017;18:70.

    PubMed  PubMed Central  Article  Google Scholar 

  71. 71.

    Ashrafian H, Li JV, Spagou K, Harling L, Masson P, Darzi A, et al. Bariatric surgery modulates circulating and cardiac metabolites. J Proteome Res. 2014;13:570–80.

    CAS  PubMed  Article  Google Scholar 

  72. 72.

    Wikoff WR, Anfora AT, Liu J, Schultz PG, Lesley SA, Peters EC, Siuzdak G. Metabolomics analysis reveals large effects of gut microflora on mammalian blood metabolites. Proc Natl Acad Sci U S A. 2009;106:3698–703.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  73. 73.

    Thorsen J, Brejnrod A, Mortensen M, Rasmussen MA, Stokholm J, Al-Soud WA, et al. Large-scale benchmarking reveals false discoveries and count transformation sensitivity in 16S rRNA gene amplicon data analysis methods used in microbiome studies. Microbiome. 2016;4:62.

    PubMed  PubMed Central  Article  Google Scholar 

  74. 74.

    Paulson JN, Stine OC, Bravo HC, Pop M. Differential abundance analysis for microbial marker-gene surveys. Nat Methods. 2013;10:1200–2.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  75. 75.


    CAS  PubMed  PubMed Central  Article  Google Scholar 

  76. 76.


    CAS  PubMed  Article  Google Scholar 

  77. 77.


    PubMed  PubMed Central  Article  CAS  Google Scholar 

  78. 78.


    PubMed  PubMed Central  Article  CAS  Google Scholar 

  79. 79.

    Jonsson V, Osterlund T, Nerman O, Kristiansson E. Statistical evaluation of methods for identification of differentially abundant genes in comparative metagenomics. BMC Genomics. 2016;17:78.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  80. 80.


    PubMed  PubMed Central  Article  Google Scholar 

  81. 81.

    White JR, Nagarajan N, Pop M. Statistical methods for detecting differentially abundant features in clinical metagenomic samples. PLoS Comput Biol. 2009;5:e1000352.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  82. 82.

    Mandal S, Van Treuren W, White RA, Eggesbo M, Knight R, Peddada SD. Analysis of composition of microbiomes: a novel method for studying microbial composition. Microb Ecol Health Dis. 2015;26:27663.

    PubMed  Google Scholar 

  83. 83.


    PubMed  PubMed Central  Article  CAS  Google Scholar 

  84. 84.


    CAS  PubMed  PubMed Central  Article  Google Scholar 

  85. 85.

    Chen EZ, Li H. A two-part mixed-effects model for analyzing longitudinal microbiome compositional data. Bioinformatics. 2016;32:2611–7.

    CAS  PubMed  Article  Google Scholar 

  86. 86.

    Zhang X, Mallick H, Tang Z, Zhang L, Cui X, Benson AK, Yi N. Negative binomial mixed models for analyzing microbiome count data. BMC Bioinformatics. 2017;18:4.

    PubMed  PubMed Central  Article  Google Scholar 

  87. 87.

    Bucci V, Tzen B, Li N, Simmons M, Tanoue T, Bogart E, et al. MDSINE: Microbial Dynamical Systems INference Engine for microbiome time-series analyses. Genome Biol. 2016;17:121.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  88. 88.


    CAS  PubMed  PubMed Central  Article  Google Scholar 

  89. 89.

    Koren O, Knights D, Gonzalez A, Waldron L, Segata N, Knight R, et al. A guide to enterotypes across the human body: meta-analysis of microbial community structures in human microbiome datasets. PLoS Comput Biol. 2013;9:e1002863.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  90. 90.

    Gajer P, Brotman RM, Bai G, Sakamoto J, Schutte UM, Zhong X, et al. Temporal dynamics of the human vaginal microbiota. Sci Transl Med. 2012;4:132ra152.

    Article  Google Scholar 

  91. 91.

    Pasolli E, Truong DT, Malik F, Waldron L, Segata N. Machine learning meta-analysis of large metagenomic datasets: tools and biological insights. PLoS Comput Biol. 2016;12:e1004977.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  92. 92.

    Anderson MJ. A new method for non‐parametric multivariate analysis of variance. Austral Ecol. 2001;26:32–46.

    Google Scholar 

  93. 93.

    Zhao N, Chen J, Carroll IM, Ringel-Kulka T, Epstein MP, Zhou H, et al. Testing in microbiome-profiling studies with MiRKAT, the Microbiome Regression-Based Kernel Association Test. Am J Hum Genet. 2015;96:797–807.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  94. 94.

    Clarke KR. Non‐parametric multivariate analyses of changes in community structure. Aust J Ecol. 1993;18:117–43.

    Article  Google Scholar 

  95. 95.

    Tang ZZ, Chen G, Alekseyenko AV. PERMANOVA-S: association test for microbial community composition that accommodates confounders and multiple distances. Bioinformatics. 2016;32:2618–25.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  96. 96.

    Randolph TW, Zhao S, Copeland W, Hullar M, Shojaie A. Kernel-Penalized regression for analysis of microbiome data. arXiv 2015;arXiv:151100297.

  97. 97.

    Jonsson V, Osterlund T, Nerman O, Kristiansson E. Variability in metagenomic count data and its influence on the identification of differentially abundant genes. J Comput Biol. 2017;24:311–26.

    CAS  PubMed  Article  Google Scholar 

  98. 98.

    Weiss S, Xu ZZ, Peddada S, Amir A, Bittinger K, Gonzalez A, et al. Normalization and microbial differential abundance strategies depend upon data characteristics. Microbiome. 2017;5:27.

    PubMed  PubMed Central  Article  Google Scholar 

  99. 99.

    Neter J, Kutner MH, Nachtsheim CJ, Wasserman W. Applied linear statistical models. Chicago: Irwin; 1996.

    Google Scholar 

  100. 100.

    Olive DJ. Multivariate linear regression. In: Olive DJ, editor. Linear regression. Cham: Springer; 2017. p. 343–87.

    Google Scholar 

  101. 101.

    Hidalgo B, Goodman M. Multivariate or multivariable regression? Am J Public Health. 2013;103:39–40.

    PubMed  PubMed Central  Article  Google Scholar 

  102. 102.

    Tsai AC. Achieving consensus on terminology describing multivariable analyses. Am J Public Health. 2013;103:e1.

    PubMed  PubMed Central  Article  Google Scholar 

  103. 103.

    Fang R, Wagner BD, Harris JK, Fillon SA. Zero-inflated negative binomial mixed model: an application to two microbial organisms important in oesophagitis. Epidemiol Infect. 2016;144:2447–55.

    CAS  PubMed  Article  Google Scholar 

  104. 104.

    Brooks JP. Challenges for case-control studies with microbiome data. Ann Epidemiol. 2016;26:336–41.

    PubMed  Article  Google Scholar 

  105. 105.

    Lozupone CA, Stombaugh J, Gonzalez A, Ackermann G, Wendel D, Vazquez-Baeza Y, et al. Meta-analyses of studies of the human microbiota. Genome Res. 2013;23:1704–14.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  106. 106.

    Bokulich NA, Rideout JR, Mercurio WG, Shiffer A, Wolfe B, Maurice CF, et al. mockrobiota: a public resource for microbiome bioinformatics benchmarking. mSystems 2016;1.

  107. 107.


    PubMed  Article  Google Scholar 

  108. 108.


    PubMed  PubMed Central  Article  CAS  Google Scholar 

  109. 109.

    Romero R, Hassan SS, Gajer P, Tarca AL, Fadrosh DW, Nikita L, et al. The composition and stability of the vaginal microbiota of normal pregnant women is different from that of non-pregnant women. Microbiome. 2014;2:4.

    PubMed  PubMed Central  Article  Google Scholar 

  110. 110.

    Ananthakrishnan AN, Luo C, Yajnik V, Khalili H, Garber JJ, Stevens BW, et al. Gut microbiome function predicts response to anti-integrin biologic therapy in inflammatory bowel diseases. Cell Host Microbe. 2017;21:603–10.

    CAS  PubMed  Article  Google Scholar 

  111. 111.


    CAS  PubMed  PubMed Central  Article  Google Scholar 

  112. 112.

    Saha JR, Butler Jr VP, Neu HC, Lindenbaum J. Digoxin-inactivating bacteria: identification in human gut flora. Science. 1983;220:325–7.

    CAS  PubMed  Article  Google Scholar 

  113. 113.


    CAS  PubMed  Article  Google Scholar 

  114. 114.

    Vetizou M, Pitt JM, Daillere R, Lepage P, Waldschmitt N, Flament C, et al. Anticancer immunotherapy by CTLA-4 blockade relies on the gut microbiota. Science. 2015;350:1079–84.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  115. 115.

    Yosef N, Regev A. Writ large: genomic dissection of the effect of cellular environment on immune response. Science. 2016;354:64–8.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  116. 116.

    Schloissnig S, Arumugam M, Sunagawa S, Mitreva M, Tap J, Zhu A, et al. Genomic variation landscape of the human gut microbiome. Nature. 2013;493:45–50.

    PubMed  Article  CAS  Google Scholar 

  117. 117.

    Franzosa EA, Huang K, Meadow JF, Gevers D, Lemon KP, Bohannan BJ, Huttenhower C. Identifying personal microbiomes using metagenomic codes. Proc Natl Acad Sci U S A. 2015;112:E2930–8.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  118. 118.

    Peet A, Kool P, Ilonen J, Knip M, Tillmann V, Group DS. Birth weight in newborn infants with different diabetes-associated HLA genotypes in three neighbouring countries: Finland, Estonia and Russian Karelia. Diabetes Metab Res Rev. 2012;28:455–61.

    CAS  PubMed  Article  Google Scholar 

  119. 119.

    Franzen O, Hu J, Bao X, Itzkowitz SH, Peter I, Bashir A. Improved OTU-picking using long-read 16S rRNA gene amplicon sequencing and generic hierarchical clustering. Microbiome. 2015;3:43.

    PubMed  PubMed Central  Article  Google Scholar 

  120. 120.


    CAS  PubMed  PubMed Central  Article  Google Scholar 

  121. 121.

    Tu Q, He Z, Zhou J. Strain/species identification in metagenomes using genome-specific markers. Nucleic Acids Res. 2014;42:e67.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  122. 122.


    PubMed  PubMed Central  Article  Google Scholar 

  123. 123.

    Ahn TH, Chai J, Pan C. Sigma: strain-level inference of genomes from metagenomic analysis for biosurveillance. Bioinformatics. 2015;31:170–7.

    CAS  PubMed  Article  Google Scholar 

  124. 124.

    Francis OE, Bendall M, Manimaran S, Hong C, Clement NL, Castro-Nallar E, et al. Pathoscope: species identification and strain attribution with unassembled sequencing data. Genome Res. 2013;23:1721–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  125. 125.

    Cleary B, Brito IL, Huang K, Gevers D, Shea T, Young S, Alm EJ. Detection of low-abundance bacterial strains in metagenomic datasets by eigengenome partitioning. Nat Biotechnol. 2015;33:1053–60.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  126. 126.


    CAS  PubMed  PubMed Central  Article  Google Scholar 

  127. 127.

    McClure R, Balasubramanian D, Sun Y, Bobrovskyy M, Sumby P, Genco CA, et al. Computational analysis of bacterial RNA-Seq data. Nucleic Acids Res. 2013;41:e140.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  128. 128.

    Ghosh S, Chan CK. Analysis of RNA-Seq data using TopHat and Cufflinks. Methods Mol Biol. 2016;1374:339–61.

    CAS  PubMed  Article  Google Scholar 

  129. 129.

    Narayanasamy S, Jarosz Y, Muller EE, Heintz-Buschart A, Herold M, Kaysen A, et al. IMP: a pipeline for reproducible reference-independent integrated metagenomic and metatranscriptomic analyses. Genome Biol. 2016;17:260.

    PubMed  PubMed Central  Article  Google Scholar 

  130. 130.

    Westreich ST, Korf I, Mills DA, Lemay DG. SAMSA: a comprehensive metatranscriptome analysis pipeline. BMC Bioinformatics. 2016;17:399.

    PubMed  PubMed Central  Article  Google Scholar 

  131. 131.

    Ni Y, Li J, Panagiotou G. COMAN: a web server for comprehensive metatranscriptomics analysis. BMC Genomics. 2016;17:622.

    PubMed  PubMed Central  Article  Google Scholar 

  132. 132.

    Leung HC, Yiu SM, Parkinson J, Chin FY. IDBA-MT: de novo assembler for metatranscriptomic data generated from next-generation sequencing technology. J Comput Biol. 2013;20:540–50.

    CAS  PubMed  Article  Google Scholar 

  133. 133.


    CAS  PubMed  PubMed Central  Article  Google Scholar 

  134. 134.

    Bose T, Haque MM, Reddy C, Mande SS. COGNIZER: a framework for functional annotation of metagenomic datasets. PLoS One. 2015;10:e0142102.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  135. 135.

    Kim J, Kim MS, Koh AY, Xie Y, Zhan X. FMAP: Functional Mapping and Analysis Pipeline for metagenomics and metatranscriptomics studies. BMC Bioinformatics. 2016;17:420.

    PubMed  PubMed Central  Article  Google Scholar 

  136. 136.


    PubMed  PubMed Central  Article  CAS  Google Scholar 

  137. 137.

    Nayfach S, Bradley PH, Wyman SK, Laurent TJ, Williams A, Eisen JA, et al. Automated and accurate estimation of gene family abundance from shotgun metagenomes. PLoS Comput Biol. 2015;11:e1004573.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  138. 138.

    Morgan XC, Kabakchiev B, Waldron L, Tyler AD, Tickle TL, Milgrom R, et al. Associations between host gene expression, the mucosal microbiome, and clinical outcome in the pelvic pouch of patients with inflammatory bowel disease. Genome Biol. 2015;16:67.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  139. 139.

    Ban Y, An L, Jiang H. Investigating microbial co-occurrence patterns based on metagenomic compositional data. Bioinformatics. 2015;31:3322–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  140. 140.

    Deng Y, Jiang YH, Yang Y, He Z, Luo F, Zhou J. Molecular ecological network analyses. BMC Bioinformatics. 2012;13:113.

    PubMed  PubMed Central  Article  Google Scholar 

  141. 141.

    Biswas S, McDonald M, Lundberg DS, Dangl JL, Jojic V. Learning microbial interaction networks from metagenomic count data. J Comput Biol. 2016;23:526–35.

    CAS  PubMed  Article  Google Scholar 

  142. 142.

    Shaw GT, Pao YY, Wang D. MetaMIS: a metagenomic microbial interaction simulator based on microbial community profiles. BMC Bioinformatics. 2016;17:488.

    PubMed  PubMed Central  Article  Google Scholar 

  143. 143.

    Shafiei M, Dunn KA, Boon E, MacDonald SM, Walsh DA, Gu H, Bielawski JP. BioMiCo: a supervised Bayesian model for inference of microbial community structure. Microbiome. 2015;3:8.

    PubMed  PubMed Central  Article  Google Scholar 

  144. 144.

    Shafiei M, Dunn KA, Chipman H, Gu H, Bielawski JP. BiomeNet: a Bayesian model for inference of metabolic divergence among microbial communities. PLoS Comput Biol. 2014;10:e1003918.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  145. 145.

    McGeachie MJ, Sordillo JE, Gibson T, Weinstock GM, Liu YY, Gold DR, et al. Longitudinal prediction of the infant gut microbiome with dynamic Bayesian networks. Sci Rep. 2016;6:20359.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

Download references


This research was supported by grants NIH U54DE023798, R01HG005220, and the STARR Cancer Consortium.

Author information




All authors wrote, read, and approved the final manuscript.

Corresponding author

Correspondence to Curtis Huttenhower.

Ethics declarations

Competing interests

The authors declare that they have no competing interests.

Publisher's Note


Rights and permissions


Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Mallick, H., Ma, S., Franzosa, E.A. et al. Experimental design and quantitative analysis of microbial community multiomics. Genome Biol 18, 228 (2017).

Download citation


  • Microbiome Profiles
  • Mock Community
  • Microbiome Data
  • Metatranscriptomics
  • Differential Abundance Analysis