Unraveling the structural basis for the unusually rich association of human leukocyte antigen DQ2.5 with class-II-associated invariant chain peptides

HLA-DM and HLA-DO interplay for the peptide editing of HLA class II in healthy tissues and leukemia

HLA class II antigen presentation is modulated by the activity of the peptide editor HLA-DM and its antagonist HLA-DO, with their interplay controlling the peptide repertoires presented by normal and malignant cells. The role of these molecules in allogeneic hematopoietic cell transplantation (alloHCT) is poorly investigated. Balanced expression of HLA-DM and HLA-DO can influence the presentation of leukemia-associated antigens and peptides targeted by alloreactive T cells, therefore affecting both anti-leukemia immunity and the potential onset of Graft versus Host Disease. We leveraged on a large collection of bulk and single cell RNA sequencing data, available at different repositories, to comprehensively review the level and distribution of HLA-DM and HLA-DO in different cell types and tissues of the human body. The resulting expression atlas will help future investigations aiming to dissect the dual role of HLA class II peptide editing in alloHCT, and their potential impact on its clinical outcome.

Structural insights into human MHC-II association with invariant chain

The loading of processed peptides on to major histocompatibility complex II (MHC-II) molecules for recognition by T cells is vital to cell-mediated adaptive immunity. As part of this process, MHC-II associates with the invariant chain (Ii) during biosynthesis in the endoplasmic reticulum to prevent premature peptide loading and to serve as a scaffold for subsequent proteolytic processing into MHC-II-CLIP. Cryo-electron microscopy structures of full-length Human Leukocyte Antigen-DR (HLA-DR) and HLA-DQ complexes associated with Ii, resolved at 3.0 to 3.1 Å, elucidate the trimeric assembly of the HLA/Ii complex and define atomic-level interactions between HLA, Ii transmembrane domains, loop domains, and class II-associated invariant chain peptides (CLIP). Together with previous structures of MHC-II peptide loading intermediates DO and DM, our findings complete the structural path governing class II antigen presentation.

A general platform for targeting MHC-II antigens via a single loop

Class-II major histocompatibility complexes (MHC-IIs) are central to the communications between CD4+ T cells and antigen presenting cells (APCs), but intrinsic structural features associated with MHC-II make it difficult to develop a general targeting system with high affinity and antigen specificity. Here, we introduce a protein platform, Targeted Recognition of Antigen-MHC Complex Reporter for MHC-II (TRACeR-II), to enable the rapid development of peptide-specific MHC-II binders. TRACeR-II has a small helical bundle scaffold and uses an unconventional mechanism to recognize antigens via a single loop. This unique antigen-recognition mechanism renders this platform highly versatile and amenable to direct structural modeling of the interactions with the antigen. We demonstrate that TRACeR-II binders can be rapidly evolved across multiple alleles, while computational protein design can produce specific binding sequences for a SARS-CoV-2 peptide of unknown complex structure. TRACeR-II sheds light on a simple and straightforward approach to address the MHC peptide targeting challenge, without relying on combinatorial selection on complementarity determining region (CDR) loops. It presents a promising basis for further exploration in immune response modulation as well as a broad range of theragnostic applications.

Analysis of CD74 Occurrence in Oncogenic Fusion Proteins

CD74 is a type II cell surface receptor found to be highly expressed in several hematological and solid cancers, due to its ability to activate pathways associated with tumor cell survival and proliferation. Over the past 16 years, CD74 has emerged as a commonly detected fusion partner in multiple oncogenic fusion proteins. Studies have found CD74 fusion proteins in a range of cancers, including lung adenocarcinoma, inflammatory breast cancer, and pediatric acute lymphoblastic leukemia. To date, there are five known CD74 fusion proteins, CD74-ROS1, CD74-NTRK1, CD74-NRG1, CD74-NRG2α, and CD74-PDGFRB, with a total of 16 different variants, each with unique genetic signatures. Importantly, the occurrence of CD74 in the formation of fusion proteins has not been well explored despite the fact that ROS1 and NRG1 families utilize CD74 as the primary partner for the formation of oncogenic fusions. Fusion proteins known to be oncogenic drivers, including those of CD74, are typically detected and targeted after standard chemotherapeutic plans fail and the disease relapses. The analysis reported herein provides insights into the early intervention of CD74 fusions and highlights the need for improved routine assessment methods so that targeted therapies can be applied while they are most effective.

