Machine learning predictions of MHC-II specificities reveal alternative binding mode of class II epitopes.

CD4+ T cells orchestrate the adaptive immune response against pathogens and cancer by recognizing epitopes presented on class II major histocompatibility complex (MHC-II) molecules. The high polymorphism of MHC-II genes represents an important hurdle toward accurate prediction and identification of CD4+ T cell epitopes. Here we collected and curated a dataset of 627,013 unique MHC-II ligands identified by mass spectrometry. This enabled us to precisely determine the binding motifs of 88 MHC-II alleles across humans, mice, cattle, and chickens. Analysis of these binding specificities combined with X-ray crystallography refined our understanding of the molecular determinants of MHC-II motifs and revealed a widespread reverse-binding mode in HLA-DP ligands. We then developed a machine-learning framework to accurately predict binding specificities and ligands of any MHC-II allele. This tool improves and expands predictions of CD4+ T cell epitopes and enables us to discover viral and bacterial epitopes following the aforementioned reverse-binding mode.

Intratumoral Tcf1+PD-1+CD8+ T Cells with Stem-like Properties Promote Tumor Control in Response to Vaccination and Checkpoint Blockade Immunotherapy.

Checkpoint blockade mediates a proliferative response of tumor-infiltrating CD8+ T lymphocytes (TILs). The origin of this response has remained elusive because chronic activation promotes terminal differentiation or exhaustion of tumor-specific T cells. Here we identified a subset of tumor-reactive TILs bearing hallmarks of exhausted cells and central memory cells, including expression of the checkpoint protein PD-1 and the transcription factor Tcf1. Tcf1+PD-1+ TILs mediated the proliferative response to immunotherapy, generating both Tcf1+PD-1+ and differentiated Tcf1-PD-1+ cells. Ablation of Tcf1+PD-1+ TILs restricted responses to immunotherapy. Tcf1 was not required for the generation of Tcf1+PD-1+ TILs but was essential for the stem-like functions of these cells. Human TCF1+PD-1+ cells were detected among tumor-reactive CD8+ T cells in the blood of melanoma patients and among TILs of primary melanomas. Thus, immune checkpoint blockade relies not on reversal of T cell exhaustion programs, but on the proliferation of a stem-like TIL subset.

Joint host-pathogen genomic analysis identifies hepatitis B virus mutations associated with human NTCP and HLA class I variation.

Evolutionary changes in the hepatitis B virus (HBV) genome could reflect its adaptation to host-induced selective pressure. Leveraging paired human exome and ultra-deep HBV genome-sequencing data from 567 affected individuals with chronic hepatitis B, we comprehensively searched for the signatures of this evolutionary process by conducting "genome-to-genome" association tests between all human genetic variants and viral mutations. We identified significant associations between an East Asian-specific missense variant in the gene encoding the HBV entry receptor NTCP (rs2296651, NTCP S267F) and mutations within the receptor-binding region of HBV preS1. Through in silico modeling and in vitro preS1-NTCP binding assays, we observed that the associated HBV mutations are in proximity to the NTCP variant when bound and together partially increase binding affinity to NTCP S267F. Furthermore, we identified significant associations between HLA-A variation and viral mutations in HLA-A-restricted T cell epitopes. We used in silico binding prediction tools to evaluate the impact of the associated HBV mutations on HLA presentation and observed that mutations that result in weaker binding affinities to their cognate HLA alleles were enriched. Overall, our results suggest the emergence of HBV escape mutations that might alter the interaction between HBV PreS1 and its cellular receptor NTCP during viral entry into hepatocytes and confirm the role of HLA class I restriction in inducing HBV epitope variations.

Contemplating immunopeptidomes to better predict them.

