Immunopeptidogenomics: Harnessing RNA-Seq to Illuminate the Dark Immunopeptidome.

Human leukocyte antigen (HLA) molecules are cell-surface glycoproteins that present peptide antigens on the cell surface for surveillance by T lymphocytes, which contemporaneously seek signs of disease. Mass spectrometric analysis allows us to identify large numbers of these peptides (the immunopeptidome) following affinity purification of solubilized HLA-peptide complexes. However, in recent years, there has been a growing awareness of the "dark side" of the immunopeptidome: unconventional peptide epitopes, including neoepitopes, which elude detection by conventional search methods because their sequences are not present in reference protein databases (DBs). Here, we establish a bioinformatics workflow to aid identification of peptides generated by noncanonical translation of mRNA or by genome variants. The workflow incorporates both standard transcriptomics software and novel computer programs to produce cell line-specific protein DBs based on three-frame translation of the transcriptome. The final protein DB also includes sequences resulting from variants determined by variant calling on the same RNA-Seq data. We then searched our experimental data against both transcriptome-based and standard DBs using PEAKS Studio (Bioinformatics Solutions, Inc). Finally, further novel software helps to compare the various result sets arising for each sample, pinpoint putative genomic origins for unconventional sequences, and highlight potential neoepitopes. We applied the workflow to study the immunopeptidome of the acute myeloid leukemia cell line THP-1, using RNA-Seq and immunopeptidome data. We confidently identified over 14,000 peptides from three replicates of purified HLA peptides derived from THP-1 cells using the conventional UniProt human proteome. Using the transcriptome-based DB generated using our workflow, we recapitulated >85% of these and also identified 1029 unconventional peptides not explained by UniProt, including 16 sequences caused by nonsynonymous variants. Our workflow, which we term "immunopeptidogenomics," can provide DBs, which include pertinent unconventional sequences and allow neoepitope discovery, without becoming too large to search. Immunopeptidogenomics is a step toward unbiased search approaches that are needed to illuminate the dark side of the immunopeptidome.

DosR proteins of Mycobacterium tuberculosis upregulate effector T cells and down regulate T regulatory cells in TB patients and their healthy contacts.

It is well established that the current problem of tuberculosis (TB) can be combated by overcoming the drawbacks of the currently available BCG vaccine. This would involve incorporation of antigens that can control TB at all stages including the dormant phase which is generally ignored. Hence, DosR regulon proteins, which are expressed in latent infection, could prove to be very good vaccine candidates as they can possibly target the silent but most predominant form of TB infection. In the present study, the immune response to two DosR proteins Rv2627 and Rv2628 has been studied in PBMCs derived from normal individuals, TB patients and healthy contacts of TB patients. It was found that these antigens were capable of stimulating a strong IFN-γ+ T cell response along with accentuation of memory T cells and other protective cytokines such as IL-2 and IL-17. At the same time these proteins decreased the frequencies of immune-suppressor regulatory T cells in in vitro stimulation of PBMC from both patients and their contacts. Considering all these facts together, we suggest Rv2627 and Rv2628 to be one of the extremely promising candidates for incorporation into a post exposure subunit vaccine against TB.

Combining ERAP1 silencing and entinostat therapy to overcome resistance to cancer immunotherapy in neuroblastoma.

