A predictive index for health status using species-level gut microbiome profiling
Taxonomic signatures of cause-specific mortality risk in human gut microbiome
Aaro Salosensaari, Ville Laitinen, … Teemu Niiranen
Gut metagenome associations with extensive digital health data in a volunteer-based Estonian microbiome cohort
Oliver Aasmets, Kertu Liis Krigul, … Elin Org
Health and disease markers correlate with gut microbiome composition across thousands of people
Ohad Manor, Chengzhen L. Dai, … Andrew T. Magis
The gut microbiome-metabolome dataset collection: a curated resource for integrative meta-analysis
Efrat Muller, Yadid M. Algavi & Elhanan Borenstein
Microbiome connections with host metabolism and habitual diet from 1,098 deeply phenotyped individuals
Francesco Asnicar, Sarah E. Berry, … Nicola Segata
Blood metabolome predicts gut microbiome α-diversity in humans
Tomasz Wilmanski, Noa Rappaport, … Nathan D. Price
Toward an improved definition of a healthy microbiome for healthy aging
Tarini Shankar Ghosh, Fergus Shanahan & Paul W. O'Toole
Temporal variability in quantitative human gut microbiome profiles and implications for clinical research
Doris Vandeputte, Lindsey De Commer, … Jeroen Raes
Gut microbiome, big data and machine learning to promote precision medicine for cancer
Giovanni Cammarota, Gianluca Ianiro, … Giampaolo Tortora
Vinod K. Gupta1,2,
Minsuk Kim ORCID: orcid.org/0000-0001-9958-03801,2,
Utpal Bakshi ORCID: orcid.org/0000-0001-8421-90311,2,
Kevin Y. Cunningham3,4,
John M. Davis III5,
Konstantinos N. Lazaridis ORCID: orcid.org/0000-0002-0437-681X6,
Heidi Nelson7,
Nicholas Chia ORCID: orcid.org/0000-0001-9652-691X1,2 &
Jaeyun Sung ORCID: orcid.org/0000-0002-3005-27981,2,5
Providing insight into one's health status from a gut microbiome sample is an important clinical goal in current human microbiome research. Herein, we introduce the Gut Microbiome Health Index (GMHI), a biologically-interpretable mathematical formula for predicting the likelihood of disease independent of the clinical diagnosis. GMHI is formulated upon 50 microbial species associated with healthy gut ecosystems. These species are identified through a multi-study, integrative analysis on 4347 human stool metagenomes from 34 published studies across healthy and 12 different nonhealthy conditions, i.e., disease or abnormal bodyweight. When demonstrated on our population-scale meta-dataset, GMHI is the most robust and consistent predictor of disease presence (or absence) compared to α-diversity indices. Validation on 679 samples from 9 additional studies results in a balanced accuracy of 73.7% in distinguishing healthy from non-healthy groups. Our findings suggest that gut taxonomic signatures can predict health status, and highlight how data sharing efforts can provide broadly applicable discoveries.
Recent advances in the field of human gut microbiome research have revealed significant associations and potential mechanistic insights regarding a vast array of complex, chronic diseases, including cancer1,2, auto-immune disease3,4,5, and metabolic syndrome6,7,8. Undoubtedly, the many microbiome studies that focused on disease contexts have been essential for elucidating underlying pathophysiological mechanisms, and for developing potential intervention strategies. As researchers uncover more details regarding which gut commensals may play a significant role in host health and disease, a promising translational application of this knowledge would be towards developing analytical tests or quantitative methods that provide indication of one's health based upon a gut microbiome snapshot9,10,11,12.
The creation of algorithm-driven markers for detecting early signs of disease prior to the occurrence of specific, diagnosable symptoms, especially from biospecimens that can be collected regularly and noninvasively, is an exciting avenue forward for personalized medicine13,14,15,16. To this end, current gut microbiome science could play a significant role in the development of stool-based tests for dynamically monitoring and predicting wellness. We can even imagine a hypothetical scenario wherein one can, through continuous monitoring, be able to detect significant changes or abnormalities in comparison to her/his normal baseline measurements; in turn, this will be the cue for additional tests or lifestyle interventions. In this sense, a proof of principle of collecting diverse longitudinal biomolecular data from human subjects, and translating the corresponding complex datasets into actionable possibilities, has been previously demonstrated by Price et al.17
Since new insights regarding the gut microbiome and human health need to be reliable and robust across a wide range of human subjects and conditions, population-level analyses can serve as an important platform for the discovery of broadly applicable principles and methodologies18,19,20,21,22. Traditionally, the collection of a substantial amount of microbiome samples (on the order of hundreds to thousands) for large cohort investigations has been undertaken at well-funded, major research centers and consortiums, mainly due to the prohibitive costs and/or lack of infrastructure for lone laboratories. However, with the recent progress in the call for broad data sharing policies and practices23, conducting analyses by crowdsourcing or repurposing data from existing published studies have already begun to play important roles in either hypothesis generation or validation in microbiome research1,24,25,26,27. As much 16s rRNA gene amplicon and shotgun metagenomic sequencing data are now readily available from multiple, independent studies conducted around the world, integration by pooling these datasets would provide a promising strategy to study health- and disease-associated signatures at large scale, as well as to gain new, holistic insights not offered by smaller, individual studies.
In this study, we introduce the Gut Microbiome Health Index (GMHI), a robust index for evaluating health status (i.e., degree of presence/absence of diagnosed disease) based on the species-level taxonomic profile of a stool shotgun metagenome (gut microbiome) sample. GMHI determines the likelihood of having a disease, independent of the clinical diagnosis; this is done so by comparing the relative abundances of two sets of microbial species associated with good and adverse health conditions, which are identified from an integrated dataset of 4347 publicly available, human stool metagenomes pooled across multiple studies encompassing various disease states. By applying GMHI to each sample in our population-scale meta-dataset, we found that GMHI distinguishes healthy from nonhealthy groups far better than ecological indices (e.g., Shannon diversity and richness) generally considered as markers for gut health and dysbiosis. Intra-study comparisons of stool metagenomes between healthy and nonhealthy phenotypes demonstrate that GMHI is the most robust and consistent predictor of health. Finally, to confirm that GMHI classification accuracy was not a result of over-fitting on the discovery cohort, we test our approach on a validation set of 679 samples from eight additional published studies and one new cohort (this study). We find that GMHI not only demonstrates strong reproducibility in stratifying healthy and nonhealthy groups, but also outperforms α-diversity indices.
A meta-dataset of integrated human stool metagenomes
An overview of our multi-study integration approach, wherein we acquired 4347 raw shotgun stool metagenomes (2636 and 1711 metagenomes from healthy and nonhealthy individuals, respectively) from 34 independently published studies, is depicted in Fig. 1a. In this study, "healthy" subjects were defined as those who were reported as not having any overt disease nor adverse symptoms at the time of the original study; alternatively, "nonhealthy" subjects were defined as those who were clinically diagnosed with a specific disease, or determined to have abnormal bodyweight based on body mass index (BMI). Accordingly, 1711 stool metagenomes from patients across 12 different disease or abnormal bodyweight conditions were pooled together into a single aggregate nonhealthy group. (Our sample selection criteria are described in "Methods." Importantly, all metagenomes were reprocessed uniformly, thereby removing a major nonbiological source of variance among different studies, as previously demonstrated28.) A description of the studies whose human stool metagenomes were collected and processed through our computational pipeline is provided in Table 1. We note that, in order to eventually identify features of the gut microbiome associated exclusively with health, it is important to be disease agnostic by considering a broad range of nonhealthy phenotypes. We provide all subjects' phenotype, age, sex, BMI, and other questionnaire measures (as provided in their respective original study) in Supplementary Data 1. Along with the additional 679 stool metagenome samples used for validation purposes (discussed below), this study provides the largest metagenomic (pooled) analysis of the human gut microbiome to date, in regards to the number of samples, phenotypes, and studies.
Fig. 1: Multi-study integration of human stool metagenomes leads to a meta-dataset of healthy and nonhealthy gut microbiomes.
a Schematic overview. A survey was conducted in PubMed and Google Scholar to search for published studies with publicly available human stool metagenome (gut microbiome) samples from healthy and nonhealthy individuals. The initial collection of stool metagenomes consisted of 7589 samples from 55 independent studies. All samples (.fastq files) were downloaded and reprocessed uniformly using identical bioinformatics methods. After quality control of sequenced reads, species-level taxonomic profiling was then performed. Studies and metagenome samples were removed based on several exclusion criterias. Finally, a total of 4347 samples (2636 and 1711 metagenomes from healthy and nonhealthy individuals, respectively) from 34 studies ranging across healthy and 12 nonhealthy phenotypes were assembled into a meta-dataset for downstream analyses. b Distribution of microbial species' prevalence across the 4347 stool metagenome samples in the meta-dataset. After removing viruses, unknown/unclassified species-level entities, and rarely observed species (i.e., detected <1% of all samples), 313 species remained for further analyses. c Principal coordinates analysis (PCoA) ordination plot based on Bray–Curtis distances shows that healthy (blue; n = 2636) and nonhealthy (orange; n = 1711) groups have significantly different distributions of gut microbiome profiles according to PERMANOVA (R2 = 0.017, P < 0.001) after adjusting for each sample's study origin. Each point corresponds to a sample. Ellipses correspond to 95% confidence regions. d In an identical PCoA plot, each color represents one of the 13 different phenotypes of health or disease. Among- and within-group dissimilarities differ only weakly (ANOSIM R = 0.21, P = 0.001).
Table 1 Human stool metagenome datasets analyzed in this study.
Importantly, we chose to integrate datasets from independent studies for two notable advantages: (i) the expansion of sample number could help to amplify the primary biological signal of interest and improve statistical power29,30; and (ii) the identified health/disease-associated signatures could encompass a wide range of heterogeneity across different sources and conditions (e.g., host genetics, geography, dietary and lifestyle patterns, age, sex, birth mode, early life exposures, medication history), thereby helping to identify robust findings despite systematic biases from batch effects or other confounding factors28,31.
After downloading, reprocessing, and performing quality filtration on all raw metagenomes, species-level taxonomic profiling was carried out using the MetaPhlAn2 pipeline32 ("Methods"). Of note, our study was mainly conducted upon species-level taxonomy information to obtain as much precise and comprehensive information about the gut microbiome as possible. A total of 1201 species were detected in at least one metagenome sample; after removing viruses, and species that were rarely observed or of unknown/unclassified identity ("Methods"), 313 species remained for further analysis (Fig. 1b and Supplementary Data 2; a phylogenetic tree showing the evolutionary relationships among these species is shown in Supplementary Fig. 1). Interestingly, six species (Bacteroides ovatus, Bacteroides uniformis, Bacteroides vulgatus, Faecalibacterium prausnitzii, Ruminococcus obeum, and Ruminococcus torques) were of high prevalence (i.e., detected in >90% of all 4347 samples).
Healthy and nonhealthy guts show species-level differences
The overall ecology of the gut microbiome has often been associated with host health8,33,34,35. Using species-level relative abundance (i.e., proportion) profiles, we examined the differences in gut microbial diversity between the healthy and nonhealthy groups. First, when using principal coordinate analysis (PCoA) ordination, we identified a significant difference between the distributions of these two groups (permutational multivariate analysis of variance (PERMANOVA), R2 = 0.02, P < 0.001; Fig. 1c). In the same PCoA plot in which the healthy and 12 nonhealthy phenotypes were presented simultaneously (Fig. 1d), we found only a weak difference among groups (analysis of similarities (ANOSIM) R = 0.21, P = 0.001).
Design rationale for a GMHI
We envision that the most intuitive way to determine how closely one's microbiome resembles that of a healthy (or nonhealthy) population is to quantify the balance between health-associated microbes relative to disease-associated microbes. Therefore, we propose an index in the form of a rational equation (and thereby yielding a dimensionless quantity) between two sets of microbial species: those that are more frequently observed in healthy compared to nonhealthy groups vs. those that are less frequently observed in healthy compared to nonhealthy groups. Next, we use our compendium of publicly available datasets, which were derived from healthy and nonhealthy human subjects, to identify these two sets of species. Finally, with these species, we tune the parameters of a predefined formula, as well as evaluate its classification accuracy. The logical rationale of each major step during the development, demonstration, and validation of our index for predicting general health status (i.e., presence/absence of diagnosed disease) from a gut microbiome sample is detailed below. In addition, a step-by-step protocol is provided in "Methods."
A prevalence-based strategy to identify health-associated microbes
We set out to identify distinct microbial species associated with healthy (H) and nonhealthy (N) groups. Here, we use a prevalence-based strategy to deal with the sparse nature of microbiome datasets. For this, we first determine pH,m and pN,m, or the prevalence of microbial species m in H and N, respectively (prevalence is defined as the proportion of samples in a given group wherein m is considered "present," i.e., relative abundance ≥ 1.0 × 10−5.) Next, for comparing the two prevalences in H and N, we apply the following two criteria: prevalence fold change \(f_m^{H,N}\) and prevalence difference \(d_m^{H,N}\), defined as \(\frac{{p_{H,m}}}{{p_{N,m}}}\) and pH,m − pN,m, respectively. A significant effect size between the two prevalences is considered to exist if both criteria satisfy (predetermined) minimum thresholds for prevalence fold change θf and prevalence difference θd (how we determine the best pair of thresholds is described below). For all detectable microbial species that simultaneously satisfy \(f_m^{H,N} \ge \theta _f\) and \(d_m^{H,N} \ge \theta _d\), we term these species observed more frequently in H (than in N) as "health-prevalent" species MH. Analogously, we identify "health-scarce" species MN, or the species observed less frequently in H (than in N), as those that satisfy \(f_m^{N,H} \ge \theta _f\) and \(d_m^{N,H} \ge \theta _d\), where \(f_m^{N,H}\) and \(d_m^{N,H}\) is defined as \(\frac{{p_{N,m}}}{{p_{H,m}}}\) and pN,m − pH,m, respectively. In this regard, the species that are eventually chosen to compose MH and MN are both dependent on θf and θd. An important strength of our prevalence-based strategy for identifying microbial associations is that it does not calculate or compare averages of measurements taken from various sources, which is challenging to justify when biological and technical heterogeneity could vary greatly across independent studies. Rather, our approach compares frequencies of a signal—on a sample-by-sample basis—between two groups, and represents a strategy more applicable to the context of integrating high-throughput data from different studies. Importantly, we chose to simultaneously test two thresholds, rather than one, in order to increase our confidence in the robustness of MH and MN, as well as to overcome biases that can occur from using only one type of threshold.
Collective abundances of two sets of microbial taxonomies
Having a strategy to identify microbial species associated with healthy (i.e., health-prevalent species MH) and nonhealthy (i.e., health-scarce species MN), we next couple these two species sets with a computational procedure that quantifies the presence/absence of diagnosed disease for any gut microbiome sample. To this end, we developed the following mathematical formula: for species of MH in sample i, their "collective abundance" \(\psi _{M_H,i}\) is defined as
$$\psi _{M_H,i} = \frac{{R_{M_H,i}}}{{\left| {M_H} \right|}}\mathop {\sum}\limits_{j \in I_{M_H}} {\left| {n_{j,i}{\mathrm{ln}}\left( {n_{j,i}} \right)} \right|},$$
where \(R_{M_H,i}\) is the richness of MH species in sample i, |MH| is the set size of MH, \(I_{M_H}\) is the index set of MH, and nj,i is the relative abundance of species j in sample i. In brief, \(\psi _{M_H,i}\) is the product of the (i) richness, that is, the numeric count of "present" taxonomies, of MH species; and (ii) the geometric mean of their relative abundances. Full details on the physical meaning and derivation of \(\psi _{M_H,i}\) are described in "Methods." For the species of MN in the same sample i, their "collective abundance" \(\psi _{M_N,i}\) can be defined analogously. Next, the collective abundances of species in sets MH and MN in sample i are compared using the ratio of \(\psi _{M_H,i}\) to \(\psi _{M_N,i}\) as
$$h_{i,M_H,M_N} = {\mathrm{log}}_{10}\left( {\frac{{\psi _{M_H,i}}}{{\psi _{M_N,i}}}} \right),$$
where \(h_{i,M_H,M_N}\) denotes the degree to which sample i portrays the collective abundance of MH to that of MN. More specifically, a positive or negative \(h_{i,M_H,M_N}\) suggests that sample i is characterized more by the microbes of MH or MN, respectively; an \(h_{i,M_H,M_N}\) of 0 indicates that there is an equal balance of both species sets. Full details on the derivation of \(h_{i,M_H,M_N}\) are provided in "Methods."
Determining health-prevalent and health-scarce species
The minimum thresholds θf and θd for prevalence fold change and prevalence difference, respectively, are used to control for the number of health-prevalent species MH and health-scarce species MN; species that simultaneously satisfy the two types of thresholds are selected to be included in one of either group. Afterwards, MH and MN is provided as input features for \(\psi _{M_H,i}\) and \(\psi _{M_N,i}\), respectively, and for the calculation of \(h_{i,M_H,M_N}\), which in turn can classify stool metagenome sample i as healthy (i.e., \(h_{i,M_H,M_N} > 0\)), nonhealthy (i.e., \(h_{i,M_H,M_N} < 0\)), or neither (i.e., \(h_{i,M_H,M_N} = 0\)). Lastly, \(h_{i,M_H,M_N}\) is tested on all 4347 stool metagenomes in our meta-dataset to find the balanced accuracy \(\chi _{M_H,M_N}\), that is, an average of the proportions of healthy and nonhealthy samples that were correctly classified, or
$$\chi _{M_H,M_N} = \frac{{P\left( {h_{i,M_H,M_N} > 0|i \in H} \right) + P\left( {h_{i,M_H,M_N} < 0|i \in N} \right)}}{2},$$
where \(P( {h_{i,M_H,M_N} > 0|i \in H} )\) is the proportion of samples in the healthy group (H) whose \(h_{i,M_H,M_N}\)s are positive, and \(P( {h_{i,M_H,M_N} < 0|i \in N})\) is the proportion of samples in the nonhealthy group (N) whose \(h_{i,M_H,M_N}\)s are negative.
We determine the final, optimal sets of health-prevalent and health-scarce species (and their corresponding prevalence thresholds) as those that result in the highest balanced accuracy \(\chi _{M_H,M_N}^{\mathrm{max}}\). This was done accordingly: after systematically testing across a range of two different thresholds (every pair of θf and θd gives different sets of MH and MN, and in turn, a different \(\chi _{M_H,M_N}\)), we found \(\chi _{M_H,M_N}^{\mathrm{max}}\) to be 69.7% when θf and θd were set to 1.4 and 10%, respectively (Supplementary Table 1). When applying the same approach for abundance profiles of all other taxonomic ranks, as well as of MetaCyc pathways, the highest accuracies found in these were as follows: Phylum, 42.1%; Class, 60.1%; Order, 62.4%; Family, 67.2%; Genus, 68.2%; and MetaCyc pathway, 59.4% (Supplementary Table 2). As evidenced by these results, taxonomic species shows the best classification accuracy. In addition, performing our method in 10-fold cross-validation using species-level abundances resulted in an accuracy of 69.6% (Supplementary Table 3), which is nearly identical to the balanced accuracy of 69.7% achieved by testing on the set of samples from which the classifier was derived. Lastly, in Supplementary Fig. 2, we show a sensitivity analysis of how the balanced accuracy χ changes with respect to the species' prevalence thresholds θf and θd.
We identified 50 microbial species that satisfy both of the aforementioned thresholds for the highest balanced accuracy; among these 50 species, 7 and 43 comprise the health-prevalent and health-scarce groups, respectively (Table 2). Interestingly, we found higher relative abundance levels of health-prevalent and health-scarce species in the healthy and nonhealthy group, respectively (Supplementary Fig. 3). Furthermore, in Supplementary Fig. 4, we show the prevalence of these species in case (i.e., nonhealthy) and/or control (i.e., healthy) for the 34 published studies upon which our stool metagenome meta-dataset was derived. Despite the heterogeneity and unevenness in prevalences across all studies, we found that, by and large, health-prevalent and health-scarce species were observed more frequently in the control and case samples, respectively. We report known associations between health-prevalent/-scarce species and health/disease in Supplementary Data 3.
Table 2 Microbial species of the Health-prevalent and Health-scarce groups.
Henceforth, we term the ratio \(h_{i,M_H,M_N}\) between these two groups of 7 Health-prevalent and 43 Health-scarce species as the GMHI. GMHI is a dimensionless metric designed to simplify the accumulation of Health-prevalent and Health-scarce species observed to be present in a microbiome sample. In practice, GMHI indicates the degree to which a subject's stool metagenome sample portrays microbial taxonomies associated with either healthy or nonhealthy. An additional summary of the design process of GMHI, enumerated in three key steps, is provided in Supplementary Note 1.
Analogous to the example mentioned above, a positive or negative GMHI allows the sample to be classified as healthy or nonhealthy, respectively; a GMHI of 0 indicates an equal balance of Health-prevalent and Health-scarce species, and thereby classified as neither. Therefore, GMHI is especially favorable in terms of the simplicity of the decision rule and the biological interpretation regarding the two sets of microbes involved in classification. Importantly, our metric can be measured on a per-sample basis, requires very little parameter tuning, and foregoes the use of qualitative assessments, for example, "low" or "high" α-diversity. Furthermore, we found no significant association between library size and GMHI (mixed-effects linear regression, P = 0.45; Supplementary Fig. 5), and that, by and large, the distributions of the index for healthy individuals do not vary much between studies (Supplementary Fig. 6).
GMHI is associated with high-density lipoprotein cholesterol
To see whether GMHI can encompass certain physiological features of health, we looked for statistical associations between GMHI and well-recognized components of physiological wellness from clinical lab tests. More specifically, we searched for correlations with GMHI and the following, as reported in their original studies: circulating blood concentrations of fasting blood glucose (from 785 subjects), triglycerides (from 915 subjects), total cholesterol (from 521 subjects), low-density lipoprotein cholesterol (LDLC; from 848 subjects), and high-density lipoprotein cholesterol (HDLC; from 841 subjects). Of note, self-reported well-being, treatment regimens, and other questionnaire data were either not provided at all or too sparsely collected to have any practical or statistical significance. When selecting for moderate correlations or better, that is, |Spearman's ρ| ≥ 0.3 (P < 0.001), we identified HDLC as the only feature that was significantly associated with GMHI (ρ = 0.34, 95% confidence interval (CI): [0.28, 0.40], P = 7.19 × 10−24; Fig. 2a). In addition, we identified significantly higher abundances of HDLC in subjects with positive GMHI compared to those with negative GMHI (Mann–Whitney U test, P = 1.22 × 10−16; Fig. 2b). This moderately positive correlation is encouraging for linking GMHI to actual health, as HDLC in the bloodstream is commonly considered as "good" cholesterol, and could be protective against heart attack and stroke, according to the American Heart Association. In relevance to this point, a recent study by Kenny et al.36 showed that cholesterol metabolism by gut microbes can influence serum cholesterol concentrations, and may thereby impact cardiovascular health. Overall, our finding demonstrates the importance of integrating clinical data with gut microbiome, and also hints at the possibility of GMHI serving as an effective and reliable predictor of cardiovascular health. In contrast, fasting blood glucose (ρ = −0.06, 95% CI: [−0.12, 0.01]), triglycerides (ρ = −0.13, 95% CI: [−0.19, −0.06]), total cholesterol (ρ = 0.15, 95% CI: [0.06, 0.23]), LDLC (ρ = 0.09, 95% CI: [0.03, 0.16]), and even age (ρ = 0.04, 95% CI: [−0.01, 0.08]) were noted to have only weak or no meaningful correlations with GMHI.
Fig. 2: GMHI is associated with high-density lipoprotein cholesterol (HDLC).
a GMHI shows a moderately positive correlation with HDLC (Spearman's ρ = 0.34, 95% CI: [0.28, 0.40], P = 7.19 × 10−24), which is a key parameter of cardiovascular health, in 841 subjects. b Significantly higher abundances of HDLC were observed in subjects with positive GMHI compared to those with negative GMHI (two-sided Mann–Whitney U test, P = 1.22 × 10−16). d Cliff's Delta. The sample size of each group, whose subjects' HDLC records were available in the original studies, is shown within parentheses. Standard box-and-whisker plots (e.g., center line, median; box limits, upper and lower quartiles; whiskers, 1.5× interquartile range; circles, outliers) are used to depict groups of numerical data.
Species-level GMHI stratifies healthy and nonhealthy groups
We calculated GMHI for each stool metagenome in our meta-dataset of 4347 samples to investigate whether the distributions of GMHI differ between healthy and nonhealthy groups. We found that the gut microbiomes in healthy have significantly higher GMHIs in comparison to gut microbiomes in nonhealthy (Mann–Whitney U test, P = 5.06 × 10−212; Cliff's Delta effect size = 0.56; Fig. 3a). (Of note, Cliff's Delta (d) is a nonparametric effect-size measure that quantifies how often one value in one distribution is higher than the values in the second distribution; it is a difference between probabilities, and thus ranges from −1 to +1.) By definition of GMHI, this result reflects the dominant influence of Health-prevalent species over Health-scarce species in the healthy group, and vice versa in the nonhealthy group.
Fig. 3: Comparisons among GMHI and other ecological metrics in stratifying healthy from nonhealthy phenotypes.
a–d Significantly higher distributions of GMHI (P = 5.06 × 10−212), Shannon diversity (P = 8.50 × 10−9), and 80% abundance coverage (P = 2.30 × 10−12) were observed in gut microbiomes of healthy than in those of nonhealthy individuals, whereas higher species richness (P = 2.30 × 10−46) was observed in nonhealthy gut microbiomes. The strongest effect size (Cliff's Delta, d) was seen with GMHI. e–h The healthy group was found to have a significantly higher distribution of GMHIs than all but one (SA) of the 12 nonhealthy phenotypes. For Shannon diversity and 80% abundance coverage, only three nonhealthy phenotypes (CD, OB, and T2D) were found to have significantly different distributions compared to healthy; both properties were higher in healthy than in CD, OB, and T2D. For species richness, 7 (ACVD, CA, CC, OB, OW, RA, and T2D) of the 12 nonhealthy phenotypes were observed to have significantly higher richness than healthy; in contrast, only CD showed significantly lower richness compared to healthy. All P values shown above the violin plots were found using the two-sided Mann–Whitney U test. *P < 0.001 in two-sided Mann–Whitney U test; n.s., not significant. The sample size of each group is shown within parentheses. ACVD atherosclerotic cardiovascular disease, CA colorectal adenoma, CC colorectal cancer, CD Crohn's disease, IGT impaired glucose tolerance, OB obesity, OW overweight, RA rheumatoid arthritis, SA symptomatic atherosclerosis, T2D type 2 diabetes, UC ulcerative colitis, UW underweight. Standard box-and-whisker plots (e.g., center line, median; box limits, upper and lower quartiles; whiskers, 1.5× interquartile range; circles, outliers) are used to depict groups of numerical data.
Next, to further identify differences between healthy and nonhealthy groups, we examined multiple measures of ecological characteristics that can be defined on a per-sample basis. For α-diversity based on the Shannon index, we found significantly higher values in healthy than in nonhealthy (Mann–Whitney U test, P = 8.50 × 10−9; Cliff's Delta = 0.10; Fig. 3b). In agreement with our results, previous investigations have also reported higher diversity in the gut microbiomes of healthy controls than in those of disease patients37,38,39. In addition, we found that the minimum number of species required to comprise at least 80% of the sample's relative abundance (henceforth called "80% abundance coverage") was significantly higher in healthy compared to nonhealthy (Mann–Whitney U test, P = 2.30 × 10−12; Cliff's Delta = 0.13; Fig. 3c). This concept, as demonstrated similarly by Kraal et al.40, has been adopted in previous studies to estimate the membership of core microbiomes41,42. Finally, species richness, which is defined as the observed number of different species, was found to be significantly lower in healthy compared to nonhealthy (Mann–Whitney U test, P = 2.30 × 10−46; Cliff's Delta = −0.26; Fig. 3d).
Finally, we investigated for differences in GMHI and in these ecological characteristics between healthy and each of the 12 phenotypes of the nonhealthy group. At the individual phenotype level, the healthy group showed significantly higher GMHI levels in all but 1 (symptomatic atherosclerosis) of the 12 different disease or abnormal bodyweight conditions (Mann–Whitney U test, P < 0.001; Fig. 3e). For Shannon diversity and 80% abundance coverage, we found that only 3 (CD, obesity, and type 2 diabetes) of the 12 nonhealthy phenotypes showed statistically significant differences (Fig. 3f, g); both properties were higher in healthy for all three comparisons. For richness, we found that 8 of the 12 nonhealthy phenotypes were significantly different compared to healthy (Fig. 3h): seven of these eight were of higher richness, whereas one (CD) was of lower richness. Taken together, our results suggest that: (i) healthy and nonhealthy gut microbiomes show distinct ecological characteristics; (ii) GMHI embodies a gut microbiome signature of wellness that is generalizable against various nonhealthy phenotypes; and (iii) GMHI can distinguish healthy from nonhealthy individuals more reliably than Shannon diversity, 80% abundance coverage, and richness.
Group proportions and Shannon diversity with respect to GMHI
For increasingly higher (more positive) and lower (more negative) values of GMHI, we observed an increasing proportion of samples from healthy and nonhealthy groups, respectively (Fig. 4a). For example, 98.2% (165 of 168) of metagenome samples with GMHIs >4.0 were from the healthy group; and 81.2% (164 of 202) of metagenome samples with GMHIs <−4.0 were of nonhealthy origin. In addition, the top 10 to 100 healthy and nonhealthy stool metagenome groups (selected based on their GMHIs) clearly clustered apart from each other in PCoA ordination (Supplementary Fig. 7), in stark contrast to the case when all samples were projected simultaneously (Fig. 1c). These observations confirm that very high (or low) collective abundance of Health-prevalent species relative to that of Health-scarce species is strongly connected to being healthy (or nonhealthy).
Fig. 4: Changes in group proportions and in Shannon diversity with respect to GMHI.
a All 4347 metagenomes were binned according to their GMHI values (x-axis). Each gray bar indicates the total number of samples in each bin (y-axis, right). Points indicate proportions (i.e., percentages) of samples in each bin corresponding to either healthy or nonhealthy individuals (y-axis, left). In bins with a positive range of GMHIs, the majority of samples classified as healthy; in contrast, samples in bins with a negative range of GMHIs mostly classified as nonhealthy. This trend was more pronounced towards bins on the far right and left. b GMHI stratifies healthy (n = 2636) and nonhealthy (n = 1711) groups more strongly compared to Shannon diversity. Each point in the scatter-plot corresponds to a metagenome sample (4347 in total). Histograms show the distribution of healthy (blue) and nonhealthy (orange) samples based on the parameter of each axis. In general, GMHI and Shannon diversity demonstrate a weak correlation (Spearman's ρ = 0.17, 95% CI: [0.14, 0.19], P = 1.7 × 10−28). The P value (H0: ρ = 0) was determined by using a t-distribution with n − 2 degrees of freedom, where n is the total number of observations.
GMHI and Shannon diversity were compared for each sample to examine their overall concordance. As shown in Fig. 4b, GMHI clearly performed much better in stratifying the healthy and nonhealthy groups compared to Shannon diversity. A small yet significant relationship was found between our metric and this conventional measure of gut health (Spearman's ρ = 0.17, 95% CI: [0.14, 0.19], P = 1.66 × 10−28). In addition, similar results were seen when GMHI was compared with 80% abundance coverage (Spearman's ρ = 0.22, 95% CI: [0.19, 0.25], P = 8.48 × 10−48) and with richness (Spearman's ρ = −0.27, 95% CI: [−0.30, −0.24], P = 4.27 × 10−74) (Supplementary Fig. 8).
Intra-study analyses favor GMHI over other ecology metrics
We next examined how well GMHI and other features of microbial ecology (i.e., Shannon diversity, 80% abundance coverage, and species richness) could distinguish healthy and nonhealthy phenotypes within individual studies. Specifically, in each of the 12 studies (out of 34 total) wherein at least 10 stool metagenome samples from both case (i.e., disease or abnormal bodyweight conditions) and control (i.e., healthy) subjects were available, we compared GMHI, Shannon diversity, 80% abundance coverage, and species richness between healthy and nonhealthy phenotype(s). By focusing on datasets from individual studies one by one, this approach not only removes a major source of batch effects, but also provides a good means to investigate the robustness of our previously observed trends (when healthy and nonhealthy samples were compared against each other in aggregate groups) across multiple, smaller studies.
We found that GMHI in healthy was significantly higher than that in any nonhealthy phenotype for 11 out of 28 case–control comparisons (Fig. 5a). For Shannon diversity and 80% abundance coverage, we found significantly higher values in healthy than in nonhealthy phenotypes for two and four case–control comparisons, respectively (Fig. 5b, c). Last, we found species richness in healthy to be significantly lower than that in nonhealthy phenotypes for three case–control comparisons (Fig. 5d). Clearly, the performance of GMHI was not perfect (and likewise for other ecological characteristics), as the expected trend from the prior pooled analysis was not replicable for all case–control comparisons within every study; overall though, GMHI strongly outperformed other microbiome ecological characteristics in distinguishing case and control.
Fig. 5: GMHI generally outperforms other microbiome ecological characteristics in distinguishing case and control across multiple study-specific comparisons.
In each of the 12 studies wherein at least 10 case (i.e., disease or abnormal bodyweight conditions) and at least 10 control (i.e., healthy) subjects were available, stool metagenomes were analyzed to compare a GMHI, b Shannon diversity, c 80% abundance coverage, and d species richness between healthy and nonhealthy phenotype(s). GMHI was found to have a significantly higher distribution in healthy for 11 (out of 28) case–control comparisons across nine different studies; Shannon diversity and 80% abundance coverage were found to have significantly higher distributions in healthy for two and four case–control comparisons (across two and four studies), respectively; and species richness was found to have a significantly lower distributions in healthy for three case–control comparisons across three different studies. Each study's phenotype sample size is shown within parentheses to the right of the phenotype abbreviation. Standard box-and-whisker plots (e.g., center line, median; box limits, upper and lower quartiles; whiskers, 1.5× interquartile range; points, samples) are used to depict groups of numerical data. The same colors in boxplots were used for the same phenotypes. P values (two-sided Mann–Whitney U test) for each study-specific comparison between healthy and nonhealthy phenotypes are shown adjacent to the boxplots accordingly: * and Ψ indicates significantly different distributions consistent with, and opposite to, respectively, the previously observed results when healthy and nonhealthy groups were compared in aggregate. * or Ψ0.01 ≤ P value < 0.05; ** or ΨΨ0.001 ≤ P value < 0.01; *** or ΨΨΨ0.0001 ≤ P value < 0.001; **** or ΨΨΨΨP value < 0.0001. ACVD atherosclerotic cardiovascular disease, CA colorectal adenoma, CC colorectal cancer, CD Crohn's disease, IGT impaired glucose tolerance, OB obesity, OW overweight, RA rheumatoid arthritis, SA symptomatic atherosclerosis, T2D type 2 diabetes, UC ulcerative colitis, UW underweight.
Analogous to the analysis above (wherein healthy was compared to each separate nonhealthy phenotype within individual studies), we compared healthy against a general nonhealthy phenotype, in which all disease samples were lumped together, when applicable. Importantly, comparisons were still made within individual studies. We found that there were statistically significant differences in GMHI between cases and controls (Mann–Whitney U test, P < 0.05; Supplementary Table 4) in 6 of the 12 studies. In contrast, we found statistically significant differences in Shannon diversity, 80% abundance coverage, and richness between cases and controls in two, three, and three (of 12) studies, respectively.
Validation of GMHI reproducibility on independent cohorts
Evaluation of any biomarker or molecular signature on independent patient samples is the gold standard for assessing its robustness15. To confirm the reproducibility of our prediction results in stratifying healthy and nonhealthy phenotypes (Fig. 3), we leveraged GMHI to predict the health status of 679 individuals whose stool metagenome samples were not part of the original formulation of GMHI. For this, we used gut microbiome data from an additional eight published studies (Supplementary Table 5), which include stool metagenomes from healthy subjects and patients with ankylosing spondylitis (AS), colorectal adenoma, colorectal cancer, Crohn's disease (CD), liver cirrhosis (LC), and nonalcoholic fatty liver disease (NAFLD). In addition, we utilized our extensive biobank of stool collections to gather our own set of samples from patients with rheumatoid arthritis (RA) ("Methods"; see Supplementary Data 4 for subject meta-data relating to both clinical and nonclinical factors). All metagenome samples in this validation dataset were pooled into one of two groups (i.e., healthy or nonhealthy), as demonstrated above.
In agreement with our results on the discovery cohort (training data), GMHIs from stool metagenomes of the healthy validation group (n = 118) were significantly higher than those of the nonhealthy validation group (n = 561) (Mann–Whitney U test, P = 3.49 × 10−28; Cliff's Delta = 0.64; Fig. 6a). In addition, the balanced accuracy resulted in 73.7%, as the classification accuracy for the healthy and nonhealthy validation group was 77.1% (91 of 118) and 70.2% (394 of 561), respectively. Notably, these results were better than the performances on the discovery cohort, wherein balanced accuracy was 69.7%, and accuracy on the healthy and nonhealthy group was 75.6% (1993 of 2636) and 63.8% (1092 of 1711), respectively.
Fig. 6: GMHI demonstrates strong reproducibility on an independent validation cohort.
The validation cohort (679 stool metagenome samples) consisted of 12 total sub-cohorts ranging across eight healthy and nonhealthy phenotypes from nine different studies. a GMHIs from stool metagenomes of the healthy group were significantly higher than those of the nonhealthy group (two-sided Mann–Whitney U test, P = 3.49 × 10−28). d Cliff's Delta. b All three healthy sub-cohorts (H1, H2, and H3) were found to have significantly higher distributions of GMHI than seven (of nine) nonhealthy sub-cohorts (AS4, CC5-I, CC5-J, CD6, LC7, NAFLD8, and RA9). No significant differences were found among H1, H2, and H3. The number in superscript adjacent to phenotype abbreviations corresponds to a particular study used in validation (see Supplementary Table 5 for study information). Standard box-and-whisker plots (e.g., center line, median; box limits, upper and lower quartiles; whiskers, 1.5× interquartile range; points, samples) are used to depict groups of numerical data. *Significantly higher distribution in healthy sub-cohort (two-sided Mann–Whitney U test, P < 0.01). The number adjacent to * indicates the healthy sub-cohort (H1, H2, or H3) to which the respective sub-cohort was compared. The sample size of each group or cohort is shown within parentheses. AS ankylosing spondylitis, CA colorectal adenoma, CC colorectal cancer, CD Crohn's disease, H healthy, LC liver cirrhosis, NAFLD, nonalcoholic fatty liver disease, RA rheumatoid arthritis.
Of note, we also compared the classification accuracy of GMHI to those of classifiers based upon the Health-prevalent species and Shannon diversity (see Supplementary Methods), and to that of a more intricate classification algorithm (Random Forests). In regards to balanced accuracies on the training data, the classifiers based upon Health-prevalent species (χ = 66.3%) and Shannon diversity (χ = 53.6%) performed comparable to, or much worse than, GMHI (χ = 69.7%); furthermore, balanced accuracy on the independent validation dataset for Health-prevalent species and Shannon diversity resulted in 59.3 and 47.0%, respectively (Supplementary Tables 6 and 7, respectively). On the other hand, the Random Forests classifier ("Methods") achieved a remarkable accuracy on the training data (χ = 98.5%). However, building complex decision rules entails the risk of over-fitting. Surely enough, this nearly perfect accuracy was mostly in part a result of outstanding over-fitting, evidenced by the poor accuracy of 52.3% (balanced accuracy) on the 679 samples of the validation cohort. In Supplementary Table 8, we provide a summary of all accuracies for classifying healthy vs. nonhealthy by the various classifiers reported in this study.
To investigate GMHI performances on the validation cohort more closely, we examined the 12 total sub-cohorts (defined per unique phenotype per individual study) ranging across eight healthy and nonhealthy phenotypes from eight additional published studies and one newly sequenced batch. As shown in Fig. 6b, all three healthy sub-cohorts were found to have significantly higher distributions of GMHI than seven (of nine) nonhealthy phenotype sub-cohorts (Mann–Whitney U test, P < 0.01; see Supplementary Table 9 for Cliff's Deltas). The classification accuracies for these three healthy sub-cohorts were 87.5% (28 of 32), 74.1% (43 of 58), and 71.4% (20 of 28); alternatively, the classification accuracies for the nonhealthy phenotype sub-cohorts were the following: 94.5% (155 of 164) for LC; 75.6% (65 of 86) for NAFLD; 73.3% (11 of 15) for CD; 67.3% (33 of 49) for RA; 55.7% (54 of 97) for AS; 37.0% (10 of 27) for CA; and 77.5% (31 of 40), 47.5% (29 of 61), and 27.3% (6 of 22) for three different cohorts of CC. Strikingly, GMHI performed well (>75.0%) in predicting adverse health for LC and NAFLD, although stool metagenomes from patients with liver disease were not part of the original discovery cohort. This finding suggests that GMHI could be applied beyond the original 12 phenotypes (of the nonhealthy group) used during the index training process. Overall, the strong reproducibility of GMHI implies that the highly diverse and complex features of gut microbiome dysbiosis implicated in pathogenesis were reasonably well captured during the dataset integration and original formulation of GMHI. Thereby, our results support previous findings by Duvallet et al.26 in regards to the presence of a generalized disease-associated gut microbial signature, which was observed to be shared across multiple studies and pathologies. Finally, from similar analyses for Shannon diversity, 80% abundance coverage, and species richness on the validation cohort, we were able to conclude that GMHI is the most accurate, robust, and clinically meaningful classifier compared to these other ecological characteristics (Supplementary Note 2 and Supplementary Fig. 9).
In this study, we present the GMHI, a simple and biologically interpretable metric to quantify the likelihood of disease presence from a gut microbiome sample. At first, we envisioned that the most intuitive way to determine how closely one's microbiome resembles that of a healthy (or nonhealthy) population is to compare the collective abundances of Health-prevalent and of Health-scarce species. By pooling massive amounts of publicly available data (4347 publicly available, shotgun metagenomic data of gut microbiomes from 34 published studies), we identified a small consortium of 50 microbial species associated with human health to serve as features for our classification model: 7 and 43 species were prevalent and scarce, respectively, in the healthy group compared to the nonhealthy group. In regards to classification accuracy, GMHI distinguished healthy from nonhealthy (as well as from individual diseases) far better than methods adopted from ecological principles (e.g., α-diversity indices), thereby paving a path forward to evaluate human (gut microbiome) health through stool metagenomic profiling. Notably, this framework can be applied to other body niches, for example, quantifying health in skin or oral microbiomes. When demonstrating the potential of GMHI on independent validation datasets, we obtained strong prediction results for healthy individuals, and for cohorts with autoimmunity and liver disease. The strong reproducibility on validation datasets suggests that sufficient dataset integration across a large population could lead to robust predictors of health. This may be due, in part, to the signature encompassing more of the heterogeneity across various sources and conditions, while amplifying signal (against noise) from the repeated phenotype characteristics.
Despite the strong reproducibility in classification accuracy demonstrated on the validation datasets, more is left to be desired in regards to achieving higher accuracies. We conjecture that misclassifications were partly because of the complex43,44,45, stochastic46,47, and highly personalized18,48 nature of gut microbiome ecologies; all of which complicate the identification of reliable signatures of health. In addition, sample collection and processing procedures, laboratory personnel, study run-dates, measurement instruments, and so forth are tremendously challenging to control in population-scale investigations. In the long term, in order to find even more accurate gut microbiome-based markers, we envision integrating larger data repositories to take into consideration a higher number of samples and sources of heterogeneity.
Several limitations of our study should be noted when interpreting our results. First, as the stool metagenomes were collected from over 40 published studies, we cannot entirely exclude experimental and technical inter-study batch effects. Our efforts to curtail batch effects include: (i) consensus preprocessing, that is, downloading all raw shotgun metagenomes and reprocessing each sample uniformly using identical bioinformatics methods; (ii) using frequencies of a signal (i.e., prevalence of "present" microbes) to identify significant associations, rather than comparing or averaging effect sizes between populations, or performing data transformations that may lead to spurious conclusions49; and (iii) validating the reproducibility of GMHI on independent datasets. Second, given our selection criteria ("Methods"), our study does not include all publicly available gut microbiome studies and samples. Certainly, more studies and samples can be taken into consideration under more relaxed criteria. Third, in an effort to be as precise as possible in describing taxonomic features of the human gut, our metagenomic analyses were performed using species-level abundances; however, microbial strains are clearly the most clinically informative and actionable unit50,51. Moreover, different strains within the same species can have significantly different associations with disease52,53,54,55, which could not be considered in our study. Nevertheless, our shotgun metagenomic approach is a significant advancement over 16s rRNA gene amplicon sequencing, which are known to be mostly limited to genus-level investigations56,57. Fourth, in the nonhealthy group, we pooled samples from only 12 phenotypes. Certainly, many more pathological states have been linked to the gut microbiome, including neurodegenerative and psychiatric disorders58,59,60. Thus, future studies will need to continuously update and expand our findings by encompassing a much broader range of conditions as new data become available. Fifth, we did not consider metagenomic functional profiles to define gut ecosystem health as demonstrated extensively by others27,61,62,63, as this too was outside the scope of our study. For microbiomes of any phenotype of interest, we posit that analyzing both taxonomic composition and functional potential are both important and complementary directions. Last, while we definitely tried to be as inclusive as possible of various geographies, ethnicities/races, and cultures, we do acknowledge that complete elimination of biases is practically impossible. Certainly, for future works, we plan to iteratively expand our application to encompass broader ranges of subjects, including those from underdeveloped countries and minority ethnicities/races, to better understand microbiome diversity and foster inclusion in microbiome research64.
Multi-study integration of human stool metagenomes
We performed exhaustive keyword searches (e.g., "gut microbiome," "metagenome," "whole-genome shotgun (WGS)") in PubMed and Google Scholar for published studies with publicly available WGS metagenome data of human stool (gut microbiome) and corresponding subject meta-data (as of March 2018). In studies wherein multiple samples were taken per individual across different time-points, we included only the first or baseline sample in the original study. We excluded studies pertaining to diet or medication interventions, or those with fewer than 10 samples. Samples from subjects who were <10 years of age were also excluded from our analysis. Last, samples that were collected from disease controls, but were not reported as healthy nor had any mentioning of diagnosed disease in the original study, were excluded from our analysis. Raw sequence files (.fastq) were downloaded from the NCBI Sequence Read Archive and European Nucleotide Archive databases (Supplementary Data 1) for the study analysis.
Reclassification of healthy samples based on reported BMI
Healthy individuals, regardless of whether they had been determined as healthy in the original studies, were considered to be part of the nonhealthy group if their reported BMI fell within the range of underweight (BMI < 18.5), overweight (BMI ≥ 25 and <30), or obese (BMI ≥ 30). Stool metagenome samples from such individuals were reclassified as underweight, overweight, or obese in our analysis.
Quality control of sequenced reads
Sequence reads were processed with the KneadData v0.5.1 quality-control pipeline (http://huttenhower.sph.harvard.edu/kneaddata), which uses Trimmomatic v0.36 and Bowtie2 v0.1 for removal of low-quality read bases and human reads, respectively. Trimmomatic v0.36 was run with parameters SLIDINGWINDOW:4:30, and Phred quality scores were thresholded at "<30." Illumina adapter sequences were removed, and trimmed nonhuman reads shorter than 60 bp in nucleotide length were discarded. Potential human contamination was filtered by removing reads that aligned to the human genome (reference genome hg19). Furthermore, stool metagenome samples of low read count after quality filtration (<1M reads) were excluded from our analysis.
Species-level taxonomic profiling
Taxonomic profiling was done using the MetaPhlAn2 v2.7.0 phylogenetic clade identification pipeline32 using default parameters. Briefly, MetaPhlAn2 classifies metagenomic reads to taxonomies based on a database (mpa_v20_m200) of clade-specific marker genes derived from ~17,000 microbial genomes (corresponding to ~13,500 bacterial and archaeal, ~3500 viral, and ~110 eukaryotic species).
Sample filtering based on taxonomic profiles
After taxonomic profiling, the following stool metagenome samples were discarded from our analysis: (i) samples composed of >5% unclassified taxonomies (100 samples); and (ii) phenotypic outliers according to a dissimilarity measure. More specifically, Bray–Curtis distances were calculated between each sample of a particular phenotype and a hypothetical sample in which the species' abundances were taken from the medians across those samples. A sample was considered as an outlier, and thereby removed from further analysis, when its dissimilarity exceeded the upper and inner fence (i.e., >1.5 times outside the interquartile range above the upper quartile and below the lower quartile) among all dissimilarities. This process removed 67 metagenome samples.
Species removal based on taxonomic profiles
As taxonomic assignment based on clade-specific marker genes may be problematic for viruses65,66, we excluded the 298 of viral origin from our analysis. Species that were labeled as either unclassified or unknown (118 species), or those of low prevalence (i.e., observed in <1% of the samples included in our meta-dataset; 472 species), were also excluded. Eventually, 313 microbial species across 4347 stool metagenome samples remained in our study for further analysis (Supplementary Data 2).
PCoA based on taxonomic profiles
The R packages "ade4" v1.7-15 and "vegan" v2.5.6 were used to perform PCoA ordination with Bray–Curtis dissimilarity as the distance measure on the stool metagenome samples, which were comprised of arcsine square root-transformed relative abundances of the aforementioned 313 microbial species identified by MetaPhlAn2. 999 permutations ("adonis2" function in the R "vegan" package v2.5.6) were performed, while random permutations were constrained within studies by using the "strata" option.
Calculation of microbiome ecological characteristics
The R package "vegan" v2.5.6 was used to calculate Shannon diversity (Shannon index) and species richness based on the species abundance profiles for each sample of our meta-dataset. To identify the 80% abundance coverage for a stool metagenome sample, the smallest number of microbial species that comprise at least 80% of the total relative abundance was identified.
Identifying microbial species more frequently observed in healthy than in nonhealthy (and vice versa):
Let pH,m and pN,m be the prevalence of microbial species m, that is, proportion of samples in a given group where m is "present" (or relative abundance ≥1.0 × 10−5), in the healthy group H and nonhealthy group N, respectively. Remark: The relative abundances for all detectable species in a microbiome (metagenome) sample sums to 1.
For m, the prevalence fold change \(f_m^{H,N}\) and prevalence difference \(d_m^{H,N}\), defined as \(\frac{{p_{H,m}}}{{p_{N,m}}}\) and pH,m − pN,m, respectively, is identified.
Let θf and θd be defined as the minimum thresholds for \(f_m^{H,N}\) and \(d_m^{H,N}\), respectively. For all detectable species in a microbiome sample, those that satisfy \(f_m^{H,N} \ge \theta _f\) and \(d_m^{H,N} \ge \theta _d\) are identified. These species are included as an element of "Health-prevalent" species MH, or the set of species more frequently observed in group H than in group N.
To identify "Health-scarce" species MN, or the set of species more frequently observed in group N than in group H, steps (b) through (c) are repeated with the following considerations:
For m, let \(f_m^{N,H}\) and \(d_m^{N,H}\) be defined as \(\frac{{p_{N,m}}}{{p_{H,m}}}\) and pN,m − pH,m, respectively.
The same thresholds θf and θd are used to identify MN. In this regard, the species that are eventually chosen to compose MH and MN are both dependent on θf and θd.
Finally, all detectable species that satisfy \(f_m^{N,H} \ge \theta _f\) and \(d_m^{N,H} \ge \theta _d\) are included in MN.
Identifying \(\psi _{M_H}\)(or \(\psi _{M_N}\)), that is, the "collective abundance" of species in MH (or MN) in a microbiome sample:
\(\psi _{M_H}\) is defined as the "collective abundance" of MH species in a microbiome sample. The calculation of \(\psi _{M_H}\) takes into consideration the following:
Species richness, that is, the numeric count of "present" species of MH.
(geometric) Mean of their relative abundances.
Basic assumptions:
\(\psi _{M_H}\) is positively correlated with \(R_{M_H}\), or the richness of MH species. Thus, (their correlation) \(\rho ( {\psi _{M_H},R_{M_H}} ) > 0\). Remark: Due to the possible large discrepancy between the cardinality (set size) of MH and that of MN, the proportion of "present" MH species is used. As such, \(R_{M_H}\) is replaced with \(\frac{{R_{M_H}}}{{\left| {M_H} \right|}}\). Thus, \(\rho ( {\psi _{M_H},\frac{{R_{M_H}}}{{| {M_H} |}}} ) > 0\).
\(\psi _{M_H}\) is positively correlated with 〈MH〉, or the mean abundance of species in MH. Thus, \(\rho ( {\psi _{M_H},\langle {M_H}\rangle }) > 0\). Remark: As it is common in microbiome data to have discrepancies between species' relative abundances to span several orders of magnitude, the geometric mean, rather than the arithmetic mean, is more appropriate to represent the mean relative abundance of MH species. More specifically, the Shannon's diversity index, which is a weighted geometric mean (by definition) and commonly applied in ecological contexts, is used. Thus, for simplicity, \(\left\langle {M_H} \right\rangle \approx \mathop {\sum}\nolimits_{j \in I_{M_H}} {| {n_{j} ln(n_{j})}|}\) is assumed, where \(I_{M_H}\) is the index set of MH, and nj is the relative abundance of species j in \(I_{M_H}\).
Given the assumptions in (b), as well as the nonnegativity of \(\frac{{R_{M_H}}}{{\left| {M_H} \right|}}\) and \(\mathop {\sum}\nolimits_{j \in I_{M_H}} {| {n_j ln ( {n_j} )} |}\), \(\psi _{M_H}\) is simply formulated as a product of the aforementioned two traits. Thus, let \(\psi _{M_H} = \frac{{R_{M_H}}}{{\left| {M_H} \right|}}\mathop {\sum}\nolimits_{j \in I_{M_H}} {| {n_{j} ln (n_{j})}|}\).
Analogously, let \(\psi _{M_N} = \frac{{R_{M_N}}}{{| {M_N}|}}\mathop {\sum}\nolimits_{j \in I_{M_N}} {| {n_jln( {n_j})}|}\)
Identifying \(h_{i,M_H,M_N}\), that is, ratio of \(\psi _{M_H}\)to \(\psi _{M_N}\), in sample i:
Formally, the log ratio of \(\psi _{M_H}\) to \(\psi _{M_N}\) in sample i can be written as
$$h_{i,M_H,M_N} = {\mathrm{log}}_{10}\left( {\frac{{\frac{{R_{M_H}}}{{\left| {M_H} \right|}}\mathop {\sum}\nolimits_{j \in I_{M_H}} {\left| {n_jln\left( {n_j} \right)} \right|} }}{{\frac{{R_{M_N}}}{{\left| {M_N} \right|}}\mathop {\sum}\nolimits_{j \in I_{M_N}} {\left| {n_jln\left( {n_j} \right)} \right|} }}} \right).$$
By definition, |MH| and |MN| is the highest richness that can be obtained by MH and MN species, respectively, in a particular microbiome sample. However, the possibility that these maximum values are rarely obtained cannot be ruled out; if so, then consequently, having a larger set size of MH (or MN) can generally result in a lower distribution of \(\frac{{R_{M_H}}}{{\left| {M_H} \right|}}\) (or \(\frac{{R_{M_N}}}{{\left| {M_N} \right|}}\)), potentially leading to biases in \(h_{i,M_H,M_N}\) when \(\left| {M_H} \right| \gg \left| {M_N} \right|\) or \(\left| {M_H} \right| \ll \left| {M_N} \right|\). Therefore, the upper limits that can be eventually used in replacement of |MH| and |MN| in Eq. (4) should reflect more of what is actually observed in real microbiome data, for example, samples ranked according to the magnitude observed between \(R_{M_H}\) and \(R_{M_N}\). In this regard, the following procedure to find alternative measures for |MH| and |MN| is used:
Identify \(R_{M_H}\) and \(R_{M_N}\) for all microbiome samples in groups H and N.
Rank order all samples consecutively by two criteria: first, by all values of \(R_{M_N}\) in ascending order (from lowest to highest); and then, by all values of \(R_{M_H}\) in descending order (from highest to lowest). This sorting strategy prioritizes having the highest possible \(R_{M_H}\) (but with the constraint of having \(R_{M_N} \approx 0\)) for the most top-ranked samples; and having the highest possible \(R_{M_N}\) (but with the constraint of having \(R_{M_H} \approx 0\)) for the most bottom-ranked samples.
Let kH be the closest integer to 1% of the number of samples in group H. As H is composed of 2636 samples, let kH be 26. Analogously, as N is composed of 1711 samples, let kN be 17.
Denote |MH|′ as the median \(R_{M_H}\) from the top kH samples, and denote |MN|′ as the median \(R_{M_N}\) from the bottom kN samples.
Replace |MH| and |MN| in Eq. (4) with |MH|′ and |MN|′, respectively.
In summary, the ratio of \(\psi _{M_H}\)to \(\psi _{M_N}\) in gut microbiome sample i can be written as
$$h_{i,M_H,M_N} = {\mathrm{log}}_{10}\left( {\frac{{\frac{{R_{M_H}}}{{\left| {M_H} \right|\prime }}\mathop {\sum}\nolimits_{j \in I_{M_H}} {\left| {n_jln\left( {n_j} \right)} \right|} }}{{\frac{{R_{M_N}}}{{\left| {M_N} \right|\prime }}\mathop {\sum}\nolimits_{j \in I_{M_N}} {\left| {n_jln\left( {n_j} \right)} \right|} }}} \right).$$
Calculating the balanced accuracy of \(h_{M_H,M_N}\):
The relative abundances of species in MH and those in MN for microbiome sample i can be provided as input features for: (i) \(\psi _{M_H}\) and \(\psi _{M_N}\), respectively; and (ii) \(h_{i,M_H,M_N}\), which in turn can classify sample i as healthy (i.e., \(h_{i,M_H,M_N} > 0\)), nonhealthy (i.e., \(h_{i,M_H,M_N} < 0\)), or neither (i.e., \(h_{i,M_H,M_N} = 0\)).
The classification accuracy or predictive performance of \(h_{i,M_H,M_N}\) is found by testing it on all samples in groups H and N, and then by finding the balanced accuracy \(\chi _{M_H,M_N}\) defined in Eq. (3).
Determining optimal sets \(M_H^\gamma\)and \(M_N^\gamma\):
The final, optimal sets of \(M_H^\gamma\) and \(M_N^\gamma\) are found by first considering a range of thresholds θf and θd. Every pair of θf and θd gives different sets of MH and MN, and in turn, different values of balanced accuracy \(\chi _{M_H,M_N}\) (see Supplementary Table 1).
The final, optimal sets of \(M_H^\gamma\) and \(M_N^\gamma\) (and their corresponding \(\theta _f^\gamma\) and \(\theta _d^\gamma\)) are determined as those that result in the highest balanced accuracy \(\chi _{M_H,M_N}^{\mathrm{max}}\).
MetaCyc pathway functional profiling of stool metagenomes
MetaCyc pathway-level relative abundances in each stool metagenome were quantified by the HUManN v2.0 pipeline63 using default parameters. The EC-filtered UniRef90 gene family database was integrated within the pipeline. Pathways that were unmapped (or unintegrated) were excluded from the analyses.
Designing a classifier based upon Random Forests
A classifier based upon a Random Forests algorithm was designed and curated in Python v3.6.4., while model implementation was performed in the "scikit-learn" Python package v0.23.1.
Stool sample collection and processing
All stool samples from patients with RA were obtained following written informed consent. The collection of biospecimens was approved by the Mayo Clinic Institutional Review Board (#14-000616). Stool samples from patients with RA were stored in their house-hold freezer (−20 °C) prior to shipment on dry ice to the Medical Genome Facility Research Core at Mayo Clinic (Rochester, MN). Once received, the samples were stored at −80 °C until DNA extraction. DNA extraction from stool samples was conducted as follows: aliquots were created from parent stool samples using a tissue punch, and the resulting child samples were then mixed with reagents from the Qiagen Power Fecal Kit. This included adding 60 μL of reagent C1 and the contents of a power bead tube (garnet beads and power bead solution). These were then vigorously vortexed to bring the sample punch into solution and centrifuged at 18,000 × G for 15 min. From there, the samples were added into a mixture of magnetic beads using a JANUS liquid handler. The samples were then run through a Chemagic MSM1 according to the manufacturer's protocol. After DNA extraction, paired-end libraries were prepared using 500 ng genomic DNA according to the manufacturer's instructions for the NEB Next Ultra Library Prep Kit (New England BioLabs). The concentration and size distribution of the completed libraries was determined using an Agilent Bioanalyzer DNA 1000 chip (Santa Clara, CA) and Qubit fluorometry (Invitrogen, Carlsbad, CA). Libraries were sequenced at 23–70 million reads per sample following Illumina's standard protocol using the Illumina cBot and HiSeq 3000/4000 PE Cluster Kit. The flow cells were sequenced as 150 × 2 paired-end reads on an Illumina HiSeq 4000 using the HiSeq 3000/4000 Sequencing Kit and HiSeq Control Software HD 3.4.0.38. Base-calling was performed using Illumina's RTA version 2.7.7.
Raw sequencing data accession IDs of all publicly available stool metagenome samples (and their corresponding studies) used in all analyses of this study are available in Supplementary Data 1 and Supplementary Data 4. Sequences for the dataset containing rheumatoid arthritis stool metagenomes used for GMHI validation have been deposited at NCBI's Sequence Read Archive (SRA) data repository (BioProject number PRJNA598446), and can be downloaded without any restrictions. The deposited sequences include .fastq files for 49 patients with rheumatoid arthritis. Measurements were taken from distinct samples. Human reads were identified and removed prior to data upload.
Code availability
R scripts demonstrating how to reproduce all findings shown in the main figures, as well as how to calculate GMHI for a given stool metagenome sample, are available at https://github.com/jaeyunsung/GMHI_2020.
Wirbel, J. et al. Meta-analysis of fecal metagenomes reveals global microbial signatures that are specific for colorectal cancer. Nat. Med. 25, 679–689 (2019).
Thomas, A. M. et al. Metagenomic analysis of colorectal cancer datasets identifies cross-cohort microbial diagnostic signatures and a link with choline degradation. Nat. Med. 25, 667–678 (2019).
Scher, J. U. & Abramson, S. B. The microbiome and rheumatoid arthritis. Nat. Rev. Rheumatol. 7, 569–578 (2011).
Jangi, S. et al. Alterations of the human gut microbiome in multiple sclerosis. Nat. Commun. 7, 12015 (2016).
CAS PubMed PubMed Central ADS Google Scholar
Lloyd-Price, J. et al. Multi-omics of the gut microbial ecosystem in inflammatory bowel diseases. Nature 569, 655–662 (2019).
Everard, A. & Cani, P. D. Diabetes, obesity and gut microbiota. Best. Pract. Res. Clin. Gastroenterol. 27, 73–83 (2013).
Mazidi, M., Rezaie, P., Kengne, A. P., Mobarhan, M. G. & Ferns, G. A. Gut microbiome and metabolic syndrome. Diabetes Metab. Syndr. 10, S150–S157 (2016).
Dao, M. C. et al. Akkermansia muciniphila and improved metabolic health during a dietary intervention in obesity: relationship with gut microbiome richness and ecology. Gut 65, 426–436 (2016).
Allaband, C. et al. Microbiome 101: studying, analyzing, and interpreting gut microbiome data for clinicians. Clin. Gastroenterol. Hepatol. 17, 218 (2019).
Staley, C., Kaiser, T. & Khoruts, A. Clinician guide to microbiome testing. Dig. Dis. Sci. 63, 3167–3177 (2018).
McBurney, M. I. et al. Establishing what constitutes a healthy human gut microbiome: state of the science, regulatory considerations, and future directions. J. Nutr. 149, 1882–1895 (2019).
Hagerty, S. L., Hutchison, K. E., Lowry, C. A. & Bryan, A. D. An empirically derived method for measuring human gut microbiome alpha diversity: demonstrated utility in predicting health-related outcomes among a human clinical sample. PLoS ONE 15, e0229204 (2020).
Hood, L. & Flores, M. A personal view on systems medicine and the emergence of proactive P4 medicine: predictive, preventive, personalized and participatory. N. Biotechnol. 29, 613–624 (2012).
Fröhlich, H. et al. From hype to reality: data science enabling personalized medicine. BMC Med. 16, 1–15 (2018).
Sung, J., Wang, Y., Chandrasekaran, S., Witten, D. M. & Price, N. D. Molecular signatures from omics data: from chaos to consensus. Biotechnol. J. 7, 946–957 (2012).
Kashyap, P. C., Chia, N., Nelson, H., Segal, E. & Elinav, E. Microbiome at the frontier of personalized medicine. Mayo Clin. Proc. 92, P15–P29 (2017).
Price, N. D. et al. A wellness study of 108 individuals using personal, dense, dynamic data clouds. Nat. Biotechnol. 35, 747–756 (2017).
Goodrich, J. K. et al. Human genetics shape the gut microbiome. Cell 159, 789–799 (2014).
Falony, G. et al. Population-level analysis of gut microbiome variation. Science 352, 560–564 (2016).
CAS PubMed ADS Google Scholar
Zhernakova, A. et al. Population-based metagenomics analysis reveals markers for gut microbiome composition and diversity. Science 352, 565–569 (2016).
McDonald, D. et al. American gut: an open platform for citizen science microbiome research. mSystems 3, e00031-18 (2018).
Goodrich, J. K. et al. Genetic determinants of the gut microbiome in UK twins. Cell Host Microbe 19, 731–743 (2016).
Amann, R. I. et al. Toward unrestricted use of public genomic data. Science 363, 350–352 (2019).
Finucane, M. M., Sharpton, T. J., Laurent, T. J. & Pollard, K. S. A taxonomic signature of obesity in the microbiome? Getting to the guts of the matter. PLoS ONE 9, e84689 (2014).
PubMed PubMed Central ADS Google Scholar
Sze, M. A. & Schloss, P. D. Looking for a signal in the noise: revisiting obesity and the microbiome. MBio 7, e01018-16 (2016).
Duvallet, C., Gibbons, S. M., Gurry, T., Irizarry, R. A. & Alm, E. J. Meta-analysis of gut microbiome studies identifies disease-specific and shared responses. Nat. Commun. 8, 1784 (2017).
Armour, C. R., Nayfach, S., Pollard, K. S. & Sharpton, T. J. A metagenomic meta-analysis reveals functional signatures of health and disease in the human gut microbiome. mSystems 4, e00332-18 (2019).
Sung, J. et al. Multi-study integration of brain cancer transcriptomes reveals organ-level molecular signatures. PLoS Comput. Biol. 9, e1003148 (2013).
Cohn, L. D. & Becker, B. J. How meta-analysis increases statistical power. Psychol. Methods 8, 243–253 (2003).
Jackson, D. & Turner, R. Power analysis for random-effects meta-analysis. Res. Synth. Methods 8, 290–302 (2017).
Ma, S. et al. Measuring the effect of inter-study variability on estimating prediction error. PLoS ONE 9, e110840 (2014).
Truong, D. T. et al. MetaPhlAn2 for enhanced metagenomic taxonomic profiling. Nat. Methods 12, 902–903 (2015).
Coyte, K. Z., Schluter, J. & Foster, K. R. The ecology of the microbiome: networks, competition, and stability. Science 350, 663–666 (2015).
Lloyd-Price, J., Abu-Ali, G. & Huttenhower, C. The healthy human microbiome. Genome Med. 8, 51 (2016).
Inkpen, S. A. Health, ecology and the microbiome. Elife 8, e47626 (2019).
Kenny, D. J. et al. Cholesterol metabolism by uncultured human gut bacteria influences host cholesterol level. Cell Host Microbe https://doi.org/10.1016/j.chom.2020.05.013 (2020).
Qin, J. et al. A metagenome-wide association study of gut microbiota in type 2 diabetes. Nature 490, 55–60 (2012).
Vogt, N. M. et al. Gut microbiome alterations in Alzheimer's disease. Sci. Rep. 7, 13537 (2017).
Zhou, Y. et al. Gut microbiota offers universal biomarkers across ethnicity in inflammatory bowel disease diagnosis and infliximab response prediction. mSystems 3, e00188-17 (2018).
Kraal, L., Abubucker, S., Kota, K., Fischbach, M. A. & Mitreva, M. The prevalence of species and strains in the human microbiome: a resource for experimental efforts. PLoS ONE 9, e97279 (2014).
Shade, A. & Handelsman, J. Beyond the Venn diagram: the hunt for a core microbiome. Environ. Microbiol. 14, 4–12 (2012).
Li, K., Bihan, M. & Methé, B. A. Analyses of the stability and core taxonomic memberships of the human microbiome. PLoS ONE 8, e63139 (2013).
Steinway, S. N., Biggs, M. B., Loughran, T. P. Jr, Papin, J. A. & Albert, R. Inference of network dynamics and metabolic interactions in the gut microbiome. PLoS Comput. Biol. 11, e1004338 (2015).
Sung, J. et al. Global metabolic interaction network of the human gut microbiota for context-specific community-scale analysis. Nat. Commun. 8, 15393 (2017).
Greenblum, S., Turnbaugh, P. J. & Borenstein, E. Metagenomic systems biology of the human gut microbiome reveals topological shifts associated with obesity and inflammatory bowel disease. Proc. Natl Acad. Sci. USA 109, 594–599 (2012).
Burns, A. R. et al. Contribution of neutral processes to the assembly of gut microbial communities in the zebrafish over host development. ISME J. 10, 655–664 (2016).
Vega, N. M. & Gore, J. Stochastic assembly produces heterogeneous communities in the Caenorhabditis elegans intestine. PLoS Biol. 15, e2000633 (2017).
Hall, A. B., Tolonen, A. C. & Xavier, R. J. Human genetic variation and the gut microbiome in disease. Nat. Rev. Genet. 18, 690–699 (2017).
Nygaard, V., Rødland, E. A. & Hovig, E. Methods that remove batch effects while retaining group differences may lead to exaggerated confidence in downstream analyses. Biostatistics 17, 29–39 (2016).
MathSciNet PubMed Google Scholar
McFarland, L. V., Evans, C. T. & Goldstein, E. J. C. Strain-specificity and disease-specificity of probiotic efficacy: a systematic review and meta-analysis. Front. Med. 5, 124 (2018).
Segata, N. On the road to strain-resolved comparative metagenomics. mSystems 3, e00190-17 (2018).
Greenblum, S., Carr, R. & Borenstein, E. Extensive strain-level copy-number variation across human gut microbiome species. Cell 160, 583–594 (2015).
Faith, J. J., Colombel, J.-F. & Gordon, J. I. Identifying strains that contribute to complex diseases through the study of microbial inheritance. Proc. Natl Acad. Sci. USA 112, 633–640 (2015).
Zhang, C. & Zhao, L. Strain-level dissection of the contribution of the gut microbiome to human metabolic disease. Genome Med. 8, 41 (2016).
Tett, A. et al. Unexplored diversity and strain-level structure of the skin microbiome associated with psoriasis. NPJ Biofilms Microbiomes 3, 14 (2017).
Poretsky, R., Rodriguez-R, L. M., Luo, C., Tsementzi, D. & Konstantinidis, K. T. Strengths and limitations of 16S rRNA gene amplicon sequencing in revealing temporal microbial community dynamics. PLoS ONE 9, e93827 (2014).
Tessler, M. et al. Large-scale differences in microbial biodiversity discovery between 16S amplicon and shotgun sequencing. Sci. Rep. 7, 6589 (2017).
Tremlett, H., Bauer, K. C., Appel-Cresswell, S., Finlay, B. B. & Waubant, E. The gut microbiome in human neurological disease: a review. Ann. Neurol. 81, 369–382 (2017).
Griffiths, J. A. & Mazmanian, S. K. Emerging evidence linking the gut microbiome to neurologic disorders. Genome Med. 10, 98 (2018).
MacQueen, G., Surette, M. & Moayyedi, P. The gut microbiota and psychiatric illness. J. Psychiatry Neurosci. 42, 75–77 (2017).
Consortium, T. H. M. P., The Human Microbiome Project Consortium. Structure, function and diversity of the healthy human microbiome. Nature 486, 207–214 (2012).
Lloyd-Price, J. et al. Strains, functions and dynamics in the expanded Human Microbiome Project. Nature 550, 61–66 (2017).
Franzosa, E. A. et al. Species-level functional profiling of metagenomes and metatranscriptomes. Nat. Methods 15, 962–968 (2018).
Porras, A. M. & Brito, I. L. The internationalization of human microbiome research. Curr. Opin. Microbiol. 50, 50–55 (2019).
Wooley, J. C., Godzik, A. & Friedberg, I. A primer on metagenomics. PLoS Comput. Biol. 6, e1000667 (2010).
Willner, D. & Hugenholtz, P. From deep sequencing to viral tagging: recent advances in viral metagenomics. Bioessays 35, 436–442 (2013).
This work was supported, in part, by the Mayo Clinic Center for Individualized Medicine (to V.K.G., M.K., U.B., K.Y.C., N.C., and J.S.), and Mark E. and Mary A. Davis to Mayo Clinic Center for Individualized Medicine (J.S.).
Microbiome Program, Center for Individualized Medicine, Mayo Clinic, Rochester, MN, 55905, USA
Vinod K. Gupta, Minsuk Kim, Utpal Bakshi, Nicholas Chia & Jaeyun Sung
Division of Surgery Research, Department of Surgery, Mayo Clinic, Rochester, MN, 55905, USA
Graduate Research Education Program (GREP), Mayo Clinic, Rochester, MN, 55905, USA
Kevin Y. Cunningham
Department of Computer Science and Engineering, University of Minnesota Twin-Cities, Minneapolis, MN, 55455, USA
Division of Rheumatology, Department of Medicine, Mayo Clinic, Rochester, MN, 55905, USA
John M. Davis III & Jaeyun Sung
Division of Gastroenterology and Hepatology, Mayo Clinic College of Medicine and Science, Rochester, MN, 55905, USA
Konstantinos N. Lazaridis
Emeritus Chair, Department of Surgery, Mayo Clinic, Rochester, MN, 55905, USA
Heidi Nelson
Vinod K. Gupta
Minsuk Kim
Utpal Bakshi
John M. Davis III
Nicholas Chia
Jaeyun Sung
V.K.G. and J.S. conceived the problem. V.K.G., M.K., and J.S. designed all analytical methodologies. V.K.G. and K.Y.C. performed the computational experiments. All authors analyzed the data. V.K.G. and J.S. wrote the manuscript, with contributions from other authors. J.M.D. is the principal investigator of the Mayo Clinic Rheumatology Biobank, from which stool samples were collected from patients with rheumatoid arthritis. All authors reviewed and approved the final manuscript.
Correspondence to Jaeyun Sung.
V.K.G. and J.S. disclose that a patent application was filed relating to the materials in this manuscript. All other authors declare no competing interests.
Peer review information Nature Communications thanks Christian Diener, Scott Olesen and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.
Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Data 1
Description of Additional Supplementary Files
Gupta, V.K., Kim, M., Bakshi, U. et al. A predictive index for health status using species-level gut microbiome profiling. Nat Commun 11, 4635 (2020). https://doi.org/10.1038/s41467-020-18476-8
Aging gut microbiota of wild macaques are equally diverse, less stable, but progressively personalized
Baptiste Sadoughi
Dominik Schneider
Microbiome (2022)
Microbial risk score for capturing microbial characteristics, integrating multi-omics data, and predicting disease risk
Chan Wang
Leopoldo N. Segal
Huilin Li
Meta-analysis defines predominant shared microbial responses in various diseases and a specific inflammatory bowel disease signal
Haya Abbas-Egbariya
Yael Haberman
Amnon Amir
Genome Biology (2022)
Alterations in the respiratory tract microbiome in COVID-19: current observations and potential significance
Carter Merenstein
Frederic D. Bushman
Ronald G. Collman
Meta-analysis of microbiome association networks reveal patterns of dysbiosis in diseased microbiomes
Tony J. Lam
Yuzhen Ye
2020 Top 50 Life and Biological Sciences Articles