Skip to main content

Prediction of response to pemetrexed in non-small-cell lung cancer with immunohistochemical phenotyping based on gene expression profiles



Palliative pemetrexed-based chemotherapy remains a standard of care treatment for the majority of patients with advanced non-squamous non-small-cell lung cancer (NSCLC). Currently, no predictive markers for pemetrexed treatment are available.


Resected tumour samples from pemetrexed-naïve NSCLC patients were collected. Gene expression profiling with respect to predicted sensitivity to pemetrexed classified predicted responders (60%) and non-responders (40%) based on differentially expressed genes encoding for pemetrexed target enzymes. Genes showing a strong correlation with these target genes were selected for measurement of corresponding protein expressions by immunohistochemical (IHC) staining. A semi-quantitative IHC scoring method was applied to construct a prediction model for response to pemetrexed. A retrospective cohort of patients with advanced NSCLC treated with first-line pemetrexed-based chemotherapy was used for external validation.


From ninety-one patients resected tumour samples were collected. The majority of patients had early or locally advanced NSCLC (96.3%). Gene expression profiling revealed five markers, which mRNA levels strongly correlated to pemetrexed target genes mRNA levels: TPX2, CPA3, EZH2, MCM2 and TOP2A. Of 63 (69%) patients IHC staining scores of these markers were obtained, which significantly differed between predicted non-responders and responders (P < 0.05). The optimized prediction model included EZH2 (OR = 0.56, 95% CI 0.35–0.90) and TPX2 (OR = 0.55, 95% CI 0.30–1.01). The model had a sensitivity of 86.8%, specificity of 63.6% and showed a good ability to distinct between responders and non-responders (C-index 0.86).

In the external study population (N = 23) the majority of patients had metastatic NSCLC (95.7%). Partial response (PR) was established in 26.1%. The sensitivity decreased drastically to 33.3%, with a specificity of 82.4% and a C-index of 0.73.


Using external validation this prediction model with IHC staining of target enzyme correlated markers showed a good discrimination, but lacked sensitivity. The role of IHC markers as response predictors for pemetrexed in clinical practice remains questionable.

Peer Review reports


In the management of advanced non-small cell lung cancer (NSCLC) systemic treatment options are rapidly expanding with the increasing use of molecular targeted agents and immunotherapy [1,2,3,4]. One of the most important therapeutic advances has been the identification of predictive molecular markers to guide patient selection for frontline treatment with these agents, like sensitizing mutations within the EGFR gene to EGFR tyrosine kinase inhibitors and protein PD-L1 overexpression to anti PD-(L)1 checkpoint inhibitors [1, 2, 5]. Despite the changing treatment landscape with increasing use of molecular targeted agents and immunotherapy, pemetrexed-based chemotherapy is still widely used as standard treatment in patients with advanced non-squamous non-small-cell lung cancer [6, 7]. Unfortunately, to date useful biomarkers predicting response to this treatment regimen are lacking.

Pemetrexed treatment shows a variable clinical efficacy apparently dependent on the histologic subtype of lung cancer. Clinical trials demonstrated efficacy of pemetrexed in nonsquamous NSCLC while efficacy was worse in squamous NSCLC and small-cell lung cancer [8,9,10]. However, tumour response to pemetrexed also differs significantly between patients with similar histology [9, 11]. In patients treated with pemetrexed monotherapy, the response rate to pemetrexed was evidently different between histological subtypes but low in both patients with squamous and non-squamous NSCLC (2.8% vs. 11.5%) [10]. In this study pemetrexed was administered as second-line treatment and patients with poor ECOG performance score were included. Although the response rate to pemetrexed was significantly higher in patients with non-squamous versus squamous NSCLC in the first-line pivotal trial [9], still more than 20% of patients with squamous NSCLC had a response to pemetrexed while the response rate was merely ~ 30% in patients with non-squamous histology [11]. These findings highlight the need for predictive molecular markers for pemetrexed-based treatment.

The main determinant of pemetrexed responsiveness is thought to be the level of expression of thymidylate synthase (TS), the primary intracellular target enzyme of pemetrexed [12, 13]. Overexpression of TS mRNA has been correlated with reduced sensitivity to pemetrexed in vitro, [14,15,16,17] and with worse clinical outcomes in patients treated with pemetrexed [18]. Moreover, the abundance of TS expression is higher in squamous cell NSCLC than in other histologic subtypes, [19, 20] which constitutes the biological hypothesis behind the superior efficacy of pemetrexed in non-squamous NSCLC. However, in clinical practice the relationship between protein expression levels of TS, measured by immunohistochemistry (IHC) methods, and the clinical efficacy of pemetrexed remains controversial [21,22,23,24,25].

Our study group earlier presented an approach to implement a more refined molecular classification of NSCLC subtypes based on gene expression profiles independent of histology [26]. Furthermore, response to pemetrexed was predicted based on expression of genes encoding different pemetrexed target enzymes including but not limited to TS, and expression signatures of correlated genes were identified. In the current study, we explore whether these differential gene expression profiles between responders and non-responders can be used to define a prediction model based on IHC scores of selected molecular markers.


