Gene Expression-Based Classification of Non-Small Cell Lung Carcinomas and Survival Prediction

Background

Current clinical therapy of non-small cell lung cancer depends on histo-pathological classification. This approach poorly predicts clinical outcome for individual patients. Gene expression profiling holds promise to improve clinical stratification, thus paving the way for individualized therapy.

Methodology and Principal Findings

A genome-wide gene expression analysis was performed on a cohort of 91 patients. We used 91 tumor- and 65 adjacent normal lung tissue samples. We defined sets of predictor genes (probe sets) with the expression profiles. The power of predictor genes was evaluated using an independent cohort of 96 non-small cell lung cancer- and 6 normal lung samples. We identified a tumor signature of 5 genes that aggregates the 156 tumor and normal samples into the expected groups. We also identified a histology signature of 75 genes, which classifies the samples in the major histological subtypes of non-small cell lung cancer. Correlation analysis identified 17 genes which showed the best association with post-surgery survival time. This signature was used for stratification of all patients in two risk groups. Kaplan-Meier survival curves show that the two groups display a significant difference in post-surgery survival time (p = 5.6E-6). The performance of the signatures was validated using a patient cohort of similar size (Duke University, n = 96). Compared to previously published prognostic signatures for NSCLC, the 17 gene signature performed well on these two cohorts.

Conclusions

The gene signatures identified are promising tools for histo-pathological classification of non-small cell lung cancer, and may improve the prediction of clinical outcome.     

The CIN4 Chromosomal Instability qPCR Classifier Defines Tumor Aneuploidy and Stratifies Outcome in Grade 2 Breast Cancer

Purpose

Quantifying chromosomal instability (CIN) has both prognostic and predictive clinical utility in breast cancer. In order to establish a robust and clinically applicable gene expression-based measure of CIN, we assessed the ability of four qPCR quantified genes selected from the 70-gene Chromosomal Instability (CIN70) expression signature to stratify outcome in patients with grade 2 breast cancer.

Methods

AURKA, FOXM1, TOP2A and TPX2 (CIN4), were selected from the CIN70 signature due to their high level of correlation with histological grade and mean CIN70 signature expression in silico. We assessed the ability of CIN4 to stratify outcome in an independent cohort of patients diagnosed between 1999 and 2002. 185 formalin-fixed, paraffin-embedded (FFPE) samples were included in the qPCR measurement of CIN4 expression. In parallel, ploidy status of tumors was assessed by flow cytometry. We investigated whether the categorical CIN4 score derived from the CIN4 signature was correlated with recurrence-free survival (RFS) and ploidy status in this cohort.

Results

We observed a significant association of tumor proliferation, defined by Ki67 and mitotic index (MI), with both CIN4 expression and aneuploidy. The CIN4 score stratified grade 2 carcinomas into good and poor prognostic cohorts (mean RFS: 83.8±4.9 and 69.4±8.2 months, respectively, p = 0.016) and its predictive power was confirmed by multivariate analysis outperforming MI and Ki67 expression.

Conclusions

The first clinically applicable qPCR derived measure of tumor aneuploidy from FFPE tissue, stratifies grade 2 tumors into good and poor prognosis groups.     

Identifying Prognostic Features by Bottom-Up Approach and Correlating to Drug Repositioning

Background

Traditionally top-down method was used to identify prognostic features in cancer research. That is to say, differentially expressed genes usually in cancer versus normal were identified to see if they possess survival prediction power. The problem is that prognostic features identified from one set of patient samples can rarely be transferred to other datasets. We apply bottom-up approach in this study: survival correlated or clinical stage correlated genes were selected first and prioritized by their network topology additionally, then a small set of features can be used as a prognostic signature.

Methods

Gene expression profiles of a cohort of 221 hepatocellular carcinoma (HCC) patients were used as a training set, ‘bottom-up’ approach was applied to discover gene-expression signatures associated with survival in both tumor and adjacent non-tumor tissues, and compared with ‘top-down’ approach. The results were validated in a second cohort of 82 patients which was used as a testing set.

