Strand Orientation Bias Detector to determine the probability of FFPE sequencing artifacts.

Formalin-fixed paraffin-embedded tissue, the most common tissue specimen stored in clinical practice, presents challenges in the analysis due to formalin-induced artifacts. Here, we present Strand Orientation Bias Detector (SOBDetector), a flexible computational platform compatible with all the common somatic SNV-calling pipelines, designed to assess the probability whether a given detected mutation is an artifact. The underlying predictor mechanism is based on the posterior distribution of a Bayesian logistic regression model trained on The Cancer Genome Atlas whole exomes. SOBDetector is a freely available cross-platform program, implemented in Java 1.8.

Transcriptomic signatures of tumors undergoing T cell attack.

BACKGROUND: Studying tumor cell-T cell interactions in the tumor microenvironment (TME) can elucidate tumor immune escape mechanisms and help predict responses to cancer immunotherapy. METHODS: We selected 14 pairs of highly tumor-reactive tumor-infiltrating lymphocytes (TILs) and autologous short-term cultured cell lines, covering four distinct tumor types, and co-cultured TILs and tumors at sub-lethal ratios in vitro to mimic the interactions occurring in the TME. We extracted gene signatures associated with a tumor-directed T cell attack based on transcriptomic data of tumor cells. RESULTS: An autologous T cell attack induced pronounced transcriptomic changes in the attacked tumor cells, partially independent of IFN-γ signaling. Transcriptomic changes were mostly independent of the tumor histological type and allowed identifying common gene expression changes, including a shared gene set of 55 transcripts influenced by T cell recognition (Tumors undergoing T cell attack, or TuTack, focused gene set). TuTack scores, calculated from tumor biopsies, predicted the clinical outcome after anti-PD-1/anti-PD-L1 therapy in multiple tumor histologies. Notably, the TuTack scores did not correlate to the tumor mutational burden, indicating that these two biomarkers measure distinct biological phenomena. CONCLUSIONS: The TuTack scores measure the effects on tumor cells of an anti-tumor immune response and represent a comprehensive method to identify immunologically responsive tumors. Our findings suggest that TuTack may allow patient selection in immunotherapy clinical trials and warrant its application in multimodal biomarker strategies.

The genomic imprint of cancer therapies helps timing the formation of metastases.

A retrospective determination of the time of metastasis formation is essential for a better understanding of the evolution of oligometastatic cancer. This study was based on the hypothesis that genomic alterations induced by cancer therapies could be used to determine the temporal order of the treatment and the formation of metastases. We analysed the whole genome sequence of a primary tumour sample and three metastatic sites derived from autopsy samples from a young never-smoker lung adenocarcinoma patient with an activating EGFR mutation. Mutation detection methods were refined to accurately detect and distinguish clonal and subclonal mutations. In comparison to a panel of samples from untreated smoker or never-smoker patients, we showed that the mutagenic effect of cisplatin treatment could be specifically detected from the base substitution mutations. Metastases that arose before or after chemotherapeutic treatment could be distinguished based on the allele frequency of cisplatin-induced dinucleotide mutations. In addition, genomic rearrangements and late amplification of the EGFR gene likely induced by afatinib treatment following the acquisition of a T790M gefitinib resistance mutation provided further evidence to tie the time of metastasis formation to treatment history. The established analysis pipeline for the detection of treatment-derived mutations allows the drawing of tumour evolutionary paths based on genomic data, showing that metastases may be seeded well before they become detectable by clinical imaging.

Dynamic classification using case-specific training cohorts outperforms static gene expression signatures in breast cancer.

The molecular diversity of breast cancer makes it impossible to identify prognostic markers that are applicable to all breast cancers. To overcome limitations of previous multigene prognostic classifiers, we propose a new dynamic predictor: instead of using a single universal training cohort and an identical list of informative genes to predict the prognosis of new cases, a case-specific predictor is developed for each test case. Gene expression data from 3,534 breast cancers with clinical annotation including relapse-free survival is analyzed. For each test case, we select a case-specific training subset including only molecularly similar cases and a case-specific predictor is generated. This method yields different training sets and different predictors for each new patient. The model performance was assessed in leave-one-out validation and also in 325 independent cases. Prognostic discrimination was high for all cases (n = 3,534, HR = 3.68, p = 1.67 E-56). The dynamic predictor showed higher overall accuracy (0.68) than genomic surrogates for Oncotype DX (0.64), Genomic Grade Index (0.61) or MammaPrint (0.47). The dynamic predictor was also effective in triple-negative cancers (n = 427, HR = 3.08, p = 0.0093) where the above classifiers all failed. Validation in independent patients yielded similar classification power (HR = 3.57). The dynamic classifier is available online at http://www.recurrenceonline.com/?q=Re_training. In summary, we developed a new method to make personalized prognostic prediction using case-specific training cohorts. The dynamic predictors outperform static models developed from single historical training cohorts and they also predict well in triple-negative cancers.

