Genome-Wide Association Analysis Identifies TYW3
Total Page:16
File Type:pdf, Size:1020Kb
Human Molecular Genetics, 2012, Vol. 21, No. 21 4774–4780 doi:10.1093/hmg/dds300 Advance Access published on July 26, 2012 Genome-wide association analysis identifies TYW3/CRYZ and NDST4 loci associated with circulating resistin levels Qibin Qi1,{, Claudia Menzaghi3,{, Shelly Smith4, Liming Liang2, Nathalie de Rekeneire7, Melissa E. Garcia6, Kurt K. Lohman4, Iva Miljkovic8, Elsa S. Strotmeyer8, Steve R. Cummings9, Alka M. Kanaya9, Frances A. Tylavsky10, Suzanne Satterfield10, Jingzhong Ding5, Eric B. Rimm1,2,11, Vincenzo Trischitta3,12, Frank B. Hu1,2,11, Yongmei Liu4,∗,{ and Lu Qi1,11,∗,{ 1Department of Nutrition and 2Department of Epidemiology, Harvard School of Public Health, Boston, MA, USA, 3Research Unit of Diabetes and Endocrine Diseases, IRCCS ‘Casa Sollievo della Sofferenza’, San Giovanni Rotondo, Italy, 4Division of Public Health Sciences and 5Department of Internal Medicine, Wake Forest School of Medicine, Downloaded from Winston-Salem, NC, USA, 6Laboratory of Epidemiology, Demography and Biometry, National Institute on Aging, Bethesda, MD, USA, 7Department of Geriatrics, Yale University School of Medicine, New Haven, CT, USA, 8Department of Epidemiology, University of Pittsburgh, Pittsburgh, PA, USA, 9Division of General Internal Medicine, University of California, San Francisco, CA, USA, 10Department of Preventive Medicine, University of Tennessee, Memphis, TN, USA, 11Channing Laboratory, Department of Medicine, Harvard Medical School, Boston, MA, USA http://hmg.oxfordjournals.org/ and 12Department of Experimental Medicine, Sapienza University of Rome, Rome, Italy Received February 2, 2012; Revised July 10, 2012; Accepted July 19, 2012 Resistin is a polypeptide hormone that was reported to be associated with insulin resistance, inflammation and risk of type 2 diabetes and cardiovascular disease. We conducted a genome-wide association (GWA) study on circulating resistin levels in individuals of European ancestry drawn from the two independent stud- by guest on August 12, 2015 ies: the Nurses’ Health Study (n 5 1590) and the Health, Aging and Body Composition Study (n 5 1658). Single-nucleotide polymorphisms (SNPs) identified in the GWA analysis were replicated in an independent cohort of Europeans: the Gargano Family Study (n 5 659). We confirmed the association with a previously known locus, the RETN gene (19p13.2), and identified two novel loci near the TYW3/CRYZ gene (1p31) and the NDST4 gene (4q25), associated with resistin levels at a genome-wide significant level, best represented by SNP rs3931020 (P 5 6.37 3 10 –12) and SNP rs13144478 (P 5 6.19 3 10218), respectively. Gene expression quantitative trait loci analyses showed a significant cis association between the SNP rs3931020 and CRYZ gene expression levels (P 5 3.68 3 1027). We also found that both of these two SNPs were significantly asso- ciated with resistin gene (RETN) mRNA levels in white blood cells from 68 subjects with type 2 diabetes (both P 5 0.02). In addition, the resistin-rising allele of the TYW3/CRYZ SNP rs3931020, but not the NDST4 SNP rs13144478, showed a consistent association with increased coronary heart disease risk [odds ratio 5 1.18 (95% CI, 1.03–1.34); P 5 0.01]. Our results suggest that genetic variants in TYW3/CRYZ and NDST4 loci may be involved in the regulation of circulating resistin levels. More studies are needed to verify the associa- tions of the SNP rs13144478 with NDST4 gene expression and resistin-related disease. ∗To whom correspondence should be addressed at: Division of Public Health Sciences, Department of Epidemiology and Prevention, Wake Forest School of Medicine, Winston-Salem, NC 27157, USA. Tel: +1 3367169503; Fax: +1 3367166427; Email: [email protected] (Y.L.); Department of Nutrition, Harvard School of Public Health, 665 Huntington Avenue, Boston, MA 02115, USA. Tel: +1 6174324116; Fax: +1 6174322435; Email: [email protected] (L.Q.). †These authors contributed equally to this work. ‡These authors contributed equally to overseeing this work. # The Author 2012. Published by Oxford University Press. All rights reserved. For Permissions, please email: [email protected] Human Molecular Genetics, 2012, Vol. 21, No. 21 4775 Table 1. Characteristics of the participants in GWAS and replication samples Study n Age (years) Females (%) BMI (kg/m2) Resistin (ng/ml) Log resistin (ng/ml) GWAS NHS (control) 892 57.2 + 6.7 100 25.6 + 4.9 18.2 + 12.2 2.76 + 0.50 NHS (case) 698 57.6 + 6.7 100 30.0 + 5.7 20.9 + 17.0 2.88 + 0.53 Health ABC 1658 73.6 + 2.9 47.1 26.6 + 4.1 19.7 + 9.9 2.90 + 0.38 Replication GFS 659 37.5 + 13.9 62.7 26.2 + 4.7 5.9 + 2.98 0.73 + 0.19 Data are means + SD or otherwise indicated. INTRODUCTION rs13144478 was imputed (MACH r2 ¼ 0.9). The SNPs rs3931020 in the TYW3/CRYZ locus (1p31) and rs13144478 Resistin is a polypeptide hormone that belongs to a family of in the NDST4 locus (4q25–q26) showed directionally consist- cysteine-rich proteins called resistin-like molecules (1). It was ent associations with resistin levels in the GFS, and the originally reported to link obesity to insulin resistance and dia- P-values for the combined results of discovery and replication betes as an adipocytokine in mice (2). Many subsequent ¼ × 212 stages reached genome-wide significance (P 6.4 10 Downloaded from human studies have yielded controversial results in the associ- 2 and 6.2 × 10 18, respectively; Table 2). Regional association ation of circulating resistin levels with insulin resistance, type plots of SNPs in the TYW3/CRYZ and NDST4 loci are dis- 2 diabetes and cardiovascular disease (3–17). In mice, resistin played in Figure 1. In addition, there were eight SNPs in a pre- is preferentially expressed in adipocytes, whereas in humans it viously known locus (12,20–23), the RETN gene (from 7630 is highly expressed in monocyte and macrophage and only to 7650 kb, on chromosome 19p13.2), with association modestly in adipose cells (18,19). Genetic factors may play http://hmg.oxfordjournals.org/ results available in our discovery data. Most of these SNPs important roles in the regulation of circulating resistin were imputed with modest imputation quality (0.2 ≤ MACH levels. It has been estimated that 70% of the variation of cir- r2 ≤ 0.7); and a previously reported SNP rs3745367 (IVS2 culating resistin levels is heritable (20). Several single- 181G/A) (20) showed strongest association with resistin nucleotide polymorphisms (SNPs) in the human resistin gene levels (P ¼ 0.007) (Supplementary Material, Table S1). (RETN) have been associated with circulating resistin levels After combining our present result and previously reported in previous candidate gene studies (12,20–23). A previous data in the GFS (20), the SNP rs3745367 showed a suggestive one-stage genome-wide association (GWA) study failed to 2 association with circulating resistin levels (P ¼ 2.66 × 10 6; indentify loci associated with circulating resisting levels at a Table 2). genome-wide significant level (24). Therefore, we undertook by guest on August 12, 2015 To further examine the positional candidate genes of the a GWA analysis for circulating resistin levels in individuals identified loci, gene expression quantitative trait loci (eQTL) of European ancestry from two discovery cohorts and an inde- analyses were performed in the Gene Expression Omnibus pendent replication study. (GEO) database with expression data from human lympho- blastoid cell lines (26). The SNP rs3931020 showed a signifi- cant cis association with CRYZ gene expression levels (P ¼ RESULTS 3.68 × 1027), and another SNP rs277369 in linkage disequi- The characteristics of participants in each cohort are shown in librium with rs3931020 (r2 ¼ 0.711) showed the most signifi- Table 1. In the discovery stage, we performed three GWA cant association (P ¼ 6.81 × 10224). The SNP rs3931020 also scans for resistin levels in a total of 3248 individuals of showed a significant but weaker cis association with TYW3 European ancestry from the Nurses’ Health Study (NHS) type gene expression (P ¼ 2.7 × 1024), compared with the CRYZ 2 diabetes controls, the NHS type 2 diabetes cases and the gene. Unfortunately, the SNP rs13144478 or correlated Health, Aging and Body Composition Study (Health ABC) SNPs near the NDST4 gene were not available in the gene ex- separately. There was no evidence of systematic bias in the dis- pression database. tribution of P-values for association tests [genomic inflation We also examined the association between the two newly factors (l) were 0.99, 0.99 and 1.02, respectively] (Supplemen- identified SNPs and RETN expression levels in white blood tary Material, Fig. S1). The results from the three GWA scans cells. Consistent with the associations with circulating resistin were combined by a meta-analysis. A graphical summary of levels, the rs3931020 C allele was significantly associated the meta-analysis results is displayed in Supplementary Mater- with higher RETN mRNA levels (P ¼ 0.023 under an additive ial, Figure S2. No associations with resistin levels reached a inheritance model, Fig. 2A), and RETN mRNA was signifi- genome-wide significance level (P ¼ 5 × 1028). cantly higher in carriers of the SNP rs13144478 T allele We took four top SNPs representing independent loci than that in non-T allele carriers (P ¼ 0.022 under a dominant having a stage 1 P , 5 × 1025 for further replication in a rep- inheritance model since only one homozygous individual for lication sample of the Gargano Family Study (GFS, n ¼ 659) the T allele was observed, Fig. 2B). Adjustments for age (Table 2). SNPs were considered to be independent if the pair- and sex did not change the observed association between the wise linkage disequilibrium r2 , 0.1 and if they were apart at SNPs and RETN mRNA expression. least 1 Mb from one another (25).