HLA Class II (DR, DQ, DP) Genes Were Separately Associated With the Progression From Seroconversion to Onset of Type 1 Diabetes Among Participants in Two Diabetes Prevention Trials (DPT-1 and TN07).

OBJECTIVE: To explore associations of HLA class II genes (HLAII) with the progression of islet autoimmunity from asymptomatic to symptomatic type 1 diabetes (T1D). RESEARCH DESIGN AND METHODS: Next-generation targeted sequencing was used to genotype eight HLAII genes (DQA1, DQB1, DRB1, DRB3, DRB4, DRB5, DPA1, DPB1) in 1,216 participants from the Diabetes Prevention Trial-1 and Randomized Diabetes Prevention Trial with Oral Insulin sponsored by TrialNet. By the linkage disequilibrium, DQA1 and DQB1 are haplotyped to form DQ haplotypes; DP and DR haplotypes are similarly constructed. Together with available clinical covariables, we applied the Cox regression model to assess HLAII immunogenic associations with the disease progression. RESULTS: First, the current investigation updated the previously reported genetic associations of DQA1*03:01-DQB1*03:02 (hazard ratio [HR] = 1.25, P = 3.50*10-3) and DQA1*03:03-DQB1*03:01 (HR = 0.56, P = 1.16*10-3), and also uncovered a risk association with DQA1*05:01-DQB1*02:01 (HR = 1.19, P = 0.041). Second, after adjusting for DQ, DPA1*02:01-DPB1*11:01 and DPA1*01:03-DPB1*03:01 were found to have opposite associations with progression (HR = 1.98 and 0.70, P = 0.021 and 6.16*10-3, respectively). Third, DRB1*03:01-DRB3*01:01 and DRB1*03:01-DRB3*02:02, sharing the DRB1*03:01, had opposite associations (HR = 0.73 and 1.44, P = 0.04 and 0.019, respectively), indicating a role of DRB3. Meanwhile, DRB1*12:01-DRB3*02:02 and DRB1*01:03 alone were found to associate with progression (HR = 2.6 and 2.32, P = 0.018 and 0.039, respectively). Fourth, through enumerating all heterodimers, it was found that both DQ and DP could exhibit associations with disease progression. CONCLUSIONS: These results suggest that HLAII polymorphisms influence progression from islet autoimmunity to T1D among at-risk subjects with islet autoantibodies.

Autoimmune susceptible HLA class II motifs facilitate the presentation of modified neoepitopes to potentially autoreactive T cells.

Rheumatoid arthritis (RA), multiple sclerosis (MS), type 1 diabetes (T1D), and celiac disease (CD), are strongly associated with susceptible HLA class II haplotypes. The peptide-binding pockets of these molecules are polymorphic, thus each HLA class II protein presents a distinct set of peptides to CD4+ T cells. Peptide diversity is increased through post-translational modifications, generating non-templated sequences that enhance HLA binding and/or T cell recognition. The high-risk HLA-DR alleles that confer susceptibility to RA are notable for their ability to accommodate citrulline, promoting responses to citrullinated self-antigens. Likewise, HLA-DQ alleles associated with T1D and CD favor the binding of deamidated peptides. In this review, we discuss structural features that promote modified self-epitope presentation, provide evidence supporting the relevance of T cell recognition of such antigens in disease processes, and make a case that interrupting the pathways that generate such epitopes and reprogramming neoepitope-specific T cells are key strategies for effective therapeutic intervention.

A Simulation Study of the Effect of Naturally Occurring Point Mutations on the SRY-DNA Complex.

Molecular dynamics (MD) simulations were conducted in order to investigate the effect of the naturally occurring point mutations of the transcription factor (TF) sex-determining region Y (SRY) on the structure and dynamics of the SRY-DNA complex. The normal SRY, along with the two mutants I13T and G40R, comprising point mutations on the SRY chain, which have been clinically identified in patients with sex developmental disorders, were modeled as DNA complexes. Our modeling work aims at elucidating atomic-level structural determinants of the aberrant SRY-DNA complexation by means of µs-long MD. The results suggest that the observed disorders brought about by the G40R-DNA and I13T-DNA may arise predominantly from the destabilization of the complex being in accord with in vitro assays found elsewhere and from modifications of the DNA bending as revealed in this study. Comparative potential of mean force computations, over a sequence of short separation distances for the three complexes, verified a higher stability of the normal SRY-DNA. Examining the way the SRY mutations modulate the SRY-DNA complex dynamics at the microscopic level is important also toward elucidating molecular determinants of function for proteins capable of binding to DNA.

