D6S265*15 marks a DRB1*15, DQB1*0602 haplotype associated with attenuated protection from type 1 diabetes mellitus

Association between alleles, haplotypes, and amino acid variations in HLA class II genes and type 1 diabetes in Kuwaiti children

Type 1 diabetes (T1D) is a complex autoimmune disorder that is highly prevalent globally. The interactions between genetic and environmental factors may trigger T1D in susceptible individuals. HLA genes play a significant role in T1D pathogenesis, and specific haplotypes are associated with an increased risk of developing the disease. Identifying risk haplotypes can greatly improve the genetic scoring for early diagnosis of T1D in difficult to rank subgroups. This study employed next-generation sequencing to evaluate the association between HLA class II alleles, haplotypes, and amino acids and T1D, by recruiting 95 children with T1D and 150 controls in the Kuwaiti population. Significant associations were identified for alleles at the HLA-DRB1, HLA-DQA1, and HLA-DQB1 loci, including DRB1*03:01:01, DQA1*05:01:01, and DQB1*02:01:01, which conferred high risk, and DRB1*11:04:01, DQA1*05:05:01, and DQB1*03:01:01, which were protective. The DRB1*03:01:01~DQA1*05:01:01~DQB1*02:01:01 haplotype was most strongly associated with the risk of developing T1D, while DRB1*11:04-DQA1*05:05-DQB1*03:01 was the only haplotype that rendered protection against T1D. We also identified 66 amino acid positions across the HLA-DRB1, HLA-DQA1, and HLA-DQB1 genes that were significantly associated with T1D, including novel associations. These results validate and extend our knowledge on the associations between HLA genes and T1D in Kuwaiti children. The identified risk alleles, haplotypes, and amino acid variations may influence disease development through effects on HLA structure and function and may allow early intervention via population-based screening efforts.

DR15-DQ6 remains dominantly protective against type 1 diabetes throughout the first five decades of life

AIMS/HYPOTHESIS

Among white European children developing type 1 diabetes, the otherwise common HLA haplotype DR15-DQ6 is rare, and highly protective. Adult-onset type 1 diabetes is now known to represent more overall cases than childhood onset, but it is not known whether DR15-DQ6 is protective in older-adult-onset type 1 diabetes. We sought to quantify DR15-DQ6 protection against type 1 diabetes as age of onset increased.

METHODS

In two independent cohorts we assessed the proportion of type 1 diabetes cases presenting through the first 50 years of life with DR15-DQ6, compared with population controls. In the After Diabetes Diagnosis Research Support System-2 (ADDRESS-2) cohort (n = 1458) clinician-diagnosed type 1 diabetes was confirmed by positivity for one or more islet-specific autoantibodies. In UK Biobank (n = 2502), we estimated type 1 diabetes incidence rates relative to baseline HLA risk for each HLA group using Poisson regression. Analyses were restricted to white Europeans and were performed in three groups according to age at type 1 diabetes onset: 0-18 years, 19-30 years and 31-50 years.

RESULTS

DR15-DQ6 was protective against type 1 diabetes through to age 50 years (OR < 1 for each age group, all p < 0.001). The following ORs for type 1 diabetes, relative to a neutral HLA genotype, were observed in ADDRESS-2: age 5-18 years OR 0.16 (95% CI 0.08, 0.31); age 19-30 years OR 0.10 (0.04, 0.23); and age 31-50 years OR 0.37 (0.21, 0.68). DR15-DQ6 also remained highly protective at all ages in UK Biobank. Without DR15-DQ6, the presence of major type 1 diabetes high-risk haplotype (either DR3-DQ2 or DR4-DQ8) was associated with increased risk of type 1 diabetes.

CONCLUSIONS/INTERPRETATION

HLA DR15-DQ6 confers dominant protection from type 1 diabetes across the first five decades of life.

Dog leucocyte antigen (DLA) class II haplotypes and risk of canine diabetes mellitus in specific dog breeds

Background

Canine diabetes mellitus (DM) is a common endocrine disease in domestic dogs. A number of pathological mechanisms are thought to contribute to the aetiopathogenesis of relative or absolute insulin deficiency, including immune-mediated destruction of pancreatic beta cells. DM risk varies considerably between different dog breeds, suggesting that genetic factors are involved and contribute susceptibility or protection. Associations of particular dog leucocyte antigen (DLA) class II haplotypes with DM have been identified, but investigations to date have only considered all breeds pooled together. The aim of this study was to analyse an expanded data set so as to identify breed-specific diabetes-associated DLA haplotypes.

