The Variable Number of Tandem Repeats Upstream of the Insulin Gene Is a Susceptibility Locus for Latent Autoimmune Diabetes in Adults
http://www.100md.com
糖尿病学杂志 2006年第6期
1 Diabetes Research Laboratories, Oxford Centre for Diabetes, Endocrinology and Metabolism, University of Oxford, Oxford, U.K.
2 Wellcome Trust Centre for Human Genetics, University of Oxford, Oxford, U.K.
3 Institute of Biomedical and Clinical Sciences, Peninsula Medical School, Exeter, U.K.
4 Centre for Diabetes and Metabolic Medicine, Bart’s and The London Queen Mary’s School of Medicine and Dentistry, London, U.K.
5 School of Clinical Medical Sciences, University of Newcastle, Newcastle, U.K.
6 Diabetes Trials Unit, Oxford Centre for Diabetes, Endocrinology and Metabolism, University of Oxford, Oxford, U.K.
HWE, Hardy-Weinburg equilibrium; LADA, latent autoimmune diabetes in adults; LD, linkage disequilibrium; PH, protective haplotype; SNP, single nucleotide polymorphism; UKPDS, U.K. Prospective Diabetes Study; VNTR, variable number of tandem repeats; VPH, very protective haplotype; W2 Repository, Warren 2 Repository; Exeter YT2D, Exeter Young-Onset type 2 Diabetes study
ABSTRACT
The etiopathological relationship between latent autoimmune diabetes in adults (LADA) and classical type 1 (insulin dependent) diabetes remains unclear. Variation at the insulin gene variable number of tandem repeats (VNTR) minisatellite influences susceptibility to type 1 diabetes, but studies in LADA have been small and inconsistent. We examined the role of insulin gene variation (using flanking variants as surrogates for VNTR subtypes) in the largest case-control study of LADA to date (400 case and 332 control subjects). Highly significant associations were identified with disease, with dominant protective effects of the T allele at eC23HphI (odds ratio [OR] 0.42 [95% CI 0.31eC0.58], P = 2.4 x 10eC8), A allele at +1,404Fnu4HI (0.50 [0.36eC0.70], P = 3.2 x 10eC5), and C allele at +3,580MspI (0.55 [0.35eC0.85], P = 0.0046). As with type 1 diabetes, the eC23HphI variant (a surrogate for the subdivision of VNTR into class I and III alleles) most clearly defined susceptibility in LADA. However, there was no association with age at diagnosis or requirement for insulin therapy 6 years postdiagnosis. This study establishes that variation within the insulin gene region does influence susceptibility to LADA, with the direction and magnitude of effect indistinguishable from that previously reported for type 1 diabetes. In conclusion, differences in VNTR-encoded susceptibility do not explain the differences in clinical presentation that distinguish classical type 1 diabetes and LADA.
Latent autoimmune diabetes in adults (LADA) and type 1 (insulin dependent) diabetes are both characterized by islet autoimmunity, indicated by the presence of circulating islet autoantibodies. However, the later age of onset (after the peak of type 1 diabetes incidence) and less aggressive clinical course (no immediate requirement for insulin therapy) of LADA result in substantial clinical overlap with type 2 diabetes (noneCinsulin dependent) (1eC3). Qualitative or quantitative differences in genetic susceptibility may underlie the marked variation in disease course and clinical presentation between LADA and type 1 diabetes.
The insulin gene region (IDDM2) accounts for 10% of familial risk for type 1 diabetes (4). Most evidence indicates that the variable number of tandem repeats (VNTR) minisatellite (which maps upstream of the translation initiation site) is the probable etiological variant, though certain flanking single nucleotide polymorphisms (SNPs) that are highly correlated with VNTR variation may also play a role (5,6). In non-African populations, much of the variation in VNTR length (repeat unit number) and composition (repeat element sequence) can be captured by dichotomizing VNTR alleles into two classes: the short class I alleles (26eC63 repeats) and longer class III alleles (141eC209 repeats) (7). In type 1 diabetes, class I alleles confer susceptibility to disease and class III, protection (8). Because of strong linkage disequilibrium (LD) between VNTR allele structure and haplotypes spanning the insulin gene, this major division in VNTR structure can be assayed using the eC23HphI SNP (9).
Further refinement of type 1 diabetes susceptibility is possible by examining the finer structure within major VNTR classes. A discrete subgroup of class I alleles, termed "IDeC," has been reported to confer protection rather than susceptibility (7). Additionally, class IIIA and IIIB alleles may differ in the degree of protection they offer, an observation reflected in the designation of the extended haplotypes on which these alleles lie as the protective haplotypes (PHs) and very protective haplotypes (VPHs), respectively (9). However, recent reexamination of these subclass associations in additional samples has cast significant doubt on these latter observations (5).
It is unclear whether these susceptibility effects of the insulin gene are evident in LADA and/or whether differences in VNTR-encoded susceptibility contribute to the clinical differences between LADA and type 1 diabetes. Published studies have yielded conflicting results (10,11). One small study from South America (n = 44 LADA and 138 control subjects) reported a significant class I association with LADA (10). However, another larger study, performed in 95 LADA case subjects and 172 control subjects from Finland, found no association of the eC23HphI site with disease (11). The present study uses the largest collection of patients yet studied to determine whether the IDDM2 locus is associated with LADA.
The LADA group (n = 400) comprised autoantibody-positive patients from three sources, the U.K. Prospective Diabetes Study (UKPDS), the Warren 2 Repository (W2 Repository), and the Exeter Young-Onset type 2 Diabetes study (Exeter YT2D); 332 normoglycemic autoantibody-negative subjects from the Diabetes in Families study served as control subjects (Table 1). All were genotyped for three SNPs known to be surrogate markers for minisatellite lineages: eC23HphI (A T), +1,404Fnu4HI (C A), and +3,580MspI (C T) (7,9).
Comparisons of genotype frequency distributions revealed borderline evidence for heterogeneity between the different case groups (eC23HphI, P = 0.045; +1,404Fnu4HI, P = 0.033; +3,580MspI, P = 0.062) (Table 2). Accordingly, results are presented for the case groups combined as well as for individual case groups with the control group. Genotype frequencies at all sites were in Hardy-Weinburg equilibrium (HWE) in the control subjects, but there were modest departures in some of the case groups (Exeter YT2D, eC23HphI, P = 0.008; W2 Repository, +1,404Fnu4HI, P = 0.01; UKPDS, +3,580MspI, P = 0.024). Only the last of these departures was in a direction consistent with the observed association effect. We attribute the other deviations to stochastic variation, since extensive regenotyping and resequencing (and failure to detect unexpected haplotypic combinations) effectively excluded genotyping error.
All three SNPs were associated with LADA in the combined group, with dominant protective effects of the T allele at eC23HphI (odds ratio [OR] 0.42 [95% CI 0.31eC0.58], P = 2.4 x 10eC8), A allele at +1,404Fnu4HI (0.50 [0.36eC0.70], P = 3.2 x 10eC5), and C allele at +3,580MspI (0.55 [0.35eC0.85], P = 0.0046). Patterns of associations were similar under the additive model (data not shown). Comparisons of the individual case groups with the control subjects (Table 3) confirmed that these associations were not an artifact of case-group heterogeneity: all case-control comparisons showed similar ORs, although sample size differences meant that not all such comparisons were significant. Step-wise logistic regression analysis in the combined LADA group versus control subjects identified the eC23HphI site (indicating VNTR class) as the primary disease-associated variant (Wald 2, P = 2.8 x 10eC8); the other sites had no additional independent effect.
Given recent discussion over the appropriate definition of LADA (12,13), we repeated these analyses, excluding those patients in the UKPDS (n = 45), W2 Repository (n = 4), and Exeter YT2D (n = 6) groups diagnosed under the age of 30 years and in whom the interval between diagnosis and insulin therapy was between 3 and 6 months. The strength of the association was unchanged in this analysis (OReC23HphI = 0.42 [95% CI 0.30eC0.58], P = 5.6 x 10eC8).
All three SNPs were, as expected, in strong LD: pairwise D' values exceeded 0.94, although r2 values ranged from 0.06 to 0.89. Given the modest departures from HWE in the case groups, the expectation-maximization algorithm was reinitialized three times to ensure convergence. This, together with presence of high LD between the SNPs, enabled accurate haplotype estimation. Four haplotypes with estimated frequencies exceeding 1% were observed, two corresponding to the subdivision of class I alleles (into IC+/ID+ and IDeC), the others to the major class III haplotypes PH and VPH (Table 4). Haplotype trend regression revealed highly significant differences in haplotype frequency distributions between the combined case subjects compared with the control subjects (P = 7.9 x 10eC6). Class I haplotypes were more frequent in the LADA than in the control subjects, in keeping with classical type 1 diabetes.
Diplotype frequencies differed significantly between case and control subjects (P < 10eC4 overall). The OR (95% CI) for VNTR diplotypes estimated under the recessive model for class I haplotypes and the dominant model for class III haplotypes were as follows: IDeC homozygotes versus all others 1.82 (1.18eC2.84), P = 0.0045; IC+/ID+ homozygotes versus all others 1.67 (1.10eC2.58), P = 0.016; PH-containing diplotypes versus all others 0.50 (0.36eC0.69), P = 2.0 x 10eC5; and VPH-containing diplotypes versus all others 0.49 (0.29eC0.80), P = 0.0028. Similar ORs were seen for the comparisons of the separate case groups against the control subjects (data not shown).
There was no association between age at diagnosis and insulin gene variation (class III positive mean ± SD 45.8 ± 10.1 years vs. class I homozygote 46.5 ± 10.5, P = 0.54). Within the UKPDS group (the only group for which such data were available—insulin requirement was determined by a protocol-driven assessment [14]), there was no association with insulin therapy requirement at 6 years postdiagnosis (class III positive 35% vs. class I homozygote 34%, P = 0.88).
This study has established that the INS gene is a susceptibility locus for autoimmune diabetes developing in adulthood as well as in children. The eC23HphI site (or, equivalently, VNTR class) appears to be the primary disease-associated variant, as in classical type 1 diabetes. We conclude that the different clinical features of LADA and type 1 diabetes cannot be explained by an absence of IDDM2-encoded genetic susceptibility in adult-onset patients (11). From the data available, there is no suggestion that the magnitude of the INS VNTR effect size in LADA differs appreciably from that seen in type 1 diabetes. In the combined LADA group, the OR (95% CI) associated with dominantly protective class III alleles was 0.42 (0.31eC0.58), a value comparable with that found in studies of classical type 1 diabetes (for example, OR 0.36 [95% CI 0.23eC0.55] and 0.22 [0.12eC0.39], as calculated from previous reports [9,15]).
Failure to detect a VNTR association with LADA in a previous study (11) may reflect several factors, including reduced power (a lower frequency of the class III allele in Finns (11,16) combined with a smaller sample size) and/or genuine differences arising through variation in genetic background or ascertainment.
Fine-structure analysis of the region demonstrated similarities in the association patterns of LADA compared with those reported for type 1 diabetes. Diplotype analyses showed that the IDeC haplotype confers similar susceptibility to LADA to other class I haplotypes and that the class III PH and VPH had similar protective effects. Both are consistent with findings from a recent reevaluation of the insulin gene in type 1 diabetes (5), although the wide CIs indicate the need for cautious interpretation. The present study had >90% power to detect the dominantly protective association of class III alleles with OR 0.5 (for = 0.05). However, considerably larger numbers of subjects would be required to detect subtle differences in predisposition associated with specific class I or class III subgroups; the largest comparison of the protective effects associated with PH and VPH in type 1 diabetes featured meta-analysis of data from 3,000 pedigrees (5).
We have made no corrections for multiple testing; application of a conservative Bonferroni correction to the lowest P value obtained (P = 2.4 x 10eC8) would still retain significance (P < 10eC4), indicating that the observed associations are likely to be real.
Class III alleles have been associated with a decreased level of transcription of INS in the thymus during development (17). It has been postulated that this difference may affect central tolerance mechanisms, promoting negative selection of autoreactive T-cells and decreasing the propensity to develop autoimmunity against insulin (17). Similar protective mechanisms are likely to operate in LADA. However, as yet unresolved differences in class I subgroup susceptibility (or, alternatively, within HLA or other type 1 diabeteseCsusceptibility loci) could, through their consequences on the development of tolerance, contribute to the differences in clinical phenotype between classical type 1 diabetes and LADA.
RESEARCH DESIGN AND METHODS
All subjects were unrelated and of exclusively British/Irish Europid origin. LADA subjects (n = 400) came from three sources of type 2 diabeteseCdiagnosed patients: the UKPDS (n = 232 from the total cohort of 5,102) (14), the W2 Repository (n = 131 from the total sample of 4,000) (18,19,20), and the Exeter YT2D study (n = 37 from the total cohort of 400) (21). LADA patients were selected using common criteria: an initial diagnosis of type 2 diabetes (age at diagnosis 25eC68 years), documented antibody-positivity for GADA and/or IA-2, and no immediate requirement for insulin (within 3 months of diagnosis). The control group comprised normoglycemic, GADA/IA-2AeCnegative subjects recruited as part of the Diabetes in Families study (n = 332). Clinical characteristics of these subjects are given in Table 1. Requirement for insulin therapy in the UKPDS was protocol defined (14).
Genotyping.
In non-African populations, low haplotype diversity in the insulin gene region enables flanking SNPs to act as surrogate markers for minisatellite lineages (7). The eC23HphI SNP (rs689) denotes the class I/III subdivision; absence of the +3,580MspI restriction site (rs2000993) distinguishes the IDeC haplotype from other class I haplotypes (IC+/ID+), and the PH/VPH subdivision is detected through the +1,404Fnu4HI site (rs3842755) (7,9). All three SNPs were genotyped using the Amplifluor method (22). Based on extensive duplicate genotyping with a variety of methods including resequencing, we estimate an overall genotyping error of <0.5%. Genotyping success rates exceeded 95% in all cohorts. Primer sequences and additional genotyping details are available from the authors.
Statistical analyses.
Deviations from HWE were sought using exact methods implemented in PEDSTATS (v6.0; University of Michigan) (23). Homogeneity (between LADA groups) and case-control association testing was performed by standard contingency table methods implemented in StatXact 6 (Cytel Software, Cambridge, MA). For each SNP, we tested for association under an additive model and assuming dominance for the protective allele (T at eC23HphI, A at +1,404Fnu4HI, and C at +3,580MspI) given previous evidence from studies of type 1 diabetes (5,15,24). Step-wise logistic regression (SPSS, v13.0 for Windows; SPSS, Chicago, IL) was used to ascertain which SNPs independently contributed to risk. Measures of pairwise LD and 3-point haplotype frequencies were estimated using the expectation-maximization algorithm implemented in HelixTree Genetics Analysis software (Golden Helix, Bozeman, MT). Haplotype frequencies were compared using haplotype trend regression in HelixTree (25). For diplotype analyses, the most probable pair of haplotypes was assigned to each individual. Though such direct assignment underestimates the associated variance, the consequences are modest: low rates of missing data and high local LD mean that >97% of assigned haplotypes had posterior probabilities >0.99. Association of age at diagnosis with genotype was assessed by linear regression modeling using SPSS. The relationship between insulin requirement 6 years postdiagnosis and genotype was determined by contingency table analysis in StatXact 6. Power calculations were performed using Quanto, v1.0 (http://hydra.usc.edu/gxe). A P value <0.05 was deemed significant.
ACKNOWLEDGMENTS
We thank Diabetes U.K. for funding this work and the Wellcome Trust for funding A.C.
We acknowledge the laboratories of Professors Gian Franco Bottazzo (Royal London Hospital, London, U.K.) and Polly Bingley (University of Bristol, Bristol, U.K.) for performing antibody testing of study subjects.
FOOTNOTES
DOI: 10.2337/db06-0089
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
REFERENCES
Tuomi T, Groop LC, Zimmet PZ, Rowley MJ, Knowles W, Mackay IR: Antibodies to glutamic acid decarboxylase reveal latent autoimmune diabetes mellitus in adults with a noneCinsulin-dependent onset of disease. Diabetes 42:359eC362, 1993
Turner R, Stratton I, Horton V, Manley S, Zimmet P, Mackay IR, Shattock M, Bottazzo GF, Holman R: UKPDS 25: autoantibodies to islet-cell cytoplasm and glutamic acid decarboxylase for prediction of insulin requirement in type 2 diabetes: UK Prospective Diabetes Study Group. Lancet 350:1288eC1293, 1997
Zimmet PZ, Tuomi T, Mackay IR, Rowley MJ, Knowles W, Cohen M, Lang DA: Latent autoimmune diabetes mellitus in adults (LADA): the role of antibodies to glutamic acid decarboxylase in diagnosis and prediction of insulin dependency. Diabet Med 11:299eC303, 1994
Todd JA: Genetic analysis of type 1 diabetes using whole genome approaches. Proc Natl Acad Sci U S A 92:8560eC8565, 1995
Barratt BJ, Payne F, Lowe CE, Hermann R, Healy BC, Harold D, Concannon P, Gharani N, McCarthy MI, Olavesen MG, McCormack R, Guja C, Ionescu-Tirgoviste C, Undlien DE, Ronningen KS, Gillespie KM, Tuomilehto-Wolf E, Tuomilehto J, Bennett ST, Clayton DG, Cordell HJ, Todd JA: Remapping the insulin gene/IDDM2 locus in type 1 diabetes. Diabetes 53:1884eC1889, 2004
Kralovicova J, Gaunt TR, Rodriguez S, Wood PJ, Day IN, Vorechovsky I: Variants in the human insulin gene that affect pre-mRNA splicing: is eC23HphI a functional single nucleotide polymorphism at IDDM2 Diabetes 55:260eC264, 2006
Stead JD, Buard J, Todd JA, Jeffreys AJ: Influence of allele lineage on the role of the insulin minisatellite in susceptibility to type 1 diabetes. Hum Mol Genet 9:2929eC2935, 2000
Bennett ST, Todd JA: Human type 1 diabetes and the insulin gene: principles of mapping polygenes. Annu Rev Genet 30:343eC370, 1996
Bennett ST, Lucassen AM, Gough SC, Powell EE, Undlien DE, Pritchard LE, Merriman ME, Kawaguchi Y, Dronsfield MJ, Pociot F, et al.: Susceptibility to human type 1 diabetes at IDDM2 is determined by tandem repeat variation at the insulin gene minisatellite locus. Nat Genet 9:284eC292, 1995
Cerrone GE, Caputo M, Lopez AP, Gonzalez C, Massa C, Cedola N, Targovnik HM, Frechtel GD: Variable number of tandem repeats of the insulin gene determines susceptibility to latent autoimmune diabetes in adults. Mol Diagn 8:43eC49, 2004
Tuomi T, Carlsson A, Li H, Isomaa B, Miettinen A, Nilsson A, Nissen M, Ehrnstrom BO, Forsen B, Snickars B, Lahti K, Forsblom C, Saloranta C, Taskinen MR, Groop LC: Clinical and genetic characteristics of type 2 diabetes with and without GAD antibodies. Diabetes 48:150eC157, 1999
Fourlanos S, Dotta F, Greenbaum CJ, Palmer JP, Rolandsson O, Colman PG, Harrison LC: Latent autoimmune diabetes in adults (LADA) should be less latent. Diabetologia 48:2206eC2212, 2005
Palmer JP, Hampe CS, Chiu H, Goel A, Brooks-Worrell BM: Is latent autoimmune diabetes in adults distinct from type 1 diabetes or just type 1 diabetes at an older age Diabetes 54 (Suppl. 2):S62eCS67, 2005
UKPDS Group: UK Prospective Diabetes Study (UKPDS). VIII. Study design, progress and performance. Diabetologia 34:877eC890, 1991
Lucassen AM, Julier C, Beressi JP, Boitard C, Froguel P, Lathrop M, Bell JI: Susceptibility to insulin dependent diabetes mellitus maps to a 4.1 kb segment of DNA spanning the insulin gene and associated VNTR. Nat Genet 4:305eC310, 1993
Undlien DE, Bennett ST, Todd JA, Akselsen HE, Ikaheimo I, Reijonen H, Knip M, Thorsby E, Ronningen KS: Insulin gene region-encoded susceptibility to IDDM maps upstream of the insulin gene. Diabetes 44:620eC625, 1995
Pugliese A, Zeller M, Fernandez A, Jr, Zalcberg LJ, Bartlett RJ, Ricordi C, Pietropaolo M, Eisenbarth GS, Bennett ST, Patel DD: The insulin gene is transcribed in the human thymus and transcription levels correlated with allelic variation at the INS VNTR-IDDM2 susceptibility locus for type 1 diabetes. Nat Genet 15:293eC297, 1997
Castleden HAJ, Shields B, Bingley PJ, Williams AJK, Sampson M, Walker M, Gibson JM, McCarthy MI, Hitman GA, Levy JC, Hattersley AT, Vaidya B, Pearson ER: GAD antibodies in probands and their relatives in a cohort clinically selected for type 2 diabetes. Diabet Med In press
Frayling TM, Walker M, McCarthy MI, Evans JC, Allen LI, Lynn S, Ayres S, Millauer B, Turner C, Turner RC, Sampson MJ, Hitman GA, Ellard S, Hattersley AT: Parent-offspring trios: a resource to facilitate the identification of type 2 diabetes genes. Diabetes 48:2475eC2479, 1999
Wiltshire S, Hattersley AT, Hitman GA, Walker M, Levy JC, Sampson M, O’Rahilly S, Frayling TM, Bell JI, Lathrop GM, Bennett A, Dhillon R, Fletcher C, Groves CJ, Jones E, Prestwich P, Simecek N, Rao PV, Wishart M, Bottazzo GF, Foxon R, Howell S, Smedley D, Cardon LR, Menzel S, McCarthy MI: A genomewide scan for loci predisposing to type 2 diabetes in a U.K. population (the Diabetes UK Warren 2 Repository): analysis of 573 pedigrees provides independent replication of a susceptibility locus on chromosome 1q. Am J Hum Genet 69:553eC569, 2001
Owen KR, Stride A, Ellard S, Hattersley AT: Etiological investigation of diabetes in young adults presenting with apparent type 2 diabetes. Diabetes Care 26:2088eC2093, 2003
Bengra C, Mifflin TE, Khripin Y, Manunta P, Williams SM, Jose PA, Felder RA: Genotyping of essential hypertension single-nucleotide polymorphisms by a homogeneous PCR method with universal energy transfer primers. Clin Chem 48:2131eC2140, 2002
Wigginton JE, Cutler DJ, Abecasis GR: A note on exact tests of Hardy-Weinberg equilibrium. Am J Hum Genet 76:887eC893, 2005
Bennett ST, Wilson AJ, Cucca F, Nerup J, Pociot F, McKinney PA, Barnett AH, Bain SC, Todd JA: IDDM2-VNTR-encoded susceptibility to type 1 diabetes: dominant protection and parental transmission of alleles of the insulin gene-linked minisatellite locus. J Autoimmun 9:415eC421, 1996
Zaykin DV, Westfall PH, Young SS, Karnoub MA, Wagner MJ, Ehm MG: Testing association of statistically inferred haplotypes with discrete and continuous traits in samples of unrelated individuals. Hum Hered 53:79eC91, 2002(Minal Desai, Eleftheria Z)
2 Wellcome Trust Centre for Human Genetics, University of Oxford, Oxford, U.K.
3 Institute of Biomedical and Clinical Sciences, Peninsula Medical School, Exeter, U.K.
4 Centre for Diabetes and Metabolic Medicine, Bart’s and The London Queen Mary’s School of Medicine and Dentistry, London, U.K.
5 School of Clinical Medical Sciences, University of Newcastle, Newcastle, U.K.
6 Diabetes Trials Unit, Oxford Centre for Diabetes, Endocrinology and Metabolism, University of Oxford, Oxford, U.K.
HWE, Hardy-Weinburg equilibrium; LADA, latent autoimmune diabetes in adults; LD, linkage disequilibrium; PH, protective haplotype; SNP, single nucleotide polymorphism; UKPDS, U.K. Prospective Diabetes Study; VNTR, variable number of tandem repeats; VPH, very protective haplotype; W2 Repository, Warren 2 Repository; Exeter YT2D, Exeter Young-Onset type 2 Diabetes study
ABSTRACT
The etiopathological relationship between latent autoimmune diabetes in adults (LADA) and classical type 1 (insulin dependent) diabetes remains unclear. Variation at the insulin gene variable number of tandem repeats (VNTR) minisatellite influences susceptibility to type 1 diabetes, but studies in LADA have been small and inconsistent. We examined the role of insulin gene variation (using flanking variants as surrogates for VNTR subtypes) in the largest case-control study of LADA to date (400 case and 332 control subjects). Highly significant associations were identified with disease, with dominant protective effects of the T allele at eC23HphI (odds ratio [OR] 0.42 [95% CI 0.31eC0.58], P = 2.4 x 10eC8), A allele at +1,404Fnu4HI (0.50 [0.36eC0.70], P = 3.2 x 10eC5), and C allele at +3,580MspI (0.55 [0.35eC0.85], P = 0.0046). As with type 1 diabetes, the eC23HphI variant (a surrogate for the subdivision of VNTR into class I and III alleles) most clearly defined susceptibility in LADA. However, there was no association with age at diagnosis or requirement for insulin therapy 6 years postdiagnosis. This study establishes that variation within the insulin gene region does influence susceptibility to LADA, with the direction and magnitude of effect indistinguishable from that previously reported for type 1 diabetes. In conclusion, differences in VNTR-encoded susceptibility do not explain the differences in clinical presentation that distinguish classical type 1 diabetes and LADA.
Latent autoimmune diabetes in adults (LADA) and type 1 (insulin dependent) diabetes are both characterized by islet autoimmunity, indicated by the presence of circulating islet autoantibodies. However, the later age of onset (after the peak of type 1 diabetes incidence) and less aggressive clinical course (no immediate requirement for insulin therapy) of LADA result in substantial clinical overlap with type 2 diabetes (noneCinsulin dependent) (1eC3). Qualitative or quantitative differences in genetic susceptibility may underlie the marked variation in disease course and clinical presentation between LADA and type 1 diabetes.
The insulin gene region (IDDM2) accounts for 10% of familial risk for type 1 diabetes (4). Most evidence indicates that the variable number of tandem repeats (VNTR) minisatellite (which maps upstream of the translation initiation site) is the probable etiological variant, though certain flanking single nucleotide polymorphisms (SNPs) that are highly correlated with VNTR variation may also play a role (5,6). In non-African populations, much of the variation in VNTR length (repeat unit number) and composition (repeat element sequence) can be captured by dichotomizing VNTR alleles into two classes: the short class I alleles (26eC63 repeats) and longer class III alleles (141eC209 repeats) (7). In type 1 diabetes, class I alleles confer susceptibility to disease and class III, protection (8). Because of strong linkage disequilibrium (LD) between VNTR allele structure and haplotypes spanning the insulin gene, this major division in VNTR structure can be assayed using the eC23HphI SNP (9).
Further refinement of type 1 diabetes susceptibility is possible by examining the finer structure within major VNTR classes. A discrete subgroup of class I alleles, termed "IDeC," has been reported to confer protection rather than susceptibility (7). Additionally, class IIIA and IIIB alleles may differ in the degree of protection they offer, an observation reflected in the designation of the extended haplotypes on which these alleles lie as the protective haplotypes (PHs) and very protective haplotypes (VPHs), respectively (9). However, recent reexamination of these subclass associations in additional samples has cast significant doubt on these latter observations (5).
It is unclear whether these susceptibility effects of the insulin gene are evident in LADA and/or whether differences in VNTR-encoded susceptibility contribute to the clinical differences between LADA and type 1 diabetes. Published studies have yielded conflicting results (10,11). One small study from South America (n = 44 LADA and 138 control subjects) reported a significant class I association with LADA (10). However, another larger study, performed in 95 LADA case subjects and 172 control subjects from Finland, found no association of the eC23HphI site with disease (11). The present study uses the largest collection of patients yet studied to determine whether the IDDM2 locus is associated with LADA.
The LADA group (n = 400) comprised autoantibody-positive patients from three sources, the U.K. Prospective Diabetes Study (UKPDS), the Warren 2 Repository (W2 Repository), and the Exeter Young-Onset type 2 Diabetes study (Exeter YT2D); 332 normoglycemic autoantibody-negative subjects from the Diabetes in Families study served as control subjects (Table 1). All were genotyped for three SNPs known to be surrogate markers for minisatellite lineages: eC23HphI (A T), +1,404Fnu4HI (C A), and +3,580MspI (C T) (7,9).
Comparisons of genotype frequency distributions revealed borderline evidence for heterogeneity between the different case groups (eC23HphI, P = 0.045; +1,404Fnu4HI, P = 0.033; +3,580MspI, P = 0.062) (Table 2). Accordingly, results are presented for the case groups combined as well as for individual case groups with the control group. Genotype frequencies at all sites were in Hardy-Weinburg equilibrium (HWE) in the control subjects, but there were modest departures in some of the case groups (Exeter YT2D, eC23HphI, P = 0.008; W2 Repository, +1,404Fnu4HI, P = 0.01; UKPDS, +3,580MspI, P = 0.024). Only the last of these departures was in a direction consistent with the observed association effect. We attribute the other deviations to stochastic variation, since extensive regenotyping and resequencing (and failure to detect unexpected haplotypic combinations) effectively excluded genotyping error.
All three SNPs were associated with LADA in the combined group, with dominant protective effects of the T allele at eC23HphI (odds ratio [OR] 0.42 [95% CI 0.31eC0.58], P = 2.4 x 10eC8), A allele at +1,404Fnu4HI (0.50 [0.36eC0.70], P = 3.2 x 10eC5), and C allele at +3,580MspI (0.55 [0.35eC0.85], P = 0.0046). Patterns of associations were similar under the additive model (data not shown). Comparisons of the individual case groups with the control subjects (Table 3) confirmed that these associations were not an artifact of case-group heterogeneity: all case-control comparisons showed similar ORs, although sample size differences meant that not all such comparisons were significant. Step-wise logistic regression analysis in the combined LADA group versus control subjects identified the eC23HphI site (indicating VNTR class) as the primary disease-associated variant (Wald 2, P = 2.8 x 10eC8); the other sites had no additional independent effect.
Given recent discussion over the appropriate definition of LADA (12,13), we repeated these analyses, excluding those patients in the UKPDS (n = 45), W2 Repository (n = 4), and Exeter YT2D (n = 6) groups diagnosed under the age of 30 years and in whom the interval between diagnosis and insulin therapy was between 3 and 6 months. The strength of the association was unchanged in this analysis (OReC23HphI = 0.42 [95% CI 0.30eC0.58], P = 5.6 x 10eC8).
All three SNPs were, as expected, in strong LD: pairwise D' values exceeded 0.94, although r2 values ranged from 0.06 to 0.89. Given the modest departures from HWE in the case groups, the expectation-maximization algorithm was reinitialized three times to ensure convergence. This, together with presence of high LD between the SNPs, enabled accurate haplotype estimation. Four haplotypes with estimated frequencies exceeding 1% were observed, two corresponding to the subdivision of class I alleles (into IC+/ID+ and IDeC), the others to the major class III haplotypes PH and VPH (Table 4). Haplotype trend regression revealed highly significant differences in haplotype frequency distributions between the combined case subjects compared with the control subjects (P = 7.9 x 10eC6). Class I haplotypes were more frequent in the LADA than in the control subjects, in keeping with classical type 1 diabetes.
Diplotype frequencies differed significantly between case and control subjects (P < 10eC4 overall). The OR (95% CI) for VNTR diplotypes estimated under the recessive model for class I haplotypes and the dominant model for class III haplotypes were as follows: IDeC homozygotes versus all others 1.82 (1.18eC2.84), P = 0.0045; IC+/ID+ homozygotes versus all others 1.67 (1.10eC2.58), P = 0.016; PH-containing diplotypes versus all others 0.50 (0.36eC0.69), P = 2.0 x 10eC5; and VPH-containing diplotypes versus all others 0.49 (0.29eC0.80), P = 0.0028. Similar ORs were seen for the comparisons of the separate case groups against the control subjects (data not shown).
There was no association between age at diagnosis and insulin gene variation (class III positive mean ± SD 45.8 ± 10.1 years vs. class I homozygote 46.5 ± 10.5, P = 0.54). Within the UKPDS group (the only group for which such data were available—insulin requirement was determined by a protocol-driven assessment [14]), there was no association with insulin therapy requirement at 6 years postdiagnosis (class III positive 35% vs. class I homozygote 34%, P = 0.88).
This study has established that the INS gene is a susceptibility locus for autoimmune diabetes developing in adulthood as well as in children. The eC23HphI site (or, equivalently, VNTR class) appears to be the primary disease-associated variant, as in classical type 1 diabetes. We conclude that the different clinical features of LADA and type 1 diabetes cannot be explained by an absence of IDDM2-encoded genetic susceptibility in adult-onset patients (11). From the data available, there is no suggestion that the magnitude of the INS VNTR effect size in LADA differs appreciably from that seen in type 1 diabetes. In the combined LADA group, the OR (95% CI) associated with dominantly protective class III alleles was 0.42 (0.31eC0.58), a value comparable with that found in studies of classical type 1 diabetes (for example, OR 0.36 [95% CI 0.23eC0.55] and 0.22 [0.12eC0.39], as calculated from previous reports [9,15]).
Failure to detect a VNTR association with LADA in a previous study (11) may reflect several factors, including reduced power (a lower frequency of the class III allele in Finns (11,16) combined with a smaller sample size) and/or genuine differences arising through variation in genetic background or ascertainment.
Fine-structure analysis of the region demonstrated similarities in the association patterns of LADA compared with those reported for type 1 diabetes. Diplotype analyses showed that the IDeC haplotype confers similar susceptibility to LADA to other class I haplotypes and that the class III PH and VPH had similar protective effects. Both are consistent with findings from a recent reevaluation of the insulin gene in type 1 diabetes (5), although the wide CIs indicate the need for cautious interpretation. The present study had >90% power to detect the dominantly protective association of class III alleles with OR 0.5 (for = 0.05). However, considerably larger numbers of subjects would be required to detect subtle differences in predisposition associated with specific class I or class III subgroups; the largest comparison of the protective effects associated with PH and VPH in type 1 diabetes featured meta-analysis of data from 3,000 pedigrees (5).
We have made no corrections for multiple testing; application of a conservative Bonferroni correction to the lowest P value obtained (P = 2.4 x 10eC8) would still retain significance (P < 10eC4), indicating that the observed associations are likely to be real.
Class III alleles have been associated with a decreased level of transcription of INS in the thymus during development (17). It has been postulated that this difference may affect central tolerance mechanisms, promoting negative selection of autoreactive T-cells and decreasing the propensity to develop autoimmunity against insulin (17). Similar protective mechanisms are likely to operate in LADA. However, as yet unresolved differences in class I subgroup susceptibility (or, alternatively, within HLA or other type 1 diabeteseCsusceptibility loci) could, through their consequences on the development of tolerance, contribute to the differences in clinical phenotype between classical type 1 diabetes and LADA.
RESEARCH DESIGN AND METHODS
All subjects were unrelated and of exclusively British/Irish Europid origin. LADA subjects (n = 400) came from three sources of type 2 diabeteseCdiagnosed patients: the UKPDS (n = 232 from the total cohort of 5,102) (14), the W2 Repository (n = 131 from the total sample of 4,000) (18,19,20), and the Exeter YT2D study (n = 37 from the total cohort of 400) (21). LADA patients were selected using common criteria: an initial diagnosis of type 2 diabetes (age at diagnosis 25eC68 years), documented antibody-positivity for GADA and/or IA-2, and no immediate requirement for insulin (within 3 months of diagnosis). The control group comprised normoglycemic, GADA/IA-2AeCnegative subjects recruited as part of the Diabetes in Families study (n = 332). Clinical characteristics of these subjects are given in Table 1. Requirement for insulin therapy in the UKPDS was protocol defined (14).
Genotyping.
In non-African populations, low haplotype diversity in the insulin gene region enables flanking SNPs to act as surrogate markers for minisatellite lineages (7). The eC23HphI SNP (rs689) denotes the class I/III subdivision; absence of the +3,580MspI restriction site (rs2000993) distinguishes the IDeC haplotype from other class I haplotypes (IC+/ID+), and the PH/VPH subdivision is detected through the +1,404Fnu4HI site (rs3842755) (7,9). All three SNPs were genotyped using the Amplifluor method (22). Based on extensive duplicate genotyping with a variety of methods including resequencing, we estimate an overall genotyping error of <0.5%. Genotyping success rates exceeded 95% in all cohorts. Primer sequences and additional genotyping details are available from the authors.
Statistical analyses.
Deviations from HWE were sought using exact methods implemented in PEDSTATS (v6.0; University of Michigan) (23). Homogeneity (between LADA groups) and case-control association testing was performed by standard contingency table methods implemented in StatXact 6 (Cytel Software, Cambridge, MA). For each SNP, we tested for association under an additive model and assuming dominance for the protective allele (T at eC23HphI, A at +1,404Fnu4HI, and C at +3,580MspI) given previous evidence from studies of type 1 diabetes (5,15,24). Step-wise logistic regression (SPSS, v13.0 for Windows; SPSS, Chicago, IL) was used to ascertain which SNPs independently contributed to risk. Measures of pairwise LD and 3-point haplotype frequencies were estimated using the expectation-maximization algorithm implemented in HelixTree Genetics Analysis software (Golden Helix, Bozeman, MT). Haplotype frequencies were compared using haplotype trend regression in HelixTree (25). For diplotype analyses, the most probable pair of haplotypes was assigned to each individual. Though such direct assignment underestimates the associated variance, the consequences are modest: low rates of missing data and high local LD mean that >97% of assigned haplotypes had posterior probabilities >0.99. Association of age at diagnosis with genotype was assessed by linear regression modeling using SPSS. The relationship between insulin requirement 6 years postdiagnosis and genotype was determined by contingency table analysis in StatXact 6. Power calculations were performed using Quanto, v1.0 (http://hydra.usc.edu/gxe). A P value <0.05 was deemed significant.
ACKNOWLEDGMENTS
We thank Diabetes U.K. for funding this work and the Wellcome Trust for funding A.C.
We acknowledge the laboratories of Professors Gian Franco Bottazzo (Royal London Hospital, London, U.K.) and Polly Bingley (University of Bristol, Bristol, U.K.) for performing antibody testing of study subjects.
FOOTNOTES
DOI: 10.2337/db06-0089
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
REFERENCES
Tuomi T, Groop LC, Zimmet PZ, Rowley MJ, Knowles W, Mackay IR: Antibodies to glutamic acid decarboxylase reveal latent autoimmune diabetes mellitus in adults with a noneCinsulin-dependent onset of disease. Diabetes 42:359eC362, 1993
Turner R, Stratton I, Horton V, Manley S, Zimmet P, Mackay IR, Shattock M, Bottazzo GF, Holman R: UKPDS 25: autoantibodies to islet-cell cytoplasm and glutamic acid decarboxylase for prediction of insulin requirement in type 2 diabetes: UK Prospective Diabetes Study Group. Lancet 350:1288eC1293, 1997
Zimmet PZ, Tuomi T, Mackay IR, Rowley MJ, Knowles W, Cohen M, Lang DA: Latent autoimmune diabetes mellitus in adults (LADA): the role of antibodies to glutamic acid decarboxylase in diagnosis and prediction of insulin dependency. Diabet Med 11:299eC303, 1994
Todd JA: Genetic analysis of type 1 diabetes using whole genome approaches. Proc Natl Acad Sci U S A 92:8560eC8565, 1995
Barratt BJ, Payne F, Lowe CE, Hermann R, Healy BC, Harold D, Concannon P, Gharani N, McCarthy MI, Olavesen MG, McCormack R, Guja C, Ionescu-Tirgoviste C, Undlien DE, Ronningen KS, Gillespie KM, Tuomilehto-Wolf E, Tuomilehto J, Bennett ST, Clayton DG, Cordell HJ, Todd JA: Remapping the insulin gene/IDDM2 locus in type 1 diabetes. Diabetes 53:1884eC1889, 2004
Kralovicova J, Gaunt TR, Rodriguez S, Wood PJ, Day IN, Vorechovsky I: Variants in the human insulin gene that affect pre-mRNA splicing: is eC23HphI a functional single nucleotide polymorphism at IDDM2 Diabetes 55:260eC264, 2006
Stead JD, Buard J, Todd JA, Jeffreys AJ: Influence of allele lineage on the role of the insulin minisatellite in susceptibility to type 1 diabetes. Hum Mol Genet 9:2929eC2935, 2000
Bennett ST, Todd JA: Human type 1 diabetes and the insulin gene: principles of mapping polygenes. Annu Rev Genet 30:343eC370, 1996
Bennett ST, Lucassen AM, Gough SC, Powell EE, Undlien DE, Pritchard LE, Merriman ME, Kawaguchi Y, Dronsfield MJ, Pociot F, et al.: Susceptibility to human type 1 diabetes at IDDM2 is determined by tandem repeat variation at the insulin gene minisatellite locus. Nat Genet 9:284eC292, 1995
Cerrone GE, Caputo M, Lopez AP, Gonzalez C, Massa C, Cedola N, Targovnik HM, Frechtel GD: Variable number of tandem repeats of the insulin gene determines susceptibility to latent autoimmune diabetes in adults. Mol Diagn 8:43eC49, 2004
Tuomi T, Carlsson A, Li H, Isomaa B, Miettinen A, Nilsson A, Nissen M, Ehrnstrom BO, Forsen B, Snickars B, Lahti K, Forsblom C, Saloranta C, Taskinen MR, Groop LC: Clinical and genetic characteristics of type 2 diabetes with and without GAD antibodies. Diabetes 48:150eC157, 1999
Fourlanos S, Dotta F, Greenbaum CJ, Palmer JP, Rolandsson O, Colman PG, Harrison LC: Latent autoimmune diabetes in adults (LADA) should be less latent. Diabetologia 48:2206eC2212, 2005
Palmer JP, Hampe CS, Chiu H, Goel A, Brooks-Worrell BM: Is latent autoimmune diabetes in adults distinct from type 1 diabetes or just type 1 diabetes at an older age Diabetes 54 (Suppl. 2):S62eCS67, 2005
UKPDS Group: UK Prospective Diabetes Study (UKPDS). VIII. Study design, progress and performance. Diabetologia 34:877eC890, 1991
Lucassen AM, Julier C, Beressi JP, Boitard C, Froguel P, Lathrop M, Bell JI: Susceptibility to insulin dependent diabetes mellitus maps to a 4.1 kb segment of DNA spanning the insulin gene and associated VNTR. Nat Genet 4:305eC310, 1993
Undlien DE, Bennett ST, Todd JA, Akselsen HE, Ikaheimo I, Reijonen H, Knip M, Thorsby E, Ronningen KS: Insulin gene region-encoded susceptibility to IDDM maps upstream of the insulin gene. Diabetes 44:620eC625, 1995
Pugliese A, Zeller M, Fernandez A, Jr, Zalcberg LJ, Bartlett RJ, Ricordi C, Pietropaolo M, Eisenbarth GS, Bennett ST, Patel DD: The insulin gene is transcribed in the human thymus and transcription levels correlated with allelic variation at the INS VNTR-IDDM2 susceptibility locus for type 1 diabetes. Nat Genet 15:293eC297, 1997
Castleden HAJ, Shields B, Bingley PJ, Williams AJK, Sampson M, Walker M, Gibson JM, McCarthy MI, Hitman GA, Levy JC, Hattersley AT, Vaidya B, Pearson ER: GAD antibodies in probands and their relatives in a cohort clinically selected for type 2 diabetes. Diabet Med In press
Frayling TM, Walker M, McCarthy MI, Evans JC, Allen LI, Lynn S, Ayres S, Millauer B, Turner C, Turner RC, Sampson MJ, Hitman GA, Ellard S, Hattersley AT: Parent-offspring trios: a resource to facilitate the identification of type 2 diabetes genes. Diabetes 48:2475eC2479, 1999
Wiltshire S, Hattersley AT, Hitman GA, Walker M, Levy JC, Sampson M, O’Rahilly S, Frayling TM, Bell JI, Lathrop GM, Bennett A, Dhillon R, Fletcher C, Groves CJ, Jones E, Prestwich P, Simecek N, Rao PV, Wishart M, Bottazzo GF, Foxon R, Howell S, Smedley D, Cardon LR, Menzel S, McCarthy MI: A genomewide scan for loci predisposing to type 2 diabetes in a U.K. population (the Diabetes UK Warren 2 Repository): analysis of 573 pedigrees provides independent replication of a susceptibility locus on chromosome 1q. Am J Hum Genet 69:553eC569, 2001
Owen KR, Stride A, Ellard S, Hattersley AT: Etiological investigation of diabetes in young adults presenting with apparent type 2 diabetes. Diabetes Care 26:2088eC2093, 2003
Bengra C, Mifflin TE, Khripin Y, Manunta P, Williams SM, Jose PA, Felder RA: Genotyping of essential hypertension single-nucleotide polymorphisms by a homogeneous PCR method with universal energy transfer primers. Clin Chem 48:2131eC2140, 2002
Wigginton JE, Cutler DJ, Abecasis GR: A note on exact tests of Hardy-Weinberg equilibrium. Am J Hum Genet 76:887eC893, 2005
Bennett ST, Wilson AJ, Cucca F, Nerup J, Pociot F, McKinney PA, Barnett AH, Bain SC, Todd JA: IDDM2-VNTR-encoded susceptibility to type 1 diabetes: dominant protection and parental transmission of alleles of the insulin gene-linked minisatellite locus. J Autoimmun 9:415eC421, 1996
Zaykin DV, Westfall PH, Young SS, Karnoub MA, Wagner MJ, Ehm MG: Testing association of statistically inferred haplotypes with discrete and continuous traits in samples of unrelated individuals. Hum Hered 53:79eC91, 2002(Minal Desai, Eleftheria Z)