EPHA5 mutation predicts the durable clinical benefit of immune checkpoint inhibitors in patients with lung adenocarcinoma

doi:10.21203/rs.3.rs-27361/v1

Download PDF

Primary research

EPHA5 mutation predicts the durable clinical benefit of immune checkpoint inhibitors in patients with lung adenocarcinoma

https://doi.org/10.21203/rs.3.rs-27361/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 06 Aug, 2020

Read the published version in Cancer Gene Therapy →

Version 1

posted

You are reading this latest preprint version

Background: Immune checkpoint inhibitor (ICI ) therapy has shown remarkable clinical benefit in lung adenocarcinoma (LUAD) patients. Genomic mutations may be applicable to predict the response to ICIs. Eph receptor A5 (EPHA5) is frequently mutated in breast cancer, lung cancer and other tumors; however, its association with outcome in patients who receive immunotherapy remains unknown.

Methods: Somatic mutation data, gene expression data and clinicopathologic information were curated from patients with LUAD who were treated with ICI therapy (MSKCC ICI-treated cohort) and patients with LUAD who did not receive ICI therapy (TCGA-LUAD cohort). The CIBERSORT and gene set enrichment analysis (GSEA) algorithms were separately used to infer the relative abundance of leukocytes and significantly enriched pathways in the defined subgroups. Tumor mutation burden (TMB), tumor neoantigen load, the fraction of copy number variants (CNVs), the mutation frequency of DNA damage repair (DDR) genes, and the prognosis of patients were compared between the EPHA5-mutated subgroup and the EPHA5 wild - type subgroup. Whole-exome sequencing (WES) data and drug sensitivity data downloaded from the Genomics of Drug Sensitivity in Cancer (GDSC) database were used to evaluate the relationship of the 50% inhibitory concentration (IC50) of multiple drugs and the EPHA5 mutation status in LUAD cells.

Results: EPHA5 mutations were associated with increased TMB, neoantigen load, levels of immune - related gene expression signatures, and enhanced tumor-infiltrating lymphocytes (TILs). Patients with EPHA5 mutations in the immunotherapy cohort achieved a longer progression - free survival (PFS) time than patients with wild-type EPHA5. Immune response pathways were among the top enriched pathways in samples with EPHA5 mutations. In addition , patients with EPHA5 mutations tended to be more sensitive to certain targeted molecular inhibitors , including serdemetan, lox2 and PF1-1.

Conclusion: Our results suggest that identifying mutations in the EPHA5 gene may provide insight into the genome-wide mutational burden and may serve as a biomarker to predict the immune response of patients with LUAD.

Cancer Biology

Lung adenocarcinoma

EPHA5

Mutation

Immune checkpoint inhibitor

Tumor mutation burden

Neoantigen load

Currently, lung cancer is still the leading cause of cancer-related deaths worldwide [1]. The majority of lung cancer patients are non-small-cell lung cancer (NSCLC) and are in an advanced stage [2]. NSCLC can be histologically classified into lung adenocarcinoma (LUAD) and squamous cell carcinoma, with LUAD being the most common subtype, comprising more than 50% of all NSCLCs [3]. At present, a therapeutic plateau has been reached, with a response rate of only approximately 20% in LUAD patients who receive traditional treatment [4]. Therefore, new selective treatments are urgently needed to improve the outcome of patients with advanced LUAD.

Immune checkpoint inhibitors (ICIs), including anti-PD-(L)1 monotherapy and anti-cytotoxic T cell lymphocyte-4 (anti-CTLA-4), have emerged as therapeutic landscape for patients with advanced cancers, including LUAD [5–8]. However, only 20%-30% of LUAD patients respond to ICIs, with most patients not achieving objective responses [9]. The mechanisms underlying variations in the response to ICIs among LUAD patients are still unclear, and potential biomarkers that are predictive of the benefit of ICIs are needed to elucidate the mechanisms of resistance in ICI nonresponders.

To date, the only Food and Drug Administration-approved predictive biomarkers of the response to ICIs in NSCLC patients are programmed death-ligand 1 (PD-L1) expression, the microsatellite instability (MSI) status, and mismatch repair deficiency [10, 11], but the sensitivity and specificity of these approaches are still modest [12]. Tumor mutation burden (TMB) [13], the intensity of CD8 + T cell infiltrates [14], copy number variants (CNVs), and the predicted neoantigen load [15] have also been proposed as distinct biomarkers of the response to ICI therapies in NSCLC patients. Moreover, many studies have found that mutations in genes associated with genomic integrity, such as EGFR and P53 [16], could cause replication stress and genomic instability, resulting in a high TMB and neoantigen burden in NSCLC patients. This raises the question of whether other mutations simultaneously affect the above factors to provide stronger predictive value for therapeutic outcomes in LUAD patients.

