Nuclear Receptor Expression Defines a Set of Prognostic Biomarkers for Lung Cancer


Background:
The identification of prognostic tumor biomarkers that also would have potential as therapeutic targets, particularly in patients with early stage disease, has been a long sought-after goal in the management and treatment of lung cancer. The nuclear receptor (NR) superfamily, which is composed of 48 transcription factors that govern complex physiologic and pathophysiologic processes, could represent a unique subset of these biomarkers. In fact, many members of this family are the targets of already identified selective receptor modulators, providing a direct link between individual tumor NR quantitation and selection of therapy. The goal of this study, which begins this overall strategy, was to investigate the association between mRNA expression of the NR superfamily and the clinical outcome for patients with lung cancer, and to test whether a tumor NR gene signature provided useful information (over available clinical data) for patients with lung cancer.

Methods and Findings:
Using quantitative real-time PCR to study NR expression in 30 microdissected non-small-cell lung cancers (NSCLCs) and their pair-matched normal lung epithelium, we found great variability in NR expression among patients' tumor and non-involved lung epithelium, found a strong association between NR expression and clinical outcome, and identified an NR gene signature from both normal and tumor tissues that predicted patient survival time and disease recurrence. The NR signature derived from the initial 30 NSCLC samples was validated in two independent microarray datasets derived from 442 and 117 resected lung adenocarcinomas. The NR gene signature was also validated in 130 squamous cell carcinomas. The prognostic signature in tumors could be distilled to expression of two NRs, short heterodimer partner and progesterone receptor, as single gene predictors of NSCLC patient survival time, including for patients with stage I disease. Of equal interest, the studies of microdissected histologically normal epithelium and matched tumors identified expression in normal (but not tumor) epithelium of NGFIB3 and mineralocorticoid receptor as single gene predictors of good prognosis.

Conclusions:
NR expression is strongly associated with clinical outcomes for patients with lung cancer, and this expression profile provides a unique prognostic signature for lung cancer patient survival time, particularly for those with early stage disease. This study highlights the potential use of NRs as a rational set of therapeutically tractable genes as theragnostic biomarkers, and specifically identifies short heterodimer partner and progesterone receptor in tumors, and NGFIB3 and MR in non-neoplastic lung epithelium, for future detailed translational study in lung cancer.

: Please see later in the article for the Editors' Summary


Published in the journal: . PLoS Med 7(12): e32767. doi:10.1371/journal.pmed.1000378
Category: Research Article
doi: 10.1371/journal.pmed.1000378

Summary

Background:
The identification of prognostic tumor biomarkers that also would have potential as therapeutic targets, particularly in patients with early stage disease, has been a long sought-after goal in the management and treatment of lung cancer. The nuclear receptor (NR) superfamily, which is composed of 48 transcription factors that govern complex physiologic and pathophysiologic processes, could represent a unique subset of these biomarkers. In fact, many members of this family are the targets of already identified selective receptor modulators, providing a direct link between individual tumor NR quantitation and selection of therapy. The goal of this study, which begins this overall strategy, was to investigate the association between mRNA expression of the NR superfamily and the clinical outcome for patients with lung cancer, and to test whether a tumor NR gene signature provided useful information (over available clinical data) for patients with lung cancer.

Methods and Findings:
Using quantitative real-time PCR to study NR expression in 30 microdissected non-small-cell lung cancers (NSCLCs) and their pair-matched normal lung epithelium, we found great variability in NR expression among patients' tumor and non-involved lung epithelium, found a strong association between NR expression and clinical outcome, and identified an NR gene signature from both normal and tumor tissues that predicted patient survival time and disease recurrence. The NR signature derived from the initial 30 NSCLC samples was validated in two independent microarray datasets derived from 442 and 117 resected lung adenocarcinomas. The NR gene signature was also validated in 130 squamous cell carcinomas. The prognostic signature in tumors could be distilled to expression of two NRs, short heterodimer partner and progesterone receptor, as single gene predictors of NSCLC patient survival time, including for patients with stage I disease. Of equal interest, the studies of microdissected histologically normal epithelium and matched tumors identified expression in normal (but not tumor) epithelium of NGFIB3 and mineralocorticoid receptor as single gene predictors of good prognosis.

Conclusions:
NR expression is strongly associated with clinical outcomes for patients with lung cancer, and this expression profile provides a unique prognostic signature for lung cancer patient survival time, particularly for those with early stage disease. This study highlights the potential use of NRs as a rational set of therapeutically tractable genes as theragnostic biomarkers, and specifically identifies short heterodimer partner and progesterone receptor in tumors, and NGFIB3 and MR in non-neoplastic lung epithelium, for future detailed translational study in lung cancer.

: Please see later in the article for the Editors' Summary

Introduction

The prevalence of lung cancer as the primary cause of cancer death in the United States has led to renewed efforts to obtain biomarker signatures that provide either prognostic or predictive information to guide therapy for individual patients (i.e., “personalized medicine”) [1][3]. Multiple genome-wide expression studies have demonstrated the usefulness of this approach for lung cancer prognosis [4][6]. However, in general, each of these studies has identified different sets of genes, even when the various studies are used to cross-validate one another, and in the majority of these studies the individual genes identified have had little impact for understanding tumor pathogenesis or as therapeutic targets. Thus, identification of prognostic biomarkers that also provide hypotheses for mechanism-based studies of carcinogenesis and offer new therapeutic targets (sometimes referred to as “theragnostics”) would be of significant benefit.

Nuclear receptors (NRs) are a large family of ligand-dependent transcription factors that respond to a number of hormonal and diet-derived lipids, including endocrine steroids, fat-soluble vitamins, fatty acids, and cholesterol metabolites [7]. NRs are also among the most successful targets of drugs approved to treat many diseases, including cancer [8]. Previously, we have shown that NR expression profiling can be used to reveal the mechanistic basis of the hierarchical transcriptional networks that govern a number of physiological processes, including development, differentiation, reproduction, circadian rhythm, and metabolism [9][13]. In the present study, we wished to investigate the potential role of the 48 members of the NR superfamily as “theragnostic” indicators in lung cancer. The strategy of examining expression of NRs, which are known therapeutic targets with defined mechanisms of action, differs from previous, open-ended genome-wide microarray studies. Thus, our goal was to determine whether one can use NR expression signatures as clinical tools for prognosis for patients with lung cancer, which also might lead to NR-selective therapies targeted at hormonal manipulation of lung cancer. As an additional aspect of this study, we wished to provide open access to our data by including a Sweave document [14][16] that contains a literate programming package to permit the full reproduction of our analysis.

