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

Integrative network analysis reveals different pathophysiological mechanisms of insulin resistance among Caucasians and African Americans

Abstract

Background

African Americans (AA) have more pronounced insulin resistance and higher insulin secretion than European Americans (Caucasians or CA) when matched for age, gender, and body mass index (BMI). We hypothesize that physiological differences (including insulin sensitivity [SI]) between CAs and AAs can be explained by co-regulated gene networks in tissues involved in glucose homeostasis.

Methods

We performed integrative gene network analyses of transcriptomic data in subcutaneous adipose tissue of 99 CA and 37 AA subjects metabolically characterized as non-diabetic, with a range of SI and BMI values.

Results

Transcripts negatively correlated with SI in only the CA or AA subjects were enriched for inflammatory response genes and integrin-signaling genes, respectively. A sub-network (module) with TYROBP as a hub enriched for genes involved in inflammatory response (corrected p = 1.7E-26) was negatively correlated with SI (r = −0.426, p = 4.95E-04) in CA subjects. SI was positively correlated with transcript modules enriched for mitochondrial metabolism in both groups. Several SI-associated co-expressed modules were enriched for genes differentially expressed between groups. Two modules involved in immune response to viral infections and function of adherens junction, are significantly correlated with SI only in CAs. Five modules involved in drug/intracellular transport and oxidoreductase activity, among other activities, are correlated with SI only in AAs. Furthermore, we identified driver genes of these race-specific SI-associated modules.

Conclusions

SI-associated transcriptional networks that were deranged predominantly in one ethnic group may explain the distinctive physiological features of glucose homeostasis among AA subjects.

Peer Review reports

Background

Ethnic differences exist in prevalence of common metabolic disorders, including type 2 diabetes (T2D) and metabolic syndrome [1]. In the United States, the recently estimated age-adjusted prevalence of T2D in adults is 7.6% among non-Hispanic Caucasians (European Americans or CA), and 14.9% among non-Hispanic African Americans (AA) [2]. Clinical (visceral adiposity and waist circumference) and biochemical (markers of inflammation) factors only partly explain these ethnic disparities [3-5]. The physiology of glucose homeostasis in African Americans is distinctive, with more pronounced insulin resistance, higher insulin secretion, and reduced hepatic insulin clearance compared to CAs [6]. These ethnic differences are significant even after accounting for age, gender, and obesity. Differences in evolutionary history have distinctly shaped the genetic (and epigenetic) architecture of AA and CA populations and may explain the ethnic differences in the stabilization points (canalization) of insulin sensitivity and other glucose homeostasis phenotypes. A recent study suggested that this ethnic stratification could be implicated in the different natural courses of T2D onset and may lead to differences in T2D prevalence [7]. Identifying molecular mechanisms that can explain the observed ethnic differences in insulin sensitivity will be helpful in developing suitable prevention, surveillance, and treatment strategies for the pathogenesis of T2D. However, the genetic cause or the physiological mechanisms for greater insulin resistance among AAs are unknown [8].

Transcript expression in human tissues is correlated with insulin sensitivity and linked to dysregulation of genes involved in the pathophysiology of obesity, insulin resistance, and T2D [9]. Lymphoblastoid cell line transcript profiles of different ethnic groups from the HapMap project also identified significant ethnic differences in expression of transcripts [10-13]. Comparisons of gene expression in subcutaneous adipose and skeletal muscle from a small cohort of CA and AA subjects by our laboratory and others indicated ethnic differences in transcript expression; some of these transcripts are associated with insulin sensitivity [14,15]. However, no analyses have identified transcriptional mechanisms associated with insulin sensitivity and related metabolic traits that predominate in one ethnic group compared to another.

Expression of transcripts involved in the same biological function tend to be co-regulated by similar factors (genetic or environmental) and can be identified as distinct network modules, where genes within a module are more highly interconnected (correlated) with each other than genes in other modules [16]. Thus, we hypothesize that distinctive physiological differences, including differences in insulin sensitivity (SI) between CAs and AAs, can be explained at least partially by gene expression network modules in tissues involved in glucose homeostasis. These modules may show characteristic patterns of insulin resistance-associated derangements predominantly in one ethnic group compared to others.

In this study, we systematically analyzed genome-wide expression data from subcutaneous adipose tissue of 99 CA and 37 AA non-diabetic subjects. A pressing task in analyzing omics data like those generated by this study is to identify and visualize a global landscape of interactomes that contribute to clinical endpoints such as disease onset and progression. Weighted gene coexpression network analysis (WGCNA) has emerged as a way to solve this problem through identification of gene modules comprised of highly interconnected genes over a gene-gene interaction heatmap [17]. WGCNA has been widely used to identify pathways and gene targets for a variety of common human diseases such as cancer [18,19], atherosclerosis [20,21], Alzheimer’s disease [16], obesity and diabetes [22,23]. A key contribution of this work is the integration of differential expression analysis, gene-trait correlation analysis, and gene network analysis to identify both the coexpressed gene modules that contribute to insulin sensitivity and the key causal regulators. Utilizing this integrative network approach, we identified transcript co-expression modules correlated with SI and other related metabolic phenotypes in each ethnic group. Several modules enriched for important biological pathway genes showed differential connectivity between ethnic groups. Our study successfully delineated transcriptional networks associated with insulin sensitivity and other glucose homeostasis traits that are deranged predominantly in one ethnic group, and may explain differences in insulin sensitivity between CAs and AAs.

Methods

Ethics statement

