Motifs of three HLA-DQ amino acid residues (α44, β57, β135) capture full association with the risk of type 1 diabetes in DQ2 and DQ8 children
figureposted on 2020-04-03, 22:52 authored by Ada AdminAda Admin, Lue Ping Zhao, George K Papadopoulos, William W. Kwok, Antonis K. Moustakas, George P. Bondinas, Helena Elding Larsson, Johnny Ludvigsson, Claude Marcus, Ulf Samuelsson, Ruihan Wang, Chul-Woo Pyo, Wyatt C. Nelson, Daniel E. Geraghty, Åke Lernmark
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-18 year-old patients (n=962) and controls (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 (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), 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 GADA and 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 TCR contacts, in the areas for CD4 binding and putative homodimer formation. The identification of three HLA-DQ AA (α44, β57, β135) conferring T1D risk should sharpen functional and translational studies.