Original Article

SLC30A8 Gene rs13266634 C/T Polymorphism in Children with Type 1 Diabetes in Tamil Nadu, India

10.4274/jcrpe.galenos.2018.2018.0195

  • Ramasamy Thirunavukkarasu
  • Arthur Joseph Asirvatham
  • Ayyappan Chitra
  • Mariakuttikan Jayalakshmi

Received Date: 28.07.2018 Accepted Date: 07.09.2018 J Clin Res Pediatr Endocrinol 2019;11(1):55-60 PMID: 30197307

Objective:

Zinc transporter 8 (ZnT8) is a multi-transmembrane protein situated in the insulin secretory granule of the islets of β-cells and is identified as a novel auto-antigen in type 1 diabetes (T1D). The gene coding for ZnT8, solute carrier family 30 member 8 (SLC30A8) is located on chromosome 8q24.11. This study aimed to identify the association of SLC30A8 rs13266634 C/T gene polymorphism with T1D in a sample of T1D children in Tamil Nadu, India.

Methods:

The family based study was conducted in 121 T1D patients and 214 of their family members as controls. The SLC30A8 gene rs13266634 C/T polymorphism was evaluated by polymerase chain reaction-restriction fragment length polymorphism.

Results:

No significant differences were observed in either allele (odds ratio: 0.92; confidence interval: 0.33-2.58; p=0.88) and genotype (CC: p=0.74; CT: p=0.82; TT: p=0.80) frequencies of rs13266634 C/T between T1D patients and controls. Transmission disequilibrium test has identified over-transmission of mutant T allele from parents to affected children (T: U=9:7) without statistical significance. Metaanalysis on the overall effects of rs13266634 C allele frequency was not different (p=0.10 and Pheterogeneity=0.99) in T1D patients as compared to the controls.

Conclusion:

The present study along with the meta-analysis does not show any substantial association of the rs13266634 C/T polymorphism with T1D development in this population.

Keywords: Type 1 diabetes, auto-antigen, polymorphisms, zinc transporter 8 autoantibody, meta-analysis

What is already known on this topic?

SLC30A8 rs13266634 C/T polymorphism in type 1 diabetes (T1D) patients from four different populations was previously reported. This gene polymorphism is associated with T1D in the German population, but not in Danish, Japanese and British populations.

What this study adds?

To our knowledge, this is the first family-based report addressing SLC30A8 gene polymorphism in South Indian patients. The present study and the meta-analysis show that the rs13266634 C/T polymorphism is not associated with type 1 diabetes in this population.


Introduction

Type 1 diabetes (T1D) is a complex, multifactorial disease caused by the selective destruction of insulin-producing pancreatic β-cells (1,2). The autoimmune destruction of pancreatic β-cells by pathogenic T cells predominately targets a number of well-known β-cell auto-antigens (3). Islet cell auto-antigens identified in T1D are Zinc transporter 8 (ZnT8), glutamic acid decarboxylase 65, tyrosine phosphatase-related molecules-2 and insulin (4). ZnT8 is a multi-transmembrane protein, belonging to the family of zinc transporters, having a role in the transport of zinc ions generated from the cytoplasm to the insulin vesicles and plays a major role in insulin maturation (5). During the process of insulin biosynthesis and secretion, frequent exocytosis of glucose stimulated insulin secretion increase the chance of ZnT8 expression on the β-cell surface (6), which further causes more ZnT8 antigen to be exposed. Once ZnT8 is exposed, it can trigger or exacerbate the production of ZnT8 autoantibodies in genetically susceptible individuals (7). Previous studies have reported autoantibodies to ZnT8 to be highly prevalent among new-onset T1D children and have suggested that they could be a marker for disease risk (8,9,10,11). The cation efflux transporter ZnT8 may influence the development of ZnT8 immunogenicity and the phenotypic features of T1D. The solute carrier family 30 member 8 (SLC30A8) gene, located in chromosome 8q24.11, encodes for the ZnT8 auto-antigen and comprises 369 amino acids (12,13). Notably, aa268-369 of the cytoplasmic domain of ZnT8, especially ZnT8-325R and ZnT8-325W, is the dominant epitope in T1D. A common non-synonymous single-nucleotide polymorphism (SNP) of SLC30A8 rs13266634 (C/T polymorphism) encodes either arginine (R) by the C allele or tryptophan (W) by the T allele at aa325 of ZnT8 (14) suggesting that rs13266634 SNP might be critical for humoral autoimmunity in T1D (11,15). Thus, the present study is based on the evidence that SLC30A8 gene polymorphism is involved in T1D development. The objective of this study was to investigate the association between rs13266634 C/T gene polymorphism and T1D among the children of Tamil Nadu and to apply these results in a meta-analysis to reveal the association between the SLC30A8 risk allele and T1D for comparison in different ethnic groups.