Methods

Sample Collection

Samples from MD Anderson Cancer Center

All tissue samples were obtained by surgical resection from patients who had provided written informed consent under approval of the institutional review boards at MD Anderson Cancer Center (MDACC). Tissues were stored at −80°C after being snap frozen in liquid nitrogen. Serial sectioning of each sample was used to histologically evaluate tumor and normal tissue for subsequent microdissection [17]. Thirty primary tumor and corresponding normal tissues (including 22 adenocarcinomas and eight squamous cell carcinomas) were selected randomly from 379 similar samples in the MDACC lung tumor collection based on stringent, predefined quality control procedures before any data analysis (see Figure 1 for overall study design). Detailed sample selection procedures are described in Text S1. Among the 30 patient samples, 17 were diagnosed with stage I (the earliest stage) disease, four with stage II disease, five with stage III disease, and four with stage IV disease. A comparative analysis of the 30 samples demonstrated that they represent an objective sampling of the MDACC lung cancer tumor collection (Figure S1; Table S1). RNAs were isolated from each sample using the Qiagen RNeasy Mini Kit (Qiagen Sciences).

Schematic of the study design for development and validation of the NR prognostic gene signature.
Fig. 1. Schematic of the study design for development and validation of the NR prognostic gene signature.
The flow chart details the design and implementation of this study. ADC, adenocarcinoma; SCC, squamous cell carcinoma.

Samples from the National Cancer Institute Director's Consortium and other datasets

Validation of the prediction model from the MDACC samples was performed using an independent dataset from a recently published National Cancer Institute (NCI) Director's Consortium study of lung cancer involving 442 resected non-small-cell lung cancers (NSCLCs) [6]. This Consortium dataset represents one of the largest microarray datasets of NSCLC samples that has been collected and studied using a common protocol. From this study, the Affymetrix U133A microarray data for the 48 NR gene expression signatures were excerpted and used as shown in Figure 1 to validate the prognostic value of NR expression. In a similar manner, the MDACC signature was validated in the NR expression microarray data taken from a second set of 117 adenocarcinomas [18]. Finally, the NR signature from the Consortium adenocarcinoma dataset was cross-validated using the 48-gene NR expression profile excerpted from the microarray of 130 squamous cell carcinomas [19].

Reverse Transcription and Quantitative PCR Assay

