Printer Friendly

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


Type 1 diabetes (T1D) is a complex, multifactorial disease caused by the selective destruction of insulin-producing pancreatic [beta]-cells (1,2). The autoimmune destruction of pancreatic [beta]-cells by pathogenic T cells predominately targets a number of well-known [beta]-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 [beta]-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.



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 [degrees]C for 5 minutes, followed by 40 cycles at 94 [degrees]C for 40 seconds and 69 [degrees]C for 45 seconds. A final extension step was carried out at 72 [degrees]C for 5 minutes. The PCR products were digested using enzyme MSp1 (Thermo Scientific, USA) incubated at 37 [degrees]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).


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 [chi square] 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).


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; [chi square]=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 ([chi square]= 0.29; [I.sup.2]=0%; [P.sub.heterogeneity]= 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.


ZnT8 is highly expressed in the pancreatic islet [beta]-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 [beta]-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.


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.


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


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, DST-PURSE, 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, Torn 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.) Elmaogullari S, Ucakturk SA, Elbeg S, Doger E, Tayfun M, Gurbuz 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 H F. 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 J P, 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, Sjogren 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 K F, 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

Ramasamy Thirunavukkarasu (1), Arthur Joseph Asirvatham (2), Ayyappan Chitra (3), Mariakuttikan Jayalakshmi (1)

(1) Madurai Kamaraj University, School of Biological Sciences, Department of Immunology, Madurai, India

(2) Government Rajaji Hospital, Clinic of Diabetology, Madurai, India

(3) Government Rajaji Hospital, Institute of Child Health and Research Centre, Madurai, India

Address for Correspondence: Mariakuttikan Jayalakshmi PhD, Madurai Kamaraj University, School of Biological Sciences, Department of Immunology, Madurai, India

E-mail: - ORCID ID:

Confict of interest: None declared

Received: 28.07.2018

Accepted: 07.09.2018

DOI: 10.4274/jcrpe.galenos.2018.2018.0195

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.
Table 1. Demographic details of the type 1 diabetes patients and

                         Type 1 diabetes patients
Details                  Male              Female

No. of Subjects (n)      70                51
Age (year)               15.8[+ or -]10.2  22.8[+ or -]10.0
Age at diagnosis (year)  15.5[+ or -]8.4   14.8[+ or -]7.9
TDDM (year)               9.5[+ or -]5.8    8.0[+ or -]5.6


No. of Subjects (n)      214
Age (year)                32.2[+ or -]15.1
Age at diagnosis (year)
TDDM (year)

TDDM: time duration of diabetes mellitus

Table 2. SLC30A8 rs13266634C/T genotypes and allele frequencies in type
1 diabetes patients and healthy controls

rs13266634C/T      T1D patients (n=121)  Controls (n=214)

CC (RR)             77 (63.6%)           140 (65.4%)
CT (WW)             37 (30.6%)            63 (29.4%)
TT (WW)              7 (5.8%)             11 (5.2%)
CC+CT (RR+RW) (1)   114 (72.2%)          203 (73.3%)
CT+TT (RW+WW) (2)   44 (27.8%)            74 (26.7%)
C (R)              191 (78.9%)           343 (80.1%)
T (W)               51 (21.1%)            85 (19.9%)

rs13266634C/T      OR (95% CI)        p-value

CC (RR)            0.92 (0.58-1.47)   0.74
CT (WW)            1.05 (0.64-1.71)   0.82
TT (WW)            1.13 (0.42-3.00)   0.8
CC+CT (RR+RW) (1)  0.88 (0.33-2.33)   0.8
CT+TT (RW+WW) (2)  1.08 (0.67-1.72)   0.74
C (R)              0.92 (0.33- 2.58)  0.88
T (W)                                 -

OR: odds ratio, CI: confdence interval, n: number in sample, T1D: type
1 diabetes
(1) Dominant model (CC+CT vs TT)
(2) Recessive model (CT+TT vs CC)

Table 3. Distribution of SLC30A8 genotype and allele among type 1
diabetes patients and controls included in the meta-analysis

Study                Arms               C      T     CC   CT   TT

                     Cases (n=1530)      2100   960  736  628  166
Brorsson et al (26)
                     Controls (n=1478)   2045   911  725  595  158
                     Cases (n=171)        198   144   63   72   36
Kawasaki et al (14)
                     Controls (n=114)     130    98
                     Cases (n=874)       1193   555  400  393   81
Gohlke et al (27)
                     Controls (n=1021)   1416   626  493  430   98
                     Cases (n=7680)     10481  4879
Raj et al (28)
                     Controls (n=7200)   9937  4463
                     Cases (n=121)        191    51   77   37    7
Present study
                     Controls (n=214)     343    85  140   63   11

Study                Population    Method

Brorsson et al (26)  Danish        Sequencing & Taqman
Kawasaki et al (14)  Japanese      PCR-RFLP
Gohlke et al (27)    German        MOLDI-TOF-MS
Raj et al (28)       British       Taq-man PCR
Present study        South Indian  PCR-RFLP

PCR-RFLP: polymerase chain reaction-restriction fragment length
COPYRIGHT 2019 Galenos Yayinevi Tic. Ltd.
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2019 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Author:Thirunavukkarasu, Ramasamy; Asirvatham, Arthur Joseph; Chitra, Ayyappan; Jayalakshmi, Mariakuttikan
Publication:Journal of Clinical Research in Pediatric Endocrinology
Article Type:Report
Geographic Code:9INDI
Date:Mar 1, 2019
Previous Article:Evaluation of Renal Function in Obese Children and Adolescents Using Serum Cystatin C Levels, Estimated Glomerular Filtration Rate Formulae and...
Next Article:Epicardial Fat Thickness in Children with Classic Congenital Adrenal Hyperplasia.

Terms of use | Privacy policy | Copyright © 2022 Farlex, Inc. | Feedback | For webmasters |