The Eph receptor tyrosine kinase (RTK) family is the largest group of surface receptors and is classified into the EphA or EphB subfamily [17]. Eph receptors have been proven to play an important role in tumor immunity [18]. For example, CD4 + T cells recognize the EphA3 epitope and then elicit selective immunoreactivity against melanoma cells [19]. In addition, CD4⁺ and CD8⁺ T cells recognize multiple epitopes of EPHA2 and enhance the immune activity of CD8⁺ T cells in cancers [20, 21]. In addition to inducing immunologic responses, Eph receptors also regulate immune cell trafficking. EphA1 and EphA4 are expressed in T cells and promote the chemotaxis and migration of CD4⁺ and CD8⁺ T cells by recruiting multiple kinase pathways [22–24]. Furthermore, a recent study showed that mutations in the EPHA5 gene impair natural killer (NK) cell-mediated cytotoxicity against NSCLC cells [25], suggesting that Eph receptors affect the tumor immune microenvironment. Thus, evaluating the roles of the Eph receptor in the response to ICIs, especially those that are driven by genetic changes, may provide additional insight into the tumor immune microenvironment.

In the present study, we first examined EPHA5, a member of the Eph/ephrin family, as a new biomarker for the response to ICIs in patients with LUAD. By comparing the gene mutation profiles of ICI responders and non-ICI responders based on the ICI-treated LUAD cohort reported by Rizvi [13], we identified a high mutation frequency of EPHA5 in ICI-responsive patients. We then analyzed the relationship of EPHA5 mutation and the prognosis of LUAD patients in both the ICI-treated LUAD cohort and the non-ICI-treated LUAD cohort. We also investigated whether EPHA5 mutation is associated with the immune microenvironment in LUAD. Our current evidence suggests that EPHA5 alterations yield promising predictive value for ICI treatment in LUAD patients.

Public datasets used in this study

An ICI-treated cohort comprising 240 NSCLC patients with a detailed clinical annotation, response data and mutation data (mutation data were obtained from targeted next-generation sequencing (NGS)) described in the cBioPortal for Cancer Genomics at http://www.cbioportal.org (240MSKCC ICI-treated cohort) was analyzed. All patients were treated with anti-PD-(L)1 monotherapy alone or anti-PD-(L)1 monotherapy combined with anti-CTLA-4. The efficacy of tumor immunotherapy was evaluated using the Response Evaluation Criteria in Solid Tumors (RECIST) version 1.1 and defined as durable clinical benefit (DCB) or no durable benefit (NDB) [13].

