HIV peptidome-wide association study reveals patient-specific epitope repertoires associated with HIV control

Immune-driven clonal cell selection at the intersection among cancer, infections, autoimmunity and senescence

Immune surveillance mechanisms play a crucial role in maintaining lifelong immune homeostasis in response to pathologic stimuli and aberrant cell states. However, their persistence, especially in the context of chronic antigenic exposure, can create a fertile ground for immune evasion. These escaping cell phenotypes, harboring a variety of genomic and transcriptomic aberrances, chiefly in human leukocyte antigen (HLA) and antigen presentation machinery genes, may survive and proliferate, featuring a scenario of clonal cell expansion with immune failure characteristics. While well characterized in solid and, to some extent, hematological malignancies, little is known about their occurrence and significance in other disease contexts. Historical literature highlights the role for escaping HLA-mediated recognition as a strategy adopted by virus to evade from the immune system, hinting at the potential for immune aberrant cell expansion in the context of chronic infections. Additionally, unmasked in idiopathic aplastic anemia as a mechanism able to rescue failing hematopoiesis, HLA clonal escape may operate in autoimmune disorders, particularly in tissues targeted by aberrant immune responses. Furthermore, senescent cell status emerging as immunogenic phenotypes stimulating T cell responses, may act as a bottleneck for the selection of such immune escaping clones, blurring the boundaries between neoplastic transformation, aging and inflammation. Here we provide a fresh overview and perspective on this immune-driven clonal cell expansion, linking pathophysiological features of neoplastic, autoimmune, infectious and senescence processes exposed to immune surveillance.

Associations of functional human leucocyte antigen class I groups with HIV viral load in a heterogeneous cohort

OBJECTIVE

Human leucocyte antigen (HLA) class I alleles are the main host genetic factors involved in controlling HIV-1 viral load (VL). Nevertheless, HLA diversity has proven a significant challenge in association studies. We assessed how accounting for binding affinities of HLA class I alleles to HIV-1 peptides facilitate association testing of HLA with HIV-1 VL in a heterogeneous cohort.

DESIGN

Cohort from the Strategic Timing of AntiRetroviral Treatment (START) study.

METHODS

We imputed HLA class I alleles from host genetic data (2546 HIV+ participants) and sampled immunopeptidomes from 2079 host-paired viral genomes (targeted amplicon sequencing). We predicted HLA class I binding affinities to HIV-1 and unspecific peptides, grouping alleles into functional clusters through consensus clustering. These functional HLA class I clusters were used to test associations with HIV VL.

RESULTS

We identified four clades totaling 30 HLA alleles accounting for 11.4% variability in VL. We highlight HLA-B∗57:01 and B∗57:03 as functionally similar but yet overrepresented in distinct ethnic groups, showing when combined a protective association with HIV+ VL (log, β -0.25; adj. P-value < 0.05). We further demonstrate only a slight power reduction when using unspecific immunopeptidomes, facilitating the use of the inferred functional HLA groups in other studies.

CONCLUSION

The outlined computational approach provides a robust and efficient way to incorporate HLA function and peptide diversity, aiding clinical association studies in heterogeneous cohorts. To facilitate access to the proposed methods and results we provide an interactive application for exploring data.

Immunogenetic Predisposition to SARS-CoV-2 Infection

Herein, we included 527 individuals from two Hospitals, Chemnitz and University-Hospital Leipzig. In total, 199 were negative for PCR and 328 were positive upon first admission. We used next generation sequencing for HLA-A, B, C, DRB1, DRB345, DQA1, DQB1, DPA1, and DPB1, and in some cases, HLA-E, F, G, and H. Furthermore, we molecularly defined 22 blood group systems comprising 26 genes and 5 platelet antigen genes. We observed a significant enrichment of homozygosity for DQA/DQB in the positive group. Within the negative subjects, HLA-B*57:01, HLA-B*55:01, DRB1*13:01, and DRB1*01:01 were enriched, and in the positive group, homozygosity for DQA/DQB, DRB1*09:01, and DRB1*15:01 was observed. DQA1*01:01, DQA1*02:01, and DQA1*01:03 were enriched in the negative group. HLA-DQB1*06:02 was enriched in the positive group, and HLA-DQB1*05:01 and HLA-DQB1*06:03 were enriched in the negative group. For the blood group systems MNS, RH, LE, FY, JK, YT, DO, and KN, enrichment was seen in both groups, depending on the antigen under observation. Homozygosity for D-positive RHD alleles, as well as the phenotypes M-N+ of the MNS blood group system and Yk(a-) of the KN system, were enriched in the positive group. All of these significances disappeared upon correction. Subjects who carried homozygous HPA-1a were more frequent in the negative group, contrasting with the finding that HPA-1ab was enriched in the positive group.