Results

Two sets of gene signatures separately identified in tumor and adjacent non-tumor tissues by bottom-up approach were developed in the training cohort. These two signatures were associated with overall survival times of HCC patients and the robustness of each was validated in the testing set, and each predictive performance was better than gene expression signatures reported previously. Moreover, genes in these two prognosis signature gave some indications for drug-repositioning on HCC. Some approved drugs targeting these markers have the alternative indications on hepatocellular carcinoma.

Conclusion

Using the bottom-up approach, we have developed two prognostic gene signatures with a limited number of genes that associated with overall survival times of patients with HCC. Furthermore, prognostic markers in these two signatures have the potential to be therapeutic targets.     

Genome-wide DNA methylation profiles in progression to in situ and invasive carcinoma of the breast with impact on gene transcription and prognosis

Background

Ductal carcinoma in situ (DCIS) of the breast is a precursor of invasive breast carcinoma. DNA methylation alterations are thought to be an early event in progression of cancer, and may prove valuable as a tool in clinical decision making and for understanding neoplastic development.

Results

We generate genome-wide DNA methylation profiles of 285 breast tissue samples representing progression of cancer, and validate methylation changes between normal and DCIS in an independent dataset of 15 normal and 40 DCIS samples. We also validate a prognostic signature on 583 breast cancer samples from The Cancer Genome Atlas. Our analysis reveals that DNA methylation profiles of DCIS are radically altered compared to normal breast tissue, involving more than 5,000 genes. Changes between DCIS and invasive breast carcinoma involve around 1,000 genes. In tumors, DNA methylation is associated with gene expression of almost 3,000 genes, including both negative and positive correlations. A prognostic signature based on methylation level of 18 CpGs is associated with survival of breast cancer patients with invasive tumors, as well as with survival of patients with DCIS and mixed lesions of DCIS and invasive breast carcinoma.

Conclusions

This work demonstrates that changes in the epigenome occur early in the neoplastic progression, provides evidence for the possible utilization of DNA methylation-based markers of progression in the clinic, and highlights the importance of epigenetic changes in carcinogenesis.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0435-x) contains supplementary material, which is available to authorized users.     

Tumor Endothelial Inflammation Predicts Clinical Outcome in Diverse Human Cancers

Background

Vascular endothelial cells contribute to the pathogenesis of numerous human diseases by actively regulating the stromal inflammatory response; however, little is known regarding the role of endothelial inflammation in the growth of human tumors and its influence on the prognosis of human cancers.

Methods

Using an experimental model of tumor necrosis factor-alpha (TNF-α)-mediated inflammation, we characterized inflammatory gene expression in immunopurified tumor-associated endothelial cells. These genes formed the basis of a multivariate molecular predictor of overall survival that was trained and validated in four types of human cancer.

Results

We report that expression of experimentally derived tumor endothelial genes distinguished pathologic tissue specimens from normal controls in several human diseases associated with chronic inflammation. We trained these genes in human cancer datasets and defined a six-gene inflammatory signature that predicted significantly reduced overall survival in breast cancer, colon cancer, lung cancer, and glioma. This endothelial-derived signature predicted outcome independently of, but cooperatively with, standard clinical and pathological prognostic factors. Consistent with these findings, conditioned culture media from human endothelial cells stimulated by pro-inflammatory cytokines accelerated the growth of human colon and breast tumors in immunodeficient mice as compared with conditioned media from untreated endothelial cells.

Conclusions

This study provides the first prognostic cancer gene signature derived from an experimental model of tumor-associated endothelial inflammation. These findings support the notion that activation of inflammatory pathways in non-malignant tumor-infiltrating endothelial cells contributes to tumor growth and progression in multiple human cancers. Importantly, these results identify endothelial-derived factors that could serve as potential targets for therapy in diverse human cancers.     

Homogeneous datasets of triple negative breast cancers enable the identification of novel prognostic and predictive signatures

Background