Training cohort

Tumour samples from pemetrexed-naïve, NSCLC patients who had undergone curative surgical resection at the Erasmus Medical Center (Rotterdam) between 1992 and 2004 were used. A detailed description of tissue collection, microarray preparation and data processing, the derivation of a gene-expression based predictive algorithm for tumour response and the identification of pemetrexed resistance-associated genes has been previously described [26]. In short, the predictive algorithm predicted tumour response based on the expression difference between internal reference genes and pemetrexed target genes TS, dihydrofolate reductase (DHFR) and glycinamide ribonucleotide formyltransferase (GARFT). Using percentile-rank based target gene expression levels relative to the internal reference genes, patients were stratified as predicted responders (±60%) and non-responders (±40%). Subsequently, significance analysis of microarray identified differentially expressed genes between these classified groups [27]. A minimized signature containing 25 genes performed optimally in predicting pemetrexed response (Additional file 1: Table S1). For the current study, we selected molecular markers from this signature if they showed a strong correlation with the gene expression of TS and if IHC stainings for these markers were commercially available. Written informed consent was obtained from all these patients. The study was conducted in accordance with the REMARK guidelines [28].

Validation cohort

In order to externally validate the model, we obtained formalin-fixed paraffin-embedded pre-treatment biopsies of a retrospective cohort of patients newly diagnosed with advanced stage (IIIB/IV) NSCLC in a large teaching hospital (Amphia hospital, Breda, the Netherlands) between January 2007 and December 2010. Patients were eligible for enrolment if they had received ≥2 cycles of platinum-combined pemetrexed chemotherapy as first-line treatment. Medical charts and radiological imaging data were reviewed to collect information regarding socio demographics, tumour histology, ECOG performance status, treatment and observed tumour response (RECIST 1.1). Patients with early stage (IA-IIB) or locally advanced (IIIA) disease, (neo)adjuvant chemotherapy, combination treatment with bevacizumab and without tissue samples from primary tumour or (lymph node) metastasis were excluded.

Tissue microarray analysis and immunohistochemistry

The tissue microarrays (TMAs) were composed of 68 of the 91 tumour tissues, in triplicate, from the Erasmus MC patient cohort used for the expression microarray analyses. TMA blocks containing 0.6 mm cores of formalin-fixed paraffin-embedded tumours were cut and antigen retrieval was performed by a 20-min incubation at 95 °C using Trisethylenediaminetetraacetic acid buffer (Klinipath, Duiven, The Netherlands). Subsequently, TMAs were stained with primary antihuman antibodies of the selected candidate markers: EZH2, TOP2A, TPX2, MCM2 and CPA3 (Additional file 2: Table S2). For each TMA multicontrol stainings were performed using a combination of tissues (liver, pancreas, tonsil, colon and appendix). The slides were stained and processed in the Ventana Benchmark ULTRA strainers, using DAB as substrate and Hematoxylin as counterstain. Tumour tissues from the validation group were equally handled, except that the tumour samples were cut in 0.4 mm instead of 0.6 mm cores for TMA blocks.

Immunohistochemical staining score

A semi-quantitative scoring method was applied to classify the intensity and quantity of IHC staining of candidate markers. The quantity score was defined as: 1: 0–30%; 2: 30–60%; 3: 60–100%. The intensity score was defined as: 0 (negative), no appreciable staining in the tumour cells; 1 (weak), barely detectable cytoplasmic/membranous or nuclear staining of tumour cells; 2 (moderate), readily appreciable staining of tumour cytoplasm/nucleus; 3 (strong), strong staining obscuring nucleus/cytoplasm of tumour cells. Multiplying quantity and intensity score yielded a total score with a range between 0 and 9. TMAs from the training group were evaluated and scored for protein expression simultaneously by K.B. and J.P.H. Samples were individually discussed until consensus was reached. For the validation group, TMA evaluation and protein staining quantification were performed independently by K.B. and S.V.

Statistical analyses

Sociodemographic and disease- and treatment-related variables were described for all patients who were included in this study and were compared between the training and validation group. We used the independent-samples t-test and the χ2-test or Fisher’s exact test for continuous and categorical variables respectively. Degree of agreement on quantity and intensity scores of the different IHC stainings was evaluated using weighted linear Cohen’s kappa scores (휅) in the validation group. Degree of agreement was determined according to widely used scale described by Landis and Koch [29]. As IHC staining scores from the training group were obtained by discussion, and thus not independently, no interobserver agreement could be calculated.

Pearson correlation coefficients (ρ) were calculated between gene expression of the candidate markers and TS in the training group and subsequently between the gene expression of those markers and their associated protein expressions. Using the described prediction algorithm of response to pemetrexed, patients from the training group were divided in predicted responders and non-responders. In the training group, we compared gene and protein expression from selected molecular markers between predicted responders and non-responders.

