Association of TP53 Mutation Status and Sex with Clinical Outcome in Non–Small Cell Lung Cancer Treated with Immune Checkpoint Inhibitors: A Retrospective Cohort Study
Article information
Abstract
Purpose
Some studies suggest that TP53 mutations are associated with the response to immune checkpoint inhibitors (ICI) in patients with non–small cell lung cancer (NSCLC) and also contribute to sex disparities in several cancers. Thus, we hypothesized that TP53 mutations might serve as sex-dependent genomic biomarkers of ICI treatment response in patients with NSCLC.
Materials and Methods
Clinical data of 100 patients with metastatic NSCLC treated with ICI monotherapy at Seoul National University Bundang Hospital (SNUBH) were retrospectively reviewed. Genomic and clinical datasets of The Cancer Genome Atlas and an ICI-treated lung cancer cohort (cBioPortal) were also analyzed.
Results
In SNUBH cohort, no statistically significant difference was observed in the median progression-free survival (PFS) according to TP53 mutation status (p=0.930); however, female patients with TP53 mutations (MT) had a significantly prolonged median PFS compared to wild-type (WT) (6.1 months in TP53 MT vs. 2.6 months in TP53 WT; p=0.021). Programmed death-ligand 1 (PD-L1) high (≥ 50%) expression was significantly enriched in female patients with TP53 MT (p=0.005). The analysis from publicly available dataset also revealed that females with NSCLC with TP53 MT showed significantly longer PFS than those with TP53 WT (p < 0.001). In The Cancer Genome Atlas analysis, expression of immune-related genes, and tumor mutation burden score in TP53 MT females were higher than in males without TP53 MT.
Conclusion
Female patients with NSCLC with TP53 mutations had high PD-L1 expression and showed favorable clinical outcomes following ICI therapy, suggesting a need for further research to explore the role of TP53 mutations for sex disparities in response to ICI therapy.
Introduction
Lung cancer is one of the leading causes of cancer-related mortality worldwide, with a 5-year survival rate of only 19%. Most patients with lung cancer (85%) suffer from non–small cell lung cancer (NSCLC) [1,2]. Immunotherapy, especially immune checkpoint inhibitor (ICI), has significantly improved clinical outcomes in patients with advanced NSCLC, achieving an outstanding duration of disease control compared to standard cytotoxic chemotherapy approaches [3-5].
Currently, only a small group (up to 20%) of patients achieve a long-standing response to ICIs, indicating the need to identify more effective biomarkers. Specific biomarkers such as programmed death-ligand 1 (PD-L1) expression, tumor mutation burden (TMB), immune signature, and mismatch repair are potential predictive factors of the benefits of ICI treatment [6]. To date, microsatellite instability status and strong PD-L1 expression in tumor cells are the only clinically recommended predictive biomarkers of ICI efficacy, as determined by polymerase chain reaction and immunohistochemistry, respectively [5,7]. Among the various known potential factors, specific genomic alterations can improve the efficacy of immunotherapy. An in-depth understanding of genetic alterations and molecular profiling based on next-generation sequencing (NGS) has dramatically changed the landscape of precision therapy for lung cancer.
The tumor suppressor gene TP53 is the most frequently observed mutated gene in cancer, resulting in the loss of tumor suppressor function or gain of oncogenic activity [8]. In addition to its role in DNA damage response, the p53 protein has also been attributed to intracellular metabolism, genetic and epigenetic stability, inflammatory remodeling, and integration with various pathways [9,10]. Therefore, understanding the relationship between p53 modification and tumor microenvironment (TME) may contribute to more effective cancer treatment. Recent studies have suggested TP53 mutations as prognostic biomarkers associated with longer survival in patients with NSCLC receiving ICI therapy [6,11]. In particular, multi-center data have shown that TP53 missense and nonsense mutations were significantly different in terms of associations with PD-L1 expression, interferon γ (IFN-γ) score, and TME composition [12]. However, the association between TP53 mutations and immunotherapy efficacy remains unclear.
Recent studies have suggested that TP53 mutations affect sex disparities in several types of cancers, including NSCLC [13], particularly causing sex-based differences in the immune features of the TME. Moreover, a meta-analysis of multiple immune features in patients with NSCLC provided a more comprehensive perspective for investigating sex bias in immunotherapy [14]. These observations emphasize the importance of sex in stratifying patients using molecular markers. While the association between TP53 mutations and immune responses is emerging, the impact of patient sex on clinical outcomes is yet to be determined. Therefore, in this study, we hypothesized that TP53 mutation might serve as a sex-dependent genomic biomarker for ICI treatment response in patients with NSCLC.
Materials and Methods
1. Patients and methods
This was a single-center retrospective study of patients with metastatic NSCLC treated with ICI therapy at Seoul National University Bundang Hospital (SNUBH), between January 2016 and July 2021. Initial samples of 142 patients were identified through a database. Eligibility criteria for participation in this study included patients aged ≥ 20 years with histologically confirmed NSCLC—metastatic or recurrent disease—and ≥ 1 measurable lesion according to the Response Evaluation Criteria in Solid Tumors (RECIST) ver. 1.1 [15], with NGS results using three different NGS assays including ver. 1 (n=43), ver. 2 (n=37) or TruSight Oncology 500 Assay (n=20). Patients were excluded if they had any systemic diseases or other malignancies that could influence the survival analysis, received only one treatment cycle of ICI, or had no information on tumor response. The final cohort included in the data analysis were 100 patients and clinical, biological, pathological, and molecular data, including tumor PD-L1 expression and TP53 mutation status assessed by NGS, were retrospectively collected from the patient’s electronic medical records (Fig. 1). PD-L1 expression was evaluated in all tumors by assessing the percentage of tumor cells with membranous staining [16].
2. Clinical outcomes
The primary endpoint was progression-free survival (PFS), defined as the time from ICI treatment initiation to disease progression or death due to any cause, whichever occurred first. Patients who exhibited no progression were censored at the last follow-up. Secondary endpoints included overall survival (OS), defined as time from ICI initiation to death due to any cause, objective response rate (ORR: complete response [CR] or partial response [PR]), and disease control rate (DCR: CR, PR or stable disease [SD] that lasted > 6 months). The radiographic best response was determined based on RECIST ver. 1.1 criteria. Patients who were surviving at the last follow-up were censored.
3. Publicly available dataset
Genomic and clinical data of a Memorial Sloan Kettering Cancer Center (MSKCC) ICI-treated cohort were obtained from the cBioPortal database [17,18]. This cohort included 240 patients with advanced NSCLC who were treated with anti–PD-L1 monotherapy or combination therapy with anti–cytotoxic T-cell lymphocyte-4, between April 2011 and January 2017. Patients with tumors molecularly profiled using MSK-IMPACT were included in the study. Demographic data, including age, smoking status, histologic subtype, ICI treatment type, number of treatment lines, and PD-L1 expression were provided; TMB and molecular data of TP53 mutation status were retrieved from the NGS data. Clinical outcomes included durable clinical benefit (DCB), which was defined as PR/SD that lasted for > 6 months, and PFS. To assess the differences in immune-related gene mutations between sex and TP53 mutation status, we obtained The Cancer Genome Atlas (TCGA)–Lung Adenocarcinoma (LUAD) and TCGA–Lung Squamous Cell Carcinoma (LUSC) cohorts, including mRNA expression profiling data and somatic mutation data after the removal of epidermal growth factor receptor (EGFR) somatic variants and anaplastic lymphoma kinase fusion cases. TCGA-LUAD and TCGA-LUSC cohorts were combined into the TCGA-NSCLC cohort for subsequent analysis. We evaluated IFNG, CXCL9, PD-L1, LAG3 genes, IFN-γ signature score, and TMB; the unit of gene expression was pan-cancer normalized log2 (fragments per kilobase of exon model per million mapped fragments+1).
4. Statistical analysis
The mean and median values of continuous variables were calculated for the population description, whereas categorical variables were described as frequencies and percentages. For baseline categorical variables, the four groups of patients with NSCLC based on TP53 mutation status and sex were compared using the chi-squared test, or Fisher’s exact test for smaller sample sizes. Kaplan-Meier analysis and the log-rank test were used to estimate the median survival times and their 95% confidence intervals (CI) and statistical significance for the four subgroups: TP53-mutated (MT) females, TP53 MT males, TP53 wild-type (WT) females, and TP53 WT males. Hazard ratios (HR) and their 95% CI were calculated using univariate analysis with Cox proportional hazard models. Variables with > 5% missing data were not included in the multivariate analysis. The variables tested in the univariate analyses of PFS were age, sex, smoking status, histological subtype, the number of lines of prior treatment line, PD-L1 expression level, ICI treatment type, and EGFR mutation status. Only variables that showed an association with survival (p < 0.10) were introduced in the multivariate models of PFS, with three additional clinical variables (age, smoking status, and PD-L1 expression level) forced into the PFS multivariate model because of their clinical interest. Survival curves were plotted using GraphPad Prism (ver. 8.0.1, GraphPad Software Inc., San Diego, CA). The Mann-Whitney U test was performed between subgroups from TCGA for gene expression analysis. The findings were considered statistically significant if the two-sided p-value was < 0.05. All statistical analyses were performed using the IBM SPSS ver. 26 (Armonk, New York, NY).
Results
1. Baseline characteristics of SNUBH cohorts
In total, 100 patients were included in this study, as shown in Table 1. At ICI initiation, 45 (45.0%) patients were over 65 years old; 63 (63.0%) were male, and 63 (63.0%) were former smokers. Seventy-four patients (74.0%) had adenocarcinoma, 15 (15.0%) had squamous cell carcinoma, and 11 (11.0%) had other histological subtypes. Most patients (n=96) had distant metastases. ICI was administered as first-line therapy in 19 patients (19.0%), while the remaining patients (n=81) had received one or more treatments before receiving ICI treatment. The majority of patients (n=88, 88.0%) received ICI monotherapy. Fifteen patients (15.0%) harbored EGFR-activating mutations and 59 patients (59.0%) harbored TP53 mutations. We compared the distribution of clinical features between TP53 mutant type and WT patients. There were no significant differences in the characteristics between patients with and without TP53 mutations according to age, performance status, histologic type, initial TNM stage, a type of treatment, number of prior treatments, PD-L1 expression level, or EGFR mutation status. However, male patients (71.2% vs. 51.2%, p=0.042) and former smokers (72.9% vs. 48.8%, p=0.014) were more frequent in the TP53 mutation group. To evaluate the sex differences in the immunotherapy response based on TP53 mutation status, we divided the patients into four groups: female patients with TP53 MT (n=17, 17.0%), female patients with TP53 WT (n=20, 20.0%), male patients with TP53 MT (n=42, 42.0%), and male patients with TP53 WT (n=21, 21.0%) (Fig. 1).
2. ORR and DCR of SNUBH cohorts
Overall, ORR was 20.0% and DCR was 43.0%. In subgroup analysis, there was no difference in ORR based on TP53 mutation status in both female patients and male patients (female: 35.3% in TP53 MT vs. 10.0% in TP53 WT, p=0.109; male: 16.7% in TP53 MT vs. 23.8% in TP53 WT, p=0.513). However, in female patients, the DCR was significantly higher in the TP53 mutated group than in the TP53 WT group (76.5% vs. 25.0%, p=0.002). In contrast, no significant difference in DCR was observed between TP53 mutated and TP53 WT male patients (33.3% vs. 52.4%, p=0.145) (Table 2).
3. Survival outcomes of SNUBH cohorts
Overall, the median follow-up was 32.1 months (95% CI, 18.1 to 46.1). Median OS was 13.6 months (95% CI, 10.0 to 17.1), and median PFS was 3.4 months (95% CI, 2.3 to 4.4). According to the TP53 mutation status, median PFS in the TP53 mutated group was 3.0 months (95% CI, 1.6 to 4.4) versus 3.9 months in the TP53 WT group (95% CI, 1.3 to 6.6; p=0.930) (Fig. 2A). In a univariate analysis for PFS, there was no significant difference in the median PFS between patients with or without TP53 mutations in all patients (S1 Table). Moreover, the multivariable analysis suggested that female sex had no significant effect on ICI benefit (HR, 0.42; 95% CI, 0.13 to 1.44; p=0.168) (S1 Table). However, when stratified by both TP53 and sex, female patients with TP53 mutation had a significantly prolonged median PFS compared to WT female patients (female with TP53 MT, 6.1 months; 95% CI, 2.2 to 10.1 and TP53 WT, 2.6 months; 95% CI, 1.5 to 3.7; p=0.021) (Fig. 2B). In multivariable analysis, female sex was independently associated with longer PFS in TP53 mutant patients (HR, 0.15; 95% CI, 0.03 to 0.82; p=0.029) (Table 3). Whereas no significant difference in median PFS was found between males with TP53 mutation and TP53 WT (male with TP53 MT, 2.2 months; 95% CI, 1.6 to 2.9 and TP53 WT, 6.6 months; 95% CI, 0.0 to 14.6 months; p=0.183) (Fig. 2C). OS was significantly higher in patients with TP53 WT compared to TP53 mutated patients (patients with TP53 MT, 12.2 months; 95% CI, 8.8 to 15.5 and TP53 WT, 25.1 months; 95% CI, 7.0 to 38.8). However, we did not find any statistically significant interaction between TP53 status and patients’ sex for OS (S2 Fig.).
4. MSKCC public data analysis
To avoid accidental bias caused by a single cohort, we further verified the association between TP53 mutations and sex in the ICI-treated cohort from the MSKCC. We used cBioPortal to collect a published cohort study of patients with NSCLC, including 240 patients receiving ICIs with available mutational and clinical data, to further explore the relationship between gene mutations and ICI prognosis. Overall, more than half of the patients were aged over 65 years old at ICI initiation, and 50.8% of the patients were female. There were 141 (58.8%) TP53 mutations and 99 (41.2%) were TP53 WT (S3 Table). The analysis showed that, when stratified by sex, female patients with TP53 mutation showed significantly longer PFS than that of TP53 WT female patients, while there was no significant difference in PFS between male patients with or without vs. female with TP53 mutations (female with TP53 MT, 5.4 months [95% CI, 3.3 to 7.5], vs. female with TP53 WT, 2.4 months [95% CI, 2.0 to 2.8]; p < 0.001 and male with TP53 MT, 3.6 months [95% CI, 2.7 to 4.6], vs. male with TP53 WT, 2.8 months [95% CI, 0.0 to 5.7]; p=0.760) (Fig. 2D). These results are consistent with the univariate results of our cohorts. Similar to our results, when stratified by sex and TP53 mutation status, a significantly higher DCB was observed in female patients with TP53 mutations than in those with TP53 WT (38.2% vs. 20.4%, p=0.003) (Table 4). Moreover, in multivariate analysis performed in the TP53 mutant group, the female sex showed an independent prognostic factor for longer PFS (HR, 0.67 [95% CI, 0.45 to 1.00]; p=0.048) (Table 5).
5. Differences in PD-L1 expression level between male and female patients from SNUBH cohorts
Comparing groups divided by sex in each group with or without TP53 mutation status, PD-L1 high (≥ 50%) expression was significantly enriched in female patients with TP53 mutation compared to those with male patients (58.8% vs. 16.7%, p=0.005), whereas there was no significant difference in PD-L1 expression between the male and female patients in TP53 WT group (p=0.161) (Table 6). Fig. 3 shows that the distribution of TP53 mutations differed significantly in terms of their association with PD-L1 expression between male and female patients (p=0.005). Among the TP53 missense mutant group, the proportion of patients with a PD-L1 score ≥ 50 was prominently higher in the female group than in the male group (p=0.033). Moreover, none of the female patients with TP53 nonsense mutations had a PD-L1 score of < 1.
6. Expression level of immune-related genes, INF-γ signature, and TMB score from TCGA
Given the potential sex differences in TP53 mutation, we explored the correlation of IFNG, CXCL9, LAG3, and PDL1 expression levels, IFN-γ signature, and TMB score stratified by sex (Sub1: female patients with TP53 mutation; Sub2: female patients without TP53 mutation; Sub3: male patients with TP53 mutation; Sub4: male patients without TP53 mutation). In the TCGA-NSCLC cohort, both female and male patients with TP53 mutations had higher TMB than TP53 WT patients (Sub1 vs. Sub2, p < 0.001; Sub3 vs. Sub4, p < 0.001). TP53 mutant male patients showed higher LAG3 expression and IFN-γ signature scores than TP53 WT male patients. However, there were no significant differences in IFNG, CXCL9, LAG3, and PDL1 gene expression levels, and IFN-γ signature score between female patients with or without TP53 mutation. We found that the expression of IFNG, CXCL9, LAG3, PDL1, and IFN-γ signature score, and TMB score in TP53 mutant female patients were higher than in male patients without TP53 mutation (Fig. 4).
Discussion
In this study, we showed that female patients with TP53 mutations had higher PD-L1 expression, increased immune-related gene expression, and favorable clinical benefits of ICI treatment. Female NSCLC patients with TP53 mutations had a significantly higher DCR and longer PFS than those with TP53 WT, while no significant difference was observed in male patients. Multivariate analysis further demonstrated that female sex was a significant prognostic factor for longer PFS in the TP53 mutation group. These findings suggest that female sex could potentially serve as a predictive biomarker for ICI therapy in TP53 mutant patients with metastatic NSCLC. Moreover, analysis of an independent dataset from the MSKCC cohort further supported our results that female NSCLC patients with TP53 mutations had a better prognosis than female patients without TP53 mutations in NSCLC.
To the best of our knowledge, this is the first study to demonstrate that TP53 mutations induce sex-dependent ICI responses. In our study, we observed a notable increase in the prevalence of PD-L1 high (≥ 50%) expression among female patients with TP53 mutations, as compared to male patients (58.8% vs. 7.0%, p=0.005). In addition, PD-L1 expression did not correlate with survival, whereas multivariate analysis showed that female sex was independently associated with better PFS in TP53 mutant patients, suggesting that PD-L1 did not represent a confounding factor. These results suggest that TP53 mutations in female patients could be involved in regulating PD-L1 expression.
TP53, a tumor suppressor, is commonly mutated in cancers. Recently, several studies have provided compelling evidence that TP53 mutations are potential prognostic biomarkers for patients with cancers such as head and neck squamous cell cancer and lung adenocarcinoma who undergo ICI treatment [6,11,19]. Dong et al. [6] retrospectively showed that in 30 metastatic NSCLC patients treated with pembrolizumab, the median PFS was significantly longer in the TP53-mutated group than in the TP53 WT group (14.5 months vs. 3.5 months, p=0.042). Additionally, Assoun et al. [11] demonstrated that TP53 gene mutations in patients with metastatic NSCLC (n=72) treated with ICI were associated with better OS, PFS, and ORR. In contrast, survival analysis of the MSK-IMPACT cohort of 350 NSCLC patients who received immunotherapy showed that NSCLC patients with TP53 truncating mutations had a significantly lower OS than those with WT TP53 (9 months vs. 14 months; p=0.019) [20]. Therefore, the association between TP53 mutations and immunotherapy efficacy remains controversial. In our study, overall patients with TP53 mutations showed no survival benefit from immunotherapy compared to those with TP53 WT. However, among TP53 mutant patients, when stratified with patients’ sex, the female patients experienced significantly longer PFS compared to male patients (6.1 months vs. 2.2 months, p=0.039). In addition, we observed the prognostic role of TP53 mutations in the context of sex from survival analysis of the MSKCC cohort of 240 patients with metastatic NSCLC who were treated with ICIs. This suggested that there might be sex biases of TP53 mutation, in particular for female NSCLC patients who were treated with immunotherapy.
Interactions between sex and p53 pathway genes have been observed in the following literature. Firstly, non-expressed mutations (NEMs) in the X chromosome are more common in females than in males, and a significant number of these NEMs are found in genes related to the p53 network, such as E3 Ubiquitin Protein Ligase 1 (HUWE1) and ATP-dependent helicase ATRX (ATRX). This suggests that females have a protective effect that limits the expression of somatic mutations [13]. Secondly, WT p53 can inhibit tumor development through various pathways [8]. However, the loss of p53 as a tumor suppressor leads to the stochastic reactivation of the inactive X chromosome (Xi) and the expression of X-linked genes from both alleles. By increasing the likelihood of stochastic transcriptional events in Xi, p53 may contribute to the sex bias observed in various cancer types [21].
Notably, immunity differs between the sexes. Females generally exhibit better adaptive and innate immune responses than males, resulting in enhanced antiviral T-cell immunity and anticancer responses [22,23]. Moreover, immune infiltration is more abundant in female tumors, and all T-cell subpopulations are significantly enriched in the TME of females [24]. In females, the TME of NSCLC is distinguished by a higher prevalence of T-cell dysfunction, increased expression of inhibitory immune checkpoint molecules, and a greater presence of immunosuppressive cells. By contrast, the TME of NSCLC in males is characterized by a lower abundance of several innate and adaptive immune cell subsets [24].
We also have analyzed the TCGA-NSCLC data for the gene expression of immune-related genes (such as IFNG, CXCL9, LAG3, and PDL1), IFN-γ signature, and TMB score considering sex and TP53 mutation status. The results showed that NSCLC with TP53 mutation had high immune-related gene expression and high TMB as noted by previous studies. However, the gene expression of immune-related genes seemed similar regardless of TP53 mutation status in female patients. TP53 mutant female patients showed higher PD-L1 expression levels than male TP53 WT patients (p=0.002). The correlation of high PD-L1 expression level and clinical benefit in TP53 mutant female patients might be a possible explanation for sex disparities in TP53 mutant NSCLC patients treated with ICI.
Our study has several limitations. First, it was retrospective in nature, with a lack of blinded independent assessment of response, although ORR and DCR were secondary endpoints. Second, PFS prolongation in female patients with TP53 mutations has not translated into a longer OS. Although we cannot directly answer this point, post-progression survival of the patients was probably affected by other prognostic factors, including subsequent treatment. Third, the limited size of our single-center series could have resulted in bias and a lack of statistical power. In the current study, a profound difference was observed in the median PFS between male patients with TP53 mutation and TP53 WT, but the p-value was not statistically significant. This could be due to the small sample size of the group. Moreover, co-mutations in TP53 and KRAS were proven to be promising biomarkers for ICI response [6]. However, the presence of co-mutations was not evaluated in this study. Future studies with larger sample sizes and prospective designs are needed to confirm these findings and evaluate the clinical implications of TP53 mutation status as a predictive biomarker for ICI therapy in patients with NSCLC.
In conclusion, we found that female NSCLC patients with TP53 mutation had high PD-L1 expression levels and showed favorable clinical outcomes after ICI therapy. The role of TP53 mutations in the sex disparities in patients with NSCLC treated with ICI needs to be further validated.
Electronic Supplementary Material
Supplementary materials are available at Cancer Research and Treatment website (https://www.e-crt.org).
Notes
Ethical Statement
This study was approved by the Institutional Review Board of the SNUBH (IRB number: B-2104-681-104), and the requirement for informed consent was waived because of the retrospective nature of this study. This study was conducted in accordance with the Declaration of Helsinki and Good Clinical Practice guidelines.
Author Contributions
Conceived and designed the analysis: Choi S, Kim SH, Lee S.
Collected the data: Choi S, Kim SH.
Contributed data or analysis tools: Choi S, Kim SH, Lee S, Seo J, Kang M, Jung EH, Kim SA, Suh KJ, Lee JY, Kim JW (Ji-Won Kim), Kim JW (Jin Won Kim), Lee JO, Kim YJ, Lee KW, Kim JH, Bang SM, Lee JS.
Performed the analysis: Choi S, Kim SH, Lee S.
Wrote the paper: Choi S, Kim SH.
Conflicts of Interest
Conflict of interest relevant to this article was not reported.
Acknowledgements
This study was supported by the SNUBH Research Fund [grant no.: 02-2020-0048].