All cDNAs were prepared for quantitative PCR (QPCR) (TaqMan method) as previously described [10]. Briefly, 2 µg of total RNA was DNAse-treated with 2 U of DNAse I in final volume of 20 µl containing 4.2 µM MgCl2. The reverse transcription reaction was performed in 100 µl final volume, followed by addition of 100 µl of DEPC-H2O. Human universal cDNAs for broadly expressed NRs or tissue-specific cDNAs for restricted-expression NRs was used to construct a standard curve of the following concentrations: 0 (i.e., no template control), 0.008, 0.04, 0.2, 1, 5, and 25 ng of 18S RNA; and 0 (i.e., no template control), 0.016, 0.08, 0.4, 2.0, 10, and 50 ng of each NR RNA. These quantities were based on the RNA concentration used for the reverse transcription reaction. A negative reverse transcription sample and a control for genomic DNA contamination were included for both 18S and NR. Per sample, 10 ng of cDNA was assayed in triplicate wells of a 384-well plate. The final forward and reverse primer concentrations used were 75 nM for 18S rRNA and 300 nM for all NRs. For this study the 48 NRs plus the two common splice variants for PPARγ (i.e., PPARγ2) and PPARδ (i.e., PPARδ2) were included in the analysis of all samples from the MDACC patient set. Primer sequences have been reported elsewhere (http://www.nursa.org/datasets.cfm?doi=10.1621/datasets.05005).

QPCR Data Analysis

Data were imported into Microsoft Excel and evaluated for PCR efficiency (e), e = 10−1/slope, where the slope was obtained from the standard curve calculated by the sequence detection system software of the ABI7900 instrument (Applied Biosystems) for the endogenous 18S reference and target NR. Relative mRNA amounts were calculated by quantity  =  e–Ct, where Ct is cycle time. The calculated quantities were averaged (avg), and the standard deviations (stdev) and coefficients of variation (CV  =  stdev/avg) were determined for the 18S and NR of each sample. Data points that showed CV >17% were considered outliers and removed. A 17% CV cutoff correlates with the maximum allowable standard deviation that can distinguish a 2-fold change with 99% confidence when samples are assayed in triplicate wells for both the endogenous reference and the gene of interest. Note that only one point per replicate may be removed. Normalized values for expression of each NR were calculated using normalized value  =  NR quantity avg/18S quantity avg. The standard deviation of the normalized value was calculated as (normalized value) × [(CV of reference)2 + (CV of gene of interest)2]½. Normalized values are represented as a bar graph. QPCR data analysis procedures have been described previously [10],[12]. The entire QPCR dataset of NR expression in normal and tumor samples from the 30-patient cohort is available in Figure S2 and online at http://www.nursa.org/datasets.cfm?doi=10.1621/datasets.05010.

Microarray Data Preprocessing

Consortium microarray raw data [3],[6] were downloaded from the NCI's caArray database and preprocessed by robust multichip average background correction and quantile–quantile normalization [20]. All gene expression values were log-transformed (on a base 2 scale). Average values were used for the different probe sets corresponding to the same gene.

Unsupervised Clustering Analysis

The hierarchical clustering algorithm [21] was used to group NR expression versus the 30-MDACC-patient cohort based on the QPCR expression profile. Gene expression values were log-transformed (on a base 2 scale) in a manner similar to the transformation of the microarray data. Euclidian distance and average link were used in the hierarchical clustering algorithm.

Supervised Classification Analysis

Supervised classification was performed using Recursive Partitioning and Regression Trees (RPART) [22][25], and was implemented using R software version 2.10.0. A detailed description of the implementation is provided in Text S1.

Survival Analysis

Overall survival time was calculated from the date of surgery until death or the last follow-up contact. Recurrence-free survival time was defined as the time interval between the date of surgery and the date of disease recurrence or death from any cause, whichever came first, or date of last follow-up evaluation. Survival curves were estimated using the product-limit method of Kaplan-Meier [26] and were compared using the log-rank test. Cox proportional-hazards analysis [27] was also performed, with survival time as the dependent variable.

Sweave Report

A Sweave document is provided () to permit others to reproduce any or all parts of our statistical analyses. Sweave is a literate programming R package that combines the source code (in R) and documentation (in LaTeX) in one file and thereby permits reproduction of published high-throughput data analysis [14][16].

Results

Association of NR Gene Expression with Survival Time

The overall schema for the design and data analysis of this study is shown in Figure 1. The expression of all 48 members (plus two splice isoforms) of the NR superfamily was investigated first in a 30-patient sample set to develop a prognostic signature. Detailed clinical data on the 30-patient cohort are given in Tables S1S3. The QPCR datasets of NR expression are shown and summarized in Figure S2 and Table S4. Raw datasets are available at http://www.nursa.org/datasets.cfm?doi  = 10.1621/datasets.05010.

Univariate analysis

We first explored whether the mRNA expression of NR genes was associated with survival outcomes for patients with lung cancer. Table S5 shows hazard ratios (HRs) and corresponding p-values for the association between each NR gene and survival outcome from the univariate Cox regression models. The mRNA expression of 37 NR genes was significantly associated with survival outcomes (p<0.05). Figure S3 shows the cumulative distribution function (CDF) of the p-values representing the associations between individual NR gene expression and survival time from univariate Cox models. These data show that the NR gene superfamily as a whole has statistically significant association with survival outcome (p<0.00001) based on a Kolmogorov-Smirnov test of the difference between the p-value distributions of NR sets and random gene sets.

Unsupervised clustering analysis

Unsupervised cluster analysis (with euclidian distance and average link) of NR expression in lung tumors revealed two distinct clusters of tissue samples (Figure 2A). Note that one tissue sample (773-ADC) did not fall into either cluster and was treated as an outlier for the clustering analysis. The two major branches of the dendrogram (cluster 1 and cluster 2) were associated with both overall survival rates (p = 0.0079) and disease recurrence rates (p = 0.0927), but no other clinical features (Table 1). Indeed, Kaplan-Meier plots for survival time and disease recurrence showed that cluster 1 and cluster 2 segregated patients into those with poor and good prognostic outcomes, respectively (p = 0.00008 for survival time; p = 0.0051 for disease recurrence) (Figure 2B and 2C). These findings suggest that the gene expression of the NR family is strongly associated with clinical outcomes for patients with lung cancer.

QPCR analysis of the NR gene expression signature in patients with lung cancer.
Fig. 2. QPCR analysis of the NR gene expression signature in patients with lung cancer.
(A) Unsupervised cluster analysis of the 30-patient MDACC lung cancer cohort using the QPCR profile of the NR superfamily. Horizontal and vertical axes represent NR and lung cancer patient clusters, respectively. (B and C) Kaplan-Meier plots showing the association of the NR gene signature with overall patient survival (B) and disease recurrence (C). p-values were obtained using the log-rank test. Red represents sample Cluster I and blue represents Cluster II, defined by an unsupervised clustering algorithm using the NR gene profiling data in (A). Circles indicate censored samples. ADC, adenocarcinoma; SCC, squamous cell carcinoma.

Tab. 1. Patient demographics summarized by unsupervised cluster analysis of lung tumors.
Patient demographics summarized by unsupervised cluster analysis of lung tumors.
Indicates p-values by t test for age and by Fisher's exact test for the other variables comparing cluster 1 and 2.

Development and Validation of a Prognostic NR Signature for Lung Cancer

Standard classification and regression tree analysis [22][25] was used to build a prognostic model based on NR gene expression in the MDACC NSCLC 30-patient cohort. We first evaluated the performance of the prognosis model using the leave-one-out cross-validation (LOOCV) method. The HR, i.e., risk of death, for the predicted high-risk versus the predicted low-risk signatures using tumor samples was 13.6; 95% confidence interval (CI), 3.07–60.92; p = 0.000014 (Figure 3A). The relatively large CI for the HR in the MDACC cohort is due to the small sample size in this dataset. These results prompted us to validate the NR signature in a larger dataset.

Kaplan-Meier plots showing the predictive power of the NR gene signature in datasets from the NCI Director's Consortium.
Fig. 3. Kaplan-Meier plots showing the predictive power of the NR gene signature in datasets from the NCI Director's Consortium.
(A) LOOCV of the recursive partitioning tree model (RPART) for the 30-sample MDACC QPCR dataset using all 48 NRs. The HR for the predicted high-risk versus the predicted low-risk signatures was 13.6; 95% CI, 3.07–60.92; p = 0.000014. (B and C) Independent validation of the 48-NR gene expression signature between the MDACC cohort and the Consortium cohort. The MDACC cohort training set (n = 30) was tested in the Consortium cohort (n = 442) (B), and vice versa (C). (D) Independent validation of the NR gene signature in the 442-sample multi-institute Consortium using RPART analysis. The microarray datasets were divided into two groups, one for the training cohort (n = 256) and the other for the testing cohort (n = 186). p-values were obtained by the log-rank test. Red and black lines represent predicted high- and low-risk groups, respectively. Circles indicate censored samples. Note that SHP expression is the single co-variable (or predictor) in the classification model that describes the survival time differences shown in (B), demonstrating that SHP is a single gene predictor that represents the entire 48-NR gene profile.

Because the majority of gene expression data now available from lung cancer samples comes from microarray expression studies, we investigated whether the NR expression profile could be validated within a completely independent dataset taken from the NCI Director's Consortium study of lung cancer involving 442 resected NSCLCs [6]. We first validated the 30-sample QPCR dataset on the 442-sample microarray dataset, and then we developed an NR signature from the microarray dataset and validated it on the QPCR data. Both directions of training and testing provided statistically significant predictive power for patient survival time (Figure 3B and 3C).

As a further validation test of the NR gene signature, we divided the 442-sample microarray data into training and testing sets, and analyzed the data using the predictive RPART model. We used the same training and testing strategy as in the genome-wide analyses of these data [6]. The training set (n = 256) included samples from the University of Michigan Cancer Center (n = 177) and the Moffitt Cancer Center (n = 79), and the testing set (n = 186) included samples from the Memorial Sloan-Kettering Cancer Center (n = 104) and the Dana-Farber Cancer Institute (n = 82). Using the NR expression profile from training data to build a predictive model yielded a HR of 2.04 (95% CI, 1.12–3.72; p = 0.018) for the predicted high-risk versus the predicted low-risk signature in testing data (Figure 3D). The higher HR value for the QPCR dataset likely reflects the greater dynamic range and quantitative nature of the QPCR assay, and the greater homogeneity of the microdissected samples.

As yet further confirmation of the prognostic NR gene signature, the NR prognosis signature developed using the MDACC dataset was validated in another dataset containing 117 adenocarcinoma samples [18]. Again, the results show that patients in the predicted low-risk group have significantly longer survival times than those in the predicted high-risk group (p = 0.0053; Figure S4). The results from this second independent dataset further confirmed the robustness of the NR prognosis signature.

Finally, we also examined whether the NR prognosis signature might be applicable to squamous cell carcinomas. To that end, the NR prognosis signature developed using the Consortium lung adenocarcinoma dataset was validated in a previously published dataset from Raponi et al. containing 130 lung squamous cell carcinomas [19]. Patients in the predicted low-risk group had significantly longer survival times than those in the predicted high-risk group (p = 0.018; Figure S5A). We then used lung squamous cell carcinomas (from the Raponi et al. dataset) to develop a NR prognosis signature and validated it in the lung adenocarcinomas (Consortium dataset). Again, the results showed that NR expression has significant predictive power (Figure S5B). Based on these results, we conclude that the NR prognosis signature may not be limited only to adenocarcinomas.

To underscore the significance of the 48-NR gene signature, we used 1,000 random gene sets, each comprising 48 transcripts, to build prediction models in the Consortium training set and then tested them in the Consortium testing set. The distribution of the association between the predicted risk groups and the survival outcomes (p-values) is shown in Figure S6; only 28 of the 1,000 random sets reached the significance level of the NR signature (p = 0.018), thereby demonstrating that the NR signature is highly selective and nonrandom. Taken together, these results strongly support the utility of the NR gene signature as a prognostic marker, even when applied and cross-validated independently by two different gene expression platforms (i.e., QPCR and microarray) and in multiple independent lung cancer datasets.

NR Signature Prediction Is Independent of Clinical Variables

We examined whether the association between the NR signature and the survival outcome was independent of clinical variables by using a multivariate Cox proportional-hazards analysis that included NR risk score, gender, age at diagnosis, use of adjuvant chemotherapy, use of adjuvant radiation therapy, and stage as the co-variables. We analyzed the Consortium testing dataset, which included samples from the Memorial Sloan-Kettering Cancer Center and the Dana-Farber Cancer Institute. The NR risk scores used in this analysis were derived from the prediction model built from the Consortium training dataset (from the University of Michigan Cancer Center and the Moffitt Cancer Center). Surprisingly, this multivariate analysis revealed that the association between NR risk scores and survival time was independent of other clinical variables (HR = 1.98; 95% CI, 1.04–3.77; p = 0.037) (Table 2). These data reveal that a significant association exists between a patient's NR profile and survival time when adjusted for other clinical variables. As expected, the correlation between tumor stage and patient survival time was also highly significant, confirming this clinical feature as a well-recognized prognostic marker used in the clinic. It is interesting to note that gender also was significantly correlated with patient survival time in our analysis (males had a higher risk of death than females).

Tab. 2. Death HRs from multivariate Cox regression analysis from two independent datasets.
Death HRs from multivariate Cox regression analysis from two independent datasets.
Memorial Sloan-Kettering Cancer Center (MSK) and Dana-Farber Cancer Institute (CAN/DF).

Next, we investigated whether combining clinical variables and NR gene expression could improve upon the prognosis based on clinical variables alone. The Consortium microarray dataset was divided into two groups, one for the training cohort and the other for the testing cohort. From these samples, the two principal components of NR gene expression and clinical variables (including gender, age, stage, and adjuvant therapy) were used to build a prediction model using the classification tree approach (Figure S7). This analysis shows that the predicted risk groups from the prognosis model based on the combination of clinical variables and NR expression have a stronger association with survival time (HR = 3.21, 95% CI, 2.01–5.12; p = 2.43×10−7) than models using the NR signature alone (HR = 1.98, 95% CI, 1.71–4.30; p = 0.037) or the clinical variables alone (HR = 2.71, p = 1.02×10−5).

Refinement of the NR Signature into Single Gene Predictors

We next explored the roles of specific NRs in the prediction models. To address this question, we further re-interrogated the classification tree model to see whether the prediction model outcomes shown in Figures 3A and 3B were due to the dominant effects of any specific NRs (see Text S1 for details). When using all 48 NR genes as input for the classification tree model, short heterodimer partner (SHP) expression was identified as the only co-variable left in the final RPART prediction model built from the 30-patient MDACC dataset. In other words, the prognosis performance of the 48-NR gene signature (shown in Figure 3A and 3B) is the same as that using SHP expression alone to build the models. Next, we removed SHP from the dataset, reanalyzed the prediction model, and found the model still had remarkable prognostic ability in the 30-patient LOOCV dataset (HR = 9.58; 95% CI, 3.00–30.6; p = 0.0000065) (Figure 4A). Furthermore, the classification tree structure revealed that when the prediction model excluded SHP, progesterone receptor (PR) was the single gene signature used. The detailed classification tree structures are shown in the Sweave document (Text S2) and plotted in Figure S8. The single gene PR signature was further validated in the testing cohort of the Consortium dataset (HR = 1.46; 95% CI, 1.12–1.90; p = 0.0048) (Figure 4B). The protective effect based on SHP and PR expression was further strengthened by univariate Cox regression modeling, which consistently showed that expression of both NRs correlated with significantly lower HRs in the microarray dataset (Figure 4C and 4D). Thus, SHP and PR represent single gene markers.

Identification of single NR gene biomarkers for lung cancer prognosis.
Fig. 4. Identification of single NR gene biomarkers for lung cancer prognosis.
(A and B) Kaplan-Meier survival plots using PR in the single gene prediction model. The MDACC cohort was tested using LOOCV (A), or it was used as a training set and independently tested in the multi-site Consortium cohort (B). For this analysis mRNA expression values for SHP were removed from the dataset in order to test the effect of other NR genes as biomarkers. In this case, PR expression is the single co-variable (or predictor) in the classification model that describes the survival time differences shown in (B), demonstrating that PR is a single gene predictor that represents the NR gene profile when SHP expression is excluded. p-values were obtained using the log-rank test. Red and black lines represent high- and low-risk groups, respectively. Circles indicate censored samples. (C and D) HRs from univariate Cox regression models for SHP and PR expression, respectively, in the MDACC and multi-site Consortium datasets.

SHP Expression Predicts Survival Time in Patients with Early Stage Lung Cancer

Since the prognosis for patients with early stage (i.e., stage I) lung cancer has substantial clinical impact on guiding therapeutic strategy, we tested whether expression of specific NRs also has predictive power to refine the prognosis of patients with stage I lung cancer. We found that the predicted risk groups using SHP as a single gene classifier were significantly associated with survival outcome for patients with stage I lung cancer in the Consortium samples (Figure 5A; p = 0.033), whereas the PR signature was marginally predictive (Figure 5B; p = 0.069).

Kaplan-Meier survival plots showing single NR gene predictors in patients with stage I lung cancer.
Fig. 5. Kaplan-Meier survival plots showing single NR gene predictors in patients with stage I lung cancer.
Predictive models for SHP (A) and PR (B) were trained in the MDACC samples and tested in the patients with stage I lung cancer in the Consortium cohort. p-values were obtained by the log-rank test. Red and black lines represent predicted high- and the low-risk groups, respectively. Circles indicate censored samples.

NR Expression in Normal Tissue Predicts Survival and Disease Recurrence

We also examined the potential prognostic value of NR expression in histologically normal lung tissue obtained from areas adjacent to the tumors of the MDACC cohort used in the above studies. When the normal tissue expression data were analyzed using the classification tree model and validated by LOOCV, the NR signature provided statistically significant predictors of both disease recurrence (HR = 4.61, 95% CI, 1.74–12.30; p = 0.00099) and overall patient survival time (HR = 2.22, 95% CI, 0.85–5.81; p = 0.094) (Figure S9). In contrast to the identification of SHP and PR in the tumor samples as key predictors, the normal lung epithelium profile classification tree structures revealed two other NRs, NGFIB3 (nerve growth factor induced gene B3, also known as NR4A1) and mineralocorticoid receptor (MR, also known as NR3C2) to be single gene predictors for survival and disease recurrence. Although the prediction models for normal tissue require further validation in an independent dataset, this analysis suggests higher expression of NGFIB3 and MR is associated with a good prognosis.

To further understand the relationship between NR expression levels in the tumor versus adjacent normal tissue, we performed a pair-wise correlation of NR gene expression in tumor and adjacent normal lung from each patient's paired tissue set. For this analysis we used both the MDACC dataset and a recently published microarray dataset from Landi et al. comparing gene expression in 135 NSCLC adenocarcinomas and adjacent normal tissue [28]. In both datasets we found that the NR expression in tumor samples correlated significantly with NR expression in adjacent normal samples (Pearson correlations were 0.87, p<0.001, and 0.92, p<0.001, respectively) (Figure S10). This analysis suggests that, in general, there is a tight correlation between tumor and normal tissue NR expression within the same patient. We next asked which of the 48 NRs might be differentially expressed between tumor and adjacent normal samples using paired t tests. Again, the t tests derived from the two datasets were consistent (Pearson correlation  = 0.7, p<0.001). Furthermore, five of the top ten NR genes that were shown to differ significantly between tumor and normal tissue were common to both datasets (i.e., AR, MR, NGFIB3, PPARγ, and RXRγ). Of note was the finding that MR and NGFIB3 were among these receptors, further emphasizing their prognostic potential.

Discussion

In this study, we investigated the expression of NRs in lung cancer and found that specifically targeting this superfamily of ligand-dependent transcription factors provided a novel prognostic biomarker. Several recent studies using genome-wide microarray experiments have proposed various sets of genetic signatures for lung cancer prognosis [4][6],[29][31]. Interestingly, the gene signatures from these studies have shared little if any overlap with one another [4],[5],[29]. In addition, because of the open-ended nature of genome-wide analyses, the signatures have provided little insight into the pathogenesis or pathophysiology of lung cancer. Moreover, these studies have yet to identify new therapeutic targets. In contrast, NRs represent a well-studied class of proteins that (1) govern complex cellular programs such as differentiation, inflammation, and metabolism [9][13], (2) are known transcriptional drivers of oncogenesis, and (3) are themselves the targets of validated drugs for many diseases including cancer [32][35]. This superfamily also includes a number of orphan receptors, many of which are currently being evaluated as potential new therapeutic targets for a number of diseases [8]. Thus, our study suggests that specifically targeting NRs may provide an alternative and clinically relevant new strategy for profiling lung tumors.

The immediate findings from this work also may have a number of important and practical implications for the use of the NR gene signature in a clinical setting. First, we demonstrated that the NR superfamily gene expression signature is an excellent predictor of both patient survival and recurrence of lung cancer that is comparable to predictors from other studies. As an example, Shedden et al. [6] compared six previously published gene signatures and showed that the HRs from these studies ranged from 0.93 to 2.30 (p-values ranged from p<0.001 to p = 0.78). In comparison, the HR for our NR signature obtained in the Consortium dataset [6] was 2.04 (p = 0.018) (Figure 3D).

The prognostic potential of NRs was validated in independent datasets, and further analysis identified tumor expression of SHP (discussed below) and PR as robust single gene predictors. The demonstration of PR as a predictive marker is in line with a previous retrospective study where PR was shown to be associated with survival time in patients with lung adenocarcinoma [36]. Expression of PR, together with estrogen receptor α (ERα), is now well established as a clinical guide to both prognostic anticipation and therapeutic intervention for breast cancer. Indeed, in thinking about the next step in our studies, the finding that certain lung cancers express specific, known therapeutic NR targets (e.g., PR, ERα, ERβ, AR, RARs, and PPARs) brings up the possibility of treating patients whose tumors express these receptors with drugs (agonists or antagonists) that target the receptors. Although it remains to be established whether NR mRNA levels will correlate with NR protein expression, several recent studies support the concept that selective receptor modulators may be effective therapeutically. For example, a number of studies have suggested anti-estrogen therapy as a lung cancer therapeutic [37][40], and in a mouse lung cancer study, the use of a PPARγ agonist had a synergistic effect in reducing tumor burden when used with cis-platinum [41]. Interestingly, treatment with progesterone has been shown to inhibit lung tumor xenograft growth in a preclinical study [36], whereas, in contrast, the Women's Health Initiative study reported increased death from lung cancer in postmenopausal women treated with both estrogen and progestins [42]. At present it is not clear whether this latter finding was due to the presence of estrogen or progestin, highlighting the need for further investigation. Nevertheless, a reasonable assumption based on the present study is that predicting responses to drugs like anti-estrogens might be accomplished by screening patients for NR expression using the methodology highlighted in this study. Evaluation of the QPCR profiles from our study revealed a high degree of patient-to-patient variability in NR expression (Figure S2), and this observation provides a strong rationale for using this approach to guide individualized treatment in the future. Similarly, our data suggest that NR profiling of individual tumors provides a clinical paradigm for identifying potential responders to NR drugs.

A second notable finding was that the NR signature could be used to refine the prognosis for patients with early stage lung cancer. Although the clinical characteristics from the MDACC dataset and the Consortium dataset are not homogeneous and the sample size is not big enough to validate the prognosis signatures for stages IA, IB, and II separately, the positive prognosis results across the two datasets indicate that the signatures may be applied across heterogeneous populations of patients. In interrogating this signature further, it was of considerable interest to find that the orphan nuclear receptor SHP is a singe gene prognostic lung cancer biomarker of early stage disease. SHP has been extensively studied for its role in liver lipid metabolism [43],[44] and as transcriptional repressor of other NRs [45]. Intriguingly, a recent report found that SHP expression was negatively associated with liver tumorigenesis in a mouse model [46]. These findings prompt further exploration into whether there is a connection between the known physiological role of SHP and lung tumorigenesis, or whether SHP has a unique pathophysiologic function in the disease pathogenesis. To that end, we note that FXR agonists, a PPARγ agonist (rosiglitazone), agents that inhibit HNF-1α action, and a number of orphan drugs are all inducers of SHP expression [47]. These compounds might be tested in preclinical models to see whether they can induce SHP expression and inhibit lung tumorigenesis or malignant behavior. Also germline mutations in SHP or polymorphisms in FXR that regulate the level of SHP expression could play a role in SHP function in lung cancer pathogenesis or behavior. Together, these findings suggest that SHP expression may not only be a prognostic biomarker, but its presence in tumors may influence the expression of other genes. Indeed, in a preliminary experiment to address this idea we analyzed genome-wide RNA expression data using the Consortium samples to test the expression association between SHP (the single gene predictor) and a number of its known target genes, including Cyclin D1, Glut4, MTP, and PGC-1α. We found that SHP expression was associated with expression of three out of the four targets: Cyclin D1, Glut4, and MTP (p<0.0001).

A third noteworthy finding from our study was the ability to predict overall survival based on NR expression in normal tissue of patients with lung cancer. Whether this may be due to a “field effect” through changes in the epithelium, to induction in a paracrine fashion by the cancer, or to some other feature of normal lung epithelium is unknown. However, this finding does suggest that interrogating the histologically normal tissue may yield insight into lung cancer oncogenesis. To that end we note that the prognostic NR signature in normal tissue is completely different than that of the adjacent tumor. In contrast to that observed in tumors, the NR signature, as distilled using RPART analysis, revealed that NGFIB3 (NR4A1) and MR are single gene biomarkers found in normal tissue for predicting disease recurrence and overall survival time, respectively. NR4A family members have been shown to be tumor suppressors in a mouse model of myeloid leukemogenesis [48]. Similarly, low expression of MR has been shown to correlate with colorectal carcinoma recurrence [49]. These studies support the notion that higher expression of NR4A1 and MR might play a protective role against lung tumor pathogenesis.

A fourth finding of our study was the independent demonstration that the NR gene signature could be tested and cross-validated using two different gene expression platforms, QPCR and microarray. Given that microarray data do not have the dynamic or quantitative properties of data generated by QPCR, the cross-validation of the NR gene signature between different platforms strongly supports the idea that the NR superfamily may be a powerful prognostic predictor that also is functionally involved in lung cancer pathophysiology. Our results also suggest that a combination of a more robust collection process (microdissection instead of tissue mass) together with more quantitative measurements (QPCR instead of microarray) may reduce variability and strengthen the data. Indeed, the 95% CIs of HRs for both SHP and PR genes from the 30-patient dataset were smaller than those from Consortium data (with a sample size of 442) (Figure 4C and 4D). HRs of the high-risk versus low-risk group, defined using unsupervised cluster results, were also higher for the 30-patient dataset (Figure 2B) than for the Consortium data (Figure S11). Thus, while labor intensive, improving sample homogeneity and the quality of the expression data is likely to provide more reliable prognostic information.

Finally, the NR expression profile provides specific, testable hypotheses on the role of the NRs in lung cancer pathogenesis. For example, blocking the function of a highly expressed tumor cell NR could inhibit tumor growth or development, while over-expressing a low-abundance tumor cell NR could test its tumor suppressive capability. The finding that non-neoplastic tissue within the vicinity of the tumor also provided an NR gene signature that was predictive for survival time may provide the basis for testing NR function in the normal lung epithelium airway field where lung tumors develop. Perhaps one of the most surprising observations from this study is that an NR signature has not appeared in the prognostic signatures obtained in any of the previous global gene expression studies. This is true in spite of the fact that, at least in the multi-site Consortium database we analyzed, excerpting just the NR expression information yielded a predictive NR gene signature that was not discovered using global gene analysis [3],[6]. Thus, our study provides a strategic rationale for using an informed candidate gene profiling approach to identify prognostic markers and to interrogate specific gene families that may play roles in the cancer biology.

Supporting Information

Attachment 1

Attachment 2

Attachment 3

Attachment 4

Attachment 5

Attachment 6

Attachment 7

Attachment 8

Attachment 9

Attachment 10

Attachment 11

Attachment 12

Attachment 13

Attachment 14

Attachment 15

Attachment 16

Attachment 17

Attachment 18


Zdroje

1. JemalA

SiegelR

WardE

HaoY

XuJ

2009 Cancer statistics, 2009. CA Cancer J Clin 59 225 249

2. SunS

SchillerJH

SpinolaM

MinnaJD

2007 New molecularly targeted therapies for lung cancer. J Clin Invest 117 2740 2750

3. XieY

MinnaJD

2008 Predicting the future for people with lung cancer. Nat Med 14 812 813

4. ChenHY

YuSL

ChenCH

ChangGC

ChenCY

2007 A five-gene signature and clinical outcome in non-small-cell lung cancer. N Engl J Med 356 11 20

5. PottiA

MukherjeeS

PetersenR

DressmanHK

BildA

2006 A genomic strategy to refine prognosis in early-stage non-small-cell lung cancer. N Engl J Med 355 570 580

6. SheddenK

TaylorJM

EnkemannSA

TsaoMS

YeatmanTJ

2008 Gene expression-based survival prediction in lung adenocarcinoma: a multi-site, blinded validation study. Nat Med 14 822 827

7. ChawlaA

RepaJJ

EvansRM

MangelsdorfDJ

2001 Nuclear receptors and lipid physiology: opening the X-files. Science 294 1866 1870

8. ShulmanAI

MangelsdorfDJ

2005 Retinoid x receptor heterodimers in the metabolic syndrome. N Engl J Med 353 604 615

9. BarishGD

DownesM

AlaynickWA

YuRT

OcampoCB

2005 A Nuclear Receptor Atlas: macrophage activation. Mol Endocrinol 19 2466 2477

10. BookoutAL

JeongY

DownesM

YuRT

EvansRM

2006 Anatomical profiling of nuclear receptor expression reveals a hierarchical transcriptional network. Cell 126 789 799

11. FuM

SunT

BookoutAL

DownesM

YuRT

2005 A Nuclear Receptor Atlas: 3T3-L1 adipogenesis. Mol Endocrinol 19 2437 2450

12. XieCQ

JeongY

FuM

BookoutAL

Garcia-BarrioMT

2009 Expression profiling of nuclear receptors in human and mouse embryonic stem cells. Mol Endocrinol 23 724 733

13. YangX

DownesM

YuRT

BookoutAL

HeW

2006 Nuclear receptor expression links the circadian clock to metabolism. Cell 126 801 810

14. CoombesKR

WangJ

BaggerlyKA

2007 Microarrays: retracing steps. Nat Med 13 1276 1277; author reply 1277–1278

15. GentlemanR

2005 Reproducible research: a bioinformatics case study. Stat Appl Genet Mol Biol 4 Article2

16. LamportL

1994 LaTeX: A document preparation system. Boston Addison Wesley

17. MaitraA

GazdarAF

WistubaII 2001 Microdissection and the study of cancer pathways. Curr Mol Med 1 153 162

18. TomidaS

TakeuchiT

ShimadaY

ArimaC

MatsuoK

2009 Relapse-related molecular signature in lung adenocarcinomas identifies patients with dismal prognosis. J Clin Oncol 27 2793 2799

19. RaponiM

ZhangY

YuJ

ChenG

LeeG

2006 Gene expression signatures for predicting prognosis of squamous cell and adenocarcinomas of the lung. Cancer Res 66 7466 7472

20. BolstadBM

IrizarryRA

AstrandM

SpeedTP

2003 A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics 19 185 193

21. EisenMB

SpellmanPT

BrownPO

BotsteinD

1998 Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci U S A 95 14863 14868

22. GarzottoM

BeerTM

HudsonRG

PetersL

HsiehYC

2005 Improved detection of prostate cancer using classification and regression tree analysis. J Clin Oncol 23 4322 4329

23. HessKR

AbbruzzeseMC

LenziR

RaberMN

AbbruzzeseJL

1999 Classification and regression tree analysis of 1000 consecutive patients with unknown primary carcinoma. Clin Cancer Res 5 3403 3410

24. KoziolJA

ZhangJY

CasianoCA

PengXX

ShiFD

2003 Recursive partitioning as an approach to selection of immune markers for tumor diagnosis. Clin Cancer Res 9 5120 5126

25. ValeraVA

WalterBA

YokoyamaN

KoyamaY

IiaiT

2007 Prognostic groups in colorectal carcinoma patients based on tumor cell proliferation and classification and regression tree (CART) survival analysis. Ann Surg Oncol 14 34 40

26. KaplanEL

MeierP

1958 Nonparametric estimation from incomplete observations. J Am Stat Assoc 53 457 481

27. CollettD

2003 Modelling survival data in medical research, 2nd edition. Boca Raton Chapman & Hall/CRC

28. LandiMT

DrachevaT

RotunnoM

FigueroaJD

LiuH

2008 Gene expression signature of cigarette smoking and its role in lung adenocarcinoma development and survival. PLoS ONE 3 e1651 doi:10.1371/journal.pone.0001651

29. BeerDG

KardiaSL

HuangCC

GiordanoTJ

LevinAM

2002 Gene-expression profiles predict survival of patients with lung adenocarcinoma. Nat Med 8 816 824

30. EndohH

TomidaS

YatabeY

KonishiH

OsadaH

2004 Prognostic model of pulmonary adenocarcinoma by expression profiling of eight genes as determined by quantitative real-time reverse transcriptase polymerase chain reaction. J Clin Oncol 22 811 819

31. LuY

LemonW

LiuPY

YiY

MorrisonC

2006 A gene expression signature predicts survival of patients with stage I non-small cell lung cancer. PLoS Med 3 e467 doi:10.1371/journal.pmed.0030467

32. JordanVC

2009 A century of deciphering the control mechanisms of sex steroid action in breast and prostate cancer: the origins of targeted therapy and chemoprevention. Cancer Res 69 1243 1254

33. MilesSA

DezubeBJ

LeeJY

KrownSE

FletcherMA

2002 Antitumor activity of oral 9-cis-retinoic acid in HIV-associated Kaposi's sarcoma. AIDS 16 421 429

34. MiyamotoH

MessingEM

ChangC

2004 Androgen deprivation therapy for prostate cancer: current status and future prospects. Prostate 61 332 353

35. RealPJ

FerrandoAA

2009 NOTCH inhibition and glucocorticoid therapy in T-cell acute lymphoblastic leukemia. Leukemia 23 1374 1377

36. IshibashiH

SuzukiT

SuzukiS

NiikawaH

LuL

2005 Progesterone receptor in non-small cell lung cancer—a potent prognostic factor and possible target for endocrine therapy. Cancer Res 65 6450 6458

37. SiegfriedJM

2006 Hormone replacement therapy and decreased lung cancer survival. J Clin Oncol 24 9 10

38. StabileLP

DavisAL

GubishCT

HopkinsTM

LuketichJD

2002 Human non-small cell lung tumors and cells derived from normal lung express both estrogen receptor alpha and beta and show biological responses to estrogen. Cancer Res 62 2141 2150

39. StabileLP

LykerJS

GubishCT

ZhangW

GrandisJR

2005 Combined targeting of the estrogen receptor and the epidermal growth factor receptor in non-small cell lung cancer shows enhanced antiproliferative effects. Cancer Res 65 1459 1470

40. TraynorAM

SchillerJH

StabileLP

KolesarJM

EickhoffJC

2008 Pilot study of gefitinib and fulvestrant in the treatment of post-menopausal women with advanced non-small cell lung cancer. Lung Cancer 64 51 59

41. GirnunGD

NaseriE

VafaiSB

QuL

SzwayaJD

2007 Synergy between PPARgamma ligands and platinum-based drugs in cancer. Cancer Cell 11 395 406

42. ChlebowskiRT

SchwartzAG

WakeleeH

AndersonGL

StefanickML

2009 Oestrogen plus progestin and lung cancer in postmenopausal women (Women's Health Initiative trial): a post-hoc analysis of a randomised controlled trial. Lancet 374 1243 1251

43. GoodwinB

JonesSA

PriceRR

WatsonMA

McKeeDD

2000 A regulatory cascade of the nuclear receptors FXR, SHP-1, and LRH-1 represses bile acid biosynthesis. Mol Cell 6 517 526

44. LuTT

MakishimaM

RepaJJ

SchoonjansK

KerrTA

2000 Molecular basis for feedback regulation of bile acid synthesis by nuclear receptors. Mol Cell 6 507 515

45. NishizawaH

YamagataK

ShimomuraI

TakahashiM

KuriyamaH

2002 Small heterodimer partner, an orphan nuclear receptor, augments peroxisome proliferator-activated receptor gamma transactivation. J Biol Chem 277 1586 1592

46. ZhangY

XuP

ParkK

ChoiY

MooreDD

2008 Orphan receptor small heterodimer partner suppresses tumorigenesis by modulating cyclin D1 expression and cellular proliferation. Hepatology 48 289 298

47. ChandaD

ParkJH

ChoiHS

2008 Molecular basis of endocrine regulation by orphan nuclear receptor Small Heterodimer Partner. Endocr J 55 253 268

48. MullicanSE

ZhangS

KonoplevaM

RuvoloV

AndreeffM

2007 Abrogation of nuclear receptors Nr4a3 and Nr4a1 leads to development of acute myeloid leukemia. Nat Med 13 730 735

49. Di FabioF

AlvaradoC

MajdanA

GologanA

VodaL

2007 Underexpression of mineralocorticoid receptor in colorectal carcinomas and association with VEGFR-2 overexpression. J Gastrointest Surg 11 1521 1528

Štítky
Interní lékařství

Článek vyšel v časopise

PLOS Medicine


2010 Číslo 12

Nejčtenější v tomto čísle

Tomuto tématu se dále věnují…


Kurzy

Zvyšte si kvalifikaci online z pohodlí domova

Diabetická neuropatie a její léčba
nový kurz
Autoři: MUDr. Michal Dubský, Ph.D.

Mnohočetný myelom: Úvodní léčba netransplantovatelných pacientů
Autoři:

Proč v prevenci zubního kazu kartáček nestačí?
Autoři: doc. MUDr. Romana Koberová – Ivančaková, CSc.

Hypersenzitivita dentinu jako příznak poškození tvrdých zubních tkání
Autoři: MUDr. Erika Lenčová, Ph.D.

Cesta pacienta s CHOPN
Autoři: doc. MUDr. Vladimír Koblížek, Ph.D.

Všechny kurzy
Kurzy Doporučená témata Časopisy
Přihlášení
Zapomenuté heslo

Nemáte účet?  Registrujte se

Zapomenuté heslo

Zadejte e-mailovou adresu se kterou jste vytvářel(a) účet, budou Vám na ni zaslány informace k nastavení nového hesla.

Přihlášení

Nemáte účet?  Registrujte se