Detailed stratified GWAS analysis for severe COVID-19 in four European populations

Given the highly variable clinical phenotype of Coronavirus disease 2019 (COVID-19), a deeper analysis of the host genetic contribution to severe COVID-19 is important to improve our understanding of underlying disease mechanisms. Here, we describe an extended genome-wide association meta-analysis of a well-characterized cohort of 3255 COVID-19 patients with respiratory failure and 12 488 population controls from Italy, Spain, Norway and Germany/Austria, including stratified analyses based on age, sex and disease severity, as well as targeted analyses of chromosome Y haplotypes, the human leukocyte antigen region and the SARS-CoV-2 peptidome. By inversion imputation, we traced a reported association at 17q21.31 to a ~0.9-Mb inversion polymorphism that creates two highly differentiated haplotypes and characterized the potential effects of the inversion in detail. Our data, together with the 5th release of summary statistics from the COVID-19 Host Genetics Initiative including non-Caucasian individuals, also identified a new locus at 19q13.33, including NAPSA, a gene which is expressed primarily in alveolar cells responsible for gas exchange in the lung.

Individual HLA heterogeneity and its implications for cellular immune evasion in cancer and beyond

Structural and functional variability of human leukocyte antigen (HLA) is the foundation for competent adaptive immune responses against pathogen and tumor antigens as it assures the breadth of the presented immune-peptidome, theoretically sustaining an efficient and diverse T cell response. This variability is presumably the result of the continuous selection by pathogens, which over the course of evolution shaped the adaptive immune system favoring the assortment of a hyper-polymorphic HLA system able to elaborate efficient immune responses. Any genetic alteration affecting this diversity may lead to pathological processes, perturbing antigen presentation capabilities, T-cell reactivity and, to some extent, natural killer cell functionality. A highly variable germline HLA genotype can convey immunogenetic protection against infections, be associated with tumor surveillance or influence response to anti-neoplastic treatments. In contrast, somatic aberrations of HLA loci, rearranging the original germline configuration, theoretically decreasing its variability, can facilitate mechanisms of immune escape that promote tumor growth and immune resistance. The purpose of the present review is to provide a unified and up-to-date overview of the pathophysiological consequences related to the perturbations of the genomic heterogeneity of HLA complexes and their impact on human diseases, with a special focus on cancer.

Challenging ancient DNA results about putative HLA protection or susceptibility to Yersina pestis

In a recent article, Immel et al. (Immel A, Key FM, Szolek A, Barquera R, Robinson MK, Harrison GF, Palmer WH, Spyrou MA, Susat J, Krause-Kyora B, et al. 2021. Analysis of genomic DNA from medieval plague victims suggests long-term effect of Yersinia pestis on human immunity genes. Mol Biol Evol. 38:4059–4076) extracted DNA from 36 individuals dead from plague in Ellwangen, Southern Germany, during the 16th century. By comparing their human leukocyte antigen (HLA) genotypes with those of 50 present-day Ellwangen inhabitants, the authors reported a significant decrease of HLA-B*51:01 and HLA-C*06:02 and a significant increase of HLA-DRB1*13:01/13:02 frequencies from ancient to modern populations. After comparing these frequencies with a larger sample of 8,862 modern Germans and performing simulations of natural selection, they concluded that these changes had been driven by natural selection. In an attempt to provide more evidence on such stimulating results, we explored the HLA frequency patterns over all of Europe, we predicted binding affinities of HLA-B/C/DRB1 alleles to 106,515 Yersinia pestis-derived peptides, and we performed forward simulations of HLA genetic profiles under neutrality. Our analyses do not sustain the conclusions of HLA protection or susceptibility to plague based on ancient DNA.

