- Research article
- Open Access
- Open Peer Review
Evidence for the association of the SLC22A4 and SLC22A5 genes with Type 1 Diabetes: a case control study
BMC Medical Geneticsvolume 7, Article number: 54 (2006)
Type 1 diabetes (T1D) is a chronic, autoimmune and multifactorial disease characterized by abnormal metabolism of carbohydrate and fat. Diminished carnitine plasma levels have been previously reported in T1D patients and carnitine increases the sensitivity of the cells to insulin. Polymorphisms in the carnitine transporters, encoded by the SLC22A4 and SLC22A5 genes, have been involved in susceptibility to two other autoimmune diseases, rheumatoid arthritis and Crohn's disease. For these reasons, we investigated for the first time the association with T1D of six single nucleotide polymorphisms (SNPs) mapping to these candidate genes: slc2F2, slc2F11, T306I, L503F, OCTN2-promoter and OCTN2-intron.
A case-control study was performed in the Spanish population with 295 T1D patients and 508 healthy control subjects. Maximum-likelihood haplotype frequencies were estimated by applying the Expectation-Maximization (EM) algorithm implemented by the Arlequin software.
When independently analyzed, one of the tested polymorphisms in the SLC22A4 gene at 1672 showed significant association with T1D in our Spanish cohort. The overall comparison of the inferred haplotypes was significantly different between patients and controls (χ2 = 10.43; p = 0.034) with one of the haplotypes showing a protective effect for T1D (rs3792876/rs1050152/rs2631367/rs274559, CCGA: OR = 0.62 (0.41–0.93); p = 0.02).
The haplotype distribution in the carnitine transporter locus seems to be significantly different between T1D patients and controls; however, additional studies in independent populations would allow to confirm the role of these genes in T1D risk.
Type 1 diabetes (T1D) is a multifactorial autoimmune T-cell-mediated disease resulting from selective destruction of the insulin producing β cells in the pancreatic islets, leading to an absolute insulin deficiency. The risk of developing T1D is determined by a complex interaction between multiple genetic and environmental factors. Although susceptibility to disease is strongly associated with alleles in the major histocompatibility complex (MHC) [1, 2], there are more than 20 putative T1D susceptibility regions identified by linkage and association studies [3, 4]. At present, several non-MHC susceptibility loci with modest genetic effects have been clearly defined. However, it is well known that many non-MHC loci predisposing to T1D remain as yet undefined .
Type 1 diabetes is a chronic degenerative disease, with altered metabolism characterized by hyperglycemia and ketoacidosis and T1D patients depend on exogenous insulin to sustain life. The role of the carnitine system in cell metabolism is mainly known in the mitochondria, where the interaction between fatty acid and glucose metabolism is fundamental for cell energy production [6, 7]. However, carnitine not only contributes to the transport of activated long-chain fatty acids into mitochondria for β-oxidation, but it also increases the sensitivity of cells to insulin . Decreased plasma carnitine levels have been reported in patients with type 2 diabetes [9–11] and some studies have investigated the carnitine status in T1D, finding similar results [12–15].
Adequate carnitine levels are required for normal lipid metabolism and are important for energy metabolism . One important component of the carnitine system is the plasma membrane carnitine transporters, named organic cation transporters (OCTN1 and OCTN2) encoded by the SLC22A4 and SLC22A5 genes, respectively. Both genes map to the cytokine gene cluster on chromosome 5q31 and show 88% homology and 77% identity in their sequences. Despite OCTN1 and OCTN2 are considered as carnitine transporters, only OCTN2 is a high-affinity human carnitine transporter, while the carnitine transport activity of OCTN1 is very low [17, 18]. In fact, a recent study has reported that the main substrate of this transporter is the ergothioneine, an intracellular antioxidant with metal ion affinity, which is transported one hundred times more efficiently than carnitine . OCTN2 is widely expressed in many adult tissues, among them in pancreas, and it participates, at least in part, in proton/organic cation antiport at the renal apical plasma membrane level .
Recent reports performed associations of some polymorphisms within the SLC22A4 and SLC22A5 genes with two other autoimmune complex diseases (rheumatoid arthritis and Crohn's disease) [20, 21]. The purpose of this study was to investigate the influence of the SLC22A4 and SLC22A5 genes in type 1 diabetes risk in the Spanish population. Six SNPs along these genes were considered good markers to map this region: slc2F2 (rs3792876) and slc2F11 (rs 2306772), which are SNPs in the SLC22A4 gene, were originally associated with RA susceptibility . Other polymorphisms within the same linkage disequilibrium (LD) block (slc2F1) and slc2F2 showed lack of association with T1D . We have studied two additional SNPs in this gene: T306I (rs272893) and L503F (rs1050152, SNP located in exon 9 of SLC22A4). The OCTN2-promoter (rs2631367) is a transversion (-207G>C) disrupting a heat shock element in the promoter region of the SLC22A5 gene and it has been described, together with L503F, as etiologic variant in Crohn disease . Finally, we analyzed an intronic SNP in the SLC22A5 gene: the OCTN2-intron (rs274559) in order to define haplotypes within these genes.
We studied 295 unrelated Spanish white T1D patients (149 men and 146 women) diagnosed according to the criteria of the American Diabetes Association (ADA) and 508 healthy controls recruited among blood donors. Both groups ethnically matched from the Madrid area. The age at onset for the T1D patients range from 1 to 55 years old (median age at onset 15 years) and all subjects were insulin-dependent at the time to study. The protocol followed the principles expressed in the Declaration of Helsinki and it was approved by the Hospital Ethics Committee.
SNPs slc2F2 (rs3792876), slc2F11 (rs2306772), T306I (rs272893), L503F (rs1050152) and OCTN2-intron (rs274559) were genotyped by TaqMan Assays on Demand under conditions recommended by manufacturer (Applied Biosystems), with identification numbers: C__3170428_10, C__3170458_1_, C__3170445_1_, C__3170459_10 and C__1173605_1, respectively. For the SLC22A5 promoter -207G>C (rs2631367), a TaqMan Assay by Design was performed.
Differences in allele or genotype frequencies for each marker were calculated by Chi-square, or Fisher's exact test when necessary. Associations were estimated by the odds ratio (OR) with 95% confidence interval (CI). Statistical analysis used Epi Info v. 6.02 (CDC Atlanta USA).
Maximum-likelihood haplotype frequencies were estimated by applying the Expectation Maximization (EM) algorithm implemented by the Arlequin software , with number of iterations set at 5000 and initial conditions at 50, with an epsilon value of 10-7. This software yields the estimated frequency of each haplotype, but it does not output the expected haplotypes for each individual. The frequency data were transformed into absolute numbers multiplying the frequencies presented in Table 3 by the total number of haplotypes in each group (patients and controls). Then, these values were introduced into contingency tables to calculate Chi-square and p-values.
Six SNPs were studied in order to check the role of SLC22A4 and SLC22A5 genes in T1D predisposition (see figure 1). Slc2F2, Slc2F11, T306I and L503F map in the SLC22A4 gene and two additional SNPs, OCTN2-promoter and OCTN2-intron, are located in the SLC22A5 gene. The analysis of the control cohort showed complete linkage disequilibrium (LD) between Slc2F11 and slc2F2, and T306I was also found in complete LD with OCTN2-intron. Both SNPs (slc2F11 and T306I) were not considered in the subsequent case-control study, because they do not supply additional information. These polymorphisms conformed to Hardy-Weinberg equilibrium.
The association of the intronic slc2F2 polymorphism in the SLC22A4 gene with rheumatoid arthritis (RA) was originally reported in a Japanese population . This SNP was even described as an etiological variant in Japan, being the mutant homozygous genotype (slc2F2*TT) associated with RA. However, this association could not be replicated either in British  or in Spanish populations . Table 1 shows the phenotype and genotype distribution of this marker in the Spanish T1D patient and control cohorts, and lack of significant association with the disease was observed. The allelic distribution was similar to that previously observed in Spanish RA patients, (minor allele frequency slc2F2*T was 7.5% in RA and 9.5% in T1D patients). We tested the carrier rate of this allele (CT+TT vs. CC) to find a possible association with T1D, but a negative result was again obtained: OR = 1.30 (0.87–1.94); p = 0.17.
Regarding the L503F (1672G>C) polymorphism, it was described by Peltekova et al.  as a functional variant associated with Crohn's disease, together with the OCTN2-promoter (-207G>C) variant. In fact, the authors found that the odds ratios conferred by the allele 1672*T, by the allele -207*C or by the TC minihaplotype were all similar. Both polymorphisms were also in strong linkage disequilibrium in our population (D' = 0.86) . However, their etiological role could not be verified in Spanish Crohn's patients and they were not described as causative polymorphisms, but as genetic markers of risk or protection haplotypes. In our diabetic population, the allele 1672*T increased disease predisposition [OR = 1.25 (1.01–1.54); p = 0.034] and the 1672*CC genotype showed a protective effect (Table 1), but no significant differences between patients and controls were found for alleles and genotypes of the OCTN2-promoter (-207G>C) and of the other SLC22A5 intronic variant (Table 2).
We continued our study by analyzing the inferred OCTN1-OCTN2 haplotypes, since haplotypes define better a DNA fragment where a susceptibility or protective gene could be located. Table 3 shows the haplotypes estimated by the Expectation-Maximization (EM) algorithm implemented by the Arlequin software with a frequency over 1%. The overall comparison of haplotypes in a 5 × 2 contingency table rendered a significant difference between patients and healthy controls (χ2 = 10.43; p = 0.034), and being a unique omnibus comparison, this result does not need ulterior correction. These results evidenced one protection CCGA [OR = 0.62 (0.41–0.93); p = 0.02] and another risk CTCA [OR = 1.23 (1.00–1.52); p = 0.05] haplotypes, although they did not withstand correction for multiple testing (Table 3). The effect of the latter is secondary to the protective haplotype, as evidenced when this one is eliminated from the comparison and no significant predisposition effect was observed [OR = 1.16 (0.93–1.44); p = 0.18]; however, when the eliminated haplotipe is CTCA, a significant protection effect is still found for the CCGA haplotype [OR = 0.67 (0.44–1.02); p = 0.047). Moreover, as other haplotypes carrying the 1672*C allele do not show any influence in T1D predisposition, the protective effect of the 1672*C allele (p = 0.034) is due to this CCGA haplotype.
To complete the study in our population, we analyzed children and adolescents with type 1 diabetes fixing the cut-off age in 15 years old (median age at onset in our cohort). No differences between young patients and either adult patients or controls were found for any isolated variant or for the inferred haplotypes (data not shown). Finally, the frequencies of both, polymorphisms and inferred haplotypes, between male and female T1D patients were similar (data not shown).
The chromosomal region 5q31 contains several genes involved in immune and inflammatory responses and the SLC22A4 and SLC22A5 genes were associated with two autoimmune diseases (rheumatoid arthritis and Crohn's disease) [20, 21]. Additionally, several genomewide scans for type 1 diabetes have identified susceptibility loci on different chromosomes, including the region 5q [27, 28]. Moreover, carnitine has been described to increase sensitivity of cells to insulin , which could hypothetically lend cells more prone to an autoimmune attack. For all these reasons we considered that these genes could be candidates to modify the susceptibility to another autoimmune disease, as type 1 diabetes.
We tested six polymorphisms in the SLC22A4 and SLC22A5 genes and in five of them no significant independent association with type 1 diabetes could be found (Tables 1 and 2). The power of this study considering a relative risk of 1.5 was 75% for Slc2F2, 86% for L503F and 99% for OCTN2-promoter and OCTN2-intron, as calculated by an UCLA Department of Statistics sofware . Therefore, we can exclude the studied SLC22A5 markers as causative candidates for T1D in Spanish patients, but in the case of Slc2F2 our results are not conclusive, albeit a recently published study does not find association of Slc2F1 and Slc2F2 SLC22A4 polymorphisms within a well-powered T1D cohort , in agreement with our data.
The overall distribution of the estimated haplotypes was significantly different between patients and healthy controls. Being ours the first study carried out in T1D, replication studies in different T1D populations would be necessary to firmly establish the role of these genes in autoimmune diabetes. We proved by using a stepwise procedure that the association with T1D of the protection haplotype CCGA is primary. The etiological polymorphisms previously found associated with increased susceptibility to RA (slc2F2*T) and to Crohn's disease (L503F*T) do not display a causative role by themselves in our population.
Low carnitine plasma levels have been found in children and adolescents with type 1 diabetes [12, 13, 15]. For this reason, we decided to conclude the study with an age-stratified analysis. No significant differences in the association of each polymorphism were observed when grouping by age at onset or gender. In the study of the inferred haplotypes again no differences were found between any group of patients and when we compared each group with controls. Therefore, it seems that the protective effect of the inferred haplotype is independent of both sex and age at disease onset.
A previous report about disequilibrium blocks at 5q31 in European-derived population  demonstrated limited haplotype diversity, concordantly with our data (5 haplotypes found out of the 24 = 16 theoretically possible). One of the defined blocks (92 kb long) corresponds to the carnitine transporter genes locus. These authors indicated that the chromosomal region 5q31 is divided into discrete blocks displaying complete linkage disequilibrium (LD). However, some degree of LD extending beyond the SLC22A4 and SLC22A5 genes could imply that most probably these genes, or other in LD with them, are responsible for the reported effect on T1D risk.
In conclusion, the overall comparison of haplotypes within the chromosomal region where the carnitine transporter genes map seems to be different between type 1 diabetes patients and healthy controls in our Spanish population. The reduced carnitine plasma levels found in both type 1 and type 2 diabetes patients [12, 13] could be explained by a mere increase in the activity of the carnitine transporters to supply the higher energetic cellular demand mainly provided by lipid metabolism in young diabetes patients. However, our results with the 1672*C allele and with the inferred haplotypes, the effect of carnitine sensitizing cells to insulin and potentially rendering them more amenable to an immune attack, and also the LD-block data described by Daly , support the role of the OCTN genes in T1D risk. Additional studies in independent populations and in both type 1 and type 2 diabetes patients will be needed to confirm the putative influence of the SLC22A4 and SLC22A5 genes in these diseases.
Major Histocompatibility Complex
Organic Cation Transporter
Single Nucleotide Polymorphism
Solute Carrier Family 22
Type 1 Diabetes
Nerup J, Platz P, Andersen OO, Christy M, Lyngsoe J, Poulsen JE, Ryder LP, Nielsen LS, Thomsen M, Svejgaard A: HL-A antigens and diabetes mellitus. Lancet. 1974, 2: 864-866. 10.1016/S0140-6736(74)91201-X.
Pociot F, McDermott MF: Genetics of type 1 diabetes mellitus. Genes Immun. 2002, 3: 235-249. 10.1038/sj.gene.6363875.
Field LL: Genetic linkage and association studies of Type I diabetes: challenges and rewards. Diabetologia. 2002, 45: 21-35. 10.1007/s125-002-8241-7.
Davies JL, Kawaguchi Y, Bennett ST, Copeman JB, Cordell HJ, Pritchard LE, Reed PW, Gough SC, Jenkins SC, Palmer SM, et al: A genome-wide search for human type 1 diabetes susceptibility genes. Nature. 1994, 371: 130-136. 10.1038/371130a0.
Anjos S, Polychronakos C: Mechanisms of genetic susceptibility to type I diabetes: beyond HLA. Mol Genet Metab. 2004, 81: 187-195. 10.1016/j.ymgme.2003.11.010.
Ramsay RR: The carnitine acyltransferases: modulators of acyl-CoA-dependent reactions. Biochem Soc Trans. 2000, 28: 182-186.
Foster DW: The role of the carnitine system in human metabolism. Ann N Y Acad Sci. 2004, 1033: 1-16. 10.1196/annals.1320.001.
Proulx F, Lacroix J, Qureshi IA, Nadeau D, Gauthier M, Lambert M: Acquired carnitine abnormalities in critically ill children. Eur J Pediatr. 1997, 156: 864-869. 10.1007/s004310050732.
Okuda Y, Kawai K, Murayama Y, Yamashita K: Postprandial changes in plasma ketone body and carnitine levels in normal and non-insulin-dependent diabetic subjects. Endocrinol Jpn. 1987, 34: 415-422.
Pregant P, Kaiser E, Schernthaner G: No effect of insulin treatment or glycemic improvement on plasma carnitine levels in type 2 diabetic patients. Clin Investig. 1993, 71: 610-612. 10.1007/BF00184484.
Tamamogullari N, Silig Y, Icagasioglu S, Atalay A: Carnitine deficiency in diabetes mellitus complications. J Diabetes Complications. 1999, 13: 251-253. 10.1016/S1056-8727(99)00052-5.
Soltesz G, Melegh B, Sandor A: The relationship between carnitine and ketone body levels in diabetic children. Acta Paediatr Scand. 1983, 72: 511-515.
Winter SC, Simon M, Zorn EM, Szabo-Aczel S, Vance WH, O'Hara T, Higashi L: Relative carnitine insufficiency in children with type I diabetes mellitus. Am J Dis Child. 1989, 143: 1337-1339.
Pregant P, Schernthaner G, Legenstein E, Lienhart L, Bruck S, Schnack C, Kaiser E: [Decreased plasma carnitine in Type I diabetes mellitus]. Klin Wochenschr. 1991, 69: 511-516. 10.1007/BF01649287.
Mamoulakis D, Galanakis E, Dionyssopoulou E, Evangeliou A, Sbyrakis S: Carnitine deficiency in children and adolescents with type 1 diabetes. J Diabetes Complications. 2004, 18: 271-274. 10.1016/S1056-8727(03)00091-6.
Akisu M, Kultursay N, Coker I, Huseyinov A: Myocardial and hepatic free carnitine concentrations in pups of diabetic female rats. Ann Nutr Metab. 2002, 46: 45-48. 10.1159/000046752.
Tamai I, Yabuuchi H, Nezu J, Sai Y, Oku A, Shimane M, Tsuji A: Cloning and characterization of a novel human pH-dependent organic cation transporter, OCTN1. FEBS Lett. 1997, 419: 107-111. 10.1016/S0014-5793(97)01441-5.
Grundemann D, Harlfinger S, Golz S, Geerts A, Lazar A, Berkels R, Jung N, Rubbert A, Schomig E: Discovery of the ergothioneine transporter. Proc Natl Acad Sci U S A. 2005, 102: 5256-5261. 10.1073/pnas.0408624102.
Tamai I, Ohashi R, Nezu J, Yabuuchi H, Oku A, Shimane M, Sai Y, Tsuji A: Molecular and functional identification of sodium ion-dependent, high affinity human carnitine transporter OCTN2. J Biol Chem. 1998, 273: 20378-20382. 10.1074/jbc.273.32.20378.
Tokuhiro S, Yamada R, Chang X, Suzuki A, Kochi Y, Sawada T, Suzuki M, Nagasaki M, Ohtsuki M, Ono M, Furukawa H, Nagashima M, Yoshino S, Mabuchi A, Sekine A, Saito S, Takahashi A, Tsunoda T, Nakamura Y, Yamamoto K: An intronic SNP in a RUNX1 binding site of SLC22A4, encoding an organic cation transporter, is associated with rheumatoid arthritis. Nat Genet. 2003, 35: 341-348. 10.1038/ng1267.
Peltekova VD, Wintle RF, Rubin LA, Amos CI, Huang Q, Gu X, Newman B, Van Oene M, Cescon D, Greenberg G, Griffiths AM, St George-Hyslop PH, Siminovitch KA: Functional variants of OCTN cation transporter genes are associated with Crohn disease. Nat Genet. 2004, 36: 471-475. 10.1038/ng1339.
Smyth DJ, Howson JM, Payne F, Maier LM, Bailey R, Holland K, Lowe CE, Cooper JD, Hulme JS, Vella A, Dahlman I, Lam AC, Nutland S, Walker NM, Twells RC, Todd JA: Analysis of polymorphisms in 16 genes in type 1 diabetes that have been associated with other immune-mediated diseases. BMC Med Genet. 2006, 7: 20-10.1186/1471-2350-7-20.
Stefan Schneider DRLE: Arlequin: A Software for Population Genetics Data Analysis. 2000, Geneva, University of Geneva, [http://anthro.unige.ch/arlequin]2.000
Barton A, Eyre S, Bowes J, Ho P, John S, Worthington J: Investigation of the SLC22A4 gene (associated with rheumatoid arthritis in a Japanese population) in a United Kingdom population of rheumatoid arthritis patients. Arthritis Rheum. 2005, 52: 752-758. 10.1002/art.20877.
Martinez A, Valdivia A, Pascual-Salcedo D, Balsa A, Fernandez-Gutierrez B, De la Concha E, Urcelay E: Role of SLC22A4, SLC22A5, and RUNX1 genes in rheumatoid arthritis. J Rheumatol. 2006, 33: 842-846.
Martinez A, Del Carmen Martin M, Mendoza JL, Taxonera C, Diaz-Rubio M, de la Concha EG, Urcelay E: Association of the organic cation transporter OCTN genes with Crohn's disease in the Spanish population. Eur J Hum Genet. 2006, 14: 222-226. 10.1038/sj.ejhg.5201529.
Nerup J, Pociot F: A genomewide scan for type 1-diabetes susceptibility in Scandinavian families: identification of new loci with evidence of interactions. Am J Hum Genet. 2001, 69: 1301-13. Epub 2001 Oct 11.. 10.1086/324341.
Lindgren CM, Widen E, Tuomi T, Li H, Almgren P, Kanninen T, Melander O, Weng J, Lehto M, Groop LC: Contribution of known and unknown susceptibility genes to early-onset diabetes in scandinavia: evidence for heterogeneity. Diabetes. 2002, 51: 1609-1617.
UCLA Department of Statistics. [http://calculatorsstatuclaedu/powercalc].
Daly MJ, Rioux JD, Schaffner SF, Hudson TJ, Lander ES: High-resolution haplotype structure in the human genome. Nat Genet. 2001, 29: 229-232. 10.1038/ng1001-229.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2350/7/54/prepub
We thank Carmen Martínez for expert technical assistance. Alfonso Martínez and Jose Luis Santiago are recipients of FIS contracts (CP04/00175 and CM05/00216). Elena Urcelay holds a "Ramón y Cajal" contract of the Spanish Science and Technology Ministry. This work was supported by grant FIS PI05/1221.
The author(s) declare that they have no competing interests.
JLS carried out the genotyping of the patients and a great part of the controls, participated in the statistical analysis and drafted the manuscript.
AM carried out a part of the genotyping of control samples and participated in the statistical analysis.
HdlC made the diagnosis and collaborated in collection of samples.
MAF participated in the genotyping and collection of samples.
MFA participated in the coordination of the study and helped to collect the DNA samples.
EgdlC coordinated the study and critically revised the manuscript.
EU conceived of the study, participated in the statistical analysis and wrote the major part of the manuscript.
All authors read and approve the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.