Association of HLA-DQ Heterodimer Residues -18ß and ß57 With Progression From Islet Autoimmunity to Diabetes in the Diabetes Prevention Trial-Type 1.

OBJECTIVE: The purpose was to test the hypothesis that the HLA-DQαß heterodimer structure is related to the progression of islet autoimmunity from asymptomatic to symptomatic type 1 diabetes (T1D). RESEARCH DESIGN AND METHODS: Next-generation targeted sequencing was used to genotype HLA-DQA1-B1 class II genes in 670 subjects in the Diabetes Prevention Trial-Type 1 (DPT-1). Coding sequences were translated into DQ α- and ß-chain amino acid residues and used in hierarchically organized haplotype (HOH) association analysis to identify motifs associated with diabetes onset. RESULTS: The opposite diabetes risks were confirmed for HLA DQA1*03:01-B1*03:02 (hazard ratio [HR] 1.36; P = 2.01 ∗ 10-3) and DQA1*03:03-B1*03:01 (HR 0.62; P = 0.037). The HOH analysis uncovered residue -18ß in the signal peptide and ß57 in the ß-chain to form six motifs. DQ*VA was associated with faster (HR 1.49; P = 6.36 ∗ 10-4) and DQ*AD with slower (HR 0.64; P = 0.020) progression to diabetes onset. VA/VA, representing DQA1*03:01-B1*03:02 (DQ8/8), had a greater HR of 1.98 (P = 2.80 ∗ 10-3). The DQ*VA motif was associated with both islet cell antibodies (P = 0.023) and insulin autoantibodies (IAAs) (P = 3.34 ∗ 10-3), while the DQ*AD motif was associated with a decreased IAA frequency (P = 0.015). Subjects with DQ*VA and DQ*AD experienced, respectively, increasing and decreasing trends of HbA1c levels throughout the follow-up. CONCLUSIONS: HLA-DQ structural motifs appear to modulate progression from islet autoimmunity to diabetes among at-risk relatives with islet autoantibodies. Residue -18ß within the signal peptide may be related to levels of protein synthesis and ß57 to stability of the peptide-DQab trimolecular complex.

The Effect of Cold Periods on the Biological Cycle of Marchalina hellenica.

Climate change is considered a major factor affecting honeybees' behavior and productivity with major consequences in both honey and agricultural production. Many research studies have expressed serious concerns about the mass losses of bee colonies and the role of bees as pollinators, while others have underlined important issues for the impact of the increase in temperature on honeybee abundance and honey yields. In the present work, we draw our attention to Marchalina hellenica, which is the most important honeydew-producing insect in Greece. A statistically significant forecasting model for the effect of cold periods in February on the life cycle of the insect is constructed, with the aid of the Cumulative Logit Model and the theory of runs. The forecasting model may help beekeepers plan the timely exploitation of honeydew secretions of pine trees, which will be beneficial for beekeepers, the rural economy, and forest protection. The new suggested model also indicates that, in view of the climate change scenarios seen in the literature, the life cycle of M. hellenica is expected to be drastically shorter.

A Localized Enantioselective Catalytic Site on Short DNA Sequences and Their Amphiphiles.

A DNA-based artificial metalloenzyme (ArM) consisting of a copper(II) complex of 4,4'-dimethyl-2,2'-bipyridine (dmbipy-Cu) bound to double-stranded DNA (dsDNA) as short as 8 base pairs with only 2 contiguous central pairs (G for guanine and C for cytosine) catalyzes the highly enantioselective Diels-Alder reaction, Michael addition, and Friedel-Crafts alkylation in water. Molecular simulations indicate that these minimal sequences provide a single site where dmbipy-Cu is groove-bound and able to function as an enantioselective catalyst. Enantioselective preference inverts when d-DNA is replaced with l-DNA. When the DNA is conjugated to a hydrophobic tail, the obtained ArMs exhibit enantioselective performance in a methanol-water mixture superior to that of non-amphiphilic dsDNA, and dsDNA-amphiphiles with more complex G•C-rich sequences.