BACKGROUND: Checkpoint immunotherapy unleashes tumor control by T cells, but it is undermined in non-immunogenic tumors, e.g. with low MHC class I expression and low neoantigen burden, such as neuroblastoma (NB). Endoplasmic reticulum aminopeptidase 1 (ERAP1) is an enzyme that trims peptides before loading on MHC class I molecules. Inhibition of ERAP1 results in the generation of new antigens able of inducing potent anti-tumor immune responses. Here, we identify a novel non-toxic combinatorial strategy based on genetic inhibition of ERAP1 and administration of the HDAC inhibitor (HDACi) entinostat that increase the immunogenicity of NB, making it responsive to PD-1 therapy. METHODS: CRISPR/Cas9-mediated gene editing was used to knockout (KO) the ERAP1 gene in 9464D NB cells derived from spontaneous tumors of TH-MYCN transgenic mice. The expression of MHC class I and PD-L1 was evaluated by flow cytometry (FC). The immunopeptidome of these cells was studied by mass spectrometry. Cocultures of splenocytes derived from 9464D bearing mice and tumor cells allowed the assessment of the effect of ERAP1 inhibition on the secretion of inflammatory cytokines and activation and migration of immune cells towards ERAP1 KO cells by FC. Tumor cell killing was evaluated by Caspase 3/7 assay and flow cytometry analysis. The effect of ERAP1 inhibition on the immune content of tumors was analyzed by FC, immunohistochemistry and multiple immunofluorescence. RESULTS: We found that inhibition of ERAP1 makes 9464D cells more susceptible to immune cell-mediated killing by increasing both the recall and activation of CD4+ and CD8+ T cells and NK cells. Treatment with entinostat induces the expression of MHC class I and PD-L1 molecules in 9464D both in vitro and in vivo. This results in pronounced changes in the immunopeptidome induced by ERAP1 inhibition, but also restrains the growth of ERAP1 KO tumors in vivo by remodelling the tumor-infiltrating T-cell compartment. Interestingly, the absence of ERAP1 in combination with entinostat and PD-1 blockade overcomes resistance to PD-1 immunotherapy and increases host survival. CONCLUSIONS: These findings demonstrate that ERAP1 inhibition combined with HDACi entinostat treatment and PD-1 blockade remodels the immune landscape of a non-immunogenic tumor such as NB, making it responsive to checkpoint immunotherapy.

Mapping the immunopeptidome of seven SARS-CoV-2 antigens across common HLA haplotypes.

Most COVID-19 vaccines elicit immunity against the SARS-CoV-2 Spike protein. However, Spike protein mutations in emerging strains and immune evasion by the SARS-CoV-2 virus demonstrates the need to develop more broadly targeting vaccines. To facilitate this, we use mass spectrometry to identify immunopeptides derived from seven relatively conserved structural and non-structural SARS-CoV-2 proteins (N, E, Nsp1/4/5/8/9). We use two different B-lymphoblastoid cell lines to map Human Leukocyte Antigen (HLA) class I and class II immunopeptidomes covering some of the prevalent HLA types across the global human population. We employ DNA plasmid transfection and direct antigen delivery approaches to sample different antigens and find 248 unique HLA class I and HLA class II bound peptides with 71 derived from N, 12 from E, 28 from Nsp1, 19 from Nsp4, 73 from Nsp8 and 45 peptides derived from Nsp9. Over half of the viral peptides are unpublished. T cell reactivity tested against 56 of the detected peptides shows CD8+ and CD4+ T cell responses against several peptides from the N, E, and Nsp9 proteins. Results from this study will aid the development of next-generation COVID vaccines targeting epitopes from across a number of SARS-CoV-2 proteins.

DNAJB1-PRKACA fusion neoantigens elicit rare endogenous T cell responses that potentiate cell therapy for fibrolamellar carcinoma.

Fibrolamellar carcinoma (FLC) is a liver tumor with a high mortality burden and few treatment options. A promising therapeutic vulnerability in FLC is its driver mutation, a conserved DNAJB1-PRKACA gene fusion that could be an ideal target neoantigen for immunotherapy. In this study, we aim to define endogenous CD8 T cell responses to this fusion in FLC patients and evaluate fusion-specific T cell receptors (TCRs) for use in cellular immunotherapies. We observe that fusion-specific CD8 T cells are rare and that FLC patient TCR repertoires lack large clusters of related TCR sequences characteristic of potent antigen-specific responses, potentially explaining why endogenous immune responses are insufficient to clear FLC tumors. Nevertheless, we define two functional fusion-specific TCRs, one of which has strong anti-tumor activity in vivo. Together, our results provide insights into the fragmented nature of neoantigen-specific repertoires in humans and indicate routes for clinical development of successful immunotherapies for FLC.

Endogenous H3.3K27M derived peptide restricted to HLA-A∗02:01 is insufficient for immune-targeting in diffuse midline glioma.

