Skip to main content
  • Research article
  • Open access
  • Published:

Gene expression profiling in whole blood identifies distinct biological pathways associated with obesity

Abstract

Background

Obesity is reaching epidemic proportions and represents a significant risk factor for cardiovascular disease, diabetes, and cancer.

Methods

To explore the relationship between increased body mass and gene expression in blood, we conducted whole-genome expression profiling of whole blood from seventeen obese and seventeen well matched lean subjects. Gene expression data was analyzed at the individual gene and pathway level and a preliminary assessment of the predictive value of blood gene expression profiles in obesity was carried out.

Results

Principal components analysis of whole-blood gene expression data from obese and lean subjects led to efficient separation of the two cohorts. Pathway analysis by gene-set enrichment demonstrated increased transcript levels for genes belonging to the "ribosome", "apoptosis" and "oxidative phosphorylation" pathways in the obese cohort, consistent with an altered metabolic state including increased protein synthesis, enhanced cell death from proinflammatory or lipotoxic stimuli, and increased energy demands. A subset of pathway-specific genes acted as efficient predictors of obese or lean class membership when used in Naive Bayes or logistic regression based classifiers.

Conclusion

This study provides a comprehensive characterization of the whole blood transcriptome in obesity and demonstrates that the investigation of gene expression profiles from whole blood can inform and illustrate the biological processes related to regulation of body mass. Additionally, the ability of pathway-related gene expression to predict class membership suggests the feasibility of a similar approach for identifying clinically useful blood-based predictors of weight loss success following dietary or surgical interventions.

Peer Review reports

Background

While excess energy intake and declining energy expenditure are clearly important contributors, individual susceptibility to obesity is also strongly influenced by genetic factors. Twin, adoption, and family studies have indicated that 40-70% of inter-individual variation in body mass index (BMI) is heritable [1, 2]. A compendium of evidence for the genetic bases of obesity have been accrued from single-gene mutation studies, Mendelian inheritance patterns, transgenic and knockout murine models, animal and human quantitative trait loci (QTL), candidate-gene association studies, and genome scan linkages and have been incorporated into the Obesity Gene Map database [3]. Also recently, a number of genome-wide association studies (GWAS) have demonstrated associations of single-nucleotide polymorphisms (SNPs) to qualitative and quantitative indices of adiposity in several populations [2, 4–10]. A combination of independent studies and meta-analysis of existing GWAS data have implicated a total of 18 genetic loci as relevant for body weight regulation to date [11].

In addition to DNA sequence variants, genetic influences are also manifested through differences in gene transcription, leading to differential messenger RNA levels. While such differences might be expected to occur in biologically relevant tissues (muscle and adipose tissue in obesity, for example), several recent studies have demonstrated an alteration in the peripheral blood transcriptome in diseases of non-hematologic origin. These include disorders such as chronic fatigue syndrome, schizophrenia and colon cancer [12–17]. Additionally, the blood transcriptome has also been found to be responsive to diverse environmental and socio-economic stimuli including ionizing radiation in cancer therapy, benzene exposure, socio-economic status, etc. [18–21]. These findings raise the intriguing possibility that blood transcriptome profiles might provide a valid biological readout for otherwise hard to study disease processes in humans and additionally generate information of high predictive and diagnostic content. In line with this argument, we postulated that differences in transcript abundance might also occur in blood from obese subjects compared to lean subjects, as a consequence of either pre-existing genetic variations, or as an adaptive response to obesity, independent of the genetic background. To test this hypothesis, we have carried out transcriptional profiling of peripheral blood from obese subjects and well-matched lean controls and conducted enrichment analysis to identify biological pathways that are preferentially associated with obesity. Our study demonstrates significant gene expression differences in blood from obese subjects compared to lean controls, particularly along the lines of differential expression of genes in key metabolic pathways regulating cell survival, protein synthesis and energy harvest. These findings are important on three levels. First, our results demonstrate the importance of blood as a biologically informative tissue in the elucidation of the obese state. Second, as differences in gene expression are often driven by sequence variants in gene regulatory regions, our study provides a mechanism for the selection of obesity-associated candidate genes for the determination of possible regulatory sequence variants. Finally, the identification of adiposity related gene expression differences in a clinically accessible tissue such as blood leads the way for the determination of biomarkers of weight regulation that could be implemented in a clinical setting.

Results

Phenotypic characterization of study subjects

Demographic and phenotypic characteristics of the subjects included in the current study are shown in Table 1. The obese and lean subjects showed statistically significant differences (p < 0.05 level) in almost all metabolic parameters tested with the exception of cholesterol, LDL-cholesterol and thyroid stimulating hormone status. Also, levels of glycated haemoglobin (HbA1c), insulin and fasting glucose were statistically significantly different but within normal clinical ranges in both groups. There were 53% (9/17) and 70% (12/17) females in the obese and lean groups, respectively. Both cohorts were closely matched for age and hormonal status (6/12 lean and 5/9 obese women were postmenopausal).

Table 1 Demographic and phenotypic characteristics of the study population.

Ascertainment of data quality

We ascertained the overall quality of the whole genome expression profiling signals by comparing the Affymetrix microarray generated expression patterns of a subset of 61 genes (with a 20% or greater change in expression between obese and lean cohorts) to expression signals generated by real-time, quantitative PCR (Taqman). The genes selected cover a range of approximately 7 logs (base 2) representing over 100-fold differences in the magnitude of gene expression on Affymetrix microarrays (average log2 signal of 4.66 for protein tyrosine phosphatase, receptor type, S to an average log2 signal of 11.35 for RAB31 gene in the obese cohort). The results are shown in Figure 1. In ~ 75% of genes tested (45/61 genes), the direction of gene expression changes between obese and lean subjects were in agreement between the Affymetrix and Taqman platforms suggesting high reproducibility of gene expression data between the two approaches. Additionally, analysis of muscle and adipose-specific marker gene expression demonstrated no evidence of contamination in the study samples (Additional File 1).

Figure 1
figure 1

Comparison of gene expression signals generated by Affymetrix microarrays and quantitative real-time PCR. Gene expression signals were generated by real-time, quantitative PCR (Taqman, black bars) and oligonucleotide microarrays (Affymetrix, white bars). Overexpression or underexpression of a gene in the obese and lean cohorts is expressed as a log ratio, to the base 2. Affymetrix and Taqman based results for each gene are shown as a stacked bar. For each gene, agreement between the results from the two platforms is indicated when both white and black bars lie on the same side of the zero (0) value on the log ratio axis; conversely, disagreement is indicated when the gray and black bars lie on opposite sides. The overall agreement between the two platforms was 85% (45/53 genes showed agreement in the direction of differential expression).

Principal components analysis of gene expression data

We performed multivariate, principal components analysis to determine whether blood gene expression signals were capable of distinguishing between the obese and lean subjects. Figure 2 shows a scatterplot representing the first two principal components based on gene expression profiles from 17 obese and 17 lean subjects. Analysis of the principal component model performance indicated that 27% of the total variance in gene expression was modelled in the first principal component (R2X) with a cross-validated prediction of 22.4%. The cross-validation results indicate that the variability captured in the first component is statistically greater than the significance limit of 2.9% (Additional file 2).

Figure 2
figure 2

Multivariate analysis of obese and lean subjects based on gene expression signals. Principal component analysis (PCA) was performed on lean and obese subjects based on 12128 Affymetrix probe-set signals. A scatterplot of the first two principal components demonstrate a general separation of the obese and lean phenotypes along the first principal component (PC1). Model parameters are as follows: Further details on the PCA model parameters are included in Supplemental Table 2.

