Purpose: Promising results in the treatment of non–small cell lung cancer (NSCLC) have been seen with agents targeting immune checkpoints, such as programmed cell death 1 (PD-1) or programmed death ligand-1 (PD-L1). However, only a select group of patients respond to these interventions. The identification of biomarkers that predict clinical benefit to immune checkpoint blockade is critical to successful clinical translation of these agents.
Methods: We conducted an integrated analysis of three independent large datasets, including The Cancer Genome Atlas of lung adenocarcinoma and two datasets from MD Anderson Cancer Center (Houston, TX), Profiling of Resistance Patterns and Oncogenic Signaling Pathways in Evaluation of Cancers of the Thorax (named PROSPECT) and Biomarker-Integrated Approaches of Targeted Therapy for Lung Cancer Elimination (named BATTLE-1). Comprehensive analysis of mRNA gene expression, reverse-phase protein array, IHC, and correlation with clinical data were performed.
Results: Epithelial–mesenchymal transition (EMT) is highly associated with an inflammatory tumor microenvironment in lung adenocarcinoma, independent of tumor mutational burden. We found immune activation coexistent with elevation of multiple targetable immune checkpoint molecules, including PD-L1, PD-L2, PD-1, TIM-3, B7-H3, BTLA, and CTLA-4, along with increases in tumor infiltration by CD4+Foxp3+ regulatory T cells in lung adenocarcinomas that displayed an EMT phenotype. Furthermore, we identify B7-H3 as a prognostic marker for NSCLC.
Conclusions: The strong association between EMT status and an inflammatory tumor microenvironment with elevation of multiple targetable immune checkpoint molecules warrants further investigation of using EMT as a predictive biomarker for immune checkpoint blockade agents and other immunotherapies in NSCLC and possibly a broad range of other cancers. Clin Cancer Res; 22(14); 3630–42. ©2016 AACR.
See related commentary by Datar and Schalper, p. 3422
This article is featured in Highlights of This Issue, p. 3417
Therapeutic agents targeting immune checkpoints, such as programmed cell death 1 (PD-1), have achieved encouraging clinical activity and been granted recent approval in the treatment of non–small cell lung cancer (NSCLC). However, it appears that only a select group of patients respond to these interventions. The identification of potential biomarkers that predict clinical benefit to immune checkpoint blockade is critical to successful clinical application of these agents. By analyzing three large independent datasets, we report that lung adenocarcinomas displaying a “mesenchymal” phenotype are associated with distinct tumor microenvironment changes, including elevated expression of multiple immune checkpoints, such as PD-1 and PD-L1, along with evidence of preexisting immunity and increases in tumor infiltration by CD3+ T cells and CD4+Foxp3+ regulatory T cells. Our data warrant further investigation of using EMT as a potential predictive biomarker to guide selection of patients who are likely to benefit from immune checkpoint blockade agents in NSCLC and possibly a range of other cancers.
Non–small cell lung cancer (NSCLC), long thought to be a “nonimmunogenic” tumor, has shown responses to immunotherapeutic approaches targeting immune checkpoints programmed cell death 1 (PD-1), programmed death ligand-1 (PD-L1), and cytotoxic T lymphocyte–associated protein 4 (CTLA-4), even in heavily pretreated patients (1–3). These findings clearly indicate the crucial role of immune checkpoint pathways in mediating immune tolerance in NSCLC. The PD-1/PD-L1 pathway has emerged as a critical inhibitory pathway that regulates T-cell response and maintains immune suppression (1, 2, 4, 5). Although anti-PD-1/PD-L1 treatment can produce durable responses, it appears to benefit only a subset of patients. Although the expression levels of PD-L1 on tumor cells and tumor-infiltrating immune cells have recently been shown to correlate with clinical response to anti-PD-1 therapy (6–8), only a subset of patients with PD-L1–expressing tumors had clinical response and others without PD-L1 staining demonstrate clinical benefit, indicating that additional factors in the tumor microenvironment exist, which define the subgroup of patients who derive benefit (6, 9).
Epithelial–mesenchymal transition (EMT) is a key process that drives cancer metastasis, drug resistance, and has been associated with poor prognosis in multiple cancers, including NSCLC (10–12). However, most studies that investigate the impact of EMT in cancer have mainly focused on metastasis and drug resistance (12, 13). The impact of EMT on reprogramming the tumor immune microenvironment is largely unknown. Our group has previously developed a robust EMT gene signature that correlates with cellular phenotypes and predicts in vitro and in vivo NSCLC resistance to EGFR or PI3K/Akt inhibitors, highlighting differential patterns of drug responsiveness for “epithelial” and “mesenchymal” NSCLC (14). We have also recently demonstrated a molecular link between EMT and intratumoral CD8+ T-cell suppression, through the regulation of PD-L1 in both animal models and human cell lines (15). To further explore the impact of EMT on tumor immune microenvironment and identify potential biomarkers for selecting patients who might preferentially benefit from PD-1/PDL-1 blockade or immunotherapies more broadly, we conducted an integrated analysis of three independent large patient datasets.
Patients and Methods
Clinical datasets and patient sample characteristics
A total of 439 patients from three clinical datasets, including The Cancer Genome Atlas (TCGA; n = 230), Profiling of Resistance Patterns and Oncogenic Signaling Pathways in Evaluation of Cancers of the Thorax (named PROSPECT hereafter; n = 152), and the Biomarker-Integrated Approaches of Targeted Therapy for Lung Cancer Elimination (named BATTLE-1 hereafter; n = 57), were included in this study but analyzed as independent datasets (16). Clinical and demographic information for 230 early-stage, surgically resected lung adenocarcinomas included in the TCGA dataset was obtained from the TCGA portal (https://tcga-data.nci.nih.gov/tcga/dataAccessMatrix.htm). The PROSPECT dataset includes tumor tissue collected from patients who underwent surgical resection of NSCLC with curative intent between 1996 and 2008 at MD Anderson Cancer Center (Houston, TX). A total of 152 tumors classified as adenocarcinoma resected from patients who had not received neoadjuvant chemotherapy or radiotherapy and had detectable tumor by H&E immunostaining were included in this study. Case selections were made by two independent pathologists who evaluated the amount of tumor tissue present in H&E slides corresponding to the paraffin blocks from surgically resected lung specimens available for this study. Cases with at least one representative tumor tissue block containing >30% viable tumor cells were selected for further immunohistochemical analysis. BATTLE-1 was a randomized, biomarker-based clinical trial for patients with recurrent or metastatic NSCLC in the second-line setting (trial registration ID: NCT00409968). Core needle biopsies were required for enrollment. Samples with detectable adenocarcinoma by H&E immunostaining were included in our analysis (n = 57). Clinical characteristics of these samples are described in Supplementary Table S1. The study was approved by the Institutional Review Boards of MD Anderson Cancer Center (Houston, TX).
mRNA expression profiling
Experimental details regarding lung adenocarcinoma from TCGA dataset, including RNA extraction from tumors, mRNA library preparation, sequencing (Illumina HiSeq platform), quality control, data processing, and quantification of gene expression, have been reported previously (17). For the PROSPECT samples, the mRNA was extracted from frozen tumor tissue corresponding to the same specimen from which the formalin-fixed, paraffin-embedded blocks were made. The frozen tumor tissue was harvested by a pathologist who examined and sampled the gross tissue specimen, and only tumor tissues were included for mRNA extraction. Array-based expression profiling of PROSPECT tumors was performed using the Illumina Human WG-6 v3 BeadChip, according to the manufacturer's instructions. BATTLE-1 adenocarcinoma was profiled using the GeneChip Human Gene 1.0 ST Array from Affymetrix. Gene expression data for the PROSPECT and BATTLE-1 dataset have been previously deposited in the GEO repository GSE42127 and GSE33072 (18, 19), respectively. For each sample, an EMT score was computed using an averaging scheme based on the mRNA expression of 76 genes previously published by our group and originally derived from the analysis of NSCLC (14, 15). The scores were calculated as the average expression level of “mesenchymal” genes minus the average expression level of “epithelial” genes. Tumor samples were then classified by EMT score as EMT-low (defined by EMT scores ≤ lowest 1/3) or EMT-high (defined by EMT scores ≥ highest 1/3) in both TCGA and PROSPECT datasets. By comparing the tumor samples with either relatively high or low EMT score, but not intermediate, we expected to increase the likelihood of identifying the immune markers associated with either “epithelial” or “mesenchymal” lung adenocarcinomas. In the BATTLE-1 dataset of metastatic tumor specimens, the sample size was much more limited (n = 57), and a narrowed spectrum of EMT scores was observed, with generally more “mesenchymal” phenotypes. This is likely due to the nature of advanced stage/metastatic disease in this patient population and the small amount of sampled tissue for each case versus resection specimens. On the basis of this narrowed distribution, a different EMT threshold was used in the BATTLE-1 dataset, defining EMT score < 0 as EMT-low (“epithelial”) and EMT score ≥ 0 as EMT-high (“mesenchymal”). The mRNA expression profiles of immune-related genes in each sample were then analyzed. The raw data files of transcriptomes were analyzed using Bioconductor R packages.
Nonsynonymous mutational burden per Mb of DNA and C>A transversion rate in TCGA tumors were calculated using the MutSig2CV algorithm as published previously (17). This algorithm takes into account recurrence of mutations, nucleotide context, gene expression, replication time, and somatic background mutation rate. Genes with a Bonferroni-corrected P < 0.025 were deemed significant. It has been previously established that transversion low samples represent tumors from lifelong never-smokers and transversion high samples as tumors from patients with 60 or more pack years (17). A linear discriminant analysis based on GC>AT frequency, GC>TA frequency, and total mutation count by using the R MASS library2 was performed based on previously published methods to classify all samples as belonging to either transversion high or transversion low categories.
Immunohistochemical staining and reverse-phase protein array analysis
For the detailed immunohistochemical and reverse-phase protein array (RPPA) analysis of PD-L1, see the Supplementary Methods section. Briefly, formalin-fixed, paraffin-embedded tumor blocks from the PROSPECT dataset were selected for immunohistochemical staining. Thirty-five tumor tissues with “mesenchymal” lung adenocarcinomas and 33 tumor tissues with “epithelial” lung adenocarcinomas, as defined by the mRNA expression, from the PROSPECT tissue bank were found to have sufficient quality material for additional analyses and were blinded to the pathologists who performed the IHC study. The markers assessed in this study included PD-L1, PD1, CD3, CD4, CD8, CD45RO, CD57, CD68, granzyme B, and FOXP3 (see Supplementary Methods for antibody details). All antibodies were detected with the Leica Bond Polymer Refine Detection Kit (Leica Microsystems), including diaminobenzidine reaction to detect the antibody labeling and hematoxylin counterstaining. The stained slides were digitally scanned using the Aperio ScanScope Turbo slide scanner (Leica Microsystems). The digital images were captured at 20× magnification.
For RPPA analysis of TCGA and PROSPECT specimens, the tumors were lysed as we have published previously (20). Five serial dilutions of each protein lysate were printed on nitrocellulose-coated slides using an Aushon Biosystems 2470 Arrayer and stained sequentially with primary and secondary antibodies in an autostainer (BioGenex), prior to signal detection using a signal amplification system and DAB-based colorimetric reaction. MicroVigene Software (VigeneTech), as well as an in-house R package, was used to assess spot intensities, and the SuperCurve method was applied to estimate protein levels in each sample. For comparisons, data were log transformed (to the base of 2) and median centered across antibodies to correct for protein loading. Rabbit polyclonal antibody to PD-L1 (cat: ab174838, Abcam) at 1:250 dilution was used to detect PD-L1 expression by RPPA. Differences in protein expression between “epithelial” and “mesenchymal” lung adenocarcinoma tumor samples were compared by t test. Pearson correlation between CD274 (PD-L1) mRNA expression and PD-L1 protein expression by RPPA were assessed, along with E-cadherin protein RPPA and PD-L1 RPPA expression. All statistical analyses were performed using R packages (version 2.10.0).
Statistical analyses were conducted using GraphPad Prism version 6.00 (GraphPad Software) or, alternatively, the R system (version 2.10.0) for statistical computing. The unpaired t test was used for comparisons between two-group means, where the data could be assumed to have been sampled from populations with normal (or approximately normal) distributions. Mann–Whitney U test was used to compare the mean ranks between two groups. All P values are two tailed, and for all analyses, P ≤ 0.05 is considered statistically significant, unless otherwise specified.
PD-1, PD-L1, and PD-L2 are significantly elevated in “mesenchymal” lung adenocarcinoma
The PD-1:PD-L1/PD-L2 pathway has been implicated in tumor escape from immune destruction in various cancers, including NSCLC (21). Although the expression of PD-L1 in NSCLC is an area of active study, whether specific subgroups of patients with NSCLC display distinct PD-L1 expression patterns is unknown. Our recent report demonstrated that “epithelial” and “mesenchymal” NSCLC have distinct patterns of drug responsiveness to EGFR and PI3K/Akt inhibitors. To study the expression of the PD-1:PD-L1/PD-L2 axes in “epithelial” and “mesenchymal” lung adenocarcinomas, mRNA expression data from TCGA and PROSPECT datasets were analyzed. As shown in Fig. 1A and B, the expression levels of PD-L1, PD-L2, and PD-1 were significantly elevated in “mesenchymal” versus “epithelial” tumors in both TCGA and PROSPECT, consistent with our recent findings (15). To confirm whether the elevated mRNA expression of PD-L1 and PD-1 correlated to protein levels, immunohistochemical staining and scoring of PD-L1 (on tumor cells) and PD-1 (on tumor-infiltrating immune cells) were performed for full tumor sections from a subset of cases in the PROSPECT dataset, for which we have sufficient tumor tissues of high quality for immunohistochemical analysis (n = 68). For these analyses, we validated the specificity of the antibody E1L3N (Cell Signaling Technology), which we employed for the immunohistochemical PD-L1 staining (Supplementary Fig. S1A–S1G). As shown in Fig. 1C and D, we found significantly higher expression of PD-L1 and PD-1 in “mesenchymal” than in “epithelial” adenocarcinoma cases (P = 0.005 and P = 0.015, respectively). We also addressed this question by using RPPA, an independent technique to detect and quantify the protein levels. We first validated the specificity of the antibody (Abcam 174838) used for RPPA versus the widely used 5H1 antibody, by Western blot analysis of cell line lysates, IHC of cell pellets, IHC of control placenta samples, and tumor samples with a range of PD-L1 expression on tumor cells and infiltrates (Supplementary Fig. S2A–S2H). The PD-L1 protein levels measured by RPPA correlated strongly with the PD-L1 mRNA levels (CD274) in both independent datasets (Fig. 1E and F). We further demonstrated a positive correlation between PD-L1 protein by RPPA and EMT score, when assessed as a continuous variable, or significantly higher levels of PD-L1 in “mesenchymal” as compared with “epithelial” tumors, when analyzed by group (Fig. 1E and F). The expression of E-cadherin, an epithelial marker that is also present on the RPPA panel, strongly negatively correlated with PD-L1 protein expression in TCGA samples (Fig. 1E) and the independent PROSPECT dataset (Fig. 1F).
A distinct tumor microenvironment immune profile with elevation of multiple immune checkpoint molecules is revealed in “mesenchymal” lung adenocarcinoma
The findings of elevated PD-1:PD-L1/PD-L2 axis in “mesenchymal” lung adenocarcinoma prompted us to investigate whether there is an association with a broader immunosuppressed phenotype beyond the PD-1:PD-L1/PD-L2 axis. To address this question, a comprehensive list of 89 immune-related genes, including costimulatory molecules, immune checkpoints, cytokines, chemokines, MHC class I and II, and genes highly expressed on dendritic cells, T cells, natural killer (NK) cells, myeloid cells, and macrophages, was generated from the literature (Supplementary Table S2; ref. 22). There is no overlap between the immune-related gene list and the previously identified EMT gene signature (19). The mRNA expression of each immune-related molecule was tested in the “mesenchymal” versus “epithelial” adenocarcinomas in TCGA and PROSPECT datasets. Strikingly, profound immune-related phenotypic changes were found in “mesenchymal” lung adenocarcinoma in contrast to much lower expression of immune-related molecules in “epithelial” lung adenocarcinoma (Fig. 2A and B). Notably, multiple immune checkpoint molecules were significantly elevated in TCGA and PROSPECT, including T-cell immunoglobulin and mucin protein-3 (TIM-3), B- and T-lymphocyte attenuator (BTLA), CTLA-4, lymphocyte activation protein 3 (LAG-3), and B7-H3 (Fig. 2C and D). Similar to PD-1, TIM-3, BTLA, CTLA-4, LAG-3, and B7-H3 all negatively regulate T-cell function through relatively unique and potentially nonoverlapping molecular mechanisms (23). Coexpression of multiple immune checkpoint molecules has been frequently found on exhausted T cells in tumors and chronic infections (23). Although the expression pattern of elevated immune checkpoint molecules was largely consistent between TCGA and PROSPECT, some differences were observed. For example, B7-H3 was significantly elevated in lung adenocarcinoma in TCGA dataset, but not in the PROSPECT dataset. By contrast, herpes virus entry mediator was elevated in PROSPECT, but not in TCGA. Although both datasets comprised patients with mainly surgically resectable disease, more patients with stage IV (4%) were included in TCGA than in PROSPECT (0%), indicating that potential differences in the patient population might contribute to the observed molecular differences. Of note, although the association between this distinct tumor microenvironment immune profile and “mesenchymal” lung adenocarcinoma was most readily observed when the data were analyzed as the top versus bottom 1/3 group of cases, we observed similar results when analyzed in either a continuous manner or categorized into EMT-low, EMT-intermediate, and EMT-high groups (Supplementary Fig. S3A–S3D).
Significant increase in tumor-infiltrating CD4+FOXP3+ regulatory T cells is found in “mesenchymal” lung adenocarcinoma
The expression of multiple immunosuppressive molecules in “mesenchymal” lung adenocarcinoma suggested the possibility of a complex immunosuppressive tumor microenvironment. To explore the cellular consequence of this observation, “mesenchymal” and “epithelial” lung adenocarcinoma specimens from the PROSPECT study were stained and scored in a blinded fashion for CD3, CD4, CD8, FOXP3, CD68, CD45RO, CD57 and, granzyme B, well-established IHC markers for T-cell subsets, regulatory T cells, macrophages, and NK cells, respectively. As shown in Fig. 3A, a robust increase in tumor-infiltrating CD3+ T cells was found in “mesenchymal” as compared with “epithelial” lung adenocarcinoma (P = 0.039). We found much higher levels of tumor infiltration by CD4+ and FOXP3+ T cells in “mesenchymal” lung adenocarcinomas (P = 0.009 and P = 0.030, respectively; Fig. 3B and C). These likely represent CD4+FOXP3+ regulatory T cells, although we are unable to confirm that interpretation as methodology for double staining cells was not used. Although there was a trend toward more tumor-infiltrating CD8+ T cells in “mesenchymal” samples, the difference was not statistically significant. As shown in Fig. 3E, we observed no significant difference in tumor infiltration of total CD68+ macrophages, suggesting that the CD4+ and FOXP3+ T cells rather than tumor-associated macrophages likely played a greater role in this EMT-associated immunosuppressive phenotype. However, additional staining to assess M1 versus M2 macrophages was not evaluated in this study. No significant difference was detected in CD45RO, CD57, or granzyme B between the “mesenchymal” and “epithelial” lung adenocarcinoma samples (data not shown).
Elevation of multiple immunostimulatory molecules and IFNγ signals are found in “mesenchymal” lung adenocarcinoma
Although PD-L1 is generally viewed as an immune inhibitory molecule, its expression has been reported to reflect an ongoing antitumor immune response (9, 24, 25). To test whether the immunosuppressive profile we found in “mesenchymal” lung adenocarcinoma is associated with the infiltrating immune response, analysis of gene expression was performed on an extensive panel of immunostimulatory molecules. As shown in Fig. 4A and B, multiple immunostimulatory molecules, including CD80, CD86, OX40L, 4-1BB, ICOS, and CD127, were significantly elevated in “mesenchymal” as compared with “epithelial” lung adenocarcinoma. In addition, as shown in Fig. 4C and D, IFNγ, IFNγ-induced protein 10 (CXCL10), and an IFNγ−inducible enzyme indoleamine 2,3-dioxygenase (IDO) were also significantly elevated in “mesenchymal” lung adenocarcinoma as compared with the “epithelial” samples. Together, these data indicate an ongoing immune response in the tumor microenvironment of lung adenocarcinomas with “mesenchymal” phenotype. The elevated IFNγ released from ongoing host immune response most likely contributed to the increased PD-L1 expression and activation of other immunosuppressive molecules observed in “mesenchymal” lung adenocarcinoma. Once again, the results from the independent TCGA and PROSPECT datasets were highly consistent with one another.
EMT is not associated with tumor mutational burden, indicating an independent factor mediating an inflammatory tumor microenvironment
Genomic mutational burden has been posited as a marker for response to immunotherapy (26, 27). To investigate whether the profound inflammatory tumor microenvironment we observed in “mesenchymal” lung adenocarcinoma is confounded by increases in mutational burden, we studied the correlation of EMT and tumor mutation frequency. Using the TCGA dataset, for which we have DNA/RNA sequencing data and can assess mutational burden, there was no correlation between tumor mutational burden, as measured by nonsynonymous mutation per Mb, and EMT status (Fig. 5A and B). By contrast, and as expected, there was a significantly higher mutational burden found in NSCLC tumors with TP53 mutation as compared with tumors with wild-type TP53 (∼45-fold higher; Fig. 5C). Smoking is known to cause a high rate of mutations, including TP53, and high frequency of cytosine to adenine (C>A) nucleotide transversions has been identified as a surrogate marker for amount of smoking exposure (17, 28). We next tested the possible impact of EMT status on mutational burden, using nucleotide transversion as a surrogate measure for the amount of tobacco exposure. Patients with high frequency of transversion (C>A) carried significantly higher mutational burden than patients with low transversion, regardless of the status of TP53 (Fig. 5D). However, the presence of both TP53 mutation and high transversion frequency led to even higher mutational burden than the combination of transversion high and wild-type TP53, indicating the additional or synergistic effects between smoking and TP53 mutation. In contrast, no impact of EMT status on mutational burden was found, regardless of smoking exposure measured by transversion frequency (Fig. 5E). These data indicate that EMT is not a surrogate for mutational burden and is an independent mediator driving changes in the tumor immune microenvironment. Interestingly, although no correlation was seen between total mutational burden and EMT, we investigated the correlation between commonly identified individual mutations in lung adenocarcinoma and EMT, which revealed a significantly higher frequency of STK11 (P = 0.001) or KEAP1 (3.9 × 10−5) mutations in “epithelial” than “mesenchymal” lung adenocarcinomas (Supplementary Table S3). None of the other studied mutations, for example, TP53, KRAS, EGFR, or BRAF, were associated with a difference in the “epithelial” or “mesenchymal” groups (Supplementary Table S3).
B7-H3 is associated with poor OS and RFS in lung adenocarcinoma
We next analyzed the correlation between immune checkpoint molecules, EMT status and overall survival (OS) or recurrence-free survival (RFS). Neither PD-L1 nor EMT status was prognostic of OS or RFS. Among all the immune checkpoint molecules, only B7-H3 was found to negatively correlate with OS and RFS (Fig. 6 and Supplementary Fig. S4 using 1/2 or 1/3 cutoff, respectively). This may indicate a potential pathway to explore in therapeutic targeting of lung adenocarcinoma with immunotherapy.
Analysis of an independent dataset of advanced lung adenocarcinoma confirms the association of EMT with inflammatory tumor microenvironment
Our analyses using the two independent datasets, TCGA and PROSPECT, demonstrated consistent conclusions in surgically resected early-stage lung adenocarcinoma samples. To test whether the association between EMT and a profound inflammatory/immunosuppressive tumor microenvironment is limited to early-stage lung adenocarcinoma or more broadly applicable to all clinical stages, we next analyzed the gene expression profile of biopsy specimens from the BATTLE-1 study (16). This dataset consists only of core needle biopsy specimens from patients with advanced, treatment-refractory NSCLC. For the sake of consistency with our TCGA and PROSPECT analyses, we only included the lung adenocarcinoma samples for this study. As shown in Supplementary Fig. S5A, we observed a similar distinction in the immune-related marker profiles. Multiple immune checkpoint molecules, including PD-L1, PD-L2, TIM3, and CTLA-4, were significantly elevated in “mesenchymal” as compared with “epithelial” lung adenocarcinoma specimens (Supplementary Fig. S5B). Furthermore, multiple immunostimulatory molecules, including CD80, CD86, OX40L, 4-1BB, CD127, IFNγ-induced protein 10, and IFNγ-inducible enzyme IDO, were also significantly elevated in “mesenchymal” compared with the “epithelial” samples (Supplementary Fig. S5C and S5D). Taken together, our data demonstrate that the association of EMT with elevation of inflammatory signals and multiple immune checkpoints is a phenomenon broadly observed across all stages of lung adenocarcinoma.
Therapies targeting the immune checkpoint molecules PD-1 and PD-L1 have achieved remarkable clinical responses in multiple types of cancers, including NSCLC. Unfortunately, durable responses have been observed in only a subset of patients with chemotherapy-refractory metastatic disease. Identification of biomarkers that predict clinical benefit to immune-based approaches is needed. Accumulating data have suggested that immune checkpoint agents are most effective in patients in whom an endogenous immune response coexists with elevation of immune checkpoints (7, 27, 29, 30). However, biomarkers to identify this subgroup of patients who carry both endogenous immune response and elevation of immune checkpoints are essentially lacking. By using integrated gene expression analysis of three independent NSCLC datasets, we demonstrated that lung adenocarcinoma with a “mesenchymal” phenotype is associated with distinct tumor microenvironment changes. It is constituted by endogenous immune activation, such as elevation of immune costimulatory molecules, IFNγ, and CXCL10, along with simultaneous elevation of multiple immune checkpoint molecules, including elevated PD-1 and PD-L1, as compared with lung adenocarcinoma with an “epithelial” phenotype. Consistent with gene analysis, high expression of PD-L1 was confirmed by IHC and RPPA in “mesenchymal” lung adenocarcinoma. Enhanced tumor infiltration by CD4+Foxp3+ regulatory T cells and CD3+ T cells was also demonstrated in patients with “mesenchymal” lung adenocarcinoma in contrast to those with an “epithelial” phenotype.
EMT, a biologic program associated with loss of cell adhesion and increased invasive behavior has been established as a major mechanism for metastasis and drug resistance in several types of epithelial cancers, including NSCLC (10, 11, 13, 31–33). Studies have demonstrated that expression of transcription factors, such as snail or neu, can induce EMT and are associated with the activation of immunosuppressive cytokines and T-lymphocyte resistance in preclinical models of melanoma, pancreatic, and breast cancer (34–37). Our study provides evidence in NSCLC that EMT is associated with much broader inflammatory changes in the tumor microenvironment, with immune activation coexistent with elevation of multiple targetable immune checkpoint molecules. These data suggest a previously underrecognized role of tumor cell EMT. That is, EMT might accelerate cancer growth and metastasis not only by direct reprogramming of cancer cells, but also by reprogramming the immune response in the local tumor microenvironment. Furthermore, our data indicate that therapies targeting immune checkpoints might have a therapeutic impact on tumor metastases and drug resistance mediated via EMT.
PD-L1 expression on tumor cells or tumor-infiltrating immune cells has recently been studied as a potential single predictive biomarker for clinical activity to anti-PD-1/PD-L1–directed therapy (7, 38). Although studies have suggested that patients with overexpression of PD-L1 have improved clinical outcomes to anti-PD-1–directed therapy, some patients with overexpression of PD-L1 do not derive benefit, and other patients with low level PD-L1 expression demonstrate robust responses, indicating that PD-L1 is not an exclusionary biomarker and suggesting that a broader measure of the tumor microenvironment is needed. PD-L1 can be upregulated through different mechanisms, such as PTEN silencing, AKT activation, or inflammatory immune responses (24, 39, 40). IFNγ, secreted by activated T cells, is known to be the primary cytokine driving PD-L1 expression. The patients who clinically respond to ipilimumab and anti-PD-1 have been found to have tumors with preexisting immunity or inflamed tumors where activated T cells are present (7, 29, 41). Gene expressions suggestive of T-cell activation, such as elevated expression of IFNγ and CXCL10, were also found in responders as compared with nonresponders (7, 41). These findings suggest that upregulation of PD-L1 expression and subsequent development of immune suppression in the tumor microenvironment were most likely driven by an ongoing immune response as an adaptive immune escape mechanism (25, 29, 42, 43). Although the association of PD-L1 and increased tumor-infiltrating immune cells is not well understood, it was shown in murine models that elevation of PD-L1 in the tumor microenvironment is dependent on the presence of CD8+ T cells secreting IFNγ (24). These data emphasize that PD-L1 expression is likely not the cause, but rather the consequence, of increased tumor infiltration by immune cells. Our recent study demonstrated that the transcription factor ZEB1, a known EMT driver, regulates the miRNA-dependent expression of PD-L1 on tumor cells and enhances the tumor response to IFNγ (15). However, the redundant mechanisms underlying the association between EMT and the tumor immune microenvironment require further study. One possibility is that the cytokine milieu generated upon EMT increases tumor infiltration by immune cells. The selective clinical response to immunotherapies targeting the PD-1/PD-L1 pathway is most likely restricted to patients whose PD-L1 expression is associated with preexisting immunity where immune activation is augmented with immune checkpoint–blocking agents (7, 24, 25, 29,39, 40, 42, 43). Biomarker panels or a collective scoring system, such as the EMT score or the recently developed PD-L1 score that incorporates infiltrating immune cells to identify this subgroup of patients who have characteristic features of high levels of PD-L1 signals of T-cell activation and an inflammatory tumor microenvironment, will help us to select the best candidate patients for immunotherapies targeting the PD-1/PD-L1 pathway (7).
Consistent with these features, our data demonstrate a strong association between an ongoing immune response and the elevation of immune checkpoints in the tumor microenvironment when lung adenocarcinomas undergo EMT. Increased gene expressions of CD80, CD86, OX40L, 4-1BB, ICOS, CD127, IFNγ, and CXCL10 were found in “mesenchymal” compared with “epithelial” tumors. This ongoing immune response likely promoted the elevation of multiple immune checkpoint molecules. One could question why tumors were not rejected if active immune cell infiltration existed in the tumor microenvironment. Several lines of evidence have demonstrated that cytotoxic CD8 T cells isolated from tumor showed functional impairment despite initial proper T-cell activation (41, 44), indicating that immunosuppressive mechanisms inhibit the T-cell function in the tumor microenvironment. As revealed in our study, in addition to the PD-1:PD-L1/PD-L2 axis, several other immune checkpoint molecules were elevated, such as TIM-3, BTLA, B7-H3, and CTLA-4. This redundant suppressive tumor environment might imply the necessity of combinatorial strategies in the design of clinical trials in future.
Although expression of PD-L1 protein or mRNA has been associated with longer survival and better clinical outcome in patients with NSCLC, such a prognostic correlation was not found in our study (9). The use of variable techniques such as RNAscope assay versus Illumina sequencing, as well as different cutoffs and processing variability, or the types of samples utilized likely explain the differences. Interestingly, B7-H3, one of the immune checkpoint molecules commonly upregulated on NSCLC, was found significantly associated with poor OS and RFS and may represent a novel target for immunotherapy in NSCLC in the future. Although all three independent datasets shared the same trend in elevated immune checkpoints, such as PD-L1, PD-L2, CTLA-4, and TIM-3, along with elevated costimulatory molecules, including CD80, CD86, OX40L, 4-1BB, and CD127, some differences were noticed in comparison of late-stage with early-stage lung adenocarcinoma samples. For example, BTLA, LAG3, PD-1, ICOS, and IFNγ were significantly elevated in “mesenchymal” as compared with “epithelial” lung adenocarcinomas in patients with early, but not late stage of disease. Although we cannot completely explain this observation, it may indicate the potential presence of different tumor immune microenvironments in advanced stage lung cancers. Further investigation of the host immune system and tumor microenvironment in advanced versus early-stage lung cancer will better elucidate the optimal strategy for personalization of immunotherapies.
Recent studies in NSCLC and melanoma revealed that the clinical response to anti-PD-1 or CTLA-4 blockade correlates with genomic mutational burden (27, 45). A high mutational load is postulated to generate non-self-antigens that can be recognized by the immune system and trigger a mutation-reactive immune response. Although both smoking and TP53 mutation are highly associated with increased tumor mutational burden in lung adenocarcinoma, our data demonstrated no association between EMT and mutational burden (46). Our study therefore indicates that EMT is an independent mediator driving inflammatory and immunosuppressive tumor microenvironment changes in addition to total mutational burden and suggests that patients with tumors carrying both high EMT score and high mutation burden are more likely to benefit from immune checkpoint therapy than those who have either alone. Interestingly, investigating the correlation between common individual mutations and EMT revealed a significantly higher frequency of STK11 and KEAP1 mutations in “epithelial” lung adenocarcinomas. This is consistent with our recent study showing that KRAS-mutant lung adenocarcinomas with concurrent STK11 or KEAP1 mutations appeared largely “immune-inert” and expressed low level of immune markers (47). Further study to investigate the underlying mechanisms behind the association of specific driver mutations and immune suppression will facilitate the implementation of personalized therapy.
In summary, our data demonstrate a strong association between EMT and an inflammatory tumor microenvironment with expression of multiple immune checkpoint molecules and immune activation. Along with using the difference between “epithelial” and “mesenchymal” groups to explore the mechanistic biology driving treatment response and resistance, further validation of potential utility of using EMT as a predictive biomarker to select patients for immune checkpoint blockade and other immunotherapies in NSCLC is needed. Because implementing an mRNA-based gene signature will be clinically challenging, other simplified testing schema will need to be devised. We are currently exploring the use of simplified methodologies that could be implemented in a clinical testing environment.
Disclosure of Potential Conflicts of Interest
J.V. Heymach reports receiving other commercial research support from AstraZeneca, Bayer, and GlaxoSmithKline; and is a consultant/advisory board member for AstraZeneca, Boerhinger Ingelheim, Exelixis, Genentech, GlaxoSmithKline, Lilly, Novartis, and Synta. No potential conflicts of interest were disclosed by the other authors.
Conception and design: Y. Lou, P. Hwu, J.V. Heymach, D.L. Gibbons
Development of methodology: L.A. Byers, J.C. Rodriguez, J.V. Heymach, D.L. Gibbons
Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): Y. Lou, E.R.P. Cuentas, W.L. Denning, L. Chen, Y.H. Fan, L.A. Byers, V.A. Papadimitrakopoulou, C. Behrens, J.C. Rodriguez, I.I. Wistuba, J.V. Heymach
Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): Y. Lou, E.R.P. Cuentas, L. Chen, L.A. Byers, J. Wang, V.A. Papadimitrakopoulou, P. Hwu, I.I. Wistuba, J.V. Heymach, L. Diao, C. Behrens
Writing, review, and/or revision of the manuscript: Y. Lou, E.R.P. Cuentas, W.L. Denning, L. Chen, L.A. Byers, V.A. Papadimitrakopoulou, J.C. Rodriguez, P. Hwu, J.V. Heymach, D.L. Gibbons
Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): Y. Lou, E.R.P. Cuentas, L. Chen, J.V. Heymach
Study supervision: J. Wang, P. Hwu, J.V. Heymach, D.L. Gibbons
This study is supported by the Department of Defense–supported PROSPECT Grant and NCI-funded Lung SPORE P50 CA070907 (to J.V. Heymach and I.I. Wistuba), Conquer Cancer Foundation of ASCO Young Investigator Award 2014 (to Y. Lou), NIH–T32 Research Training in Academic Medical Oncology (to Y. Lou), NIH R01 grant R01CA1668484 (to J.V. Heymach), LUNGevity Foundation Research Award (to D.L. Gibbons and L.A. Byers), Uniting Against Lung Cancer/Lung Cancer Research Foundation (to D.L. Gibbons), Rexanna's Foundation for Fighting Lung Cancer (to D.L. Gibbons), The Stading Lung Cancer Research Fund (to J.V. Heymach), MD Anderson Cancer Center Physician Scientist Award (to D.L. Gibbons), K08-CA151651 (D.L. Gibbons), and CPRIT RP150405 (D.L. Gibbons). D.L. Gibbons is an R. Lee Clark Fellow of the University of Texas MD Anderson Cancer Center, supported by the Jeane F. Shelby Scholarship Fund. This work was also supported by the generous philanthropic contributions to The University of Texas MD Anderson Lung Cancer Moon Shots Program.
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.