HLA allele-calling using multi-ancestry whole-exome sequencing from the UK Biobank identifies 129 novel associations in 11 autoimmune diseases

The human leukocyte antigen (HLA) region on chromosome 6 is strongly associated with many immune-mediated and infection-related diseases. Due to its highly polymorphic nature and complex linkage disequilibrium patterns, traditional genetic association studies of single nucleotide polymorphisms do not perform well in this region. Instead, the field has adopted the assessment of the association of HLA alleles (i.e., entire HLA gene haplotypes) with disease. Often based on genotyping arrays, these association studies impute HLA alleles, decreasing accuracy and thus statistical power for rare alleles and in non-European ancestries. Here, we use whole-exome sequencing (WES) from 454,824 UK Biobank (UKB) participants to directly call HLA alleles using the HLA-HD algorithm. We show this method is more accurate than imputing HLA alleles and harness the improved statistical power to identify 360 associations for 11 auto-immune phenotypes (at least 129 likely novel), leading to better insights into the specific coding polymorphisms that underlie these diseases. We show that HLA alleles with synonymous variants, often overlooked in HLA studies, can significantly influence these phenotypes. Lastly, we show that HLA sequencing may improve polygenic risk scores accuracy across ancestries. These findings allow better characterization of the role of the HLA region in human disease.

HLA-II immunopeptidome profiling and deep learning reveal features of antigenicity to inform antigen discovery

CD4+ T cell responses are exquisitely antigen specific and directed toward peptide epitopes displayed by human leukocyte antigen class II (HLA-II) on antigen-presenting cells. Underrepresentation of diverse alleles in ligand databases and an incomplete understanding of factors affecting antigen presentation in vivo have limited progress in defining principles of peptide immunogenicity. Here, we employed monoallelic immunopeptidomics to identify 358,024 HLA-II binders, with a particular focus on HLA-DQ and HLA-DP. We uncovered peptide-binding patterns across a spectrum of binding affinities and enrichment of structural antigen features. These aspects underpinned the development of context-aware predictor of T cell antigens (CAPTAn), a deep learning model that predicts peptide antigens based on their affinity to HLA-II and full sequence of their source proteins. CAPTAn was instrumental in discovering prevalent T cell epitopes from bacteria in the human microbiome and a pan-variant epitope from SARS-CoV-2. Together CAPTAn and associated datasets present a resource for antigen discovery and the unraveling genetic associations of HLA alleles with immunopathologies.

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.

Multiple autoimmunity and epitope spreading in monozygotic twins

We report clinical, serologic, and immunogenetic studies of a set of monozygotic male twin patients who develop autoimmune thyroiditis and vitiligo associated with the HLA-DRB1*04-DQB1*03:02 and HLA-DRB1*03-DQB1*0201 haplotypes. The patients had detectable anti-thyroid and anti-melanocyte autoantibodies. A critical review is presented regarding the role of MHC II molecules linked to clinical manifestations of various autoimmune diseases displayed in a single patient, as is the case in the twin patients reported here.

•

Multiple autoimmunity is a clinicopathological issue that is not well understood.

•

Monozygotic twins with thyroiditis, vitiligo, HLA-DRB1*04-DQB1*03:02 and HLA-DRB1*03-DQBI*0201 haplotypes.

•

Multiple autoantibodies related with intermolecular epitope spreading. Epitope handling by MHC proteins probably related with multiple autoimmunity.

Differential Frequencies of HLA-DRB1, DQA1, and DQB1 Alleles and Haplotypes Are Observed in the Arbovirus-Related Neurological Syndromes

BACKGROUND

We took advantage of the 2015-2016 Brazilian arbovirus outbreak (Zika [ZIKV]/dengue/chikungunya viruses) associated with neurological complications to type HLA-DRB1/DQA1/DQB1 variants in patients exhibiting neurological complications and in bone marrow donors from the same endemic geographical region.

METHODS

DRB1/DQA1/DQB1 loci were typed using sequence-specific oligonucleotides. In silico studies were performed using X-ray resolved dimer constructions.

RESULTS