Identification of differentially expressed genes

Genes showing differential expression between the obese and lean subjects were identified via the Comparative Marker Selection module in GenePattern [22], using the signal-to-noise algorithm for ranking genes. A permutation testing was performed to compute the significance (nominal p-value) of the rank assigned to each gene. A false discovery rate (FDR) was also calculated to control for multiple testing. A total of 12127 probesets were detected above background (set to 50 units) among which 374 probesetes were overexpressed (2-fold or greater) and 75 probesets were underexpressed (2-fold or greater) in the obese samples compared to the leans. The results of the differential gene analysis are presented in Additional Files 3 and 4. Inspection of the gene list showed that a majority of the genes upregulated in the obese subjects were genes known to be selectively expressed in erythrocytes/reticulocytes. These included genes such as carbonic anhydrase, ferrochelatase, synuclein, glycophorin B, etc. This finding is consistent with previous observations of higher red blood cell counts (hematocrit) in obesity [23–26] and provides evidence for the expansion of transcriptionally active reticulocytes in obesity. Conversely, several genes related to immune function showed reduced expression in the obese subjects.

Pathway analysis of gene expression difference between lean and obese subjects

The transcriptome data was next subjected to bioinformatic pathway analysis by the Gene Set Enrichment Analysis (GSEA) algorithm [27]. The values for the GSEA algorithmic parameters used in the current study are indicated in Additional File 5 and details about the GSEA algorithm have been explained in Materials and Methods. Pathway analysis was conducted either with the Kyoto Encyclopedia for Genes and Genomes (KEGG) metabolic pathway database [28], or a user-created custom database consisting of pathways drawn from several sources (Additional File 6). Pathways were evaluated by their normalized enrichment score (NES), nominal p-values (permuted) and false discovery rates, as described in [29].

KEGG Pathway analysis

Enrichment analysis of gene expression profiles against KEGG pathways identified 5 pathways at p permuted < 0.05 level (Additional File 7). Notable among them were the 'apoptosis', 'ribosome', and 'oxidative phosphorylation' pathways. The pathway enrichment plots and expression profiles of a subset of genes contributing significantly to the enrichment of these 3 pathways are collectively shown in Figure 3. A number of genes, including apoptotic protease activating factor 1, baculoviral IAP repeat containing 2, caspase 7, Fas, interleukin 1 beta, interleukin 1 receptor associated kinase 4, etc. contributed to the core enrichment of the 'apoptosis' pathway in the obese subjects (Additional File 8). Enrichment of the ribosome pathway was effected by coordinate upregulation of several ribosomal protein genes (ribosomal protein L31, S7, S24, L35, L7 for example). Several genes involved in the mitochondrial process of electron transfer and ATP synthesis demonstrated increased expression in the obese cohort leading to a significant enrichment of the 'oxidative phosphorylation' pathway in this group. Some of the genes contributing to core enrichment of this pathway included cytochrome c oxidase subunits 6C, 7B and 7C, NADH-coenzyme Q reductase, NADH deyhdrogenase beta subcomplex 3, etc.

Figure 3
figure 3

Gene-set enrichment analysis. Gene-set enrichment analysis against the KEGG database for differentially enriched pathways in whole blood between obese and lean subjects. Enrichment plots for the 3 pathways upregulated in the obese cohort are shown on the left side with the relative gene positions indicated by the straight lines (line plot) under each graph. Lines clustered to the left represent higher ranked genes in the ranked list. Expression profiles for a subset of genes (shaded in yellow in the line plots) contributing to core enrichment for each pathway are shown to the right as a heatmap. The heatmap compares subject-level gene expression in both obese and lean subjects. Gene expression is normalized for each row. Lower levels of expression are represented in shades of blue and higher expression in red.

Custom Pathway analysis

In addition to investigating pathway enrichment based on the KEGG database, we also subjected a set of 'custom' pathways to analysis by GSEA (Additional File 6). GSEA analysis of the custom pathways identified 2 pathways as significantly upregulated in the obese, at a nominal p-value < 5% and FDR < 5%. These were the 'electron transport chain pathway' and the 'erythrocyte/reticulocytespecific_affytechnote' pathways (Additional File 9). The 'electron transport chain pathway' (National Cancer Institute Pathway Interaction Database [30]) is a subset of the KEGG 'oxidative phosphorylation' pathway. The 'erythrocyte/reticulocytespecific_affytechnote' pathway consists of genes reported to be selectively enriched for expression in erythrocytes/reticulocytes (Affymetrix, [31, 32]). Identification of this gene-set as an obesity-upregulated pathway further supports our earlier observation of increased expression of individual erythrocyte/reticulocyte specific genes in the obese subjects. Details are provided in Additional Files 10 and 11.

Effects of gender on pathway enrichment

Since our study cohort contained both male and female subjects, the contribution of gender to pathway enrichment was investigated. To determine whether pathway ranks were influenced by gender, we carried out independent gene-set enrichment analyses on subgroups comprised of female or male subjects only. We compared the relative ranks of the KEGG pathways in the three analyses as an indication of their sensitivity to gender. 'Apoptosis' was ranked 7th, 8th and 3rd and 'oxidative phosphorylation' was ranked 10th, 12th and 18th for All subjects, Females and Males respectively. The 'ribosome' pathway was the top ranked pathway for All subjects and Females analysis, but was ranked 27th in the analysis involving the Males. We repeated the same subgroup analyses on the custom pathway set and in all cases the 'electron transport chain pathway' and the 'erythrocyte/reticulocytespecific_affytechnote' pathways remained the top 2 ranked pathways for all groups tested. Details are provided in Additional File 12.

Effect of cell populations on pathway enrichment

Since whole-blood consists of a mixture of various cell types, we investigated the relation between the observed enrichment in "ribosome", "apoptosis" and "oxidative phosphorylation" pathways in the obese and enrichment of reticulocytes/erythrocytes in obese subjects as previously reported [23–26]. We scaled the gene expression data independently by the expression of 2 erythrocyte-specific transcripts, hemoglobin D (HBD) and erythrocyte membrane protein, band 2 (EMPB2) and subjected the scaled data to gene-set enrichment analysis. Of the original 3 pathways found to be enriched in the obese subjects, the "ribosome" pathway was still the top differentially expressed pathway with both unscaled and scaled data. However, the "apoptosis" and "oxidative phosphorylation" pathways were no longer significantly enriched, with either of the scaled datasets. Pathway enrichment results with scaled data are provided in Additional File 13.

Class prediction via blood gene expression

