Search

HLA-DQA1 and -DQB1 genes have significant and potentially causal associations with autoimmune type 1 diabetes (T1D). To follow up on the earlier analysis on high-risk HLA-DQ2.5 and DQ8.1, the current analysis uncovers seven residues (αa1, α157, α196, β9, β30, β57, and β70) that are resistant to T1D among subjects with DQ4-, 5-, 6-, and 7-resistant DQ haplotypes. These 7 residues form 13 common motifs: 6 motifs are significantly resistant, 6 motifs have modest or no associations (P values >0.05), and 1 motif has 7 copies observed among control subjects only. The motifs “DAAFYDG,” “DAAYHDG,” and “DAAYYDR” have significant resistance to T1D (odds ratios [ORs] 0.03, 0.25, and 0.18; P = 6.11 × 10−24, 3.54 × 10−15, and 1.03 × 10−21, respectively). Remarkably, a change of a single residue from the motif “DAAYHDG” to “DAAYHSG” (D to S at β57) alters the resistance potential, from resistant motif (OR 0.15; P = 3.54 × 10−15) to a neutral motif (P = 0.183), the change of which was significant (Fisher P value = 0.0065). The extended set of linked residues associated with T1D resistance and unique to each cluster of HLA-DQ haplotypes represents facets of all known features and functions of these molecules: antigenic peptide binding, peptide–MHC class II complex stability, β167-169 RGD loop, T-cell receptor binding, formation of homodimer of α-β heterodimers, and cholesterol binding in the cell membrane rafts. Identification of these residues is a novel understanding of resistant DQ associations with T1D. Our analyses endow potential molecular approaches to identify immunological mechanisms that control disease susceptibility or resistance to provide novel targets for immunotherapeutic strategies.

HLA-DQA1 and -DQB1 are strongly associated with type 1 diabetes (T1D), and DQ8.1 and DQ2.5 are major risk haplotypes. Next-generation targeted sequencing of HLA-DQA1 and -DQB1 in Swedish newly diagnosed 1- to 18 year-old patients (n = 962) and control subjects (n = 636) was used to construct abbreviated DQ haplotypes, converted into amino acid (AA) residues, and assessed for their associations with T1D. A hierarchically organized haplotype (HOH) association analysis allowed 45 unique DQ haplotypes to be categorized into seven clusters. The DQ8/9 cluster included two DQ8.1 risk and the DQ9 resistant haplotypes, and the DQ2 cluster included the DQ2.5 risk and DQ2.2 resistant haplotypes. Within each cluster, HOH found residues α44Q (odds ratio [OR] 3.29, P = 2.38 * 10−85) and β57A (OR 3.44, P = 3.80 * 10−84) to be associated with T1D in the DQ8/9 cluster representing all ten residues (α22, α23, α44, α49, α51, α53, α54, α73, α184, β57) due to complete linkage disequilibrium (LD) of α44 with eight such residues. Within the DQ2 cluster and due to LD, HOH analysis found α44C and β135D to share the risk for T1D (OR 2.10, P = 1.96 * 10−20). The motif “QAD” of α44, β57, and β135 captured the T1D risk association of DQ8.1 (OR 3.44, P = 3.80 * 10−84), and the corresponding motif “CAD” captured the risk association of DQ2.5 (OR 2.10, P = 1.96 * 10−20). Two risk associations were related to GAD65 autoantibody (GADA) and IA-2 autoantibody (IA-2A) but in opposite directions. CAD was positively associated with GADA (OR 1.56, P = 6.35 * 10−8) but negatively with IA-2A (OR 0.59, P = 6.55 * 10−11). QAD was negatively associated with GADA (OR 0.88; P = 3.70 * 10−3) but positively with IA-2A (OR 1.64; P = 2.40 * 10−14), despite a single difference at α44. The residues are found in and around anchor pockets 1 and 9, as potential T-cell receptor contacts, in the areas for CD4 binding and putative homodimer formation. The identification of three HLA-DQ AAs (α44, β57, β135) conferring T1D risk should sharpen functional and translational studies.

Next-generation targeted sequencing of HLA-DRB1 and HLA-DRB3, -DRB4, and -DRB5 (abbreviated as DRB345) provides high resolution of functional variant positions to investigate their associations with type 1 diabetes risk and with autoantibodies against insulin (IAA), GAD65 (GADA), IA-2 (IA-2A), and ZnT8 (ZnT8A). To overcome exceptional DR sequence complexity as a result of high polymorphisms and extended linkage disequilibrium among the DR loci, we applied a novel recursive organizer (ROR) to discover disease-associated amino acid residues. ROR distills disease-associated DR sequences and identifies 11 residues of DRB1, sequences of which retain all significant associations observed by DR genes. Furthermore, all 11 residues locate under/adjoining the peptide-binding groove of DRB1, suggesting a plausible functional mechanism through peptide binding. The 15 residues of DRB345, located respectively in the β49–55 homodimerization patch and on the face of the molecule shown to interact with and bind to the accessory molecule CD4, retain their significant disease associations. Further ROR analysis of DR associations with autoantibodies finds that DRB1 residues significantly associated with ZnT8A and DRB345 residues with GADA. The strongest association is between four residues (β14, β25, β71, and β73) and IA-2A, in which the sequence ERKA confers a risk association (odds ratio 2.15, P = 10−18), and another sequence, ERKG, confers a protective association (odds ratio 0.59, P = 10−11), despite a difference of only one amino acid. Because motifs of identified residues capture potentially causal DR associations with type 1 diabetes, this list of residuals is expected to include corresponding causal residues in this study population.

The possible contribution of HLA-DRB3, -DRB4, and -DRB5 alleles to type 1 diabetes risk and to insulin autoantibody (IAA), GAD65 (GAD autoantibody [GADA]), IA-2 antigen (IA-2A), or ZnT8 against either of the three amino acid variants R, W, or Q at position 325 (ZnT8RA, ZnT8WA, and ZnT8QA, respectively) at clinical diagnosis is unclear. Next-generation sequencing (NGS) was used to determine all DRB alleles in consecutively diagnosed patients ages 1–18 years with islet autoantibody–positive type 1 diabetes (n = 970) and control subjects (n = 448). DRB3, DRB4, or DRB5 alleles were tested for an association with the risk of DRB1 for autoantibodies, type 1 diabetes, or both. The association between type 1 diabetes and DRB1*03:01:01 was affected by DRB3*01:01:02 and DRB3*02:02:01. These DRB3 alleles were associated positively with GADA but negatively with ZnT8WA, IA-2A, and IAA. The negative association between type 1 diabetes and DRB1*13:01:01 was affected by DRB3*01:01:02 to increase the risk and by DRB3*02:02:01 to maintain a negative association. DRB4*01:03:01 was strongly associated with type 1 diabetes (P = 10−36), yet its association was extensively affected by DRB1 alleles from protective (DRB1*04:03:01) to high (DRB1*04:01:01) risk, but its association with DRB1*04:05:01 decreased the risk. HLA-DRB3, -DRB4, and -DRB5 affect type 1 diabetes risk and islet autoantibodies. HLA typing with NGS should prove useful to select participants for prevention or intervention trials.