Chemometric Study of Fatty Acid Composition of Virgin Olive Oil from Four Widespread Greek Cultivars.

Virgin olive oil (VOO) is one of the key components of the Mediterranean diet owing to the presence of monounsaturated fatty acids and various bioactive compounds. These beneficial traits, which are usually associated with the cultivar genotype, are highlighting the demand of identifying characteristics of olive oil that will ensure its authenticity. In this work, the fatty acid (FA) composition of 199 VOO samples from Koroneiki, Megaritiki, Amfissis, and Manaki cultivars was determined and studied by chemometrics. Olive cultivar greatly influenced the FA composition, namely, oleic acid (from 75.36% for Amfissis to 65.81% for Megaritiki) and linoleic acid (from 13.35% for Manaki to 6.70% for Koroneiki). Spearman's rho correlation coefficients revealed differences and similarities among the olive oil cultivars. The use of the forward stepwise algorithm identified the FAs arachidonic acid, gadoleic acid, linoleic acid, α-linolenic acid, palmitoleic acid, and palmitic acid as the most significant for the differentiation of samples. The application of linear and quadratic cross-validation discriminant analysis resulted in the correct classification of 100.00% and 99.37% of samples, respectively. The findings demonstrated the special characteristics of the VOO samples derived from the four cultivars and their successful botanical differentiation based on FA composition.

Salinity Stress Alters the Secondary Metabolic Profile of M. sativa, M. arborea and Their Hybrid (Alborea).

Increased soil salinity, and therefore accumulation of ions, is one of the major abiotic stresses of cultivated plants that negatively affect their growth and yield. Among Medicago species, only Medicago truncatula, which is a model plant, has been extensively studied, while research regarding salinity responses of two important forage legumes of Medicago sativa (M. sativa) and Medicago arborea (M. arborea) has been limited. In the present work, differences between M. arborea, M. sativa and their hybrid Alborea were studied regarding growth parameters and metabolomic responses. The entries were subjected to three different treatments: (1) no NaCl application (control plants), (2) continuous application of 100 mM NaCl (acute stress) and (3) gradual application of NaCl at concentrations of 50-75-150 mM by increasing NaCl concentration every 10 days. According to the results, M. arborea maintained steady growth in all three treatments and appeared to be more resistant to salinity. Furthermore, results clearly demonstrated that M. arborea presented a different metabolic profile from that of M. sativa and their hybrid. In general, it was found that under acute and gradual stress, M. sativa overexpressed saponins in the shoots while M. arborea overexpressed saponins in the roots, which is the part of the plant where most of the saponins are produced and overexpressed. Alborea did not perform well, as more metabolites were downregulated than upregulated when subjected to salinity stress. Finally, saponins and hydroxycinnamic acids were key players of increased salinity tolerance.

The KAG motif of HLA-DRB1 (ß71, ß74, ß86) predicts seroconversion and development of type 1 diabetes.

BACKGROUND: HLA-DR4, a common antigen of HLA-DRB1, has multiple subtypes that are strongly associated with risk of type 1 diabetes (T1D); however, some are risk neutral or resistant. The pathobiological mechanism of HLA-DR4 subtypes remains to be elucidated. METHODS: We used a population-based case-control study of T1D (962 patients and 636 controls) to decipher genetic associations of HLA-DR4 subtypes and specific residues with susceptibility to T1D. Using a birth cohort of 7865 children with periodically measured islet autoantibodies (GADA, IAA or IA-2A), we proposed to validate discovered genetic associations with a totally different study design and time-to-seroconversions prior to clinical onset of T1D. A novel analytic strategy hierarchically organized the HLA-DRB1 alleles by sequence similarity and identified critical amino acid residues by minimizing local genomic architecture and higher-order interactions. FINDINGS: Three amino acid residues of HLA-DRB1 (ß71, ß74, ß86) were found to be predictive of T1D risk in the population-based study. The "KAG" motif, corresponding to HLA-DRB1×04:01, was most strongly associated with T1D risk ([O]dds [R]atio=3.64, p = 3.19 × 10-64). Three less frequent motifs ("EAV", OR = 2.55, p = 0.025; "RAG", OR = 1.93, p = 0.043; and "RAV", OR = 1.56, p = 0.003) were associated with T1D risk, while two motifs ("REG" and "REV") were equally protective (OR = 0.11, p = 4.23 × 10-4). In an independent birth cohort of HLA-DR3 and HLA-DR4 subjects, those having the "KAG" motif had increased risk for time-to-seroconversion (Hazard Ratio = 1.74, p = 6.51 × 10-14) after adjusting potential confounders. INTERPRETATIONS: DNA sequence variation in HLA-DRB1 at positions ß71, ß74, and ß86 are non-conservative (ß74 A→E, ß71 E vs K vs R and ß86 G vs V). They result in substantial differences in peptide antigen anchor pocket preferences at p1, p4 and potentially neighboring regions such as pocket p7. Differential peptide antigen binding is likely to be affected. These sequence substitutions may account for most of the HLA-DR4 contribution to T1D risk as illustrated in two HLA-peptide model complexes of the T1D autoantigens preproinsulin and GAD65. FUNDING: National Institute of Diabetes and Digestive and Kidney Diseases and the Swedish Child Diabetes Foundation and the Swedish Research Council.