Are tumor-associated carbohydrates the missing link between the gut microbiome and response to immune checkpoint inhibitor treatment in cancer?

Immune checkpoint inhibitor therapy has dramatically improved survival in a significant subset of patients with several solid tumor types. Increasing the number of patients benefitting from this form of therapy is an important translational research goal. Correlations between the composition of the gut microbiome and response to immune checkpoint inhibitor therapy raised the possibility that direct modulation of the gut microbiome may significantly improve the clinical benefit of this treatment. Several lines of observations suggest that tumor-associated carbohydrates, including those recognized as blood group-related glycolipid antigens, such as the Forssman antigen, may be some of the key factors behind this clinical correlation. Such antigens are expressed in human cancer, humans often produce antibodies against those, and they can induce antibody directed cellular cytotoxicity. Importantly, these antibodies are often induced by antigens present in microbes of the gut. If identified, these antibodies could be boosted by appropriate vaccination techniques and thus enhance anti-tumor immunity with minimal side effects.

Histopathology and proteomics are synergistic for High-Grade Serous Ovarian Cancer platinum response prediction.

Patients with High-Grade Serous Ovarian Cancer (HGSOC) exhibit varied responses to treatment, with 20-30% showing de novo resistance to platinum-based chemotherapy. While hematoxylin-eosin (H&E) pathological slides are used for routine diagnosis of cancer type, they may also contain diagnostically useful information about treatment response. Our study demonstrates that combining H&E-stained Whole Slide Images (WSIs) with proteomic signatures using a multimodal deep learning framework significantly improves the prediction of platinum response in both discovery and validation cohorts. This method outperforms the Homologous Recombination Deficiency (HRD) score in predicting platinum response and overall patient survival. The study sets new performance benchmarks and explores the intersection of histology and proteomics, highlighting phenotypes related to treatment response pathways, including homologous recombination, DNA damage response, nucleotide synthesis, apoptosis, and ER stress. This integrative approach has the potential to improve personalized treatment and provide insights into the therapeutic vulnerabilities of HGSOC.

Mutational signature-based identification of DNA repair deficient gastroesophageal adenocarcinomas for therapeutic targeting.

Homologous recombination (HR) and nucleotide excision repair (NER) are the two most frequently disabled DNA repair pathways in cancer. HR-deficient breast, ovarian, pancreatic and prostate cancers respond well to platinum chemotherapy and PARP inhibitors. However, the frequency of HR deficiency in gastric and esophageal adenocarcinoma (GEA) still lacks diagnostic and functional validation. Using whole exome and genome sequencing data, we found that a significant subset of GEA, but very few colorectal adenocarcinomas, show evidence of HR deficiency by mutational signature analysis (HRD score). High HRD gastric cancer cell lines demonstrated functional HR deficiency by RAD51 foci assay and increased sensitivity to platinum chemotherapy and PARP inhibitors. Of clinical relevance, analysis of three different GEA patient cohorts demonstrated that platinum treated HR deficient cancers had better outcomes. A gastric cancer cell line with strong sensitivity to cisplatin showed HR proficiency but exhibited NER deficiency by two photoproduct repair assays. Single-cell RNA-sequencing revealed that, in addition to inducing apoptosis, cisplatin treatment triggered ferroptosis in a NER-deficient gastric cancer, validated by intracellular GSH assay. Overall, our study provides preclinical evidence that a subset of GEAs harbor genomic features of HR and NER deficiency and may therefore benefit from platinum chemotherapy and PARP inhibitors.

Increased frequency of CHD1 deletions in prostate cancers of African American men is associated with rapid disease progression without inducing homologous recombination deficiency.