Diffuse midline glioma (DMG) is a childhood brain tumor with an extremely poor prognosis. Chimeric antigen receptor (CAR) T cell therapy has recently demonstrated some success in DMG, but there may a need to target multiple tumor-specific targets to avoid antigen escape. We developed a second-generation CAR targeting an HLA-A∗02:01 restricted histone 3K27M epitope in DMG, the target of previous peptide vaccination and T cell receptor-mimics. These CAR T cells demonstrated specific, titratable, binding to cells pulsed with the H3.3K27M peptide. However, we were unable to observe scFv binding, CAR T cell activation, or cytotoxic function against H3.3K27M+ patient-derived models. Despite using sensitive immunopeptidomics, we could not detect the H3.3K27M26-35-HLA-A∗02:01 peptide on these patient-derived models. Interestingly, other non-mutated peptides from DMG were detected bound to HLA-A∗02:01 and other class I molecules, including a novel HLA-A3-restricted peptide encompassing the K27M mutation and overlapping with the H3 K27M26-35-HLA-A∗02:01 peptide. These results suggest that targeting the H3 K27M26-35 mutation in context of HLA-A∗02:01 may not be a feasible immunotherapy strategy because of its lack of presentation. These findings should inform future investigations and clinical trials in DMG.

Immunolyser: A web-based computational pipeline for analysing and mining immunopeptidomic data.

Immunopeptidomics has made tremendous contributions to our understanding of antigen processing and presentation, by identifying and quantifying antigenic peptides presented on the cell surface by Major Histocompatibility Complex (MHC) molecules. Large and complex immunopeptidomics datasets can now be routinely generated using Liquid Chromatography-Mass Spectrometry techniques. The analysis of this data - often consisting of multiple replicates/conditions - rarely follows a standard data processing pipeline, hindering the reproducibility and depth of analysis of immunopeptidomic data. Here, we present Immunolyser, an automated pipeline designed to facilitate computational analysis of immunopeptidomic data with a minimal initial setup. Immunolyser brings together routine analyses, including peptide length distribution, peptide motif analysis, sequence clustering, peptide-MHC binding affinity prediction, and source protein analysis. Immunolyser provides a user-friendly and interactive interface via its webserver and is freely available for academic purposes at https://immunolyser.erc.monash.edu/. The open-access source code can be downloaded at our GitHub repository: https://github.com/prmunday/Immunolyser. We anticipate that Immunolyser will serve as a prominent computational pipeline to facilitate effortless and reproducible analysis of immunopeptidomic data.

A combined immunopeptidomics, proteomics, and cell surface proteomics approach to identify immunotherapy targets for diffuse intrinsic pontine glioma.

Introduction: Diffuse intrinsic pontine glioma (DIPG), recently reclassified as a subtype of diffuse midline glioma, is a highly aggressive brainstem tumor affecting children and young adults, with no cure and a median survival of only 9 months. Conventional treatments are ineffective, highlighting the need for alternative therapeutic strategies such as cellular immunotherapy. However, identifying unique and tumor-specific cell surface antigens to target with chimeric antigen receptor (CAR) or T-cell receptor (TCR) therapies is challenging. Methods: In this study, a multi-omics approach was used to interrogate patient-derived DIPG cell lines and to identify potential targets for immunotherapy. Results: Through immunopeptidomics, a range of targetable peptide antigens from cancer testis and tumor-associated antigens as well as peptides derived from human endogenous retroviral elements were identified. Proteomics analysis also revealed upregulation of potential drug targets and cell surface proteins such as Cluster of differentiation 27 (CD276) B7 homolog 3 protein (B7H3), Interleukin 13 alpha receptor 2 (IL-13Rα2), Human Epidermal Growth Factor Receptor 3 (HER2), Ephrin Type-A Receptor 2 (EphA2), and Ephrin Type-A Receptor 3 (EphA3). Discussion: The results of this study provide a valuable resource for the scientific community to accelerate immunotherapeutic approaches for DIPG. Identifying potential targets for CAR and TCR therapies could open up new avenues for treating this devastating disease.

Multi-Omic Analysis of Two Common P53 Mutations: Proteins Regulated by Mutated P53 as Potential Targets for Immunotherapy.

