Next Article in Journal
SpoVG Is Necessary for Sporulation in Bacillus anthracis
Previous Article in Journal
Investigation of Heterologously Expressed Glucose-6-Phosphate Dehydrogenase Genes in a Yeast zwf1 Deletion
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Microorganisms
Volume 8
Issue 4
10.3390/microorganisms8040547
microorganisms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Investigation of Adiposity Measures and Operational Taxonomic unit (OTU) Data Transformation Procedures in Stool Samples from a German Cohort Study Using Machine Learning Algorithms

by
Martina Troll
1,2,*,
Stefan Brandmaier
1,2,
Sandra Reitmeier
3,4,
Jonathan Adam
1,2,
Sapna Sharma
1,2,
Alice Sommer
2,5,
Marie-Abèle Bind
5,
Klaus Neuhaus
3,
Thomas Clavel
3,6,
Jerzy Adamski
7,8,9,
Dirk Haller
3,4,
Annette Peters
2,10 and
Harald Grallert
1,2,11,*
1
Research Unit of Molecular Epidemiology, Helmholtz Zentrum München, 85764 Neuherberg, Germany
2
Institute of Epidemiology, Helmholtz Zentrum München, 85764 Neuherberg, Germany
3
ZIEL Institute for Food & Health, Technical University of Munich, 85354 Freising-Weihenstephan, Germany
4
Chair of Nutrition and Immunology, Technical University of Munich, 85354 Freising-Weihenstephan, Germany
5
Department of Statistics, Faculty of Arts and Sciences, Harvard University, Cambridge, MA 02138-2901, USA
6
Functional Microbiome Research Group, Institute of Medical Microbiology, RWTH University Hospital, 52074 Aachen, Germany
7
Research Unit Molecular Endocrinology and Metabolism, Helmholtz Zentrum München, 85764 Neuherberg, Germany
8
Chair of Experimental Genetics, Technical University of Munich, 85350 Freising-Weihenstephan, Germany
9
Department of Biochemistry, Yong Loo Lin School of Medicine, National University of Singapore, Singapore 117597, Singapore
10
Chair of Epidemiology, Faculty of Medicine, Ludwig-Maximilians-University München, 81377 Munich, Germany
11
German Center for Diabetes Research (DZD), 85764 Neuherberg, Germany
 
add Show full affiliation list
 
remove Hide full affiliation list
*
Authors to whom correspondence should be addressed.
Microorganisms 2020, 8(4), 547; https://doi.org/10.3390/microorganisms8040547
Submission received: 9 March 2020 / Revised: 1 April 2020 / Accepted: 7 April 2020 / Published: 10 April 2020
(This article belongs to the Section Systems Microbiology)
Browse Figures
Versions Notes

Abstract

:
The analysis of the gut microbiome with respect to health care prevention and diagnostic purposes is increasingly the focus of current research. We analyzed around 2000 stool samples from the KORA (Cooperative Health Research in the Region of Augsburg) cohort using high-throughput 16S rRNA gene amplicon sequencing representing a total microbial diversity of 2089 operational taxonomic units (OTUs). We evaluated the combination of three different components to assess the reflection of obesity related to microbiota profiles: (i) four prediction methods (i.e., partial least squares (PLS), support vector machine regression (SVMReg), random forest (RF), and M5Rules); (ii) five OTU data transformation approaches (i.e., no transformation, relative abundance without and with log-transformation, as well as centered and isometric log-ratio transformations); and (iii) predictions from nine measurements of obesity (i.e., body mass index, three measures of body shape, and five measures of body composition). Our results showed a substantial impact of all three components. The applications of SVMReg and PLS in combination with logarithmic data transformations resulted in considerably predictive models for waist circumference-related endpoints. These combinations were at best able to explain almost 40% of the variance in obesity measurements based on stool microbiota data (i.e., OTUs) only. A reduced loss in predictive performance was seen after sex-stratification in waist–height ratio compared to other waist-related measurements. Moreover, our analysis showed that the contribution of OTUs less prevalent and abundant is minor concerning the predictive power of our models.

1. Introduction

As a consequence of modern lifestyles, excessive body weight has become an important health burden, with a prevalence that has greatly surpassed underweight on a global scale [1]. Understanding biological mechanisms associated with obesity in combination with changes in the daily habits and living conditions of people will help address this challenge of current health systems and societies. One of the biological factors influencing obesity is the gut microbiome [2], which has also been linked to comorbidities, such as diabetes [3], and to even further pathophysiological processes, such as neurological and mental disorders [4]. These associations seem not surprising as extensive research has shown many axes between the gut microbiome and metabolically relevant organs such as the liver, adipose tissue, and even the brain [5]. Moreover, the gut microbiota has its role in the mediation and conversion of external input from the environment, including medication [6,7], food, and energy metabolism [5,8,9]. These findings introduced the gut microbiome as a key player for the comprehension of and possibly novel intervention strategies for the obese phenotype. However, there is no clear consensus on how the association between obesity and the gut microbiome is shaped [5,10,11]. For improved understanding of the interactions and derivation of correct conclusions, it is important to dissect data processing strategies and their impact on study outcomes [12,13].
Human gut microbiota are commonly analyzed using high-throughput 16S rRNA gene amplicon analysis of feces, an easily accessible and noninvasive sampling procedure for health assessments. Its ubiquity in bacteria and the combination of highly conserved and variable regions within the 16S rRNA gene makes it suitable for compositional analyses of the bacterial gut content [14]. During data processing, sequenced DNA fragments are usually clustered into operational taxonomic units (OTUs) based on sequence similarity above a determined threshold, though other analysis strategies are also emerging [15]. Despite the ability to rapidly obtain a comprehensive overview of microbiota profiles using amplicon sequencing, the resulting data is of inherent compositional nature due to technical instruments being limited in their maximum measurement capacities. The occurrence of OTUs is therefore expressed in relative abundances with inherent dependencies among them. In effect, an increase in counts measured of one OTU will cause a decrease in available measurement slots for other OTUs. However, it is assumed that the resulting counts reflect the underlying distribution in native samples to a high degree [12]. Further, many OTUs are found only in one or few samples, causing sparsity in the OTU table, i.e., great amounts of zeros, which poses an additional challenge in OTU data analysis [13].
Consequently, standard statistical procedures are likely not able to account for these constraints and lead to suboptimal data evaluation. Log-ratio transformations, in particular centered log-ratio (CLR) and isometric log-ratio (ILR), as preprocessing methods for OTU data have been discussed in coping with some of these issues, as they are invariable, whether they have been calculated from the measured count data or the unknown complete original sample. In the case of the CLR, the OTUs are displayed in relationship to the geometric mean of the sample, while ILR uses orthonormal bases [12,16].
In the present study, we evaluated the suitability of different methods to capture the complexity of a recently generated OTU data set in approx. 2000 individuals from the cross-sectional KORA (Cooperative Health Research in the Region of Augsburg) FF4 study (2013/14). Specifically, we investigated the practical consequences of data transformations in light of machine learning prediction. Moreover, we compared nine measures of obesity and how well they could be represented by the fecal microbiota profiles.