All participants provided written informed consent under protocols originally approved (IRB #53028) by the University of Arkansas for Medical Sciences (UAMS). Our current study was approved by the Institutional Review Board at Wake Forest School of Medicine (IRB#00011083).

Study subjects

We used subcutaneous adipose tissue samples from 136 non-diabetic individuals. These subjects were recruited previously; detailed methods for recruitment, physical examination, physiological experiments, and obtaining biopsies have been published [24]. In brief, Caucasian (European-American) or African-American men and women who were generally in good health, between 19 and 60 years of age, and had a BMI between 19 and 45 kg/m2 were recruited. All participants had a screening visit, at which time height, weight, and waist and hip circumference were measured; body fat was determined by dual x-ray absorptiometry (DXA) scan (Hologic QDR-4500); fasting blood samples for lipid measurements were taken; and a standard 75-g oral glucose tolerance test (OGTT) was done, with measurement of glucose and insulin at baseline and 30- minute intervals for 2 hours. OGTT data were analyzed by homeostatic model assessment [25]. Based on their OGTT results, individuals who did not have diabetes were selected for this study.

At the second visit, adipose biopsies were obtained from subjects under fasting conditions. Insulin modified (0.03 U/kg) frequently sampled intravenous glucose tolerance tests (FSIGT) were performed as described previously [14,24]. We used the MINMOD Millennium program to analyze FSIGT data to determine insulin sensitivity (SI) and acute insulin response (AIRG) [26]. Biopsies were obtained from abdominal subcutaneous fat near the umbilicus under local anesthesia (Lidocaine) using a Bergstrom needle. Each biopsy is made by opening the Bergstrom needle and applying suction using a 120 ml syringe. Samples were rinsed immediately in sterile normal saline, quick-frozen in liquid nitrogen, and stored at −80°C for further use. Demographic characteristics of our cohort are shown in Table 1.

Table 1 Demographics of study population

Laboratory measurements

Insulin was measured by the UAMS Clinical Research Center core laboratory using an immuno-chemiluminometric assay (Invitron Limited, Monmouth, UK). Plasma glucose was measured by glucose oxidase methods at LabCorp, Inc. (Burlington, NC). Plasma triglyceride, total cholesterol, and HDL cholesterol concentrations were measured directly (at LabCorp) by enzymatic colorimetric methods using an automated Hitachi Cobas C system and assay kits from Roche. LDL cholesterol and VLDL cholesterol concentrations were calculated indirectly by the Friedewald equation.

RNA extraction

Total RNA was isolated from whole adipose tissue using the RNAeasy Lipid Tissue Mini kit (QIAGEN Inc-USA, Valencia, CA). The quantity and quality of the isolated total RNA samples were determined by ultraviolet spectrophotometry (Nanodrop, Thermo Scientific, Pittsburgh, PA) and electrophoresis (Experion nucleic acid analyzer, BioRad Laboratories, Inc., Hercules, CA), respectively. High-quality RNA with RIN (RNA integrity number) >8 was used for genome-wide transcriptome analysis.

Microarray studies

Genome-wide transcriptome analysis and initial array processing were done at the Center for Human Genomics Core Laboratory (Wake Forest School of Medicine) using HumanHT-12 v4 Expression BeadChip (Illumina, San Diego, CA) whole genome gene expression arrays according to the vendor-recommended standard protocol [24]. In brief, biotin-labeled cRNA was prepared from 500 ng of total RNA using the Illumina TotalPrep RNA amplification kit, and 750 ng of biotinylated cRNA was hybridized to the chip in a hybridization oven at 58°C for 18 hours. Chips were processed through a series of posthybridization washes, blocked in blocking buffer, and moved to streptavidin-Cy3-containing detection buffer. Chips were then washed, dried, and scanned in the Illumina BeadArray Reader. All RNA samples were processed using the same lot of chips and reagents, and all experiments were performed within a month to avoid technical variability. Raw expression intensity was background subtracted and normalized by the average normalization algorithm as implemented in GenomeStudio Gene Expression Module v1.0 application software (Illumina). Normalized data were used for further analysis. Data for 17,434 probes annotated as an Entrez gene and with a detection p-value of ≤0.01 in 90% of the samples were used for final analysis. All expression data has been deposited in Gene Expression Omnibus (GEO; http://www.ncbi.nlm.nih.gov/projects/geo/) under accession number GSE65221.

Real-time qPCR validation

Total RNA (1 μg) was reverse transcribed using a Qiagen QuantiTect reverse transcription kit (QIAGEN Inc., Valencia, CA) according to the manufacturer’s protocol. Transcript-specific oligonucleotide primers for quantitative real-time PCR (qRT-PCR) were designed to capture splice variants captured by the probe sequences in the array [14,24]. Selected transcripts were measured by qRT-PCR using Power SYBR green chemistry (Applied Biosystems, Inc., Foster City, CA) and normalized to the expression of human ribosomal protein, large P0 (RPLP0 or 36B4) gene [24]. The standard curves were generated for absolute quantification using pooled cDNA from the samples assayed.

Statistical and bioinformatic data analysis

To identify differential expression of individual transcripts between CA and AA subjects, normalized gene-expression data from adipose tissue were first adjusted using linear regressions for age and gender (expression ~ age + gender). The residuals following regression were analyzed by Student’s T-test to identify transcripts with differential expression between CA and AA samples. Differentially expressed transcripts were defined as those with a ≥1.2 fold change and a p value ≤ 0.05. The residual expression values after covariate adjustments (for age and gender) were used for further testing, including network construction and gene-expression physiological phenotypic correlations. We used Spearman’s correlation coefficient to identify correlation of each transcript with SI. False discovery rates (FDR) were also calculated; a transcript was considered associated with SI if FDR ≤0.05.

The CA and AA data sets from all 17,434 expressed probes were independently processed through the weighted gene co-expression network analysis (WGCNA) [17,27]. The weighted network analysis began with a matrix of the Pearson correlations between all gene pairs, and then converted the correlation matrix into an unsigned adjacency matrix using a power function, so that the resulting adjacency matrix, i.e., the weighted coexpression network, is approximately scale-free. To explore the modular structures of the co-expression network, the adjacency matrix was further transformed into a topological overlap matrix (TOM) [17]. Because topological overlap between two genes reflects both their direct interaction and their indirect interactions through all other genes in the network, this approach helps create more cohesive and biologically more meaningful modules. To identify modules of highly coregulated genes, we used average linkage hierarchical clustering to group genes based on the topological overlap of their connectivity, followed by a dynamic cut-tree algorithm to dynamically cut clustering dendrogram branches into gene modules. To distinguish between modules, each module was assigned a unique color identifier. The gene with the highest intramodular connectivity was considered the “hub gene” (most connected gene).

To examine how each gene module was related to gluco-metabolic traits (e.g., SI) we first performed principal component analysis for each module, and then computed module-trait correlation between the first principal component (module eigengene) and each trait [16,28]. The significance (p value) and false discovery rate (FDR) of each correlation were also calculated. The FDR was estimated through random permutation of sample names of the trait data. A module was associated with a trait if the FDR was ≤0.05. Overrepresentation of canonical pathways and biological processes in modules was measured through Fisher’s exact test and corrected for numbers of modules and functional categories tested. Ingenuity Pathway Analysis (IPA) also was used to reveal enrichment of canonical pathways within selected modules.

To quantify differences in adipose tissue transcript network organization between CA and AA subjects, we employed a modular differential connectivity (MDC) metric [16]. In brief, MDC is the ratios of the connectivity of all gene pairs in a module from CA subjects, to that of the same gene pairs from AA subjects. MDC is a continuous measure ranging from 0 to infinity. MDC > 1 indicates gain of connectivity or enhanced co-regulation between genes, whereas MDC < 1 indicates loss of connectivity or reduced co-regulation between genes. The statistical significance of the MDC metrics was computed through the FDR procedure. The significance or FDR of the MDC statistic can be accessed by permuting the data underlying the two networks. We estimated FDR based both on shuffled samples (i.e., networks with nonrandom nodes but random connections) and shuffled gene labels (i.e., networks with random nodes but nonrandom connections), and then selected the larger value as the final FDR estimate.

To further identify key regulators (drivers) genes of the modules identified by WGCNA, we applied the key driver analysis [16] to the module-based unweighted coexpression networks derived from ARACNE (Algorithm for the Reconstruction of Accurate Cellular Networks) \( \overline{\mu}+\sigma \left(\mu \right) \) are nominated as key regulators [16]. This criterion identified genes with number of NHNNs significantly above the corresponding average value.

We also tested conservation of SI-associated modules between ethnic groups by analyzing overlap of module members and enrichment of modules for common or ethnic-specific SI-associated genes. Fisher’s exact test was used to quantitate significant overlap and enrichment. Modules in CAs that show significant overlap (FET_P ≤ 3.1E-06, i.e., Bonferroni corrected p value of 0.05) with a module in AAs and are enriched for SI-associated genes in both ethnicities were considered conserved, while modules enriched for genes associated with SI in only one ethnic group were considered ethnic-specific.

We used two different methods to control for multiple testing issues: (1) the highly conservative Bonferroni correction method, to constrain the study-wise significance level; and (2) an empirical false discovery rate (FDR) method, which constrains the overall rate of false positive events. For instance, for the enrichment of functional categories in modules, we corrected for the number of modules and Gene Ontology processes and other categories tested. For the enrichment of expression clinical trait correlation signatures in modules, we used a combination of FDR <0.05 and nominal P < 0.05.

Results

Adipose tissue transcripts show ethnic differences in expression

We identified significant differential expression (fold change 1.2 and p < 0.05) of 655 transcript probes between CA and AA subjects. These differentially expressed probes included 225 up-regulated and 193 down-regulated mRNA transcripts (Additional file 1: Table S1) of RefSeq genes with known function. Genes with >2-fold higher expression in AA subjects included crystallin-beta 2 (CRYBB2), splicing factor 1 (SF1), and C-type lectin domain family 4, member A (CLEC4A). Genes with >2 fold lower expression included diabetes-associated protein in insulin-sensitive tissue (DAPIT/USMG5), macrophage receptor with collagenous structure (MARCO), and dehydrogenase/reductase member 4 like 1 (DHRS4L1) (Figure 1A). We also validated differential expression of MARCO by qRT-PCR (Figure 1B) and detected a positive correlation of this gene with BMI in CA subjects (Figure 1C). Annotation of these differentially expressed genes using DAVID [30] indicated overrepresentation of these genes in important functional categories (Additional file 1: Table S2). Highly represented categories included genes involved in signal transduction and cell communication, endoplasmic reticulum, regulation of blood pressure, and the carboxylic acid biosynthetic process. In IPA analyses, enrichment of genes involved in glutathione-mediated detoxification and NRF2-mediated oxidative stress response pathways were among those differentially expressed in AA subjects.

Figure 1
figure 1

Genes showing strong differential expression in adipose tissue of African Americans (AA) compared to Caucasians (CA). A) Genome-wide expression microarray identified 2-fold higher expression of the CRYBB2, SF1, GSTT2B, and CLEA4A gene in AA subjects, while expression of DHRS4L1, MARCO (Macrophage receptor with collagenous structure), and USMG5 was lower; B) Technical validation by qRT-PCR confirmed lower expression of MARCO gene in adipose tissue of AAs compared to CAs. MARCO expression is shown after 36B4 normalization and ln transformation. In A and B the box represents the interquartile range, which contains 50% of the values. The whiskers are lines that extend the box to the highest and lowest values, excluding outliers. A line across the box indicates the median; and C)MARCO expression is positively correlated with BMI in 99 non-diabetic CA subjects.