Methods

The 12 most highly represented breeds in the UK Canine Diabetes Register were selected for study. DLA-typing data from 646 diabetic dogs and 912 breed-matched non-diabetic controls were analysed to enable breed-specific analysis of the DLA. Dogs were genotyped for allelic variation at DLA-DRB1, -DQA1, -DQB1 loci using DNA sequence-based typing. Genotypes from all three loci were combined to reveal three-locus DLA class II haplotypes, which were evaluated for statistical associations with DM. This was performed for each breed individually and for all breeds pooled together.

Results

Five dog breeds were identified as having one or more DLA haplotype associated with DM susceptibility or protection. Four DM-associated haplotypes were identified in the Cocker Spaniel breed, of which one haplotype was shared with Border Terriers. In the three breeds known to be at highest risk of DM included in the study (Samoyed, Tibetan Terrier and Cairn Terrier), no DLA haplotypes were found to be associated with DM.

Conclusions

Novel DLA associations with DM in specific dog breeds provide further evidence that immune response genes contribute susceptibility to this disease in some cases. It is also apparent that DLA may not be contributing obvious or strong risk for DM in some breeds, including the seven breeds analysed for which no associations were identified.

Early prediction of autoimmune (type 1) diabetes

Underlying type 1 diabetes is a genetic aetiology dominated by the influence of specific HLA haplotypes involving primarily the class II DR-DQ region. In genetically predisposed children with the DR4-DQ8 haplotype, exogenous factors, yet to be identified, are thought to trigger an autoimmune reaction against insulin, signalled by insulin autoantibodies as the first autoantibody to appear. In children with the DR3-DQ2 haplotype, the triggering reaction is primarily against GAD signalled by GAD autoantibodies (GADA) as the first-appearing autoantibody. The incidence rate of insulin autoantibodies as the first-appearing autoantibody peaks during the first years of life and declines thereafter. The incidence rate of GADA as the first-appearing autoantibody peaks later but does not decline. The first autoantibody may variably be followed, in an apparently non-HLA-associated pathogenesis, by a second, third or fourth autoantibody. Although not all persons with a single type of autoantibody progress to diabetes, the presence of multiple autoantibodies seems invariably to be followed by loss of functional beta cell mass and eventually by dysglycaemia and symptoms. Infiltration of mononuclear cells in and around the islets appears to be a late phenomenon appearing in the multiple-autoantibody-positive with dysglycaemia. As our understanding of the aetiology and pathogenesis of type 1 diabetes advances, the improved capability for early prediction should guide new strategies for the prevention of type 1 diabetes.

Genetic risk analysis of a patient with fulminant autoimmune type 1 diabetes mellitus secondary to combination ipilimumab and nivolumab immunotherapy

BACKGROUND

Checkpoint inhibitor immunotherapy is becoming an effective treatment modality for an increasing number of malignancies. As a result, autoinflammatory side-effects are also being observed more commonly in the clinic. We are currently unable to predict which patients will develop more severe toxicities associated with these treatment regimens.

CASE PRESENTATION

We present a patient with stage IV melanoma that developed rapid onset autoimmune type 1 diabetes (T1D) in response to combination ipilimumab and nivolumab immunotherapy. At the time of the patient's presentation with diabetes ketoacidosis, a confirmed anti-GAD antibody seroconversion was noted. Longer-term follow-up of this patient has demonstrated a durable complete response based on PET CT imaging along with a persistently undetectable C-peptide level. Single nucleotide polymorphism gene sequencing and HLA risk allele analysis has revealed the patient to lack any established genetic predisposition to the development of autoimmune T1D.

CONCLUSIONS

While larger studies are necessary to better understand the role of genetic risk factors for the development of autoimmune toxicities in those patients undergoing checkpoint inhibitor immunotherapy, these results suggest that pre-screening patients for known T1D risk alleles may not be indicated. Additional investigation is needed to determine whether an approach such as T cell receptor clonotypic analysis to identify the presence of autoreactive T cell clones may be an effective approach for predicting which patients are at risk for the development of autoinflammatory toxicities while undergoing checkpoint inhibitor immunotherapy.