Evolutionary pressures rendered by animal husbandry practices for avian influenza viruses to adapt to humans

Commercial poultry operations produce and crowd billions of birds every year, which is a source of inexpensive animal protein. Commercial poultry is intensely bred for desirable production traits, and currently presents very low variability at the major histocompatibility complex. This situation dampens the advantages conferred by the MHC’s high genetic variability, and crowding generates immunosuppressive stress. We address the proteins of influenza A viruses directly and indirectly involved in host specificities. We discuss how mutants with increased virulence and/or altered host specificity may arise if few class I alleles are the sole selective pressure on avian viruses circulating in immunocompromised poultry. This hypothesis is testable with peptidomics of MHC ligands. Breeding strategies for commercial poultry can easily and inexpensively include high variability of MHC as a trait of interest, to help save billions of dollars as a disease burden caused by influenza and decrease the risk of selecting highly virulent strains.

The similarity of class II HLA genotypes defines patterns of autoreactivity in idiopathic bone marrow failure disorders

Idiopathic aplastic anemia (IAA) is a rare autoimmune bone marrow failure (BMF) disorder initiated by a human leukocyte antigen (HLA)-restricted T-cell response to unknown antigens. As in other autoimmune disorders, the predilection for certain HLA profiles seems to represent an etiologic factor; however, the structure-function patterns involved in the self-presentation in this disease remain unclear. Herein, we analyzed the molecular landscape of HLA complexes of a cohort of 300 IAA patients and almost 3000 healthy and disease controls by deeply dissecting their genotypic configurations, functional divergence, self-antigen binding capabilities, and T-cell receptor (TCR) repertoire specificities. Specifically, analysis of the evolutionary divergence of HLA genotypes (HED) showed that IAA patients carried class II HLA molecules whose antigen-binding sites were characterized by a high level of structural homology, only partially explained by specific risk allele profiles. This pattern implies reduced HLA binding capabilities, confirmed by binding analysis of hematopoietic stem cell (HSC)-derived self-peptides. IAA phenotype was associated with the enrichment in a few amino acids at specific positions within the peptide-binding groove of DRB1 molecules, affecting the interface HLA-antigen-TCR β and potentially constituting the basis of T-cell dysfunction and autoreactivity. When analyzing associations with clinical outcomes, low HED was associated with risk of malignant progression and worse survival, underlying reduced tumor surveillance in clearing potential neoantigens derived from mechanisms of clonal hematopoiesis. Our data shed light on the immunogenetic risk associated with IAA etiology and clonal evolution and on general pathophysiological mechanisms potentially involved in other autoimmune disorders.

Chickens as a simple system for scientific discovery: The example of the MHC

Chickens have played many roles in human societies over thousands of years, most recently as an important model species for scientific discovery, particularly for embryology, virology and immunology. In the last few decades, biomedical models like mice have become the most important model organism for understanding the mechanisms of disease, but for the study of outbred populations, they have many limitations. Research on humans directly addresses many questions about disease, but frank experiments into mechanisms are limited by practicality and ethics. For research into all levels of disease simultaneously, chickens combine many of the advantages of humans and of mice, and could provide an independent, integrated and overarching system to validate and/or challenge the dogmas that have arisen from current biomedical research. Moreover, some important systems are simpler in chickens than in typical mammals. An example is the major histocompatibility complex (MHC) that encodes the classical MHC molecules, which play crucial roles in the innate and adaptive immune systems. Compared to the large and complex MHCs of typical mammals, the chicken MHC is compact and simple, with single dominantly-expressed MHC molecules that can determine the response to infectious pathogens. As a result, some fundamental principles have been easier to discover in chickens, with the importance of generalist and specialist MHC alleles being the latest example.

Immunodominant T-cell epitopes from the SARS-CoV-2 spike antigen reveal robust pre-existing T-cell immunity in unexposed individuals