We analyzed genomic data derived from the prostate cancer of African and European American men in order to identify differences that may contribute to racial disparity of outcome and that could also define novel therapeutic strategies. In addition to analyzing patient derived next generation sequencing data, we performed FISH based confirmatory studies of Chromodomain helicase DNA-binding protein 1 (CHD1) loss on prostate cancer tissue microarrays. We created CRISPR edited, CHD1 deficient prostate cancer cell lines for genomic, drug sensitivity and functional homologous recombination (HR) activity analysis. We found that subclonal deletion of CHD1 is nearly three times as frequent in prostate tumors of African American men than in men of European ancestry and it associates with rapid disease progression. We further showed that CHD1 deletion is not associated with homologous recombination deficiency associated mutational signatures in prostate cancer. In prostate cancer cell line models CHD1 deletion did not induce HR deficiency as detected by RAD51 foci formation assay or mutational signatures, which was consistent with the moderate increase of olaparib sensitivity. CHD1 deficient prostate cancer cells, however, showed higher sensitivity to talazoparib. CHD1 loss may contribute to worse outcome of prostate cancer in African American men. A deeper understanding of the interaction between CHD1 loss and PARP inhibitor sensitivity will be needed to determine the optimal use of targeted agents such as talazoparib in the context of castration resistant prostate cancer.

[Simplified, low-cost gene expression profiling for the prediction of outcome in breast cancer based on routine histologic specimens]. / Egyszerusített, alacsony költségu génexpressziós teszt az emlodaganatos betegségek kimenetelének elorejelzésére rutin patológiai blokkok felhasználásával.

BACKGROUND: Grade 2 breast carcinomas do not form a uniform prognostic group. AIM: To extend the number of patients and the investigated genes of a previously identified prognostic signature described by the authors that reflect chromosomal instability in order to refine characterization of grade 2 breast cancers and identify driver genes. METHODS: Using publicly available databases, the authors selected 9 target and 3 housekeeping genes that are capable to divide grade 2 breast carcinomas into prognostic groups. Gene expression was investigated by polymerase chain reaction in 249 formalin-fixed, paraffin-embedded breast tumors. The results were correlated with relapse-free survival. RESULTS: Histologically grade 2 carcinomas were split into good and a poor prognosis groups. Centroid-based ranking showed that 3 genes, FOXM1, TOP2A and CLDN4 were able to separate the good and poor prognostic groups of grade 2 breast carcinomas. CONCLUSION: Using appropriately selected control genes, a limited set of genes is able to split prognostic groups of breast carcinomas independently from their grade.

Intratumoral T-cell and B-cell receptor architecture associates with distinct immune tumor microenvironment features and clinical outcomes of anti-PD-1/L1 immunotherapy.

BACKGROUND: Effective cooperation between B-cells and T-cells within the tumor microenvironment may lead to the regression of established tumors. B-cells and T-cells can recognize tumor antigens with exquisite specificity via their receptor complexes. Nevertheless, whether a diverse intratumoral B-cells and T-cell receptor (BCR, TCR) repertoire affects the tumor immune microenvironment (TIME) and clinical outcomes in patients treated with immunotherapy is unclear. METHODS: We extracted information on BCR and TCR repertoire diversity from large clinical datasets and measured the association between immune receptor diversity features, the TIME, and clinical outcomes of patients treated with anti-PD-1/PD-L1 immunotherapy. RESULTS: In multiple tumor types, an increasingly diverse TCR repertoire was strongly associated with a highly activated TIME, while BCR diversity was more associated with antibody responses but not with the overall B-cell infiltration nor with measures related to intratumoral CD8+T cell activity. Neither TCR nor BCR diversity was independent prognostic biomarkers of survival across multiple cancer types. However, both TCR and BCR diversity improved the performance of predictive models combined with established biomarkers of response to immunotherapy. CONCLUSION: Overall, these data indicate a currently unexplored immunological role of intratumoral B-cells associated with BCR diversity and antibody responses but independent of classical anticancer T-cells intratumoral activities.

Activation of invasion by oncogenic reprogramming of cholesterol metabolism via increased NPC1 expression and macropinocytosis.