We identified significant correlations of 1,625 transcripts with insulin sensitivity (SI) in both CA and AA subjects (812 and 813 transcripts showed positive and negative correlations, respectively with SI). Some transcripts showed significant correlation (FDR ≤ 0.05) with SI only in CA (n = 2,797 transcripts) or only in AA subjects (n = 1,983 transcripts) (Figure 2). These transcripts were enriched for important gene ontology categories and pathways (Additional file 1: Table S3 and Figure 2). Transcripts negatively correlated with SI only in CA subjects (n = 1,367 transcripts) were enriched for inflammatory response pathway genes (75 genes, corrected p-value = 6.90E-09). Transcripts negatively correlated with SI only in AA subjects (n = 932 transcripts) were enriched for integrin signaling pathway genes (24 genes, corrected p-value = 8.90E-06).

Figure 2
figure 2

Insulin sensitivity (S I )-associated genes common and distinct for ethnic groups are enriched for important functional category and biological pathways. Venn diagram showing overlap between Caucasian and African-American subjects for genes correlated (FDR ≤ 0.05) with SI in subcutaneous adipose tissue. See Additional file 1: Table S3 for further details.

Co-expression modules in adipose are associated with SI and related metabolic phenotypes in Caucasians

Network distance coupled with hierarchical clustering and dynamic tree cutting methods implemented in WGCNA identified 93 co-expression networks (modules) in adipose tissue of CA subjects (module sizes from 1,957 to 10 transcripts). Modules are denoted by colors (Figure 3A). Four modules in CA subjects included >1,000 probes, the largest being the “turquoise” module, with 1,957 transcript probes. Transcripts that could not be included in any module (3,777 of 17,434 probes) were classified as “gray”. Module membership of each transcript, simple adjacency matrix-based connectivity, and TOM-based connectivity measures of those transcripts within each module are shown in Additional file 1: Table S4.

Figure 3
figure 3

Weighted gene co-expression analysis of adipose tissue transcript in adipose tissue from Caucasian (CA) subjects identified gene networks (modules) associated with gluco-metabolic phenotypes. A) Topological Overlap Matrix (TOM) plot of the adipose tissue expression network connections in 99 CA subjects. Light color represents low topological overlap and progressively darker red color represents higher overlap. Each module is assigned by unique color; these are shown along the left side and the top, B) Correlation heatmap of 93 module eigengenes and 12 gluco-metabolic traits, C) A Circos plot showing module names (Ring 1), module color representation (Ring 2), correlations between module eigengenes and insulin sensitivity (SI) (Ring 3), with the bar height representing absolute value of correlation coefficient (range = [0.0012,0.53]) and the bar color for the correlation sign (red – positive and green – negative), module specificity to CA (high bar – CA-specific and low bar – not CA-specific) (Ring 4), significance of the enrichment for up-regulated genes in AAs with the bar height representing –log10(FET p value) (Ring 5) (range = [2.2e-39,1]), significance of the enrichment for down-regulated genes in AAs with the bar height representing –log10(FET p value) (Ring 6) (range = [5.8e-6,1]), and significance of the best enrichment for the modules in AA with the bar height representing –log10(FET p value) (Ring 7) (range = [0, 0.015]). The module names and color representations match those in Figure 3A. The modules are sorted based on the correlations between the module eigengenes and SI.

Correlation analysis of module eigengenes (principal components of the module) with gluco- and cardio-metabolic traits of CA subjects identified some strong module-trait correlations. As shown in Figure 3B, some modules were correlated with multiple phenotypes (also see Additional file 1: Table S5). The clustering of SI and triglyceride (TG) indicates a significant overlap of modules correlated with these two traits in CA subjects. Modules showing significant positive correlation with SI showed significant negative correlation with TG in CA subjects. Module eigengene of 27 modules showed significant correlation with SI (p ≤ 0.05 and FDR ≤ 0.05) (Figure 3C). Eight large modules (>25 gene members) showed strong correlation (positive or negative correlation ≥0.4) with either FSIGT-derived (SI and AIRg) [26] or OGTT-derived (HOMA-IR and Matsuda index) [25] glucose homeostasis traits (Table 2). The “yellow” module showed strong negative correlation (r = −0.426, P = 4.95E-04) with SI and was strongly enriched for genes involved in inflammatory response (corrected p = 1.7E-26). The TYROBP gene (TYRO protein tyrosine kinase binding protein) was the most connected hub gene of this co-expression network module based on TOM-based intramodular connectivity, consistent with the finding of the prominent role of TYROBP in other complex diseases [16,31].

Table 2 Adipose tissue co-expression modules associated with insulin sensitivity (S I ) or other glucose homeostasis traits in Caucasians