The identification of T-cell epitopes is key for a complete molecular understanding of immune recognition mechanisms in infectious diseases, autoimmunity and cancer. T-cell epitopes further provide targets for personalized vaccines and T-cell therapy, with several therapeutic applications in cancer immunotherapy and elsewhere. T-cell epitopes consist of short peptides displayed on Major Histocompatibility Complex (MHC) molecules. The recent advances in mass spectrometry (MS) based technologies to profile the ensemble of peptides displayed on MHC molecules - the so-called immunopeptidome - had a major impact on our understanding of antigen presentation and MHC ligands. On the one hand, these techniques enabled researchers to directly identify hundreds of thousands of peptides presented on MHC molecules, including some that elicited T-cell recognition. On the other hand, the data collected in these experiments revealed fundamental properties of antigen presentation pathways and significantly improved our ability to predict naturally presented MHC ligands and T-cell epitopes across the wide spectrum of MHC alleles found in human and other organisms. Here we review recent computational developments to analyze experimentally determined immunopeptidomes and harness these data to improve our understanding of antigen presentation and MHC binding specificities, as well as our ability to predict MHC ligands. We further discuss the strengths and limitations of the latest approaches to move beyond predictions of antigen presentation and tackle the challenges of predicting TCR recognition and immunogenicity.

Building and analyzing metacells in single-cell genomics data.