The DQA1*01, DQA1*05, DQB1*02, or DQB1*06 genotypes/haplotypes and DQA1/DQB1 haplotypes that encode the putative DQA1/DQB1 dimers were overrepresented in the whole group of patients and in patients exhibiting peripheral neurological spectrum disorders (PSD) or encephalitis spectrum disorders (ESD). The DRB1*04, DRB1*13, and DQA1*03 allele groups protected against arbovirus neurological manifestation, being underrepresented in whole group of patients and ESD and PSD groups. Genetic and in silico studies revealed that DQA1/DQB1 dimers (1) were primarily associated with susceptibility to arbovirus infections; (2) can bind to a broad range of ZIKV peptides (235 of 1878 peptides, primarily prM and NS2A); and (3) exhibited hydrophilic and highly positively charged grooves when compared to the DRA1/DRB1 cleft. The protective dimer (DRA1/DRB1*04) bound a limited number of ZIKV peptides (40 of 1878 peptides, primarily prM).

CONCLUSION

Protective haplotypes may recognize arbovirus peptides more specifically than susceptible haplotypes.

The Crystal Structure of the MHC Class I (MHC-I) Molecule in the Green Anole Lizard Demonstrates the Unique MHC-I System in Reptiles

The reptile MHC class I (MCH-I) and MHC class II proteins are the key molecules in the immune system; however, their structure has not been investigated. The crystal structure of green anole lizard peptide-MHC-I-β2m (pMHC-I or pAnca-UA*0101) was determined in the current study. Subsequently, the features of pAnca-UA*0101 were analyzed and compared with the characteristics of pMHC-I of four classes of vertebrates. The amino acid sequence identities between Anca-UA*0101 and MHC-I from other species are <50%; however, the differences between the species were reflected in the topological structure. Significant characteristics of pAnca-UA*0101 include a specific flip of ∼88° and an upward shift adjacent to the C terminus of the α1- and α2-helical regions, respectively. Additionally, the lizard MHC-I molecule has an insertion of 2 aa (VE) at positions 55 and 56. The pushing force from 55-56VE triggers the flip of the α1 helix. Mutagenesis experiments confirmed that the 55-56VE insertion in the α1 helix enhances the stability of pAnca-UA*0101. The peptide presentation profile and motif of pAnca-UA*0101 were confirmed. Based on these results, the proteins of three reptile lizard viruses were used for the screening and confirmation of the candidate epitopes. These data enhance our understanding of the systematic differences between five classes of vertebrates at the gene and protein levels, the formation of the pMHC-I complex, and the evolution of the MHC-I system.

Comparison of HLA ligand elution data and binding predictions reveals varying prediction performance for the multiple motifs recognized by HLA-DQ2.5

Binding prediction tools are commonly used to identify peptides presented on MHC class II molecules. Recently, a wealth of data in the form of naturally eluted ligands has become available and discrepancies between ligand elution data and binding predictions have been reported. Quantitative metrics for such comparisons are currently lacking. In this study, we assessed how efficiently MHC class II binding predictions can identify naturally eluted peptides, and investigated instances with discrepancies between the two methods in detail. We found that, in general, MHC class II eluted ligands are predicted to bind to their reported restriction element with high affinity. But, for several studies reporting an increased number of ligands that were not predicted to bind, we found that the reported MHC restriction was ambiguous. Additional analyses determined that most of the ligands predicted to not bind, are predicted to bind other co‐expressed MHC class II molecules. For selected alleles, we addressed discrepancies between elution data and binding predictions by experimental measurements and found that predicted and measured affinities correlate well. For DQA1*05:01/DQB1*02:01 (DQ2.5) however, binding predictions did miss several peptides that were determined experimentally to be binders. For these peptides and several known DQ2.5 binders, we determined key residues for conferring DQ2.5 binding capacity, which revealed that DQ2.5 utilizes two different binding motifs, of which only one is predicted effectively. These findings have important implications for the interpretation of ligand elution data and for the improvement of MHC class II binding predictions.

HLA-DO Modulates the Diversity of the MHC-II Self-peptidome