The “pink” and “light yellow” modules, which are both enriched for mitochondrion genes are positively correlated with SI with r = 0.411 (P = 3.53E-04) and 0.413 (P = 7.3E-06), respectively. Two nuclear-encoded mitochondrial genes, SDHB (Succinate dehydrogenase complex, subunit B) and DCI (Dodecenoyl-Coenzyme A delta isomerase or ECI1), were the hub genes for these two modules. The “light yellow” module included genes involved in carboxylic acid metabolism (corrected p = 5.7E-06), while the “pink” module was enriched for branched-chain amino acid (valine, leucine, and isoleucine) degradation genes (corrected p = 2.70E-14, see Additional file 1: Table S6).

Co-expression modules in adipose are associated with SI and related metabolic phenotypes in African Americans

WGCNA of adipose tissue transcripts of 37 African American subjects using the same parameters as described above for Caucasians (beta = 5, scale R^2 = 0.81, trunc. R^2 = 0.99, hiercutoff =0.99 and minModuleSize =10, Additional file 2: Figure S1) identified 173 co-expression modules (module size ranges from 1980 to 10 probes, Figure 4A, Additional file 1: Table S7). Three modules in AA subjects included >1,000 probes, the largest being the “turquoise” module, with 1980 transcript probes. In AA subjects, clustering of module-trait correlations indicated the strongest overlap between SI- and AIRg-associated modules (Figure 4B). SI was correlated with the eigengene of 34 modules (p ≤ 0.05 and FDR ≤ 0.05) (Figure 4C and Additional file 1: Table S8). Sixteen large modules showed strong correlations (positive or negative correlation ≥0.4) with either FSIGT- or OGTT-derived glucose homeostasis traits in AA subjects (Table 3). Those with the strongest negative correlations with SI were “midnight blue” (r = −0.58, p = 4.050E-04), “green” (r = −0.56, p = 0.001), and “thistle” (r = −0.52, p = 0.015). These modules were enriched for genes involved in function of nucleolus (p = 3.7E-07), protein refolding (p = 0.000032), and antioxidant activity (p = 0.0014), respectively (Additional file 1: Table S9). Annexin A1 (ANXA1 or Lipocortin I), a gene with phospholipase A2 inhibitory activity that belongs to a family of Ca(2+)-dependent phospholipid binding proteins, showed a strong negative correlation with insulin sensitivity (r = −0.66, p <0.00001) and the hub gene of the “midnight blue” module. The “black” module showed the strongest positive correlation with SI (r = 0.61, p = 0.001) in adipose tissue of AA subjects. This module showed significant overlap (p ≤ 3.1E-06) with “tan” (57 genes), “light yellow” (33 genes), and “pink” (41genes) modules identified in CA subjects (Figure 5). Similar to the “light yellow” module in CA subjects, the “black” module in AA subjects was enriched for genes involved in mitochondrial carboxylic acid metabolism (corrected p = 2.80E-13); the nuclear encoded mitochondrial gene DCI appeared as the hub gene for this module.

Figure 4
figure 4

Weighted gene co-expression analysis of adipose tissue transcript in adipose tissue identified distinct pattern of association between gene networks (modules) and gluco-metabolic phenotypes in African American subjects. A) TOM plot of the adipose tissue indicates distinct expression network connections in 37 AA subjects. Light color represents low topological overlap and progressively darker red color represents higher overlap. Each module is assigned by unique color; these are shown along the left side and the top, B) Correlation heatmap of 173 module eigengenes and 12 gluco-metabolic traits, C) A Circos plot showing module names (Ring 1), module color representation (Ring 2), correlations between module eigengenes and insulin sensitivity (SI) (Ring 3), with the bar height representing absolute value of correlation coefficient (range = [0.0033,0.61]) and the bar color for the correlation sign (red – positive and green – negative), module specificity to AA (high bar – AA-specific and low bar – not AA-specific) (Ring 4), significance of the enrichment for up-regulated genes in AAs with the bar height representing–log10(FET p value) (Ring 5) (range = [8.2e-19,1]), significance of the enrichment for down-regulated genes in AAs with the bar height representing –log10(FET p value) (Ring 6) (range = [1.8e-7,1]), and significance of the best enrichment for the modules in CA with the bar height representing–log10(FET p value) (Ring 7) (range = [0, 0.13]). The module names and color representations match those in Figure 4A. The modules are sorted based on the correlations between the module eigengenes and SI.

Table 3 Adipose tissue co-expression modules associated with insulin sensitivity (S I ) or other glucose homeostasis traits in African Americans
Figure 5
figure 5

Overlap of insulin sensitivity (S I ) associated modules in Caucasians and African Americans. Fisher’s exact test (FET) analysis was used to identify overlap among gene members of SI-associated modules in CA and AA subjects. SI-associated modules in CAs that showed significant overlap (FET_P value < 3.1E-06) with at least one SI-associated module in AAs, or vice versa, are shown. P ≥ 0.05, P ≤ 3.1E-06, P ≤ 1.0E-10 and P ≤ 1.0E-20 are marked in gray, pink, brown, and red, respectively. Enrichment of these modules for biological pathways or processes is shown.

To validate our findings, we compared our co-expression modules with the coexpression network modules from an independent human adipose gene expression data (called deCODE data) from 640 Caucasian subjects [23]. Twelve modules in the AA network and nineteen modules in the CA network significantly (Bonferroni corrected FET p-value < 0.05) overlap with the network modules from the deCODE data. Notably, 5 of the top 10 CA modules associated with SI identified in our study significantly overlap with the modules from the deCODE data, and the yellow (inflammatory response), pink (mitochondrion) and peru (steroid metabolism) modules are among the top five most conserved modules, with FET p values as 3.6 × 10−151 (5.6 fold enrichment), 6.6 × 10−58 (17.5 fold enrichment) and 1.4 × 10−28 (13 fold enrichment), respectively. On the other hand, 4 of the top 10 AA modules associated with SI significantly overlapped with the gene modules from the deCODE data, but to a lesser degree (Bonferroni corrected FET p-values are between 8.2 × 10−3 and 2.4 × 10−8). This comparison basically validates our findings about both the AA and CA subjects, although the results from the CA data are more profoundly replicated in another independent study in CA subjects. However, clinical and metabolic phenotype data are not publicly available for deCODE study subjects to validate our other findings, including correlation between modules and phenotypes.

Some metabolic trait-associated modules show ethnic differences

We tested the enrichment of modules for genes differentially expressed in AA compared to CA subjects. Nineteen co-expression modules in adipose of CA subjects were enriched for differentially expressed genes (enrichment p ≤0.05). Among these modules, five showed significant correlations with SI (Additional file 1: Table S10A). The “red” module in CA subjects included the largest number of genes (50 genes, p =3.17E-26) upregulated in African Americans, correlated with AIRG (r = 0.276, p = 0.0476), and was enriched for genes involved in cytoskeletal protein binding and adherens junction (enrichment p-value <3.0E-10). Similarly, 18 transcript modules in adipose of AA subjects were enriched for genes differentially expressed in African Americans compared to Caucasians, and five were significantly correlated with SI (Additional file 1: Table S10B). The “magenta” network module in AA adipose tissue was enriched for genes differentially up-regulated in AA compared to CA subjects (p = 0.0084, includes 11 up-regulated transcripts), showed a significant positive correlation with AIRG (r = 0.51; p = 0.006), and was enriched for genes involved in unsaturated fatty acid biosynthesis (corrected p = 0.0037) and glucose metabolism (corrected p = 0.007).

