- Research article
- Open Access
Allelic variants in the PHTF1-PTPN22, C12orf30 and CD226 regions as candidate susceptibility factors for the type 1 diabetes in the Estonian population
BMC Medical Genetics volume 11, Article number: 11 (2010)
Type 1 diabetes is a multifactorial disease with a strong genetic component. The aim of the study was to assess the impact of single nucleotide polymorphisms (SNPs) in several genes as susceptible markers in the risk of type 1 diabetes in the Estonian population. Methods: The rs6679677 (1p13), rs17696736 (12q24) and rs763361 (18q22) were genotyped in a total of 230 controls and 154 type 1 diabetes patients of Estonian origin.
The rs6679677 A (OR = 2.13, 95%CI = 1.48-3.08, p = 0.00001), rs17696736 G (OR = 1.53, 95%CI = 1.14-2.04, p = 0.0046) and rs763361 T (OR = 1.48, 95%CI = 1.11-1.98, p = 0.0084) alleles were associated with risk of type 1 diabetes.
The current study supports the rs6679677 (PHTF1-PTPN22), rs17696736 (C12orf30) and rs763361 (CD226) SNPs as susceptibility factors for type 1 diabetes outside the major histocompatibility region (MHC) region. The full study had 80% or above to detect an odds ratio of 1.8 under the assumption of an additive model at type 1 error rate, α = 0.05.
The genetic predisposition to type 1 diabetes is strongly linked to MHC region  but non - MHC locus associated genes are also considered to predispose in the risk of type 1 diabetes. Evidence strongly supports the insulin gene , the CTLA-4 locus  and the PTPN22 gene  as important candidates outside the MHC region. Meanwhile, the established genetic associations with type 1 diabetes explain roughly half of the genetic risk for type 1 diabetes, indicating that other loci exist. In favor of it, genome wide association studies have recently reported new chromosomal regions and novel candidates in the risk of type 1 diabetes [5, 6]. Of note, within the CD226 gene (18q22) the rs763361 SNP was found to be associated with type 1 diabetes . The CD226 gene is expressed on the majority of immune cells including natural killer (NK) cells and T cells mediating their activation and differentiation . Moreover, the rs17696736 and the rs6679677 in the chromosomal regions of 12q24 and 1q13 respectively, were also reported to be associated with type 1 diabetes . In particular, the rs17696736 SNP (C12orf30 gene) is located within a large linkage disequilibrium (LD) block that contains several genes of possible function relevance to type 1 diabetes . The rs6679677 SNP is an intergenic SNP located downstream of the round spermatid basic protein 1 (RSBN1) and in the promoter site of the putative homeodomain transcription factor 1 (PHTF1) gene. The complexity of type 1 diabetes, the heterogeneity of allele frequencies across populations and the difficulties to localize the true causative variants within a LD block recommends the replication analysis of reported SNPs in order to establish their impact on individual ethnic groups. In the current study we sought to replicate recently reported associations from the genome wide association (GWA) studies in order to evaluate their impact as possible candidates in the risk of type 1 diabetes in Estonian population.
An ethnically homogenous population of Estonian origin including 230 controls (mean age 45.9 ± 14.5 years, 139 females) and 154 type 1 diabetes patients (mean age at diagnosis 22.0 ± 14.3 years, 77 females) were enrolled into a case-control study. Estonia is situated in the north of Europe, on the eastern shores of the Baltic Sea and Estonians belong to the Finno-Ugric ethnic group. Of note, our population cohort is originated from Saaremaa County in Estonia. Saaremaa is the biggest island of Estonia where the population has been relatively permanent over generations compared to the major mainland regions in Estonia. All patients with type 1 diabetes were diagnosed according to the criteria established by the Expert Committee on the Diagnosis and Classification of Diabetes Mellitus . The Ethics Committee of the University of Tartu approved the study, and informed consent was obtained from all participants.
The genotyping of rs6679677 (PHTF1-PTPN22), rs17696736 (C12orf30) and rs763361 (CD226) SNPs was carried out using the SNP genotyping assay (Applied Biosystems, Foster City, CA). PCR reactions were run on the ABI 7000 instrument (Applied Biosystems) using the following cycling parameters: after the first step at 95°C for 10 min, 40 cycles of denaturation at 92°C for 15 s and extension at 60°C for 1 min and the genotypes were analysed using the allelic discrimination function of the instrument. A 100% genotyping success rate was obtained for all studied SNPs.
The R (version 2.4.1) statistical computer program (The R Foundation for Statistical Computing, Boston, MA, http://www.r-project.org) was employed to obtain odds ratio (OR) values as well as Wald's confidence intervals (CI) for alleles and genotypes. Hardy-Weinberg equilibrium (HWE) was performed using the SHEsis software . Permutation analysis was performed using the Haploview 3.2 version software . 10000 permutations were carried out to estimate the significance of the results, correcting for the three loci tested. Power calculations were performed using the PS software . A p value of < 0.05 was considered significant for all analyses.
The rs6679677, rs17696736 and rs763361 SNPs were genotyped in all studied subjects and the analysis of genotypes distribution showed no significant deviation from the HWE in type 1 diabetes or in the control group (p > 0.05). A significant higher frequency of the rs6679677 AA (OR = 4.23, 95%CI = 1.56-11.52, p = 0.0047), rs17696736 GG (OR = 2.32, 95%CI = 1.25-4.28, p = 0.0074) and rs763361 TT (OR = 2.29, 95%CI = 1.25-4.18, p = 0.0071) genotype was observed in the group of type 1 diabetes patients when compared to controls (table 1). Moreover, the distribution frequency of the rs6679677 A (OR = 2.13, 95%CI = 1.48-3.08, p = 0.00001), rs17696736 G (OR = 1.53, 95%CI = 1.14-2.04, p = 0.0046) and rs763361 T (OR = 1.48, 95%CI = 1.11-1.98, p = 0.0084) allele was significantly higher in type 1 diabetes when compared with the control group (table 1). Power calculation analysis showed that given the frequencies and genotype odds ratio observed in our population for the homozygotes (high risk alleles), our case control study had 85%, 89% and 90% power (at α = 0.05), to detect significant differences for the rs6679677, rs1769676 and rs763361 markers, respectively. Of note, all observed type 1 diabetes associations with the rs6679677 A, rs17696736 G and rs763361 T alleles remained significant after permutation analysis (10000 permutations, giving p = 0.0004, p = 0.016 and p = 0.034, respectively).
In the present study we attempted to assess the impact of novel reported findings from genome wide type 1 diabetes association studies in a population of Estonian origin. In favor of it, the association of the rs6679677 SNP in susceptibility of type 1 diabetes has been confirmed; thus making it an additional candidate within the 1q13 region. Notably, the rs6679677 SNP is located in the same haplotype block with the rs2476601 SNP of the PTPN22 gene. Having previously demonstrated a significant association of T allele of rs2476601 SNP in type 1 diabetes , we sought to assess the hypothesis that the rs2476601 SNP might not be the sole causal variant in the region. However, according to our previous findings  and based on the genetic data alone, we were unable to distinguish between the two studied SNPs since they are in perfect LD (R2 = 1.00, D' = 1.00). This is in correlation with the results of recent studies in which an association of the rs6679677 SNP with type 1 diabetes has also been reported and the effects of the two SNPs could not be distinguished due to tight LD [14, 15]. Moreover, through using the Genomatix web tool for SNP analysis http://www.genomatix.de, as Qu et al. demonstrated , we revealed that the substitution of cytosine (C) to adenine (A) in the rs6679677 locus predicts the loss of two (CAAT, CDP) and the generation of one (Sox-5) transcription binding sites. In addition we agree with Qu et al.  that the region around PTPN22 may harbor additional disease related alleles not yet found.
Furthermore, we confirmed that the rs763361 SNP in the CD226 gene was associated with susceptibility to type 1 diabetes. The rs763361 SNP is predicted to be a functional SNP, residing in the final 3' of the exon 7 and alters the exon splicing silencer (ESS) sequence . The disease-susceptible minor allele A, which encodes Ser307, maintains the wild type ESS sequence and compared to the resistant Gly307 allele would have a higher silencer function that reduces the productive splicing of exons 6 and 7. Notably, the CD226 molecule is expressed on a variety of hematopoietic cells and interaction with its ligands poliovirus receptor (PVR, CD155) and Nectin-2 (CD112) results in a variety of cellular responses of innate and adaptive immunity . Recently, new studies have provided additional evidence and support the role of rs763361 SNP in the risk of autoimmune thyroid disease , rheumatoid arthritis, multiple sclerosis [18, 19] and Wegener's granulomatosis .
Moreover, we showed that the rs17696736 SNP in the C12orf30 gene also contribute to susceptibility of type 1 diabetes. Interestingly, the C120rf30 gene maps into a region of extensive linkage disequilibrium including several genes that represent functional candidates of type 1 diabetes-susceptibility because of their presumed roles in immune signaling. Remarkably, the SH2B3 (SH2B adaptor protein 3), TRAFD1 (TRAF-type zinc finger domain containing 1) and PTPN11 (protein tyrosine phosphatase, non-receptor type 11) map on the same region with the C12orf30 gene. Therefore, for this region in particular, extensive resequencing, further genotyping and targeted functional studies will be essential steps in identifying which gene, or genes, are causal.
One common issue in the study of complex diseases, including type 1 diabetes, is the limited sample size, resulting in adequate power to detect association or the risk of false positive results. Typically, the full study had 80% power or above to detect an odds ratio of 1.8 under the assumption of an additive model at type 1 error rate, α = 0.05. Further, take in consideration the sample size of our study, permutation analysis has been applied in order to evaluate the impact of our findings. Of note, all the reported associations remained significant after applied permutation analysis.
Noticeably, further genotyping, resequencing and fine mapping could assist in clarifying the true causative genes within a strong LD block and could reveal new susceptible allelic variants for type 1 diabetes. The potential of the latter approaches to detect new associations have recently been demonstrated by Nejentsev et al . Moreover, we should also take into account that due to allelic heterogeneity, there is a variation in risk allele frequencies between populations, which may account for differences in disease prevalence between ethnic groups.
In the current study, we have demonstrated a significant association between recently reported SNPs and type 1 diabetes in Estonian population. To confirm our findings, further associated studies are needed to elucidate the functional role of these SNPs as likely mediators for autoimmunity and as potential therapeutic targets for antidiabetic drugs. The probability that the newly reported loci with type 1 diabetes may predispose other autoimmune diseases should also be examined.
Horn GT, Bugawan TL, Long CM, Erlich HA: Allelic sequence variation of the HLA-DQ loci: relationship to serology and to insulin-dependent diabetes susceptibility. Proc Natl Acad Sci USA. 1988, 85: 6012-6016. 10.1073/pnas.85.16.6012.
Bennett ST, Lucassen AM, Gough SCL, Powell EE, Undlien DE, Pritchard LE, Merriman ME, Kawaguchi Y, Dronsfield MJ, Pociot F, Nerup J, Bouzekri N, Cambon-Thomsen A, Rønningen KS, Barnett AH, Bain SC, Todd JA: Susceptibility to human type 1 diabetes at IDDM2 is determined by tandem repeat variation at the insulin gene minisatellite locus. Nat Genet. 1995, 9: 284-292. 10.1038/ng0395-284.
Ueda H, Howson JM, Esposito L, Heward J, Snook H, Chamberlain G, Rainbow DB, Hunter KM, Smith AN, Di Genova G, Herr MH, Dahlman I, Payne F, Smyth D, Lowe C, Twells RC, Howlett S, Healy B, Nutland S, Rance HE, Everett V, Smink LJ, Lam AC, Cordell HJ, Walker NM, Bordin C, Hulme J, Motzo C, Cucca F, Hess JF, Metzker ML, Rogers J, Gregory S, Allahabadia A, Nithiyananthan R, Tuomilehto-Wolf E, Tuomilehto J, Bingley P, Gillespie KM, Undlien DE, Rønningen KS, Guja C, Ionescu-Tîrgovişte C, Savage DA, Maxwell AP, Carson DJ, Patterson CC, Franklyn JA, Clayton DG, Peterson LB, Wicker LS, Todd JA, Gough SC: Association of the T-cell regulatory gene CTLA4 with susceptibility to autoimmune disease. Nature. 2003, 423: 506-511. 10.1038/nature01621.
Bottini N, Musumeci L, Alonso A, Rahmouni S, Nika K, Rostamkhani M, MacMurray J, Meloni GF, Lucarelli P, Pellecchia M, Eisenbarth GS, Comings D, Mustelin T: A functional variant of lymphoid tyrosine phosphatase is associated with type I diabetes. Nat Genet. 2004, 36: 337-338. 10.1038/ng1323.
Todd JA, Walker NM, Cooper JD, Smyth DJ, Downes K, Plagnol V, Bailey R, Nejentsev S, Field SF, Payne F, Lowe CE, Szeszko JS, Hafler JP, Zeitels L, Yang JH, Vella A, Nutland S, Stevens HE, Schuilenburg H, Coleman G, Maisuria M, Meadows W, Smink LJ, Healy B, Burren OS, Lam AA, Ovington NR, Allen J, Adlem E, Leung HT, Wallace C, Howson JM, Guja C, Ionescu-Tîrgovişte C, Genetics of Type 1 Diabetes in Finland, Simmonds MJ, Heward JM, Gough SC, Wellcome Trust Case Control Consortium, Dunger DB, Wicker LS, Clayton DG: Robust associations of four new chromosome regions from genome-wide analyses of type 1 diabetes. Nat Genet. 2007, 39: 857-864. 10.1038/ng2068.
Barrett JC, Clayton DG, Concannon P, Akolkar B, Cooper JD, Erlich HA, Julier C, Morahan G, Nerup J, Nierras C, Plagnol V, Pociot F, Schuilenburg H, Smyth DJ, Stevens H, Todd JA, Walker NM, Rich SS, The Type 1 Diabetes Genetics Consortium: Genome-wide association study and meta-analysis find that over 40 loci affect risk of type 1 diabetes. Nat Genet. 2009,
Shibuya A, Campbell D, Hannum C, Yssel H, Franz-Bacon K, McClanahan T, Kitamura T, Nicholl J, Sutherland G, Lanier L: DNAM-1, a novel adhesion molecule involved in the cytolytic function of T lymphocytes. Immunity. 1996, 4: 573-581. 10.1016/S1074-7613(00)70060-4.
Wellcome Trust Case Control Consortium: Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls. Nature. 2008, 447: 661-678. 10.1038/nature05911.
Expert Committee on the diagnosis and classification of diabetes mellitus: Report of the Expert Committee on the diagnosis and classification of diabetes mellitus. Diabetes Care. 2003, 26 (Suppl 1): 5-20.
Shi YY, He L: SHEsis, a powerful software platform for analyses of linkage disequilibrium, haplotype construction, and genetic association at polymorphism loci. Cell Res. 2005, 15: 97-98. 10.1038/sj.cr.7290286.
Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005, 21: 263-265. 10.1093/bioinformatics/bth457.
Purcell S, Cherny SS, Sham PC: Genetic power calculator: design of linkage and association genetic mapping studies of complex traits. Bioinformatics. 2003, 19: 149-150. 10.1093/bioinformatics/19.1.149.
Douroudis K, Prans E, Haller K, Nemvalts V, Rajasalu T, Tillmann V, Kisand K, Uibo R: Protein tyrosine phosphatase non-receptor type 22 gene variants at position 1858 are associated with type 1 and type 2 diabetes in Estonian population. Tissue Antigens. 2008, 72: 425-430. 10.1111/j.1399-0039.2008.01115.x.
Smyth DJ, Cooper JD, Howson JM, Walker NM, Plagnol V, Stevens H, Clayton DG, Todd JA: PTPN22 Trp620 explains the association of chromosome 1p13 with type 1 diabetes and shows a statistical interaction with HLA class II genotypes. Diabetes. 2008, 57: 1739-1737. 10.2337/db07-1131.
Qu H, Tessier M, Hudson T, Polychronakos C: Confirmation of the association of the R620W polymorphism in the protein tyrosine phosphatase PTPN22 with type 1 diabetes in a family based study. J Med Genet. 2005, 42: 266-270. 10.1136/jmg.2004.026971.
Bottino C, Castriconi R, Pende D, Rivera P, Nanni M, Carnemolla B, Cantoni C, Grassi J, Marcenaro S, Reymond n, Vitale m, Moretta L, Lopez M, Moretta A: Identification of PVR (CD155) and Nectin-2 (CD112) as cell surface ligands for the human DNAM-1 (CD226) activating molecule. J Exp Med. 2003, 198: 557-567. 10.1084/jem.20030788.
Hafler JP, Maier LM, Cooper JD, Plagnol V, Hinks A, Simmonds MJ, Stevens HE, Walker NM, Healy B, Howson JM, Maisuria M, Duley S, Coleman G, Gough SC, International Multiple Sclerosis Genetics Consortium (IMSGC), Worthington J, Kuchroo VK, Wicker LS, Todd JA: CD226 Gly307Ser association with multiple autoimmune diseases. Genes Immun. 2009, 10: 5-10. 10.1038/gene.2008.82.
International Multiple Sclerosis Genetics Consortium (IMSGC): The expanding genetic overlap between multiple sclerosis and type 1 diabetes. Genes Immun. 2009, 10: 11-14. 10.1038/gene.2008.83.
Wieczorek S, Hoffjan S, Chan A, Rey L, Harper L, Fricke H, Holle JU, Gross WL, Epplen JT, Lamprecht P: Novel association of the CD226 (DNAM-1) Gly307Ser polymorphism in Wegener's granulomatosis and confirmation for multiple sclerosis in German patients. Genes Immun. 2009, 10: 591-595. 10.1038/gene.2009.44.
Nejentsev S, Walker N, Riches D, Egholm M, Todd JA: Rare variants of IFIH1, a gene implicated in antiviral responses, protect against type 1 diabetes. Science. 2009, 324: 387-389. 10.1126/science.1167728.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2350/11/11/prepub
We thank Dr. M. Möls from the Institute of Mathematical Statistics, University of Tartu for advice in statistical analysis. This study was supported by the Estonian Science Foundation grant no. 7749, partly by INNOVE SA and by the European Union through the European Regional Development Fund.
The authors declare that they have no competing interests.
KD performed the genotyping, the statistical analysis, interpreted the data and drafted the manuscript. KK, TR and VN participated in materials' and data collection. RU participated in the initiation and coordination of the study. All authors read and approved the final manuscript.
About this article
Cite this article
Douroudis, K., Kisand, K., Nemvalts, V. et al. Allelic variants in the PHTF1-PTPN22, C12orf30 and CD226 regions as candidate susceptibility factors for the type 1 diabetes in the Estonian population. BMC Med Genet 11, 11 (2010). https://doi.org/10.1186/1471-2350-11-11
- CD226 Gene
- Linkage Disequilibrium Block
- Permutation Analysis
- PTPN22 Gene
- C12orf30 Gene