Cancer cells are dependent on cholesterol, and they possess strictly controlled cholesterol homeostasis mechanisms. These allow them to smoothly switch between cholesterol synthesis and uptake to fulfill their needs and to adapt environmental changes. Here we describe a mechanism of how cancer cells employ oncogenic growth factor signaling to promote uptake and utilization of extracellular cholesterol via Myeloid Zinc Finger 1 (MZF1)-mediated Niemann Pick C1 (NPC1) expression and upregulated macropinocytosis. Expression of p95ErbB2, highly oncogenic, standard-treatment resistant form of ErbB2 mobilizes lysosomes and activates EGFR, invasion and macropinocytosis. This is connected to a metabolic shift from cholesterol synthesis to uptake due to macropinocytosis-enabled flow of extracellular cholesterol. NPC1 increase facilitates extracellular cholesterol uptake and is necessary for the invasion of ErbB2 expressing breast cancer spheroids and ovarian cancer organoids, indicating a regulatory role for NPC1 in the process. The ability to obtain cholesterol as a byproduct of increased macropinocytosis allows cancer cells to direct the resources needed for the energy-consuming cholesterol synthesis towards other activities such as invasion. These results demonstrate that macropinocytosis is not only an alternative energy source for cancer cells but also an efficient way to provide building material, such as cholesterol, for its macromolecules and membranes.

Prognostic value of circulating proteins in patients undergoing surgery for pancreatic cancer.

BACKGROUND: Pancreatic ductal adenocarcinoma (PDAC) is the fourth leading cause of cancer death. Less than 20% of patients are diagnosed with resectable disease. Identifying truly resectable disease is challenging because 20%-40% of the patients subjected to resection are found to have advanced disease during surgery. The aim of our study was to identify panels of circulating proteins that could be used to distinguish patients with unresectable PDAC from patients with resectable PDAC and to identify prognostic signatures for both groups. METHODS: We measured 92 circulating immuno-oncology-related proteins using the proximity extension assay from Olink Proteomics in 273 patients eligible for surgery for PDAC. Two bioinformaticians worked independently of one another on the same data. LASSO and Ridge regression were used in the statistical analyses. RESULTS: One protein index for determining resectability had an AUC value of 0.66. Several indices for prognosis had AUC values between 0.50 and 0.75 and were therefore not better than existing prognostic markers. DISCUSSION: Our study did not reveal any new high-performing protein panels that could be used to identify patients with inoperable PDAC before surgery. The panel of 92 proteins investigated has previously been found to be applicable for diagnostic use in patients with PDAC, but it does not seem to warrant further investigation regarding resectability in the subgroup of patients with PDAC referred to surgery.

Replicative stress in gastroesophageal adenocarcinoma is associated with chromosomal instability and sensitivity to DNA damage response inhibitors.

Gastroesophageal adenocarcinoma (GEA) is an aggressive, often lethal, malignancy that displays marked chromosomal instability (CIN). To understand adaptive responses that enable CIN, we analyzed paired normal, premalignant, and malignant gastric lesions from human specimens and a carcinogen-induced mouse model, observing activation of replication stress, DNA damage response (DDR), and cell cycle regulator p21 in neoplastic progression. In GEA cell lines, expression of DDR markers correlated with ploidy abnormalities, including high-level focal amplifications and whole-genome duplication (WGD). Moreover, high expression of DNA damage marker H2AX correlated with CIN, WGD, and inferior patient survival. By developing and implementing a composite diagnostic score that incorporates TP53 mutation status, ploidy abnormalities, and H2AX expression, among other genomic information, we can identify GEA cell lines with enhanced sensitivity to DDR pathway inhibitors targeting Chk1/2 and Wee1. Anti-tumor properties were further augmented in combination with irinotecan (SN38) but not gemcitabine chemotherapy. These results implicate specific DDR biomarkers and ploidy abnormalities as diagnostic proxy that may predict premalignant progression and response to DDR pathway inhibitors.

Replicative stress in gastroesophageal cancer is associated with chromosomal instability and sensitivity to DNA damage response inhibitors.

Gastroesophageal adenocarcinoma (GEA) is an aggressive malignancy with chromosomal instability (CIN). To understand adaptive responses enabling DNA damage response (DDR) and CIN, we analyzed matched normal, premalignant, and malignant gastric lesions from human specimens and a carcinogen-induced mouse model, observing activation of replication stress, DDR, and p21 in neoplastic progression. In GEA cell lines, expression of DDR markers correlated with ploidy abnormalities, such as number of high-level focal amplifications and whole-genome duplication (WGD). Integrating TP53 status, ploidy abnormalities, and DDR markers into a compositive score helped predict GEA cell lines with enhanced sensitivity to Chk1/2 and Wee1 inhibition, either alone or combined with irinotecan (SN38). We demonstrate that Chk1/2 or Wee1 inhibition combined with SN38/irinotecan shows greater anti-tumor activity in human gastric cancer organoids and an in vivo xenograft mouse model. These findings indicate that specific DDR biomarkers and ploidy abnormalities may predict premalignant progression and response to DDR pathway inhibitors.