Logistic regression with dependent variable predicted tumour response by gene expression signature was applied to the training cohort to build a prediction model with the IHC staining scores of selected molecular markers as independent variables. Optimized model derivation was performed using purposeful selection by stepwise in- and exclusion of molecular markers [30]. Univariable logistic regression identified molecular markers associated with predicted tumour response. We specified a priori that molecular markers with p < 0.2 on univariable analysis would be candidate variables for multivariable logistic regression model. In the iterative process of variable selection, variables were removed from the model if they were non-significant and not a confounder. We used backward selection with a p-value < 0.05 and/or change of effect size of (an)other included variable(s) > 20% to remain in the model. The fit of a reduced model versus full model was compared with the likelihood ratio test (LRT), following a chi-squared distribution. Subsequently, the model was externally validated in the validation cohort. In both cohorts, the model performance of the derived model was assessed by examining the predictive classification accuracy and discriminatory ability (C-index). A C-index of 1.0 would indicate perfect discrimination, whereas a C-index of 0.5 indicates total absence of discrimination. All statistical analyses were performed with the use of SPSS, version 24.0 (IBM Corporation, Armonk, NY).


The selection of patients in both training and validation cohort is depicted in Fig. 1. Of the 91 surgically resected samples of the primary tumour, 68 (74.7%) samples were suitable for further processing into TMAs. Since samples of five patients could not be used due to insufficient TMA material, the training cohort ultimately consisted of 63 patients whose samples were prepared with additional IHC stainings. For validation, 44 of 142 (31%) patients who received pemetrexed had advanced stage NSCLC treated with ≥2 cycles first-line treatment with platinum-combined pemetrexed, excluding combinations with bevacizumab. Of these patients, 18 (40.9%) were excluded from further analysis because their diagnosis was cytology-based and no histologic biopsy was obtained. Histology samples were retrieved from primary tumour (N = 13) or lymph nodes or distant metastases (N = 13). Three patients had insufficient tumour material available for additional IHC staining and therefore the validation cohort consisted of 23 patients.

Fig. 1
figure 1

Flowchart of training and validation population. Abbreviations: NSCLC, non-small-cell lung cancer; TMA, tissue micro array; IHC, immunohistochemical

Patient characteristics

Patient and treatment characteristics are shown in Table 1. In the training cohort 46 patients (73.0%) were male compared to 10 (43.5%) in the validation cohort. Eighteen patients (28.6%) had squamous NSCLC in the training cohort while in the validation cohort only patients were included with non-squamous histology. The majority of patients (95.2%) in the training group had early stage NSCLC opposed to all patients with advanced disease stage in the validation group. No data were available with regard to the ECOG performance score of the patients in the training cohort. In the validation cohort, 21.7% of the patients had a performance score of 2. Corresponding to the differences in disease stage between the cohorts, all patients in the training cohort underwent surgical resection in contrast to palliative pemetrexed-based chemotherapy in the validation cohort. Patients in the training cohort had worse overall survival compared to the validation cohort (4.5 months vs. 28 months).

Table 1 Characteristics of patients in the training population and the validation population

Selection molecular markers

Of the 25-genes containing optimized gene expression signature predicting response to pemetrexed in the training group, five molecular markers were selected based on their correlation with the gene expression level of target genes (TS, DHFR, GARFT) and the commercial availability of corresponding IHC stainings: Enhancer of zeste homolog 2 (EZH2), Topoisomerase II Alpha (TOP2A), Microtubule Nucleation Factor (TPX2), Carboxypeptidase A3 (CPA3) and Minichromosome Maintenance Complex Component 2 (MCM2). All markers showed a positive correlation to the mRNA level of TS (EZH2, ρ = 0.732; MCM2, ρ = 0.804; TOP2A, ρ = 0.814; TPX2, ρ = 0.825), except for CPA3 which was negatively correlated (ρ = − 0.467) (Additional file 4: Figure S1). The correlation of gene mRNA level with their corresponding IHC staining score had a range between 0.303 (CPA3) and 0.578 (EZH2) (Additional file 5: Figure S2).

The IHC stainings of the same markers were applied to the TMAs of the samples of patients in the validation cohort. The strength of agreement between the observers with regard to the intensity score ranged between 휅 = 0.515 (CPA3) and 휅 = 1 (MCM2), and with regard to the quantity score between 휅 = 0.547 (TPX2) and 휅 = 0.851 (CPA3). Weighted kappa values of IHC staining scores are outlined in Additional file 3: Table S3.

Of all selected markers, both mRNA levels and IHC staining scores were significantly higher in predicted non-responders than responders in the training group, except for CPA which mRNA level and IHC staining score was significantly lower in non-responders compared to responders. These results are depicted in Fig. 2.

Fig. 2
figure 2