Nine residues in HLA-DQ molecules determine with susceptibility and resistance to type 1 diabetes among young children in Sweden.

HLA-DQ molecules account over 50% genetic risk of type 1 diabetes (T1D), but little is known about associated residues. Through next generation targeted sequencing technology and deep learning of DQ residue sequences, the aim was to uncover critical residues and their motifs associated with T1D. Our analysis uncovered (αa1, α44, α157, α196) and (ß9, ß30, ß57, ß70, ß135) on the HLA-DQ molecule. Their motifs captured all known susceptibility and resistant T1D associations. Three motifs, "DCAA-YSARD" (OR = 2.10, p = 1.96*10-20), "DQAA-YYARD" (OR = 3.34, 2.69*10-72) and "DQDA-YYARD" (OR = 3.71, 1.53*10-6) corresponding to DQ2.5 and DQ8.1 (the latter two motifs) associated with susceptibility. Ten motifs were significantly associated with resistance to T1D. Collectively, homozygous DQ risk motifs accounted for 43% of DQ-T1D risk, while homozygous DQ resistant motifs accounted for 25% protection to DQ-T1D risk. Of the identified nine residues five were within or near anchoring pockets of the antigenic peptide (α44, ß9, ß30, ß57 and ß70), one was the N-terminal of the alpha chain (αa1), one in the CD4-binding region (ß135), one in the putative cognate TCR-induced αß homodimerization process (α157), and one in the intra-membrane domain of the alpha chain (α196). Finding these critical residues should allow investigations of fundamental properties of host immunity that underlie tolerance to self and organ-specific autoimmunity.

Discrimination of botanical origin of olive oil from selected Greek cultivars by SPME-GC-MS and ATR-FTIR spectroscopy combined with chemometrics.

BACKGROUND: Consumers today wish to know the botanical origin of the olive oil they purchase. The objective of the present study was the development of robust chemometric models based on gas chromatography-mass spectrometry (GC-MS) and attenuated total reflectance-Fourier transform infrared spectroscopy (ATR-FTIR) for the purpose of botanical differentiation of three commercial Greek olive oil cultivars. RESULTS: Using the solid-phase microextraction technique (SPME), volatile compounds (VC) were obtained and analyzed by GC-MS. Five hydrocarbons and one ester were selected by the forward stepwise algorithm, which best discriminated the olive oil samples. From ATR-FTIR analysis, the spectral regions chosen from the forward stepwise algorithm were associated with CO stretching vibration of the esters of triglycerides and the CH bending vibrations of the CH2 aliphatic group and double bonds. Application of the supervised methods of linear and quadratic discriminant cross-validation analysis, based on VC data, provided a correct classification score of 97.4% and 100.0%, respectively. Corresponding statistical analyses were used in the mid-infrared spectra, by which 96.1% of samples were discriminated correctly. CONCLUSION: ATR-FTIR and SPME-GC-MS techniques in conjunction with the appropriate feature selection algorithm and classification methods proved to be powerful tools for the authentication of Greek olive oil. The proposed methodology could be used in an industrial setting for determination of the botanical origin of Greek olive oil. © 2020 Society of Chemical Industry.