We next examined whether biological pathways implicated from gene-set enrichment analysis of the current study could provide a set of mechanism-based gene predictors that would be capable of predicting obese and lean subjects with high accuracy. We created an initial, inclusive set containing all genes (features) belonging to the ribosome, apoptosis or oxidative phosphorylation pathways (183 genes). Since this list was also likely to contain redundant and non-informative genes, we applied two independent feature selection algorithms to identify a smaller set of genes that would be capable of distinguishing between the obese and lean phenotypes with high success rates, based on the metrics specific to the two algorithms used (described in detail in Materials and Methods). A search for overlapping genes scoring high in both algorithms (ranked within the top 20 genes in both) resulted in a set of 11 genes. The logged gene expression signals from the full (183) and filtered (11) gene-sets were then used as inputs into four different classifiers representing distinct algorithmic approaches to classification and prediction. These included the Naive Bayes, Logistic Regression, Random Forests and ZeroR classifiers. A full description of the classifiers is presented in Materials and Methods and Additional File 14. Each classifier was first trained on a randomly selected 66% of the samples and then used to predict the class for the remaining 33% samples. The process was repeated 100 times for each classifier. Classifier performance was evaluated by four parameters (true positive, true negative, false positive and false negative rates). A description of the performance evaluators can be found in Additional File 14. The classifier ZeroR simply predicts the same class for all instances and was used as a baseline classifier. Any classifier should perform significantly better than ZeroR in order to be considered useful. Table 2 compares the performance of the four classifiers with either the full gene-set (183 genes) or the filtered set (11 genes). For each of the four performance evaluators, we plotted the average and standard deviation values for the four parameters over the 100 iterations. Overall, the Naïve Bayes and logistic regression classifiers performed better than the decision-tree based classifier (Random Forests) and all three classifiers performed significantly better than ZeroR. A comparison of the classifier results with the full (183) or filtered (11) gene-set inputs showed that both inputs had similar true positive and false negative rates. Both Naïve Bayes and logistic regression classifiers displayed high sensitivity as indicated by true positive rates close to 1.0. These two classifiers also demonstrated lower false-positive rates with the filtered gene set compared to the full gene set. Additionally, the filtered gene set classifiers displayed higher specificities (true negative rates) compared to the full gene set based classifiers. Based on these results, we found the 11-gene based Naïve Bayes or logistic regression based classifiers to perform better compared to the 183-gene classifiers for predicting class membership. The identities of the 11 genes are shown in Table 3 and appear to be primarily composed of genes from the oxidative phosphorylation and apoptosis pathways.

Table 2 Classification of lean and obese subjects
Table 3 Identity of genes constituting the 11-gene classifier.

Discussion

Our study demonstrates significant gene expression differences in whole blood from age-matched obese and lean subjects of Northern European White genetic ancestry. These differences further lead to the identification of differentially enriched biological pathways in obesity and lead to an increased appreciation and understanding of genomic changes in whole blood related to body mass expansion. The current study is not designed to resolve whether the observed transcriptional differences are causal or caused, i.e. whether the differences in gene expression are related to the development of obesity or reflect an adaptive mechanism in response to increased body mass. Although blood is usually not considered to be a target organ for obesity, certain observations are pertinent. First, the physiological role of blood as a sentinel tissue and a systemic integrator of tissue and organ-level perturbations could lead to adaptive responses in response to major metabolic perturbations such as excessive build-up of body mass and the attendant increases in the demand for nutrient and oxygen transport. Secondly, the chronic low-grade tissue inflammation observed in obesity [33] is expected to have a direct effect on circulating leukocytes, including immune dysfunction and apoptosis. Finally, macrophages in blood share many functional and antigenic properties with preadipocytes and adipocytes and transcriptome profiles of preadipocytes are reportedly closer to the macrophages than to adipocytes [34]. In this context, our study provides the first detailed investigation of the blood transcriptome in relation to obesity and provides evidence in favor of its dynamic involvement in the process. It is important to note here that the between-group differences in gene expression were usually small and there was considerable heterogeneity in individual gene expression values among subjects in the obese or lean categories. However, the between-group variation exceeded the within-group variation for several genes leading to statistically significant differences between the groups. Additionally, as demonstrated by principal components analysis, blood gene expression profiles were able to distinguish lean subjects from obese subjects even when the subject classes were not exposed a priori (unsupervised clustering). Since gene expression measures were used as input for the PCA analysis, these results suggest that the differences in blood transcript levels between obese and lean subjects were significant and informative enough to cause a separation between the two classes.

The application of pathway analysis provided additional information and insight into the biological processes that are differentially regulated in obese and lean blood samples. Some of the pathways with increased component transcript abundances included the "ribosome", "apoptosis" and "oxidative phosphorylation" pathways. Upregulation of the ribosomal pathway in the obese subjects was due to an increased expression of several ribosomal protein-encoding genes, indicative of enhanced protein synthesis in blood cells, possibly as a consequence of enhanced metabolic demands in the obese state. This observation is consistent with a recent report that links ribosomal RNA synthesis to cellular energy supply through activation of the AMP-activated protein kinase [35]. The presence of increased apoptosis in the obese phenotype has also been well documented in animal and human cell culture models. For example, increased cardiomyocyte apoptosis has been reported in leptin-deficient ob/ob mice and leptin-resistant db/db mice [36]. Prolonged exposure to free fatty acids also have pro-apoptotic effects on human pancreatic islets [37] and circulating cytokines, such as tumor necrosis factor alpha (TNF-α) have been reported to induce apoptosis in cultured human preadipocytes and adipocytes [38]. Our findings now provide evidence for activation of a similar apoptotic program in blood from obese subjects. While the current study does not allow us to pinpoint the cause of the enhanced apoptosis, we speculate that obesity-associated chronic inflammation [39, 40] or lipotoxicity are contributing factors. Finally, the observed upregulation of the 'oxidative phosphorylation' pathway in obese subjects is consistent with a response to increased energy demands in obese subjects. Functional and gene expression studies have previously indicated impairment in oxidative phosphorylation and mitochondrial function in subjects with type 2 diabetes compared to controls [29, 41, 42]. Our findings are consistent with Takamura et al., who demonstrated an upregulation of oxidative phosphorylation genes in the livers of obese, type 2 diabetic patients compared to non-obese diabetics [43]. More interestingly, our findings now point to a similar involvement of energy-harvesting mechanisms in obese blood and provide further evidence in favor of a role for mitochondrial dysfunction in obesity [44, 45]. A gender-based sub-analysis demonstrated relative stability of the "apoptosis" and "oxidative phosphorylation" pathway ranks in both genders; in contrast, the "ribosome" pathway differed significantly in rank between females and males, suggesting a gender-specific effect (Additional File 7). Since a majority of genes upregulated in the obese subjects are highly expressed in erythrocytes and reticulocytes, we scaled the gene expression data independently by the expression of two erythrocyte-specific transcripts, hemoglobin D (HBD) and erythrocyte membrane protein, band 2 (EMPB2) and subjected the scaled data to gene-set enrichment analysis. Of the three pathways found to be differentially upregulated in the obese subjects, the "ribosome" pathway remained the top differentially expressed pathway (with the scaled data) whereas the "apoptosis" and "oxidative phosphorylation" pathways were no longer significantly enriched, with either of the scaled datasets. These findings suggest that an increase in erythrocyte/reticulocyte numbers in the obese (differential hematocrit) is a possible explanatory mechanism for the observed increase in transcript levels for "apoptosis" and "oxidative phosphorylation" in the obese subjects. The results for the "ribosome" pathway, in contrast, suggest a significant upregulation of the transcripts for the component genes of this pathway in the obese subjects, even after adjustment for erythrocyte-specific gene expression. We note one caveat to the scaling approach used here for investigating cell number effects. Since the same amount of cRNA was used from each sample for hybridization, the relative enrichment of cell types is expected to have a real effect on gene expression only for genes that are differentially expressed among the cell types (e.g. hemoglobin transcripts that are expressed only in reticulocytes and not lymphocytes). For genes expressed at comparable levels across cell types, the differential cell type representation should not have an effect on expression unless there is a true upregulation or downregulation of these genes between the two groups (although the cellular origin for the differential expression may not be known). Scaling the gene expression data by the expression of reticulocyte/erythrocyte specific genes cannot distinguish between the above two mechanisms of enhanced gene expression and can lead to potentially incorrect conclusions. However, our results clearly demonstrate that inter-individual variations in hematocrit, especially between obese and lean subjects, may affect interpretation of expression data and should be considered as an important co-variate in future studies.

