OpenCustomDB: Integration of Unannotated Open Reading Frames and Genetic Variants to Generate More Comprehensive Customized Protein Databases.

peer reviewedProteomic diversity in biological samples can be characterized by mass spectrometry (MS)-based proteomics using customized protein databases generated from sets of transcripts previously detected by RNA-seq. This diversity has only been increased by the recent discovery that many translated alternative open reading frames rest unannotated at unsuspected locations of mRNAs and ncRNAs. These novel protein products, termed alternative proteins, have been left out of all previous custom database generation tools. Consequently, genetic variations that impact alternative open reading frames and variant peptides from their translated proteins are not detectable with current computational workflows. To fill this gap, we present OpenCustomDB, a bioinformatics tool that uses sample-specific RNaseq data to identify genomic variants in canonical and alternative open reading frames, allowing for more than one coding region per transcript. In a test reanalysis of a cohort of 16 patients with acute myeloid leukemia, 5666 peptides from alternative proteins were detected, including 201 variant peptides. We also observed that a significant fraction of peptide-spectrum matches previously assigned to peptides from canonical proteins got better scores when reassigned to peptides from alternative proteins. Custom protein libraries that include sample-specific sequence variations of all possible open reading frames are promising contributions to the development of proteomics and precision medicine. The raw and processed proteomics data presented in this study can be found in PRIDE repository with accession number PXD029240.Canada Research Chairs Progra

Autophagy degrades immunogenic endogenous retroelements induced by 5-azacytidine in acute myeloid leukemia

The hypomethylating agent 5-azacytidine (AZA) is the first-line therapy for acute myeloid leukemia (AML) patients unfit for intensive chemotherapy. Evidence suggests that the anti-tumor effect of AZA results partly from T-cell cytotoxic responses against MHC-I-associated peptides (MAPs) whose expression is induced by hypomethylation. Through a proteogenomic approach, we analyzed the impact of AZA on the transcriptome and MAP repertoire of four AML cell lines and validated salient findings in the transcriptome of 437 primary AML samples. We demonstrate that AZA caused pleiotropic changes in AML cells via perturbation of transcription, translation, and protein degradation. Overall, 1,364 MAPs were upregulated in AZA-treated cells, including several cancer-testis antigens. Increased MAP abundance was due to the upregulation of corresponding transcripts in a minority of cases and post-translational events in most cases. Furthermore, AZA-induced hypomethylation increased the abundance of numerous transcripts, of which 38% were endogenous retroelements (EREs). Upregulated ERE transcripts triggered innate immune responses but were degraded by autophagy and not processed into MAPs. Autophagy resulted from the formation of protein aggregates caused by AZA-dependent inhibition of DNMT2, a tRNA-methyl transferase enzyme. We found that autophagy inhibition had a synergistic effect with AZA on AML cell proliferation and survival, increased ERE levels and triggered pro-inflammatory responses. Finally, autophagy gene signatures were associated with a lower abundance of CD8+T-cell markers in AML patients expressing high levels of EREs. Altogether, this work demonstrates that the impact of AZA is regulated at several levels and suggests that inhibiting autophagy could improve the immune recognition of AML blasts in patients

Autophagy degrades immunogenic endogenous retroelements induced by 5-azacytidine in acute myeloid leukemia.

peer reviewedThe hypomethylating agent 5-azacytidine (AZA) is the first-line treatment for AML patients unfit for intensive chemotherapy. The effect of AZA results in part from T-cell cytotoxic responses against MHC-I-associated peptides (MAPs) deriving from hypermethylated genomic regions such as cancer-testis antigens (CTAs), or endogenous retroelements (EREs). However, evidence supporting higher ERE MAPs presentation after AZA treatment is lacking. Therefore, using proteogenomics, we examined the impact of AZA on the repertoire of MAPs and their source transcripts. AZA-treated AML upregulated both CTA and ERE transcripts, but only CTA MAPs were presented at greater levels. Upregulated ERE transcripts triggered innate immune responses against double-stranded RNAs but were degraded by autophagy, and not processed into MAPs. Autophagy resulted from the formation of protein aggregates caused by AZA-dependent inhibition of DNMT2. Autophagy inhibition had an additive effect with AZA on AML cell proliferation and survival, increased ERE levels, increased pro-inflammatory responses, and generated immunogenic tumor-specific ERE-derived MAPs. Finally, autophagy was associated with a lower abundance of CD8+ T-cell markers in AML patients expressing high levels of EREs. This work demonstrates that AZA-induced EREs are degraded by autophagy and shows that inhibiting autophagy can improve the immune recognition of AML blasts in treated patients

