Association of a genetic variant of the ZPR1 zinc finger gene with type 2 diabetes mellitus
- Authors:
- Published online on: November 7, 2014 https://doi.org/10.3892/br.2014.379
- Pages: 88-92
Abstract
Introduction
Type 2 diabetes mellitus (DM) is a major public health issue that affects over one billion people worldwide (1). Type 2 DM is a complex disease that involves genetic and environmental factors and their interactions (2). Due to the high prevalence of type 2 DM, identifying the genes or genetic loci associated with the risk or protection of type 2 DM is important for understanding the mechanisms underlying the disease and for benefiting the patients with personalized prevention and treatment programs.
Genome-wide association studies (GWASs) and subsequent meta-analyses have identified >56 susceptibility loci for type 2 DM (3–11). However, these susceptibility loci have been identified predominantly in Caucasian populations. Differences in allele frequencies and the effect of size among different ethnicity groups yielded the discovery of new loci in different populations (12). Although several single-nucleotide polymorphisms (SNPs) have been identified as susceptibility loci for type 2 DM in Japanese individuals (10,11), the genes that confer susceptibility to this condition remain to be identified definitively.
Previous GWASs have identified various loci and genes that confer susceptibility to coronary heart disease (CHD) for Caucasian populations (13,14). As type 2 DM is an important risk factor for CHD, we hypothesized that certain polymorphisms may contribute to the genetic susceptibility to CHD through affecting the susceptibility to type 2 DM. The purpose of the present study was to examine a possible association of type 2 DM in Japanese individuals with 29 SNPs identified as susceptibility loci for CHD by meta-analyses of the GWASs.
Subjects and methods
Study population
The study population comprised of 3,757 Japanese individuals (1,444 subjects with type 2 DM and 2,313 controls) who either visited outpatient clinics or were admitted to the participating hospitals (Gifu Prefectural General Medical Center, Gifu; Gifu Prefectural Tajimi Hospital, Tajimi; Japanese Red Cross Nagoya First Hospital, Nagoya; Inabe General Hospital, Inabe; Hirosaki University Hospital, Reimeikyo Rehabilitation Hospital and Hirosaki Stroke Center, Hirosaki, Japan) between 2002 and 2012 due to various symptoms or for an annual health checkup. Written informed consent was obtained from all the participants and the Institutional Review Board of each participating hospital approved the study.
Type 2 DM is defined according to the criteria of the World Health Organization, as described previously (15,16). Subjects with type 2 DM had a fasting plasma glucose level of ≥6.93 mmol/l (126 mg/dl), a blood glycosylated hemoglobin (hemoglobin A1c) content of ≥6.5%, or were taking anti-diabetic medication. Individuals with type 1 DM, maturity onset diabetes of the young, DM associated with mitochondrial diseases or single-gene disorders, pancreatic diseases or other metabolic or endocrinological diseases were excluded from the study. Individuals on medication that may cause secondary DM were also excluded. The control individuals had a fasting plasma glucose level of <6.05 mmol/l (110 mg/dl), a blood hemoglobin A1c content of <6.2% and had no history of DM or of receiving anti-diabetic medication.
Selection and genotyping of polymorphisms
SNPs that were recently identified as susceptibility loci for CHD in Caucasian populations were searched for by meta-analyses of GWASs (13,14). These SNPs were examined with the dbSNP database (National Center for Biotechnology Information; http://www.ncbi.nlm.nih.gov/SNP/) to find SNPs with a minor allele frequency of >0.015 in a Japanese population. Finally, 29 SNPs (data not shown) were selected and the association with type 2 DM was examined. Wild-type and variant alleles of the SNPs were determined from the original sources.
Venous blood (7 ml) was collected into tubes containing 50 mmol/l ethylenediaminetetraacetic acid (disodium salt) and genomic DNA was isolated with a kit (Genomix; Talent, Trieste, Italy). Genotypes of the 29 SNPs were determined at G&G Science (Fukushima, Japan) by a method that combines polymerase chain reaction and sequence-specific oligonucleotide probes with suspension array technology (Luminex Corporation, Austin, TX, USA). The overall call rate of genotyping of 29 SNPs was 99%. The detailed genotyping methodology was performed as described previously (17).
Statistical analysis
The χ2 test was used to compare the categorical variables, whereas the Mann-Whitney U test was used for analysis of the quantitative data. Allele frequencies of each SNP were compared between subjects with type 2 DM and controls by the χ2 test. A false discovery rate (FDR) was calculated to compensate for multiple comparisons of genotypes, and FDR≤0.05 was considered to indicate a statistical significance for association. Multivariable logistic regression analysis was performed with type 2 DM as a dependent variable and age, gender (0, women; 1, men), body mass index (BMI) and the genotype of SNP as independent variables. The SNP was assessed according to dominant (the combined group of heterozygotes and variant homozygotes verses wild-type homozygotes), recessive (variant homozygotes verses the combined group of wild-type homozygotes and heterozygotes) and two additive [additive 1 (heterozygotes verses wild-type homozygotes) and additive 2 (variant homozygotes verses wild-type homozygotes)] genetic models. As fasting plasma glucose level and blood hemoglobin A1c content were not normally distributed (P<0.01 by the Kolmogorov-Smirnov Lilliefors test), these parameters were compared among genotypes by the non-parametric Kruskal-Wallis test. Statistical analysis was performed with JMP version 11 and JMP Genomics version 6.0 software (SAS Institute, Cary, NC, USA).
Results
Clinical characteristics of the study subjects
The clinical characteristics of the study subjects are shown in Table I. Age, the frequency of males, BMI, the prevalence of smoking, myocardial infarction, dyslipidemia and hypertension, as well as serum concentrations of triglycerides and creatinine, were higher, whereas the serum concentrations of high-density lipoprotein (HDL) cholesterol were lower in subjects with type 2 DM compared to controls.
Associations of SNPs to type 2 DM
Allele frequencies were compared between subjects with type 2 DM and controls by the χ2 test and five SNPs with P<0.05 are shown in Table II. Among these SNPs, rs964184 (C→G) of the ZPR1 zinc finger gene (ZPR1) was significantly (FDR≤0.05) associated with the prevalence of type 2 DM. The genotype distributions of five SNPs were in Hardy-Weinberg equilibrium (P>0.05) among subjects with type 2 DM and controls.
![]() | Table IIComparison of the single-nucleotide polymorphism (P<0.05) allele frequencies by the χ2 test between subjects with type 2 diabetes mellitus (DM) and controls. |
Multivariable logistic regression analysis with adjustment for age, gender and BMI revealed that rs964184 of ZPR1 was significantly associated with type 2 DM in the dominant and additive 1 and 2 models, with the minor G allele representing a risk factor for this condition (Table III). As hypertriglyceridemia is an important risk factor for type 2 DM, additional multivariable logistic regression analysis was performed with adjustment for serum triglycerides concentrations or hypertriglyridemia (serum concentration of triglycerides ≥1.65 mmol/l or taking anti-dyslipidemic medication) in addition to age, gender and BMI (Table III). rs964184 was also significantly associated with type 2 DM in the dominant and additive 1 models in this analysis.
![]() | Table IIIMultivariable logistic regression analysis of rs964184 of ZPR1 zinc finger gene and type 2 diabetes mellitus with additional adjustments to age, gender and BMI. |
Associations of rs964184 to fasting plasma glucose level and blood hemoglobin A1c content
Finally, the associations of rs964184 genotypes to fasting plasma glucose level and blood hemoglobin A1c content were examined by the Kruskal-Wallis test (Table IV). rs964184 was significantly associated with the two parameters and the G allele was associated with the increases in fasting plasma glucose level and in blood hemoglobin A1c content.
![]() | Table IVAssociation of rs964184 of ZPR1 zinc finger gene to fasting plasma glucose level and blood hemoglobin A1c content as determined by the Kruskal-Wallis test. |
Discussion
The associations of 29 SNPs identified as susceptibility loci for CHD by meta-analyses of GWASs to type 2 DM were examined and it was observed that rs964184 of ZPR1 was significantly associated with type 2 DM in Japanese individuals. The prevalence of type 2 DM, fasting plasma glucose level and blood hemoglobin A1c content were increased by 18.0, 6.7 and 7.4%, respectively, for individuals with the GG genotype of rs964184 compared to those with the CC genotype.
rs964184 is located in the intron region of ZPR1 at chromosome 11q23.3. ZPR1 is an essential regulatory protein for cell proliferation and signal transduction and may have multiple physiological functions (18,19). The most relevant transcription factor that binds to the promoter region of ZPR1 is peroxisome proliferator-activated receptor γ, which plays an important role in insulin sensitivity and obesity (20,21). The promoter region of ZPR1 is also bound by hepatocyte nuclear factor 4α, which activates a variety of genes involved in glucose, fatty acid and cholesterol metabolism (22).
ZPR1 is located ~1.6 kb upstream of the APOA5-A4-C3-A1 gene complex. Previous studies have shown that several polymorphisms in or near APOA5 are significantly associated with serum triglycerides concentrations (23–26). rs964184 of ZPR1 has been associated with serum triglycerides and this may be attributable to linkage disequilibrium with functional SNPs in APOA5, which influence metabolism of chylomicrons, very-low-density lipoprotein and HDL (27). As an increase in serum triglycerides concentration is an important risk factor for type 2 DM (28), a multivariable logistic regression analysis was performed with adjustment for serum triglycerides levels or hypertriglyceridemia in addition to age, gender and BMI. There was a significant association of rs964184 with type 2 DM in this analysis, indicating that the association was independent, at least in part, of serum triglycerides levels in the study. The previous GWASs suggested that APOA5 polymorphisms may also play an important role in the development of type 2 DM (29,30). A subgroup analysis by ethnicity of a meta-analysis revealed a significant association of the −1131T→C polymorphism of APOA5 with type 2 DM in Asian populations (31). This observation may support the hypothesis that rs964184 of ZPR1 is associated with type 2 DM through the interaction with APOA5 in Japanese individuals.
Although the contribution of rs964184 to the increased susceptibility to type 2 DM was examined in several GWASs mainly with Caucasians, the significant association was not detected (32,33). The reason for the discrepancy between the previous studies and the present results remains unclear. The variations in the minor G allele frequencies due to the ethnic differences may be, at least in part, responsible for this discrepancy. The frequencies of the CG and GG genotypes of rs9645184 were 23.7 and 2.4%, respectively, in Caucasian populations (32), whereas in the present study population they were 39.1 and 7.2%, respectively. The G allele of rs964184 was therefore higher in the present population (26.8%) compared to the Caucasian population (13–14%) (23,33,34). In addition, the prevalence of type 2 DM in the present population was 38.2%, which was more than twice that reported previously (32). The higher frequency of the G allele and the higher prevalence of type 2 DM in the present population compared to those in the previous studies (33,34) may increase the statistical power to detect the association of rs964184 with type 2 DM.
In conclusion, the results indicate that rs964184 (C→G) of ZPR1 may be a susceptibility locus for type 2 DM in Japanese individuals. Validation of these findings is required in other independent subject panels or ethnic groups.
Acknowledgements
The present study was supported by a Collaborative Research Grant from the Gifu Prefectural General Medical Center (no. H24-26 to Y.Y.) and a Grant-in-Aid for Scientific Research from the Ministry of Education, Culture, Sports, Science and Technology of Japan (no. 24590746 to Y.Y.).
References
Unwin N, Gan D and Whiting D: The IDF Diabetes Atlas: proving evidence, raising awareness and promoting action. Diabetes Res Clin Pract. 87:2–3. 2010. View Article : Google Scholar | |
Hu FB: Globalization of diabetes: the role of diet, lifestyle, and genes. Diabetes Care. 34:1249–1257. 2011. View Article : Google Scholar : PubMed/NCBI | |
McCarthy MI: Genomics, type 2 diabetes, and obesity. N Engl J Med. 363:2339–2350. 2010. View Article : Google Scholar : PubMed/NCBI | |
Shu XO, Long J, Cai Q, et al: Identification of new genetic risk variants for type 2 diabetes. PLoS Genet. 6:e10011272010. View Article : Google Scholar : PubMed/NCBI | |
Kooner JS, Saleheen D, Sim X, et al: Genome-wide association study in individuals of South Asian ancestry identifies six new type 2 diabetes susceptibility loci. Nat Genet. 43:984–989. 2011. View Article : Google Scholar : PubMed/NCBI | |
Cho YS, Chen CH, Hu C, et al: Meta-analysis of genome-wide association studies identifies eight new loci for type 2 diabetes in East Asians. Nat Genet. 44:67–72. 2011. View Article : Google Scholar : PubMed/NCBI | |
Prokopenko I, McCarthy MI and Lindgren CM: Type 2 diabetes: new genes, new understanding. Trends Genet. 24:613–621. 2008. View Article : Google Scholar : PubMed/NCBI | |
Tabassum R, Chauhan G, Dwivedi OP, et al: Genome-wide association study for type 2 diabetes in Indians identifies a new susceptibility locus at 2q21. Diabetes. 62:977–986. 2013. View Article : Google Scholar : PubMed/NCBI | |
Sladek R, Rocheteau G, Rung J, et al: A genome-wide association study identifies novel risk loci for type 2 diabetes. Nature. 445:881–885. 2007. View Article : Google Scholar : PubMed/NCBI | |
Yasuda K, Miyake K, Horikawa Y, et al: Variants in KCNQ1 are associated with susceptibility to type 2 diabetes mellitus. Nat Genet. 40:1092–1097. 2008. View Article : Google Scholar : PubMed/NCBI | |
Yamauchi T, Hara K, Maeda S, et al: A genome-wide association study in the Japanese population identifies susceptibility loci for type 2 diabetes at UBE2E2 and C2CD4A-C2CD4B. Nat Genet. 42:864–868. 2010. View Article : Google Scholar : PubMed/NCBI | |
McCarthy MI: The importance of global studies of the genetics of type 2 diabetes. Diabetes Metab J. 35:91–100. 2011. View Article : Google Scholar : PubMed/NCBI | |
Schunkert H, Konig IR, Kathiresan S, et al: Large-scale association analysis identifies 13 new susceptibility loci for coronary artery disease. Nat Genet. 43:333–338. 2011. View Article : Google Scholar : PubMed/NCBI | |
CARDIoGRAMplusC4D Consortium, . Deloukas P, Kanoni S, et al: Large-scale association analysis identifies new risk loci for coronary artery disease. Nat Genet. 45:25–33. 2013. View Article : Google Scholar : PubMed/NCBI | |
Expert Committee on the Diagnosis. Classification of Diabetes Mellitus: Report of the expert committee on the diagnosis and classification of diabetes mellitus. Diabetes Care. 26 (Suppl 1):S5–S20. 2003. View Article : Google Scholar : PubMed/NCBI | |
Kuzuya T, Nakagawa S, Satoh J, et al: Report of the Committee on the classification and diagnosis criteria of diabetes mellitus. Diabetes Res Clin Pract. 55:65–85. 2002. View Article : Google Scholar : PubMed/NCBI | |
Itoh Y, Mizuki N, Shimada T, et al: High-throughput DNA typing of HLA-A, -B, -C, and -DRB1 loci by a PCR-SSOP-Luminex method in the Japanese population. Immunogenetics. 57:717–729. 2005. View Article : Google Scholar : PubMed/NCBI | |
Galcheva-Gargova Z, Konstantinov KN, Wu IH, et al: Binding of zinc finger protein ZPR1 to the epidermal growth factor receptor. Science. 272:1797–1802. 1996. View Article : Google Scholar : PubMed/NCBI | |
Galcheva-Gargova Z, Gangwani L, Konstantinov KN, et al: The cytoplasmic zinc finger protein ZPR1 accumulates in the nucleolus of proliferating cells. Mol Biol Cell. 9:2963–2971. 1998. View Article : Google Scholar : PubMed/NCBI | |
Mangelsdorf DJ, Thummel C, Beato M, et al: The nuclear receptor superfamily: the second decade. Cell. 83:835–839. 1995. View Article : Google Scholar : PubMed/NCBI | |
Corton JC, Anderson SP, Stauber A, et al: Central role of peroxisome proliferator-activated receptors in the actions of peroxisome proliferators. Annu Rev Pharmacol Toxicol. 40:491–518. 2000. View Article : Google Scholar : PubMed/NCBI | |
Sladek FM, Zhong WM, Lai E, et al: Liver-enriched transcription factor HNF-4 is a novel member of the steroid hormone receptor superfamily. Genes Dev. 4:2353–2365. 1990. View Article : Google Scholar : PubMed/NCBI | |
Teslovich TM, Musunuru K, Smith AV, et al: Biological, clinical and population relevance of 95 loci for blood lipids. Nature. 466:707–713. 2010. View Article : Google Scholar : PubMed/NCBI | |
Kathiresan S, Willer CJ, Peloso GM, et al: Common variants at 30 loci contribute to polygenic dyslipidemia. Nat Genet. 41:56–65. 2009. View Article : Google Scholar : PubMed/NCBI | |
Willer CJ, Sanna S, Jackson AU, et al: Newly identified loci that influence lipid concentrations and risk of coronary artery disease. Nat Genet. 40:161–169. 2008. View Article : Google Scholar | |
Waterworth DM, Ricketts SL, Song K, et al: Genetic variants influencing circulating lipid levels and risk of coronary artery disease. Arterioscler Thromb Vasc Biol. 30:2264–2276. 2010. View Article : Google Scholar : PubMed/NCBI | |
O'Brien PJ, Alborn WE, Sloan JH, et al: The novel apolipoprotein A5 is present in human serum, is associated with VLDL, HDL, and chylomicrons, and circulates at very low concentrations compared with other apolipoproteins. Clin Chem. 51:351–359. 2005. | |
Bernstein RM, Davis BM, Olefsky JM, et al: Hepatic insulin responsiveness in patients with endogenous hypertriglyceridaemia. Diabetologia. 14:249–253. 1978. View Article : Google Scholar : PubMed/NCBI | |
Talmud PJ, Cooper JA, Hattori H, et al: The apolipoprotein A-V genotype and plasma apolipoprotein A-V and triglyceride levels: prospective risk of type 2 diabetes. Results from the Northwick Park Heart Study II. Diabetologia. 49:2337–2340. 2006. View Article : Google Scholar : PubMed/NCBI | |
Soter MO, Gomes KB, Fernandes AP, et al: −1131T>C and SW19 polymorphisms in APOA5 gene and lipid levels in type 2 diabetic patients. Mol Biol Rep. 39:7541–7548. 2012. View Article : Google Scholar | |
Yin YW, Sun QQ, Wang PJ, et al: Genetic polymorphism of apolipoprotein A5 gene and susceptibility to type 2 diabetes mellitus: a meta-analysis of 15,137 subjects. PLoS One. 9:e891672014. View Article : Google Scholar : PubMed/NCBI | |
Kristiansson K, Perola M, Tikkanen E, et al: Genome-wide screen for metabolic syndrome susceptibility loci reveals strong lipid gene contribution but no evidence for common genetic basis for clustering of metabolic syndrome traits. Circ Cardiovasc Genet. 5:242–249. 2012. View Article : Google Scholar | |
van de Woestijne AP, van der Graaf Y, de Bakker PI, et al: Rs964184 (APOA5-A4-C3-A1) is related to elevated plasma triglyceride levels, but not to an increased risk for vascular events in patient with clinically manifest vascular disease. PLoS One. 9:e1010822014. | |
Johansen CT, Wang J, Lanktree MB, et al: Excess of rare variants in genes identified by genome-wide association study of hypertriglyceridemia. Nat Genet. 42:684–687. 2010. View Article : Google Scholar : PubMed/NCBI |