Coexpression Network Analysis in Abdominal and Gluteal Adipose Tissue Reveals Regulatory Genetic Loci for Metabolic Syndrome and Related Phenotypes

Metabolic Syndrome (MetS) is highly prevalent and has considerable public health impact, but its underlying genetic factors remain elusive. To identify gene networks involved in MetS, we conducted whole-genome expression and genotype profiling on abdominal (ABD) and gluteal (GLU) adipose tissue, and whole blood (WB), from 29 MetS cases and 44 controls. Co-expression network analysis for each tissue independently identified nine, six, and zero MetS–associated modules of coexpressed genes in ABD, GLU, and WB, respectively. Of 8,992 probesets expressed in ABD or GLU, 685 (7.6%) were expressed in ABD and 51 (0.6%) in GLU only. Differential eigengene network analysis of 8,256 shared probesets detected 22 shared modules with high preservation across adipose depots (DABD-GLU = 0.89), seven of which were associated with MetS (FDR P<0.01). The strongest associated module, significantly enriched for immune response–related processes, contained 94/620 (15%) genes with inter-depot differences. In an independent cohort of 145/141 twins with ABD and WB longitudinal expression data, median variability in ABD due to familiality was greater for MetS–associated versus un-associated modules (ABD: 0.48 versus 0.18, P = 0.08; GLU: 0.54 versus 0.20, P = 7.8×10−4). Cis-eQTL analysis of probesets associated with MetS (FDR P<0.01) and/or inter-depot differences (FDR P<0.01) provided evidence for 32 eQTLs. Corresponding eSNPs were tested for association with MetS–related phenotypes in two GWAS of >100,000 individuals; rs10282458, affecting expression of RARRES2 (encoding chemerin), was associated with body mass index (BMI) (P = 6.0×10−4); and rs2395185, affecting inter-depot differences of HLA-DRB1 expression, was associated with high-density lipoprotein (P = 8.7×10−4) and BMI–adjusted waist-to-hip ratio (P = 2.4×10−4). Since many genes and their interactions influence complex traits such as MetS, integrated analysis of genotypes and coexpression networks across multiple tissues relevant to clinical traits is an efficient strategy to identify novel associations

Summary of <i>cis</i> eQTL results for selected expression probes significantly associated with MetS and showing differential ABD-GLU expression in MolOBB network and single-gene association analyses (FDR P<0.01).

<p>Independent association between the <i>cis</i> eQTL SNPs and metabolic traits was assessed in two GWA datasets.</p><p>DE = differential expression;</p><p>*genomic control adjusted pvalues from GIANT consortium;</p><p>**genomic control adjusted pvalues from <a href="http://www.plosgenetics.org/article/info:doi/10.1371/journal.pgen.1002505#pgen.1002505-Teslovich1" target="_blank">[4]</a>.</p

Sources of gene expression variation in different tissues.

<p>A) 626 probesets associated with MetS in MolOBB ABD B) 205 probesets associated with MetS in MolOBB GLU. Variances are decomposed in MolTWIN ABD and WB.</p

Modules for which eigengenes were significantly correlated (FDR p<0.01) with MetS in ABD (N = 9) and GLU (N = 6).

<p>FDR corrected pvalues for the associations with MetS and six quantitative metabolic traits are shown.</p><p>*DE = differentially expressed between adipose depots;</p><p>**Diast bp = diastolic blood pressure;</p><p>***Syst bp = systolic blood pressure.</p

MetS-associated modules found in the different fat depots.

<p>Visualization of the ABD brown (A), GLU darkgreen (B), consensus yellow in ABD (C) and consensus yellow in GLU (D) modules, respectively. For each module the top 150 pairwise correlations (intramodular connectivities) are shown. Genes with the top 10 highest ranked module membership are displayed in larger circles.</p

Biological Processes GO terms were significantly enriched (FDR P<0.01) in 15 modules associated with MetS in ABD and GLU.

<p>*P value = Fisher Exact Test;</p><p>**FE = Fold Enrichment.</p

Characteristics of the participants included in ABD gene expression analyses (with both nonmissing expression and phenotype data).

<p>Values are means ± standard deviation for each quantitative variable.</p><p>*Individuals (N = 7) with treatment for lipid abnormalities or hypertension were assigned as having reduced HDL or raised blood pressure. T2D is defined as fasting blood glucose >7.0 mmol/l or antiglycemic treatment).</p

Gene expression changes related to MetS in ABD and GLU.

<p>A,B. Heatmap of scaled expression values of 893 and 335 differentially expressed genes between MetS cases and controls in ABD (A) and GLU (B) samples. The dendogram depicts hierarchical clustering of the differentially expressed genes. The bottom bars show black boxes for MetS, the presence of the MetS components (reduced HDL, raised TG, raised fasting glucose, raised blood pressure) and gender (females). C. The top 10 genes differentially expressed in ABD and the top 10 genes differentially expressed in GLU. The horizontal bars displays −log₁₀ Pvalues for differential expression between MetS cases and controls in ABD and GLU and between depots.</p