Current prognostic gene signatures for breast cancer mainly reflect proliferation status and have limited value in triple-negative (TNBC) cancers. The identification of prognostic signatures from TNBC cohorts was limited in the past due to small sample sizes.

Methodology/Principal Findings

We assembled all currently publically available TNBC gene expression datasets generated on Affymetrix gene chips. Inter-laboratory variation was minimized by filtering methods for both samples and genes. Supervised analysis was performed to identify prognostic signatures from 394 cases which were subsequently tested on an independent validation cohort (n = 261 cases).

Conclusions/Significance

Using two distinct false discovery rate thresholds, 25% and <3.5%, a larger (n = 264 probesets) and a smaller (n = 26 probesets) prognostic gene sets were identified and used as prognostic predictors. Most of these genes were positively associated with poor prognosis and correlated to metagenes for inflammation and angiogenesis. No correlation to other previously published prognostic signatures (recurrence score, genomic grade index, 70-gene signature, wound response signature, 7-gene immune response module, stroma derived prognostic predictor, and a medullary like signature) was observed. In multivariate analyses in the validation cohort the two signatures showed hazard ratios of 4.03 (95% confidence interval [CI] 1.71–9.48; P = 0.001) and 4.08 (95% CI 1.79–9.28; P = 0.001), respectively. The 10-year event-free survival was 70% for the good risk and 20% for the high risk group. The 26-gene signatures had modest predictive value (AUC = 0.588) to predict response to neoadjuvant chemotherapy, however, the combination of a B-cell metagene with the prognostic signatures increased its response predictive value. We identified a 264-gene prognostic signature for TNBC which is unrelated to previously known prognostic signatures.     

Molecular Characteristics and Metastasis Predictor Genes of Triple-Negative Breast Cancer: A Clinical Study of Triple-Negative Breast Carcinomas

Background

Triple-negative breast cancer is a subtype of breast cancer with aggressive tumor behavior and distinct disease etiology. Due to the lack of an effective targeted medicine, treatment options for triple-negative breast cancer are few and recurrence rates are high. Although various multi-gene prognostic markers have been proposed for the prediction of breast cancer outcome, most of them were proven clinically useful only for estrogen receptor-positive breast cancers. Reliable identification of triple-negative patients with a favorable prognosis is not yet possible.

Methodology/Principal Findings

Clinicopathological information and microarray data from 157 invasive breast carcinomas were collected at National Taiwan University Hospital from 1995 to 2008. Gene expression data of 51 triple-negative and 106 luminal breast cancers were generated by oligonucleotide microarrays. Hierarchical clustering analysis revealed that the majority (94%) of triple-negative breast cancers were tightly clustered together carrying strong basal-like characteristics. A 45-gene prognostic signature giving 98% predictive accuracy in distant recurrence of our triple-negative patients was determined using the receiver operating characteristic analysis and leave-one-out cross validation. External validation of the prognostic signature in an independent microarray dataset of 59 early-stage triple-negative patients also obtained statistical significance (hazard ratio 2.29, 95% confidence interval (CI) 1.04–5.06, Cox P = 0.04), outperforming five other published breast cancer prognostic signatures. The 45-gene signature identified in this study revealed that TGF-β signaling of immune/inflammatory regulation may play an important role in distant metastatic invasion of triple-negative breast cancer.

Conclusions/Significance

Gene expression data and recurrence information of triple-negative breast cancer were collected and analyzed in this study. A novel set of 45-gene signature was found to be statistically predictive in disease recurrence of triple-negative breast cancer. The 45-gene signature, if further validated, may be a clinically useful tool in risk assessment of distant recurrence for early-stage triple-negative patients.     

A network-based gene expression signature informs prognosis and treatment for colorectal cancer patients

Background

Several studies have reported gene expression signatures that predict recurrence risk in stage II and III colorectal cancer (CRC) patients with minimal gene membership overlap and undefined biological relevance. The goal of this study was to investigate biological themes underlying these signatures, to infer genes of potential mechanistic importance to the CRC recurrence phenotype and to test whether accurate prognostic models can be developed using mechanistically important genes.