We used the modular differential connectivity (MDC) metric [16] to quantify differences in adipose tissue transcript network organization between CA and AA subjects. This analysis suggested significant differential connectivity (FDR < 0.2) in 59 CA adipose tissue-derived modules (Additional file 1: Table S11), almost all showing gain of connectivity (MDC > 1). Among these modules, 13 showed significant correlation with SI in CA subjects and were enriched for genes in functional categories important for insulin resistance (Table 4). Similar analyses of AA adipose tissue-derived networks identified MDC in 11 modules; 10 had high gain of connectivity (MDC> > 1) and 1 had loss of connectivity (MDC < 1). Only one of these modules was correlated with insulin sensitivity (Additional file 1: Table S12).

Table 4 Insulin sensitivity associated co-expression modules in adipose tissue of Caucasian subjects shows differential connectivity

We identified 16 and 15 conserved SI-associated network modules in CAs and AAs, respectively (FDR ≤ 0.05) (Table 5A and B, Figures 3C and 4C). These modules showed significant overlap (FET_p value ≤3.1E-06) of module gene members in both ethnicities, and were significantly enriched for genes correlated with SI in both ethnicities. We also identified five and nine ethnic-specific SI-associated network modules, respectively, in CAs and AAs (FDR ≤ 0.05) (Table 6). These modules were significantly enriched for genes correlated with SI in only one ethnicity (Figures 3C and 4C). A subset of modules (three in CA and five in AA) also showed no significant overlap (FET_p value > 3.1E-06) of gene members between ethnicities, which indicates the presence of highly ethnic-specific mechanisms associated with insulin sensitivity.

Table 5 Insulin sensitivity associated co-expression modules conserved in adipose tissue of both Caucasian and African American subjects
Table 6 Ethnicity specific Insulin sensitivity associated co-expression modules

Among these ethnicity-specific SI-associated modules, “bisque” is most significantly enriched for transcripts negatively correlated with SI in CAs (p =1.85E-16). This module, which has a strong differential connectivity (MDC = 11.16), is enriched for the interferon signaling pathway (in IPA analysis B-H P-value = 5.01E-20, Figure 6A, B). Several interferon-inducible genes involved in immune response to viral infection (including OAS2, ISG15, STAT1, and ADAR) were identified as “drivers” in this sub-network (Figure 6C) based on the ARACNE [29] and key driver analysis [16,31]. The pro-inflammatory T-cell-derived cytokine interferon gamma (IFNγ) promotes the macrophage M1 phenotype. IFNγ also modulates adipocyte metabolic functions. McGillicuddy et al. (2009) demonstrated that IFNγ may play an important role in T-cell modulation of diet-induced obesity, insulin resistance, and type 2 diabetes via activation of the adipocyte JAK-STAT pathway [32]. IFNγ attenuates insulin signaling, lipid storage, and differentiation in human adipocytes, an effect likely mediated via sustained STAT1 activation [33]. Indeed, 19 genes of this module were associated with metabolic and cardiovascular disorders based on the Genetic Association Database (http://geneticassociationdb.nih.gov/). Specifically, IRF7 and NLRC5 are associated with type 2 diabetes, and SP110 and ZBP1 are associated with waist circumference and body fat distribution, respectively.

Figure 6
figure 6

The “bisque” module is a Caucasian-specific insulin sensitivity (S I )-associated network module. A) Enrichment of the “bisque” module genes for canonical pathways based on IPA analysis, B) Genes in this module are involved in the interferon signaling pathway, and are negatively correlated with SI (shown in green color), and C) Network plots of the “bisque” module in adipose of CA subjects. The red nodes are “drivers”; the size of a node is proportional to the number of nodes in the node's 2-layer neighborhood.

Among the 9 AA-specific and SI-associated modules, the “light yellow” module is most enriched for the transcripts positively correlated with SI in AAs (p = 6.95E-09). The top driver gene of this sub-network (Figure 7A), YWHAE (14-3-3-epsilon protein), may play a role in the regulation of SI by modulating the interaction between the insulin receptor and IRS1. This module was enriched for oxidoreductase activity (p = 2.90E-04), and was marginally enriched for genes in ERK5 (p = 0.007) and vascular endothelial growth factor (VEGF) signaling pathways (p = 0.018) (Figure 7B) on IPA analysis. Both pathways include the YWHAE gene. The top interaction network based on the IPA database also includes YWHAE (14-3-3) (Figure 7C). We identified 35 genes from this module in the genetic association database, including genes associated with type 2 diabetes, obesity, and insulin resistance (i.e. PNPLA2, SOD2, ZNF716, VAMP3, SHC1, IL18, DLC1, CYFIP1, and ASPSCR1).

Figure 7
figure 7

The “light yellow” module is an African American-specific insulin sensitivity (S I )-associated network module. A) Network plots of the “light yellow” module in adipose of AA subjects. The red nodes are “drivers”; the size of a node is proportional to the number of nodes in the node's 2-layer neighborhood, B) Enrichment of module genes for canonical pathways in IPA analysis, C) The most significant network among the module members based on IPA interaction network analysis involves the YWHAE (14-3-3) gene. Genes positively or negatively correlated with SI in this interaction network are shown in red and green, respectively.

Although the “bisque” and “light yellow” modules significantly overlap with others in the AA and CA networks (Table 6), less than 26% of their members are in common with their counterparts. Therefore, they are considered as ethnicity-specific. ARACNE analysis identified seven Caucasian specific driver genes (ISG15, ADAR, EPSTI1, IRF7, HERC6, MX1 and IFI44L) for the “bisque” module. Similarly, among the 24 drivers of the “light yellow” module in AA, eight (YWHAE, LLPH, TNFRSF14, CCDC91, FLJ38717, LOC196752, COX19, and SLC44A4) are also the drivers of the same set of genes in CA. However, the remaining 16 drivers (CITED4, UBA7, SIRT2, LOC100132323, RFNG, LOC440345, KIAA1683, ZBTB40, C10orf6, PATE2, LOC158301, LOC389517, LOC641972, SNX26/ARHGAP33, LOC92017, and SFRS5/SRSF5) of this module appear to be AA-specific.

Discussion

The primary goal of our study was to measure transcriptional modulation related to SI and related metabolic traits in adipose tissue, and to examine how these might relate to molecular mechanisms of insulin resistance in African Americans compared to European Americans. Our analyses revealed distinct ethnic-specific, SI-associated configurations and biological functions of these co-expression modules. Our study provides critical insight into the transcriptional mechanisms and molecular processes underlying insulin sensitivity that could be implicated in T2D onset and may explain differences in T2D prevalence between CA and AA individuals.