Several recent publications have reported on the successful application of gene expression signatures as classifiers or predictors of phenotypic class, disease progression and therapeutic prognosis, primarily in the area of diagnosis and treatment of several types of cancers [16, 46–48]. However, the biological mechanisms linking the predictive genes to the outcomes being predicted are not always clear. This lack of mechanism has often been criticized as a barrier to the clinical utility of the gene predictors. One solution to the problem is to choose gene predictors from biological pathways associated a priori with the phenotype or outcome of interest. This approach was pursued in this study and led to the identification of an 11-gene based classifier that could distinguish and predict obese and lean subjects with high accuracy. Our motivation for this exercise was to provide proof-of-concept data to test if blood gene expression patterns can have predictive value in the context of obesity. While such prediction is not necessarily required for distinguishing obesity from leanness, blood based gene biomarkers can significantly advance the clinical management of obesity by, for example, allowing the prediction of weight loss success from diet or bariatric surgery.

One potential downstream application of differential gene expression analysis in whole-blood is the selection of candidate genes with possible regulatory polymorphisms (single nucleotide polymorphisms in promoter regions, for example) that associate with obesity and help explain the observed differences in expression. Comprehensive sequencing of the regulatory regions of such candidate genes are expected to yield additional insights into the genetics of obesity such as the identification of expression QTLs (eQTLs). While a direct subject-level association of gene regulatory polymorphisms to gene expression levels is outside the scope of the current work, we conducted a preliminary analysis of the existence of putative regulatory variants in the 11 gene predictors identified in our analysis. Based on data from the NCBI dbSNP database (Build 131), several genes contained common sequence variants near the 5'-end of the gene spanning a region 2000 bases upstream of the start codon (SNPs rs2515192 and rs3019164 for ATP6V1C1, rs1317775 and rs1318199 for BIRC2, rs11709092 for PRKAR2A, etc.). It is reasonable to speculate that a subset of these upstream sequence variants could influence transcription.

Our study relied on whole-blood collected in PAXgene tubes instead of peripheral blood mononuclear cells (PBMCs), consistent with our ultimate goal of identifying clinically relevant and useful predictors of weight loss success. This procedure, however, has the disadvantage of investigating a relatively heterogeneous cell population where noise could mask gene expression differences in specific cell types. PBMC's, consisting of lymphocytes and monocytes provide a consistent and homogeneous sample for transcriptome analysis. However, the extra fractionation procedure for PBMCs requires a prolonged period before RNA stabilization leading to significant ex vivo changes in gene expression profiling [49]. Additionally, compared to whole blood, several cell types including neutrophils, basophils, eosinophils, platelets, reticulocytes and erythrocytes are depleted in PBMCs which lead to loss of important transcription information. On the other hand, PAX samples show a decrease in the number of expressed genes and lower gene expression values with higher variability compared to the PBMCs [50], primarily due to the high abundance of globin transcripts that constitute over 70% of whole blood mRNA [51]. However, the PAXgene system employs an easy way to collect, store, transport and stabilize RNA from whole blood and based on our overall goals, was the method of choice for our analysis. In this context, the ability of gene expression signatures from biologically relevant pathways to accurately classify and predict obese and lean classes, as observed in this study, provides further validation of our approach and suggests future suitability of the PAXgene based whole blood transcriptome for yielding clinically usable biomarkers related to weight regulation. Additional sensitivity could be obtained in future studies via selective reduction of the globin transcript from whole blood RNA samples [52, 53].

There are the following limitations to the current study. First, since the study employed whole blood, the relative contribution of the number and transcriptional programs in specific cell types towards the observed gene expression differences cannot be clearly delineated. Second, the relatively small sample sizes reduced the power for detection of subtle differences in expression. Also, due to small sample numbers, we had to rely on cross-validation methods for calculation of prediction errors instead of testing candidate predictors on new samples. The possibility of over-fitting cannot, therefore, be entirely ruled out.

Conclusions

Gene expression profiling in whole blood demonstrated significant differences in transcript levels that were capable of separating obese and lean phenotypes in multivariate analysis. Gene-set enrichment analysis further identified differences in biological pathways relating to cell survival, protein synthesis and energy harvest between the obese and lean groups. A subset of genes responsible for pathway enrichment also acted as efficient predictors of phenotype (obese or lean) when their expression signatures were used as inputs to Naive Bayes or logistic regression based classifiers. Together, our study is the first to investigate the information content in whole blood in relation to obesity. Our findings demonstrate that the investigation of gene expression profiles from whole blood can inform and illustrate the biological processes related to regulation of body mass. Additionally, the ability of pathway-related gene expression to predict class membership suggests the feasibility of a similar approach to identify blood-based robust predictors of weight loss success in response to dietary and surgical interventions.

Methods

Study Subjects

Twenty consecutive obese subjects enrolled in the Ottawa Hospital Weight Management Program at the Ottawa Hospital, Ottawa, with a body mass index (BMI) of 30-50 kg/m2, were recruited for study. All subjects were of Northern European White genetic ancestry. Patients were excluded on the basis of medical conditions possibly affecting whole blood gene expression, including out of normal range thyroid indices (TSH, free T3) at week 1 or week 13, diabetes mellitus treated with insulin or oral hypoglycemic agents, cigarette smoking, congestive heart failure, obstructive sleep apnea, active malignancy. Patients treated with weight-altering medications including tricyclic antidepressants, paroxetine, mirtazepine, lithium, valproate, gabapentin and typical and atypical antipsychotics, fluoxetine in doses greater than 20mg, bupropion, topiramate, systemic glucocorticoids and weight management drugs were also excluded. Blood samples were collected at baseline prior to initiation of weight loss therapy. Twenty lean subjects from the same genetic ancestry (Northern European White), with a BMI ≤ the 10th percentile for age and sex and no prior history of having had a BMI> 25th percentile for more than a 2 year consecutive period, were recruited from the Ottawa community. Lean subjects were excluded if they had any medical conditions affecting weight gain such as hyperthyroidism, anorexia nervosa, bulimia, major depression, or malabsorption syndromes. BMI for obese and lean subjects was categorized according to the population percentiles for age and sex using the Canadian Heart Health Survey data for subjects over the age of 18 years (data on file; Health Canada). The study protocol was approved by the Human Research Ethics Committees of the Ottawa Hospital and the University of Ottawa Heart Institute and informed consent was obtained from all participants prior to their enrolling into the program.

Sample preparation for transcriptome analysis

2.5 ml of fasting whole blood was drawn from study subjects by standard venipuncture and directly transferred to PAXgene blood RNA tubes (Qiagen, Santa Clara, CA). PAXgene tubes were processed at designated times after phlebotomy by the PAXgene protocol. Isolation of total RNA was accomplished according to the manufacturer's instructions. Prior to further processing, RNA quality was ascertained by electropherograms on the Agilent 2100 Bioanalyzer. Extracted RNA from all samples was stored -70°C until processed for microarray hybridizations.

Microarray hybridization and data analysis

Hybridization of 100 nanograms of labeled cRNA from each sample was carried out on Affymetrix GeneChip® Human Genome U133 Plus 2.0 Arrays according to the manufacturer's instructions. Microarray data was deposited in the Gene Expression Omnibus data repository (accession number GSE18897). Gene expression signals were generated from hybridized and scanned Affymetrix arrays by the GC-RMA algorithm [54]. Probesets with a normalized average expression level of less than 50 units in all of the tested groups were eliminated from further analysis. Significance of differential gene expression was ascertained via the signal-to-noise algorithm from the GenePattern Comparative Marker Selection module [22], employing a permutation-based t-test and false discovery rate (FDR) control. The Signal-to-Noise feature selection method is a variation of the more commonly used t-test statistic and looks at the difference of the means in each of the classes scaled by the sum of the standard deviations: Sx = (μ0-μ1)/(σ0 + σ1) where μ0 is the mean of class 0 and σ0 is the standard deviation of class 0. The Signal-to-Noise statistic penalizes genes that have higher variance in each class more than those genes that have a high variance in one class and a low variance in another.

