- Research article
- Open Access
- Open Peer Review
A comprehensive analysis of adiponectin QTLs using SNP association, SNP cis-effects on peripheral blood gene expression and gene expression correlation identified novel metabolic syndrome (MetS) genes with potential role in carcinogenesis and systemic inflammation
BMC Medical Genomicsvolume 6, Article number: 14 (2013)
Metabolic syndrome (MetS) is an aberration associated with increased risk for cancer and inflammation. Adiponectin, an adipocyte-produced abundant protein hormone, has countering effect on the diabetogenic and atherogenic components of MetS. Plasma levels of adiponectin are negatively correlated with onset of cancer and cancer patient mortality. We previously performed microsatellite linkage analyses using adiponectin as a surrogate marker and revealed two QTLs on chr5 (5p14) and chr14 (14q13).
Using individuals from 85 extended families that contributed to the linkage and who were measured for 42 clinical and biologic MetS phenotypes, we tested QTL-based SNP associations, peripheral white blood cell (PWBC) gene expression, and the effects of cis-acting SNPs on gene expression to discover genomic elements that could affect the pathophysiology and complications of MetS.
Adiponectin levels were found to be highly intercorrelated phenotypically with the majority of MetS traits. QTL-specific haplotype-tagging SNPs associated with MetS phenotypes were annotated to 14 genes whose function could influence MetS biology as well as oncogenesis or inflammation. These were mechanistically categorized into four groups: cell-cell adhesion and mobility, signal transduction, transcription and protein sorting. Four genes were highly prioritized: cadherin 18 (CDH18), myosin X (MYO10), anchor protein 6 of AMPK (AKAP6), and neuronal PAS domain protein 3 (NPAS3). PWBC expression was detectable only for the following genes with multi-organ or with multi-function properties: NPAS3, MARCH6, MYO10 and FBXL7. Strong evidence of cis-effects on the expression of MYO10 in PWBC was found with SNPs clustered near the gene’s transcription start site. MYO10 expression in PWBC was marginally correlated with body composition (p= 0.065) and adipokine levels in the periphery (p = 0.064). Variants of genes AKAP6, NPAS3, MARCH6 and FBXL7 have been previously reported to be associated with insulin resistance, inflammatory markers or adiposity studies using genome-wide approaches whereas associations of CDH18 and MYO10 with MetS traits have not been reported before.
Adiponectin QTLs-based SNP association and mRNA expression identified genes that could mediate the association between MetS and cancer or inflammation.
MetS is a cluster of phenotypes characterized by preferential deposition of fat in the abdominal/visceral region, insulin resistance, dyslipidemia, increased blood pressure and increased plasma levels of adipokines/cytokines . Its complications include increased blood pressure/hypertension, glucose intolerance and type 2 diabetes. MetS has been associated with several forms of cancer including breast, endometrial, cervical, ovarian, esophageal, colon, rectal, pancreatic, hepatic, biliary and kidney cancers . Overweight/obesity by itself accounts for 14% of cancer deaths in men and 20% in women . Childhood cancer survivors are also known to present features of the MetS in later life .
Adiponectin is an adipocyte-produced kinin encoded by the gene apM1 (adipose most abundant transcript 1) at chr3q27 . Adiponectin circulates in two forms, the low molecular weight (LMW) dimers, and the high molecular weight (HMW) oligomers. Tissue specificity of adiponectin n-mers is determined by their relative affinity for receptors (AdipoR1 and R2). AdipoR1 has a higher affinity for HMW. It is ubiquitously expressed but highly enriched in skeletal muscle. AdipoR2 can bind both forms and is mainly expressed in the liver . Adiponectin has been shown to enhance insulin sensitivity and to exert anti-diabetogenic and anti-atherogenic functions as well as anti-inflammatory and anti-angiogenic activities [6–8]. Its plasma levels are inversely correlated with increased risk for obesity-related malignancies .
To explore possible genetic mechanisms linking MetS pathways with severe adverse outcomes such as development of cancer, we previously conducted a genome-wide linkage study using microsatellite markers at approximately 10 centiMorgan (cM) intervals in 1,100 individuals from 170 nuclear families (out of 85 extended families) of predominantly Northern European ancestry. This study identified two highly significant quantitative trait loci (QTLs) on chr5p14 and chr14q13 . In the present study, we tested for genes harbored within these QTLs that could account for the association between MetS and carcinogenesis or systemic inflammation. For this purpose, we evaluated 1,137 individuals from the 85 extended families that contributed to the original linkage signals and who were phenotyped for 42 clinical and/or biological MetS traits (described in Results). We then performed QTL-focused SNP association and PWBC target gene expression analyses to identify those genes predisposing for this association.
Subjects and phenotypes
The present study cohort included 1,137 individuals from 85 extended families. Details of their recruitment and ascertainment traits procedures have been previously described [11, 12]. Recruitment was initiated via an obese proband (BMI ≥ 30) with minimal availability of one obese and one never-obese (BMI ≤ 27) siblings and at least one, preferably both, parents. The clinical phenotypic components, which included weight, height, BMI, waist circumference (WC), hip circumference (HC), waist to hip ratio (WHR), fasting glucose (FG), fasting insulin (FI), insulin to glucose ratio (IGR), homeostasis model assessment (HOMA), plasma triglycerides (TG), total cholesterol (TC), LDL-cholesterol (LDL-c) and calculated LDL-cholesterol levels (cal. LDL-c), HDL-cholesterol (HDL-c), systolic and diastolic blood pressure (sBP and dBP) and pulse, as described . The biological traits included total fat mass both in kilogram and percentage (Fatkg and Fatpct) and total fat free mass both in kilograms and as percentage of total body mass (Leankg and Leanpct) by DXA (Dual-emission X-ray absorptiometry) ; total abdominal, visceral and subcutaneous fat size (TAF, VF and SubQF) by computed tomography (CT) scans of the fourth lumbar spine ; respiratory quotient (RQ) and resting energy expenditure (REE) measured in resting subjects using a Deltatrac Indirect Calorimeter (Sensor Medics, VIASYS Healthcare, Conshohocken, PA) after a 10 hr fast; insulin/glucose responsiveness indices: insulin sensitivity (SI), glucose effectiveness (SG), acute insulin response to glucose (AIRG) and disposition index (DI) by Minimal Model ; lipids/lipoprotein sizing [HDL median diameter (HMED), LDL-cholesterol median diameter (LMEDn), LDL-cholesterol dominant peak diameter (LDLppd) and apoB-containing non-HDL median diameter (BMED) which includes VLDL, ILDL, LPα and LDL] measured by polyacrylamide gradient gel electrophoresis ; circulating levels of adiponectin, leptin by a double antibody, equilibrium RIA (Millipore Corporation, Billerica, MA); and TNF-alpha, interleukin-1beta (IL-1β) and interleukin-6 (IL-6)] measured as described . All study procedures for adults and children were approved by the Institutional Review Boards of the Medical College of Wisconsin and Children’s Hospital of Wisconsin, respectively. Informed consent was obtained from the participating subjects or guardians of participating minors.
SNP genotyping and data cleaning
Genomic DNA was extracted and prepared from whole blood using commercial kits (Puregene, Minneapolis, MN). Genome-wide SNP genotyping was performed using Affymetrix Genome-Wide Human SNP 6.0 arrays and SNP calls were generated by Genotype Console 3.2. Individuals with fewer than 95% of all available markers called were excluded. 869,222 autosomal SNPs were prepared by Preswalk and checked for Mendelian consistency with SimWalk2. A SNP was eliminated if: 1) fewer than 95% of the cohort were typed successfully; 2) the SNP was monoallelic; 3) the SNP had more than two alleles; 4) fewer than five copies of the SNP existed in the current study cohort. Hardy-Weinberg equilibrium (HWE) was tested for each SNP using SOLAR ; SNPs with excessive deviation from HWE (p < 10-8) were excluded. Missing SNP data were imputed with MERLIN .
Genome-wide transcriptional profiles of a subset of the SNP genotyping cohort (369 individuals from 55 nuclear families) were obtained as previously described  with modifications. Briefly, for each individual 2.5 ml blood was collected into a PAXgene® Blood RNA Tube (BD, Franklin Lakes, NJ) following an overnight fast. Total RNA was isolated from each tube using the PAXgene Blood RNA Kit (Qiagen, Valencia, CA) and anti-sense RNA (aRNA) was synthesized using the MessageAmp II-Biotin aRNA kit (Ambion, Austin, TX). A total of 1.5 μg aRNA was hybridized to Illumina HumanWG-6 version 2 or version 3 chips (Illumina, San Diego, CA) and expression detected on the Illumina® BeadArray™ 500GX Reader. Illumina GenomeStudio software (version 2010.3) was used for preliminary data analysis with standard background normalization. The data discussed in this publication have been deposited in NCBI’s Gene Expression Omnibus  and are accessible through GEO Series accession number GSE45987 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE45987).
Measured genotype analysis
In SOLAR, each SNP genotype was converted to a covariate measure equal to 0, 1 or 2 copies of the minor allele (or, for missing genotypes, a weighted covariate based on imputation). These covariates were included in variance-components mixed models for measured genotype analyses  versus null models that incorporated the random effect of kinship and fixed effects such as age, age2, sex and their interactions. Individual scores from a principal components analysis of representative SNPs were also included to correct for possible population stratification .
Within chromosomal regions showing the strongest prior evidence of linkage (as logarithm of odds, LOD, score) with MetS traits in our MRC-OB cohort, we selected all available SNPs on the Affymetrics 6.0 array that mapped within a 1-LOD confidence interval of the maximum LOD score in each region. Each SNP covariate was tested independently in a 1 degree of freedom likelihood ratio test. To take into account the linkage disequilibrium between SNPs in each region, we calculated the effective number of independent SNPs (Neff) using the method of Moskvina & Schmidt . Critical p-values were calculated using Bonferroni/Šidák correction for each linkage region based on its Neff.
Gene expression analysis
Microarray data were available in two batches, one based on Version 2 arrays (48,701 probes, 307 samples) and the other on Version 3 (48,803 probes, 230 samples). To guard against possible batch effects and probe differences, each batch was analyzed separately: The number of probe transcripts detectable at p≤0.05 by BeadStudio software was counted, a false discovery rate (FDR) was computed across all probes, and transcripts detectable at 5% FDR were retained. Expression levels were log2 transformed and inverse-quartile normalized. Transformed and normalized expression levels for probes that mapped to the 1-LOD QTL regions were tested for association with phenotypes of interest in models that included the random effect of kinship. Gene-centric p-values were calculated by combining independent p-values from the two microarray batches and multiple probes using Stouffer’s weighted Z-score method  implemented in R.
Intercorrelation between plasma adiponectin levels and MetS phenotypes
Based on their attributes in different biological pathways, we categorized our 42 measured MetS phenotypes into five groups: body composition, insulin responsiveness, lipids/lipoprotein profiles, cardiovascular performance, and adipokine and cytokines (Table 1). Means ± SD of the 42 MetS phenotypes and levels of their heritability have been previously described . We found that plasma adioponectin levels were highly heritable (additive heritability, h2 = 0.48) in our study cohort. Intercorrelation analysis of leading MetS component phenotypes revealed significant phenotypic correlations to varied degrees (See the Additional file 1: Table S1). Plasma adiponectin is also significantly correlated with the majority of tested phenotypes (29 out of 41, nominal p<0.05). In particular, it is highly correlated with the MetS-defining [27, 28] phenotypes WC (correlation p-value=1.86×10-22), FG (7.8×10-4), plasma TG (1.3×10-8), plasma HDL-c (5.2×10-12), sBP (1.55×10-6) and dBP (7.1×10-3). Overall, levels of adiponectin inversely correlate with a phenotypes representing adverse adipose distribution (weight, BMI, WC, HC, WHR, Fatkg, Fatpct and SubQF, VF and TAF), elevated fasting glucose levels and insulin resistance indices (FG, FI, IGR and HOMA), dyslipidemia (TG), cardiovascular malfunction (sBP, dBP and pulse) and leptin. On the contrary, adiponectin levels correlate positively with several phenotypes indexing leanness (Leanpct), insulin sensitivity (SI and SG) and HDL-c. These results indicate that adiponectin is a reliable and significant phenotypic surrogate of both the clinical and biological traits of MetS in our cohort.
SNP genotype/phenotype association
Consistent with common practice, we defined confidence regions for the previously defined adiponectin linkage peaks as the physical positions adjacent to the peak showing linkage evidence of at least LOD score - 1: 5p14 [peak LOD=4.1; chr5: 9,792,000-23,021,100bp (NCBI36/hg18)] and 14q13 [peak LOD=3.2; chr14: 23,131,000-36,761,868bp (NCBI36/hg18)]. Haplotype-tagging SNPs located within these regions were tested for statistical association with each of the 42 MetS phenotypes. To account for multiple testing within QTLs while accounting for linkage disequilibrium (LD) among SNPs, QTL-specific significance thresholds were calculated based upon Bonferroni/Šidák correction for the effective number of independent SNPs at 5p14 (nSNPs =3540, neffective= 2429.61, pα=0.05= 2.11×10-5) and 14q13 (nSNPs=3794, neffective =2760.37, pα=0.05=1.86×10-5). SNPs that passed region-specific significance thresholds for respective traits and/or the highest nominally significant SNP of each trait for each trait are shown in Tables 2 and 3.
Fine-mapping of this region identified QTL-wide significant SNP associations with body composition phenotypes, including weight and BMI (Figure 1A and 1B). As shown in Table 2, 7 SNPs led by rs10061119 are mapped to an intronic region of the gene encoding cadherin 18 (CDH18). In addition, SNPs annotated near myosin X (MYO10), showed nominally significant associations with plasma adiponectin (Table 2 and Figure 1C). The second strongest SNP association cluster with adiponectin (Figure 1C, pointed by a red arrow) is intergenic and is approximately 400kb 5’ to the transcription start site of the cadherin-18 gene CDH18 (Figure 1C, framed in green).
SNPs mapped to CDH18 (cadherin 18) were also nominally associated with other MetS phenotypes including WC, HC and RQ. SNP rs404639 mapped to an intron of CDH12, another member of the cadherin family, and was associated with plasma leptin. SNP rs6894869 was significantly associated with lipid particle size phenotype BMED and mapped to an intron of CTNND2, which functions in the same pathway as cadherins, plays a crucial role in CNS development  and is a partner of the Alzheimer disease-associated gene presenilin-1 . MYO10 was also nominally associated with Leankg, LDLppd and IL-6. For the other MetS traits the strongest associations were with SNPs near several other genes with biological relevance to cancer. SNPs near TAG, a tumor antigen gene, were associated with HMED (high density lipid particle sizing), REE (resting energy expenditure), and insulin phenotypes FG and HOMA. SNPs mapped to MARCH11 and FBXL7, two ubiquitin ligases involved in the pro-inflammation pathways, were associated with total fat mass, FI, SI and pro-inflammation marker IL-1beta (MARCH11), and with lean body mass (FBXL7).
This QTL region has a high density of genes with regulatory functions including signal transduction, transcription and post-transcription processing. Our fine-mapping identified SNPs associated with each of the 42 MetS phenotypes. Table 3 summarizes the SNPs most highly associated with each trait, the attributes of the variants, and the gene ontology of their annotated genes. Variants of NPAS3, a neuronal transcription factor thought to be involved in brain tumor suppression, were found to be associated with 11 of the MetS phenotypes (Fatkg, Leankg, SubQF, TAF, REE, TG, TC, LDL-c, pulse, IL1b and IL-6), suggesting a pleitropic effect (Table 3). Two of these associations (TG and TC) were borderline significant, marginally reaching the regional threshold (Table 3 and Figure 2A). SNPs near BRMS1L, a breast cancer suppressor gene, was associated with LMEDn and LDLppd. SNPs of EGLN3, a regulator of transcription factor HIF that affects apoptosis in hemangioblastoma and clear cell renal cancer, were associated with IGR and LDL-c. SNPs of NOVA1, a post-transcription processing of the GnRH paraneoplastic antigen, were associated at nominal significance with several phenotypes including weight, BMI, WC, FG, FI, HOMA, SI, AIR, and sBP.
Prominent SNP associations in this region were mapped to genes for signal transducers include kinases and kinase-related proteins AKAP6 and PRKD1 and the phosphotase PPP2R3C. SNPs of the gene encoding AKAP6, an anchor protein for AMPK, which may connect MetS with cancer, were associated with HC, RQ and REE/leanmass. SNPs of PRKD1, a kinase that phosphorylates and therefore activates E-cadherin functions, were associated with WHR, Fatpct, Leanpct, and TNF-alpha. SNPs of PPP2R3C, a protein phosphatase, were associated with DI. SNP rs9322942 of SRP54, the gene for a signal recognition particle protein, is associated with adiponectin. Figure 2 shows the SNP association patterns with phenotypes RQ and TG at this QTL.
The cis-SNP effects on gene expression
To test the function of the interrogated SNPs of the two adiponectin QTLs, we evaluated the effect of each variant on the levels of transcripts produced from the genes located within these two regions. We surveyed transcripts that were robustly detectable after filtering through our quality control procedures (see Methods). We then assessed the association between QTL-specific SNPs and each transcript annotated to the same locus.
Figure 3 shows the effects of SNPs of the chr5 locus on expression of MYO10. We found strong signals clustering in proximity to the transcription start site (TSS) of MYO10. SNPs exhibiting the strongest cis-effects (lowest p-value=5.43×10-6) are in high linkage disequilibrium (data not shown) and span a 17 kb region that has been shown previously to be highly important for transcription initiation [31–34].
We further analyzed the relationship between MYO10 gene expression patterns in PWBC and key MetS phenotypes for which we had found evidence of SNP-phenotype associations (Table 2). Table 4 shows the levels of the correlations between the level of MYO10 gene expression and MetS traits. We found positive correlations of MYO10 expression with measures of adiposity. MYO10 expression correlated negatively with plasma levels of adiponectin (although this did not reach statistical significance) and was positively correlated with plasma leptin at
Summary of candidate gene prioritization
Table 5 summarizes all the biologically relevant cancer or pro-inflammation genes, annotated by the SNP associations and gene expression correlations with MetS phenotypes under the adiponectin QTL peaks of 5p14 and 14q13. Figures 1 and 2 show SNP/phenotype association plots of SNPs located within these QTLs that we consider to be of the highest priority based on their biological relevance, in relation to the phenotypes from which the top candidate gene was identified.
As previously seen in other studies [36, 37] plasma adiponectin was significantly correlated with several of the clinical outcome and biologic precursor phenotypes of MetS in our cohort. It correlated with BMI, total body fat mass, visceral vs. subcutaneous fat distribution, insulin to glucose response indices, levels of plasma lipids/lipoprotein density profiles, and circulating levels of leptin and cytokines (Table 1). As expected, the key clinical traits of MetS are also inter-correlated with each other in our subjects (Additional file 1: Table S1). Using adiponectin level as a surrogate for MetS, we previously identified two adiponectin-linked QTLs (5p14 and 14q13) in our study population. We now have used SNP association and the cis-effects of local SNPs on PWBC expression to discover MetS-related genes within these two QTLs. By doing so, we were able to identify multiple novel genes that have significant associations with individual MetS phenotypes. These genes can be mechanistically grouped into four categories (Table 5): (1) cell-cell adhesion and mobility (CDH18, CDH12, CTNND2 and MYO10) (2) signal transduction (AKAP, PRKD1 and PPP2R3C); (3) transcription (NPAS3, BRMS1, Nova-1, EGLN3 and TAG); and (4) protein sorting (MARCH11 and FBXL7). It is interesting that many of these genes are known to play a role in cancer biology. The link between MetS and increased cancer risk is well known  thus our prioritized genes provide a possible explanation for the observed association between MetS and cancer.
Our strongest SNPs associations on the chr5 QTL clustered around genes in the cell-cell adhesion pathways (Table 2 and Figure 1). Cell to cell adhesion is a key event in major cellular processes including proliferation, mobility, differentiation and cell death. It is also a basic mechanism facilitating the paracrine communication between cells. CDH18 is a cadherin that is expressed in multiple tissues, most prominently in the CNS  and is one of several cadherin genes located in this region of chromosome 5. The CDH18 gene was most strongly associated with weight, BMI, waist circumference and RQ in our cohort, indicating a possible role in the development of body composition. Mechanistically this could be mediated via a role in CNS but since cell adhesion is a phenomenon that is integral to the regional expansion of adipose tissue it is possible that the gene may play a direct role in the preferential deposition of fat into the visceral abdominal region, a fundamental phenotype of MetS. Aberrations leading to alteration of cell adhesion are linked with several tumor types. N-cadherins specifically have been reported to play a role in the formation, growth, invasion  and migration of malignant tumors in various settings, including Crohn’s disease  and prostate cancer [42, 43]. Deletions in the region of the chromosome 5 cadherin cluster have been specifically associated with risk of some malignancies . Thus it is possible that the CDH18 gene and other nearby genes identified in cohort (including CTNND2, CDH12 and MYO10) play a mechanistic role in the observed connection between MetS and cancer risk.
Our strongest SNP association signals for plasma adiponectin levels clustered near the gene MYO10 on chr5. MYO10 protein is a member of the myosin family and has a role in trafficking adhesion molecules including integrins and adherins . It is essential for the initiation, stability and formation of filopodia, a “finger-like” cellular protrusion that can sense environmental cues and a key structure for cell motility including cancer cells . In a recent study of breast cancer with poor prognosis, expression of MYO10 is significantly increased in patients with breast cancer. The expression of MYO10 also differentially clustered with clinocopathological markers including p53 mutation, estrogen receptor (ER), tumor grade and patient survival . In our analysis, variants in MYO10 also associated with lean body mass, LDL diameter and IL-6, thus reflecting pleiotropy in its role in MetS, possibly via modulation of adiponectin levels. No data has shown that MYO10 directly transports CDH18, however our data suggest that these two genes might work together in a biological pathway that predisposes individuals with MetS for carcinogenesis or high risk of metastasis. Our evaluation of the function of these tagging SNPs (or SNPs in LD with these) suggested strong cis-effects of variants near the TSS of MYO10. This region thus may harbor causal SNP(s) that directly affect binding of transcription factors. Results from human cell lines recently published by the ENCODE project has demonstrated that this 17 kb region is marked with transcription-activating histone modifications, is highly accessible to protein binding, and can be bound by multiple transcription factors in chromatin immunoprecipitation assays (Figure 3) [33, 34]. Our approach therefore is able to efficiently identify genetic elements that can function in influencing genes important for our traits of interest.
At the QTL on chr14, SNP associations clustered near genes that function as signal transducers like AKAP6 and transcription factors like NPAS3 (Table 3 and Figure 2). The protein encoded by AKAP6 is an anchor protein for the well-characterized protein kinase A (PKA), also known as AMP-activated protein kinase (AMPK). AKAPs compartmentalize AMPK, thereby regulating AMPK’s biological functions. AKAPs by themselves may also function as integrators of AMPK-related signal transduction pathways . AMPK is a key regulator sensing cellular energy status and systemic energy balance  and mediates the effects of adipose tissue kinins including adiponectin and leptin in regulating body weight, glucose and lipid homeostasis . It may also play a role in cancer via several processes including arresting of cell cycle progression, triggering inhibition of protein synthesis, and cell growth required in cell proliferation . AMPK is also known to overcome growth factor signaling from a variety of stimuli mediated by the proto-oncogenes Akt and ERK. The AMPK signaling network also contains a number of tumor suppressor genes including p53, which is a universal tumor suppressor . AMPK is thus one of the most prominent biological molecules proposed to connect metabolic pathways of MetS to cancer .
NPAS3, associated with several lipid and cytokine phenotypes in our cohort, encodes a neuronal basic helix-loop-helix (Per, Arnt, Sim) domain transcription factor, belonging to a family that has been reported to be crucial for neuronal development and maintenance  and exerts a regulatory role in neurogenesis and brain glucose metabolism . Brain tissue of an NPAS3 knockout mouse model shows altered levels of the glycolysis metabolites, the pentose phosphate shunt and Krebs cycle components, and aberrant behavior in mice. Glucose sensing by the brain is known to have effects on feeding, behavior, weight, weight control and associated metabolic consequences, as in MetS. NPAS3 can also modulate cell cycle proliferation, apoptosis, migration and cell invasion and exhibits features of a tumor suppressor. Absence of NPAS3 expression in astrocytomas is a negative prognostic marker for survival .
Increasing evidence has suggested that peripheral blood can be used to reveal the relative differences in gene expression among individuals [54, 55]. Obtaining gene expression information using peripheral blood is simple and minimally invasive, making it possible to assay a large number of our subjects. Furthermore, obesity can be considered as a low-grade inflammation state , due to long-term low-level immune response triggered by accumulation of macrophages in visceral adipose tissue, chronic insulin resistance, adiponkine and cytokine production and malfunctional cardiovascular performance [57, 58]. Profiling gene expression in peripheral blood therefore may capture changes in the expression of genes important for components of MetS such as visceral adiposity. In our analysis, we found some significant signals of SNP associations near pro-inflammatory genes. MARCH11 (associated with fat, FI, SI and IL-1beta in our cohort) and FBXL7 (associated with WHR in our cohort) belong to the family of ubiquitin ligases that are engaged in the process of protein sorting and transport of deformed ones from the trans-Golgi network [59, 60]. MARCH11 expression was shown in the Gene Expression Omnibus (GEO) database to be higher in patients with inflammatory dilated cardiomyopathy. FBXL7, an F-box protein family member, is specific for sorting of phosphorylated proteins and their degradation. Thus MARCH11 and FBXL7 may be involved in inflammation, a characteristic phenotype of MetS.
We do acknowledge that PWBC can reveal only a portion of the total gene expression information and that profiling it in other focal tissues like adipose and liver might reveal complementary genomic information. Further exploration of the expression levels of the other genes associated with MetS (and known to play a role in cancer biology) in other target tissues might help to reveal their gene expression correlations with MetS traits.
Our study has identified several novel as well as previously reported genes whose SNP variants are associated with phenotypes that are important components of MetS. e.g. genome-wide study of the genetic basis of insulin resistance reported that a variant of AKAP6 was associated with BMI-adjusted fasting insulin (5×10-7) in a meta-analysis cohort of European descent . NPAS3 was found in the current study to be associated with a multitude of phenotypes including two inflammatory markers IL-1β and IL-6. Interestingly, in a recent report of a genome-wide search for associations between SNPs and biomarkers of systemic inflammation before and after an anti-inflammatory drug treatment, a variant of NPAS3 was found to be one of the two most significantly associated gene loci with levels of C-reactive protein , and studies of visceral adiposity and BMI identified FBXL7, MARCH6 and PRKD1 as candidate loci [63–65].
In the present study of a family cohort of Northern European origin, we conducted a comprehensive analysis using both SNP associations and gene expression correlations on two previously identified adiponectin QTLs and identified 14 novel genes associated with various metabolic syndrome traits. We found that many of these genes are involved in fundamental cellular functions that could provide a mechanistic link between the development of metabolic syndrome and cancer biology as well as systemic inflammation. Further studies are needed to elucidate the roles these genes play in the development of metabolic syndrome and its relationship with cancer biology. Eventually this could lead to the design of novel diagnostic markers and/or pharmaceutical agents for the prevention and treatment of MetS patients who are at higher risks for developing cancer and systemic inflammation.
Day C: Metabolic syndrome, or What you will: definitions and epidemiology. Diab Vasc Dis Res. 2007, 4 (1): 32-38. 10.3132/dvdr.2007.003.
Russo A, Autelitano M, Bisanti L: Metabolic syndrome and cancer risk. Eur J Cancer. 2008, 44 (2): 293-297. 10.1016/j.ejca.2007.11.005.
Talvensaari KK, Lanning M, Tapanainen P, Knip M: Long-term survivors of childhood cancer have an increased risk of manifesting the metabolic syndrome. J Clin Endocrinol Metab. 1996, 81 (8): 3051-3055. 10.1210/jc.81.8.3051.
Chandran M, Phillips SA, Ciaraldi T, Henry RR: Adiponectin: more than just another fat cell hormone. Diab Care. 2003, 26: 2442-2450. 10.2337/diacare.26.8.2442.
Pajvani UB, Hawkins M, Combs TP, Rajala MW, Doebber T, Berger JP, Wagner JA, Wu M, Knopps A, Xiang AH, Utzschneider KM, Kahn SE, Olefsky JM, Buchanan TA, Scherer PE: Complex distribution, not absolute amount of adiponectin, correlates with thiazolidinedione-mediated improvement in insulin sensitivity. J Biol Chem. 2004, 279: 12152-12162.
Ouchi N, Walsh K: Adiponectin as an anti-inflammatory factor. Clin Chim Acta. 2007, 380 (1–2): 24-30.
Kubota N, Terauchi Y, Yamauchi T, Kubota T, Moroi M, Matsui J, Eto K, Yamashita T, Kamon J, Satoh H, Yano W, Froguel P, Nagai R, Kimura S, Kadowaki T, Noda T: Disruption of adiponectin causes insulin resistance and neointimal formation. J Biol Chem. 2002, 277 (29): 25863-25866. 10.1074/jbc.C200251200.
Magge SN, Stettler N, Koren D, Levitt Katz LE, Gallagher PR, Mohler ER, Rader DJ: Adiponectin is associated with favorable lipoprotein profile, independent of BMI and insulin resistance, in adolescents. J Clin Endocrinol Metab. 2011, 96 (5): 1549-1554. 10.1210/jc.2010-2364.
Kelesidis I, Kelesidis T, Mantzoros CS: Adiponectin and cancer: a systematic review. Br J Cancer. 2006, 94 (9): 1221-1225. 10.1038/sj.bjc.6603051.
Comuzzie AG, Funahashi T, Sonnenberg G, Martin LJ, Jacob HJ, Black AE, Maas D, Takahashi M, Kihara S, Tanaka S, Matsuzawa Y, Blangero J, Cohen D, Kissebah A: The genetic basis of plasma variation in adiponectin, a global endophenotype for obesity and the metabolic syndrome. J Clin Endocrinol Metab. 2001, 86 (9): 4321-4325. 10.1210/jc.86.9.4321.
Kissebah AH, Sonnenberg GE, Myklebust J, Goldstein M, Broman K, James RG, Marks JA, Krakower GR, Jacob HJ, Weber J, Martin L, Blangero J, Comuzzie AG: Quantitative trait loci on chromosomes 3 and 17 influence phenotypes of the metabolic syndrome. Proc Natl Acad Sci U S A. 2000, 97 (26): 14478-14483. 10.1073/pnas.97.26.14478.
Sonnenberg GE, Krakower GR, Martin LJ, Olivier M, Kwitek AE, Comuzzie AG, Blangero J, Kissebah AH: Genetic determinants of obesity-related lipid traits. J Lipid Res. 2004, 45 (4): 610-615. 10.1194/jlr.M300474-JLR200.
Svendsen OL, Haarbo J, Heitmann BL, Gotfredsen A, Christiansen C: Measurement of body fat in elderly subjects by dual-energy x-ray absorptiometry, bioelectrical impedance, and anthropometry. Am J Clin Nutr. 1991, 53 (5): 1117-1123.
Peiris AN, Hennes MI, Evans DJ, Wilson CR, Lee MB, Kissebah AH: Relationship of anthropometric measurements of body fat distribution to metabolic profile in premenopausal women. Acta Med Scand Suppl. 1988, 723: 179-188.
Bergman RN: Toward physiological understanding of glucose tolerance. Minimal-model approach. Diabetes. 1989, 38: 1512-1527. 10.2337/diabetes.38.12.1512.
Rainwater DL, Moore PH, Shelledy WR, Dyer TD, Slifer SH: Characterization of a composite gradient gel for the electrophoretic separation of lipoproteins. J Lipid Res. 1997, 38 (6): 1261-1266.
Dupont NC, Wang K, Wadhwa PD, Culhane JF, Nelson EL: Validation and comparison of luminex multiplex cytokine analysis kits with ELISA: determinations of a panel of nine cytokines in clinical sample culture supernatants. J Reprod Immunol. 2005, 66 (2): 175-191. 10.1016/j.jri.2005.03.005.
Almasy L, Blangero J: Multipoint quantitative-trait linkage analysis in general pedigrees. Am J Hum Genet. 1998, 62 (5): 1198-1211. 10.1086/301844.
Burdick JT, Chen W-M, Abecasis GR, Cheung VG: In silico method for inferring missing genotypes in pedigrees. Nat Genet. 2006, 38: 1002-1004. 10.1038/ng1863.
Göring HH, Curran JE, Johnson MP, Dyer TD, Charlesworth J, Cole SA, Jowett JBM, Abraham LJ, Rainwater DL, Comuzzie AG, Mahaney MC, Almasy L, MacCluer JW, Kissebah AH, Collier GR, Moses EK, Blangero J: Discovery of expression QTLs using large-scale transcriptional profiling in human lymphocytes. Nat Genet. 2007, 39 (10): 1208-1216. 10.1038/ng2119.
Edgar R, Domrachev M, Lash AE: Gene expression omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002, 30 (1): 207-210. 10.1093/nar/30.1.207.
Boerwinkle E, Chakraborty R, Sing CF: The use of measured genotype information in the analysis of quantitative phenotypes in man. I. Models and analytical methods. Ann Hum Genet. 1986, 50: 181-194. 10.1111/j.1469-1809.1986.tb01037.x.
Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, Reich D: Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet. 2006, 38: 904-909. 10.1038/ng1847.
Moskvina V, Schmidt KM: On multiple testing correction in genome-wide association studies. Genet Epidemiol. 2008, 32: 567-573. 10.1002/gepi.20331.
Whitlock MC: Combining probability from independent tests: the weighted Z-method is superior to Fisher’s approach. J Evol Biol. 2005, 18: 1368-1373. 10.1111/j.1420-9101.2005.00917.x.
Zhang Y, Kent JW, Olivier M, Ali O, Broeckel U, Abdou RM, Dyer TD, Comuzzie AE, Curran JE, Carless MA, Rainwater DL, Goring HHH, Blangero J, Kissebah AH: QTL-based association analyses reveal novel genes influencing pleiotropy of Metabolic Syndrome (MetS). Obesity. 2013, accepted
May Executive Summary of The Third Report of The National Cholesterol Education Program (NCEP): Expert panel on detection, evaluation, and treatment of high blood cholesterol in adults (adult treatment panel III). JAMA. 2001, 285 (19): 2486-2497. 10.1001/jama.285.19.2486.
Alberti KG, Eckel RH, Grundy SM, Zimmet PZ, Cleeman JI, Donato KA, Fruchart JC, James WP, Loria CM, Smith SC, International Diabetes Federation Task Force on Epidemiology and Prevention; Hational Heart, Lung, and Blood Institute; American Heart Association; World Heart Federation; International Atherosclerosis Society; International Association for the Study of Obesity: Harmonizing the metabolic syndrome: a joint interim statement of the international diabetes federation task force on epidemiology and prevention; national heart, lung, and blood institute; American heart association; world heart federation; international atherosclerosis society; and international association for the study of obesity. Circulation. 2009, 120 (16): 1640-1645. 10.1161/CIRCULATIONAHA.109.192644.
Ho C, Zhou J, Medina M, Goto T, Jacobson M, Bhide PG, Kosik KS: Delta-catenin is a nervous system-specific adherens junction protein which undergoes dynamic relocalization during development. J Comp Neurol. 2000, 420 (2): 261-276. 10.1002/(SICI)1096-9861(20000501)420:2<261::AID-CNE8>3.0.CO;2-Q.
Zhou J, Liyanage U, Medina M, Ho C, Simmons AD, Lovett M, Kosik KS: Presenilin 1 interaction in the brain with a novel member of the Armadillo family. Neuroreport. 1997, 8 (8): 2085-2090. 10.1097/00001756-199705260-00054.
Kent WJ, Sugnet CW, Furey TS, Roskin KM, Pringle TH, Zahler AM, Haussler D: The human genome browser at UCSC. Genome Res. 2002, 12 (6): 996-1006.
UCSC Genome Database. http://genome.ucsc.edu/.
Myers RM, Stamatoyannopoulos J, Snyder M, Dunham I, Hardison RC, Bernstein BE, Gingeras TR, Kent WJ, Birney E, ENCODE Project Consortium: A user’s guide to the encyclopedia of DNA elements (ENCODE). PLoS Biol. 2011, 9 (4): e1001046-10.1371/journal.pbio.1001046.
Rosenbloom KR, Dreszer TR, Long JC, Malladi VS, Sloan CA, Raney BJ, Cline MS, Karolchik D, Barber GP, Clawson H, Diekhans M, Fujita PA, Goldman M, Gravell RC, Harte RA, Hinrichs AS, Kirkup VM, Kuhn RM, Learned K, Maddren M, Meyer LR, Pohl A, Rhead B, Wong MC, Zweig AS, Haussler D, Kent WJ: ENCODE whole-genome data in the UCSC genome browser: update 2012. Nucleic Acids Res. 2011, 40 (Database issue): 1-6.
Pruim RJ, Welch RP, Sanna S, Teslovich TM, Chines PS, Gliedt TP, Boehnke M, Abecasis GR, Willer CJ: LocusZoom: Regional visualization of genome-wide association scan results. Bioinformatics. 2010, 26 (18): 2336-2337. 10.1093/bioinformatics/btq419.
Koh SB, Yoon J, Kim JY, Yoo BS, Lee SH, Park JK, Choe KH: Relationships between serum adiponectin with metabolic syndrome and components of metabolic syndrome in non-diabetic Koreans: ARIRANG study. Yonsei Med J. 2011, 52 (2): 234-241. 10.3349/ymj.2011.52.2.234.
Patel DA, Srinivasan SR, Xu JH, Chen W, Berenson GS: Adiponectin and its correlates of cardiovascular risk in young adults: the Bogalusa heart study. Metabolism. 2006, 55 (11): 1551-1557. 10.1016/j.metabol.2006.06.028.
Esposito K, Chiodini P, Colao A, Lenzi A, Giugliano D: Metabolic syndrome and risk of cancer: a systematic review and meta-analysis. Diabetes Care. 2012, 35 (11): 2402-2411. 10.2337/dc12-0336.
Gene Expression Atlas at the European Bioinformatics Institute. http://www.ebi.ac.uk/gxa/gene/ENSG00000145526?ef=cell_type.
Hazan RB, Phillips GR, Qiao RF, Norton L, Aaronson SA: Exogenous expression of N-cadherin in breast cancer cells induces cell migration, invasion, and metastasis. J Cell Biol. 2000, 148 (4): 779-790. 10.1083/jcb.148.4.779.
Burke JP, Cunningham MF, Sweeney C, Docherty NG, O’Connell PR: N-cadherin is overexpressed in Crohn‘s stricture fibroblasts and promotes intestinal fibroblast migration. Inflamm Bowel Dis. 2011, 17 (8): 1665-1673. 10.1002/ibd.21543.
Gravdal K, Halvorsen OJ, Haukaas SA, Akslen LA: A switch from E-cadherin to N-cadherin expression indicates epithelial to mesenchymal transition and is of strong and independent importance for the progress of prostate cancer. Clin Cancer Res. 2007, 13 (23): 7003-7011. 10.1158/1078-0432.CCR-07-1263.
Tanaka H, Kono E, Tran CP, Miyazaki H, Yamashiro J, Shimomura T, Fazli L, Wada R, Huang J, Vessella RL, An J, Horvath S, Gleave M, Rettig MB, Wainberg ZA, Reiter RE: Monoclonal antibody targeting of N-cadherin inhibits prostate cancer growth, metastasis and castration resistance. Nat Med. 2010, 16 (12): 1414-1420. 10.1038/nm.2236.
Kajikawa Y, Ikeda M, Takemoto S, Tomoda J, Ohmaru N, Kusachi S: Association of circulating levels of leptin and adiponectin with metabolic syndrome and coronary heart disease in patients with various coronary risk factors. Int Heart J. 2011, 52 (1): 17-22. 10.1536/ihj.52.17.
Almagro S, Durmort C, Chervin-Pétinot A, Heyraud S, Dubois M, Lambert O, Maillefaud C, Hewat E, Schaal JP, Huber P, Gulino-Debrac D: The motor protein myosin-X transports VE-cadherin along filopodia to allow the formation of early endothelial cell-cell contacts. Mol Cell Biol. 2010, 30 (7): 1703-1717. 10.1128/MCB.01226-09.
Arjonen A, Kaukonen R, Ivaska J: Filopodia and adhesion in cancer cell motility. Cell Adh Migr. 2011, 5 (5): 421-430. 10.4161/cam.5.5.17723.
Michel JJ, Scott JD: AKAP mediated signal transduction. Annu Rev Pharmacol Toxicol. 2002, 42: 235-257. 10.1146/annurev.pharmtox.42.083101.135801.
Kahn BB, Alquier T, Carling D, Hardie DG: AMP-activated protein kinase: ancient energy gauge provides clues to modern understanding of metabolism. Cell Metab. 2005, 1 (1): 15-25. 10.1016/j.cmet.2004.12.003.
Carretero J, Medina PP, Blanco R, Smit L, Tang M, Roncador G, Maestre L, Conde E, Lopez-Rios F, Clevers HC, Sanchez-Cespedes M: Dysfunctional AMPK activity, signalling through mTOR and survival in response to energetic stress in LKB1-deficient lung cancer. Oncogene. 2007, 26 (11): 1616-1625. 10.1038/sj.onc.1209951.
Motoshima H, Goldstein BJ, Igata M, Araki E: AMPK and cell proliferation–AMPK as a therapeutic target for atherosclerosis and cancer. J Physiol. 2006, 574: 63-71. 10.1113/jphysiol.2006.108324.
Rutter J, Reick M, Wu LC, McKnight SL: Regulation of clock and NPAS2 DNA binding by the redox state of NAD cofactors. Science. 2001, 293 (5529): 510-514. 10.1126/science.1060698.
Sha L, Macintyre L, Machell JA, Kelly MP, Porteous DJ, Brandon NJ, Muir WJ, Blackwood DH, Watson DG, Clapcote SJ, Pickard BS: Transcriptional regulation of neurodevelopmental and metabolic pathways by NPAS3. Mol Psychiatry. 2012, 17: 267-279. 10.1038/mp.2011.73.
Moreira F, Kiehl TR, So K, Ajeawung NF, Honculada C, Gould P, Pieper RO, Kamnasaran D: NPAS3 demonstrates features of a tumor suppressive role in driving the progression of Astrocytomas. Am J Pathol. 2011, 179 (1): 462-476. 10.1016/j.ajpath.2011.03.044.
Charlesworth JC, Curran JE, Johnson MP, Göring HHH, Dyer TD, Diego VP, Kent JW, Mahaney MC, Almasy L, MacCluer JW, Moses EK, Blangero J: Transcriptomic epidemiology of smoking: the effect of smoking on gene expression in lymphocytes. BMC Medical Genomics. 2010, 3: 29-10.1186/1755-8794-3-29.
Yamaoka M, Maeda N, Nakamura S, Kashine S, Nakagawa Y, Hiuge-Shimizu A, Okita K, Imagawa A, Matsuzawa Y, Matsubara K, Funahashi T, Shimomura I: A pilot investigation of visceral fat adiposity and gene expression profile in peripheral blood cells. PLoS One. 2012, 7 (10): e47377-10.1371/journal.pone.0047377.
Scarpellini E, Tack J: Obesity and metabolic syndrome: an inflammatory condition. Dig Dis. 2012, 30 (2): 148-153. 10.1159/000336664.
Hotamisligil GS: Endoplasmic reticulum stress and the inflammatory basis of metabolic disease. Cell. 2010, 140: 900-917. 10.1016/j.cell.2010.02.034.
Neels JG, Olefsky JM: Inflamed fat: what starts the fire?. J Clin Invest. 2006, 116: 33-35.
Morokuma Y, Nakamura N, Kato A, Notoya M, Yamamoto Y, Sakai Y, Fukuda H, Yamashina S, Hirata Y, Hirose S: MARCH-XI, a novel transmembrane ubiquitin ligase implicated in ubiquitin-dependent protein sorting in developing spermatids. J Biol Chem. 2007, 282 (34): 24806-24815. 10.1074/jbc.M700414200.
Cenciarelli C, Chiaur DS, Guardavaccaro D, Parks W, Vidal M, Pagano M: Identification of a family of human F-box proteins. Curr Biol. 1999, 9 (20): 1177-1179. 10.1016/S0960-9822(00)80020-2.
Manning AK, Hivert MF, Scott RA, Grimsby JL, Bouatia-Naji N, Chen H, Rybin D, Liu CT, Bielak LF, Prokopenko I, Amin N, Barnes D, Cadby G, Hottenga JJ, Ingelsson E, Jackson AU, Johnson T, Kanoni S, Ladenvall C, Lagou V, Lahti J, Lecoeur C, Liu Y, Martinez-Larrad MT, Montasser ME, Navarro P, Perry JR, Rasmussen-Torvik LJ, Salo P, Sattar N: A genome-wide approach accounting for body mass index identifies genetic variants influencing fasting glycemic traits and insulin resistance. Nat Genet. 2012, 44 (6): 659-669. 10.1038/ng.2274.
Aslibekyan S, Kabagambe EK, Irvin MR, Straka RJ, Borecki IB, Tiwari HK, Tsai MY, Hopkins PN, Shen J, Lai CQ, Ordovas JM, Arnett DK: A genome-wide association study of inflammatory biomarker changes in response to fenofibrate treatment in the genetics of lipid lowering drug and diet network. Pharmacogenet Genomics. 2012, 22 (3): 191-197. 10.1097/FPC.0b013e32834fdd41.
Fox CS, Liu Y, White CC, Feitosa M, Smith AV, Heard-Costa N, Lohman K, Johnson AD, Foster MC, Greenawalt DM, Griffin P, Ding J, Newman AB, Tylavsky F, Miljkovic I, Kritchevsky SB, Launer L, Garcia M, Eiriksdottir G, Carr JJ, Gudnason V, Harris TB, Cupples LA, Borecki IB, GIANT Consortium; MAGIC Consortium; GLGC Consortium: Genome-wide association for abdominal subcutaneous and visceral adipose reveals a novel locus for visceral fat in women. PLoS Genet. 2012, 8 (5): e1002695-10.1371/journal.pgen.1002695.
Wang KS, Liu X, Zheng S, Zeng M, Pan Y, Callahan K: A novel locus for body mass index on 5p15.2: a meta-analysis of two genome-wide association studies. Gene. 2012, 500 (1): 80-84. 10.1016/j.gene.2012.03.046.
Speliotes EK, Willer CJ, Berndt SI, Monda KL, Thorleifsson G, Jackson AU, Lango Allen H, Lindgren CM, Luan J, Mägi R, Randall JC, Vedantam S, Winkler TW, Qi L, Workalemahu T, Heid IM, Steinthorsdottir V, Stringham HM, Weedon MN, Wheeler E, Wood AR, Ferreira T, Weyant RJ, Segrè AV, Estrada K, Liang L, Nemesh J, Park JH, Gustafsson S, Kilpeläinen TO: Association analyses of 249,796 individuals reveal 18 new loci associated with body mass index. Nat Genet. 2010, 42 (11): 937-948. 10.1038/ng.686.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1755-8794/6/14/prepub
Dr. Ahmed Kissebah, the senior author of this work, was fully involved in the preparation of this report up to his death on 17 May 2012. We are saddened by his loss and deeply respectful of his long and distinguished career in obesity research.
We acknowledge the contribution of Roland James (Medical College of Wisconsin), the project manager of TOPS Center for Obesity and Metabolic Research, for his extensive effort in recruitment, phenotyping and database management; Jacqueline Marks (Medical College of Wisconsin) for her supervision of all biochemical procedures; Melena Zelembaba and Regina Cole (Medical College of Wisconsin) for their technical assistance in the Affymetrix chips genotyping procedures; and Ruth Gielow (TOPS Club, Inc.) and her assistants in the recruitment process. Finally, we are most grateful to all the members of TOPS Club, Inc. and their families who volunteered for this study. This work is supported by grants from the NIH (DK071895-03 and DK65598-01) to A.H.K and by TOPS (Take Off Pounds Sensibly) Club, Inc.
The authors declare that they have no competing interests.
Y.Z. analyzed data, contributed to discussion, wrote and revised the manuscript. J.W.K. performed statistical genetic analyses, contributed to discussion and co-wrote the manuscript. M.O. provided laboratory assistance in genotyping, contributed to discussion and data analysis. O.A. assisted in the phenotyping procedures, contributed to discussion and revision of the manuscript. D.C. assisted in data analysis and contributed to revision of the manuscript. R.M.A. assisted in the phenotyping procedures. U.B. facilitated the genotyping procedures. J.E.C. and M.A.C. did cytokine and transcriptome procedures. T.D.D. prepared the genotyping raw data for genetic analysis procedures. A.C. contributed to the original linkage analysis. H.H.H.G. assisted in the transcriptome data analysis. D.L.R. did the lipoprotein particle sizing. J.B. supervised and directed the overall statistical genetic analyses and contributed to discussion. A.H.K designed the study, analyzed data, contributed to discussion and revised the manuscript. All authors read and approved the final manuscript.