Boxplots of gene and protein expression levels of selected markers in predicted responders and non-responders to pemetrexed in the training group. Boxplots (Tukey) representing medians and interquartile ranges of relative mRNA expression and IHC staining scores respectively. Whiskers represent minimum and maximum 1.5 interquartile range and dots are outliers. Abbreviations: IHC, immunohistochemical; EZH2, Enhancer of zeste homolog; TOP2A, Topoisomerase II; TPX2, Microtubule Nucleation Factor; CPA3, Carboxypeptidase A3; MCM2, Minichromosome Maintenance Complex Component 2. Abbreviations: IHC, immunohistochemical; EZH2, Enhancer of zeste homolog; TOP2A, Topoisomerase II Alpha; TPX2, Microtubule Nucleation Factor; CPA3, Carboxypeptidase A3; MCM2, Minichromosome Maintenance Complex Component 2

Model derivation

The model coefficients and odds ratios (OR) with corresponding 95% confidence intervals (CI) of the prediction model with dependent variable tumour response to pemetrexed using univariable and multivariable logistic regression analyses are presented in Table 2. Univariable analyses of the relationship between the IHC staining scores of the selected markers and the gene expression based predicted tumour response to pemetrexed were performed using the training cohort. Higher IHC staining scores of all markers were significantly associated with worse predicted tumour response, except for CPA3 which repeatedly showed the reverse association compared to the other markers. Using multivariable analysis, only a higher IHC staining score of EZH2 was significantly related with a worse predicted tumour response to pemetrexed (OR 0.56, 95% CI 0.35–0.90; P = 0.015). The staining scores of all other markers failed to demonstrate a significant association on multivariable analysis. Although IHC staining score of TPX2 was not significantly associated with tumour response in the multivariable model (OR 0.55, 95% CI 0.30–1.01; P = 0.054), this variable was still included in the final optimized prediction model. Removal of this variable led to a significantly reduced model fit (P < 0.001).

Table 2 Prediction model derivation to predict tumor response using IHC staining scores of selected molecular markers in training group (N = 63)

Model performance and validation

In Table 3 the different test characteristics describing the performance of the model in the training and validation cohort are shown. In the training cohort, 38 patients were predicted responders (63.3%) by gene expression profiling and 86.8% (33 of 38) were correctly classified responders by the prediction model (sensitivity 86.8%, 95% CI 71.9–95.6). Fourteen patients were predicted non-responders (36.7%) by gene expression profiling while 63.6% (14 of 22) were correctly classified as non-responders by the model (specificity 63.6%, 95% CI 40.7–82.8).

Table 3 Conditional and post-test probability performance of the IHC based prediction model in the training and validation cohort

In the validation cohort, the same classification by the prediction model was applied, however an actual tumour response was obtained as these patients were treated with pemetrexed. The response rate was 26.1% and therefore the prevalence of response was substantially lower than in the training cohort. Of the six patients who experienced a partial response, two patients were correctly classified by the prediction model resulting in a sensitivity of 33.3% (95% CI 4.3–77.7). Herewith, the sensitivity in this cohort is significantly worse than the sensitivity in the training cohort (Fisher’s exact test, P = 0.011). The positive predictive value (PPV) also decreased substantially. If we classified both patients with a partial response and stable disease as responders (56.5%), performance characteristics of the model declined dramatically (Table 3).

The ROC curve showed a C-index of 0.86 (95% CI 0.77–0.96) in the training cohort, representing a good discriminatory performance. The C-index decreased to 0.73 (95% CI 0.52–0.93) if the prediction model was applied to the validation cohort (Fig. 3).

Fig. 3
figure 3

Receiver operating characteristic curve showing model performance of two-protein expression prediction model in training and validation population. Diagonal line reflects total absence of discrimination (AUC = 0.5). Abbreviations: AUC, area under the curve


Currently, the profit of and need for molecular markers to select therapy for individual patients is increasingly recognized. The last several years, treatment of advanced NSCLC has become obviously more complex and therefore tools to choose therapies that are most likely to benefit patients are required. Indeed the registration of new therapeutic agents as frontline therapy is accompanied by selective markers. For patients with EGFR mutation-positive, ALK rearrangement-positive or ROS1 rearrangement positive tumours first-line molecular targeted tyrosine kinase inhibitors are recommended, while patients with high PD-L1 expression in the tumour are suitable for pembrolizumab first-line treatment [6]. Despite the wide implementation of molecular markers, the administration of chemotherapy in NSCLC patients is still solely based on histology even though its capacity to predict response has been proved to be suboptimal.

In an earlier study, we demonstrated an approach to predict response to pemetrexed based on the use of gene expression profiles [26]. Samples of different histological subtypes including squamous NSCLC were used. Prediction of response to pemetrexed based on gene expression profiling of its target enzymes, failed to show the expected disadvantage for the squamous cell histological subtype and these results therefore challenge the restricted use of pemetrexed in non-squamous NSCLC. In the present study, we developed a prediction model using immunohistochemistry scores of selected candidate genes (EZH2 and TPX2) from the gene expression signature predicting pemetrexed response. In the training group, the use of the model resulted in good performance characteristics of the model in the sense of a high sensitivity and PPV. The model also showed the ability to discriminate well between responders and non-responders. Unfortunately, the results of the obtained model could not be validated when applied to an external cohort of patients treated with pemetrexed.

Although it is hypothesized that superior efficacy of pemetrexed in non-squamous over squamous NSCLC is related to the level of TS expression, multiple clinical studies failed to demonstrate the association between TS protein expression and clinical outcomes [22, 25]. This can be ascribed firstly to the fact that current semi-quantitative IHC methods might lack sensitivity required to measure protein expression of TS opposed to quantitative analysis of mRNA expression. And secondly, the biological significance of TS for pemetrexed responsiveness might be of less importance than other molecular processes, e.g. gene expression or amplification of other (target) genes. To overcome these limitations in the present study, we carefully selected markers related to mRNA gene expression of TS but also to other target enzymes DHFR and GARFT.

In the literature, both EZH2 and TPX2 have been previously linked to survival in NSCLC. TPX2 is involved in key steps during mitotic events and increased expression has been associated with poor overall survival in NSCLC [31, 32]. EZH2 epigenetically silences multiple genes involved in cell differentiation, growth and invasion. It is often overexpressed in NSCLC promoting cancer progression and a more aggressive tumour behaviour [33, 34]. Downregulation of EZH2 has been associated with higher expression of oestrogen receptor and increased sensitivity to tamoxifen in advanced breast cancer patients [35]. Similarly, one can speculate that EZH2 might change the expression of genes related to responsiveness of pemetrexed through its ability to silence other genes. Although lack of a control arm precludes discrimination between a prognostic or predictive factor, we purposefully focused on radiological response rather than survival.

It is crucial to predict treatment effects for individual patients in order to avoid unnecessary toxicities and to offer alternative treatment options, with few false positives [36]. The clinical value of the derived model is probably limited as the sensitivity was poor in the validation group. Moreover, the low PPV makes the classifier not useful for clinical decision making, as many patients who are predicted responders will then actually undergo potentially harmful treatment with low chance of tumour response. The failure of the prediction model to adequately perform in the external patient population might be ascribed to an insufficient sensitivity of used IHC assays to measure significant differences in protein expression or a discrepancy between protein expression and protein activity. Additionally, spatiotemporal heterogeneity might have led to different intrinsic tumour properties in the validation group as these patients had advanced disease and in half of the cases tumour samples were obtained from lymph node or distant metastases. Finally, other factors might influence pemetrexed activity such as cell transport and intracellular formation of polyglutamate metabolites [37]. Our study was limited in the number of patients included, especially in the external validation.

cohort. We recognize that the differences in histology between the training and validation cohort is a major shortcoming. It was impossible to include patients with non-squamous NSCLC to the validation cohort, as selection was treatment-based and pemetrexed is only recommended in patients with nonsquamous NSCLC. For ALK and ROS1-rearrangement positive adenocarcinoma patients might experience more benefit to pemetrexed-based chemotherapy [38, 39], molecular characteristics would have been desirable. Unfortunately, those data were not available in our cohort, but the high number of smokers profoundly reduces the chance of rearrangements. Although response rates were in accordance with the literature, patients in the validation cohort experienced a very poor median overall survival of only 4.5 months. This can be partially explained by the presence of a substantial group of patients (> 20%) with a poor ECOG performance score and suboptimal treatment with carboplatin instead of cisplatin combination. Whether these patients appropriately represent the population with advanced NSCLC is therefore highly questionable and we cannot exclude their genetic profile to be different. However, given the results we do not expect that expanding the number of samples will lead to a clinically useful biomarker.


There remains an unmet need to identify biomarkers to select patients for standard pemetrexed-based treatment. Prediction of pemetrexed responsiveness with IHC stainings of markers correlated to TS and other target enzymes could not be validated using external validation. Future research focusing on metabolomics, pharmacokinetics and pharmacogenetics might offer new insights into tailoring therapy. Until a well-validated biomarker is identified, histology should remain the standard to select advanced NSCLC patients eligible for treatment with pemetrexed.



Confidence interval


Carboxypeptidase A3


Dihydrofolate reductase


Enhancer of zeste homolog 2


Glycinamide ribonucleotide formyltransferase




Likelihood ratio test


Minichromosome Maintenance Complex Component 2


Non-small-cell lung cancer


Odds ratio


Positive predictive value


Tissue microarrays


Topoisomerase II Alpha


Microtubule Nucleation Factor