2. Materials and Methods

2.1. Study Population

The analyses were based on participants’ data from the KORA (Cooperative Health Research in the Region of Augsburg) FF4 study. As the 2nd population-based follow-up examination of the KORA S4 study (1999/2001), it was conducted in 2013/2014 with 2279 individuals in southern Germany. The investigations were carried out in accordance with the Declaration of Helsinki, including written informed consent of all participants. All study methods were approved by the ethics committee of the Bavarian Chamber of Physicians, Munich (FF4: EC No. 06068, 25 October 2012) [17,18]. The informed consent given by KORA study participants does not cover data posting in public databases. However, data are available upon request from KORA-gen (http://epi.helmholtz-muenchen.de/kora-gen/) by means of a project agreement. Requests should be sent to [email protected] and are subject to approval by the KORA Board.

2.2. Microbiota Profiling by 16S rRNA Amplicon Sequencing

The microbiome measurement procedures via high-throughput 16S rRNA gene sequencing have been described in detail elsewhere [19]. In brief, of all 2272 participants from the KORA FF4 study, individuals with missing consent were excluded (n = 197). The remaining participants were asked to collect a stool sample in a collection tube containing DNA stabilizer (Stratec DNA Stool Stabilizer, Cat. no. 1038111100), which were finally stored at −80 °C. After DNA extraction, the V3-V4 regions of the 16S rRNA gene were amplified via 2-step PCR (primers 341F-ovh and 785r-ovh) and sequenced on Illumina HiSeq. To guarantee high-quality output, samples with low read counts were resequenced on an Illumina MiSeq. Raw sequence reads were processed using the IMNGS platform [20] that integrates a UPARSE-based OTU clustering approach at 97% sequence similarity level [21]. Only OTUs occurring at a relative abundance ≥0.25% in at least one sample across the entire cohort study were kept to avoid the analysis of spurious taxa [22]. For technical quality control, samples with an insufficient amount of DNA content (<12 ng/μL) or reads (<7000) were excluded as well (n = 42). In total, 2089 OTUs were identified for the complete study with usable fecal samples.

2.3. Exclusion Criteria

Of the available 2033 individuals with fecal microbiome measurements, participants with antibiotic medication intake (n = 41), malignant neoplasms of digestive organs according to WHO ICD-10 codes C15-C26 (n = 27), and underweight with a BMI < 18.5 kg/m2 according to WHO criteria (n = 8) were excluded. All individuals with missing information in one of the adiposity measures were additionally excluded (n = 35), resulting in a study population of 1923 individuals.

2.4. Adiposity Measures

Nine adiposity measures were examined in this study: body mass index (BMI; kg/m2), waist circumference (WC; cm), waist–hip ratio (WHR), waist–height ratio (WHtR), body adiposity index (BAI; according to the formula by Reference [23]: (hip circumference (cm))/((height (cm) × 0.01)1.5) − 18), fat mass index (FMI; kg/m2), lean body mass index (LBMI; kg/m2), appendicular muscle mass index (AMMI; kg/m2), and body fat (BF; %). Body composition was measured via bioelectrical impedance analysis (BIA 2000-S device, Data Input, Pöcking, Germany) with Kyle’s equations [24,25,26,27]. To enable comparability of results, all adiposity measures were scaled to a mean of 0 and standard deviation of 1 for each (sub-)set for machine learning, separately.

2.5. Machine Learning Algorithms

In this study, four different machine learning algorithms, i.e., support vector machine regression (SVMReg) with normalized poly kernel (NPK), random forest (RF), and M5Rules as well as partial least squares (PLS) with different numbers of principal components, 20 (PLS20) and 4 (PLS4), were compared. All prediction models were calculated in Weka version 3.8.3, except for PLS in Weka version 3.6.15 [28,29], using default options except the following parametrizations: M5rules: minNumInstances 4.0; SVM: filter type no normalization/standardization; and NPK: E2.0. From the Weka summary statistics, the correlation coefficients (CC) were used for model evaluation. They represent the Pearson correlation between the actual and the predicted values of the outcome (i.e., obesity measurement) from 10-fold cross-validation.

2.6. Statistical Analysis

The remaining analyses were conducted in R version 3.6.0 [30]. For each (sub-)set of the data, data transformations were executed separately and each transformation was calculated sample-wise. Relative abundance (RA) was calculated with each sample sum to 100%. Logarithm (Log) refers to the natural logarithm. Centered log-ratio (CLR) and isometric log-ratio (ILR) were generated with the compositions package [31] in R. Raw count data was included for further analyses as the control value. Alpha diversity as Simpson effective counts according to Jost [32] was calculated via the vegan package in R [33]. For comparison of prediction rates with and without sparsely occurring OTUs, across-sample prevalence and abundance were calculated for the complete OTU data set and respective OTUs were excluded in a stepwise fashion. For prevalence, OTUs present in less than x% of the samples were excluded; for abundance, OTUs with a relative abundance equal to or less than x% were excluded. The term prevalence was used in this study to describe the number of times an OTU was nonzero in the samples.

3. Results and Discussion

3.1. Benchmark Data on Study Population and Adiposity Measures

First, we examined the features of nine adiposity measures within our study population. On average, study participants were older than the German median age [34] as well as overweight according to the World Health Organization (WHO) guidelines [35] (Table 1). Three different components of adiposity assessment were examined, including representation of overall obesity (i.e., BMI), body shape (i.e., waist circumference, waist–hip ratio, and waist–height ratio), as well as body composition measurements with its two subcomponents of muscle-related measures (i.e., lean body mass index and appendicular muscle mass index) and fat-related measures (i.e., fat mass index, body fat percentage, and body adiposity index). Each measure exhibited a large range, approximating a comprehensive coverage of possible forms for each phenotype, except underweight (Figure 1A and Table S1). For example, the lowest BMI measurement in our data was 18.6 kg/m2 (low normal weight) and the highest was at 62.2 kg/m2 (severely obese).
On the topic of sex differences in adiposity measures, three aspects can be noted from the present data (Figure 1A): (i) The distribution of most measures related to body shape as well as body composition showed distinct sex-specific differences in their distribution. (ii) In certain measures, women showed dual peaks in their distribution of obesity, particularly in waist-related measurements. This was likely due to the simultaneous existence of android and gynoid fat distributions within adult female populations. The former might in part be due to menopause in a slightly older study population. This is in contrast to adult men with a predominant android phenotype [36,37]. (iii) On the other hand, accounting for height was able to reduce sex-related biases in some indices, such as in BMI and waist–height ratio.
Additionally, three aspects on the relation of adiposity measures among each other are set out from the comparison of Spearman correlations (Figure 1B): (i) The same two indices, BMI and waist–height ratio, both showed reasonable correlations to all other adiposity measurements, ranging from 0.51 to 0.90 for BMI and 0.49 to 0.91 for waist–height ratio. (ii) Closely related measurements, i.e., within one (sub-)component of adiposity assessment, demonstrated particularly high correlations: e.g., lean body mass index and appendicular muscle mass index (0.99), fat mass index and body fat percentage (0.91), or waist circumference and waist–height ratio (0.91). (iii) Furthermore, measurements from one closely related group sometimes exhibited negative correlations to measurements from another group, e.g., appendicular muscle mass index and body fat percentage (−0.25) or waist–hip ratio and body adiposity index (−0.10).

3.2. Adiposity Measurement Selection

We trained prediction models for obesity using adiposity measurements as outcome and OTU data as predictors. We calculated models for all possible combinations of body measurements, preprocessing procedures, and machine learning algorithms (Table 2). For evaluation, we used correlation coefficients (CCs) from cross-validation between actual and predicted values of the outcome (see Section 2, Materials and Methods), with a potential value of 1 at its best and of 0 at its worst (as a negative CC would indicate a model, which predicts the opposite of the measured values).
Our results revealed the lowest prediction rates for measurements related to body fat (i.e., body adiposity index, fat mass index, and body fat percentage). CCs for these measurements varied from 0.04 (body fat percentage in combination with relative abundance and M5Rules) to 0.28 (fat mass index and body fat percentage with support vector machine regression (SVMReg) and log-ratio transformations). By comparison, the highest prediction rates were achieved with measures related to central fat accumulation (i.e., waist circumference, waist–hip ratio, and waist–height ratio). These measurements were able to predict the variation in the phenotypes up to almost 40%. Additionally, our results showed that measurements related to muscle mass (i.e., appendicular muscle mass index and lean body mass index) were also well predicted by fecal OTU data. Specifically, the explained variation in muscle-related measurements of our study reached 36% and 35% for lean body mass index and appendicular muscle mass index, respectively.
These findings may lead to several indications. (i) Our data suggested a link between the gut microbiota profiles and central fat rather than body fat content in general. These findings are in line with previous studies that identified a stronger relationship of central obesity measurements than to other obesity measurements such as BMI. In particular, visceral fat has shown associations with the gut microbiome and this could be the reason behind the higher correlation towards waist measurements compared to other parameters. For example, Zierer et al. speculated that the association of the gut microbiome to fatty acid metabolism “may be better reflected by visceral fat measures” [38] as well as an interrelation with food to influence visceral fat [38,39,40,41]. (ii) The prediction capacity in muscle measurements supported recently published experiments in mice that demonstrated effects of the gut microbiota on the quantity, mass, and functional characteristics of skeletal muscles [42]. (iii) Taken together, these findings might suggest a link of gut microbiota composition and cardiometabolic health. Central obesity, for example, high waist–height ratio, has been identified as a risk factor for cardiovascular disease mortality [43], whereas muscle mass seems to protect from cardiovascular disease [44]. Additionally, previous studies already have shown associations of the gut microbiome with cardiometabolic risk factors [45] and are supported by mechanistic studies in mice, suggesting a role for gut bacteria in dietary energy harvest and fatty acid metabolism [41]. These findings should convey further investigations, examining the role of the gut microbiome in cardiometabolic health. Among other things, this observational study could be embedded in a hypothetical randomized experiment, as suggested by Bind and Rubin [46].

3.3. Comparison of OTU Data Transformation Approaches

Because of ongoing debates on optimal data normalization, transformation, and general preprocessing [13,47], we aimed to evaluate different methods within our data set and in combination with our target outcomes, i.e., obesity phenotypes (Table 2). Logarithmic transformations (i.e., sample-wise relative abundance with log-transformation, CLR, and ILR) generally improved the prediction performance independent of outcome or machine learning algorithm, except for random forest (RF). The best results were achieved with waist measurements, in particular in combination with SVMReg and CLR (optimal CCs for each waist measurement between 0.36 and 0.39) and PLS4 in all log-transformations (0.35–0.37). This was underlined when comparing RA values with and without a log-transformation. This ensured a distinction between different methods and the improvement resulting from employing a logarithmic transformation as such [48]. These results supported the assumption of logarithmic transformations as a useful tool in microbiota data analysis and from a machine learning perspective in particular.
Comparing the different transformation methods that include a logarithmic transformation, the results suggested only minor differences between the methods in terms of practical application to OTU data. Although CLR and ILR achieved the overall best results in combination with waist circumference, they did not lead to a big improvement in performance in comparison to log-transformed RAs. For our subsequent analyses with SVMReg, the easier interpretation of CLR values as a relation to the geometric mean made it preferable to ILR.
Log-transformations have been critically discussed for their potential to increase the impact of less prevalent and abundant OTUs and vice versa for dominant OTUs in certain contexts [47]. Our analysis conversely suggests an improvement in performance with logarithmic transformations in most cases, with no emphasis on sparsely occurring OTUs in the prediction of waist–height ratio (Table 2, Table S2, and Figure 2). Moreover, log-transformations are able to manage skewedness and to limit the large variation in raw count data and already lead to major improvements in predictions in our data set. Therefore, we assumed that difficulties for many algorithms in coping with said variation, and in particular SVMReg and PLS4, benefited from log transformations. For example, for waist–height ratio predictions, the CC improved from 0.21 (SVMReg) and 0.23 (PLS4) in raw count data up to 0.37 (SVMReg) and 0.36 (PLS4) in logarithmically transformed data. In contrast, an algorithm based on randomness like random forest (RF) is able to work comparatively well with raw count data and RAs and even almost consistently better than with logarithmically transformed data. In the same example as before, RF achieved a correlation of 0.31 in raw count data but of 0.25 in CLR.

3.4. Implications from Machine Learning Algorithms

The juxtaposition of machine learning algorithms contained mixed messages. On the one hand, algorithms without assumptions on data structure, such as SVMReg and RF, showed good overall performance, which may underline the importance of accounting for the complex nature of OTU data. For instance, OTUs are clusters of 16S rRNA gene sequences using a similarity threshold without possible functional differentiation. Thus, an OTU may contain different strains of the same species or species which are taxonomically very close. Consequently, such groupings may contain strains even with opposite effects on a phenotype [49,50]. In addition, OTUs as representation of bacteria in the human gut likely do not exist as isolated entities but have many interrelations among each other [51,52]. In particular SVMReg, which transforms into higher dimensional feature spaces via the use of kernels, seemed well-suited to account for such interdependencies of OTUs.
On the other hand, linear models, here exemplified by PLS, performed similarly well on the OTU data if used with a low number of principal components. This was in particular the case after logarithmic transformation, which lead to a normality-like distribution optimal for a linear approach. In this light, our data supports a linear dependency between logarithmically transformed OTU data and adiposity endpoints in our data set.
The results concerning the number of principal components also raised the question of data overfitting. A high number of components (PLS20) consistently lead to worse performance, while a low number (PLS4) represented one of the best performing methods used. Therefore, the number of parameters contained in a linear model should be critically assessed and kept to a necessary minimum.

3.5. Effect of Sex-Stratification on Predictive Performance

In a next step, we examined prediction rates of abdominal measurements in a sex-stratified approach using two of the best performing combinations (i.e., SVMReg with CLR and PLS4 with log-transformed RAs). Of the three abdominal measurements compared in our study, waist–height ratio had the lowest decrease in predictive performance in the sex-stratified data sets in both methods. With PLS4, the CCs in waist–height ratio decreased from 0.36 in the total population to 0.28 and 0.34 in men and women, respectively. By comparison, with waist circumference the CCs dropped from 0.37 to 0.27 (men) and 0.31 (women). This indicates that the sex-specific variation in waist circumference-related measures could be reasonably addressed by accounting for height (Figure 1A and Table 3). Therefore, waist–height ratio combined the advantages of BMI, i.e., height adjustment and reasonable correlations to other adiposity measurements (Figure 1B), with the improved prediction capacity of waist measurements.
However, generally, it could be observed that the sex stratification reduced prediction capabilities. To attribute this effect to the smaller sample size seems implausible, as the sensitivity of machine learning, in contrast to p-values, is unlikely to be impaired in population sizes that vary in the magnitude of almost 1000 to 2000 (Table 1 and Table 3). Still, a data split could cause particularly sex-specific OTUs to be less prominent in either of the strata and, hence, considered to be less informative by the algorithm than in the complete data set. This should warrant specific attention to sex differences and appropriate analysis approaches for future studies.
However, the effects on the predictive performance were also distributed unevenly among the sexes, whereby the phenotypes are consistently better predicted in women. One reason for this might be the slightly narrowed distribution of waist–height ratio values in men compared to women (Figure 1A). Another factor might be the, on average, increased ecosystem diversity in women (Figure S1). In a previous publication on the dietary habits of the KORA FF4 population, Breuninger et al. showed that women consume a more fiber-dense diet than men [53] and, hence, a rather favorable nutrition for the gut microbiome linked to increased gut bacterial diversity.

3.6. Dependency of Prediction Capacity on Prevalent and Abundant OTUs

In order to analyze the contribution of low abundant/prevalent OTUs to the predictive performance of our models, we excluded OTUs in a stepwise fashion according to their across-sample prevalence and abundance in the data set. For this purpose, we used the same two combinations as methods (i.e., SVMReg with CLR and PLS4 with log-transformed RAs) to predict waist–height ratio. However, our results suggested that the prediction of waist–height ratio was mainly based on frequent and abundant OTUs (Figure 2 and Table S2). We observed the same prediction capacity after the exclusion of OTUs which occurred with a prevalence of less than 30% of the samples, which make up more than 72% of all OTUs in the data. Even after the exclusion of 98% of all OTUs, the top 52 entities were able to explain almost a quarter of the variation in waist–height ratio in the population as estimated by PLS4. Similar approaches have been implemented in previous studies, e.g., using a threshold of 25% [38,39]. Within our analyses of waist–height ratio, the same threshold would merely lead to a decrease in prediction capacity by approximately 0.01–0.03 with SVMReg and no decrease in PLS4. The same could be observed for across-sample relative abundance. An exclusion of OTUs with an abundance of 0.01% or less, which excluded close to 50% of OTUs, led to almost no (−0.01, SVMReg) to no decrease (PLS4) in correlation coefficients. In a scenario where the quality of prediction is mainly determined by prevalent and abundant OTUs, it seems plausible that the effects of pseudocounts, an added value to offset zeros before log-transformation, was only minimal as shown in our analyses (Table S3). This is relevant to the criticism of log-ratio transformations, in that they need coping mechanisms for the zero-inflation in OTU data sets and of which pseudocounts are one possible choice [13].

3.7. Limitations

This study has several limitations. First, the study population reflected a middle-aged to elderly part of the population (38–88 years, Table S1). Therefore, the findings cannot be extrapolated to all age groups. Moreover, causal relationships cannot be established in this kind of study, and further validation from experimental data is needed. In addition, the tests were restricted to a dataset from one cohort study, and therefore, these results on their own cannot be generalized. The cross-sectional nature of our study only permits the assessment of correlations of a specific point in time for the participants. Recent changes in weight cannot be excluded but should average out with sufficient sample size.

4. Conclusions

In summary, we applied a selection of preexisting methods to compare the relation between adiposity measures and OTU data for the KORA FF4 study specifically. We identified waist-related obesity measures to be reflected best by our OTU data. In particular waist–height ratio, which showed less bias in sex-stratification compared to the total population, seemed to be the most appropriate choice. Our analyses showed that almost 40% of the variance in abdominal obesity could be reflected by composition of fecal microbiota. The predictive ability is mostly guaranteed by prevalent and abundant OTUs.

Supplementary Materials

The following are available online at https://www.mdpi.com/2076-2607/8/4/547/s1, Figure S1: Alpha diversity as effective number of species by Simpson index in men and women; Table S1: Extended benchmark data on total study population (n = 1923); Table S2a. Comparison of levels of prevalence in waist–height ratio; Table S2b. Comparison of levels of relative abundance; Table S3. Comparison of pseudocounts; Table S4: OTUs.

Author Contributions

Conceptualization, M.T., S.B., A.P., and H.G.; methodology, M.T., S.B., and T.C.; software, M.T. and S.B.; validation, M.T. and J.A. (Jonathan Adam); formal analysis, M.T.; investigation, M.T.; resources, S.R., K.N., T.C., D.H., A.P., and H.G.; data curation, S.R. and T.C.; writing—original draft preparation, M.T. and S.B.; writing—review and editing, M.T., S.B., S.R., J.A. (Jonathan Adam), S.S., A.S., M.-A.B., K.N., T.C., J.A. (Jerzy Adamski), D.H., A.P., and H.G.; visualization, M.T.; supervision, A.P. and H.G.; project administration, H.G.; funding acquisition, T.C., D.H., and A.P. All authors have read and agreed to the published version of the manuscript.

Funding

The KORA study was initiated and financed by the Helmholtz Zentrum München—German Research Center for Environmental Health, which is funded by the German Federal Ministry of Education and Research (BMBF) and by the State of Bavaria. Furthermore, KORA research was supported within the Munich Center of Health Sciences (MC-Health), Ludwig-Maximilians-Universität, as part of LMUinnovativ. The Technical University of Munich provided funding for the ZIEL Institute for Food & Health and Caroline Ziegler provided technical support for sample preparation and 16S rRNA gene amplicon sequencing. Microbiota profiling of KORA samples was supported by enable Kompetenzcluster der Ernährungsforschung (No. 01EA1409A) and the European Union Joint Programming Initiative DINAMIC (No. 2815ERA04E, 2815ERA11E). T.C. received also funding from the German Research Foundation (DFG): Project-ID 403224013—SFB 1382. Research reported in this publication was also supported by the John Harvard Distinguished Science Fellow Program within the FAS Division of Science of Harvard University and by the Office of the Director, National Institutes of Health under award number DP5OD021412. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.

Acknowledgments

We thank all KORA participants and technical assistants without whose contributions this study could not have been realized.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. World Health Organization Malnutrition. Key Facts. 2018. Available online: https://www.who.int/news-room/fact-sheets/detail/malnutrition (accessed on 20 March 2020).
  2. Khan, M.J.; Gerasimidis, K.; Edwards, C.A.; Shaikh, M.G. Role of Gut Microbiota in the Aetiology of Obesity: Proposed Mechanisms and Review of the Literature. J. Obes. 2016, 2016, 1–27. [Google Scholar] [CrossRef] [Green Version]
  3. Hartstra, A.V.; Bouter, K.E.C.; Bäckhed, F.; Nieuwdorp, M. Insights into the Role of the Microbiome in Obesity and Type 2 Diabetes. Diabetes Care 2015, 38, 159–165. [Google Scholar] [CrossRef] [Green Version]
  4. Fung, T.C.; Olson, C.A.; Hsiao, E.Y. Interactions between the microbiota, immune and nervous systems in health and disease. Nat. Neurosci. 2017, 20, 145–155. [Google Scholar] [CrossRef]
  5. Schroeder, B.O.; Bäckhed, F. Signals from the gut microbiota to distant organs in physiology and disease. Nat. Med. 2016, 22, 1079–1089. [Google Scholar] [CrossRef]
  6. Jackson, M.A.; Verdi, S.; Maxan, M.-E.; Shin, C.M.; Zierer, J.; Bowyer, R.C.E.; Martin, T.; Williams, F.M.K.; Menni, C.; Bell, J.T.; et al. Gut microbiota associations with common diseases and prescription medications in a population-based cohort. Nat. Commun. 2018, 9, 2655. [Google Scholar] [CrossRef] [Green Version]
  7. Wu, H.; Esteve, E.; Tremaroli, V.; Khan, M.T.; Caesar, R.; Mannerås-Holm, L.; Ståhlman, M.; Olsson, L.M.; Serino, M.; Planas-Fèlix, M.; et al. Metformin alters the gut microbiome of individuals with treatment-naive type 2 diabetes, contributing to the therapeutic effects of the drug. Nat. Med. 2017, 23, 850–858. [Google Scholar] [CrossRef]
  8. Tseng, C.-H.; Wu, C.-Y. The gut microbiome in obesity. J. Formos. Med. Assoc. 2019, 118, S3–S9. [Google Scholar] [CrossRef]
  9. Bohan, R.; Tianyu, X.; Tiantian, Z.; Ruonan, F.; Hongtao, H.; Qiong, W.; Chao, S. Gut microbiota: a potential manipulator for host adipose tissue and energy metabolism. J. Nutr. Biochem. 2019, 64, 206–217. [Google Scholar] [CrossRef]
  10. 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 2014, 9, e84689. [Google Scholar] [CrossRef] [Green Version]
  11. Sze, M.A.; Schloss, P.D. Looking for a Signal in the Noise: Revisiting Obesity and the Microbiome. mBio 2016, 7, e01018-16. [Google Scholar] [CrossRef] [Green Version]
  12. Gloor, G.B.; Macklaim, J.M.; Pawlowsky-Glahn, V.; Egozcue, J.J. Microbiome Datasets Are Compositional: And This Is Not Optional. Front. Microbiol. 2017, 8, 2224. [Google Scholar] [CrossRef] [Green Version]
  13. Weiss, S.; Xu, Z.Z.; Peddada, S.; Amir, A.; Bittinger, K.; Gonzalez, A.; Lozupone, C.; Zaneveld, J.R.; Vázquez-Baeza, Y.; Birmingham, A.; et al. Normalization and microbial differential abundance strategies depend upon data characteristics. Microbiome 2017, 5, 27. [Google Scholar] [CrossRef] [Green Version]
  14. Morgan, X.C.; Huttenhower, C. Chapter 12: Human Microbiome Analysis. PLoS Comput. Biol. 2012, 8, e1002808. [Google Scholar] [CrossRef] [Green Version]
  15. Callahan, B.J.; McMurdie, P.J.; Rosen, M.J.; Han, A.W.; Johnson, A.J.A.; Holmes, S.P. DADA2: High-resolution sample inference from Illumina amplicon data. Nat. Methods 2016, 13, 581–583. [Google Scholar] [CrossRef] [Green Version]
  16. Egozcue, J.J. Isometric Logratio Transformations for Compositional Data Analysis. Math. Geol. 2003, 35, 279–300. [Google Scholar] [CrossRef]
  17. Holle, R.; Happich, M.; Löwel, H.; Wichmann, H.E.; MONICA/KORA Study Group. KORA—A research platform for population based health research. Gesundheitswesen 2005, 67, S19–S25. [Google Scholar] [CrossRef] [Green Version]
  18. Laxy, M.; Knoll, G.; Schunk, M.; Meisinger, C.; Huth, C.; Holle, R. Quality of Diabetes Care in Germany Improved from 2000 to 2007 to 2014, but Improvements Diminished since 2007. Evidence from the Population-Based KORA Studies. PLoS ONE 2016, 11, e0164704. [Google Scholar] [CrossRef]
  19. Reitmeier, S.; Kießling, S.; Clavel, T.; List, M.; Almeida, E.L.; Ghosh, T.S.; Neuhaus, K.; Grallert, H.; Troll, M.; Rathmann, W.; et al. Arrhythmic gut microbiome signatures for risk profiling of Type-2 Diabetes. bioRxiv 2019. [Google Scholar] [CrossRef]
  20. Lagkouvardos, I.; Joseph, D.; Kapfhammer, M.; Giritli, S.; Horn, M.; Haller, D.; Clavel, T. IMNGS: A comprehensive open resource of processed 16S rRNA microbial profiles for ecology and diversity studies. Sci. Rep. 2016, 6, 33721. [Google Scholar] [CrossRef]
  21. Edgar, R.C. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 2013, 10, 996–998. [Google Scholar] [CrossRef]
  22. Reitmeier, S.; Hitch, T.C.; Fikas, N.; Hausmann, B.; Ramer-Tait, A.E.; Neuhaus, K.; Berry, D.; Haller, D.; Lagkouvardos, I.; Clavel, T. Handling of spurious sequences affects the outcome of high-throughput 16S rRNA gene amplicon profiling. Review 2020. [Google Scholar] [CrossRef]
  23. Bergman, R.N.; Stefanovski, D.; Buchanan, T.A.; Sumner, A.E.; Reynolds, J.C.; Sebring, N.G.; Xiang, A.H.; Watanabe, R.M. A Better Index of Body Adiposity. Obesity 2011, 19, 1083–1089. [Google Scholar] [CrossRef]
  24. Markus, M.R.P.; Rospleszcz, S.; Ittermann, T.; Baumeister, S.E.; Schipf, S.; Siewert-Markus, U.; Lorbeer, R.; Storz, C.; Ptushkina, V.; Peters, A.; et al. Glucose and insulin levels are associated with arterial stiffness and concentric remodeling of the heart. Cardiovasc. Diabetol. 2019, 18, 145. [Google Scholar] [CrossRef]
  25. Kyle, U.G.; Genton, L.; Karsegard, L.; Slosman, D.O.; Pichard, C. Single prediction equation for bioelectrical impedance analysis in adults aged 20--94 years. Nutrition 2001, 17, 248–253. [Google Scholar] [CrossRef]
  26. Kyle, U.G.; Schutz, Y.; Dupertuis, Y.M.; Pichard, C. Body composition interpretation. Contributions of the fat-free mass index and the body fat mass index. Nutrition 2003, 19, 597–604. [Google Scholar] [CrossRef]
  27. Kyle, U.G.; Genton, L.; Hans, D.; Pichard, C. Validation of a bioelectrical impedance analysis equation to predict appendicular skeletal muscle mass (ASMM). Clin. Nutr. 2003, 22, 537–543. [Google Scholar] [CrossRef]
  28. Frank, E.; Hall, M.A.; Witten, I.H. Online Appendix for Data mining: Practical Machine Learning Tools and Techniques. 2016. Available online: https://waikato.github.io/weka-wiki/citing_weka/ (accessed on 20 March 2020).
  29. Hall, M.; Frank, E.; Holmes, G.; Pfahringer, B.; Reutemann, P.; Witten, I.H. The WEKA Data Mining Software: An Update. SIGKDD Explor. 2009, 11, 10–18. [Google Scholar] [CrossRef]
  30. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2019; Available online: https://www.R-project.org/ (accessed on 21 March 2020).
  31. van den Boogaart, K.G.; Tolosana-Delgado, R.; Bren, M. compositions: Compositional Data Analysis. 2018. Available online: https://CRAN.R-project.org/package=compositions (accessed on 21 March 2020).
  32. Jost, L. Entropy and diversity. Oikos 2006, 113, 363–375. [Google Scholar] [CrossRef]
  33. Oksanen, J.; Blanchet, F.G.; Friendly, M.; Kindt, R.; Legendre, P.; McGlinn, D.; Minchin, P.R.; O’Hara, R.B.; Simpson, G.L.; Solymos, P.; et al. vegan: Community Ecology Package. 2019. Available online: https://CRAN.R-project.org/package=vegan (accessed on 21 March 2020).
  34. United Nations, Department of Economic and Social Affairs, Population Division. World Population Ageing 2017 – Highlights; United Nations: New York, NY, USA, 2017; Available online: https://www.un.org/en/development/desa/population/publications/pdf/ageing/WPA2017_Highlights.pdf (accessed on 21 March 2020).
  35. World Health Organization. Obesity: Preventing and Managing the Global Epidemic: Report of a WHO Consultation (WHO Technical Report Series 894); World Health Organization: Geneva, Switzerland, 2000; p. 9 (Table 2.1). Available online: https://www.who.int/nutrition/publications/obesity/WHO_TRS_894/en/ (accessed on 19 March 2020).
  36. Ley, C.J.; Lees, B.; Stevenson, J.C. Sex- and menopause-associated changes in body-fat distribution. Am. J. Clin. Nutr. 1992, 55, 950–954. [Google Scholar] [CrossRef] [Green Version]
  37. Kirchengast, S.; Gruber, D.; Sator, M.; Hartmann, B.; Knogler, W.; Huber, J. Menopause-associated differences in female fat patterning estimated by dual-energy X-ray absorptiometry. Ann. Hum. Biol. 1997, 24, 45–54. [Google Scholar] [CrossRef]
  38. Zierer, J.; Jackson, M.A.; Kastenmüller, G.; Mangino, M.; Long, T.; Telenti, A.; Mohney, R.P.; Small, K.S.; Bell, J.T.; Steves, C.J.; et al. The fecal metabolome as a functional readout of the gut microbiome. Nat. Genet. 2018, 50, 790–795. [Google Scholar] [CrossRef]
  39. Beaumont, M.; Goodrich, J.K.; Jackson, M.A.; Yet, I.; Davenport, E.R.; Vieira-Silva, S.; Debelius, J.; Pallister, T.; Mangino, M.; Raes, J.; et al. Heritable components of the human fecal microbiome are associated with visceral fat. Genome Biol. 2016, 17, 189. [Google Scholar] [CrossRef] [Green Version]
  40. Le Roy, C.I.; Bowyer, R.C.E.; Castillo-Fernandez, J.E.; Pallister, T.; Menni, C.; Steves, C.J.; Berry, S.E.; Spector, T.D.; Bell, J.T. Dissecting the role of the gut microbiota and diet on visceral fat mass accumulation. Sci. Rep. 2019, 9, 9758. [Google Scholar] [CrossRef]
  41. Arora, T.; Bäckhed, F. The gut microbiota and metabolic disease: current understanding and future perspectives. J. Intern. Med. 2016, 280, 339–349. [Google Scholar] [CrossRef]
  42. Lahiri, S.; Kim, H.; Garcia-Perez, I.; Reza, M.M.; Martin, K.A.; Kundu, P.; Cox, L.M.; Selkrig, J.; Posma, J.M.; Zhang, H.; et al. The gut microbiota influences skeletal muscle mass and function in mice. Sci. Transl. Med. 2019, 11, eaan5662. [Google Scholar] [CrossRef] [Green Version]
  43. Rost, S.; Freuer, D.; Peters, A.; Thorand, B.; Holle, R.; Linseisen, J.; Meisinger, C. New indexes of body fat distribution and sex-specific risk of total and cause-specific mortality: a prospective cohort study. BMC Public Health 2018, 18, 427. [Google Scholar] [CrossRef] [Green Version]
  44. Tyrovolas, S.; Panagiotakos, D.; Georgousopoulou, E.; Chrysohoou, C.; Tousoulis, D.; Haro, J.M.; Pitsavos, C. Skeletal muscle mass in relation to 10 year cardiovascular disease incidence among middle aged and older adults: the ATTICA study. J. Epidemiol. Community Health 2020, 74, 26–31. [Google Scholar] [CrossRef] [Green Version]
  45. Fei, N.; Bernabé, B.P.; Lie, L.; Baghdan, D.; Bedu-Addo, K.; Plange-Rhule, J.; Forrester, T.E.; Lambert, E.V.; Bovet, P.; Gottel, N.; et al. The human microbiota is associated with cardiometabolic risk across the epidemiologic transition. PLoS ONE 2019, 14, e0215262. [Google Scholar] [CrossRef] [Green Version]
  46. Bind, M.-A.C.; Rubin, D.B. Bridging observational studies and randomized experiments by embedding the former in the latter. Stat. Methods Med. Res. 2019, 28, 1958–1978. [Google Scholar] [CrossRef] [Green Version]
  47. McKnight, D.T.; Huerlimann, R.; Bower, D.S.; Schwarzkopf, L.; Alford, R.A.; Zenger, K.R. Methods for normalizing microbiome data: An ecological perspective. Methods Ecol. Evol. 2019, 10, 389–400. [Google Scholar] [CrossRef]
  48. Costea, P.I.; Zeller, G.; Sunagawa, S.; Bork, P. A fair comparison. Nat. Methods 2014, 11, 359. [Google Scholar] [CrossRef]
  49. Nguyen, N.-P.; Warnow, T.; Pop, M.; White, B. A perspective on 16S rRNA operational taxonomic unit clustering using sequence similarity. NPJ Biofilms Microbiomes 2016, 2, 16004. [Google Scholar] [CrossRef] [Green Version]
  50. Human Microbiome Project Consortium Structure, function and diversity of the healthy human microbiome. Nature 2012, 486, 207–214. [CrossRef] [Green Version]
  51. Venturelli, O.S.; Carr, A.C.; Fisher, G.; Hsu, R.H.; Lau, R.; Bowen, B.P.; Hromada, S.; Northen, T.; Arkin, A.P. Deciphering microbial interactions in synthetic human gut microbiome communities. Mol. Syst. Biol. 2018, 14, e8157. [Google Scholar] [CrossRef]
  52. Faust, K.; Sathirapongsasuti, J.F.; Izard, J.; Segata, N.; Gevers, D.; Raes, J.; Huttenhower, C. Microbial Co-occurrence Relationships in the Human Microbiome. PLoS Comput. Biol. 2012, 8, e1002606. [Google Scholar] [CrossRef]
  53. Breuninger, T.A.; Riedl, A.; Wawro, N.; Rathmann, W.; Strauch, K.; Quante, A.; Peters, A.; Thorand, B.; Meisinger, C.; Linseisen, J. Differential associations between diet and prediabetes or diabetes in the KORA FF4 study. J. Nutr. Sci. 2018, 7, e34. [Google Scholar] [CrossRef] [Green Version]
Microorganisms 08 00547 g001 550
Figure 1. Comparison of adiposity measures: (A) Sex-stratified distribution of adiposity measures and (B) Spearman correlation in study population (n = 1923). Adiposity measurements were z-transformed (mean = 0, standard deviation = 1). AMMI: appendicular muscle mass index; BAI: body adiposity index; BF (%): body fat percentage; BMI: body mass index; FMI: fat mass index; LBMI: lean body mass index; WC: waist circumference; WHR: waist–hip ratio; WHtR: waist–height ratio.
Figure 1. Comparison of adiposity measures: (A) Sex-stratified distribution of adiposity measures and (B) Spearman correlation in study population (n = 1923). Adiposity measurements were z-transformed (mean = 0, standard deviation = 1). AMMI: appendicular muscle mass index; BAI: body adiposity index; BF (%): body fat percentage; BMI: body mass index; FMI: fat mass index; LBMI: lean body mass index; WC: waist circumference; WHR: waist–hip ratio; WHtR: waist–height ratio.
Microorganisms 08 00547 g001
Microorganisms 08 00547 g002 550
Figure 2. Dependency of prediction capacity on prevalent and abundant OTUs in waist–height ratio, implemented with the combinations SVMReg NPK with CLR, and PLS4 with RA+Log: In each step, (A) OTUs present in less than x% of the samples were excluded (prevalence); this means for example in the step 80%, that OTUs with a prevalence of less than 80% were excluded for the calculations. (B) OTUs with an across-sample relative abundance of x% or less were excluded (abundance).
Figure 2. Dependency of prediction capacity on prevalent and abundant OTUs in waist–height ratio, implemented with the combinations SVMReg NPK with CLR, and PLS4 with RA+Log: In each step, (A) OTUs present in less than x% of the samples were excluded (prevalence); this means for example in the step 80%, that OTUs with a prevalence of less than 80% were excluded for the calculations. (B) OTUs with an across-sample relative abundance of x% or less were excluded (abundance).
Microorganisms 08 00547 g002
Table
Table 1. Benchmark data on the KORA (Cooperative Health Research in the Region of Augsburg) FF4 study population (mean (standard deviation)).
Table 1. Benchmark data on the KORA (Cooperative Health Research in the Region of Augsburg) FF4 study population (mean (standard deviation)).
CharacteristicsOverall (n = 1923)Men (n = 936)Women (n = 987)
Age (years)60.0 (12.1)60.3 (12.3)59.7 (11.9)
Body mass index (kg/m2)27.9 (5.0)28.3 (4.5)27.5 (5.4)
Waist circumference (cm)97.0 (14.2)102.8 (12.4)91.5 (13.6)
Waist–hip ratio0.91 (0.09)0.96 (0.07)0.85 (0.07)
Waist–height ratio 0.58 (0.08)0.59 (0.07)0.56 (0.09)
Body adiposity index31.0 (6.0)27.8 (4.0)34.0 (5.9)
Fat mass index (kg/m2)9.4 (3.4)8.1 (2.8)10.6 (3.5)
Lean body mass index (kg/m2)18.5 (2.6)20.1 (2.1)16.9 (2.1)
Appendicular muscle mass index (kg/m2)7.7 (1.3)8.6 (1.0)6.8 (1.0)
Body fat (%)32.9 (7.1)28.1 (5.3)37.5 (5.5)
Table
Table 2. Comparison of adiposity measures, transformation methods of operational taxonomic unit (OTU) count data, and machine learning algorithms in the study population (n = 1923): Each method includes 10-fold cross-validation, and for each logarithmic transformation, one pseudocount was added to the raw counts. CC: correlation coefficient; CLR: centered log-ratio transformation; ILR: isometric log-ratio transformation; NPK: normalized poly kernel; PLS: partial least squares regression; RF: random forest; RA: relative abundance; RMSE: root mean squared error; SVMReg: support vector machine regression.
Table 2. Comparison of adiposity measures, transformation methods of operational taxonomic unit (OTU) count data, and machine learning algorithms in the study population (n = 1923): Each method includes 10-fold cross-validation, and for each logarithmic transformation, one pseudocount was added to the raw counts. CC: correlation coefficient; CLR: centered log-ratio transformation; ILR: isometric log-ratio transformation; NPK: normalized poly kernel; PLS: partial least squares regression; RF: random forest; RA: relative abundance; RMSE: root mean squared error; SVMReg: support vector machine regression.
Data TransformationMachine Learning Algorithm
SVMReg NPKRFM5RulesPLS20PLS4
CCRMSECCRMSECCRMSECCRMSECCRMSE
Body mass index
Raw counts0.181.010.260.970.191.060.171.240.210.99
RA (100%)0.181.010.260.970.231.020.211.110.220.99
RA + Log0.300.970.250.970.261.010.221.220.320.96
Raw counts + CLR0.330.950.210.980.251.010.221.230.310.97
Raw counts + ILR0.330.950.190.980.280.990.221.230.310.97
Waist circumference
Raw counts0.240.980.320.950.221.070.221.190.250.98
RA (100%)0.240.980.330.950.231.080.271.090.280.97
RA + Log0.360.940.290.960.340.970.271.180.370.94
Raw counts + CLR0.390.930.280.960.340.970.271.190.370.94
Raw counts + ILR0.390.930.280.960.300.990.271.190.370.94
Waist-hip ratio
Raw counts0.230.990.310.960.181.060.211.160.230.98
RA (100%)0.230.990.290.960.251.020.241.110.260.97
RA + Log0.340.940.260.970.290.990.281.160.350.95
Raw counts + CLR0.360.940.250.970.281.000.281.170.350.95
Raw counts + ILR0.360.940.250.970.281.000.281.170.350.95
Waist-height ratio
Raw counts0.211.000.310.950.221.060.201.230.230.98
RA (100%)0.211.000.300.960.231.030.251.090.260.98
RA + Log0.330.950.270.960.300.990.241.210.360.94
Raw counts + CLR0.370.940.250.970.291.000.241.220.360.94
Raw counts + ILR0.370.940.260.970.281.000.241.220.360.95
Body adiposity index
Raw counts0.131.020.140.990.131.070.101.280.131.00
RA (100%)0.131.020.140.990.091.150.121.150.131.01
RA + Log0.230.990.140.990.181.030.131.290.230.99
Raw counts + CLR0.240.990.120.990.151.060.131.300.241.00
Raw counts + ILR0.240.990.120.990.131.060.131.300.241.00
Fat mass index
Raw counts0.141.020.170.990.111.070.121.210.161.00
RA (100%)0.141.020.180.980.161.050.141.140.161.01
RA + Log0.260.980.150.990.181.040.171.260.260.98
Raw counts + CLR0.280.970.170.990.191.030.171.270.260.99
Raw counts + ILR0.280.970.160.990.171.050.171.270.260.99
Lean body mass index
Raw counts0.250.990.320.950.181.120.191.260.240.98
RA (100%)0.250.990.330.950.251.020.251.100.270.97
RA + Log0.340.940.270.960.300.990.281.180.330.96
Raw counts + CLR0.360.940.260.970.281.000.281.190.330.96
Raw counts + ILR0.360.940.280.960.291.000.281.190.330.96
Appendicular muscle mass index
Raw counts0.250.990.320.950.221.070.201.240.230.98
RA (100%)0.250.990.290.960.211.060.241.100.270.97
RA + Log0.340.940.260.970.290.990.281.170.330.96
Raw counts + CLR0.350.940.290.960.291.000.281.180.330.96
Raw counts + ILR0.350.940.280.960.320.980.281.180.330.96
Body fat percentage
Raw counts0.141.020.190.980.081.120.101.190.151.00
RA (100%)0.141.020.180.980.041.240.121.160.151.00
RA + Log0.230.970.160.990.181.040.181.250.260.99
Raw counts + CLR0.280.970.140.990.151.050.181.250.251.00
Raw counts + ILR0.280.970.120.990.191.040.181.250.261.00
Table
Table 3. Comparison of predictive performance in study population and sex-stratified subsets: Each (sub-)set was scaled (adiposity measure) and transformed (OTUs) separately. CLR: centered log-ratio; NPK: normalized poly kernel; RA: relative abundance (100%); PLS4: partial least squares with four principal components; SVMReg: support vector machine regression.
Table 3. Comparison of predictive performance in study population and sex-stratified subsets: Each (sub-)set was scaled (adiposity measure) and transformed (OTUs) separately. CLR: centered log-ratio; NPK: normalized poly kernel; RA: relative abundance (100%); PLS4: partial least squares with four principal components; SVMReg: support vector machine regression.
Abdominal Adiposity MeasuresTotal Population (n = 1923)Men (n = 936)Women (n = 987)
CCRMSECCRMSECCRMSE
SVMReg NPK and CLR
Waist circumference0.390.930.240.980.320.96
Waist–hip ratio0.360.940.240.980.260.98
Waist–height ratio0.370.940.260.970.340.95
PLS4 and RA + Log
Waist circumference0.370.940.270.990.310.98
Waist–hip ratio0.350.950.241.010.271.00
Waist–height ratio0.360.940.280.990.340.97

Share and Cite

MDPI and ACS Style

Troll, M.; Brandmaier, S.; Reitmeier, S.; Adam, J.; Sharma, S.; Sommer, A.; Bind, M.-A.; Neuhaus, K.; Clavel, T.; Adamski, J.; et al. Investigation of Adiposity Measures and Operational Taxonomic unit (OTU) Data Transformation Procedures in Stool Samples from a German Cohort Study Using Machine Learning Algorithms. Microorganisms 2020, 8, 547. https://doi.org/10.3390/microorganisms8040547

AMA Style

Troll M, Brandmaier S, Reitmeier S, Adam J, Sharma S, Sommer A, Bind M-A, Neuhaus K, Clavel T, Adamski J, et al. Investigation of Adiposity Measures and Operational Taxonomic unit (OTU) Data Transformation Procedures in Stool Samples from a German Cohort Study Using Machine Learning Algorithms. Microorganisms. 2020; 8(4):547. https://doi.org/10.3390/microorganisms8040547

Chicago/Turabian Style

Troll, Martina, Stefan Brandmaier, Sandra Reitmeier, Jonathan Adam, Sapna Sharma, Alice Sommer, Marie-Abèle Bind, Klaus Neuhaus, Thomas Clavel, Jerzy Adamski, and et al. 2020. "Investigation of Adiposity Measures and Operational Taxonomic unit (OTU) Data Transformation Procedures in Stool Samples from a German Cohort Study Using Machine Learning Algorithms" Microorganisms 8, no. 4: 547. https://doi.org/10.3390/microorganisms8040547

APA Style

Troll, M., Brandmaier, S., Reitmeier, S., Adam, J., Sharma, S., Sommer, A., Bind, M. -A., Neuhaus, K., Clavel, T., Adamski, J., Haller, D., Peters, A., & Grallert, H. (2020). Investigation of Adiposity Measures and Operational Taxonomic unit (OTU) Data Transformation Procedures in Stool Samples from a German Cohort Study Using Machine Learning Algorithms. Microorganisms, 8(4), 547. https://doi.org/10.3390/microorganisms8040547

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop