Multi-omic signature of body weight change: results from a population-based cohort study
- Simone Wahl†1, 2, 3Email author,
- Susanne Vogt†1,
- Ferdinand Stückler4,
- Jan Krumsiek4,
- Jörg Bartel4,
- Tim Kacprowski5,
- Katharina Schramm6, 7,
- Maren Carstensen8, 9,
- Wolfgang Rathmann10,
- Michael Roden8, 9, 11,
- Carolin Jourdan1,
- Antti J Kangas12,
- Pasi Soininen12, 13,
- Mika Ala-Korpela12, 13, 14, 15,
- Ute Nöthlings16,
- Heiner Boeing17,
- Fabian J Theis4, 18,
- Christa Meisinger1,
- Melanie Waldenberger1, 2,
- Karsten Suhre19, 20,
- Georg Homuth5,
- Christian Gieger1, 2,
- Gabi Kastenmüller19,
- Thomas Illig1, 2, 21,
- Jakob Linseisen1,
- Annette Peters1, 3, 22,
- Holger Prokisch6, 7,
- Christian Herder8, 9,
- Barbara Thorand†1, 3 and
- Harald Grallert†1, 2, 3
© Wahl et al.; licensee BioMed Central. 2015
Received: 9 September 2014
Accepted: 20 January 2015
Published: 9 March 2015
Excess body weight is a major risk factor for cardiometabolic diseases. The complex molecular mechanisms of body weight change-induced metabolic perturbations are not fully understood. Specifically, in-depth molecular characterization of long-term body weight change in the general population is lacking. Here, we pursued a multi-omic approach to comprehensively study metabolic consequences of body weight change during a seven-year follow-up in a large prospective study.
We used data from the population-based Cooperative Health Research in the Region of Augsburg (KORA) S4/F4 cohort. At follow-up (F4), two-platform serum metabolomics and whole blood gene expression measurements were obtained for 1,631 and 689 participants, respectively. Using weighted correlation network analysis, omics data were clustered into modules of closely connected molecules, followed by the formation of a partial correlation network from the modules. Association of the omics modules with previous annual percentage weight change was then determined using linear models. In addition, we performed pathway enrichment analyses, stability analyses, and assessed the relation of the omics modules with clinical traits.
Four metabolite and two gene expression modules were significantly and stably associated with body weight change (P-values ranging from 1.9 × 10−4 to 1.2 × 10−24). The four metabolite modules covered major branches of metabolism, with VLDL, LDL and large HDL subclasses, triglycerides, branched-chain amino acids and markers of energy metabolism among the main representative molecules. One gene expression module suggests a role of weight change in red blood cell development. The other gene expression module largely overlaps with the lipid-leukocyte (LL) module previously reported to interact with serum metabolites, for which we identify additional co-expressed genes. The omics modules were interrelated and showed cross-sectional associations with clinical traits. Moreover, weight gain and weight loss showed largely opposing associations with the omics modules.
Long-term weight change in the general population globally associates with serum metabolite concentrations. An integrated metabolomics and transcriptomics approach improved the understanding of molecular mechanisms underlying the association of weight gain with changes in lipid and amino acid metabolism, insulin sensitivity, mitochondrial function as well as blood cell development and function.
KeywordsMetabolomics Transcriptomics Weight change Obesity Molecular epidemiology Bioinformatics
With an estimated 671 million obese individuals worldwide in 2013 , obesity has reached epidemic proportions. Considering the manifold health problems associated with excess body weight, including cardiovascular disease and type 2 diabetes, obesity poses a serious public health problem . Understanding the mechanisms by which excess body weight contributes to cardiometabolic risk is a prerequisite for advances in therapeutic approaches. Despite extensive research, however, the complex molecular basis of body weight-related metabolic perturbations is not fully understood.
Advances in the field of high-throughput omics technologies, including metabolomics and transcriptomics, offer the opportunity to simultaneously measure hundreds or thousands of molecules, for example, metabolites and gene transcripts, thereby allowing a deeper characterization of obesity-related pathomechanisms on a molecular level . In recent years, a number of cross-sectional efforts suggested a relationship between obesity and the human blood metabolome (for example, [4-6]) and transcriptome (for example, [7,8]), which extend to different tissues such as adipose tissue . In addition, weight loss upon behavioral intervention was associated with changes in the blood metabolome [5,10], suggesting that the observed obesity-related molecular signatures are at least in part reversible.
However, the effect of long-term body weight change on the human blood metabolome and transcriptome in the general population – rather than under clinical settings – is less well explored. Few prospective studies have investigated the association of body weight change with concentrations of a larger set of metabolites in healthy subjects and these are restricted to a panel of lipoprotein subclasses [11,12]. In addition, although multi-omic approaches have been fruitful in different applications to enhance the understanding of complex molecular pathways (for example, [13-15]), the potential of integrating multiple omics techniques has rarely been used in the study of weight change-associated metabolic effects in humans .
Here, we used data from Cooperative Health Research in the Region of Augsburg (KORA) S4/F4, which constitutes a large phenotypically and molecularly well-characterized population-based cohort. We aimed to characterize associations of body weight change over a seven-year follow-up period with serum metabolomics and whole blood transcriptomics data assessed at follow-up, to determine distinct groups of molecules associated with weight change using weighted correlation network analysis (WGCNA), to study the interrelation of these groups using partial correlation networks, to investigate the stability of the findings in relevant subgroups and upon additional multivariable adjustment, for example, of subjects with weight gain versus weight reduction, and to determine the relation of the identified omics signatures with clinical traits.
Written informed consent was obtained from all participants. The KORA studies have been approved by the ethics committee of the Bavarian Medical Association.
KORA (Cooperative Health Research in the Augsburg Region) is a research platform of independent population-based health surveys and subsequent follow-up examinations of community-dwelling adults living in the region of Augsburg in Southern Germany. Study design, sampling method and data collection have been described in detail elsewhere . The KORA S4 survey (1999 to 2001) comprised 4,261 participants, 25 to 74 years old . Of these, 3,080 subjects participated in the follow-up examination KORA F4 (2006 to 2008). The present study is based on a subsample of 1,658 participants of KORA S4/F4 with metabolomics data from two platforms available. Gene expression data were available for a subsample of 703 subjects, 62 to 77 years old in F4. Women who were pregnant at the time of the examinations or during the follow-up period were excluded prior to analysis.
Anthropometric measurements and interviews
In both examinations, body weight, height, waist and hip circumference as well as systolic and diastolic blood pressure were measured using standard protocols as described elsewhere . Information on lifestyle factors and comorbidities was collected in a structured interview by trained interviewers. Intake of medication within seven days prior to examination was recorded with the IDOM-Software . We combined information obtained in KORA S4 and F4 to determine changes in the variables during the follow-up period. Change in lifestyle factors, comorbidities and medication was defined as described in Additional file 1.
At the follow-up examination, blood samples were collected during study center visits between 8 a.m. and 11 a.m., after subjects were instructed to fast overnight for at least eight hours. Whole blood was collected using PAXgene Blood RNA tubes (BD, Heidelberg, Germany) and stored at −80°C until analysis. Red blood cell (RBC) count, hematocrit, mean corpuscular hemoglobin (MCH), mean corpuscular haemoglobin concentration (MCHC) as well as mean cell volume of erythrocytes (MCV) were measured with the impedance and the cyanmethemoglobin method using the LH 750 Hematology Analyzer (Beckman Coulter, Brea, CA, USA). Glycated hemoglobin (HbA1c) was determined using the high performance liquid chromatography method (HA 8160, Menarini, Florence, Italy). For serum collection, blood was drawn into serum gel S-Monovette tubes (Sarstedt, Nümbrecht, Germany), gently inverted two to three times and rested for 30 minutes at room temperature to obtain complete coagulation. The material was then centrifuged for 10 minutes (2,750 g at 15°C). Serum was aliquoted into synthetic straws which were kept for a maximum of six hours at 4°C before storage at −80°C until analysis. Fasting glucose levels (FGlc) and glucose levels two hours post challenge during an oral glucose tolerance test (2hGlc) were assessed with the hexokinase method (GLU Flex), low (LDL) and high density lipoprotein (HDL) cholesterol levels with the CHOD-PAP method (LDL: ALDL Flex, HDL: AHDL Flex), and fasting triglyceride (TG) levels with the GPO-PAP method (TGL Flex; all assays from Dade Behring, Eschborn, Germany). C-reactive protein was measured with nephelometry on a BN II using reagents from Siemens (Eschborn, Germany).
Metabolite detection and quantification were performed on two different platforms, a commercial mass spectrometry (MS)-based platform at the company Metabolon Inc. (Durham, NC, USA) (N = 1,768) [20,21] and a nuclear magnetic resonance (NMR) spectroscopy-based platform (N = 1,788) .
The analytical platform developed by Metabolon is based on two ultrahigh-performance liquid chromatography/tandem mass spectrometry (UHPLC/MS/MS2) injections and one gas chromatography/mass spectrometry (GC/MS) injection per sample. The two UHPLC injections were optimized for basic and acidic species, respectively. More detail is given in  and . The platform provides relative quantification for a total of 517 compounds, 325 of which could be identified based on a standard library of MS/MS spectra. The identified small molecules cover a large number of metabolite classes, including fatty acids, ketone bodies, glycerophospholipids, sphingolipids, acylcarnitines, bile acid metabolites, amino acids, peptides, carbohydrates, xenobiotics, vitamins and nucleotide metabolites. The full list of metabolites with pathway annotations is given in Table S1 in Additional file 2.
Technical effects were controlled for by dividing the metabolite concentration values by the median value of samples measured on the same day for each metabolite. See Figure S1 in Additional file 3 for the distribution of the variability (expressed as coefficient of variation (CV)) of the day medians across the metabolites. Metabolite medians varied with a median CV of 0.494 (range 0.307 to 1.220). In addition, outlier values of >4 standard deviations from the mean of the respective metabolite on the log10 scale were set to missing. Finally, 81 metabolites (42 identified, 39 unidentified) with more than 50% missing values, and 5 observations with more than 20% missings were excluded, leaving a total of 1,763 observations of 434 metabolites (281 identified, 153 unidentified) for analysis.
NMR spectroscopy platform
As another technique to obtain quantitative information on metabolic compounds, an NMR spectroscopy platform was used. The precise experimental methodology has been described elsewhere [13,22]. Computational strategies of metabolite identification and quantification from the NMR spectra are described by Inouye et al. . A total of 130 metabolite concentrations and derived measures were obtained from the NMR spectroscopy platform (see Table S1 in Additional file 2 for a list of metabolites including full names and pathway annotations).
Preprocessing of NMR data was similar to Metabolon data in terms of outlier exclusion and detection rate thresholds. None of the metabolite traits had more than 50% missing values; however, four observations with more than 20% missings were excluded from the data set leaving a total of 1,784 observations for analysis.
Combined analysis of multiple metabolomics platforms
For 1,658 subjects, both Metabolon and NMR data were available (see Additional file 1: Table S2 in Additional file 2 and Figure S2 in Additional file 3 for information on the concordance of measurements for metabolites measured on both platforms). [M] and [N] is added to the metabolite names in this work to indicate measurement on the Metabolon and the NMR spectroscopy platform, respectively. Metabolites of both platforms were assigned to super-pathways and sub-pathways in accordance with the pathways proposed by Metabolon on the basis of Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Thereby, super-pathways largely represent chemical classes analogous to the top-level pathway definitions in KEGG, whereas sub-pathways try to classify the metabolites according to their role in metabolism. As an example, the metabolites acetoacetate and 3-hydroxybutyrate (BHBA) were assigned to the super-pathway ‘Lipids’ and to the sub-pathway ‘Ketone bodies’, according to their role in ketogenesis.
Gene expression measurement
RNA isolation and gene expression measurement have been described in detail elsewhere [23,24]. Gene expression data were exported to the statistical environment R, version 2.14.2 . Data were quantile normalized using the R package lumi, version 2.8.0, from the Bioconductor platform .
Multiple imputation of missing values
Data from both metabolomics platforms contained a large number of missing values. These were treated with multiple imputation by chained equations (MICE) using the R package mice, version 2.21 . Few missing values in phenotypic variables were imputed in the same imputation process. Twenty imputed data sets were generated with 10 iterations each. See Additional file 1: Table S3 in Additional file 2, and Figures S3 to S7 in Additional file 3 for descriptives of missing values in the data set, for considerations on missing value handling, for the precise imputation settings and for imputation diagnostics.
Weighted correlation network analysis (WGCNA)
Previous studies have shown that clusters of related genes may be more reproducibly associated with a phenotype or disease than single genes  and that testing groups of metabolites instead of single metabolites improved power in a genome-wide association study . Furthermore, high correlations were observed within groups of molecules in the metabolomics and transcriptomics data sets in this study. Thus, our strategy was – in addition to providing single metabolite/transcript associations – to cluster metabolomics and transcriptomics data prior to testing the association with body weight change, in order to obtain joint association signals of groups of correlated metabolites/transcripts. Specifically, metabolite and transcript levels determined at follow-up were clustered by means of WGCNA using the R package WGCNA, version 1.34 [30,31].
Transcripts were pre-selected based on their association with metabolites to reduce the number of transcripts to those related to the blood metabolome, aiming to improve the power and specificity to observe ΔBW-related transcripts that are relevant with regard to metabolic processes. Prior to analysis, gene expression data were log2 transformed and adjusted for microarray-specific technical variables (RNA integrity number, amplification plate indicators as well as sample storage time)  by determining the residuals from linear models of gene expression on these variables. Then, association with metabolites was determined using linear models with transformed metabolite as response and adjusted transcript as covariate, adjusting for age and sex, and a linear model additionally adjusted for body weight and ΔBW to avoid the selection of transcripts related to metabolites due to their common association with these variables. A total of 2,537 transcripts with at least a suggestive association (P <10−5) in both models were selected.
WGCNA was jointly applied to the 411 serum metabolites and derived measures (281 identified metabolites from the Metabolon platform and 130 metabolites and derived measures from the NMR spectroscopy platform, see Table S1 in Additional file 2; N = 1,631) and separately to the 2,537 metabolite-related transcripts (N = 689). Briefly, the topological overlap matrix (TOM) was derived from the Pearson’s correlations among the variables using the signed topological overlap dissimilarity measure . Thereby, a soft-thresholding power of 13 (metabolite network) and 8 (transcript network) was selected based on the scale-free topology criterion (Figure S8 in Additional file 3). The TOM was then subjected to hierarchical clustering with distance between clusters defined through average linkage, and distinct clusters (‘modules’) were derived using a dynamic tree cutting algorithm with a minimum module size of five molecules , followed by merging closely correlated modules at a dendrogram height of 0.25 (see Figure S9 in Additional file 3 for the resulting cluster dendrograms). For each module, the module eigengene (ME), a module-representative variable, also interpretable as the ‘center’ of the module, was derived as the first principal component of a principal component analysis on the scaled matrix of molecules corresponding to the respective module. Strength of module membership, that is, the contribution of the molecules to the module, was computed as the correlation of each molecule with the respective ME.
Modeling association with body weight change
Linear models were used to determine the association of body weight change with the metabolite and gene expression modules. Thereby, body weight change was defined as annual percentage body weight change: ΔBW = 100% * ((weight (F4) - weight (S4))/weight (S4))/follow-up years, where weight gain was coded as positive weight change, and weight loss as negative weight change. We decided to study ΔBW rather than change in body mass index (ΔBMI), since in adults, ΔBW might be considered a more appropriate measure of weight change. When looking at change over time in adults, the division by squared height (as BMI = BW/(body height)2) does not improve accuracy but might rather add noise due to measurement error of height measurement, and causes an unwanted increase in BMI due to height shrinkage during aging. However, ΔBW and ΔBMI correlate highly (correlation coefficient 0.98), so that results should not be very different.
Each ME was modeled as response variable, and ΔBW, sex, age and body weight at baseline were modeled as covariates. The effects of ΔBW on MEs were tested using Wald tests, and P-values were corrected for multiple testing using the Bonferroni method (at P <0.05/(number of modules)). Association of single metabolites and metabolite-related transcripts with ΔBW were assessed in a similar way, using significance levels of P <0.05/411 and P <0.05/2,537, respectively.
The chosen modeling strategy is restricted to finding linear associations between ΔBW and MEs/molecules across the complete ΔBW range, comprising both subjects with weight loss (negative ΔBW) and weight gain (positive ΔBW). Assuming that weight loss and weight gain might not show strictly opposing metabolic effects, the analysis was repeated stratified to the groups with weight loss and weight gain. Formally, the above mentioned model was extended by subgroup index as covariate and by interaction terms of the subgroup index with ΔBW and the other covariates. Similarly, subgroup analyses were performed in obese (N = 426; defined as BMI >30) and non-obese (N = 1,205) subjects, in central obese (N = 522; defined as waist-hip ratio >1 in men and waist-hip ratio >0.85 in women) and not central obese (N = 1,109) subjects, in men (N = 828) and women (N = 803), as well as in subjects >55 years old (N = 855) and >55 years old (N = 776) at baseline.
Furthermore, as an explorative approach towards a biological explanation of the observed associations, we studied their sensitivity through additional adjustment for three groups of variables. The first model was adjusted for changes in lifestyle factors during follow-up, including change in physical activity, smoking, alcohol consumption and sleeping behavior as well as nutrition habits at baseline. The second model was adjusted for incident comorbidities during follow-up, including diabetes, cancer, myocardial infarction and stroke. The third model was adjusted for changes in medication, including beta blocker, metformin, other anti-diabetic medication, systemic corticosteroids, oral contraceptives and antidepressants. [See Additional file 1 for definitions of these variables].
Finally, cross-sectional associations of modules with binary and log-transformed continuous clinical traits in KORA F4 (metabolic syndrome, myocardial infarction (MI), stroke, HDL cholesterol, LDL cholesterol, FGlc, 2hGlc, HbA1c, systolic and diastolic blood pressure, C-reactive protein) were investigated by means of logistic and linear regression models, respectively, adjusted for age, sex, body weight, lipid-lowering medication (overall, statin, fibrates), antihypertensive medication, antidiabetic medication as well as systemic corticoid intake.
Pathway enrichment analyses
To formally investigate whether the identified metabolite modules were enriched for specific biological pathways (super- and sub-pathways as described above, Table S1 in Additional file 2), weighted enrichment analyses were performed as applied before in different contexts [33,34]. Briefly, for each pathway c and module m, the enrichment statistics Scm was defined as the sum of module membership measures across all metabolites assigned to the respective pathway, whereas metabolites from other modules were assigned zero weight. Pathway assignment of all metabolites was randomly permuted 100,000 times, and enrichment statistics Scm(perm) were computed. Permutation P-values were then defined as the number of enrichment statistics Scm(perm) larger than the original Scm.
For the gene expression modules, we explored enrichment for gene ontology (GO) terms using the R packages GO.db, version 2.9.0, AnnotationDbi, version 1.22.6, and org.Hs.eg.db, version 2.9.0. Furthermore, the commercial software Ingenuity Pathway Analysis (IPA) was applied to identify enriched canonical pathways (IPA build version 312825 M, content version 18841524, release date: 24 June 2014; analysis date: 4 July 2014 ). The reference set was restricted to genes represented on the IlluminaHT-12 v3 BeadChip, and only human annotations were considered. In case multiple probes mapped to one gene, the probe exhibiting the largest module membership was considered for downstream analyses. Pathway analyses were performed with IPA’s Core Analysis module. Both GO and IPA enrichment analyses are based on Fisher’s exact test.
Construction of a multi-omics network
A partial correlation network was constructed from the modules significantly associated with ΔBW as described in detail . For each pair of modules, the partial correlation coefficient of the respective MEs was calculated as the Pearson’s correlation coefficient of the residuals of the two MEs with regard to all other MEs, as well as sex, age, body weight and ΔBW. Since strong interrelationships among the modules might result in spurious negative partial correlations, pairwise marginal correlation (that is, the Pearson’s correlation, uncorrected for any other variables) was taken as a prerequisite.
Effects on inter- and intra-module connectivity within the multi-omics network were studied as follows: inter- and intra-module connectivity was defined as the correlation of the MEs between two modules, and the average module membership strength (definition see above) across the module members of a certain module, respectively. Significance of the difference in inter-/intra-module connectivity between the groups of weight gain and weight loss was determined through permutation testing (where weight change status was randomly shuffled).
All statistical analyses were performed in R, version 3.0.1 .
Results and discussion
Characteristics of the study population
Metabolomics data ( N = 1,631)
Combined metabolomics and transcriptomics data ( N = 689)
Body weight (kg), baseline
Body weight (kg), follow-up
Δ Body weight (%)
Δ Body weight/year (%)
BMI (kg/m2), baseline
BMI (kg/m2), follow-up
Age (years), baseline
Age (years), follow-up
Relative frequency (%)
Relative frequency (%)
Weight change direction (reduction/gain)
Body weight change is globally associated with the metabolite profile
Four of the eight MetMs were significantly associated with annual percentage body weight change (ΔBW), in linear models adjusted for age, sex and baseline body weight (positive associations for MetM1, P = 1.2 × 10−24, MetM3, P = 2.2 × 10−4, and MetM4, P = 7.3 × 10−17; negative association for MetM5, P = 1.7 × 10−14; all significant after Bonferroni correction for 27 modules). These four modules comprised a total of 147 metabolites. Together, these metabolites covered major branches of metabolism captured by the metabolomics platforms, including lipid metabolism, amino acids and peptides, carbohydrate metabolism, cofactors and vitamins, and energy metabolism (Figure S10 in Additional file 3). This suggests a global association of body weight change with the serum metabolome.
These results demonstrate that ΔBW strongly associates with lipoprotein constituents (see also Figure S11 in Additional file 3), amino acids and peptides, as well as metabolites of energy metabolism, and that clustering helped to reveal pathways jointly and strongly associated with ΔBW.
The metabolic signature associated with body weight change is consistent with known pathophysiology of obesity
Overall, the metabolite signature associated with ΔBW concurred with known aspects of the pathophysiology of weight change and obesity. The associations of ΔBW with lipoprotein subclasses (positive association with VLDL, LDL and S-HDL subclasses, negative association with larger HDL particles and with HDL and LDL particle size; Figure 3 and Figure S11 in Additional file 3) are in agreement with the observations of two smaller prospective studies that analyzed the effect of weight change over similar time periods (9 and 6.5 years, respectively) on lipoprotein subclasses [11,12]. Specifically, ΔBW was positively associated with increases in VLDL and LDL subclasses, and with decreases in L-HDL, whereas S-HDL behaved oppositely . ΔBW was also negatively related to LDL and HDL particle sizes [11,12]. The clustering of S_HDL_TG [N] (TG in S-HDL) within the TG/VLDL module in our study was also in agreement with its close correlation with VLDL subclasses in Inouye et al. , where S-HDL behaved differently from larger HDL subclasses with regard to metabolite-transcript associations.
Mechanisms by which body weight increase gives rise to the described changes may include an increased release of free fatty acids from adipose tissue, triggering hepatic TG and VLDL production  and increasing the activity of hepatic lipase . Hepatic lipase is involved in the exchange of TGs from VLDL against cholesterol esters from HDL, thereby promoting the production of small dense LDL. Together with phospholipid transfer protein (PLTP) and cholesterol ester transfer protein (CETP), which also show increased levels upon obesity , hepatic lipase is centrally involved in regulating HDL particle size. Interestingly, this was reflected in oppositional associations of genetic variants in the respective genes LIPC, PLTP and CETP with small versus large HDL subclasses .
ΔBW was strongly associated with amino acid concentrations, most prominently BCAAs, phenylalanine, tyrosine and glutamate. The increase of these amino acids in obesity has long been known , and has also been observed in more recent studies (for example, ). The underlying mechanism might be an impaired catabolism of BCAA upon obesity . Experimental studies show that BCAAs inhibit the insulin receptor substrate via the mTOR/p70S6K/S6K pathway . Accordingly, in the study by Newgard et al. , addition of BCAAs to a high-fat diet in rats promoted the development of insulin resistance. Recently, BCAAs, phenylalanine and tyrosine were shown to associate with future insulin resistance , future type 2 diabetes  and prevalent metabolic syndrome . In our study, the BCAA module associated positively with markers of insulin resistance (HbA1c: P = 7.2 × 10−10; 2hGlc: P = 9.0 × 10−9) and metabolic syndrome prevalence (P = 3.8 × 10−5) (Figure 4).
The metabolic signature associated with body weight change points towards mitochondrial dysfunction
The positive association of ΔBW with several metabolites in the TG/VLDL and BCAA modules is supportive of a link between body weight gain and mitochondrial dysfunction. Mitochondria are important organelles in the regulation of metabolism . As a main function, they produce adenosine triphosphate (ATP) from carbohydrates, fats and proteins via the tricarboxylic acid cycle, which is supplied with pyruvate (from glycolysis), acetyl-CoA (from β-oxidation) and amino acid metabolites (from protein catabolism) [49,50]. In states of insufficient oxygen supply or mitochondrial dysfunction, pyruvate is not imported into the mitochondria but instead converted to lactate via lactic acid fermentation and to alanine via transamination. Obesity is associated with decreased fatty acid β-oxidation so that obese individuals are more dependent on the glycolytic pathway for ATP production , resulting in an increased pyruvate production. At the same time, obesity is associated with diminished mitochondrial biogenesis , reduced mitochondrial size and diminished respiratory chain activity . These effects might explain increased circulating levels of pyruvate, lactate and alanine upon weight gain in our study. Concentrations of these metabolites have previously been shown to be elevated in obesity , and lactate levels increased upon weight gain during chemotherapy in early breast cancer patients . Furthermore, concentrations of pyruvate, lactate and alanine were found to be predictive of future glucose intolerance, independent of body mass .
All three metabolites clustered in the TG/VLDL module, together with VLDL subclasses, fatty acids, propionylcarnitine and phenylalanine. The blood concentrations of fatty acids and circulating acylcarnitines (among them propionylcarnitine) were found to be increased in obese subjects . Elevated levels of these metabolites are also thought to be linked to mitochondrial dysfunction . In addition, aluminum-induced mitochondrial dysfunction was found to promote VLDL secretion in human hepatocytes , suggesting a role of mitochondrial dysfunction as a further link between ΔBW and dyslipidemia as well as cardiometabolic disease.
Furthermore, two C5 acylcarnitines, isovalerylcarnitine and 2 − methylbutyroylcarnitine, BCAAs, tyrosine, tryptophan and creatinine clustered in the BCAA module. Short-chain acylcarnitines (C4 and C5) are products of BCAA catabolism. As catabolism of BCAAs is performed in the mitochondrial matrix, elevated BCAA levels could be an indicator for mitochondrial dysfunction . A recent retrospective study found primary mitochondrial respiratory chain disease to be significantly associated with elevated BCAA levels . Moreover, respiratory chain inhibition in cultured muscle cells was found to be significantly associated with a reduced uptake of several other amino acids, among them phenylalanine, tyrosine and tryptophan (together with a greater secretion of lactate, alanine and creatine) . Together, metabolites of the TG/VLDL and the BCAA module provide a link between mitochondrial dysfunction and body weight gain already in a non-obese state.
Individual metabolite associations with body weight change provide further interesting insights
As a result of the minimum module size of five molecules chosen in WGCNA, metabolites reflecting biological pathways that are represented by fewer than five correlated metabolites on the metabolomics platforms are less likely to cluster in modules sharing association with ΔBW. Since these might, nevertheless, be interesting, we provide the overall single metabolite association results in Table S4 in Additional file 2. These include the positive association of ΔBW with the tryptophan metabolites hydroxytryptophan [M] and kynurenine [M], which are successors of tryptophan in the serotonin and niacin biosynthesis pathways, respectively, and negative association with serotonin (5HT) [M]. They also include the positive association with the xenobiotics caffeine [M] and piperine [M]. Increased concentrations in these substances may result either from an increased consumption of coffee or caffeinated drinks, and herbs or spices, or from a lower degradation or excretion. Furthermore, we observed a negative association of ΔBW with quinate [M] and catechol sulfate [M], a positive association with bradykinin, des-arg(9) [M], the active metabolite of the vasodilating peptide hormone bradykinin, and a positive association with the metabolites N1-methyl-3-pyridon-4-carboxamide [M] and N1-methyladenosine [M] from the nucleotide super-pathway.
Body weight change associates with a gene expression module related to erythrocyte development and a lipid-leukocyte-like module
Two of the 19 gene expression modules (GenMs) were significantly associated with ΔBW in linear models adjusted for age, sex and baseline body weight (positive association for GenM6, P = 3.8 × 10−12; negative association for GenM14, P = 1.9 × 10−4).
Consequently, we hypothesized that the transcripts in the RBC module are reflective of erythrocyte development, since immature RBCs, reticulocytes, contain remnant mRNA  which is depleted during erythrocyte maturation. An increased hematopoiesis upon diet-induced obesity in rats has been observed, putatively through action of leptin in the bone marrow , whereas glycosylated hemoglobin shows an inverse relationship with erythrocyte survival . A shift towards a larger proportion of immature RBCs upon weight gain would be consistent with these observations. In agreement with our findings, a small transcriptomics study of obesity reported the majority of the genes significantly upregulated in obesity to be related to reticulocytes .
The proportion of immature reticulocytes was shown to be elevated in cardiac disease patients . In a transcriptomics study of participants from the Framingham heart study, a cluster of transcripts specific to CD71+ early erythroid cells was significantly upregulated in coronary heart disease (CHD) cases compared to controls . This cluster, comprising 126 transcripts, showed a large overlap with the RBC module identified in our study, with nine of the above-mentioned core transcripts being contained in this cluster. Furthermore, Zhang et al. observed an increased RBC count and increased hemoglobin levels in obese but not in non-obese CHD patients as compared to healthy controls . Also, the obese patients were more likely to have acute coronary syndrome, which the authors attribute to a potential role of RBCs in the development of plaque instability. Our results suggest that even in non-obese subjects, weight change is related to RBC development.
Another, smaller, gene expression module (GenM14, comprising 17 transcripts) was negatively associated with ΔBW (Figure 5, Table S6 in Additional file 2). GenM14 contained all 11 transcripts of the ‘lipid-leukocyte (LL) module’ previously described as a leukocyte gene expression module strongly related to blood lipids  and a large number of serum metabolites including lipoprotein subclasses, lipids, glycoproteins and amino acids . We therefore termed it the ‘LL-like module.’ The core of the LL-like module (module membership strength >0.8) comprised the LL genes HDC, GATA2, SLC45A3, MS4A2 and SPRYD5. Inouye et al.  discussed this module as being involved in basophil and mast cell-related immune response and allergy. For instance, the core gene HDC codes for a protein converting histidine to histamine, which is secreted by basophils and mast cells in response to IgE sensitization. Accordingly, when we applied Ingenuity pathway analysis to our data, ‘FC Epsilon RI Signaling’ (P = 1.1 × 10−4), ‘Histamine Biosynthesis’ (P = 9.1 × 10−4) and ‘Airway Inflammation in Asthma’ (P = 3.7 × 10−3) were the top three canonical pathways (Table S8 in Additional file 2). Six further transcripts, IL4, TRIM49L1, TEX101, EPAS1/HIF2-α, CCNA1 and CAV2, were co-expressed with the LL module genes, although being less strongly correlated with the module center (module membership strengths ranging from 0.48 to 0.54), suggesting that they might share functionality with the LL module genes. Indeed, IL4 codes for the cytokine interleukin 4 which has long been known to induce differentiation of naïve T cells to Th2 cells that play a role in allergen response and which is secreted by basophils as a reaction to allergens . EPAS1/HIF2-α encodes a component of the hypoxia inducible transcription factor (HIF), which regulates responses to reduced oxygen and for which also a role in regulating inflammation  and energy balance  has been reported. Of note, the association of the LL-like module with ΔBW as well as with the TG/VLDL module was negative. Although these findings are in line with the negative association between the LL module and VLDL metabolites reported by Inouye et al. , they are contradictory to an analysis by Gonen et al., in which VLDL was found to trigger the release of histamine from human basophils . Furthermore, obesity is a risk factor for asthma and weight gain was found to increase the risk of developing airway hyperresponsiveness . It remains to be determined how these results fit with our observation of decreased expression of genes related to basophil/mast cell level or function being associated with weight gain.
Neither of the ΔBW-related GenM’s seemed to comprise genes with a well-established relationship to lipid metabolism, as might be expected after preselecting metabolite-related transcripts. Exemplarily, we looked up the genes LIPC, CETP and PLTP discussed above within the context of lipoprotein metabolism, as well as ABCG1 which has been discussed together with CETP as a strongly upregulated transcript in adipose tissue in response to diet-induced weight loss . Whereas ABCG1 transcripts tended to show a negative association with ΔBW (best P = 6.7 × 10−5 for transcript ILMN_2329927, which clustered in GenM1), transcripts of the other three genes were not related to either ΔBW or metabolites. These results could have been expected considering the tissue origin of these proteins. In line with these findings, the strong obesity-related changes in adipose tissue gene expression were weakly represented by blood cell transcriptomics in the study by Emilsson et al. .
Interrelation of the omics modules and sensitivity analyses
To investigate the interrelatedness of the identified ΔBW-related metabolite and gene expression modules conditional on all other modules and the above-mentioned covariates, we constructed a partial correlation network from the MEs (Figure 1; a larger illustration is provided in Figure S12 in Additional file 3). The six ΔBW-related modules were interrelated. The strongest positive partial correlation was observed between the TG/VLDL module and the LDL/IDL module (P = 2.3 × 10−54). In addition, the TG/VLDL module showed a strong negative correlation with the HDL module (P = 1.8 × 10−72) and with the LL-like module (P = 6.8 × 10−29).
In addition, the effect of ΔBW on inter- and intra-module connectivity of the network elements was investigated, since previous studies suggested sensitivity of metabolic network topology towards external factors [13,69] (see Methods). We did not observe any significant differences in network connectivity between the groups with weight gain and weight loss (all P >0.01).
The generally opposing associations of weight loss and weight gain with blood metabolite concentrations are in line with the studies by Mäntyselkä et al.  and Naganuma et al. , where the majority of ΔBW associations with lipoprotein measures were linear across the weight change range, and weight loss and weight gain showed opposing effects. Interestingly, the effect of weight loss (≥5% across 6.5 years) versus stable weight on VLDL subclasses and L-HDL was stronger in absolute terms than the effect of weight gain (≥5%) versus stable weight . These findings are in accordance with the stronger associations for the TG/VLDL module and the HDL module observed for weight loss in our study. Although larger studies in subjects with a larger range of weight change might have more power to differentially investigate the effects of weight loss versus weight gain, our results suggest that differences are not large and that, in general, weight loss is capable of reversing the effects of weight gain on the blood metabolome and transcriptome. Accordingly, it was shown in a randomized controlled trial that normalization of obesity led to a reversal of an unfavorable LDL subclass pattern .
Second, further subgroup analyses were performed, assuming that the weight change effect might depend on (central) obesity, on sex and on age (Figure 6, Figure S13 in Additional file 3). Again, no significant subgroup-specific effects were observed, although, as already mentioned above, this study might not provide sufficient power to study effect modification.
Body weight change over a period of seven years might be due to several reasons, including changes in lifestyle, the occurrence of diseases and changes in medication. For these reasons, we investigated the sensitivity of the observed associations with ΔBW towards adjustment for changes in lifestyle factors, for disease incidence and, finally, for changes in medication in three separate models (Table S9 in Additional file 2, Figure S14 in Additional file 3). None of the three models showed a change in effect sizes across the modules, indicating that the observed associations were primarily due to the change in body weight per se rather than the mechanisms that might have facilitated weight change. Note, however, that the majority of the variables reflecting changes in lifestyle, disease and medication were obtained from interviews and might, therefore, have insufficient accuracy. Also, nutrition was only obtained from the baseline time point so that the effect of changes could not be investigated.
Together, these results suggest that the metabolite-gene network and its relation to weight change reflect a largely stable system.
Several extensions of our study seem worthwhile. First, it would be interesting to obtain a higher resolution of body weight measurements during follow-up, as well as of metabolomics and gene expression measurements, to decipher the longitudinal sequence of metabolic changes and to study the metabolic processes related to weight cycling. An important limitation of this study is that metabolomics and transcriptomics data were not available from the baseline time point from all platforms and for all subjects. The present observational study does not allow conclusions on the effect directions underlying the observed associations of body weight change with metabolite or transcript levels. In addition, extending whole blood transcriptomics to different tissues seems extremely promising, considering that blood might only weakly reflect weight-related transcriptional changes in tissues  and that blood metabolites originate from different tissues. However, gene expression signatures identified in blood will be of large practical relevance since blood is most easily accessible also in a clinical setting. Also, in the context of weight change and its metabolic consequences, integrating metabolomics and blood cell transcriptomics is relevant from the perspective that blood cells may interact with blood substances in the etiology of atherosclerotic events .
Through the integration of two-platform serum metabolomic and whole blood transcriptomic data and the formation of modules of closely connected molecules, we obtained a comprehensive characterization of the metabolomic and transcriptomic signature of body weight change over a seven-year period in a large population-based cohort. Weight gain and weight loss were strongly and opposingly associated with the blood metabolome, with VLDL, LDL and large HDL subclasses, TGs, BCAAs and markers of energy metabolism as core molecules of the four metabolite modules. These associations point towards the development of dyslipidemia, disturbed amino acid metabolism as well as mitochondrial dysfunction upon weight gain. Two weight change-related gene expression modules pinpoint reticulocytes and immune cells (mast cells, basophils) as blood cell types putatively playing a role in body weight-related blood metabolism. Metabolite and gene expression modules were associated with clinical phenotypes, suggesting a role in linking excess body weight with metabolic and cardiovascular comorbidities. Our findings also support the hypothesis that clustering omics data prior to analyzing associations with a phenotype has increased power to identify biologically relevant pathways [28,29]. The LL-like module was found to be associated with weight change, although none of the contributing genes showed a univariate association with weight change that would have passed significance after correction for multiple comparisons.
Together, our study provides evidence for a largely reversible effect of long-term body weight gain in the general population on an integrated blood metabolomic and transcriptomic network. This improves the knowledge on molecular processes elicited by weight change and potentially linking it to comorbidities.
This work has been supported by the “Kompetenznetz Adipositas (Competence Network Obesity)” funded by the Federal Ministry of Education and Research (FKZ: 01GI1121B) and the German Center for Diabetes Research (DZD e.V.). FS was supported by a PhD student fellowships from the “Studienstiftung des Deutschen Volkes”. FJT and JB are supported by an ERC starting grant (LatentCauses). JK is supported by a grant from the German Helmholtz Postdoctoral Programme. PS has been supported by the Academy of Finland. MAK has been supported by the Sigrid Juselius Foundation and the Strategic Research Funding from the University of Oulu.
The KORA research platform (KORA, Cooperative Research in the Region of Augsburg) 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. We are indebted to all study participants. Furthermore, we are grateful to the field staff in Augsburg and Munich who were involved in the conduct of the KORA studies, and the staff of the Genome Analysis Center at the Helmholtz Zentrum München involved in the omics measurements.
- Ng M, Fleming T, Robinson M, Thomson B, Graetz N, Margono C, et al. Global, regional, and national prevalence of overweight and obesity in children and adults during 1980–2013: a systematic analysis for the Global Burden of Disease Study 2013. Lancet. 2014;384:766–81.View ArticlePubMedPubMed CentralGoogle Scholar
- Guh DP, Zhang W, Bansback N, Amarsi Z, Birmingham CL, Anis AH. The incidence of co-morbidities related to obesity and overweight: a systematic review and meta-analysis. BMC Public Health. 2009;9:88.View ArticlePubMedPubMed CentralGoogle Scholar
- Cornelis MC, Hu FB. Systems epidemiology: a new direction in nutrition and metabolic disease research. Curr Nutr Rep. 2013;2:225–35.View ArticleGoogle Scholar
- Newgard CB, An J, Bain JR, Muehlbauer MJ, Stevens RD, Lien LF, et al. A branched-chain amino acid-related metabolic signature that differentiates obese and lean humans and contributes to insulin resistance. Cell Metab. 2009;9:311–26.View ArticlePubMedPubMed CentralGoogle Scholar
- Oberbach A, Bluher M, Wirth H, Till H, Kovacs P, Kullnick Y, et al. Combined proteomic and metabolomic profiling of serum reveals association of the complement system with obesity and identifies novel markers of body fat mass changes. J Proteome Res. 2011;10:4769–88.View ArticlePubMedGoogle Scholar
- Wahl S, Yu Z, Kleber M, Singmann P, Holzapfel C, He Y, et al. Childhood obesity is associated with changes in the serum metabolite profile. Obes Facts. 2012;5:660–70.View ArticlePubMedGoogle Scholar
- Emilsson V, Thorleifsson G, Zhang B, Leonardson AS, Zink F, Zhu J, et al. Genetics of gene expression and its effect on disease. Nature. 2008;452:423–8.View ArticlePubMedGoogle Scholar
- Ghosh S, Dent R, Harper ME, Gorman SA, Stuart JS, McPherson R. Gene expression profiling in whole blood identifies distinct biological pathways associated with obesity. BMC Med Genomics. 2010;3:56.View ArticlePubMedPubMed CentralGoogle Scholar
- Pietilainen KH, Naukkarinen J, Rissanen A, Saharinen J, Ellonen P, Keranen H, et al. Global transcript profiles of fat in monozygotic twins discordant for BMI: pathways behind acquired obesity. PLoS Med. 2008;5:e51.View ArticlePubMedPubMed CentralGoogle Scholar
- Reinehr T, Wolters B, Knop C, Lass N, Hellmuth C, Harder U, et al. Changes in the serum metabolite profile in obese children with weight loss. Eur J Nutr. In press.
- Naganuma R, Sakurai M, Miura K, Yoshita K, Morikawa Y, Kido T, et al. Relation of long-term body weight change to change in lipoprotein particle size in Japanese men and women: the INTERMAP Toyama Study. Atherosclerosis. 2009;206:282–6.View ArticlePubMedGoogle Scholar
- Mantyselka P, Kautiainen H, Saltevo J, Wurtz P, Soininen P, Kangas AJ, et al. Weight change and lipoprotein particle concentration and particle size: a cohort study with 6.5-year follow-up. Atherosclerosis. 2012;223:239–43.View ArticlePubMedGoogle Scholar
- Inouye M, Kettunen J, Soininen P, Silander K, Ripatti S, Kumpula LS, et al. Metabonomic, transcriptomic, and genomic variation of a population cohort. Mol Syst Biol. 2010;6:441.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhou M, Wang S, Zhao A, Wang K, Fan Z, Yang H, et al. Transcriptomic and metabonomic profiling reveal synergistic effects of quercetin and resveratrol supplementation in high fat diet fed mice. J Proteome Res. 2012;11:4961–71.View ArticlePubMedGoogle Scholar
- Xu T, Holzapfel C, Dong X, Bader E, Yu Z, Prehn C, et al. Effects of smoking and smoking cessation on human serum metabolite profile: results from the KORA cohort study. BMC Med. 2013;11:60.View ArticlePubMedPubMed CentralGoogle Scholar
- Holle R, Happich M, Lowel H, Wichmann HE, MONICA KORA Study Group: KORA--a research platform for population based health research. Gesundheitswesen. 2005; Suppl 1:19–25.
- Wolfenstetter SB, Menn P, Holle R, Mielck A, Meisinger C, von Lengerke T. Body weight changes and outpatient medical care utilisation: results of the MONICA/KORA cohorts S3/F3 and S4/F4. Psychosoc Med. 2012;9:Doc09. doi:10.3205/psm000087.
- Rathmann W, Haastert B, Icks A, Lowel H, Meisinger C, Holle R, et al. High prevalence of undiagnosed diabetes mellitus in Southern Germany: target populations for efficient screening. The KORA survey 2000. Diabetologia. 2003;46:182–9.PubMedGoogle Scholar
- Mühlberger N, Behrend C, Stark R, Holle R. Datenbankgestützte Online-Erfassung von Arzneimitteln im Rahmen gesundheitswissenschaftlicher Studien - Erfahrungen mit der IDOM-Software. Informatik, Biometrie und Epidemiologie in Medizin und Biologie. 2003;34:601–11.Google Scholar
- DeHaven CD, Evans AM, Dai H, Lawton KA. Organization of GC/MS and LC/MS metabolomics data into chemical libraries. J Cheminform. 2010;2:9.View ArticlePubMedPubMed CentralGoogle Scholar
- Evans AM, DeHaven CD, Barrett T, Mitchell M, Milgram E. Integrated, nontargeted ultrahigh performance liquid chromatography/electrospray ionization tandem mass spectrometry platform for the identification and relative quantification of the small-molecule complement of biological systems. Anal Chem. 2009;81:6656–67.View ArticlePubMedGoogle Scholar
- Soininen P, Kangas AJ, Wurtz P, Tukiainen T, Tynkkynen T, Laatikainen R, et al. High-throughput serum NMR metabonomics for cost-effective holistic studies on systemic metabolism. Analyst. 2009;134:1781–5.View ArticlePubMedGoogle Scholar
- Doring A, Gieger C, Mehta D, Gohlke H, Prokisch H, Coassin S, et al. SLC2A9 influences uric acid concentrations with pronounced sex-specific effects. Nat Genet. 2008;40:430–6.View ArticlePubMedGoogle Scholar
- Schurmann C, Heim K, Schillert A, Blankenberg S, Carstensen M, Dorr M, et al. Analyzing illumina gene expression microarray data from different tissues: methodological aspects of data analysis in the metaxpress consortium. PLoS One. 2012;7:e50938.View ArticlePubMedPubMed CentralGoogle Scholar
- R Core Team. R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. 2014. http://www.R-project.org.
- Bioconductor: Open Source Software for Bioinformatics. http://www.bioconductor.org.
- van Buuren S, Groothuis-Oudshoorn K. mice: Multivariate Imputation by Chained Equations in R. J Stat Softw. 2011;45:1-67.
- Chuang HY, Lee E, Liu YT, Lee D, Ideker T. Network-based classification of breast cancer metastasis. Mol Syst Biol. 2007;3:140.View ArticlePubMedPubMed CentralGoogle Scholar
- Inouye M, Ripatti S, Kettunen J, Lyytikainen LP, Oksala N, Laurila PP, et al. Novel Loci for metabolic networks and multi-tissue expression studies reveal genes for atherosclerosis. PLoS Genet. 2012;8:e1002907.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang B, Horvath S. A general framework for weighted gene co-expression network analysis. Stat Appl Genet Mol Biol. 2005;4:Article 17. Epub 2005 Aug 12.
- Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008;9:559.View ArticlePubMedPubMed CentralGoogle Scholar
- Langfelder P, Zhang B, Horvath S. Defining clusters from a hierarchical cluster tree: the Dynamic Tree Cut package for R. Bioinformatics. 2008;24:719–20.View ArticlePubMedGoogle Scholar
- Krumsiek J, Suhre K, Illig T, Adamski J, Theis FJ. Bayesian independent component analysis recovers pathway signatures from blood metabolomics data. J Proteome Res. 2012;11:4120–31.View ArticlePubMedGoogle Scholar
- Wahl S, Krug S, Then C, Kirchhofer A, Kastenmüller G, Brand T, et al. Comparative analysis of plasma metabolomics response to metabolic challenge tests in healthy subjects and influence of the FTO obesity risk allele. Metabolomics. 2013;10:386–401.Google Scholar
- Ingenuity. http://www.ingenuity.com.
- Krumsiek J, Suhre K, Illig T, Adamski J, Theis FJ. Gaussian graphical modeling reconstructs pathway reactions from high-throughput metabolomics data. BMC Syst Biol. 2011;5:21.View ArticlePubMedPubMed CentralGoogle Scholar
- Klop B, Elte JW, Cabezas MC. Dyslipidemia in obesity: mechanisms and potential targets. Nutrients. 2013;5:1218–40.View ArticlePubMedPubMed CentralGoogle Scholar
- Brunzell J, Hokanson J. Dyslipidemia of central obesity and insulin resistance. Diabetes Care. 1999;22:C10–3.View ArticlePubMedGoogle Scholar
- Tzotzas T, Desrumaux C, Lagrost L. Plasma phospholipid transfer protein (PLTP): review of an emerging cardiometabolic risk factor. Obes Rev. 2009;10:403–11.View ArticlePubMedGoogle Scholar
- Tukiainen T, Kettunen J, Kangas AJ, Lyytikainen LP, Soininen P, Sarin AP, et al. Detailed metabolic and genetic characterization reveals new associations for 30 known lipid loci. Hum Mol Genet. 2012;21:1444–55.View ArticlePubMedGoogle Scholar
- Freedman DS, Otvos JD, Jeyarajah EJ, Barboriak JJ, Anderson AJ, Walker JA. Relation of lipoprotein subclasses as measured by proton nuclear magnetic resonance spectroscopy to coronary artery disease. Arterioscler Thromb Vasc Biol. 1998;18:1046–53.View ArticlePubMedGoogle Scholar
- Musunuru K, Orho-Melander M, Caulfield MP, Li S, Salameh WA, Reitz RE, et al. Ion mobility analysis of lipoprotein subfractions identifies three independent axes of cardiovascular risk. Arterioscler Thromb Vasc Biol. 2009;29:1975–80.View ArticlePubMedPubMed CentralGoogle Scholar
- Mora S, Otvos JD, Rosenson RS, Pradhan A, Buring JE, Ridker PM. Lipoprotein particle size and concentration by nuclear magnetic resonance and incident type 2 diabetes in women. Diabetes. 2010;59:1153–60.View ArticlePubMedPubMed CentralGoogle Scholar
- Felig P, Marliss E, Cahill GJ. Plasma amino acid levels and insulin secretion in obesity. N Engl J Med. 1969;281:811–6.View ArticlePubMedGoogle Scholar
- Lu J, Xie G, Jia W, Jia W. Insulin resistance and the metabolism of branched-chain amino acids. Front Med. 2013;7:53–9.View ArticlePubMedGoogle Scholar
- Wurtz P, Tiainen M, Makinen VP, Kangas AJ, Soininen P, Saltevo J, et al. Circulating metabolite predictors of glycemia in middle-aged men and women. Diabetes Care. 2012;35:1749–56.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang TJ, Larson MG, Vasan RS, Cheng S, Rhee EP, McCabe E, et al. Metabolite profiles and the risk of developing diabetes. Nat Med. 2011;17:448–53.View ArticlePubMedPubMed CentralGoogle Scholar
- Wiklund PK, Pekkala S, Autio R, Munukka E, Xu L, Saltevo J, et al. Serum metabolic profiles in overweight and obese women with and without metabolic syndrome. Diabetol Metab Syndr. 2014;6:40.View ArticlePubMedPubMed CentralGoogle Scholar
- Demine S, Reddy N, Renard P, Raes M, Arnould T. Unraveling biochemical pathways affected by mitochondrial dysfunctions using metabolomic approaches. Metabolites. 2014;4:831–78.View ArticlePubMedPubMed CentralGoogle Scholar
- Rogge MM. The role of impaired mitochondrial lipid oxidation in obesity. Biol Res Nurs. 2009;10:356–73.View ArticlePubMedGoogle Scholar
- Kelley D, He J, Menshikova E, Ritov V. Dysfunction of mitochondria in human skeletal muscle in type 2 diabetes. Diabetes. 2002;51:2944–50.View ArticlePubMedGoogle Scholar
- Keun HC, Sidhu J, Pchejetski D, Lewis JS, Marconell H, Patterson M, et al. Serum molecular signatures of weight change during early breast cancer chemotherapy. Clin Cancer Res. 2009;15:6716–23.View ArticlePubMedGoogle Scholar
- Mailloux R, Lemire J, Appanna V. Aluminum-induced mitochondrial dysfunction leads to lipid accumulation in human hepatocytes: a link to obesity. Cell Physiol Biochem. 2007;20:627–38.View ArticlePubMedGoogle Scholar
- Clarke C, Xiao R, Place E, Zhang Z, Sondheimer N, Bennett M, et al. Mitochondrial respiratory chain disease discrimination by retrospective cohort analysis of blood metabolites. Mol Genet Metab. 2013;110:145–52.View ArticlePubMedGoogle Scholar
- Shaham O, Slate NG, Goldberger O, Xu Q, Ramanathan A, Souza AL, et al. A plasma signature of human mitochondrial disease revealed through metabolic profiling of spent media from cultured muscle cells. Proc Natl Acad Sci U S A. 2010;107:1571–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Goh SH, Josleyn M, Lee YT, Danner RL, Gherman RB, Cam MC, et al. The human reticulocyte transcriptome. Physiol Genomics. 2007;30:172–8.View ArticlePubMedGoogle Scholar
- Trottier MD, Naaz A, Li Y, Fraker PJ. Enhancement of hematopoiesis and lymphopoiesis in diet-induced obese mice. Proc Natl Acad Sci U S A. 2012;109:7622–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Virtue M, Furne J, Nuttall F, Levitt M. Relationship between GHb concentration and erythrocyte survival determined from breath carbon monoxide concentration. Diabetes Care. 2004;27:931–5.View ArticlePubMedGoogle Scholar
- Kendall RG, Mellors I, Hardy J, Mardle B. Patients with pulmonary and cardiac disease show an elevated proportion of immature reticulocytes. Clin Lab Haematol. 2001;23:27–31.View ArticlePubMedGoogle Scholar
- Joehanes R, Ying S, Huan T, Johnson AD, Raghavachari N, Wang R, et al. Gene expression signatures of coronary heart disease. Arterioscler Thromb Vasc Biol. 2013;33:1418–26.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang Y, Ma AQ, Gong M, Lu Q, Lu M, Tian G. Red blood cell level is increased in obese but not in non-obese patients with coronary heart disease. J Geriatr Cardiol. 2010;7:143–6.Google Scholar
- Inouye M, Silander K, Hamalainen E, Salomaa V, Harald K, Jousilahti P, et al. An immune response network associated with blood lipid levels. PLoS Genet. 2010;6:e1001113.View ArticlePubMedPubMed CentralGoogle Scholar
- Sokol CL, Barton GM, Farr AG, Medzhitov R. A mechanism for the initiation of allergen-induced T helper type 2 responses. Nat Immunol. 2008;9:310–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Imtiyaz HZ, Williams EP, Hickey MM, Patel SA, Durham AC, Yuan LJ, et al. Hypoxia-inducible factor 2alpha regulates macrophage function in mouse models of acute and tumor inflammation. J Clin Invest. 2010;120:2699–714.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang H, Zhang G, Gonzalez FJ, Park SM, Cai D. Hypoxia-inducible factor directs POMC gene to mediate hypothalamic glucose sensing and energy balance regulation. PLoS Biol. 2011;9:e1001112.View ArticlePubMedPubMed CentralGoogle Scholar
- Gonen B, O’Donnell P, Post T, Quinn T, Schulman E. Very low density lipoproteins (VLDL) trigger the release of histamine from human basophils. Biochim Biophys Acta. 1987;917:418–24.View ArticlePubMedGoogle Scholar
- Shore SA. Obesity, airway hyperresponsiveness, and inflammation. J Appl Physiol. 2010;108:735–43.View ArticlePubMedGoogle Scholar
- Johansson LE, Danielsson AP, Parikh H, Klintenberg M, Norstrom F, Groop L, et al. Differential gene expression in adipose tissue from obese human subjects during weight loss and weight maintenance. Am J Clin Nutr. 2012;96:196–207.View ArticlePubMedGoogle Scholar
- Valcarcel B, Ebbels TM, Kangas AJ, Soininen P, Elliot P, Ala-Korpela M, et al. Genome metabolome integrated network analysis to uncover connections between genetic variants and complex traits: an application to obesity. J R Soc Interface. 2014;11:20130908.View ArticlePubMedPubMed CentralGoogle Scholar
- Siri-Tarino PW, Williams PT, Fernstrom HS, Rawlings RS, Krauss RM. Reversal of small, dense LDL subclass phenotype by normalization of adiposity. Obesity. 2009;17:1768–75.View ArticlePubMedPubMed CentralGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.