Next-Generation HLA Sequence Analysis Uncovers Seven HLA-DQ Amino Acid Residues and Six Motifs Resistant to Childhood Type 1 Diabetes.

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.

Binding Dynamics of siRNA with Selected Lipopeptides: A Computer-Aided Study of the Effect of Lipopeptides' Functional Groups and Stereoisomerism.

The engineering issues pertaining to nanoparticle systems toward targeted gene therapies have not been fully probed. Recent experiments have identified specific structural characteristics of a novel class of lipopeptides (LP) that may lead to potent nanocarriers intended as RNAi therapeutics, albeit the molecular mechanism that underlies their performance remains unexplored. We conducted molecular dynamics simulations in atomistic detail coupled with free energy computations to study the dynamics and thermodynamics of an acrylate- and an epoxide-derived LP, members of the aforesaid class, upon their binding to siRNA in aqueous solution aiming at examining structure-potency relations. We found that the entropic part of the free energy of binding predominates; moreover, the first LP class tends to disrupt the Watson-Crick base pairing of siRNA, whereas the latter leaves the double helix intact. Moreover, the identified tug-of-war effect between LP-water and LP-siRNA hydrogen bonding in the supramolecular complex can underpin synthesis routes toward tuning the association dynamics. Our simulations on two diastereomers of the epoxide-derived LP showed significant structural and energetics differences upon binding, as a result of steric effects imposed by the different absolute configurations at their chiral centers. These findings may serve as crucial design parameters toward modulating the interplay between complex stability and ease of releasing the nucleic acid drug into the cell.

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.

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.

A modified flow cytometry method for objective estimation of human CD4+ regulatory T cells (CD4+ Tregs) in peripheral blood, via CD4/CD25/CD45RO/FoxP3 labeling.

BACKGROUND: Several methods exist for flow-cytometric estimation of human peripheral blood CD4+ T regulatory cells (CD4+ Tregs). METHODS: We report our experience with the estimation of human CD4+ Tregs via three different characterizations using flow cytometry (CD25high FoxP3+ , CD25high CD127low/- FoxP3+ , and CD4+ CD25high/int CD45ROFoxP3+ ) in normal subjects. We have used these methods on the control populations from two studies (32 and 36 subjects, respectively), the latter two methods retrospectively on the subjects of the first study. The six CD4+ T cell fractions obtained by the third method were differentially colored to ascertain the distribution of these cell fractions in the CD25/FoxP3, CD45RO/FoxP3, and CD25/CD127 dot plots from CD4/CD25/CD45RO/FoxP3 and CD4/CD25/CD45RO/CD127 panels. RESULTS: Each approach gives significantly different estimates of Tregs (expressed as percentage of CD4+ T cells), with the second almost invariably yielding higher percentages than the other two. Only the third approach can distinguish among effector and naïve Tregs and FoxP3+ non-Tregs. Analysis of CD25/CD127 dot plots reveals that Treg delineation via the widely used definition of CD4+ CD25high CD127low/- cells unavoidably yields a mixture of nearly all effector and most of naïve Tregs, as well as FoxP3+ non-Tregs plus other cells. Delineation of effector/naïve Tregs and FoxP3+ non-Tregs is possible via CD45RO/CD25 dot plots but not by CD45RO/FoxP3 counterparts (as done previously) because of overlapping FoxP3 intensities among Tregs and non-Tregs. CONCLUSION: Our comparison shows that CD4/CD25/CD45RO/FoxP3 panels are an objective means of estimating effector and naïve Tregs via colored dot plots, aiding thus in Treg delineation in health and detecting aberrations in disease.

Discriminative T cell recognition of cross-reactive islet-antigens is associated with HLA-DQ8 transdimer-mediated autoimmune diabetes.