Methods and Findings

We investigated eight published CRC gene expression signatures and found no functional convergence in Gene Ontology enrichment analysis. Using a random walk-based approach, we integrated these signatures and publicly available somatic mutation data on a protein-protein interaction network and inferred 487 genes that were plausible candidate molecular underpinnings for the CRC recurrence phenotype. We named the list of 487 genes a NEM signature because it integrated information from Network, Expression, and Mutation. The signature showed significant enrichment in four biological processes closely related to cancer pathophysiology and provided good coverage of known oncogenes, tumor suppressors, and CRC-related signaling pathways. A NEM signature-based Survival Support Vector Machine prognostic model was trained using a microarray gene expression dataset and tested on an independent dataset. The model-based scores showed a 75.7% concordance with the real survival data and separated patients into two groups with significantly different relapse-free survival (p = 0.002). Similar results were obtained with reversed training and testing datasets (p = 0.007). Furthermore, adjuvant chemotherapy was significantly associated with prolonged survival of the high-risk patients (p = 0.006), but not beneficial to the low-risk patients (p = 0.491).

Conclusions

The NEM signature not only reflects CRC biology but also informs patient prognosis and treatment response. Thus, the network-based data integration method provides a convergence between biological relevance and clinical usefulness in gene signature development.     

Distinct Gene Expression Signatures in Lynch Syndrome and Familial Colorectal Cancer Type X

Introduction

Heredity is estimated to cause at least 20% of colorectal cancer. The hereditary nonpolyposis colorectal cancer subset is divided into Lynch syndrome and familial colorectal cancer type X (FCCTX) based on presence of mismatch repair (MMR) gene defects.

Purpose

We addressed the gene expression signatures in colorectal cancer linked to Lynch syndrome and FCCTX with the aim to identify candidate genes and to map signaling pathways relevant in hereditary colorectal carcinogenesis.

Experimental design

The 18 k whole-genome c-DNA-mediated annealing, selection, extension, and ligation (WG-DASL) assay was applied to 123 colorectal cancers, including 39 Lynch syndrome tumors and 37 FCCTX tumors. Target genes were technically validated using real-time quantitative RT-PCR (qRT-PCR) and the expression signature was validated in independent datasets.

Results

Colorectal cancers linked to Lynch syndrome and FCCTX showed distinct gene expression profiles, which by significance analysis of microarrays (SAM) differed by 2188 genes. Functional pathways involved were related to G-protein coupled receptor signaling, oxidative phosphorylation, and cell cycle function and mitosis. qRT-PCR verified altered expression of the selected genes NDUFA9, AXIN2, MYC, DNA2 and H2AFZ. Application of the 2188-gene signature to independent datasets showed strong correlation to MMR status.

Conclusion

Distinct genetic profiles and deregulation of different canonical pathways apply to Lynch syndrome and FCCTX and key targets herein may be relevant to pursue for refined diagnostic and therapeutic strategies in hereditary colorectal cancer.     

Integrated Analysis of Gene Expression Profiles Associated with Response of Platinum/Paclitaxel-Based Treatment in Epithelial Ovarian Cancer

Purpose

This study aims to explore gene expression signatures and serum biomarkers to predict intrinsic chemoresistance in epithelial ovarian cancer (EOC).

Patients and Methods

Gene expression profiling data of 322 high-grade EOC cases between 2009 and 2010 in The Cancer Genome Atlas project (TCGA) were used to develop and validate gene expression signatures that could discriminate different responses to first-line platinum/paclitaxel-based treatments. A gene regulation network was then built to further identify hub genes responsible for differential gene expression between the complete response (CR) group and the progressive disease (PD) group. Further, to find more robust serum biomarkers for clinical application, we integrated our gene signatures and gene signatures reported previously to identify secretory protein-encoding genes by searching the DAVID database. In the end, gene-drug interaction network was constructed by searching Comparative Toxicogenomics Database (CTD) and literature.