Methods


Subjects

The study subjects comprised 121 T1D patients from the Department of Diabetology, Government Rajaji Hospital in Madurai, Tamil Nadu, India, along with 214 their first degree relatives (120 parents and 94 siblings) as controls. All patients were evaluated by clinical history and routine laboratory tests. The patients met the revised criteria of the American Diabetes Association (ADA) for the screening of T1D (16). Genomic DNA was extracted from 5 mL of peripheral blood sample by salting out method (17).

Ethic board consent for the study was approved by the Institutional Ethics Committees of Govt. Rajaji Hospital (Ref. No. 23339/E4/3/10) and Madurai Kamaraj University (MKU/IRB/11/11) and consented in writing by the participants.


Genotype Analysis

Subjects were genotyped for rs13266634 C/T polymorphism of SLC30A8 gene by polymerase chain reaction (PCR)-restriction fragment length polymorphism (18,19). The region surrounding the polymorphism was amplified with the following primers: Forward, 5’-GGACAGAAAGAGTTCCCATAGCG-3’; Reverse, 5’-ATAGCAGCATGTTTGAAGGTGGC-3’. PCR was performed at 95 °C for 5 minutes, followed by 40 cycles at 94 °C for 40 seconds and 69 °C for 45 seconds. A final extension step was carried out at 72 °C for 5 minutes. The PCR products were digested using enzyme MSp1 (Thermo Scientific, USA) incubated at 37 °C for 4 hours and visualized on 2% agarose gel. In the wild-type genotype (CC) the fragments obtained were of 234 and 195 bp. In the heterozygote genotype (CT), three fragments were detected of 429, 234 and 195 bp. Only one fragment of 429 bp was identified in the homozygote genotype (TT).


Meta-analysis

An extensive literature search was done to examine the association between T1D and SLC30A8 gene. The original data were collected from the following electronic databases: PubMed, Elsevier, Science Direct, Web of Science and Google Scholar with key words “Zinc transporter protein member 8, ZnT8, SLC30A8 gene polymorphism, SLC30A8 or SLC30A8 variant, combined with autoimmunity, autoimmune diabetes, T1D mellitus”. All searches were done independently by more than two research investigators. The following inclusion criteria were applied: 1) studies should be case-controlled; and 2) all patients should meet the diagnostic criteria for T1D according to the ADA. Studies were excluded if they did not report on genotype frequency or if they had insufficient data.


Statistical Analysis

The obtained clinical data were subjected to Student t-test and χ2 test after segregating the data based on age, number and sex of the subjects. Odds ratio (OR) and their p-values were calculated by logistic regression, which was performed using STATA 14v software (STATA Corporation, College road, TX, USA). In addition, the transmission/disequilibrium test (TDT) was employed to detect preferential transmission from heterozygous parents to affected offspring (20). The TDT analysis was done by Haploview 4.2v. software (Broad Institute, Cambridge, MA, USA). The level of significance was set at p<0.05. Heterogeneity evaluation was performed by the Cochran’s Q-test (21) and p<0.10 was considered statistically significant. If not significant, OR and 95% confident interval (CI) was estimated by fixed effect model (22), otherwise the random effect model was used (23). Heterogeneity of the data was quantified using the I2 test (24). I2 value of 25%, 50% and 75% were nominally considered low, moderate and high estimates, respectively. Funnel plot and Egger’s linear regression test was used for the analysis of publication bias (25). Meta-analysis was performed with Rev Man 5.0v. software (RevMan 5.0, The Cochrane Collaboration, Oxford, UK).


Results