Presentation of antigenic peptides on MHC-II molecules is essential for tolerance to self and for initiation of immune responses against foreign antigens. DO (HLA-DO in humans, H2-O in mice) is a nonclassical MHC-II protein that has been implicated in control of autoimmunity and regulation of neutralizing antibody responses to viruses. These effects likely are related to a role of DO in selecting MHC-II epitopes, but previous studies examining the effect of DO on presentation of selected CD4 T cell epitopes have been contradictory. To understand how DO modulates MHC-II antigen presentation, we characterized the full spectrum of peptides presented by MHC-II molecules expressed by DO-sufficient and DO-deficient antigen-presenting cells in vivo and in vitro using quantitative mass spectrometry approaches. We found that DO controlled the diversity of the presented peptide repertoire, with a subset of peptides presented only when DO was expressed. Antigen-presenting cells express another nonclassical MHC-II protein, DM, which acts as a peptide editor by preferentially catalyzing the exchange of less stable MHC-II peptide complexes, and which is inhibited when bound to DO. Peptides presented uniquely in the presence of DO were sensitive to DM-mediated exchange, suggesting that decreased DM editing was responsible for the increased diversity. DO-deficient mice mounted CD4 T cell responses against wild-type antigen-presenting cells, but not vice versa, indicating that DO-dependent alterations in the MHC-II peptidome could be recognized by circulating T cells. These data suggest that cell-specific and regulated expression of HLA-DO serves to fine-tune MHC-II peptidomes, in order to enhance self-tolerance to a wide spectrum of epitopes while allowing focused presentation of immunodominant epitopes during an immune response.

What to do with HLA-DO/H-2O two decades later?

The main objective of antigen processing is to orchestrate the selection of immunodominant epitopes for recognition by CD4 T cells. To achieve this, MHC class II molecules have evolved with a flexible peptide-binding groove in need of a bound peptide. Newly synthesized MHC-II molecules bind a class II invariant chain (Ii) upon synthesis and are shuttled to a specialized compartment, where they encounter exogenous antigens. Ii serves multiple functions, one of which is to maintain the shape of the MHC-II groove so that it can readily bind exogenous antigens upon dissociation of the Ii peptide in MHC- II compartment. MIIC contains processing enzymes, one or both accessory molecules, HLA-DM/H2-M (DM) and HLA-DO/H2-O (DO), and optimal denaturing conditions. In a process known as "editing," DM facilitates the dissociation of the invariant chain peptide, CLIP, for exchange with exogenous antigens. Despite the availability of mechanistic insights into DM functions, understanding how DO contributes to epitope selection has proven to be more challenging. The current dogma assumes that DO inhibits DM, whereas an opposing model suggests that DO fine-tunes the epitope selection process. Understanding which of these, or potentially other models of DO function is important, as DO variants have been linked to autoimmunity, cancer, and the generation of broadly neutralizing antibodies to viruses. This review therefore attempts to evaluate experimental evidence in support of these hypotheses, with an emphasis on the less discussed model, and to explore intriguing questions about the importance of DO in biology.

Class II MHC antigen processing in immune tolerance and inflammation

Presentation of peptide antigens by MHC-II proteins is prerequisite to effective CD4 T cell tolerance to self and to recognition of foreign antigens. Antigen uptake and processing pathways as well as expression of the peptide exchange factors HLA-DM and HLA-DO differ among the various professional and non-professional antigen-presenting cells and are modulated by cell developmental state and activation. Recent studies have highlighted the importance of these cell-specific factors in controlling the source and breadth of peptides presented by MHC-II under different conditions. During inflammation, increased presentation of selected self-peptides has implications for maintenance of peripheral tolerance and autoimmunity.

Dynamics Determine Signaling in a Multicomponent System Associated with Rheumatoid Arthritis

Strategies that target multiple components are usually required for treatment of diseases originating from complex biological systems. The multicomponent system consisting of the DR4 major histocompatibility complex type II molecule, the glycopeptide CII259-273 from type II collagen, and a T-cell receptor is associated with development of rheumatoid arthritis (RA). We introduced non-native amino acids and amide bond isosteres into CII259-273 and investigated the effect on binding to DR4 and the subsequent T-cell response. Molecular dynamics simulations revealed that complexes between DR4 and derivatives of CII259-273 were highly dynamic. Signaling in the overall multicomponent system was found to depend on formation of an appropriate number of dynamic intramolecular hydrogen bonds between DR4 and CII259-273, together with the positioning of the galactose moiety of CII259-273 in the DR4 binding groove. Interestingly, the system tolerated modifications at several positions in CII259-273, indicating opportunities to use analogues to increase our understanding of how rheumatoid arthritis develops and for evaluation as vaccines to treat RA.

Genetic variants at the 16p13 locus confer risk for eosinophilic esophagitis