We identified significant differential expression of 655 transcript probes between CA and AA subjects, which represent 418 genes with known function. A significant enrichment of these genes in important functional categories, including the glutathione-mediated detoxification and oxidative stress response pathways, was observed. Several isoforms of glutathione S-transferases (GSTT2, GSTT2B, GSTM3, and GSTM5) showed higher expression, while microsomal glutathione S-transferase 1 (MGST1) and superoxide dismutase 2 (SOD2) showed lower expression in AA subjects compared to CA subjects. The increased GST transcript expression we observed may indicate a stronger oxidative stress response capacity (in response to environmental factors) [34] in African Americans. However, further molecular genetic studies will be required to define its role in relation to insulin resistance and obesity.

Intriguingly, AA subjects show lower expression of a transcript splice form of USMG5 (DAPIT). DAPIT is a phylogenetically conserved ATP synthase (mitochondrial and lysosomal) that was down-regulated at the mRNA level in skeletal muscle of streptozotocin (STZ)-treated diabetic rats [35]. Its protein expression was up-regulated in epididymal adipose tissue of STZ-treated mice [36]. However, this gene was not associated with insulin sensitivity or BMI in our study. Compared to CA subjects, AA subjects also show significantly lower expression of the MARCO gene, which was positively correlated with BMI in CA subjects. The MARCO gene is a scavenger receptor expressed on macrophages and dendritic cells. Macrophages from MARCO-deficient mice exhibited lower interleukin-12 (IL-12) production in responses to stimulation with lipopolysaccharide and interferon-γ. MARCO −/− mice have severely impaired ability to clear bacteria from the lungs and increased mortality during bacterial infection, possibly due to impaired Th1 polarization [37,38]. Further studies will be required to reveal the species-specific and tissue-specific role of the DAPIT and MARCO genes.

We detected associations between several transcript modules, SI, and related metabolic phenotypes. Analysis to identify enrichment of these trait-associated models for functional categories of genes revealed molecular mechanisms that may cause insulin resistance in both ethnic groups, and indicated some pathophysiological mechanisms that are restricted to one ethnic group. The “black” module in adipose of AA subjects overlapped significantly with the “light yellow” module identified in CA subjects. The DCI (or ECI1) was the hub gene for both modules. This gene encodes a key mitochondrial enzyme involved in beta-oxidation of unsaturated fatty acids [39]. These modules also showed strong enrichment for genes involved in mitochondrial carboxylic acid and branched chain amino acid metabolism. Positive correlations of these modules enriched for mitochondrial function with insulin sensitivity, and negative correlations with BMI in both AA and CA subjects, indicate that a common molecular mechanism may be involved in obesity-associated insulin resistance.

Previous studies had indicated a role of inflammatory and immune response pathway genes in SI and obesity. However, these studies were restricted to Caucasian subjects. In line with these findings, we identified negative correlations between SI and the “yellow” module in adipose of CA subjects (r = −0.426, p = 4.95E-04) and a strong positive correlation of this module with BMI (r = 0.607, p = 4.03E-09). This module was strongly enriched for genes involved in inflammatory response (corrected p = 1.7E-26). IPA indicated enrichment of genes involved in Fcγ receptor-mediated phagocytosis in macrophages and monocytes (p = 2.51E-12). The TYROBP is the hub gene for this transcript network module. This gene encodes a transmembrane signaling polypeptide containing an immunoreceptor tyrosine-based activation motif in its cytoplasmic domain. TYROBP is involved in signaling events that mediate crosstalk between dendritic cells and natural killer cells. Recently, it was identified as a central regulator of a gene expression network in brain tissue involved in late-onset Alzheimer’s disease [16,40]. The “yellow” module in CA subjects showed the strongest overlap with the “turquoise” module of AA subjects (452 genes, p = 5.01E-169). This “turquoise” module showed enrichment for immune cell activation and the module eigengene showed significant correlation with SI (r = −0.315, FDR =2.84E-02) in AA subjects. We found enrichment of inflammatory response pathway genes (p-value = 6.90E-09) among genes that negatively correlated with SI only in CA subjects. Interestingly, a recent study in an animal model suggested that local pro-inflammatory response in adipose tissue is adaptive, required for proper adipose tissue remodeling and expansion to enable safe storage of excess nutrient and metabolic homeostasis [41].

Many adipose co-expression modules were enriched for genes differentially expressed between the two ethnicities or had differential connectivity (MDC) among members of the network module. Notably, the “bisque” module showed strong differential connectivity (MDC = 11.16) and was highly enriched for transcripts negatively correlated with SI only in CAs. Driver genes of this sub-network (OAS2, ISG15, STAT1, and ADAR) are involved in immune response to viral infections, and expression of these genes is inducible by interferon. A recent personal-omics analysis showed transcriptional modulation of genes in pathways involved in glucose regulation of insulin secretion occurred shortly after viral infection, and was predicted as involved in elevation of glucose and HbA1c in a Caucasian subject [42]. A recent in vitro analysis showed significant ADAR genotype (and haplotype)-dependent modulation of IFNγ-producing cells (or IFNγ-Elispot responses upon viral infection) in CAs but not in AAs [43]. We showed that ADAR is one of the 7 CA-specific driver genes in this module. IRF7, another CA-specific driver of this module, was modulated by lipopolysaccharide, interferon, and virus-responsive cis-regulatory genetic variants in monocytes and dendritic cells [44,45]. IRF7 modulates expression of several downstream genes in this pathway as a master trans-regulator. The “light yellow” module was enriched for transcripts correlated with SI only in AAs. The top driver gene of this sub-network, YWHAE (14-3-3-epsilon protein), may play a role in the regulation of SI by modulating the interaction between the insulin receptor and IRS1. However, the functions of many other genes in this sub-network are not clear.

Our study cohort included far more Caucasian subjects (N = 99) than African American subjects (N = 37). Unequal racial distribution in our cohort was influenced by the population demographics of Arkansas (recruitment site) and volunteer participation bias. However, to our knowledge, this is the largest single cohort of adipose tissue gene expression related to glucose homeostasis that includes phenotypically well-characterized subjects of both races. While sample size is an important factor for all analyses (e.g. differential expression, gene-trait correlation and network construction), the integrative network analysis was developed in this study to minimize the impact of a relatively small sample size, to obtain a more comprehensive understanding of differentially expressed genes and the network context under which individual genes operate. The highly significant overlap between the results from the two ethnic groups studied here and an independent but much larger study (the deCODE study) in Caucasian subjects suggests that the integrative network analysis approach was appropriate in balancing for the impact of the small sample size.

Conclusions

By comparing co-expression networks in adipose tissue of metabolically well-characterized CA and AA subjects, we found derangements in transcriptional networks of genes involved in mitochondrial energy metabolism pathways that lead to insulin resistance in both ethnicities. We also found evidence that other molecular mechanisms involved in modulating glucose homeostasis are restricted predominantly to subjects from one ethnic origin. Distinct genetic (and epigenetic) regulatory architecture, including differences in the frequency of expression regulatory alleles among populations, may determine the structure of co-expression networks observed in adipose tissues of these subjects [46,47]. Further studies will be required to identify how genetic and epigenetic factors determine the structure of co-expression networks in adipose tissue that modulate glucose homeostasis and related physiological traits.