The COVID-19 pandemic has revealed a range of disease phenotypes in infected patients with asymptomatic, mild, or severe clinical outcomes, but the mechanisms that determine such variable outcomes remain unresolved. In this study, we identified immunodominant CD8 T-cell epitopes in the spike antigen using a novel TCR-binding algorithm. The predicted epitopes induced robust T-cell activation in unexposed donors demonstrating pre-existing CD4 and CD8 T-cell immunity to SARS-CoV-2 antigen. The T-cell reactivity to the predicted epitopes was higher than the Spike-S1 and S2 peptide pools in the unexposed donors. A key finding of our study is that pre-existing T-cell immunity to SARS-CoV-2 is contributed by TCRs that recognize common viral antigens such as Influenza and CMV, even though the viral epitopes lack sequence identity to the SARS-CoV-2 epitopes. This finding is in contrast to multiple published studies in which pre-existing T-cell immunity is suggested to arise from shared epitopes between SARS-CoV-2 and other common cold-causing coronaviruses. However, our findings suggest that SARS-CoV-2 reactive T-cells are likely to be present in many individuals because of prior exposure to flu and CMV viruses.

Unique pathogen peptidomes facilitate pathogen-specific selection and specialization of MHC alleles

A key component of pathogen-specific adaptive immunity in vertebrates is the presentation of pathogen-derived antigenic peptides by major histocompatibility complex (MHC) molecules. The excessive polymorphism observed at MHC genes is widely presumed to result from the need to recognize diverse pathogens, a process called pathogen-driven balancing selection. This process assumes that pathogens differ in their peptidomes—the pool of short peptides derived from the pathogen’s proteome—so that different pathogens select for different MHC variants with distinct peptide-binding properties. Here, we tested this assumption in a comprehensive data set of 51.9 Mio peptides, derived from the peptidomes of 36 representative human pathogens. Strikingly, we found that 39.7% of the 630 pairwise comparisons among pathogens yielded not a single shared peptide and only 1.8% of pathogen pairs shared more than 1% of their peptides. Indeed, 98.8% of all peptides were unique to a single pathogen species. Using computational binding prediction to characterize the binding specificities of 321 common human MHC class-I variants, we investigated quantitative differences among MHC variants with regard to binding peptides from distinct pathogens. Our analysis showed signatures of specialization toward specific pathogens especially by MHC variants with narrow peptide-binding repertoires. This supports the hypothesis that such fastidious MHC variants might be maintained in the population because they provide an advantage against particular pathogens. Overall, our results establish a key selection factor for the excessive allelic diversity at MHC genes observed in natural populations and illuminate the evolution of variable peptide-binding repertoires among MHC variants.

Human immunotypes impose selection on viral genotypes through viral epitope specificity

BACKGROUND

Understanding the genetic interplay between human hosts and infectious pathogens is crucial for how we interpret virulence factors. Here, we tested for associations between HIV and host genetics, and interactive genetic effects on viral load (VL) in HIV+ ART-naive clinical trial participants.

METHODS

HIV genomes were sequenced and the encoded amino acid (AA) variants were associated with VL, human single nucleotide polymorphisms (SNPs) and imputed HLA alleles, using generalized linear models with Bonferroni correction.

RESULTS

Human (388,501 SNPs) and HIV (3,010 variants) genetic data was available for 2,122 persons. Four HIV variants were associated with VL (p-values<1.66×10 -5). Twelve HIV variants were associated with a range of 1-512 human SNPs (p-value<4.28×10 -11). We found 46 associations between HLA alleles and HIV variants (p-values<1.29×10 -7). We found HIV variants and immunotypes when analyzed separately, were associated with lower VL, whereas the opposite was true when analyzed in concert. Epitope binding prediction showed HLA alleles to be weaker binders of associated HIV AA variants relative to alternative variants on the same position.

CONCLUSIONS

Our results show the importance of immunotype specificity on viral antigenic determinants, and the identified genetic interplay puts emphasis that viral and human genetics should be studied in the context of each other.

Does Antigen Glycosylation Impact the HIV-Specific T Cell Immunity?