Nucleotide excision repair deficiency is a targetable therapeutic vulnerability in clear cell renal cell carcinoma.

Purpose: Due to a demonstrated lack of DNA repair deficiencies, clear cell renal cell carcinoma (ccRCC) has not benefitted from targeted synthetic lethality-based therapies. We investigated whether nucleotide excision repair (NER) deficiency is present in an identifiable subset of ccRCC cases that would render those tumors sensitive to therapy targeting this specific DNA repair pathway aberration. Experimental Design: We used functional assays that detect UV-induced 6-4 pyrimidine-pyrimidone photoproducts to quantify NER deficiency in ccRCC cell lines. We also measured sensitivity to irofulven, an experimental cancer therapeutic agent that specifically targets cells with inactivated transcription-coupled nucleotide excision repair (TC-NER). In order to detect NER deficiency in clinical biopsies, we assessed whole exome sequencing data for the presence of an NER deficiency associated mutational signature previously identified in ERCC2 mutant bladder cancer. Results: Functional assays showed NER deficiency in ccRCC cells. Irofulven sensitivity increased in some cell lines. Prostaglandin reductase 1 (PTGR1), which activates irofulven, was also associated with this sensitivity. Next generation sequencing data of the cell lines showed NER deficiency-associated mutational signatures. A significant subset of ccRCC patients had the same signature and high PTGR1 expression. Conclusions: ccRCC cell line based analysis showed that NER deficiency is likely present in this cancer type. Approximately 10% of ccRCC patients in the TCGA cohort showed mutational signatures consistent with ERCC2 inactivation associated NER deficiency and also substantial levels of PTGR1 expression. These patients may be responsive to irofulven, a previously abandoned anticancer agent that has minimal activity in NER-proficient cells.

CRISPR screens reveal genetic determinants of PARP inhibitor sensitivity and resistance in prostate cancer.

Prostate cancer harboring BRCA1/2 mutations are often exceptionally sensitive to PARP inhibitors. However, genomic alterations in other DNA damage response genes have not been consistently predictive of clinical response to PARP inhibition. Here, we perform genome-wide CRISPR-Cas9 knockout screens in BRCA1/2-proficient prostate cancer cells and identify previously unknown genes whose loss has a profound impact on PARP inhibitor response. Specifically, MMS22L deletion, frequently observed (up to 14%) in prostate cancer, renders cells hypersensitive to PARP inhibitors by disrupting RAD51 loading required for homologous recombination repair, although this response is TP53-dependent. Unexpectedly, loss of CHEK2 confers resistance rather than sensitivity to PARP inhibition through increased expression of BRCA2, a target of CHEK2-TP53-E2F7-mediated transcriptional repression. Combined PARP and ATR inhibition overcomes PARP inhibitor resistance caused by CHEK2 loss. Our findings may inform the use of PARP inhibitors beyond BRCA1/2-deficient tumors and support reevaluation of current biomarkers for PARP inhibition in prostate cancer.

Nucleotide excision repair deficiency is a targetable therapeutic vulnerability in clear cell renal cell carcinoma.

Due to a demonstrated lack of DNA repair deficiencies, clear cell renal cell carcinoma (ccRCC) has not benefitted from targeted synthetic lethality-based therapies. We investigated whether nucleotide excision repair (NER) deficiency is present in an identifiable subset of ccRCC cases that would render those tumors sensitive to therapy targeting this specific DNA repair pathway aberration. We used functional assays that detect UV-induced 6-4 pyrimidine-pyrimidone photoproducts to quantify NER deficiency in ccRCC cell lines. We also measured sensitivity to irofulven, an experimental cancer therapeutic agent that specifically targets cells with inactivated transcription-coupled nucleotide excision repair (TC-NER). In order to detect NER deficiency in clinical biopsies, we assessed whole exome sequencing data for the presence of an NER deficiency associated mutational signature previously identified in ERCC2 mutant bladder cancer. Functional assays showed NER deficiency in ccRCC cells. Some cell lines showed irofulven sensitivity at a concentration that is well tolerated by patients. Prostaglandin reductase 1 (PTGR1), which activates irofulven, was also associated with this sensitivity. Next generation sequencing data of the cell lines showed NER deficiency-associated mutational signatures. A significant subset of ccRCC patients had the same signature and high PTGR1 expression. ccRCC cell line-based analysis showed that NER deficiency is likely present in this cancer type. Approximately 10% of ccRCC patients in the TCGA cohort showed mutational signatures consistent with ERCC2 inactivation associated NER deficiency and also substantial levels of PTGR1 expression. These patients may be responsive to irofulven, a previously abandoned anticancer agent that has minimal activity in NER-proficient cells.