References

  1. Bonds DE, Zaccaro DJ, Karter AJ, Selby JV, Saad M, Goff Jr DC. Ethnic and racial differences in diabetes care: The Insulin Resistance Atherosclerosis Study. Diabetes Care. 2003;26:1040–6.

    Article  PubMed  Google Scholar 

  2. Cowie CC, Rust KF, Byrd-Holt DD, Gregg EW, Ford ES, Geiss LS, et al. Prevalence of diabetes and high risk for diabetes using A1C criteria in the U.S. population in 1988–2006. Diabetes Care. 2010;33:562–8.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Alderete TL, Toledo-Corral CM, Goran MI. Metabolic basis of ethnic differences in diabetes risk in overweight and obese youth. Curr Diab Rep. 2014;14:455–0455.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Araneta MR, Barrett-Connor E. Ethnic differences in visceral adipose tissue and type 2 diabetes: Filipino, African-American, and white women. Obes Res. 2005;13:1458–65.

    Article  PubMed  Google Scholar 

  5. Hyatt TC, Phadke RP, Hunter GR, Bush NC, Munoz AJ, Gower BA. Insulin sensitivity in African-American and white women: association with inflammation. Obesity (Silver Spring). 2009;17:276–82.

    Article  CAS  Google Scholar 

  6. Haffner SM, D'Agostino R, Saad MF, Rewers M, Mykkanen L, Selby J, et al. Increased insulin resistance and insulin secretion in nondiabetic African-Americans and Hispanics compared with non-Hispanic whites. The Insulin Resistance Atherosclerosis Study. Diabetes. 1996;45:742–8.

    Article  CAS  PubMed  Google Scholar 

  7. Kodama K, Tojjar D, Yamada S, Toda K, Patel CJ, Butte AJ. Ethnic differences in the relationship between insulin sensitivity and insulin response: a systematic review and meta-analysis. Diabetes Care. 2013;36:1789–96.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. McCormack S, Grant SF. Genetics of obesity and type 2 diabetes in African Americans. J Obes. 2013;19:396416.

    Google Scholar 

  9. Keller MP, Attie AD. Physiological insights gained from gene expression analysis in obesity and diabetes. Annu Rev Nutr. 2010;30:341–64.

    Article  CAS  PubMed  Google Scholar 

  10. Spielman RS, Bastone LA, Burdick JT, Morley M, Ewens WJ, Cheung VG. Common genetic variants account for differences in gene expression among ethnic groups. Nat Genet. 2007;39:226–31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Storey JD, Madeoy J, Strout JL, Wurfel M, Ronald J, Akey JM. Gene-expression variation within and among human populations. Am J Hum Genet. 2007;80:502–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Stranger BE, Nica AC, Forrest MS, Dimas A, Bird CP, Beazley C, et al. Population genomics of human gene expression. Nat Genet. 2007;39:1217–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Zhang W, Duan S, Kistner EO, Bleibel WK, Huang RS, Clark TA, et al. Evaluation of genetic variation contributing to differences in gene expression between populations. Am J Hum Genet. 2008;82:631–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Elbein SC, Kern PA, Rasouli N, Yao-Borengasser A, Sharma NK, Das SK. Global gene expression profiles of subcutaneous adipose and muscle from glucose-tolerant, insulin-sensitive, and insulin-resistant individuals matched for BMI. Diabetes. 2011;60:1019–29.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Smith LM, Yao-Borengasser A, Starks T, Tripputi M, Kern PA, Rasouli N. Insulin resistance in African-American and Caucasian women: differences in lipotoxicity, adipokines, and gene expression in adipose tissue and muscle. J Clin Endocrinol Metab. 2010;95:4441–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Zhang B, Gaiteri C, Bodea LG, Wang Z, McElwee J, Podtelezhnikov AA, et al. Integrated systems approach identifies genetic nodes and networks in late-onset Alzheimer's disease. Cell. 2013;153:707–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Zhang B, Horvath S. A general framework for weighted gene co-expression network analysis. Stat Appl Genet Mol Biol. 2005, Volume 4, Issue 1, ISSN (Online) 1544-6115, doi:10.2202/1544-6115.1128, August 2005.

  18. Tran LM, Zhang B, Zhang Z, Zhang C, Xie T, Lamb JR, et al. Inferring causal genomic alterations in breast cancer using gene expression data. BMC Syst Biol. 2011;5:121–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Horvath S, Zhang B, Carlson M, Lu KV, Zhu S, Felciano RM, et al. Analysis of oncogenic signaling networks in glioblastoma identifies ASPM as a molecular target. Proc Natl Acad Sci U S A. 2006;103:17402–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Gargalovic PS, Imura M, Zhang B, Gharavi NM, Clark MJ, Pagnon J, et al. Identification of inflammatory gene modules based on variations of human endothelial cell responses to oxidized lipids. Proc Natl Acad Sci U S A. 2006;103:12741–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Huan T, Zhang B, Wang Z, Joehanes R, Zhu J, Johnson AD, et al. A systems biology framework identifies molecular underpinnings of coronary heart disease. Arterioscler Thromb Vasc Biol. 2013;33:1427–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Chen Y, Zhu J, Lum PY, Yang X, Pinto S, MacNeil DJ, et al. Variations in DNA elucidate molecular networks that cause disease. Nature. 2008;452:429–35.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. 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.

    Article  CAS  PubMed  Google Scholar 

  24. Das SK, Sharma NK, Hasstedt SJ, Mondal AK, Ma L, Langberg KA, et al. An integrative genomics approach identifies activation of thioredoxin/thioredoxin reductase-1-mediated oxidative stress defense pathway and inhibition of angiogenesis in obese nondiabetic human subjects. J Clin Endocrinol Metab. 2011;96:E1308–13.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Matsuda M, DeFronzo RA. Insulin sensitivity indices obtained from oral glucose tolerance testing: comparison with the euglycemic insulin clamp. Diabetes Care. 1999;22:1462–70.

    Article  CAS  PubMed  Google Scholar 

  26. Boston RC, Stefanovski D, Moate PJ, Sumner AE, Watanabe RM, Bergman RN. MINMOD Millennium: a computer program to calculate glucose effectiveness and insulin sensitivity from the frequently sampled intravenous glucose tolerance test. Diabetes Technol Ther. 2003;5:1003–15.

    Article  CAS  PubMed  Google Scholar 

  27. Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008;9:559.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Haas BE, Horvath S, Pietilainen KH, Cantor RM, Nikkola E, Weissglas-Volkov D, et al. Adipose co-expression networks across Finns and Mexicans identify novel triglyceride-associated genes. BMC Med Genomics. 2012;5:61–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Margolin AA, Nemenman I, Basso K, Wiggins C, Stolovitzky G, Dalla FR, et al. ARACNE: an algorithm for the reconstruction of gene regulatory networks in a mammalian cellular context. BMC Bioinformatics. 2006;7 Suppl 1:S7.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57.

    Article  CAS  Google Scholar 

  31. Wang IM, Zhang B, Yang X, Zhu J, Stepaniants S, Zhang C, et al. Systems analysis of eleven rodent disease models reveals an inflammatome signature and key drivers. Mol Syst Biol. 2012;8:594.

    Article  PubMed  PubMed Central  Google Scholar 

  32. McGillicuddy FC, Chiquoine EH, Hinkle CC, Kim RJ, Shah R, Roche HM, et al. Interferon gamma attenuates insulin signaling, lipid storage, and differentiation in human adipocytes via activation of the JAK/STAT pathway. J Biol Chem. 2009;284:31936–44.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Zhao GN, Jiang DS, Li H. Interferon regulatory factors: at the crossroads of immunity, metabolism, and disease. Biochim Biophys Acta. 2014;1852(2):365–78.

    Article  PubMed  Google Scholar 

  34. Board PG, Menon D. Glutathione transferases, regulators of cellular metabolism and physiology. Biochim Biophys Acta. 1830;2013:3267–88.

    Google Scholar 

  35. Paivarinne H, Kainulainen H. DAPIT, a novel protein down-regulated in insulin-sensitive tissues in streptozotocin-induced diabetes. Acta Diabetol. 2001;38:83–6.

    Article  CAS  PubMed  Google Scholar 

  36. Kontro H, Hulmi JJ, Rahkila P, Kainulainen H. Cellular and tissue expression of DAPIT, a phylogenetically conserved peptide. Eur J Histochem. 2012;56:e18.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Jozefowski S, Arredouani M, Sulahian T, Kobzik L. Disparate regulation and function of the class A scavenger receptors SR-AI/II and MARCO. J Immunol. 2005;175:8032–41.

    Article  CAS  PubMed  Google Scholar 

  38. Jozefowski S. The role of the class A scavenger receptors, SR-A and MARCO, in the immune system. Part 2. Contribution to recognition and phagocytosis of pathogens as well as induction of immune response. Postepy Hig Med Dosw (Online). 2012;66:120–31.

    Google Scholar 

  39. Janssen U, Fink T, Lichter P, Stoffel W. Human mitochondrial 3,2-trans-enoyl-CoA isomerase (DCI): gene structure and localization to chromosome 16p13.3. Genomics. 1994;23:223–8.

    Article  CAS  PubMed  Google Scholar 

  40. De Felice FG, Ferreira ST. Inflammation, defective insulin signaling, and mitochondrial dysfunction as common molecular denominators connecting type 2 diabetes to Alzheimer disease. Diabetes. 2014;63:2262–72.

    Article  PubMed  Google Scholar 

  41. Wernstedt AI, Tao C, Morley TS, Wang QA, Delgado-Lopez F, Wang ZV, et al. Adipocyte inflammation is essential for healthy adipose tissue expansion and remodeling. Cell Metab. 2014;20:103–18.

    Article  Google Scholar 

  42. Chen R, Mias GI, Li-Pook-Than J, Jiang L, Lam HY, Chen R, et al. Personal omics profiling reveals dynamic molecular and medical phenotypes. Cell. 2012;148:1293–307.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Haralambieva IH, Ovsyannikova IG, Umlauf BJ, Vierkant RA, Shane PV, Jacobson RM, et al. Genetic polymorphisms in host antiviral genes: associations with humoral and cellular immunity to measles vaccine. Vaccine. 2011;29:8988–97.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Fairfax BP, Humburg P, Makino S, Naranbhai V, Wong D, Lau E, et al. Innate immune activity conditions the effect of regulatory variants upon monocyte gene expression. Science. 2014;343:1246949.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Lee MN, Ye C, Villani AC, Raj T, Li W, Eisenhaure TM, et al. Common genetic variants modulate pathogen-sensing responses in human dendritic cells. Science. 2014;343:1246980.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Das SK, Sharma NK. Expression quantitative trait analyses to identify causal genetic variants for type 2 diabetes susceptibility. World J Diabetes. 2014;5:97–114.

    PubMed  PubMed Central  Google Scholar 

  47. Kang HP, Yang X, Chen R, Zhang B, Corona E, Schadt EE, et al. Integration of disease-specific single nucleotide polymorphisms, expression quantitative trait loci and coexpression networks reveal novel candidate genes for type 2 diabetes. Diabetologia. 2012;55:2205–13.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