It is largely unknown how post-translational protein modifications, including glycosylation, impacts recognition of self and non-self T cell epitopes presented by HLA molecules. Data in the literature indicate that O- and N-linked glycosylation can survive epitope processing and influence antigen presentation and T cell recognition. In this perspective, we hypothesize that glycosylation of viral proteins and processed epitopes contribute to the T cell response to HIV. Although there is some evidence for T cell responses to glycosylated epitopes (glyco-epitopes) during viral infections in the literature, this aspect has been largely neglected for HIV. To explore the role of glyco-epitope specific T cell responses in HIV infection we conducted in silico and ex vivo immune studies in individuals with chronic HIV infection. We found that in silico viral protein segments with potentially glycosylable epitopes were less frequently targeted by T cells. Ex vivo synthetically added glycosylation moieties generally masked T cell recognition of HIV derived peptides. Nonetheless, in some cases, addition of simple glycosylation moieties produced neo-epitopes that were recognized by T cells from HIV infected individuals. Herein, we discuss the potential importance of these observations and compare limitations of the employed technology with new methodologies that may have the potential to provide a more accurate assessment of glyco-epitope specific T cell immunity. Overall, this perspective is aimed to support future research on T cells recognizing glycosylated epitopes in order to expand our understanding on how glycosylation of viral proteins could alter host T cell immunity against viral infections.

HIV-1 and human genetic variation

Over the past four decades, research on the natural history of HIV infection has described how HIV wreaks havoc on human immunity and causes AIDS. HIV host genomic research, which aims to understand how human genetic variation affects our response to HIV infection, has progressed from early candidate gene studies to recent multi-omic efforts, benefiting from spectacular advances in sequencing technology and data science. In addition to invading cells and co-opting the host machinery for replication, HIV also stably integrates into our own genome. The study of the complex interactions between the human and retroviral genomes has improved our understanding of pathogenic mechanisms and suggested novel preventive and therapeutic approaches against HIV infection.

Minding the gap in HIV host genetics: opportunities and challenges

Genome-wide association studies (GWAS) have been successful in identifying and confirming novel genetic variants that are associated with diverse HIV phenotypes. However, these studies have predominantly focused on European cohorts. HLA molecules have been consistently associated with HIV outcomes, some of which have been found to be population specific, underscoring the need for diversity in GWAS. Recently, there has been a concerted effort to address this gap that leads to health care (disease prevention, diagnosis, treatment) disparities with marginal improvement. As precision medicine becomes more utilized, non-European individuals will be more and more disadvantaged, as the genetic variants identified in genomic research based on European populations may not accurately reflect that of non-European individuals. Leveraging pre-existing, large, multiethnic cohorts, such as the UK Biobank, 23andMe, and the National Institute of Health's All of Us Research Program, can contribute in raising genomic research in non-European populations and ultimately lead to better health outcomes.

HLA Heterozygote Advantage against HIV-1 Is Driven by Quantitative and Qualitative Differences in HLA Allele-Specific Peptide Presentation

Pathogen-mediated balancing selection is regarded as a key driver of host immunogenetic diversity. A hallmark for balancing selection in humans is the heterozygote advantage at genes of the human leukocyte antigen (HLA), resulting in improved HIV-1 control. However, the actual mechanism of the observed heterozygote advantage is still elusive. HLA heterozygotes may present a broader array of antigenic viral peptides to immune cells, possibly resulting in a more efficient cytotoxic T-cell response. Alternatively, heterozygosity may simply increase the chance to carry the most protective HLA alleles, as individual HLA alleles are known to differ substantially in their association with HIV-1 control. Here, we used data from 6,311 HIV-1-infected individuals to explore the relative contribution of quantitative and qualitative aspects of peptide presentation in HLA heterozygote advantage against HIV. Screening the entire HIV-1 proteome, we observed that heterozygous individuals exhibited a broader array of HIV-1 peptides presented by their HLA class I alleles. In addition, viral load was negatively correlated with the breadth of the HIV-1 peptide repertoire bound by an individual's HLA variants, particularly at HLA-B. This suggests that heterozygote advantage at HLA-B is at least in part mediated by quantitative peptide presentation. We also observed higher HIV-1 sequence diversity among HLA-B heterozygous individuals, suggesting stronger evolutionary pressure from HLA heterozygosity. However, HLA heterozygotes were also more likely to carry certain HLA alleles, including the highly protective HLA-B*57:01 variant, indicating that HLA heterozygote advantage ultimately results from a combination of quantitative and qualitative effects in antigen presentation.