The p53 protein is mutated in more than 50% of human cancers. Mutated p53 proteins not only lose their normal function but often acquire novel oncogenic functions, a phenomenon termed mutant p53 gain-of-function. Mutant p53 has been shown to affect the transcription of a range of genes, as well as protein-protein interactions with transcription factors and other effectors; however, no one has intensively investigated and identified these proteins, or their MHC presented epitopes, from the viewpoint of their ability to act as targets for immunotherapeutic interventions. We investigated the molecular changes that occurred after the TP53 null osteosarcoma cells, SaOS-2, were transfected with one of two conformational p53-mutants, either R175H or R273H. We then examined the phenotypic and functional changes using macroscopic observations, proliferation, gene expression and proteomics alongside immunopeptidome profiling of peptide antigen presentation in the context of major histocompatibility complex (MHC) class I molecules. We identified several candidate proteins in both TP53 mutant cell lines with differential expression when compared to the TP53 null vector control, SaOS-V. Quantitative SWATH proteomics combined with immune-peptidome analysis of the class-I eluted peptides identified several epitopes presented on pMHC and in silico analysis shortlisted which antigens were expressed in a range of cancerous but not adjacent healthy tissues. Out of all the candidates, KLC1 and TOP2A showed high levels of expression in every tumor type examined. From these proteins, three A2 and four pan HLA-A epitopes were identified in both R175H and R273H from TOP2A. We have now provided a short list of future immunotherapy targets for the treatment of cancers harboring mutated TP53.

Defective Proinsulin Handling Modulates the MHC I Bound Peptidome and Activates the Inflammasome in ß-Cells.

How immune tolerance is lost to pancreatic ß-cell peptides triggering autoimmune type 1 diabetes is enigmatic. We have shown that loss of the proinsulin chaperone glucose-regulated protein (GRP) 94 from the endoplasmic reticulum (ER) leads to mishandling of proinsulin, ER stress, and activation of the immunoproteasome. We hypothesize that inadequate ER proinsulin folding capacity relative to biosynthetic need may lead to an altered ß-cell major histocompatibility complex (MHC) class-I bound peptidome and inflammasome activation, sensitizing ß-cells to immune attack. We used INS-1E cells with or without GRP94 knockout (KO), or in the presence or absence of GRP94 inhibitor PU-WS13 (GRP94i, 20 µM), or exposed to proinflammatory cytokines interleukin (IL)-1ß or interferon gamma (IFNγ) (15 pg/mL and 10 ng/mL, respectively) for 24 h. RT1.A (rat MHC I) expression was evaluated using flow cytometry. The total RT1.A-bound peptidome analysis was performed on cell lysates fractionated by reverse-phase high-performance liquid chromatography (RP-HPLC), followed by liquid chromatography coupled with tandem mass spectrometry (LC-MS/MS). The nucleotide-binding oligomerization domain, leucine rich repeat and pyrin domain containing protein (NLRP1), nuclear factor of kappa light polypeptide gene enhancer in B-cells inhibitor alpha (IκBα), and (pro) IL-1ß expression and secretion were investigated by Western blotting. GRP94 KO increased RT1.A expression in ß-cells, as did cytokine exposure compared to relevant controls. Immunopeptidome analysis showed increased RT1.A-bound peptide repertoire in GRP94 KO/i cells as well as in the cells exposed to cytokines. The GRP94 KO/cytokine exposure groups showed partial overlap in their peptide repertoire. Notably, proinsulin-derived peptide diversity increased among the total RT1.A peptidome in GRP94 KO/i along with cytokines exposure. NLRP1 expression was upregulated in GRP94 deficient cells along with decreased IκBα content while proIL-1ß cellular levels declined, coupled with increased secretion of mature IL-1ß. Our results suggest that limiting ß-cell proinsulin chaperoning enhances RT1.A expression alters the MHC-I peptidome including proinsulin peptides and activates inflammatory pathways, suggesting that stress associated with impeding proinsulin handling may sensitize ß-cells to immune-attack.

Isolation of HLA Bound Peptides by Immunoaffinity Capture and Identification by Mass Spectrometry.