Results

A 349-gene predictive model and an 18-gene model independent of key clinical features with high accuracy were developed for prediction of chemoresistance in EOC. Among them, ten important hub genes and six critical signaling pathways were identified to have important implications in chemotherapeutic response. Further, ten potential serum biomarkers were identified for predicting chemoresistance in EOC. Finally, we suggested some drugs for individualized treatment.

Conclusion

We have developed the predictive models and serum biomarkers for platinum/paclitaxel response and established the new approach to discover potential serum biomarkers from gene expression profiles. The potential drugs that target hub genes are also suggested.     

Wnt Pathway Activity in Breast Cancer Sub-Types and Stem-Like Cells

Introduction

Wnt signalling has been implicated in stem cell regulation however its role in breast cancer stem cell regulation remains unclear.

Methods

We used a panel of normal and breast cancer cell lines to assess Wnt pathway gene and protein expression, and for the investigation of Wnt signalling within stem cell-enriched populations, mRNA and protein expression was analysed after the selection of anoikis-resistant cells. Finally, cell lines and patient-derived samples were used to investigate Wnt pathway effects on stem cell activity in vitro.

Results

Wnt pathway signalling increased in cancer compared to normal breast and in both cell lines and patient samples, expression of Wnt pathway genes correlated with estrogen receptor (ER) expression. Furthermore, specific Wnt pathway genes were predictive for recurrence within subtypes of breast cancer. Canonical Wnt pathway genes were increased in breast cancer stem cell-enriched populations in comparison to normal breast stem cell-enriched populations. Furthermore in cell lines, the ligand Wnt3a increased whilst the inhibitor DKK1 reduced mammosphere formation with the greatest inhibitory effects observed in ER+ve breast cancer cell lines. In patient-derived metastatic breast cancer samples, only ER-ve mammospheres were responsive to the ligand Wnt3a. However, the inhibitor DKK1 efficiently inhibited both ER+ve and ER-ve breast cancer but not normal mammosphere formation, suggesting that the Wnt pathway is aberrantly activated in breast cancer mammospheres.

Conclusions

Collectively, these data highlight differential Wnt signalling in breast cancer subtypes and activity in patient-derived metastatic cancer stem-like cells indicating a potential for Wnt-targeted treatment in breast cancers.     

Prognostic breast cancer signature identified from 3D culture model accurately predicts clinical outcome across independent datasets

Background

One of the major tenets in breast cancer research is that early detection is vital for patient survival by increasing treatment options. To that end, we have previously used a novel unsupervised approach to identify a set of genes whose expression predicts prognosis of breast cancer patients. The predictive genes were selected in a well-defined three dimensional (3D) cell culture model of non-malignant human mammary epithelial cell morphogenesis as down-regulated during breast epithelial cell acinar formation and cell cycle arrest. Here we examine the ability of this gene signature (3D-signature) to predict prognosis in three independent breast cancer microarray datasets having 295, 286, and 118 samples, respectively.

Methods and Findings

Our results show that the 3D-signature accurately predicts prognosis in three unrelated patient datasets. At 10 years, the probability of positive outcome was 52, 51, and 47 percent in the group with a poor-prognosis signature and 91, 75, and 71 percent in the group with a good-prognosis signature for the three datasets, respectively (Kaplan-Meier survival analysis, p<0.05). Hazard ratios for poor outcome were 5.5 (95% CI 3.0 to 12.2, p<0.0001), 2.4 (95% CI 1.6 to 3.6, p<0.0001) and 1.9 (95% CI 1.1 to 3.2, p = 0.016) and remained significant for the two larger datasets when corrected for estrogen receptor (ER) status. Hence the 3D-signature accurately predicts breast cancer outcome in both ER-positive and ER-negative tumors, though individual genes differed in their prognostic ability in the two subtypes. Genes that were prognostic in ER+ patients are AURKA, CEP55, RRM2, EPHA2, FGFBP1, and VRK1, while genes prognostic in ER− patients include ACTB, FOXM1 and SERPINE2 (Kaplan-Meier p<0.05). Multivariable Cox regression analysis in the largest dataset showed that the 3D-signature was a strong independent factor in predicting breast cancer outcome.