Eosinophilic esophagitis (EoE) is a chronic inflammatory disease of the esophagus triggered by immune hypersensitivity to food. Herein, we tested whether genetic risk factors for known, non-allergic, immune-mediated diseases, particularly those involving autoimmunity, were associated with EoE risk. We used the high-density Immunochip platform, encoding 200,000 genetic variants for major auto-immune disease. Accordingly, 1214 subjects with EoE of European ancestry and 3734 population controls were genotyped and assessed using data directly generated or imputed from the previously published GWAS. We found lack of association of EoE with the genetic variants in the major histocompatibility complex (MHC) class I, II, and III genes and nearly all other loci using a highly powered study design with dense genotyping throughout the locus. Importantly, we identified an EoE risk locus at 16p13 with genome-wide significance (P_combined=2.05 × 10⁻⁹, odds ratio = 0.76−0.81). This region is known to encode for the genes CLEC16A, DEXI, and CIITI, which are expressed in immune cells and esophageal epithelial cells. Suggestive EoE risk were also seen 5q23 (intergenic) and 7p15 (JAZF1). Overall, we have identified an additional EoE risk locus at 16p13 and highlight a shared and unique genetic etiology of EoE with a spectrum of immune-associated diseases.

Quantification of HLA-DM-Dependent Major Histocompatibility Complex of Class II Immunopeptidomes by the Peptide Landscape Antigenic Epitope Alignment Utility

The major histocompatibility complex of class II (MHCII) immunopeptidome represents the repertoire of antigenic peptides with the potential to activate CD4⁺ T cells. An understanding of how the relative abundance of specific antigenic epitopes affects the outcome of T cell responses is an important aspect of adaptive immunity and offers a venue to more rationally tailor T cell activation in the context of disease. Recent advances in mass spectrometric instrumentation, computational power, labeling strategies, and software analysis have enabled an increasing number of stratified studies on HLA ligandomes, in the context of both basic and translational research. A key challenge in the case of MHCII immunopeptidomes, often determined for different samples at distinct conditions, is to derive quantitative information on consensus epitopes from antigenic peptides of variable lengths. Here, we present the design and benchmarking of a new algorithm [peptide landscape antigenic epitope alignment utility (PLAtEAU)] allowing the identification and label-free quantification (LFQ) of shared consensus epitopes arising from series of nested peptides. The algorithm simplifies the complexity of the dataset while allowing the identification of nested peptides within relatively short segments of protein sequences. Moreover, we apply this algorithm to the comparison of the ligandomes of cell lines with two different expression levels of the peptide-exchange catalyst HLA-DM. Direct comparison of LFQ intensities determined at the peptide level is inconclusive, as most of the peptides are not significantly enriched due to poor sampling. Applying the PLAtEAU algorithm for grouping of the peptides into consensus epitopes shows that more than half of the total number of epitopes is preferentially and significantly enriched for each condition. This simplification and deconvolution of the complex and ambiguous peptide-level dataset highlights the value of the PLAtEAU algorithm in facilitating robust and accessible quantitative analysis of immunopeptidomes across cellular contexts. In silico analysis of the peptides enriched for each HLA-DM expression conditions suggests a higher affinity of the pool of peptides isolated from the high DM expression samples. Interestingly, our analysis reveals that while for certain autoimmune-relevant epitopes their presentation increases upon DM expression others are clearly edited out from the peptidome.

A TCRα framework-centered codon shapes a biased T cell repertoire through direct MHC and CDR3β interactions

Selection of biased T cell receptor (TCR) repertoires across individuals is seen in both infectious diseases and autoimmunity, but the underlying molecular basis leading to these shared repertoires remains unclear. Celiac disease (CD) occurs primarily in HLA-DQ2.5+ individuals and is characterized by a CD4+ T cell response against gluten epitopes dominated by DQ2.5-glia-α1a and DQ2.5-glia-α2. The DQ2.5-glia-α2 response recruits a highly biased TCR repertoire composed of TRAV26-1 paired with TRBV7-2 harboring a semipublic CDR3β loop. We aimed to unravel the molecular basis for this signature. By variable gene segment exchange, directed mutagenesis, and cellular T cell activation studies, we found that TRBV7-3 can substitute for TRBV7-2, as both can contain the canonical CDR3β loop. Furthermore, we identified a pivotal germline-encoded MHC recognition motif centered on framework residue Y40 in TRAV26-1 engaging both DQB1*02 and the canonical CDR3β. This allowed prediction of expanded DQ2.5-glia-α2-reactive TCR repertoires, which were confirmed by single-cell sorting and TCR sequencing from CD patient samples. Our data refine our understanding of how HLA-dependent biased TCR repertoires are selected in the periphery due to germline-encoded residues.