Pathway analysis

Bioinformatic pathway analysis was conducted with the Gene Set Enrichment Analysis (GSEA) software package [27, 55]. GSEA is a computational method to detect statistically significant, concordant differences in a priori defined gene sets (pathways) between two biological states. GSEA accomplishes this task by calculating a weighted Kolmogorov-Smirnov statistic, adjusted for gene-set size (known as the Normalized Enrichment Score, NES) for each gene-set, based on the over-representation of members of a gene-set towards the top or bottom of a list of genes ranked by the strength of their correlation (positive or negative) to one of the two phenotypes. The statistical significance of NES score is estimated by a permutation test based on random shuffling of the phenotype or tag (gene) labels. GSEA addresses the problem of multiple testing (testing hundreds of gene-sets simultaneously) by calculating a false-discovery rate and a family-wise error rate on the ES p-values.

Quantitative real time polymerase chain reaction (RT-PCR)

Whole blood was collected in PAXgene™ blood tubes (Qiagen, Santa Clara, CA) and total RNA was extracted using the PAXgene™ blood kit. All RNA was treated with DNase I to remove genomic DNA contamination. The RNA was converted to cDNA in a 96-well microtiter plate on an ABI PRISM 7700 Sequence Detector System (Applied Biosystems, Foster City, CA) using the Applied Biosystems High Capacity cDNA archive kit. Gene expression was conducted on the Applied Biosystems 7900 using TaqMan® RT-PCR technology. A global median absolute deviation (MAD) was computed from the gene expression values by taking the median deviation for each set of technical replicates, using either the Ct values or log2 calculated abundances. Outliers were defined as having more than five times the global MAD. Following technical and biological outlier identification the data was normalized using reference housekeeper genes. The mean Ct value of all reference genes across all samples ("global mean Ct") was subtracted from the mean Ct value of all reference genes within each sample ("sample reference mean") to determine a normalization factor for each sample. The normalization factor for a given sample was then subtracted from its Ct value resulting in a normalized Ct. All Ct values were then converted to log2 abundances.

Class Prediction from gene expression

Class prediction (obese or lean) from gene expression data was carried out through the WEKA Explorer and WEKA Experimenter applications. First, 183 genes belonging to the 3 obese-upregulated pathways (ribosome, apoptosis and oxidative phosphorylation) were used to identify a subset of maximally informative features (genes) for classifier testing while removing irrelevant or redundant features that could negatively impact algorithm performance. Feature selection was accomplished by two independent 'filtering-based' algorithms (Information Gain and Cfs Subset Evaluator) and using 10-fold cross validation for each method [56, 57]. We did not use 'wrapper-based' feature-selection because we wanted the selected features to be independent of classification algorithms [58]. Both procedures resulted in a list of genes that were then ranked based on their importance in each feature selection method. From these ranked lists, we selected a total of 11 genes that were ranked within the top 20 genes in both lists. Gene expression signals for these 11 genes were then used as input in 4 different classifiers (Naïve Bayes, Logistic Regression, Random Forests and ZeroR) representing 4 different algorithmic approaches (Bayesian, regression, decision trees and rule-based, respectively) which were independently tested for predictive performance (Additional File 13) [59, 60]. Classifer-specific parameters were kept at the defaults provided in WEKA Experimenter. Each classifier used 66% of the samples for training (from a total of 34 obese plus lean subjects) and 33% for testing (chosen at random for each round) for a total of 100 iterations. For each classifier, the true positive rate, true negative rate, false positive rate, and false negative rates were calculated (average plus standard deviation over 100 iterations) and the values used to compare individual classifiers for their predictive performance.