Abstract 2993: Unmutated tumor antigens are abundant and contribute to tumor control in melanoma

peer reviewedAbstract Recognition of MHC-I-associated tumor antigens (TAs) by CD8+ T cells is central to antitumor immunity. Owing to the elevated tumor mutational burden (TMB) in melanoma, the marked efficacy of immune checkpoint blockade (ICB) has been attributed to the recognition of mutated TAs. However, recent reports showed that response to ICB in melanomas with low TMB is associated with CD8+ T-cell reactivity against melanocyte lineage-associated antigens (LSAs). Here, we systematically evaluated the contribution of all TA classes, i.e., mutated and unmutated, canonical and non-canonical, to the antigenic landscape of melanoma. We characterized the TAs from melanoma biopsies and patient-derived cell lines using proteogenomics. Out of 79450 MHC-I-associated peptides (MAPs) identified from 19 samples, we found 557 unmutated TAs classified as tumor-specific (TSA), tumor-associated (TAA), or LSAs. These TAs most often derived from annotated open-reading frames, followed by ncRNAs and intergenic regions. By contrast, only 6 MAPs were mutated and tumor-specific, which could be partially explained by a decreased expression of mutations within MAP-generating genomic regions. While the number of unmutated TAs with predicted presentation (TApres) in melanoma patients was similar between responders and non-responders pre-ICB, non-responders showed marks of inefficient antigen presentation. In consequence, only responders lost TApres upon treatment, in tandem with an expansion in tumor-infiltrating lymphocytes. These results reveal a previously underappreciated contribution of unmutated TAs to tumor control in melanoma and suggest that enhancing their recognition could improve the ICB efficacy in non-responders. Citation Format: Anca Apavaloaei, Qingchuan Zhao, Leslie Hesnard, Krystel Vincent, Marie-Pierre Hardy, Chantal Durette, Joël Lanoix, Jean-Philippe Laverdure, Jean-David Larouche, Maria Virginia Ruiz Cuevas, Grégory Ehx, Sébastien Lemieux, Pierre Thibault, Claude Perreault. Unmutated tumor antigens are abundant and contribute to tumor control in melanoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 2993

The tumor-specific antigen landscape of acute myeloid leukemia

peer reviewedIn theory, tumor-specific antigens (TSAs) represent ideal targets for cancer immunotherapy, but none have been identified in acute myeloid leukemia (AML) so far. We have recently designed a proteogenomic approach allowing the identification of MHC class I-associated peptides deriving from any reading frame of all genomic regions (not only from exons). Using this approach, we could identify 82 TSAs from 18 primary human AML samples. All these TSAs derived from unmutated transcripts which were aberrantly expressed by allegedly non-coding regions of the genome (that would have been missed by conventional exome-based approaches). In agreement with previous reports showing that high intron retention is a key feature in AML, we found that introns were the main source of TSAs (54.8%). The TSAs were presented by 24 different HLA allotypes, together covering 98,04% of the world population. The examination of the expression of the RNAs coding for the TSAs showed that most of them (68%) were expressed by at least 50% of AML samples in the “The Cancer Genome Atlas” (TCGA) cohort while all of them were either expressed at very low levels or not expressed at all in normal tissues (from the “Genotype-Tissue Expression” (GTEx) project). We have yet to discover the specific epigenetic changes leading to TSA expression in cancer cells. Meanwhile, we propose that vaccination against ab- errantly expressed TSAs could be used to treat a vast majority of patients, with minimal risks of collateral damage to healthy normal tissues

Induced pluripotent stem cells display a distinct set of MHC I-associated peptides shared by human cancers.

peer reviewedPrevious reports showed that mouse vaccination with pluripotent stem cells (PSCs) induces durable anti-tumor immune responses via T cell recognition of some elusive oncofetal epitopes. We characterize the MHC I-associated peptide (MAP) repertoire of human induced PSCs (iPSCs) using proteogenomics. Our analyses reveal a set of 46 pluripotency-associated MAPs (paMAPs) absent from the transcriptome of normal tissues and adult stem cells but expressed in PSCs and multiple adult cancers. These paMAPs derive from coding and allegedly non-coding (48%) transcripts involved in pluripotency maintenance, and their expression in The Cancer Genome Atlas samples correlates with source gene hypomethylation and genomic aberrations common across cancer types. We find that several of these paMAPs were immunogenic. However, paMAP expression in tumors coincides with activation of pathways instrumental in immune evasion (WNT, TGF-β, and CDK4/6). We propose that currently available inhibitors of these pathways could synergize with immune targeting of paMAPs for the treatment of poorly differentiated cancers

Breast cancer immunopeptidomes contain numerous shared tumor antigens

peer reviewedBackground Hormone-receptor-positive breast cancer (HR+) is an immunologically cold cancer that has not benefited from advances in immunotherapy. In contrast, triple-negative breast cancer (TNBC) displays high levels of leukocytic infiltration and responds to immune checkpoint inhibitors. CD8 T cells, the main effectors of anti-cancer responses, recognize MHC I-associated peptides (MAPs). Our work aimed to characterize the repertoire of MAPs presented by HR+ and TNBC tumors. Methods Using a proteogenomic approach relying on mass spectrometry, we identified 57 094 unique MAPs in 26 primary breast cancer samples (14 HR+, 12 TNBC). Results MAP source genes showed a high overlap between both subtypes (>70%). We identified 25 tumor-specific antigens (TSAs) derived from various genomic regions, of which 24 were unmutated. TSAs were mainly identified in TNBC samples (70%) and were more highly shared among TCGA TNBC than HR+ samples. In the TNBC TCGA cohort, the predicted number of TSAs positively correlated with leukocytic infiltration (p<0.05) and overall survival (p<0.05, figure 1), suggesting that these TSAs are immunogenic in vivo. We also identified 49 overexpressed tumor-associated antigens (TAAs), some of which derived from cancer-associated fibroblasts. FEST assays confirmed the in vitro immunogenicity of our TSAs and TAAs. Conclusions Well-defined antigens were identified in both subtypes of breast cancer and represent attractive targets for cancer immunotherapy. The higher prevalence and immunogenicity of TSAs in TNBC tumors provide a molecular rationale for the responsiveness of TNBC to immune checkpoint inhibitors

Breast cancer immunopeptidomes contain numerous shared tumor antigens.

peer reviewedHormone receptor-positive breast cancer (HR+) is immunologically cold and has not benefited from advances in immunotherapy. In contrast, subsets of triple-negative breast cancer (TNBC) display high leukocytic infiltration and respond to checkpoint blockade. CD8+ T cells, the main effectors of anticancer responses, recognize MHC I-associated peptides (MAPs). Our work aimed to characterize the repertoire of MAPs presented by HR+ and TNBC tumors. Using mass spectrometry, we identified 57,094 unique MAPs in 26 primary breast cancer samples. MAP source genes highly overlapped between both subtypes. We identified 25 tumor-specific antigens (TSAs) mainly deriving from aberrantly expressed regions. TSAs were most frequently identified in TNBC samples and were more shared among The Cancer Genome Atlas (TCGA) database TNBC than HR+ samples. In the TNBC cohort, the predicted number of TSAs positively correlated with leukocytic infiltration and overall survival, supporting their immunogenicity in vivo. We detected 49 tumor-associated antigens (TAAs), some of which derived from cancer-associated fibroblasts. Functional expansion of specific T cell assays confirmed the in vitro immunogenicity of several TSAs and TAAs. Our study identified attractive targets for cancer immunotherapy in both breast cancer subtypes. The higher prevalence of TSAs in TNBC tumors provides a rationale for their responsiveness to checkpoint blockade