A biallelic multiple nucleotide length polymorphism explains functional causality at 5p15.33 prostate cancer risk locus.

To date, single-nucleotide polymorphisms (SNPs) have been the most intensively investigated class of polymorphisms in genome wide associations studies (GWAS), however, other classes such as insertion-deletion or multiple nucleotide length polymorphism (MNLPs) may also confer disease risk. Multiple reports have shown that the 5p15.33 prostate cancer risk region is a particularly strong expression quantitative trait locus (eQTL) for Iroquois Homeobox 4 (IRX4) transcripts. Here, we demonstrate using epigenome and genome editing that a biallelic (21 and 47 base pairs (bp)) MNLP is the causal variant regulating IRX4 transcript levels. In LNCaP prostate cancer cells (homozygous for the 21 bp short allele), a single copy knock-in of the 47 bp long allele potently alters the chromatin state, enabling de novo functional binding of the androgen receptor (AR) associated with increased chromatin accessibility, Histone 3 lysine 27 acetylation (H3K27ac), and ~3-fold upregulation of IRX4 expression. We further show that an MNLP is amongst the strongest candidate susceptibility variants at two additional prostate cancer risk loci. We estimated that at least 5% of prostate cancer risk loci could be explained by functional non-SNP causal variants, which may have broader implications for other cancers GWAS. More generally, our results underscore the importance of investigating other classes of inherited variation as causal mediators of human traits.

ATM deficiency confers specific therapeutic vulnerabilities in bladder cancer.

Ataxia-telangiectasia mutated (ATM) plays a central role in the cellular response to DNA damage and ATM alterations are common in several tumor types including bladder cancer. However, the specific impact of ATM alterations on therapy response in bladder cancer is uncertain. Here, we combine preclinical modeling and clinical analyses to comprehensively define the impact of ATM alterations on bladder cancer. We show that ATM loss is sufficient to increase sensitivity to DNA-damaging agents including cisplatin and radiation. Furthermore, ATM loss drives sensitivity to DNA repair-targeted agents including poly(ADP-ribose) polymerase (PARP) and Ataxia telangiectasia and Rad3 related (ATR) inhibitors. ATM loss alters the immune microenvironment and improves anti-PD1 response in preclinical bladder models but is not associated with improved anti-PD1/PD-L1 response in clinical cohorts. Last, we show that ATM expression by immunohistochemistry is strongly correlated with response to chemoradiotherapy. Together, these data define a potential role for ATM as a predictive biomarker in bladder cancer.

ALK peptide vaccination restores the immunogenicity of ALK-rearranged non-small cell lung cancer.

Anaplastic lymphoma kinase (ALK)-rearranged non-small cell lung cancer (NSCLC) is treated with ALK tyrosine kinase inhibitors (TKIs), but the lack of activity of immune checkpoint inhibitors (ICIs) is poorly understood. Here, we identified immunogenic ALK peptides to show that ICIs induced rejection of ALK+ tumors in the flank but not in the lung. A single-peptide vaccination restored priming of ALK-specific CD8+ T cells, eradicated lung tumors in combination with ALK TKIs and prevented metastatic dissemination of tumors to the brain. The poor response of ALK+ NSCLC to ICIs was due to ineffective CD8+ T cell priming against ALK antigens and is circumvented through specific vaccination. Finally, we identified human ALK peptides displayed by HLA-A*02:01 and HLA-B*07:02 molecules. These peptides were immunogenic in HLA-transgenic mice and were recognized by CD8+ T cells from individuals with NSCLC, paving the way for the development of a clinical vaccine to treat ALK+ NSCLC.