References

  1. Pausova Z, Gossard F, Gaudet D, Tremblay J, Kotchen TA, Cowley AW, Hamet P: Heritability estimates of obesity measures in siblings with and without hypertension. Hypertension. 2001, 38 (1): 41-47.

    Article  CAS  PubMed  Google Scholar 

  2. Wardle J, Carnell S, Haworth CM, Plomin R: Evidence for a strong genetic influence on childhood adiposity despite the force of the obesogenic environment. Am J Clin Nutr. 2008, 87 (2): 398-404.

    CAS  PubMed  Google Scholar 

  3. Rankinen T, Zuberi A, Chagnon YC, Weisnagel SJ, Argyropoulos G, Walts B, Perusse L, Bouchard C: The human obesity gene map: the 2005 update. Obesity (Silver Spring). 2006, 14 (4): 529-644. 10.1038/oby.2006.71.

    Article  Google Scholar 

  4. Willer CJ, Speliotes EK, Loos RJ, Li S, Lindgren CM, Heid IM, Berndt SI, Elliott AL, Jackson AU, Lamina C, et al: Six new loci associated with body mass index highlight a neuronal influence on body weight regulation. Nat Genet. 2009, 41 (1): 25-34. 10.1038/ng.287.

    Article  CAS  PubMed  Google Scholar 

  5. Thorleifsson G, Walters GB, Gudbjartsson DF, Steinthorsdottir V, Sulem P, Helgadottir A, Styrkarsdottir U, Gretarsdottir S, Thorlacius S, Jonsdottir I, et al: Genome-wide association yields new sequence variants at seven loci that associate with measures of obesity. Nat Genet. 2009, 41 (1): 18-24. 10.1038/ng.274.

    Article  CAS  PubMed  Google Scholar 

  6. Meyre D, Delplanque J, Chevre JC, Lecoeur C, Lobbens S, Gallina S, Durand E, Vatin V, Degraeve F, Proenca C, et al: Genome-wide association study for early-onset and morbid adult obesity identifies three new risk loci in European populations. Nat Genet. 2009, 41 (2): 157-159. 10.1038/ng.301.

    Article  CAS  PubMed  Google Scholar 

  7. Loos RJ, Lindgren CM, Li S, Wheeler E, Zhao JH, Prokopenko I, Inouye M, Freathy RM, Attwood AP, Beckmann JS, et al: Common variants near MC4R are associated with fat mass, weight and risk of obesity. Nat Genet. 2008, 40 (6): 768-775. 10.1038/ng.140.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Scuteri A, Sanna S, Wei-Min C, Uda M, Albai G, Strait J, Najjar S, Nagaraja R, Orrð M, Usala G, et al: Genome-Wide Association Scan Shows Genetic Variants in the FTO Gene Are Associated with Obesity-Related Traits. PLoS Genetics. 2007, Public Library of Science, 3: e115-1210. 10.1371/journal.pgen.0030115.

    Google Scholar 

  9. Hinney A, Nguyen TT, Scherag A, Friedel S, Bronner G, Muller TD, Grallert H, Illig T, Wichmann HE, Rief W, et al: Genome wide association (GWA) study for early onset extreme obesity supports the role of fat mass and obesity associated gene (FTO) variants. PLoS ONE. 2007, 2 (12): e1361-10.1371/journal.pone.0001361.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Frayling TM, Timpson NJ, Weedon MN, Zeggini E, Freathy RM, Lindgren CM, Perry JRB, Elliott KS, Lango H, Rayner NW, et al: A Common Variant in the FTO Gene Is Associated with Body Mass Index and Predisposes to Childhood and Adult Obesity. Science. 2007, 316 (5826): 889-894. 10.1126/science.1141634.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Scherag A, Dina C, Hinney A, Vatin V, Scherag S, Vogel CI, Muller TD, Grallert H, Wichmann HE, Balkau B, et al: Two new Loci for body-weight regulation identified in a joint analysis of genome-wide association studies for early-onset extreme obesity in French and german study groups. PLoS Genet. 6 (4): e1000916-10.1371/journal.pgen.1000916.

  12. Glatt SJ, Everall IP, Kremen WS, Corbeil J, Šášik R, Khanlou N, Han M, Liew C-C, Tsuang MT: Comparative gene expression analysis of blood and brain provides concurrent validation of SELENBP1 up-regulation in schizophrenia. Proceedings of the National Academy of Sciences of the United States of America. 2005, 102 (43): 15533-15538. 10.1073/pnas.0507666102.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Vernon SD, Unger ER, Dimulescu IM, Rajeevan M, Reeves WC: Utility of the blood for gene expression profiling and biomarker discovery in chronic fatigue syndrome. Dis Markers. 2002, 18 (4): 193-199.

    Article  CAS  PubMed  Google Scholar 

  14. Batliwalla FM, Baechler EC, Xiao X, Li W, Balasubramanian S, Khalili H, Damle A, Ortmann WA, Perrone A, Kantor AB, et al: Peripheral blood gene expression profiling in rheumatoid arthritis. Genes Immun. 2005, 6 (5): 388-397. 10.1038/sj.gene.6364209.

    Article  CAS  PubMed  Google Scholar 

  15. Batliwalla FM, Li W, Ritchlin CT, Xiao X, Brenner M, Laragione T, Shao T, Durham R, Kemshetti S, Schwarz E, et al: Microarray analyses of peripheral blood cells identifies unique gene expression signature in psoriatic arthritis. Mol Med. 2005, 11 (1-12): 21-29. 10.2119/2006-00003.Gulko.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. DePrimo S, Wong L, Khatry D, Nicholas S, Manning W, Smolich B, O'Farrell A-M, Cherrington J: Expression profiling of blood samples from an SU5416 Phase III metastatic colorectal cancer clinical trial: a novel strategy for biomarker identification. BMC Cancer. 2003, 3 (1): 3-10.1186/1471-2407-3-3.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Tsuang MT, Nossova N, Yager T, Tsuang M-M, Guo S-C, Shyu KG, Glatt SJ, Liew CC: Assessing the validity of blood-based gene expression profiles for the classification of schizophrenia and bipolar disorder: A preliminary report. American Journal of Medical Genetics Part B: Neuropsychiatric Genetics. 2005, 133B (1): 1-5. 10.1002/ajmg.b.30161.

    Article  Google Scholar 

  18. Tusher V, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci USA. 2001, 98: 5116-5121. 10.1073/pnas.091062498.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Amundson SA, Grace MB, McLeland CB, Epperly MW, Yeager A, Zhan Q, Greenberger JS, Fornace AJ: Human In vivo Radiation-Induced Biomarkers: Gene Expression Changes in Radiotherapy Patients. Cancer Res. 2004, 64 (18): 6368-6371. 10.1158/0008-5472.CAN-04-1883.

    Article  CAS  PubMed  Google Scholar 

  20. Forrest MS, Lan Q, Hubbard AE, Zhang L, Vermeulen R, Zhao X, Li G, Wu YY, Shen M, Yin S, et al: Discovery of novel biomarkers by microarray analysis of peripheral blood mononuclear cell gene expression in benzene-exposed workers. Environ Health Perspect. 2005, 113 (6): 801-807. 10.1289/ehp.7635.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Miller GE, Chen E, Fok AK, Walker H, Lim A, Nicholls EF, Cole S, Kobor MS: Low early-life social class leaves a biological residue manifested by decreased glucocorticoid and increased proinflammatory signaling. Proceedings of the National Academy of Sciences. 2009, 106 (34): 14716-14721. 10.1073/pnas.0902971106.

    Article  CAS  Google Scholar 

  22. Gould J, Getz G, Monti S, Reich M, Mesirov JP: Comparative gene marker selection suite. Bioinformatics. 2006, 22 (15): 1924-1925. 10.1093/bioinformatics/btl196.

    Article  CAS  PubMed  Google Scholar 

  23. Kalofoutis A, Jullien G, Koutselinis A, Miras C: Erythrocyte 2,3-diphosphoglycerate as related to diabetes and obesity. Clin Chem. 1975, 21 (10): 1414-1415.

    CAS  PubMed  Google Scholar 

  24. Snyder LR, Born S, Lechner AJ: Blood oxygen affinity in high-and low-altitude populations of the deer mouse. Respir Physiol. 1982, 48 (1): 89-105. 10.1016/0034-5687(82)90052-4.

    Article  CAS  PubMed  Google Scholar 

  25. Rao GM, Morghom LO: Erythrocyte count and hemoglobin levels in diabetic women. Eur J Appl Physiol Occup Physiol. 1984, 52 (3): 272-275. 10.1007/BF01015208.

    Article  CAS  PubMed  Google Scholar 

  26. Wysocki M, Krotkiewski M, Braide M, Bagge U: Hemorheological disturbances, metabolic parameters and blood pressure in different types of obesity. Atherosclerosis. 1991, 88 (1): 21-28. 10.1016/0021-9150(91)90253-Y.

    Article  CAS  PubMed  Google Scholar 

  27. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, et al: Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005, 102 (43): 15545-15550. 10.1073/pnas.0506580102.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Ogata H, Goto S, Sato K, Fujibuchi W, Bono H, Kanehisa M: KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 1999, 27 (1): 29-34. 10.1093/nar/27.1.29.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Mootha VK, Lindgren CM, Eriksson KF, Subramanian A, Sihag S, Lehar J, Puigserver P, Carlsson E, Ridderstrale M, Laurila E, et al: PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat Genet. 2003, 34 (3): 267-273. 10.1038/ng1180.

    Article  CAS  PubMed  Google Scholar 

  30. Schaefer CF, Anthony K, Krupa S, Buchoff J, Day M, Hannay T, Buetow KH: PID: the Pathway Interaction Database. Nucleic Acids Res. 2009, D674-679. 10.1093/nar/gkn653. 37 Database

  31. Technote A: Signature genes for different blood subtypes. Affymetrix Technote. [http://media.affymetrix.com/support/technical/technotes/bloodappendix_technote.pdf]

  32. Whitney AR, Diehn M, Popper SJ, Alizadeh AA, Boldrick JC, Relman DA, Brown PO: Individuality and variation in gene expression patterns in human blood. Proc Natl Acad Sci USA. 2003, 100 (4): 1896-1901. 10.1073/pnas.252784499.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Solinas G, Vilcu C, Neels JG, Bandyopadhyay GK, Luo JL, Naugler W, Grivennikov S, Wynshaw-Boris A, Scadeng M, Olefsky JM, et al: JNK1 in hematopoietically derived cells contributes to diet-induced inflammation and insulin resistance without affecting obesity. Cell Metab. 2007, 6 (5): 386-397. 10.1016/j.cmet.2007.09.011.

    Article  CAS  PubMed  Google Scholar 

  34. Charriere G, Cousin B, Arnaud E, Andre M, Bacou F, Penicaud L, Casteilla L: Preadipocyte conversion to macrophage. Evidence of plasticity. J Biol Chem. 2003, 278 (11): 9850-9855. 10.1074/jbc.M210811200.

    Article  CAS  PubMed  Google Scholar 

  35. Hoppe S, Bierhoff H, Cado I, Weber A, Tiebe M, Grummt I, Voit R: AMP-activated protein kinase adapts rRNA synthesis to cellular energy supply. Proc Natl Acad Sci USA. 2009, 106 (42): 17781-17786. 10.1073/pnas.0909873106.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Barouch LA, Gao D, Chen L, Miller KL, Xu W, Phan AC, Kittleson MM, Minhas KM, Berkowitz DE, Wei C, et al: Cardiac Myocyte Apoptosis Is Associated With Increased DNA Damage and Decreased Survival in Murine Models of Obesity. Circ Res. 2006, 98 (1): 119-124. 10.1161/01.RES.0000199348.10580.1d.

    Article  CAS  PubMed  Google Scholar 

  37. Lupi R, Dotta F, Marselli L, Del Guerra S, Masini M, Santangelo C, Patané G, Boggi U, Piro S, Anello M, et al: Prolonged Exposure to Free Fatty Acids Has Cytostatic and Pro-Apoptotic Effects on Human Pancreatic Islets. Diabetes. 2002, 51 (5): 1437-1442. 10.2337/diabetes.51.5.1437.

    Article  CAS  PubMed  Google Scholar 

  38. Zhang HH, Kumar S, Barnett AH, Eggo MC: Dexamethasone Inhibits Tumor Necrosis Factor-{{alpha}}-Induced Apoptosis and Interleukin-1{beta} Release in Human Subcutaneous Adipocytes and Preadipocytes. J Clin Endocrinol Metab. 2001, 86 (6): 2817-2825. 10.1210/jc.86.6.2817.

    CAS  PubMed  Google Scholar 

  39. Dandona P, Aljada A, Bandyopadhyay A: Inflammation: the link between insulin resistance, obesity and diabetes. Trends in Immunology. 2004, 25 (1): 4-7. 10.1016/j.it.2003.10.013.

    Article  CAS  PubMed  Google Scholar 

  40. Xu H, Barnes GT, Yang Q, Tan G, Yang D, Chou CJ, Sole J, Nichols A, Ross JS, Tartaglia LA, et al: Chronic inflammation in fat plays a crucial role in the development of obesity-related insulin resistance. J Clin Invest. 2003, 112 (12): 1821-1830.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Kelley DE, He J, Menshikova EV, Ritov VB: Dysfunction of mitochondria in human skeletal muscle in type 2 diabetes. Diabetes. 2002, 51 (10): 2944-2950. 10.2337/diabetes.51.10.2944.

    Article  CAS  PubMed  Google Scholar 

  42. Patti ME, Butte AJ, Crunkhorn S, Cusi K, Berria R, Kashyap S, Miyazaki Y, Kohane I, Costello M, Saccone R, et al: Coordinated reduction of genes of oxidative metabolism in humans with insulin resistance and diabetes: Potential role of PGC1 and NRF1. Proc Natl Acad Sci USA. 2003, 100 (14): 8466-8471. 10.1073/pnas.1032913100.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Takamura T, Misu H, Matsuzawa-Nagata N, Sakurai M, Ota T, Shimizu A, Kurita S, Takeshita Y, Ando H, Honda M, et al: Obesity upregulates genes involved in oxidative phosphorylation in livers of diabetic patients. Obesity (Silver Spring). 2008, 16 (12): 2601-2609. 10.1038/oby.2008.419.

    Article  CAS  Google Scholar 

  44. Wallace DC: A MITOCHONDRIAL PARADIGM OF METABOLIC AND DEGENERATIVE DISEASES, AGING, AND CANCER: A Dawn for Evolutionary Medicine. Annual Review of Genetics. 2005, 39 (1): 359-407. 10.1146/annurev.genet.39.110304.095751.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Harper ME, Dent R, Monemdjou S, Bezaire V, Van Wyck L, Wells G, Kavaslar GN, Gauthier A, Tesson F, McPherson R: Decreased mitochondrial proton leak and reduced expression of uncoupling protein 3 in skeletal muscle of obese diet-resistant women. Diabetes. 2002, 51 (8): 2459-2466. 10.2337/diabetes.51.8.2459.

    Article  CAS  PubMed  Google Scholar 

  46. Abeel T, Helleputte T, Van de Peer Y, Dupont P, Saeys Y: Robust biomarker identification for cancer diagnosis with ensemble feature selection methods. Bioinformatics. 26 (3): 392-398. 10.1093/bioinformatics/btp630.

  47. Golub T, Slonim D, Tamayo P, Huard C, Gaasenbeek M, Mesirov J, Coller H, Loh M, Downing J, Caligiuri M, et al: Molecular classification of cancer: class discovery and class prediction by gene expression monitoring. Science. 1999, 286: 531-537. 10.1126/science.286.5439.531.

    Article  CAS  PubMed  Google Scholar 

  48. Dhanasekaran SM, Barrette TR, Ghosh D, Shah R, Varambally S, Kurachi K, Pienta KJ, Rubin MA, Chinnaiyan AM: Delineation of prognostic biomarkers in prostate cancer. Nature. 2001, 412 (6849): 822-826. 10.1038/35090585.

    Article  CAS  PubMed  Google Scholar 

  49. Pahl A, Brune K: Gene expression changes in blood after phlebotomy: implications for gene expression profiling. Blood. 2002, 100 (3): 1094-1095. 10.1182/blood-2002-03-0813.

    Article  CAS  PubMed  Google Scholar 

  50. Min J, Barrett A, Watts T, Pettersson F, Lockstone H, Lindgren C, Taylor J, Allen M, Zondervan K, McCarthy M: Variability of gene expression profiles in human blood and lymphoblastoid cell lines. BMC Genomics. 11 (1): 96-10.1186/1471-2164-11-96.

  51. Raghavachari N, Xu X, Munson PJ, Gladwin MT: Characterization of Whole Blood Gene Expression Profiles as a Sequel to Globin mRNA Reduction in Patients with Sickle Cell Disease. PLoS ONE. 2009, 4 (8): e6484-10.1371/journal.pone.0006484.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Wright C, Bergstrom D, Dai H, Marton M, Morris M, Tokiwa G, Wang Y, Fare T: Characterization of globin RNA interference in gene expression profiling of whole-blood samples. Clin Chem. 2008, 54 (2): 396-405. 10.1373/clinchem.2007.093419.

    Article  CAS  PubMed  Google Scholar 

  53. Liu J, Walter E, Stenger D, Thach D: Effects of globin mRNA reduction methods on gene expression profiles from whole blood. J Mol Diagn. 2006, 8 (5): 551-558. 10.2353/jmoldx.2006.060021.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Irizarry RA, Wu Z, Jaffee HA: Comparison of Affymetrix GeneChip expression measures. Bioinformatics. 2006, 22 (7): 789-794. 10.1093/bioinformatics/btk046.

    Article  CAS  PubMed  Google Scholar 

  55. Subramanian A, Kuehn H, Gould J, Tamayo P, Mesirov JP: GSEA-P: a desktop application for Gene Set Enrichment Analysis. Bioinformatics. 2007, 23 (23): 3251-3253. 10.1093/bioinformatics/btm369.

    Article  CAS  PubMed  Google Scholar 

  56. Indra Devi M, Rajaram R, Selvakuberan K: Generating best features for web page classification. Webology. 2008, 5 (1): Article 52.

    Google Scholar 

  57. Cover T, Thomas J: Elements of Information Theory. 1991, New York: Wiley, full_text.

    Book  Google Scholar 

  58. Williams N, Zander S, Armitage G: A preliminary performance comparison of five machine learning algorithms for practical IP traffic flow classification. ACM SIGCOMM Computer Communication Review. 2006, 36 (5): 7-15. 10.1145/1163593.1163596.

    Article  Google Scholar 

  59. Breiman L: Random Forests. Machine Learning. 2001, 455 (1): 5-32. 10.1023/A:1010933404324.

    Article  Google Scholar 

  60. Del Rio A, Gasteiger J: Simple method for the prediction of the separation of racemates with high-performance liquid chromatography on Whelk-O1 chiral stationary phase. Journal of Chromatography A. 2008, 1185 (1): 49-58. 10.1016/j.chroma.2008.01.034.

    Article  CAS  PubMed  Google Scholar 

Pre-publication history

Download references

Acknowledgements

This work was conducted with a grant support from GlaxoSmithKline. Part of the study was supported by NIH grants NHLBI-5R25HL059868-10 and NIDDK-1R21DK088319-01 (Ghosh) and a grant from the Heart & Stroke Foundation of Ontario (NA-5413; McPherson, Dent and Harper).

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Sujoy Ghosh or Ruth McPherson.

Additional information

Competing interests

The authors SG, SG(Gorman), and JS are former or current employees of GlaxoSmithKline and have equity in the company. The other authors (RD, RM, MEH) have no disclosures.

Authors' contributions

SG carried out the experimental design, data analysis, interpretation and drafted the manuscript. RD provided phenotype information and samples for transcriptome analysis and edited the manuscript. MEH participated in experimental design and interpretation of microarray data and manuscript editing. SG(Gorman) performed the Taqman analysis. JS coordinated the sample management, RNA isolation and microarray hybridizations. RM had overall oversight of the study and helped prepare the final version of the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

12920_2010_191_MOESM1_ESM.DOC

Additional file 1: Comparison of expression of adipocyte and muscle specific genes in whole blood samples utilized in the current study. Data for adipocyte-specific and muscle-specific gene expression was obtained from microarray data available on 79 different tissues from the Genomics Institute of Novartis Research Foundation http://www.gnf.org. Relative expression in whole blood was also obtained from the same source. The average expression and standard deviation in adipocyte-specific and muscle-specific gene expression observed in whole blood samples used in the current study (DOC 49 KB)

12920_2010_191_MOESM2_ESM.DOC

Additional file 2: PCA model output from multivariate analysis on obese and lean subjects based on whole blood gene expression signals. Analysis of performance of the PCA model separating obese from lean subjects based on blood gene expression signals. (DOC 34 KB)

12920_2010_191_MOESM3_ESM.DOC

Additional file 3: Top 100 differentially expressed genes in whole blood from obese and lean subjects. Using the GenePattern algorithm http://www.broadinstitute.org/cancer/software/genepattern/ a list of the top 50 upregulated and top 50 downregulated genes in obese and lean samples was generated and plotted on a heat-map for visualization. Higher expression levels are indicated in red and lower expression levels are indicated in blue. Genes (rows) are indicated by their Affymetrix probeset identifiers and samples (columns) are indicated by their obese or lean categories. (DOC 78 KB)

12920_2010_191_MOESM4_ESM.XLS

Additional file 4: Differential gene expression analysis between obese and lean subjects from blood transcriptome data. Identification of top 200 differentially expressed genes between Obese and Lean subjects using the Comparative Marker Selection module in GenePattern. Results for top 100 upregulated and top 100 downregulated genes (Obese vs. Lean) are shown. (XLS 64 KB)

12920_2010_191_MOESM5_ESM.DOC

Additional file 5: Description and value ranges of the parameters used in gene-set enrichment analysis (GSEA) in the present study. For detailed explanation of parameters and acceptable value ranges, please see additional documentation at http://www.broadinstitute.org/gsea/doc/GSEAUserGuideFrame.html. (DOC 24 KB)

12920_2010_191_MOESM6_ESM.XLS

Additional file 6: Custom pathway database used in the GSEA studies. Column 1 indicates the pathway name, column 2 indicates the source for the pathway information. Subsequent columns represent the gene symbols for the genes constituting the pathway. (XLS 369 KB)

12920_2010_191_MOESM7_ESM.XLS

Additional file 7: List of pathways determined to be upregulated in the Obese subjects compared to the Lean subjects and vice versa. These results were obtained by querying the pathways in the KEGG database (> = 10 and < = 200 gene members) and using the GSEA algorithm. (XLS 61 KB)

12920_2010_191_MOESM8_ESM.XLS

Additional file 8: List of gene members of the 3 pathways identified by GSEA as upregulated in obese subjects. List of genes in the oxidative phosphorylation, ribosome and apoptosis signaling pathways that were expressed in blood from obese and lean subjects. The column 'Core Enrichment' describes whether a gene contributed significantly towards the enrichment of the respective pathway in the gene set enrichment analysis. (XLS 42 KB)

12920_2010_191_MOESM9_ESM.XLS

Additional file 9: List of pathways determined to be upregulated in the Obese subjects compared to the Lean subjects and vice versa. These results were obtained by querying the pathways in a user-defined custom database (> = 10 and < = 200 gene members) and using the GSEA algorithm. (XLS 29 KB)

12920_2010_191_MOESM10_ESM.XLS

Additional file 10: Expression patterns of genes reported to be enriched in erythrocyte/reticulocyte fraction among lean and obese subjects. The average log2 expression signal in obese and lean subjects, differential expression (log ratio) and statistical significance of the differences in expression in the two groups are indicated. (XLS 26 KB)

12920_2010_191_MOESM11_ESM.DOC

Additional file 11: Overexpression of erythrocyte/reticulocyte enriched genes in obese blood samples. A scatter plot of log average expression of genes (x-axis) versus the differences of log expression between the obese and lean cohorts (y-axis) was created (also known as a MA plot). Each gene is indicated by a gray dot with the exception of genes reported to be enriched in erythrocytes/reticulocytes (compared to other blood cell types) which are shown as black pluses. A value of 0 on the y-axis signifies no differences in gene expression between the lean and obese cohorts. (DOC 85 KB)

12920_2010_191_MOESM12_ESM.PDF

Additional file 12: Gene set enrichment analysis between obese and lean subjects considering male only or female only cohorts. Pathways are ranked in descending order of their enrichment for each comparison. Results with males only are shown first, followed by the results with females only. (PDF 74 KB)

12920_2010_191_MOESM13_ESM.PDF

Additional file 13: Gene set enrichment analysis between obese and lean subjects after scaling of gene expression data by expression levels of erythrocyte membrane protein band 2 (EMPB2) and hemoglobin D (HBD) genes respectively. Pathways are ranked in descending order of enrichment. Top part refers to results obtained after scaling with EMPB2; bottom part shows the results following scaling with HBD. (PDF 76 KB)

12920_2010_191_MOESM14_ESM.DOC

Additional file 14: Description of classifiers and classifier performance evaluators used in the study. Brief descriptions of the Naïve Bayes, Logistic Regression, Random Forests and ZeroR classifiers along with the feature selection algorithms used (Information Gain and Cfs Subset evaluator) are given. Mathematical formulas for true and false positive and negative rates (classifier evaluation metrics) are also provided. (DOC 34 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Ghosh, S., Dent, R., Harper, ME. et al. Gene expression profiling in whole blood identifies distinct biological pathways associated with obesity. BMC Med Genomics 3, 56 (2010). https://doi.org/10.1186/1755-8794-3-56

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1755-8794-3-56

Keywords