This article describes the purification of HLA-bound peptides and their subsequent sequencing by mass spectrometry. These methods can be used for both HLA class I and class II molecules and can be adapted to different species depending on the availability of specific antibodies. Peptides can be successfully isolated from a variety of sample types, including in vitro cultured cells and primary tissues. The method involves the affinity capture of HLA-peptide complexes and separation of peptides from HLA heavy chains, followed by tailored interrogation by mass spectrometry to take into account the non-tryptic nature of endogenously derived HLA-bound peptides. © 2021 Wiley Periodicals LLC. Basic Protocol 1: Preparation of immunoaffinity column Alternate Protocol 1: Preparation of microscale immunoaffinity column Basic Protocol 2: Generation of cell lysate and HLA immunoaffinity purification Alternate Protocol 2: Microscale immunoaffinity purification Basic Protocol 3: Separation of HLA peptides by reverse-phase high-performance liquid chromatography (RP-HPLC) Alternate Protocol 3: Isolation of HLA peptides using molecular weight cutoff (MWCO) filter Basic Protocol 4: Mass spectrometry and data analysis.

In-depth mining of the immunopeptidome of an acute myeloid leukemia cell line using complementary ligand enrichment and data acquisition strategies.

The identification of T cell epitopes derived from tumour specific antigens remains a significant challenge for the development of peptide-based vaccines and immunotherapies. The use of mass spectrometry-based approaches (immunopeptidomics) can provide powerful new avenues for the identification of such epitopes. In this study we report the use of complementary peptide antigen enrichment methods and a comprehensive mass spectrometric acquisition strategy to provide in-depth immunopeptidome data for the THP-1 cell line, a cell line used widely as a model of human leukaemia. To accomplish this, we combined robust experimental workflows that incorporated ultrafiltration or off-line reversed phase chromatography to enrich peptide ligand as well as a multifaceted data acquisition strategy using an Orbitrap Fusion LC-MS instrument. Using the combined datasets from the two ligand enrichment methods we gained significant depth in immunopeptidome coverage by identifying a total of 41,816 HLA class I peptides from THP-1 cells, including a significant number of peptides derived from different oncogenes or over expressed proteins associated with cancer. The physicochemical properties of the HLA-bound peptides dictated their recovery using the two ligand enrichment approaches and their distribution across the different precursor charge states considered in the data acquisition strategy. The data highlight the complementarity of the two enrichment procedures, and in cases where sample is not limiting, suggest that the combination of both approaches will yield the most comprehensive immunopeptidome information.

Mycobacterium tuberculosis MymA is a TLR2 agonist that activate macrophages and a TH1 response.

Cell wall of Mycobacterium tuberculosis (M.tb) is a major source of immunogenic proteins that can be tested as vaccine candidates. MymA (Rv3083), a 55 kDa M.tb flavin containing monooxygenase, is involved in modification of mycolic acids during acidic shock following M.tb internalization in macrophage. In this study, we have investigated the role of this cell wall associated protein in activation of macrophages by toll like receptor (TLRs) engagement and subsequent signaling. Our results showed that MymA stimulation of THP1 cells and human monocyte derived macrophages (MDM) lead to upregulation of TLR2 and co-stimulatory molecules CD40, CD80, CD86 and HLA-DR. This upregulation is partially reduced by TLR2 blocking antibodies. The activation of macrophage following MymA stimulation also resulted in release of proinflammatory cytokines, TNF-α and IL-12. Moreover, MymA also polarized the immune response towards TH1 as shown by an increased IFN-γ level in the supernatant of stimulated peripheral blood mononuclear cells (PBMC). In consensus with the TLR2 signaling involving MyD88 and NF-κB, we also observed several fold increase in mRNA for TLR2, MyD88 and NF-κB on MymA induction of THP-1 and MDM by qRT-PCR. The increased production of NF-κB following recognition of MymA by TLR2 was further confirmed by HEK-TLR2 reporter cell line colorimetric assay. In conclusion, immunological evaluation revealed that MymA is a TLR2 agonist that upregulates signaling via MyD88 and NF-κB in macrophages to stimulate the release of proinflammatory cytokines. The MymA protein should be investigated further for expression in recombinant BCG as a pre-exposure vaccine or as a post-exposure subunit vaccine candidate.

Characterization of MymA protein as a flavin-containing monooxygenase and as a target of isoniazid.