Human leukocyte antigen (HLA)-DQ8 transdimer (HLA-DQA1*0501/DQB1*0302) confers exceptionally high risk in autoimmune diabetes. However, little is known about HLA-DQ8 transdimer-restricted CD4 T cell recognition, an event crucial for triggering HLA-DQ8 transdimer-specific anti-islet immunity. Here, we report a high degree of epitope overlap and T cell promiscuity between susceptible HLA-DQ8 and HLA-DQ8 transdimer. Despite preservation of putative residues for T cell receptor (TCR) contact, stronger disease-associated responses to cross-reactive, immunodominant islet epitopes are elicited by HLA-DQ8 transdimer. Mutagenesis at the α chain of HLA-DQ8 transdimer in complex with the disease-relevant GAD65250-266 peptide and in silico analysis reveal the DQ α52 residue located within the N-terminal edge of the peptide-binding cleft for the enhanced T cell reactivity, altering avidity and biophysical affinity between TCR and HLA-peptide complexes. Accordingly, a structurally promiscuous but nondegenerate TCR-HLA-peptide interface is pivotal for HLA-DQ8 transdimer-mediated autoimmune diabetes.

Eleven Amino Acids of HLA-DRB1 and Fifteen Amino Acids of HLA-DRB3, 4, and 5 Include Potentially Causal Residues Responsible for the Risk of Childhood Type 1 Diabetes.

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.

Epitope Stealing as a Mechanism of Dominant Protection by HLA-DQ6 in Type 1 Diabetes.

The heterozygous DQ2/8 (DQA1*05:01-DQB1*02:01/DQA1*03:01-DQB1*03:02) genotype confers the highest risk in type 1 diabetes (T1D), whereas the DQ6/8 (DQA1*02:01-DQB1*06:02/DQA1*03:01-DQB1*03:02) genotype is protective. The mechanism of dominant protection by DQ6 (DQB1*06:02) is unknown. We tested the hypothesis that DQ6 interferes with peptide binding to DQ8 by competition for islet epitope ("epitope stealing") by analysis of the islet ligandome presented by HLA-DQ6/8 and -DQ8/8 on dendritic cells pulsed with islet autoantigens preproinsulin (PPI), GAD65, and IA-2, followed by competition assays using a newly established "epitope-stealing" HLA/peptide-binding assay. HLA-DQ ligandome analysis revealed a distinct DQ6 peptide-binding motif compared with the susceptible DQ2/8 molecules. PPI and IA-2 peptides were identified from DQ6, of DQ6/8 heterozygous dendritic cells, but no DQ8 islet peptides were retrieved. Insulin B6-23, a highly immunogenic CD4 T-cell epitope in patients with T1D, bound to both DQ6 and DQ8. Yet, binding of InsB6-23 to DQ8 was prevented by DQ6. We obtained first functional evidence of a mechanism of dominant protection from disease, in which HLA molecules associated with protection bind islet epitopes in a different, competing, HLA-binding register, leading to "epitope stealing" and conceivably diverting the immune response from islet epitopes presented by disease-susceptible HLA molecules in the absence of protective HLA.

DRB4*01:01 Has a Distinct Motif and Presents a Proinsulin Epitope That Is Recognized in Subjects with Type 1 Diabetes.

DRB4*01:01 (DRB4) is a secondary HLA-DR product that is part of the high-risk DR4/DQ8 haplotype that is associated with type 1 diabetes (T1D). DRB4 shares considerable homology with HLA-DR4 alleles that predispose to autoimmunity, including DRB1*04:01 and DRB1*04:04. However, the DRB4 protein sequence includes distinct residues that would be expected to alter the characteristics of its binding pockets. To identify high-affinity peptides that are recognized in the context of DRB4, we used an HLA class II tetramer-based approach to identify epitopes within multiple viral Ags. We applied a similar approach to identify antigenic sequences within glutamic acid decarboxylase 65 and pre-proinsulin that are recognized in the context of DRB4. Seven sequences were immunogenic, eliciting high-affinity T cell responses in DRB4+ subjects. DRB1*04:01-restricted responses toward many of these peptides have been previously described, but responses to a novel pre-proinsulin 9-28 peptide were commonly observed in subjects with T1D. Furthermore, T cells that recognized this peptide in the context of DRB4 were present at significantly higher frequencies in patients with T1D than in healthy controls, implicating this as a disease-relevant specificity that may contribute to the breakdown of ß cell tolerance in genetically susceptible individuals. We then deduced a DRB4 motif and confirmed its key features through structural modeling. This modeling suggested that the core epitope within the pre-proinsulin 9-28 peptide has a somewhat unusual binding motif, with tryptophan in the fourth binding pocket of DRB4, perhaps influencing the availability of this complex for T cell selection.