Thymidylate synthase


  1. Mok TS, Wu Y-L, Thongprasert S, Yang C-H, Chu D-T, Saijo N, et al. Gefitinib or carboplatin–paclitaxel in pulmonary adenocarcinoma. N Engl J Med. 2009;361:947–57.

    Article  CAS  PubMed  Google Scholar 

  2. Shaw AT, Kim D-W, Nakagawa K, Seto T, Crinó L, Ahn M-J, et al. Crizotinib versus chemotherapy in advanced ALK -positive lung Cancer. N Engl J Med. 2013;368:2385–94.

    Article  CAS  PubMed  Google Scholar 

  3. Reck M, Kaiser R, Mellemgaard A, Douillard J-Y, Orlov S, Krzakowski M, et al. Docetaxel plus nintedanib versus docetaxel plus placebo in patients with previously treated non-small-cell lung cancer (LUME-lung 1): a phase 3, double-blind, randomised controlled trial. Lancet Oncol. 2014;15:143–55.

    Article  CAS  PubMed  Google Scholar 

  4. Borghaei H, Paz-Ares L, Horn L, Spigel DR, Steins M, Ready NE, et al. Nivolumab versus docetaxel in advanced nonsquamous non–small-cell lung Cancer. N Engl J Med. 2015;373:1627–39.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Reck M, Rodríguez-Abreu D, Robinson AG, Hui R, Csőszi T, Fülöp A, et al. Pembrolizumab versus chemotherapy for PD-L1–positive non–small-cell lung Cancer. N Engl J Med. 2016;375:1823–33.

    Article  CAS  PubMed  Google Scholar 

  6. Hanna N, Johnson D, Temin S, Baker S, Brahmer J, Ellis PM, et al. Systemic therapy for stage IV non-small-cell lung Cancer: American Society of Clinical Oncology clinical practice guideline update. J Clin Oncol. 2017;JCO2017746065.

    Article  CAS  PubMed  Google Scholar 

  7. Novello S, Barlesi F, Califano R, Cufer T, Ekman S, Levra MG, et al. Metastatic non-small-cell lung cancer: ESMO clinical practice guidelines for diagnosis, treatment and follow-up†. Ann Oncol. 2016;27(suppl_5):v1–27.

    Article  CAS  PubMed  Google Scholar 

  8. M a S, Smit EF, Lorigan P, Konduri K, Reck M, Szczesna A, et al. Phase III study of pemetrexed plus carboplatin compared with etoposide plus carboplatin in chemotherapy-naive patients with extensive-stage small-cell lung cancer. J Clin Oncol. 2009;27:4787–92.

    Article  CAS  Google Scholar 

  9. Scagliotti GV, Parikh P, von Pawel J, Biesma B, Vansteenkiste J, Manegold C, et al. Phase III study comparing cisplatin plus gemcitabine with cisplatin plus pemetrexed in chemotherapy-naive patients with advanced-stage non-small-cell lung cancer. J Clin Oncol. 2008;26:3543–51.

    Article  CAS  PubMed  Google Scholar 

  10. Hanna N, F a S, Fossella FV, Pereira JR, De Marinis F, von Pawel J, et al. Randomized phase III trial of pemetrexed versus docetaxel in patients with non-small-cell lung cancer previously treated with chemotherapy. J Clin Oncol. 2004;22:1589–97.

    Article  CAS  PubMed  Google Scholar 

  11. Scagliotti G, Hanna N, Fossella F, Sugarman K, Blatter J, Peterson P, et al. The differential efficacy of pemetrexed according to NSCLC histology: a review of two phase III studies. Oncologist. 2009;14:253–63.

    Article  CAS  PubMed  Google Scholar 

  12. Shih C, Chen VJ, Gossett LS, Gates SB, MacKellar WC, Habeck LL, et al. LY231514, a pyrrolo [2,3-d]pyrimidine-based antifolate that inhibits multiple folate-requiring enzymes. Cancer Res. 1997;57:1116–23 Accessed 2 Oct 2017.

    CAS  PubMed  Google Scholar 

  13. Chattopadhyay S, Moran RG, Goldman ID. Pemetrexed: biochemical and cellular pharmacology, mechanisms, and clinical applications. Mol Cancer Ther. 2007;6:404–17.

    Article  CAS  PubMed  Google Scholar 

  14. Hanauske A-R, Eismann U, Oberschmidt O, Pospisil H, Hoffmann S, Hanauske-Abel H, et al. In vitro chemosensitivity of freshly explanted tumor cells to pemetrexed is correlated with target gene expression. Investig New Drugs. 2007;25:417–23.

    Article  CAS  Google Scholar 

  15. Sigmond J, Backus HHJ, Wouters D, Temmink OH, Jansen G, Peters GJ. Induction of resistance to the multitargeted antifolate Pemetrexed (ALIMTA) in WiDr human colon cancer cells is associated with thymidylate synthase overexpression. Biochem Pharmacol. 2003;66:431–8.

    Article  CAS  PubMed  Google Scholar 

  16. Takezawa K, Okamoto I, Okamoto W, Takeda M, Sakai K, Tsukioka S, et al. Thymidylate synthase as a determinant of pemetrexed sensitivity in non-small cell lung cancer. Br J Cancer. 2011;104:1594–601.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Ozasa H, Oguri T, Uemura T, Miyazaki M, Maeno K, Sato S, et al. Significance of thymidylate synthase for resistance to pemetrexed in lung cancer. Cancer Sci. 2010;101:161–6.

    Article  CAS  PubMed  Google Scholar 

  18. Liu Y, Yin T-J, Zhou R, Zhou S, Fan L, Zhang R-G. Expression of thymidylate synthase predicts clinical outcomes of pemetrexed-containing chemotherapy for non-small-cell lung cancer: a systemic review and meta-analysis. Cancer Chemother Pharmacol. 2013;72:1125–32.

    Article  CAS  PubMed  Google Scholar 

  19. Ceppi P, Volante M, Saviozzi S, Rapa I, Novello S, Cambieri A, et al. Squamous cell carcinoma of the lung compared with other histotypes shows higher messenger RNA and protein levels for thymidylate synthase. Cancer. 2006;107:1589–96.

    Article  CAS  PubMed  Google Scholar 

  20. Monica V, Scagliotti GV, Ceppi P, Righi L, Cambieri A, Lo Iacono M, et al. Differential thymidylate synthase expression in different variants of large-cell carcinoma of the lung. Clin Cancer Res. 2009;15:7547–52.

    Article  CAS  PubMed  Google Scholar 

  21. Christoph DC, Asuncion BR, Hassan B, Tran C, Maltzman JD, O’Shannessy DJ, et al. Significance of folate receptor alpha and thymidylate synthase protein expression in patients with non–small-cell lung Cancer treated with Pemetrexed. J Thorac Oncol. 2013;8:19–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Wynes MW, Konopa K, Singh S, Reyna-Asuncion B, Ranger-Moore J, Sternau A, et al. Thymidylate synthase protein expression by IHC and gene copy number by SISH correlate and show great variability in non–small cell lung Cancer. J Thorac Oncol. 2012;7:982–92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Sun J-M, Han J, Ahn JS, Park K, Ahn M-J. Significance of thymidylate synthase and thyroid transcription factor 1 expression in patients with nonsquamous non-small cell lung cancer treated with pemetrexed-based chemotherapy. J Thorac Oncol. 2011;6:1392–9.

    Article  PubMed  Google Scholar 

  24. Nicolson MC, Fennell DA, Ferry D, O’Byrne K, Shah R, Potter V, et al. Thymidylate synthase expression and outcome of patients receiving Pemetrexed for advanced nonsquamous non–small-cell lung Cancer in a prospective blinded assessment phase II clinical trial. J Thorac Oncol. 2013;8:930–9.

    Article  CAS  PubMed  Google Scholar 

  25. Shimizu T, Nakanishi Y, Nakagawa Y, Tsujino I, Takahashi N, Nemoto N, et al. Association between expression of thymidylate synthase, dihydrofolate reductase, and glycinamide ribonucleotide formyltransferase and efficacy of pemetrexed in advanced non-small cell lung cancer. Anticancer Res. 2012;32:4589–96 Accessed 18 Sep 2017.

    CAS  PubMed  Google Scholar 

  26. Hou J, Lambers M, den Hamer B, M a d B, Hoogsteden HC, Grosveld F, et al. Expression profiling-based subtyping identifies novel non-small cell lung cancer subgroups and implicates putative resistance to pemetrexed therapy. J Thorac Oncol. 2012;7:105–14.

    Article  CAS  PubMed  Google Scholar 

  27. Tusher VG, Tibshirani R, Chu G. Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci U S A. 2001;98:5116–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. McShane LM, Altman DG, Sauerbrei W, Taube SE, Gion M, Clark GM. Reporting recommendations for tumor marker prognostic studies (REMARK). JNCI J Natl Cancer Inst. 2005;97:1180–4.

    Article  CAS  PubMed  Google Scholar 

  29. Landis JR, Koch GG. The measurement of observer agreement for categorical data. Biometrics. 1977;33:159–74 Accessed 3 Oct 2017.

    Article  CAS  PubMed  Google Scholar 

  30. Hosmer DW, Lemeshow S, Sturdivant RX. Applied logistic regression, 3rd edition: John Wiley and Sons Ltd, New York, April 2013.,+3rd+Edition-p-9780470582473. Accessed 13 May 2018.

  31. Schneider MA, Christopoulos P, Muley T, Warth A, Klingmueller U, Thomas M, et al. AURKA, DLGAP5, TPX2, KIF11 and CKAP5: five specific mitosis-associated genes correlate with poor prognosis for non-small cell lung cancer patients. Int J Oncol. 2017;50:365–72.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Martens-de Kemp SR, Nagel R, Stigter-van Walsum M, van der Meulen IH, van Beusechem VW, Braakhuis BJM, et al. Functional genetic screens identify genes essential for tumor cell survival in head and neck and lung cancer. Clin Cancer Res. 2013;19:1994–2003.

    Article  CAS  PubMed  Google Scholar 

  33. Behrens C, Solis LM, Lin H, Yuan P, Tang X, Kadara H, et al. EZH2 protein expression associates with the early pathogenesis, tumor progression, and prognosis of non-small cell lung carcinoma. Clin Cancer Res. 2013;19:6556–65.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Huqun R, Ishikawa R, Zhang J, Miyazawa H, Goto Y, Shimizu Y, et al. Enhancer of zeste homolog 2 is a novel prognostic biomarker in nonsmall cell lung cancer. Cancer. 2012;118:1599–606.

    Article  CAS  PubMed  Google Scholar 

  35. Reijm EA, Timmermans AM, Look MP, Meijer-van Gelder ME, Stobbe CK, van Deurzen CHM, et al. High protein expression of EZH2 is related to unfavorable outcome to tamoxifen in metastatic breast cancer. Ann Oncol. 2014;25:2185–90.

    Article  CAS  PubMed  Google Scholar 

  36. van Kruijsdijk RCM, Visseren FLJ, Boni L, Groen HJM, Dingemans AMC, Aerts JGJV, et al. Pemetrexed plus carboplatin versus pemetrexed in pretreated patients with advanced non-squamous non-small-cell lung cancer: treating the right patients based on individualized treatment effect prediction. Ann Oncol. 2016;27:1280–6.

    Article  PubMed  Google Scholar 

  37. Gonen N, Assaraf YG. Antifolates in cancer therapy: structure, activity and mechanisms of drug resistance. Drug Resist Updat. 2012;15:183–210.

    Article  CAS  PubMed  Google Scholar 

  38. Shaw a T, Varghese a M, Solomon BJ, Costa DB, Novello S, Mino-Kenudson M, et al. Pemetrexed-based chemotherapy in patients with advanced, ALK-positive non-small cell lung cancer. Ann Oncol. 2013;24:59–66.

    Article  PubMed  Google Scholar 

  39. Park S, Ahn B-C, Lim SW, Sun J-M, Kim HR, Hong MH, et al. Characteristics and outcome of ROS1-positive non-small cell lung Cancer patients in routine clinical practice. J Thorac Oncol. 2018;13:1373–82.

    Article  PubMed  Google Scholar 