Tuberculosis is a global health problem especially with the emergence of drug-resistant Mycobacterium tuberculosis strains, creating an urgent need to identify new drug targets. The mycobacterial cell wall is an attractive target for chemotherapeutic agents. Gene products of mymA operon are known to be required for the maintenance of cell wall and play an important role in persistence, thus making them important drug targets. This study was undertaken to biochemically characterize the MymA as a flavin-containing monooxygenase (FMO). Our results established its enzymatic activity in vitro and found that the mycobacterial FMO requires NADPH and FAD as cofactors, similar to other characterized bacterial FMOs. The enzyme follows Michaelis-Menten kinetics to catalyze substrates such as trimethylamine and thiourea. We also propose that MymA could be one of the targets of the antituberculosis drug, isoniazid (INH), which is a cell wall inhibitor. Molecular docking studies revealed that INH targeted NADPH-binding site of the MymA. Further, experimental validation revealed that INH inhibits MymA with the IC50 value of 4.9 µm. Thus, this study characterizes for the first time that MymA is a mycobacterial FMO, which may be a target of INH.

Analysis of the DosR regulon genes to select cytotoxic T lymphocyte epitope specific vaccine candidates using a reverse vaccinology approach.

OBJECTIVE/BACKGROUND: There is an urgent need for a more effective vaccine against Mycobacterium tuberculosis (Mtb). Although CD4+ T cells play a central role in host immunity to Mtb, recent evidence suggests a critical role of CD8+ T cells in combating Mtb. In the present study, we have predicted HLA antigen class I binding peptides of DosR operon using an in-silico approach. This method is useful as an initial computational filtration of probable epitopes based on their binding ability and antigenicity. METHODS: CD8+ epitopes were predicted by software NetMHC 3.4 and BIMAS. Self-peptides were found and excluded by indigenously developed Perl script. Antigenicity of promiscuous peptides was predicted using a VaxiJen server. The top VaxiJen scoring antigenic peptides were docked to globally relevant HLA allele using CABS dock and Hex program. RESULTS: A total of 1436 overlapping nonamer peptides were generated which gave 46 promiscuous epitopes, 25 were predicted to be antigenic. Rv2627 epitope "SAFRPPLV" which gave the highest Vaxijen score of 1.9157 and showed binding to all the three HLA loci. The top VaxiJen scoring antigenic peptides were docked and had significant interactions with residues of the HLA class I molecule indicating them to be good cytotoxic T lymphocyte epitopes. CONCLUSION: Our study has generated several promiscuous antigenic peptides capable of binding to major histocompatibility complex class I with high affinity. These epitopes can become part of a postexposure multivalent subunit vaccine upon experimental validation.

Predicting promiscuous antigenic T cell epitopes of Mycobacterium tuberculosis mymA operon proteins binding to MHC Class I and Class II molecules.

Limited efficacy of Bacillus Calmette-Guérin vaccine has raised the need to explore other immunogenic candidates to develop an effective vaccine against Mycobacterium tuberculosis (Mtb). Both CD4+ and CD8+ T cells play a critical role in host immunity to Mtb. Infection of macrophages with Mtb results in upregulation of mymA operon genes thereby suggesting their importance as immune targets. In the present study, after exclusion of self-peptides mymA operon proteins of Mtb were analyzed in silico for the presence of Human Leukocyte Antigen (HLA) Class I and Class II binding peptides using Bioinformatics and molecular analysis section, NetMHC 3.4, ProPred and Immune epitope database software. Out of 56 promiscuous epitopes obtained, 41 epitopes were predicted to be antigenic for MHC Class I. In MHC Class II, out of 336 promiscuous epitopes obtained, 142 epitopes were predicted to be antigenic. The comparative bioinformatics analysis of mymA operon proteins found Rv3083 to be the best vaccine candidate. Molecular docking was performed with the most antigenic peptides of Rv3083 (LASGAASVV with alleles HLA-B51:01, HAATSGTLI with HLA-A02, IVTATGLNI and EKIHYGLKVNTA with HLA-DRB1_01:01) to study the structural basis for recognition of peptides by various HLA molecules. The software binding prediction was validated by the obtained molecular docking score of peptide-HLA complex. These peptides can be further investigated for their immunological relevance in patients of tuberculosis using major histocompatibility complex tetramer approach.