HLA-DRB1*15:01-DQA1*01:02-DQB1*06:02 Haplotype Protects Autoantibody-Positive Relatives From Type 1 Diabetes Throughout the Stages of Disease Progression

The HLA-DRB1*15:01-DQA1*01:02-DQB1*06:02 haplotype is linked to protection from the development of type 1 diabetes (T1D). However, it is not known at which stages in the natural history of T1D development this haplotype affords protection. We examined a cohort of 3,358 autoantibody-positive relatives of T1D patients in the Pathway to Prevention (PTP) Study of the Type 1 Diabetes TrialNet. The PTP study examines risk factors for T1D and disease progression in relatives. HLA typing revealed that 155 relatives carried this protective haplotype. A comparison with 60 autoantibody-negative relatives suggested protection from autoantibody development. Moreover, the relatives with DRB1*15:01-DQA1*01:02-DQB1*06:02 less frequently expressed autoantibodies associated with higher T1D risk, were less likely to have multiple autoantibodies at baseline, and rarely converted from single to multiple autoantibody positivity on follow-up. These relatives also had lower frequencies of metabolic abnormalities at baseline and exhibited no overall metabolic worsening on follow-up. Ultimately, they had a very low 5-year cumulative incidence of T1D. In conclusion, the protective influence of DRB1*15:01-DQA1*01:02-DQB1*06:02 spans from autoantibody development through all stages of progression, and relatives with this allele only rarely develop T1D.

Genetics of the HLA region in the prediction of type 1 diabetes

Type 1 diabetes (T1D) is one of the most widely studied complex genetic disorders, and the genes in HLA are reported to account for approximately 40-50% of the familial aggregation of T1D. The major genetic determinants of this disease are polymorphisms of class II HLA genes encoding DQ and DR. The DR-DQ haplotypes conferring the highest risk are DRB1*03:01-DQA1*05:01-DQB1*02:01 (abbreviated "DR3") and DRB1*04:01/02/04/05/08-DQA1*03:01-DQB1*03:02/04 (or DQB1*02; abbreviated "DR4"). The risk is much higher for the heterozygote formed by these two haplotypes (OR = 16.59; 95% CI, 13.7-20.1) than for either of the homozygotes (DR3/DR3, OR = 6.32; 95% CI, 5.12-7.80; DR4/DR4, OR = 5.68; 95% CI, 3.91). In addition, some haplotypes confer strong protection from disease, such as DRB1*15:01-DQA1*01:02-DQB1*06:02 (abbreviated "DR2"; OR = 0.03; 95% CI, 0.01-0.07). After adjusting for the genetic correlation with DR and DQ, significant associations can be seen for HLA class II DPB1 alleles, in particular, DPB1*04:02, DPB1*03:01, and DPB1*02:02. Outside of the class II region, the strongest susceptibility is conferred by class I allele B*39:06 (OR =10.31; 95% CI, 4.21-25.1) and other HLA-B alleles. In addition, several loci in the class III region are reported to be associated with T1D, as are some loci telomeric to class I. Not surprisingly, current approaches for the prediction of T1D in screening studies take advantage of genotyping HLA-DR and HLA-DQ loci, which is then combined with family history and screening for autoantibodies directed against islet-cell antigens. Inclusion of additional moderate HLA risk haplotypes may help identify the majority of children with T1D before the onset of the disease.

Genetic variation within the HLA class III influences T1D susceptibility conferred by high-risk HLA haplotypes