Somatic mutation data, gene expression data (Illumina HiSeq, RNA-Seq) and overall survival (OS) data from the TCGA-LUAD cohort were downloaded from the Genomic Data Commons (https://portal.gdc.cancer.gov/) by using the R package TCGAbiolinks [26], while disease-free survival (DFS) data were downloaded from cBioPortal [27].

The GDSC-LUAD dataset contains somatic mutation data and drug sensitivity data on multiple LUAD cell lines and was downloaded from the Genomics of Drug Sensitivity in Cancer (GDSC) database (https://www.cancerrxgene.org) [28].

Comparison of the prognosis of LUAD patients with mutant EPHA5 and those with wild-type EPHA5

A subset of ICI-treated LUAD patients (n = 187, from the 240MSKCC ICI-treated cohort) with annotated clinical records was divided into two subgroups, EPHA5-MT (mutant type) and EPHA5-WT (wild-type), and PFS, defined as the date the patient began immunotherapy to the date of progression, was assessed by Kaplan-Meier (KM) analysis and compared between the two subgroups. We also calculated the DFS and OS of 494 non-ICI-treated patients from the TCGA-LUAD cohort depending on the status of EPHA5 mutations by KM analysis.

EPHA5 Mutation And Tumor Immunogenicity

The somatic mutation data from the ICI-treated LUAD cohort were analyzed by targeted NGS (MSK-IMPACT) as described by Hira Rizvi. In the TCGA-LUAD and GDSC-LUAD cohorts, TMB was determined by using the nonsynonymous mutations as the raw mutation count and dividing by 38 Mb [29], while the TMB score was then applied to the ICI-treated LUAD cohort. We used the ComplexHeatmap [30] package in R to visualize the genetic mutations in the top 20 genes with the highest alteration frequencies and clinical characteristics of patients in the ICI-treated LUAD and TCGA-LUAD cohorts. The maftools [31] package in R was used to visualize EPHA5 mutations in the ICI-treated LUAD and TCGA-LUAD cohorts.

Tumor-infiltrating lymphocyte (TIL) analysis and T cell-inflamed gene expression (RNA) profiling

The CIBERSORT algorithm [32] (http://cibersort.stanford.edu/) was used to compare the abundance of 22 TILs between the EPHA5-MT and EPHA5-WT subgroups in the TCGA-LUAD cohort. In addition, differences in the mRNA expression of T cell-inflamed genes associated with antigen presentation, adaptive immune resistance, chemokine expression, and cytolytic activity between the EPHA5-WT and EPHA5-WT subgroups in the TCGA-LUAD cohort were also compared [33]. The differences in mRNA expression were calculated as follows: p-value < 0.05 and log₂(fold change) > 0.58 or < (-0.58).

Drug Sensitivity Analysis Of The EPHA5-WT And EPHA5-MT Subgroups

According to the mutation status of EPHA5, the cell lines corresponding to code "LUAD" in the GDSC database were divided into two groups, and then the two groups of cells were compared for their sensitivity to compounds including chemotherapy and molecular-targeted drugs.

CNV Analysis

The Affymetrix SNP 6.0 microarray data of the TCGA-LUAD cohorts were downloaded from the Broad GDAC Firehose (http://gdac.broadinstitute.org/). We then used GenePattern [34] (https://cloud.genepattern.org/gp/pages/index.jsf) to analyze the downloaded CNV segments with GISTIC 2.0. We set the confidence level to 0.99 and did not exclude the X chromosome before analysis. All other GISTIC 2.0 parameters were default values.

Pathway Enrichment Analysis

The gene expression profiles of the TCGA-LUAD cohort were downloaded and partitioned into two groups according to the mutation status of EPHA5. The differentially expressed genes between the two groups were obtained by using the R package edgeR [35] and then fed into the R-function in the clusterProfiler [36] package for gene set enrichment analysis (GSEA). A gene set with a nominal p-value less than 0.05 was considered significantly enriched in Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) and Reactome pathways.

EPHA5 Mutation And DNA Damage Repair (DDR) Gene Mutation

The gene sets associated with the DDR pathways were obtained from the Molecular Signatures Database (MSigDB) of the Broad Institute [37]. The numbers of nonsynonymous mutations in DDR pathway-related genes were compared between the EPHA5-WT subgroup and the EPHA5-MT subgroup in the ICI-treated LUAD cohort, TCGA-LUAD cohort, and GDSC-LUAD cohort.

Statistical analysis

The comparisons of TMB, neoantigen load, immune cell abundance, immune-related gene expression, age, pack years, MSI score and number of mutations in DDR pathway genes between the EPHA5-WT and EPHA5-MT subgroups were performed using the Mann-Whitney U test. Differences in the distributions of baseline characteristics, including gene mutation status, sex, smoking history, treatment stage, and response, between subgroups were investigated using the chi-square test or Fisher’s exact test. Clinical outcomes, including OS, PFS, and DFS, were calculated using KM analysis and compared using the log-rank test. All p values are two-sided, with p ≤ 0.05 taken as significant. All statistical analyses were performed using R software ver. 3.6.1 (R Foundation, Vienna, Austria). The ggpubr package (https://CRAN.R-project.org/package=ggpubr) was used to generate boxplots for visualization. The maftools package [31] was used to visualize CNV, and the false discovery rate (FDR) was set to 0.05.

Association of the EPHA5 mutation status with clinical characteristics of patients with LUAD

To identify the mutated genes that might be related to the response to ICIs in advanced LUAD patients, we first investigated genes with different mutation frequencies between patients who achieved a DCB and those who achieved NDB among all patients who received ICIs in the ICI-treated LUAD cohort. The gene mutation pattern showed that EPHA5 was one of the top 20 frequently mutated genes in patients who achieved DCB compared to patients who achieved NDB. Since EPHA5 encodes an RTK involved in cancer development, and EPHA5 was found to be mutated in multiple cancers. We thus speculated that EPHA5 may be related to the response to ICIs. To verify our hypothesis, we initially analyzed the mutation characteristics of EPHA5 based on both the ICI-treated LUAD and TCGA-LUAD datasets. The results revealed that the mutation sites of EPHA5 spanned the entire gene, including the ligand binding site and receptor activity sites, which might affect the protein structure and function of EPHA5 (Fig. 1a). We then divided the two cohorts into the EPHA5-WT and EPHA5-MT subgroups and analyzed the correlation between the EPHA5 mutation status and clinical characteristics of patients. As shown in Fig. 1b, c, analysis of the two datasets showed that a history of tobacco exposure exhibited the strongest association with the EPHA5 mutation status among all clinical variables considered. Interestingly, EPHA5 mutation was also associated with TMB and neoantigen load. EPHA5 mutation also showed a positive correlation with the mutation frequency of several common mutations found in LUAD, such as missense mutations in TP53, KEAP1 and LRP1B, which affect the efficacy of immunotherapies, suggesting that EPHA5 mutations may be related to the response to ICIs in patients with LUAD.

EPHA5 mutation is associated with the clinical benefit of ICIs in LUAD patients

Since TMB and neoantigen load serve as predictive markers of the response to ICIs, we first aimed to determine whether EPHA5 mutation in LUAD patients affects the whole TMB and alters tumor antigenicity. By analyzing the ICI-treated LUAD and non-ICI-treated TCGA-LUAD datasets, we determined that the EPHA5-MT subgroups had a higher TMB and neoantigen load (Fig. 2a), suggesting a potential relationship between mutated EPHA5 and the response to ICIs in LUAD patients. We then assessed the association of EPHA5 mutation with the prognosis of LUAD patients. The results showed that patients in the EPHA5-MT subgroup experienced longer PFS than those in the EPHA5-WT subgroup when analyzing the ICI-treated LUAD cohorts (p = 0.039) (Fig. 2b); however, there was no significant difference in either DFS or OS between the EPHA5-WT subgroup and the EPHA5-MT subgroup in the non-ICl-treated TCGA-LUAD cohort (p>0.05) (Fig. 2c). This finding suggests that LUAD patients with EPHA5 mutations benefit from ICI treatment but not non-ICI treatment.

EPHA5 Mutation Is Associated With High CNV Counts

The acquisition of tumor aneuploidy correlates with markers of immune evasion and is associated with reduced immune-mediated cytotoxic function, and high CNV counts have been confirmed to be associated with poor clinical outcomes regarding NSCLC immunotherapy. Therefore, we explored the associations of the EPHA5 mutation status with CNVs at the chromosome arm level in the TCGA-LUAD cohorts. As shown in Fig. 3, we found that CNV counts were relatively lower in the EPHA5-MT subgroup than in the EPHA5-WT subgroup at 1q21.2, 1q22, 8q24.21, 9p21.3, and 14q13.3. However, CNV counts could not be validated in the ICI-treated LUAD dataset due to limited genomic coverage of the targeted panel.

EPHA5 Mutation Is Associated With Activated Antitumor Immunity

TILs are important markers for predicting the response to ICI treatment. We thus continued to survey the relationships between EPHA5 mutation and common immune infiltrates in the LUAD-TCGA cohort. Our results showed that tumor-infiltrating T lymphocytes, including CD8 + T cells, CD4 + resting memory T cells and macrophage M1 cells, were generally more abundant in the EPHA5-MT tumors than in the EPHA5-WT tumors, as evaluated by the expression levels of immune cell marker genes and CIBERSORT analysis (Fig. 4a, b). The immune signature analysis also revealed significantly increased expression levels of chemokines and cytolytic activity-associated gene signatures (CCL5, CXCL9, CD8A, CD8B, and GZMB) in EPHA5 mutant tumors compared to wild-type EPHA5 tumors. We also found that the levels of LAG3, an immune checkpoint marker, were higher in the EPHA5 mutation subgroups (Fig. 4c).

These results indicate that EPHA5 mutation is strongly associated with a relatively hot immune microenvironment, which firmly supports its ability to predict the response to ICIs in LUAD patients.

EPHA5 mutation is associated with high mutation frequencies of DDR pathway genes

Considering that DDR gene alterations are associated with the response to ICIs in tumors, we further examined the association of the EPHA5 mutation status and DDR gene mutations in both the ICI-treated LUAD cohort and the non-ICI-treated TCGA-LUAD cohort. We observed a trend towards an enrichment of DDR gene alterations in the EPHA5 mutation subgroups. Although this trend was not observed in the base excision repair (BER) pathway in the ICI-treated LUAD cohort, eight other pathways (i.e., homologous recombination (HR), mismatch repair (MMR), Fanconi anemia (FA), nucleotide excision repair (NER), nonhomologous end joining (NHEJ), DDR, DNA double-strand breaks (DSBs) and single-strand breaks (SSBs)) demonstrated enrichment for mutations in the EPHA5-MT subgroups versus the EPHA5-WT subgroups (Fig. 5).

EPHA5 Mutation And Drug Selection Of LUAD Cells

In addition to examining the role of EPHA5 mutation in disease progression, we also investigated the role of EPHA5 mutation in chemotherapy and molecular-targeted drug therapy in LUAD patients. By analyzing the GDSC database, we found that the mutation status of EPHA5 had no effect on cell sensitivity to common chemotherapy drugs, such as cisplatin and irinotecan (Fig. 6a). Interestingly, the EPHA5-WT subgroup and the EPHA5-MT subgroup had different sensitivities to 37 molecular-targeted drugs, of which serdemetan, lox2 and PF1-1 were relatively sensitive in the EPHA5-MT subgroup, while the other 34 drugs were more resistant in the EPHA5-MT subgroup (Fig. 6b and Fig. S1). This result may provide a reference for the treatment choice in LUAD patients with EPHA5 mutations.

Pathway Enrichment Analysis Of EPHA5 Mutation

Finally, we investigated whether EPHA5 mutation is involved in immune-related pathways using the LUAD-TCGA cohorts. GSEA showed that the insulin receptor, negative regulation of Notch signaling, MHC II antigen presentation, regulation of p53, and NK cell-mediated cytotoxicity pathways, which enhance antitumor immune activity in cancers, were significantly enriched in with the EPHA5-MT subgroup, while the positive regulation of the ERK1 and ERK2 cascade, PI3K cascade, FGFR2 function, and NF-kappaB signaling pathways were significantly enriched in the EPHA5-WT subgroup. In addition, consistent with previous studies showing that microRNAs (miRNAs) may play a role in regulating tumor immunity through the direct regulation of genes involved in immune activation [38], our results revealed a significant enrichment of gene silencing by miRNAs in the EPHA5-MT subgroup (Fig. 7). These results suggest that EPHA5 mutation might affect tumor immunity in LUAD patients by influencing multiple pathways.

In this study, we observed that EPHA5 mutation was enriched in patients who respond to ICIs and is strongly predictive of DCB among LUAD patients treated with ICIs. Based on the analysis of the ICI-treated LUAD dataset and the non-ICI-treated TCGA-LUAD dataset, as well as the GDSC-LUAD dataset, EPHA5 mutation is significantly associated with enhanced tumor immunogenicity, activated antitumor immunity, a high mutation burden of DDR genes, and low CNV counts at multiple chromosomal arms. We also determined that LUAD patients with EPHA5 mutations are more sensitive to certain molecular-targeted drugs, such as serdemetan, lox2 and PF1-1. Furthermore, samples with EPHA5 mutations were characterized by the upregulation of signaling pathways involved in a hot immune microenvironment, such as MHC II antigen presentation and NK cell-mediated cytotoxicity. These results suggest that EPHA5 mutation may be a predictor of the immune response and promote the immune response in patients with LUAD.

The Eph/ephrin family is the largest known family of RTKs in mammals [39]. The Eph/ephrin family of receptors consists of 9 EphA members and 5 EphB members classified according to sequence homology [40]. There is evidence that Eph receptors may mediate immune cell activation [18, 41], but data regarding the potential correlation between genomic alterations of Eph/ephrin family genes and the response to ICIs in LUAD patients are not available. Our present study represents the first report on the association between the response to ICIs and specific mutated genes belonging to the Eph/ephrin family in LUAD patients. Among the multiple tumor-related genes identified, we found that EPHA5 mutation was strongly associated with prolonged PFS in patients who received ICI treatment. These findings add great value to the robust link between the Eph/ephrin family and immunotherapy.

Previous studies have observed that EPHA5 is frequently mutated in several cancers, including NSCLC [25, 42, 43]. Mutations in the EPHA5 gene could impair NK cell-mediated cytotoxicity against cancer cells [25] and predict the overall prognosis of cancer patients [44]. In addition to the effects on gene mutations, changes in EPHA5 transcription levels also affect tumor progression (e.g., enhanced drug resistance, invasion and metastatic ability of cancer cells) [45–47]. In our present study, we focused on the role of genomic alterations of EPHA5 in the response to ICIs. We examined the clinical relevance of EPHA5 alterations in the MSKCC and TCGA databases. We observed that mutated EPHA5 was beneficial for LUAD patients who received ICIs. However, for patients who did not receive ICIs, mutated EPHA5 was not significantly associated with prognosis, suggesting that EPHA5 mutation may be related to the tumor immune environment. The results of the subsequent analysis were consistent with our expectations. EPHA5 mutations are associated with a high tumor TMB, neoantigen load, T lymphocyte infiltration, and expression of immune-related signatures. EPHA5 mutations are also related to CNV in certain chromosome arms and mutations in genes related to the DDR signaling pathways. Through GSEA, we further revealed the potential mechanisms by which mutant EPHA5 regulates the immune response in LUAD. This evidence suggests that mutations in a single gene may be used to predict the response to ICIs in LUAD patients.

This retrospective analysis has several limitations. First, the relatively small number of patients treated with ICIs limits the statistical power of our analysis. Second, the ICI-treated cohort includes patients who received anti-PD-1 as well as anti-PD-L1 [anti-PD-(L)1]-based combination therapy. A difference in treatment may lead to a difference in response. Third, since a single targeted NGS panel was used in the ICI-treated LUAD dataset and gene expression data are lacking, this dataset could not be used to assess the correlation of EPHA5 mutations with immune cell infiltration, CNV counts and related pathways. However, we assessed the correlation of EPHA5 mutations with immune characteristics in LUAD by analyzing the non-ICI-treated TCGA-LUAD dataset. Although our study provides some interesting findings, it should still be tested in a larger and more uniformly treated cohort.

Our study provides solid evidence that EPHA5 mutation is associated with a high objective response rate, good DCB, and prolonged PFS in patients receiving ICIs. Therefore, EPHA5 mutation can serve as a novel predictive biomarker for the response to ICIs in LUAD patients. Further exploration of molecular mechanisms and prospective clinical trials are warranted.

ICI：Immune checkpoint inhibitor; LUAD：lung adenocarcinoma; EPHA5: Eph receptor A5; GSEA: gene set enrichment analysis; TMB: Tumor mutation burden; CNV: copy number variant; DDR: DNA damage repair; WES: Whole-exome sequencing; GDSC: Drug Sensitivity in Cancer; IC50: 50% inhibitory concentration; TILs: tumor-infiltrating lymphocytes; PFS: progression-free survival; NSCLC: non-small-cell lung cancer; anti-CTLA-4: anti-cytotoxic T cell lymphocyte-4; PD-L1: programmed death-ligand 1; RTK: Eph receptor tyrosine kinase; NK: natural killer; NGS: targeted next-generation sequencing; DCB: durable clinical benefit; NDB: no durable benefit; OS: overall survival; DFS: disease-free survival; GDSC: Genomics of Drug Sensitivity in Cancer; KM: Kaplan-Meier; GO: Gene Ontology; KEGG: Kyoto Encyclopedia of Genes and Genomes; TCGA: The Cancer Genome Atlas; BER: base excision repair; HR: homologous recombination; MMR: mismatch repair; FA: Fanconi anemia; NER: nucleotide excision repair; NHEJ: nonhomologous end joining; DSBs: DNA double-strand breaks; SSBs: single-strand breaks; miRNAs: microRNAs;

Acknowledgements

Not applicable.

Authors’ contributions

Conceptualization and methodology: AL and PL; Formal analysis and data curation: WH and WZ; Validation and investigation: TW, LG and JZ; Writing—original draft preparation and Writing—review and editing: WH, AL and PL; Approval of final manuscript: all authors. All authors read and approved the final manuscript.

Funding

No funding was received.

Availability of data and materials

The analyzed data sets generated during the present study are available from the corresponding author on reasonable request.

Ethics approval and consent to participate

Not applicable.

Patient consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424.
Kumarakulasinghe NB, van Zanwijk N, Soo RA. Molecular targeted therapy in the treatment of advanced stage non-small cell lung cancer (NSCLC). Respirology. 2015;20(3):370–8.
Travis W, Brambilla E, Noguchi M, Nicholson A, Geisinger K, Yatabe Y, et al. International association for the study of lung cancer/american thoracic society/european respiratory society international multidisciplinary classification of lung adenocarcinoma. Journal of thoracic oncology: official publication of the International Association for the Study of Lung Cancer. 2011;6(2):244–85.
Liu K, Guo J, Liu K, Fan P, Zeng Y, Xu C, et al. Integrative analysis reveals distinct subtypes with therapeutic implications in KRAS-mutant lung adenocarcinoma. EBioMedicine. 2018;36:196–208.
Hugo W, Zaretsky JM, Sun L, Song C, Moreno BH, Hu-Lieskovan S, et al. Genomic and Transcriptomic Features of Response to Anti-PD-1 Therapy in Metastatic Melanoma. Cell. 2017;168(3):542.
Sarfaty M, Leshno M, Gordon N, Moore A, Neiman V, Rosenbaum E, et al. Cost Effectiveness of Nivolumab in Advanced Renal Cell Carcinoma. Eur Urol. 2018;73(4):628–34.
Seiwert TY, Burtness B, Mehra R, Weiss J, Berger R, Eder JP, et al. Safety and clinical activity of pembrolizumab for treatment of recurrent or metastatic squamous cell carcinoma of the head and neck (KEYNOTE-012): an open-label, multicentre, phase 1b trial. Lancet Oncol. 2016;17(7):956–65.
Assi HI, Kamphorst AO, Moukalled NM, Ramalingam SS. Immune checkpoint inhibitors in advanced non-small cell lung cancer. Cancer. 2018;124(2):248–61.
Carbognin L, Pilotto S, Milella M, Vaccaro V, Brunelli M, Calio A, et al. Differential Activity of Nivolumab, Pembrolizumab and MPDL3280A according to the Tumor Expression of Programmed Death-Ligand-1 (PD-L1): Sensitivity Analysis of Trials in Melanoma, Lung and Genitourinary Cancers. PLoS One. 2015;10(6):e0130142.
Reck M, Rodriguez-Abreu D, Robinson AG, Hui R, Csoszi T, Fulop A, et al. Pembrolizumab versus Chemotherapy for PD-L1-Positive Non-Small-Cell Lung Cancer. N Engl J Med. 2016;375(19):1823–33.
Le DT, Durham JN, Smith KN, Wang H, Bartlett BR, Aulakh LK, et al. Mismatch repair deficiency predicts response of solid tumors to PD-1 blockade. Science. 2017;357(6349):409–13.
Taube JM, Klein A, Brahmer JR, Xu H, Pan X, Kim JH, et al. Association of PD-1, PD-1 ligands, and other features of the tumor immune microenvironment with response to anti-PD-1 therapy. Clin Cancer Res. 2014;20(19):5064–74.
Rizvi H, Sanchez-Vega F, La K, Chatila W, Jonsson P, Halpenny D, et al. Molecular Determinants of Response to Anti-Programmed Cell Death (PD)-1 and Anti-Programmed Death-Ligand 1 (PD-L1) Blockade in Patients With Non-Small-Cell Lung Cancer Profiled With Targeted Next-Generation Sequencing. J Clin Oncol. 2018;36(7):633–41.
Tang H, Wang Y, Chlewicki LK, Zhang Y, Guo J, Liang W, et al. Facilitating T Cell Infiltration in Tumor Microenvironment Overcomes Resistance to PD-L1 Blockade. Cancer Cell. 2016;29(3):285–96.
Lauss M, Donia M, Harbst K, Andersen R, Mitra S, Rosengren F, et al. Mutational and putative neoantigen load predict clinical benefit of adoptive T cell therapy in melanoma. Nat Commun. 2017;8(1):1738.
Dong ZY, Zhong WZ, Zhang XC, Su J, Xie Z, Liu SY, et al. Potential Predictive Value of TP53 and KRAS Mutation Status for Response to PD-1 Blockade Immunotherapy in Lung Adenocarcinoma. Clin Cancer Res. 2017;23(12):3012–24.
Gale NW, Holland SJ, Valenzuela DM, Flenniken A, Pan L, Ryan TE, et al. Eph receptors and ligands comprise two major specificity subclasses and are reciprocally compartmentalized during embryogenesis. Neuron. 1996;17(1):9–19.
Shiuan E, Chen J. Eph Receptor Tyrosine Kinases in Tumor Immunity. Cancer Res. 2016;76(22):6452–7.
Chiari R, Hames G, Stroobant V, Texier C, Maillere B, Boon T, et al. Identification of a tumor-specific shared antigen derived from an Eph receptor and presented to CD4 T cells on HLA class II molecules. Cancer Res. 2000;60(17):4855–63.
Tatsumi T, Herrem CJ, Olson WC, Finke JH, Bukowski RM, Kinch MS, et al. Disease stage variation in CD4 + and CD8 + T-cell reactivity to the receptor tyrosine kinase EphA2 in patients with renal cell carcinoma. Cancer Res. 2003;63(15):4481–9.
Alves PM, Faure O, Graff-Dubois S, Gross DA, Cornet S, Chouaib S, et al. EphA2 as target of anticancer immunotherapy: identification of HLA-A*0201-restricted epitopes. Cancer Res. 2003;63(23):8476–80.
Aasheim HC, Delabie J, Finne EF. Ephrin-A1 binding to CD4 + T lymphocytes stimulates migration and induces tyrosine phosphorylation of PYK2. Blood. 2005;105(7):2869–76.
Hjorthaug HS, Aasheim HC. Ephrin-A1 stimulates migration of CD8 + CCR7 + T lymphocytes. Eur J Immunol. 2007;37(8):2326–36.
Holen HL, Nustad K, Aasheim HC. Activation of EphA receptors on CD4 + CD45RO + memory cells stimulates migration. J Leukoc Biol. 2010;87(6):1059–68.
Zhang J, Zhang Z, Song W, Liu J, Zhang Y. EPHA5 mutation impairs natural killer cell-mediated cytotoxicity against non-small lung cancer cells and promotes cancer cell migration and invasion. Mol Cell Probes. 2020:101566.
Colaprico A, Silva TC, Olsen C, Garofano L, Cava C, Garolini D, et al. TCGAbiolinks: an R/Bioconductor package for integrative analysis of TCGA data. Nucleic Acids Res. 2016;44(8):e71.
Cerami E, Gao J, Dogrusoz U, Gross BE, Sumer SO, Aksoy BA, et al. The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer Discov. 2012;2(5):401–4.
Yang W, Soares J, Greninger P, Edelman E, Lightfoot H, Forbes S, et al. Genomics of Drug Sensitivity in Cancer (GDSC): a resource for therapeutic biomarker discovery in cancer cells. Nucleic acids research. 2013;41:D955-61.
Chalmers ZR, Connelly CF, Fabrizio D, Gay L, Ali SM, Ennis R, et al. Analysis of 100,000 human cancer genomes reveals the landscape of tumor mutational burden. Genome Med. 2017;9(1):34.
Gu Z, Eils R, Schlesner M. Complex heatmaps reveal patterns and correlations in multidimensional genomic data. Bioinformatics. 2016;32(18):2847–9.
Mayakonda A, Lin DC, Assenov Y, Plass C, Koeffler HP. Maftools: efficient and comprehensive analysis of somatic variants in cancer. Genome Res. 2018;28(11):1747–56.
Newman AM, Liu CL, Green MR, Gentles AJ, Feng W, Xu Y, et al. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. 2015;12(5):453–7.
Thorsson V, Gibbs DL, Brown SD, Wolf D, Bortone DS, Ou Yang TH, et al. The Immune Landscape of Cancer. Immunity. 2018;48(4):812 – 30.e14.
Reich M, Liefeld T, Gould J, Lerner J, Tamayo P, Mesirov J. GenePattern 2.0. Nature genetics. 2006; 38(5):500–1.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.
Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. Omics. 2012;16(5):284–7.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102(43):15545–50.
Vishnubalaji R, Shaath H, Elango R, Alajez NM. Noncoding RNAs as potential mediators of resistance to cancer immunotherapy. Semin Cancer Biol. 2019.
Tuzi NL, Gullick WJ. eph, the largest known family of putative growth factor receptors. Br J Cancer. 1994;69(3):417–21.
Pasquale EB. The Eph family of receptors. Curr Opin Cell Biol. 1997;9(5):608–15.
Darling T, Lamb T. Emerging Roles for Eph Receptors and Ephrin Ligands in Immunity. Frontiers in immunology. 2019;10:1473.
Jeong HM, Kim RN, Kwon MJ, Oh E, Han J, Lee SK, et al. Targeted exome sequencing of Korean triple-negative breast cancer reveals homozygous deletions associated with poor prognosis of adjuvant chemotherapy-treated patients. Oncotarget. 2017;8(37):61538–50.
Quinn AM, Hickson N, Adaway M, Priest L, Jaeger E, Udar N, et al. Diagnostic Mutation Profiling and Validation of Non-Small-Cell Lung Cancer Small Biopsy Samples using a High Throughput Platform. J Thorac Oncol. 2015;10(5):784–92.
Qian J, Nie W, Lu J, Zhang L, Zhang Y, Zhang B, et al. Racial differences in characteristics and prognoses between Asian and white patients with nonsmall cell lung cancer receiving atezolizumab: An ancillary analysis of the POPLAR and OAK studies. Int J Cancer. 2019.
Li Y, Chu J, Feng W, Yang M, Zhang Y, Zhang Y, et al. EPHA5 mediates trastuzumab resistance in HER2-positive breast cancers through regulating cancer stem cell-like properties. FASEB journal: official publication of the Federation of American Societies for Experimental Biology. 2019;33(4):4851–65.
Gu S, Feng J, Jin Q, Wang W, Zhang S. Reduced expression of EphA5 is associated with lymph node metastasis, advanced TNM stage, and poor prognosis in colorectal carcinoma. Histology histopathology. 2017;32(5):491–7.
Staquicini F, Qian M, Salameh A, Dobroff A, Edwards J, Cimino D, et al. Receptor tyrosine kinase EphA5 is a functional molecular target in human lung cancer. J Biol Chem. 2015;290(12):7345–59.

Figure S1. EPHA5 mutation and drug selection of LUAD cells. Boxplot showing molecular-targeted drugs with significant 50% inhibitory concentration (IC50) values in the EPHA5-WT and EPHA5-MT subgroups.

FigureS1.tif

Download PDF

Journal Publication

published 06 Aug, 2020

Read the published version in Cancer Gene Therapy →

Version 1

posted

You are reading this latest preprint version

EPHA5 mutation predicts the durable clinical benefit of immune checkpoint inhibitors in patients with lung adenocarcinoma

Status:

Journal Publication

Version 1

Abstract

Figures

Background

Methods

Public datasets used in this study

EPHA5 Mutation And Tumor Immunogenicity

Drug Sensitivity Analysis Of The EPHA5-WT And EPHA5-MT Subgroups

CNV Analysis

Pathway Enrichment Analysis

EPHA5 Mutation And DNA Damage Repair (DDR) Gene Mutation

Statistical analysis

Results

Association of the EPHA5 mutation status with clinical characteristics of patients with LUAD

EPHA5 Mutation Is Associated With High CNV Counts

EPHA5 Mutation Is Associated With Activated Antitumor Immunity

EPHA5 Mutation And Drug Selection Of LUAD Cells

Pathway Enrichment Analysis Of EPHA5 Mutation

Discussion

Conclusion

Abbreviations

Declarations

References

Supplementary Figure Caption

Supplementary Files

Status:

Journal Publication

Version 1