The advent of high-throughput single-cell genomics technologies has fundamentally transformed biological sciences. Currently, millions of cells from complex biological tissues can be phenotypically profiled across multiple modalities. The scaling of computational methods to analyze and visualize such data is a constant challenge, and tools need to be regularly updated, if not redesigned, to cope with ever-growing numbers of cells. Over the last few years, metacells have been introduced to reduce the size and complexity of single-cell genomics data while preserving biologically relevant information and improving interpretability. Here, we review recent studies that capitalize on the concept of metacells-and the many variants in nomenclature that have been used. We further outline how and when metacells should (or should not) be used to analyze single-cell genomics data and what should be considered when analyzing such data at the metacell level. To facilitate the exploration of metacells, we provide a comprehensive tutorial on the construction and analysis of metacells from single-cell RNA-seq data ( https://github.com/GfellerLab/MetacellAnalysisTutorial ) as well as a fully integrated pipeline to rapidly build, visualize and evaluate metacells with different methods ( https://github.com/GfellerLab/MetacellAnalysisToolkit ).

The MHC Motif Atlas: a database of MHC binding specificities and ligands.

The highly polymorphic Major Histocompatibility Complex (MHC) genes are responsible for the binding and cell surface presentation of pathogen or cancer specific T-cell epitopes. This process is fundamental for eliciting T-cell recognition of infected or malignant cells. Epitopes displayed on MHC molecules further provide therapeutic targets for personalized cancer vaccines or adoptive T-cell therapy. To help visualizing, analyzing and comparing the different binding specificities of MHC molecules, we developed the MHC Motif Atlas (http://mhcmotifatlas.org/). This database contains information about thousands of class I and class II MHC molecules, including binding motifs, peptide length distributions, motifs of phosphorylated ligands, multiple specificities or links to X-ray crystallography structures. The database further enables users to download curated datasets of MHC ligands. By combining intuitive visualization of the main binding properties of MHC molecules together with access to more than a million ligands, the MHC Motif Atlas provides a central resource to analyze and interpret the binding specificities of MHC molecules.

Metacells untangle large and complex single-cell transcriptome networks.

BACKGROUND: Single-cell RNA sequencing (scRNA-seq) technologies offer unique opportunities for exploring heterogeneous cell populations. However, in-depth single-cell transcriptomic characterization of complex tissues often requires profiling tens to hundreds of thousands of cells. Such large numbers of cells represent an important hurdle for downstream analyses, interpretation and visualization. RESULTS: We develop a framework called SuperCell to merge highly similar cells into metacells and perform standard scRNA-seq data analyses at the metacell level. Our systematic benchmarking demonstrates that metacells not only preserve but often improve the results of downstream analyses including visualization, clustering, differential expression, cell type annotation, gene correlation, imputation, RNA velocity and data integration. By capitalizing on the redundancy inherent to scRNA-seq data, metacells significantly facilitate and accelerate the construction and interpretation of single-cell atlases, as demonstrated by the integration of 1.46 million cells from COVID-19 patients in less than two hours on a standard desktop. CONCLUSIONS: SuperCell is a framework to build and analyze metacells in a way that efficiently preserves the results of scRNA-seq data analyses while significantly accelerating and facilitating them.

Structural dissimilarity from self drives neoepitope escape from immune tolerance.

T-cell recognition of peptides incorporating nonsynonymous mutations, or neoepitopes, is a cornerstone of tumor immunity and forms the basis of new immunotherapy approaches including personalized cancer vaccines. Yet as they are derived from self-peptides, the means through which immunogenic neoepitopes overcome immune self-tolerance are often unclear. Here we show that a point mutation in a non-major histocompatibility complex anchor position induces structural and dynamic changes in an immunologically active ovarian cancer neoepitope. The changes pre-organize the peptide into a conformation optimal for recognition by a neoepitope-specific T-cell receptor, allowing the receptor to bind the neoepitope with high affinity and deliver potent T-cell signals. Our results emphasize the importance of structural and physical changes relative to self in neoepitope immunogenicity. Considered broadly, these findings can help explain some of the difficulties in identifying immunogenic neoepitopes from sequence alone and provide guidance for developing novel, neoepitope-based personalized therapies.

Mass Spectrometry Based Immunopeptidomics Leads to Robust Predictions of Phosphorylated HLA Class I Ligands.

The presentation of peptides on class I human leukocyte antigen (HLA-I) molecules plays a central role in immune recognition of infected or malignant cells. In cancer, non-self HLA-I ligands can arise from many different alterations, including non-synonymous mutations, gene fusion, cancer-specific alternative mRNA splicing or aberrant post-translational modifications. Identifying HLA-I ligands remains a challenging task that requires either heavy experimental work for in vivo identification or optimized bioinformatics tools for accurate predictions. To date, no HLA-I ligand predictor includes post-translational modifications. To fill this gap, we curated phosphorylated HLA-I ligands from several immunopeptidomics studies (including six newly measured samples) covering 72 HLA-I alleles and retrieved a total of 2,066 unique phosphorylated peptides. We then expanded our motif deconvolution tool to identify precise binding motifs of phosphorylated HLA-I ligands. Our results reveal a clear enrichment of phosphorylated peptides among HLA-C ligands and demonstrate a prevalent role of both HLA-I motifs and kinase motifs on the presentation of phosphorylated peptides. These data further enabled us to develop and validate the first predictor of interactions between HLA-I molecules and phosphorylated peptides.

The C-terminal extension landscape of naturally presented HLA-I ligands.

HLA-I molecules play a central role in antigen presentation. They typically bind 9- to 12-mer peptides, and their canonical binding mode involves anchor residues at the second and last positions of their ligands. To investigate potential noncanonical binding modes, we collected in-depth and accurate HLA peptidomics datasets covering 54 HLA-I alleles and developed algorithms to analyze these data. Our results reveal frequent (442 unique peptides) and statistically significant C-terminal extensions for at least eight alleles, including the common HLA-A03:01, HLA-A31:01, and HLA-A68:01. High resolution crystal structure of HLA-A68:01 with such a ligand uncovers structural changes taking place to accommodate C-terminal extensions and helps unraveling sequence and structural properties predictive of the presence of these extensions. Scanning viral proteomes with the C-terminal extension motifs identifies many putative epitopes and we demonstrate direct recognition by human CD8+ T cells of a 10-mer epitope from cytomegalovirus predicted to follow the C-terminal extension binding mode.

Single-cell transcriptome analysis of fish immune cells provides insight into the evolution of vertebrate immune cell types.

The immune system of vertebrate species consists of many different cell types that have distinct functional roles and are subject to different evolutionary pressures. Here, we first analyzed conservation of genes specific for all major immune cell types in human and mouse. Our results revealed higher gene turnover and faster evolution of trans-membrane proteins in NK cells compared with other immune cell types, and especially T cells, but similar conservation of nuclear and cytoplasmic protein coding genes. To validate these findings in a distant vertebrate species, we used single-cell RNA sequencing of lck:GFP cells in zebrafish and obtained the first transcriptome of specific immune cell types in a nonmammalian species. Unsupervised clustering and single-cell TCR locus reconstruction identified three cell populations, T cells, a novel type of NK-like cells, and a smaller population of myeloid-like cells. Differential expression analysis uncovered new immune-cell-specific genes, including novel immunoglobulin-like receptors, and neofunctionalization of recently duplicated paralogs. Evolutionary analyses confirmed the higher gene turnover of trans-membrane proteins in NK cells compared with T cells in fish species, suggesting that this is a general property of immune cell types across all vertebrates.

The Length Distribution and Multiple Specificity of Naturally Presented HLA-I Ligands.

HLA-I molecules bind short peptides and present them for recognition by CD8+ T cells. The length of HLA-I ligands typically ranges from 8 to 12 aa, but variability is observed across different HLA-I alleles. In this study we collected recent in-depth HLA peptidomics data, including 12 newly generated HLA peptidomes (31,896 unique peptides) from human meningioma samples, to analyze the peptide length distribution and multiple specificity across 84 different HLA-I alleles. We observed a clear clustering of HLA-I alleles with distinct peptide length distributions, which enabled us to study the structural basis of peptide length distributions and predict peptide length distributions from HLA-I sequences. We further identified multiple specificity in several HLA-I molecules and validated these observations with binding assays. Explicitly modeling peptide length distribution and multiple specificity improved predictions of naturally presented HLA-I ligands, as demonstrated in an independent benchmarking based on the new human meningioma samples.

High-throughput and Sensitive Immunopeptidomics Platform Reveals Profound Interferonγ-Mediated Remodeling of the Human Leukocyte Antigen (HLA) Ligandome.

Comprehensive knowledge of the human leukocyte antigen (HLA) class-I and class-II peptides presented to T-cells is crucial for designing innovative therapeutics against cancer and other diseases. However methodologies for their purification for mass-spectrometry analysis have been a major limitation. We designed a novel high-throughput, reproducible and sensitive method for sequential immuno-affinity purification of HLA-I and -II peptides from up to 96 samples in a plate format, suitable for both cell lines and tissues. Our methodology drastically reduces sample-handling and can be completed within five hours. We challenged our methodology by extracting HLA peptides from multiple replicates of tissues (n = 7) and cell lines (n = 21, 108 cells per replicate), which resulted in unprecedented depth, sensitivity and high reproducibility (Pearson correlations up to 0.98 and 0.97 for HLA-I and HLA-II). Because of the method's achieved sensitivity, even single measurements of peptides purified from 107 B-cells resulted in the identification of more than 1700 HLA-I and 2200 HLA-II peptides. We demonstrate the feasibility of performing drug-screening by using ovarian cancer cells treated with interferon gamma (IFNγ). Our analysis revealed an augmented presentation of chymotryptic-like and longer ligands associated with IFNγ induced changes of the antigen processing and presentation machinery. This straightforward method is applicable for basic and clinical applications.

Estimating the Contribution of Proteasomal Spliced Peptides to the HLA-I Ligandome.

Spliced peptides are short protein fragments spliced together in the proteasome by peptide bond formation. True estimation of the contribution of proteasome-spliced peptides (PSPs) to the global human leukocyte antigen (HLA) ligandome is critical. A recent study suggested that PSPs contribute up to 30% of the HLA ligandome. We performed a thorough reanalysis of the reported results using multiple computational tools and various validation steps and concluded that only a fraction of the proposed PSPs passes the quality filters. To better estimate the actual number of PSPs, we present an alternative workflow. We performed de novo sequencing of the HLA-peptide spectra and discarded all de novo sequences found in the UniProt database. We checked whether the remaining de novo sequences could match spliced peptides from human proteins. The spliced sequences were appended to the UniProt fasta file, which was searched by two search tools at a false discovery rate (FDR) of 1%. We find that 2-6% of the HLA ligandome could be explained as spliced protein fragments. The majority of these potential PSPs have good peptide-spectrum match properties and are predicted to bind the respective HLA molecules. However, it remains to be shown how many of these potential PSPs actually originate from proteasomal splicing events.

Polymorphic sites preferentially avoid co-evolving residues in MHC class I proteins.

Major histocompatibility complex class I (MHC-I) molecules are critical to adaptive immune defence mechanisms in vertebrate species and are encoded by highly polymorphic genes. Polymorphic sites are located close to the ligand-binding groove and entail MHC-I alleles with distinct binding specificities. Some efforts have been made to investigate the relationship between polymorphism and protein stability. However, less is known about the relationship between polymorphism and MHC-I co-evolutionary constraints. Using Direct Coupling Analysis (DCA) we found that co-evolution analysis accurately pinpoints structural contacts, although the protein family is restricted to vertebrates and comprises less than five hundred species, and that the co-evolutionary signal is mainly driven by inter-species changes, and not intra-species polymorphism. Moreover, we show that polymorphic sites in human preferentially avoid co-evolving residues, as well as residues involved in protein stability. These results suggest that sites displaying high polymorphism may have been selected during vertebrates' evolution to avoid co-evolutionary constraints and thereby maximize their mutability.

The pathogen-related yeast protein Pry1, a member of the CAP protein superfamily, is a fatty acid-binding protein.

Members of the CAP superfamily (cysteine-rich secretory proteins, antigen 5, and pathogenesis-related 1 proteins), also known as SCP superfamily (sperm-coating proteins), have been implicated in many physiological processes, including immune defenses, venom toxicity, and sperm maturation. Their mode of action, however, remains poorly understood. Three proteins of the CAP superfamily, Pry1, -2, and -3 (pathogen related in yeast), are encoded in the Saccharomyces cerevisiae genome. We have shown previously that Pry1 binds cholesterol in vitro and that Pry function is required for sterol secretion in yeast cells, indicating that members of this superfamily may generally bind sterols or related small hydrophobic compounds. On the other hand, tablysin-15, a CAP protein from the horsefly Tabanus yao, has been shown to bind leukotrienes and free fatty acids in vitro Therefore, here we assessed whether the yeast Pry1 protein binds fatty acids. Computational modeling and site-directed mutagenesis indicated that the mode of fatty acid binding is conserved between tablysin-15 and Pry1. Pry1 bound fatty acids with micromolar affinity in vitro, and its function was essential for fatty acid export in cells lacking the acyl-CoA synthetases Faa1 and Faa4. Fatty acid binding of Pry1 is independent of its capacity to bind sterols, and the two sterol- and fatty acid-binding sites are nonoverlapping. These results indicate that some CAP family members, such as Pry1, can bind different lipids, particularly sterols and fatty acids, at distinct binding sites, suggesting that the CAP domain may serve as a stable, secreted protein domain that can accommodate multiple ligand-binding sites.

Deciphering HLA-I motifs across HLA peptidomes improves neo-antigen predictions and identifies allostery regulating HLA specificity.

The precise identification of Human Leukocyte Antigen class I (HLA-I) binding motifs plays a central role in our ability to understand and predict (neo-)antigen presentation in infectious diseases and cancer. Here, by exploiting co-occurrence of HLA-I alleles across ten newly generated as well as forty public HLA peptidomics datasets comprising more than 115,000 unique peptides, we show that we can rapidly and accurately identify many HLA-I binding motifs and map them to their corresponding alleles without any a priori knowledge of HLA-I binding specificity. Our approach recapitulates and refines known motifs for 43 of the most frequent alleles, uncovers new motifs for 9 alleles that up to now had less than five known ligands and provides a scalable framework to incorporate additional HLA peptidomics studies in the future. The refined motifs improve neo-antigen and cancer testis antigen predictions, indicating that unbiased HLA peptidomics data are ideal for in silico predictions of neo-antigens from tumor exome sequencing data. The new motifs further reveal distant modulation of the binding specificity at P2 for some HLA-I alleles by residues in the HLA-I binding site but outside of the B-pocket and we unravel the underlying mechanisms by protein structure analysis, mutagenesis and in vitro binding assays.

Unsupervised HLA Peptidome Deconvolution Improves Ligand Prediction Accuracy and Predicts Cooperative Effects in Peptide-HLA Interactions.

Ag presentation on HLA molecules plays a central role in infectious diseases and tumor immunology. To date, large-scale identification of (neo-)Ags from DNA sequencing data has mainly relied on predictions. In parallel, mass spectrometry analysis of HLA peptidome is increasingly performed to directly detect peptides presented on HLA molecules. In this study, we use a novel unsupervised approach to assign mass spectrometry-based HLA peptidomics data to their cognate HLA molecules. We show that incorporation of deconvoluted HLA peptidomics data in ligand prediction algorithms can improve their accuracy for HLA alleles with few ligands in existing databases. The results of our computational analysis of large datasets of naturally processed HLA peptides, together with experimental validation and protein structure analysis, further reveal how HLA-binding motifs change with peptide length and predict new cooperative effects between distant residues in HLA-B07:02 ligands.

Identification of innate lymphoid cells in single-cell RNA-Seq data.

Innate lymphoid cells (ILCs) consist of natural killer (NK) cells and non-cytotoxic ILCs that are broadly classified into ILC1, ILC2, and ILC3 subtypes. These cells recently emerged as important early effectors of innate immunity for their roles in tissue homeostasis and inflammation. Over the last few years, ILCs have been extensively studied in mouse and human at the functional and molecular level, including gene expression profiling. However, sorting ILCs with flow cytometry for gene expression analysis is a delicate and time-consuming process. Here we propose and validate a novel framework for studying ILCs at the transcriptomic level using single-cell RNA-Seq data. Our approach combines unsupervised clustering and a new cell type classifier trained on mouse ILC gene expression data. We show that this approach can accurately identify different ILCs, especially ILC2 cells, in human lymphocyte single-cell RNA-Seq data. Our new model relies only on genes conserved across vertebrates, thereby making it in principle applicable in any vertebrate species. Considering the rapid increase in throughput of single-cell RNA-Seq technology, our work provides a computational framework for studying ILC2 cells in single-cell transcriptomic data and may help exploring their conservation in distant vertebrate species.

Protein homology reveals new targets for bioactive small molecules.

MOTIVATION: The functional impact of small molecules is increasingly being assessed in different eukaryotic species through large-scale phenotypic screening initiatives. Identifying the targets of these molecules is crucial to mechanistically understand their function and uncover new therapeutically relevant modes of action. However, despite extensive work carried out in model organisms and human, it is still unclear to what extent one can use information obtained in one species to make predictions in other species. RESULTS: Here, for the first time, we explore and validate at a large scale the use of protein homology relationships to predict the targets of small molecules across different species. Our results show that exploiting target homology can significantly improve the predictions, especially for molecules experimentally tested in other species. Interestingly, when considering separately orthology and paralogy relationships, we observe that mapping small molecule interactions among orthologs improves prediction accuracy, while including paralogs does not improve and even sometimes worsens the prediction accuracy. Overall, our results provide a novel approach to integrate chemical screening results across multiple species and highlight the promises and remaining challenges of using protein homology for small molecule target identification. AVAILABILITY AND IMPLEMENTATION: Homology-based predictions can be tested on our website http://www.swisstargetprediction.ch. CONTACT: david.gfeller@unil.ch or vincent.zoete@isb-sib.ch. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.