Human leukocyte antigen (HLA) class II DRB1 and DQB1 represent the major type I diabetes (T1D) genetic susceptibility loci; however, other genes in the HLA region are also involved in T1D risk. We analyzed 1411 pedigrees (2865 affected individuals) from the type I diabetes genetics consortium genotyped for HLA classical loci and for 12 single-nucleotide polymorphisms (SNPs) in the class III region previously shown to be associated with T1D in a subset of 886 pedigrees. Using the transmission disequilibrium test, we compared the proportion of SNP alleles transmitted from within the high-risk DR3 and DR4 haplotypes to affected offspring. Markers rs4151659 (mapping to CFB) and rs7762619 (mapping 5' of LTA) were the most strongly associated with T1D on DR3 (P=1.2 x 10(-9) and P=2 x 10(-12), respectively) and DR4 (P=4 x 10(-15) and P=8 x 10(-8), respectively) haplotypes. They remained significantly associated after stratifying individuals in analyses for B*1801, A*0101-B*0801, DPB1*0301, DPB1*0202, DPB1*0401 or DPB1*0402. Rs7762619 and rs4151659 are in strong linkage disequilibrium (LD) (r(2)=0.82) with each other, but a joint analysis showed that the association for each SNP was not solely because of LD. Our data support a role for more than one locus in the class III region contributing to risk of T1D.

Several loci in the HLA class III region are associated with T1D risk after adjusting for DRB1-DQB1

AIM

Several studies have indicated that genes in the human leucocyte antigen (HLA) region additional to the HLA class II DRB1-DQB1 contribute to type 1 diabetes (T1D) susceptibility. The aim of this study was to assess if markers in the class III Major Histocompatibility Complex (MHC) region are associated with T1D after accounting for linkage disequilibrium (LD) with DRB1-DQB1.

METHODS

We investigated 356 single nucleotide polymorphisms (SNPs) in the class III region covering 1.1 megabases in two subsets of data: 289 Human Biological Data Interchange (HBDI) Caucasian families and 597 additional Caucasian families collected by the Type 1 Diabetes Genetics Consortium (T1DGC). Analysis conditioning on DRB1-DQB1 was performed using the overall conditional genotype method.

RESULTS

Thirteen SNPs replicated in both subsets of the data and showed evidence of an additional effect on disease risk. Although some of the SNPs are in tight LD with each other, at least six of the associations were not because of LD with other class III markers. The strongest association within class III markers was with rs2395106 that maps 5' to the NOTCH4 gene, which has also been implicated in susceptibility to rheumatoid arthritis. The second association was with rs707915 mapping to the MSH5 gene, in a block of six markers significantly associated with T1D after adjusting for LD with DR-DQ. In total, six-independent associations within class III were observed although results were not adjusted for LD with class I.

CONCLUSIONS

Our data confirm that the class III region is involved in T1D susceptibility and suggest that more than one gene in the region is involved.

Extreme genetic risk for type 1A diabetes in the post-genome era

A series of genes and loci influencing the genetic risk of type 1A (immune-mediated) diabetes are now well characterized. These include genes of the major histocompatibility complex (MHC), polymorphisms 5' of the insulin gene, and PTPN22, as well as more recently defined loci from genome-wide association studies. By far the major determinants of risk for type 1A diabetes are genes within or linked to the MHC and in particular alleles of class II genes (HLA-DR, DQ, and DP). There is evidence that MHC class I alleles contribute and there are additional MHC-linked influences such that for a major subset of relatives of patients there is a risk as high as 80% for siblings, and for the general population a risk as high as 20% can be defined at birth just by analyzing the MHC. We believe the search for additional MHC loci will require analysis of the remarkable long-range identity (up to 9 million base pairs) of extended MHC haplotypes. Current prediction algorithms will likely be greatly improved for the general population when the additional contributing loci of the MHC are defined.

HLA DR-DQ haplotypes and genotypes and type 1 diabetes risk: analysis of the type 1 diabetes genetics consortium families

OBJECTIVE

The Type 1 Diabetes Genetics Consortium has collected type 1 diabetic families worldwide for genetic analysis. The major genetic determinants of type 1 diabetes are alleles at the HLA-DRB1 and DQB1 loci, with both susceptible and protective DR-DQ haplotypes present in all human populations. The aim of this study is to estimate the risk conferred by specific DR-DQ haplotypes and genotypes.

RESEARCH DESIGN AND METHODS

Six hundred and seven Caucasian families and 38 Asian families were typed at high resolution for the DRB1, DQA1, and DQB1 loci. The association analysis was performed by comparing the frequency of DR-DQ haplotypes among the chromosomes transmitted to an affected child with the frequency of chromosomes not transmitted to any affected child.

RESULTS

A number of susceptible, neutral, and protective DR-DQ haplotypes have been identified, and a statistically significant hierarchy of type 1 diabetes risk has been established. The most susceptible haplotypes are the DRB1*0301-DQA1*0501-DQB1*0201 (odds ratio [OR] 3.64) and the DRB1*0405-DQA1*0301-DQB1*0302, DRB1*0401-DQA1*0301-DQB*0302, and DRB1*0402-DQA1*0301-DQB1*0302 haplotypes (ORs 11.37, 8.39, and 3.63), followed by the DRB1*0404-DQA1*0301-DQB1*0302 (OR 1.59) and the DRB1*0801-DQB1*0401-DQB1*0402 (OR 1.25) haplotypes. The most protective haplotypes are DRB1*1501-DQA1*0102-DQB1*0602 (OR 0.03), DRB1*1401-DQA1*0101-DQB1*0503 (OR 0.02), and DRB1*0701-DQA1*0201-DQB1*0303 (OR 0.02).

CONCLUSIONS

Specific combinations of alleles at the DRB1, DQA1, and DQB1 loci determine the extent of haplotypic risk. The comparison of closely related DR-DQ haplotype pairs with different type 1 diabetes risks allowed identification of specific amino acid positions critical in determining disease susceptibility. These data also indicate that the risk associated with specific HLA haplotypes can be influenced by the genotype context and that the trans-complementing heterodimer encoded by DQA1*0501 and DQB1*0302 confers very high risk.

Searching for additional disease loci in a genomic region

Our aim is to review methods to optimize detection of all disease genes in a genetic region. As a starting point, we assume there is sufficient evidence from linkage and/or association studies, based on significance levels or replication studies, for the involvement in disease risk of the genetic region under study. For closely linked markers, there will often be multiple associations with disease, and linkage analyses identify a region rather than the specific disease-predisposing gene. Hence, the first task is to identify the primary (major) disease-predisposing gene or genes in a genetic region, and single nucleotide polymorphisms thereof, that is, how to distinguish true associations from those that are just due to linkage disequilibrium with the actual disease-predisposing variants. Then, how do we detect additional disease genes in this genetic region? These two issues are of course very closely interrelated. No existing programs, either individually or in aggregate, can handle the magnitude and complexity of the analyses needed using currently available methods. Further, even with modern computers, one cannot study every possible combination of genetic markers and their haplotypes across the genome, or even within a genetic region. Although we must rely heavily on computers, in the final analysis of multiple effects in a genetic region and/or interaction or independent effects between unlinked genes, manipulation of the data by the individual investigator will play a crucial role. We recommend a multistrategy approach using a variety of complementary methods described below.

Relative predispositional effects of HLA class II DRB1-DQB1 haplotypes and genotypes on type 1 diabetes: a meta-analysis

The direct involvement of the human leukocyte antigen class II DR-DQ genes in type 1 diabetes (T1D) is well established, and these genes display a complex hierarchy of risk effects at the genotype and haplotype levels. We investigated, using data from 38 studies, whether the DR-DQ haplotypes and genotypes show the same relative predispositional effects across populations and ethnic groups. Significant differences in risk within a population were considered, as well as comparisons across populations using the patient/control (P/C) ratio. Within a population, the ratio of the P/C ratios for two different genotypes or haplotypes is a function only of the absolute penetrance values, allowing ranking of risk effects. Categories of consistent predisposing, intermediate ('neutral'), and protective haplotypes were identified and found to correlate with disease prevalence and the marked ethnic differences in DRB1-DQB1 frequencies. Specific effects were identified, for example for predisposing haplotypes, there was a statistically significant and consistent hierarchy for DR4 DQB1*0302s: DRB1*0405 =*0401 =*0402 > *0404 > *0403, with DRB1*0301 DQB1*0200 (DR3) being significantly less predisposing than DRB1*0402 and more than DRB1*0404. The predisposing DRB1*0401 DQB1*0302 haplotype was relatively increased compared with the protective haplotype DRB1*0401 DQB1*0301 in heterozygotes with DR3 compared with heterozygotes with DRB1*0101 DQB1*0501 (DR1). Our results show that meta-analyses and use of the P/C ratio and rankings thereof can be valuable in determining T1D risk factors at the haplotype and amino acid residue levels.

Sequence variation within the major histocompatibility complex subregion centromeric of HLA class II in type 1 diabetes

The extended major histocompatibility complex (xMHC) has been studied intensively with regard to type 1 diabetes (T1D) predisposition. So far, little attention has been given to the subregion centromeric of MHC class II. We selected five single nucleotide polymorphisms in genes with potential immune-related functions in the genomic regions of death-domain-associated protein 6 (DAXX, apoptosis associated), TAP-binding protein (TAPBP, human leukocyte antigen class I loading) and retinoic acid receptor beta (RXRB, vitamin D receptor function) that may bear relevance to the pathogenesis of T1D. A total of 277 unrelated individuals with juvenile-onset T1D and 286 control subjects were genotyped using sequence-specific priming-polymerase chain reaction. The genotype and allelic frequencies of the markers tested were not significantly different between patients and control subjects. Subsequent haplotype analysis showed six DAXX-TAPBP-RXRB haplotypic configurations. No difference was observed between patients and control cohorts when stratified for T1D high-risk DQ2-DR17 and DQ8-DR4 haplotypes. However, the distribution of these haplotypes affected T1D susceptibility encoded by the intermediate risk haplotypes DQ5-DR1 and DQ2-DR7 by increasing and decreasing susceptibility, respectively. We propose that studying genetic variants in the xMHC may be particularly rewarding to define disease pathways in patients displaying intermediate risk DQ-DR haplotypes.

14th International HLA and Immunogenetics Workshop: report on the HLA component of type 1 diabetes

The type 1 diabetes (T1D) component of the 13th International Histocompatibility Workshop (IHW) obtained microsatellite (msat) and human leukocyte antigen (HLA)-DR/DQ data on case/control and family samples through an international collaboration. The aim was to detect the effects of susceptibility loci on the HLA complex independent of the primary determinants in the class II region (HLA-DR/DQ). As part of the activity of the 14th International HLA and Immunogenetics Workshop (14th IHIWS), a T1D workshop was held to present analyses of the 13th IHW data and to discuss the current status of knowledge about the genetics of T1D. These data are now available online through dbMHC, a web-based resource established by the National Center for Biotechnology. Continuing work since the 13th IHW has resulted in published work showing heterogeneity of DR3 haplotypes in data sets from the 13th IHW and Human Biological Data Interchange (HBDI). In addition, we identified markers that define DRB1*1501 DQB1*0602 haplotypes conferring reduced protection from diabetes in a Swedish 13th IHW data set. Further analyses of the 13th IHW data set not only showed some significant results but also demonstrated extensive heterogeneity reminiscent of non-HLA genes. The haplotype analysis in HBDI families identified two msats with significant effects on susceptibility and statistically significant age of onset effects at class III markers that are not because of linkage disequilibrium, with class I alleles known to affect age of onset. The above studies underscore the importance of refining our understanding of susceptibility associated with genes in the HLA complex.

Genetic determinants of type 1 diabetes across populations

T1D results from autoimmune-mediated destruction of the pancreatic beta cells, a process that is conditioned by multiple genes and environmental factors. The main genetic determinants map to the major histocompatibility complex (MHC), and in particular DR and DQ, although, genes outside the MHC contribute, including the insulin gene, PTPN22, and CTLA-4. There are remarkable differences in genetic susceptibility to T1D between populations. We believe this variation reflects differing frequencies of diabetes causative and protective alleles and haplotypes, and thus remains a major genetic influence linked to the MHC region not accounted for by DR and DQ alleles. In this article, we discuss global variations in genetic susceptibility to T1D in view of current genetic understanding.

An immunologic homunculus for type 1 diabetes

Autoimmune diseases such as the diabetes that develops in NOD mice depend on immunologic recognition of specific autoantigens, but recognition can result in a pathogenic or protective T cell response. A study by Du et al. in this issue of the JCI demonstrates that TGF-beta signaling by T cells recognizing the insulin peptide B:9-23 is essential for such protection and that this inhibitory cytokine functions in both a paracrine and an autocrine manner (see the related article beginning on page 1360). We propose that the insulin peptide B:9-23 and a conserved TCR motif form an "immunologic homunculus" underlying the relatively common targeting of insulin by T cells that, as demonstrated by the study of Du and coworkers, results in a protective T cell response, or diabetes, as shown by other investigators, for related T cell receptors.