The demographic details of the T1D subjects and controls are given in Table 1. There was no significant differences observed in allele (OR=0.92; CI=0.33-2.58; p=0.88) and genotype (CC: OR=0.92; CI=0.58-1.47; p=0.74; CT: OR=1.05; CI=0.64-1.71; p=0.82; TT: OR=1.13; CI=0.42-3.00; p=0.80) frequencies of rs13266634 C/T between T1D patients and controls, respectively (Table 2). Upon analysis of 30 parent-offspring trios (one affected child and both parents) of the study cohort, TDT analysis identified over-transmission of mutant T allele of rs13266634 C/T polymorphism from parents to affected children (T: U=9:7; MAF=0.194; χ2=0.25; p=0.61) without statistical significance.

Meta-analysis of the data via literature survey was able to retrieve 18 studies. Of these, nine were excluded after screening the abstracts, review and irrelevant subject matter. Three studies did not provide comprehensive information. Two studies were not considered as they provided insufficient genotype frequencies. The remaining four studies (14,26,27,28) associated with rs13266634 C/T polymorphism in the SLC30A8 gene of T1D, which met the required criteria, were included in the present meta-analysis. Along with the present study, a total of five eligible studies with a total of 10,376 T1D patients and 10,027 control subjects were included in the meta-analysis.

Characteristics of the said studies and the distribution of rs13266634 C/T genotypes and alleles in T1D patients and controls are given in Table 3. Overall effects of rs13266634 C allele frequency in T1D patients (OR=0.97; CI=0.92-1.01; p=0.10) based on pooled analysis were not different from the controls (Table 4). There was no evidence of virtual asymmetry (χ2= 0.29; I2=0%; Pheterogeneity= 0.99) which indicated that no publication bias crept in the meta-analysis (Figure 1).

In the Forest plot the area of squares, horizontal lines and diamond shows the weight of specific study, confidence intervals and the summary of fixed-effects OR, respectively (Table 4).

In the Funnel plot the open circle represents various studies considered for this plot correlation (Figure 1). No evidence of publication bias was found.


Discussion

ZnT8 is highly expressed in the pancreatic islet β-cells and recognized as one of the four major auto-antigens in T1D patients. It has been observed that autoantibodies are generated against ZnT8 prior to the onset of disease. It is known that rs13266634 C/T SNP is responsible for the autoimmune response to ZnT8 (12). The rs13266634 C/T plays a susceptibility role in the presence of impaired, autoimmunity-mediated β-cell dysfunction which leads to T1D development (13). Studies of the role of rs13266634 C/T polymorphism in T1D among a global population are scanty. This work appears to be the first family based TDT analysis on rs13266634 SNP with its allele transmission from parents to offspring. As for TDT results, the present study documents over-transmission of mutant T allele of rs13266634 in T1D. In a case control scenario, the present study indicates that there is a lack of association of rs13266634 C/T polymorphism to T1D. A few earlier studies also lent support to this contention in the Danish, Japanese and British populations (14,26,28). However, a German study indicates a higher occurrence of the C allele and CC genotype of rs13266634 C/T polymorphism in early onset of T1D patients compared to controls (27). A recent study revealed that an adjacent locus of rs2466293 in the SLC30A8 gene seems to predispose to the risk of T1D in individuals of non-European descent (29).

Until now, several publications have investigated the correlation of rs13266634 C/T polymorphisms with T1D (14,26,27,28). However, the results remain inconclusive. In order to reach a more concrete opinion of this contentious matter, a meta-analysis was performed with expanded sample size, aiming to explore the relationship of polymorphism at rs13266634 C/T of the SLC30A8 gene with susceptibility to T1D. However the result of the meta-analysis indicated that the C allele conferred no risk in the development of T1D. Nevertheless, we should point out that one of the previous meta-analyses on T2D revealed that the rs13266634 C/T polymorphism is significantly associated with impaired glucose tolerance (30).


Study Limitations

The study is limited by a relatively small number of subjects. Varied studies from different ethnicities with large sample size are required to conclusively confirm the role of rs13266634 C/T polymorphism in T1D.


Conclusion

This result demonstrates that the allele, genotype, genetic models and allele transmission of rs13266634 C/T polymorphism are not strongly associated with T1D in the children of a Tamil Nadu population. The meta-analysis also indicates that the rs13266634 C/T polymorphism was not associated with T1D.


Acknowledgements

This work was supported by the grants of the Tamil Nadu State Council for Science and Technology, Govt. of Tamil Nadu. India.

Ethics

Ethics Committee Approval: Ethic board consent for the study was approved by the Institutional Ethics Committees of Govt. Rajaji Hospital (Ref. No. 23339/E4/3/10) and Madurai Kamaraj University (MKU/IRB/11/11).