This work was partly supported by the grants R01 DK039311 (to S.D.) from the National Institutes of Health (NIH)/ National Institute of Diabetes and Kidney Disease (NIDDK), R01 DK090111 (to S.D.) from NIH/ NIDDK, R01AG046170 (to B.Z.) from the NIH/National Institute on Aging (NIA), R01CA163772 (to B.Z.) from NIH/National Cancer Institute (NCI), and U01AI111598-01 (to B.Z.) from NIH/National Institute of Allergy and Infectious Diseases (NIAID). Microarray based gene expression data generation was supported by WFSM intramural faculty development fund. We thank Lijun Ma (Wake Forest School of Medicine) for technical assistance in using GenomeStudio and Yongzhong Zhao (Icahn School of Medicine at Mount Sinai) for submitting genome-wide expression data in GEO. We also thank Karen Klein (Biomedical Research Services Administration, Wake Forest University Health Sciences) and Ethel Kouba (Internal Medicine, Endocrinology) for critical reading and editing of our manuscript. We thank the late Steven C. Elbein, MD for recruitment and metabolic characterization of the study subjects analyzed in this study.

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Swapan Kumar Das or Bin Zhang.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

SKD designed the project, researched and analyzed data, and wrote, reviewed, and edited manuscript. NKS contributed data and reviewed the manuscript. BZ researched and analyzed data, reviewed and edited the manuscript. All authors read and approved the final manuscript.

Additional files

Additional file 1: Table S1.

Genes showing differential expression in adipose tissue of African Americans compared to Caucasians. Table S2. Enrichment of functional category and biological pathway among genes showing differential expression in adipose tissue of African Americans compared to Caucasians. Table S3. Enrichment of functional category and biological pathway among genes correlated with insulin sensitivity (SI) in adipose tissue. Table S4. Weighted gene co-expression network analysis (WGCNA) identified 93 co-expression networks (modules) in adipose tissue of Caucasian subjects. Table S5. Correlation of module eigengenes with gluco- and cardio- metabolic phenotypes in Caucasian subjects. Table S6. Co-expression modules in adipose tissue of Caucasians are enriched for genes in functional categories and biological pathways. Table S7. Weighted gene co-expression network analysis (WGCNA) identified 173 co-expression networks (modules) in adipose tissue of African American subjects. Table S8. Correlation of module eigengenes with gluco- and cardio- metabolic phenotypes in African American subjects. Table S9. Co-expression modules in adipose tissue of African Americans are enriched for genes in functional categories and biological pathways. Table S10. A) Adipose tissue derived transcript modules in CA are enriched for genes differentially expressed in African Americans. B) Adipose tissue derived transcript modules in AA are enriched for genes differentially expressed in African Americans. Table S11. Co-expression modules in adipose tissue of Caucasians show differential connectivity. Table S12. Co-expression modules in adipose tissue of African Americans show differential connectivity.

Additional file 2: Figure S1.

Scale free topology plot of weighted adipose tissue co-expression network in (A) CA subjects (N = 99) and (B) AA subjects (N = 37) constructed with the power adjacency function power (s, beta = 5). Network satisfies a scale free topology approximately (black linear regression line R2 = 0.81. But a better fit is provided by an exponentially truncated power law (red line, R2 = 0.99).

Rights and permissions

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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Das, S.K., Sharma, N.K. & Zhang, B. Integrative network analysis reveals different pathophysiological mechanisms of insulin resistance among Caucasians and African Americans. BMC Med Genomics 8, 4 (2015). https://doi.org/10.1186/s12920-015-0078-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12920-015-0078-0

Keywords