Download references


Not applicable


This study was supported by ZonMw, the Netherlands [grant number 152001017]. The funding was used to perform the immunohistochemical stainings on the TMAs of the patients in the validation cohort. The funding body had no role in the design of this study and collection, analysis and interpretation of the data or in writing the manuscript.

Availability of data and materials

The microarray data of the training cohort are available at the Gene Expression Omnibus of the National Center for Biotechnology Information (GSE19188).

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

Author information

Authors and Affiliations



JA, BS and SV formulated the study design. SV, JH, KB, LV, JPH and SP performed the data acquisition. Scoring of the immunohistochemical stainings was done by SV, KB and JPH. SV and BS completed the statistical analyses. All authors were involved in manuscript preparation and editing. All authors read and approved the final manuscript.

Corresponding author

Correspondence to J. G. J. V. Aerts.

Ethics declarations

Ethics approval and consent to participate

Written informed consent was obtained from all participants in the training cohort. Tissues were studied under an anonymous tissue protocol approved by the medical ethical committee of Erasmus University Medical Center. As data in the validation cohort were collected retrospectively, no approval by a medical research and ethics committee was necessary according to Dutch guidelines (

Consent for publication

Not applicable

Competing interests

Joachim Aerts has a consultant/advisory role with Eli Lilly and Company, Roche, Bristol-Myers Squibb, MSD and Boehringer Ingelheim. The other authors declare no conflicts of interest.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Table S1. Minimized signature for prediction of pemetrexed response. (DOCX 66 kb)

Additional file 2:

Table S2. Antibodies used for immunohistochemical analyses. (DOCX 33 kb)

Additional file 3:

Table S3. Interobserver agreement of the IHC staining score with regard to tumour quantity and intensity of staining. (DOCX 39 kb)

Additional file 4:

Figure S1. Scatter plots of gene expression levels of selected molecular markers and TS gene expression. Dot plots showing correlations between relative mRNA expression of TS and mRNA expression of TOP2A, MCM2, EZH2, CPA3, TPX2. Abbreviations: IHC, immunohistochemical; EZH2, Enhancer of zeste homolog; TOP2A, Topoisomerase II; TPX2, Microtubule Nucleation Factor; CPA3, Carboxypeptidase A3; MCM2, Minichromosome Maintenance Complex Component 2. (EPS 180 kb)

Additional file 5:

Figure S2. Scatter plots of gene expression levels of selected molecular markers with their protein expression level and the associated correlations. Dot plots showing correlations between relative mRNA expression and IHC staining score of TOP2A, MCM2, EZH2, CPA3, TPX2. Abbreviations: IHC, immunohistochemical; EZH2, Enhancer of zeste homolog; TOP2A, Topoisomerase II; TPX2, Microtubule Nucleation Factor; CPA3, Carboxypeptidase A3; MCM2, Minichromosome Maintenance Complex Component 2. (EPS 176 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Visser, S., Hou, J., Bezemer, K. et al. Prediction of response to pemetrexed in non-small-cell lung cancer with immunohistochemical phenotyping based on gene expression profiles. BMC Cancer 19, 440 (2019).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: