MHC class I–associated peptides (MAPs) define the immune self for CD8⁺ T lymphocytes and are key targets of cancer immunosurveillance. Here, the goals of our work were to determine whether the entire set of protein-coding genes could generate MAPs and whether specific features influence the ability of discrete genes to generate MAPs. Using proteogenomics, we have identified 25,270 MAPs isolated from the B lymphocytes of 18 individuals who collectively expressed 27 high-frequency HLA-A,B allotypes. The entire MAP repertoire presented by these 27 allotypes covered only 10% of the exomic sequences expressed in B lymphocytes. Indeed, 41% of expressed protein-coding genes generated no MAPs, while 59% of genes generated up to 64 MAPs, often derived from adjacent regions and presented by different allotypes. We next identified several features of transcripts and proteins associated with efficient MAP production. From these data, we built a logistic regression model that predicts with good accuracy whether a gene generates MAPs. Our results show preferential selection of MAPs from a limited repertoire of proteins with distinctive features. The notion that the MHC class I immunopeptidome presents only a small fraction of the protein-coding genome for monitoring by the immune system has profound implications in autoimmunity and cancer immunology.

Colorectal cancer is the second leading cause of cancer death worldwide, and the incidence of this disease is expected to increase as global socioeconomic changes occur. Immune checkpoint inhibition therapy is effective in treating a minority of colorectal cancer tumors; however, microsatellite stable tumors do not respond well to this treatment. Emerging cancer immunotherapeutic strategies aim to activate a cytotoxic T cell response against tumor-specific antigens, presented exclusively at the cell surface of cancer cells. These antigens are rare and are most effectively identified with a mass spectrometry–based approach, which allows the direct sampling and sequencing of these peptides. Although the few tumor-specific antigens identified to date are derived from coding regions of the genome, recent findings indicate that a large proportion of tumor-specific antigens originate from allegedly noncoding regions. Here, we employed a novel proteogenomic approach to identify tumor antigens in a collection of colorectal cancer–derived cell lines and biopsy samples consisting of matched tumor and normal adjacent tissue. The generation of personalized cancer databases paired with mass spectrometry analyses permitted the identification of more than 30,000 unique MHC I–associated peptides. We identified 19 tumor-specific antigens in both microsatellite stable and unstable tumors, over two-thirds of which were derived from noncoding regions. Many of these peptides were derived from source genes known to be involved in colorectal cancer progression, suggesting that antigens from these genes could have therapeutic potential in a wide range of tumors. These findings could benefit the development of T cell–based vaccines, in which T cells are primed against these antigens to target and eradicate tumors. Such a vaccine could be used in tandem with existing immune checkpoint inhibition therapies, to bridge the gap in treatment efficacy across subtypes of colorectal cancer with varying prognoses. Data are available via ProteomeXchange with identifier PXD028309.

Hormone 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.

Epithelial ovarian cancer (EOC) has not significantly benefited from advances in immunotherapy, mainly because of the lack of well-defined actionable antigen targets. Using proteogenomic analyses of primary EOC tumors, we previously identified 91 aberrantly expressed tumor-specific antigens (TSAs) originating from unmutated genomic sequences. Most of these TSAs derive from non-exonic regions, and their expression results from cancer-specific epigenetic changes. The present study aimed to evaluate the immunogenicity of 48 TSAs selected according to two criteria: presentation by highly prevalent HLA allotypes and expression in a significant fraction of EOC tumors. Using targeted mass spectrometry analyses, we found that pulsing with synthetic TSA peptides leads to a high-level presentation on dendritic cells. TSA abundance correlated with the predicted binding affinity to the HLA allotype. We stimulated naïve CD8 T cells from healthy blood donors with TSA-pulsed dendritic cells and assessed their expansion with two assays: MHC-peptide tetramer staining and TCR Vβ CDR3 sequencing. We report that these TSAs can expand sizeable populations of CD8 T cells and, therefore, represent attractive targets for EOC immunotherapy.

Treatment with immune checkpoint inhibitors in colorectal cancer (CRC) has largely benefited patients with microsatellite instability-high (MSI-H) and not the larger proportion of patient with microsatellite-stable (MSS) tumors. This clinical dichotomy has fueled the view that high mutational burden is the dominant driver of tumor immunogenicity and that MSS CRC fails to respond because it is "antigen poor". To directly test this premise and define the origins of presented tumor antigens, we integrated HLA class I immunopeptidomics and matched RNA-seq from 26 primary CRC tumors spanning MSI-H and MSS subtypes. Using patient-specific canonical and cancer-specific proteogenomic databases, we identified 115,292 unique MHC-associated peptides (MAPs) across 61 HLA alleles, with a mean of 9,292 MAPs per tumor and no significant difference in MAP counts between MSI-H and MSS tumors. In toto, we identified 266 tumor antigens, all coded by unmutated genomic sequences, comprising 70 aberrantly expressed tumor-specific antigens (aeTSAs) and 196 tumor-associated antigens (TAAs). In our cohort, MSS tumors presented more TAAs and a comparable number of aeTSAs per tumor relative to MSI-H tumors. In TCGA-COAD stratified analyses (483 tumors), MSS tumors yielded more presentable aeTSAs and TAAs per patient than MSI-H tumors. Across both subtypes, aeTSAs arose predominantly from intronic translation, UTR usage, retroelement activation, and germline-like transcription, including recurrent aeTSAs from PIWIL1, L1TD1, and endogenous retroviral loci. Together, these data demonstrate that MSS CRC is not antigen poor and highlight non-canonical translation as a major, previously underappreciated contributor to the CRC immunopeptidome.

MHC-I-associated peptides (MAPs) derive from selective yet highly diverse genomic regions, including allegedly non-protein-coding sequences, such as endogenous retroelements (EREs). Quantifying canonical (exonic) and non-canonical MAPs-encoding RNA expression in malignant and benign cells is critical for identifying tumor antigens (TAs) but represents a challenge for immunologists. We present BamQuery, a computational tool attributing an exhaustive RNA expression to MAPs of any origin (exon, intron, UTR, intergenic) from bulk and single-cell RNA-sequencing data. We show that non-canonical MAPs (including TAs) can derive from multiple different genomic regions (up to 35,343 for EREs), abundantly expressed in normal tissues. We also show that supposedly tumor-specific mutated MAPs, viral MAPs, and MAPs derived from proteasomal splicing can arise from different unmutated non-canonical genomic regions. The genome-wide approach of BamQuery allows comprehensive mapping of all MAPs in healthy and cancer tissues. BamQuery can also help predict MAP immunogenicity and identify safe and actionable TAs.