Informed Consent: All parents were informed about the purpose of the study, and a signed consent for study participation was obtained.

Peer-review: Internally peer-reviewed.

Authorship Contributions

Concept: Mariakuttikan Jayalakshmi, Design: Mariakuttikan Jayalakshmi, Ramasamy Thirunavukkarasu, Data Collection or Processing: Ramasamy Thirunavukkarasu, Arthur Joseph Asirvatham, Ayyappan Chitra, Analysis or Interpretation: Ramasamy Thirunavukkarasu, Mariakuttikan Jayalakshmi, Literature Research: Ramasamy Thirunavukkarasu, Writing: Ramasamy Thirunavukkarasu, Mariakuttikan Jayalakshmi.

Financial Disclosure: This study was supported by the grants of the Tamil Nadu State Council for Science and Technology (ref. No.TNSCST/S&T Proj/ MS/ VR/2011-2012 by Government of Tamil Nadu and UGC - NRCBS, CAS, DSTPURSE, and UGC - Meritorious fellowship by Government of India.

  1. Gillespie KM. Type 1 diabetes: pathogenesis and prevention. CMAJ 2006;2:165-170.
  2. Pihoker C, Gilliam LK, Hampe CS, Lernmark A. Autoantibodies in diabetes. Diabetes 2005;54(Suppl 2):52-61.
  3. Lieberman SM, DiLorenzo TP. A comprehensive guide to antibody and T-cell responses in type 1 diabetes. Tissue Antigens 2003;62:359-377.
  4. Arvan P, Pietropaolo M, Ostrov D, Rhodes CJ. Islet autoantigens: Structure, function, localization, and regulation. Cold Spring Harb Perspect Med 2012;2.
  5. Kambe T, Yamaguchi-Iwai Y, Sasaki R, Nagao M. Overview of mammalian zinc transporters. Cell Mol Life Sci 2004;61:49-68.
  6. Chimienti F, Devergnas S, Favier A, Seve M. Identification and cloning of a beta-cell-specific zinc transporter, ZnT-8, localized into insulin secretory granules. Diabetes 2004;53:2330-2337.
  7. Yi B, Huang G, Zhou Z. Different role of zinc transporter 8 between type 1 diabetes mellitus and type 2 diabetes mellitus. J Diabetes Investig 2016;7:459-465. Epub 2016 Jan 9.
  8. Vaziri-Sani F, Oak S, Radtke J, Lernmark K, Lynch K, Agardh CD, Cilio CM, Lethagen AL, Ortqvist E, Landin-Olsson M, Törn C, Hampe CS. ZnT8 autoantibody titers in type 1 diabetes patients decline rapidly after clinical onset. Autoimmunity 2010;43:598-606. Epub 2010 Mar 19
  9. Yang L, Luo S, Huang G, Peng J, Li X, Yan X, Lin J, Wenzlau JM, Davidson HW, Hutton JC, Zhou Z. The diagnostic value of zinc transporter 8 autoantibody (ZnT8A) for type 1 diabetes in Chinese. Diabetes Metab Res Rev 2010;26:579-584.
  10. Yi B, Huang G, Zhou ZG. Current and Future Clinical Applications of Zinc Transporter-8 in Type 1 Diabetes Mellitus. Chin Med J (Engl) 2015;128:2387-2394.
  11. Elmaoğulları S, Uçaktürk SA, Elbeg Ş, Döğer E, Tayfun M, Gürbüz F, Bideci A. Prevalence of ZnT8 antibody in Turkish children and adolescents with new onset type 1 diabetes. J Clin Res Pediatr Endocrinol 2018;10:108-112. Epub 2017 Sep 25
  12. Wenzlau JM, Juhl K, Yu L, Moua O, Sarkar SA, Gottlieb P, Rewers M, Eisenbarth GS, Jensen J, Davidson HW, Hutton JC. The cation efflux transporter ZnT8 (Slc30A8) is a major autoantigen in human type 1 diabetes. Proc Natl Acad Sci U S A 2007;104:17040-17045. Epub 2007 Oct 17
  13. Wenzlau JM, Liu Y, Yu L, Moua O, Fowler KT, Rangasamy S, Walters J, Eisenbarth GS, Davidson HW, Hutton JC. A common nonsynonymous single nucleotide polymorphism in the SLC30A8 gene determines ZnT8 autoantibody specificity in type 1 diabetes. Diabetes 2008;57:2693-2697. Epub 2008 Jun 30
  14. Kawasaki E, Uga M, Nakamura K, Kuriya G, Satoh T, Fujishima K, Ozaki M, Abiru N, Yamasaki H, Wenzlau JM, Davidson HW, Hutton JC, Eguchi K. Association between anti‐ZnT8 autoantibody specificities and SLC30A8 Arg325Trp variant in Japanese patients with type 1 diabetes. Diabetologia 2008;51:2299-2302. Epub 2008 Oct 11
  15. Wenzlau JM, Moua O, Liu Y, Eisenbarth GS, Hutton JC, Davidson HW. Identification of a major humoral epitope in Slc30A8 (ZnT8). Ann N Y Acad Sci 2008;1150:252-255.
  16. American Diabetes Association. Standards of Medical Care in Diabetes - 2010. Diabetes Care 2010;33(Suppl 1):11-61.
  17. Miller SA, Dykes DD, Polesky HF. A simple salting out procedure for extracting DNA from human nucleated cells. Nucleic Acids Res 1988;16:1215.
  18. Feghih H, Khatami SR, Azarpira N, Foroughmand AM. SLC30A8 gene polymorphism (rs13266634 C/T) and type 2 diabetes mellitus in south Iranian population. Mol Biol Rep 2014;41:2709-2715. Epub 2014 Jan 22
  19. Huang Q, Yin JY, Dai XP, Wu J, Chen X, Deng CS, Yu M, Gong ZC, Zhou HH, Liu ZQ. Association analysis of SLC30A8 rs13266634 and rs16889462 polymorphisms with type 2 diabetes mellitus and repaglinide response in Chinese patients. Eur J Clin Pharmacol 2010;66:1207-1215. Epub 2010 Aug 31
  20. Spielman RS, McGinnis RE, Ewens WJ. Transmission test for linkage disequilibrium: the insulin gene region and insulin-dependent diabetes mellitus (IDDM). Am J Hum Genet 1993;52:506-516.
  21. Cochran WG. The combination of estimates from different experiments. Biometrics 1954;10:101-129.
  22. Mantel N, Haenszel W. Statistical aspects of the analysis of data from retrospective studies of disease. J Natl Cancer Inst 1959;22:719-748.
  23. DerSimonian R, Laird N. Meta-analysis in clinical trials. Control Clin Trials 1986;7:177-188.
  24. Higgins JP, Thompson SG. Quantifying heterogeneity in a meta-analysis. Stat Med 2002;21:1539-1558.
  25. Egger M, Davey Smith G, Schneider M, Minder C. Bias in meta-analysis detected by a simple, graphical test. BMJ 1997;315:629-634.
  26. Brorsson C, Bergholdt R, Sjögren M, Eising S, Sorensen KM, Hougaard DM, Orho-Melander M, Groop L, Pociot F. A non-synonymous variant in SLC30A8 is not associated with type 1 diabetes in the Danish population. Mol Genet Metab 2008;94:386-388. Epub 2008 Apr 8
  27. Gohlke H, Ferrari U, Koczwara K, Bonifacio E, Illig T, Ziegler AG. SLC30A8 (ZnT8) polymorphism is associated with young age at type 1 diabetes onset. Rev Diabet Stud 2008;5:25-27. Epub 2008 May 10
  28. Raj SM, Howson JM, Walker NM, Cooper JD, Smyth DJ, Field SF, Stevens, Todd JA. No association of multiple type 2 diabetes loci with type 1 diabetes. Diabetologia 2009;52:2109-2116. Epub 2009 May 20
  29. Gomes KF, Semzezem C, Batista R, Fukui RT, Santos AS, Correia MR, Passos-Bueno MR, Silva ME. Importance of Zinc Transporter 8 Autoantibody in the Diagnosis of Type 1 Diabetes in Latin Americans. Sci Rep 2017;16:7:207.
  30. Xu K, Zha M, Wu X, Yu Z, Yu R, Xu X, Chen H, Yang T. Association between rs13266634 C/T polymorphisms of solute carrier family 30 member 8 (SLC30A8) and type 2 diabetes, impaired glucose tolerance, type 1 diabetes--a meta-analysis. Diabetes Res Clin Pract 2011;91:195-202. Epub 2010 Dec 4