Conclusions

The 3D-signature accurately predicts breast cancer outcome across multiple datasets and holds prognostic value for both ER-positive and ER-negative breast cancer. The signature was selected using a novel biological approach and hence holds promise to represent the key biological processes of breast cancer.     

Combining gene signatures improves prediction of breast cancer survival

Background

Several gene sets for prediction of breast cancer survival have been derived from whole-genome mRNA expression profiles. Here, we develop a statistical framework to explore whether combination of the information from such sets may improve prediction of recurrence and breast cancer specific death in early-stage breast cancers. Microarray data from two clinically similar cohorts of breast cancer patients are used as training (n = 123) and test set (n = 81), respectively. Gene sets from eleven previously published gene signatures are included in the study.

Principal Findings

To investigate the relationship between breast cancer survival and gene expression on a particular gene set, a Cox proportional hazards model is applied using partial likelihood regression with an L2 penalty to avoid overfitting and using cross-validation to determine the penalty weight. The fitted models are applied to an independent test set to obtain a predicted risk for each individual and each gene set. Hierarchical clustering of the test individuals on the basis of the vector of predicted risks results in two clusters with distinct clinical characteristics in terms of the distribution of molecular subtypes, ER, PR status, TP53 mutation status and histological grade category, and associated with significantly different survival probabilities (recurrence: p = 0.005; breast cancer death: p = 0.014). Finally, principal components analysis of the gene signatures is used to derive combined predictors used to fit a new Cox model. This model classifies test individuals into two risk groups with distinct survival characteristics (recurrence: p = 0.003; breast cancer death: p = 0.001). The latter classifier outperforms all the individual gene signatures, as well as Cox models based on traditional clinical parameters and the Adjuvant! Online for survival prediction.

Conclusion

Combining the predictive strength of multiple gene signatures improves prediction of breast cancer survival. The presented methodology is broadly applicable to breast cancer risk assessment using any new identified gene set.     

Expression Profiling of Mitochondrial Voltage-Dependent Anion Channel-1 Associated Genes Predicts Recurrence-Free Survival in Human Carcinomas

Background

Mitochondrial voltage-dependent anion channels (VDACs) play a key role in mitochondria-mediated apoptosis. Both in vivo and in vitro evidences indicate that VDACs are actively involved in tumor progression. Specifically, VDAC-1, one member of the VDAC family, was thought to be a potential anti-cancer therapeutic target. Our previous study demonstrated that the human gene VDAC1 (encoding the VDAC-1 isoform) was significantly up-regulated in lung tumor tissue compared with normal tissue. Also, we found a significant positive correlation between the gene expression of VDAC1 and histological grade in breast cancer. However, the prognostic power of VDAC1 and its associated genes in human cancers is largely unknown.

Methods

We systematically analyzed the expression pattern of VDAC1 and its interacting genes in breast, colon, liver, lung, pancreatic, and thyroid cancers. The genes differentially expressed between normal and tumor tissues in human carcinomas were identified.

Results

The expression level of VDAC1 was uniformly up-regulated in tumor tissue compared with normal tissue in breast, colon, liver, lung, pancreatic, and thyroid cancers. Forty-four VDAC1 interacting genes were identified as being commonly differentially expressed between normal and tumor tissues in human carcinomas. We designated VDAC1 and the 44 dysregulated interacting genes as the VDAC1 associated gene signature (VAG). We demonstrate that the VAG signature is a robust prognostic biomarker to predict recurrence-free survival in breast, colon, and lung cancers, and is independent of standard clinical and pathological prognostic factors.

Conclusions

VAG represents a promising prognostic biomarker in human cancers, which may enhance prediction accuracy in identifying patients at higher risk for recurrence. Future therapies aimed specifically at VDAC1 associated genes may lead to novel agents in the treatment of cancer.