Steven Eschrich

Gene expression-based survival prediction in lung adenocarcinoma: A multi-site, blinded validation study

Although prognostic gene expression signatures for survival in early-stage lung cancer have been proposed, for clinical application, it is critical to establish their performance across different subject populations and in different laboratories. Here we report a large, training–testing, multi-site, blinded validation study to characterize the performance of several prognostic models based on gene expression for 442 lung adenocarcinomas. The hypotheses proposed examined whether microarray measurements of gene expression either alone or combined with basic clinical covariates (stage, age, sex) could be used to predict overall survival in lung cancer subjects. Several models examined produced risk scores that substantially correlated with actual subject outcome. Most methods performed better with clinical data, supporting the combined use of clinical and molecular information when building prognostic models for early-stage lung cancer. This study also provides the largest available set of microarray data with extensive pathological and clinical annotation for lung adenocarcinomas.

Persistent Activation of Stat3 Signaling Induces Survivin Gene Expression and Confers Resistance to Apoptosis in Human Breast Cancer Cells

Purpose: Signal transducer and activator of transcription 3 (Stat3) protein is persistently activated in breast cancer and promotes tumor cell survival. To gain a better understanding of the role of constitutive Stat3 signaling in breast cancer progression, we evaluated the expression profile of potential Stat3-regulated genes that may confer resistance to apoptosis. Experimental Design: Stat3 signaling was blocked with antisense oligonucleotides in human MDA-MB-435s breast cancer cells and Affymetrix GeneChip microarray analysis was done. The candidate Stat3 target gene Survivin was further evaluated in molecular assays using cultured breast cancer cells and immunohistochemistry of breast tumor specimens. Results: Survivin, a member of the inhibitor of apoptosis protein family, was identified as a potential Stat3-regulated gene by microarray analysis. This was confirmed in Survivin gene promoter studies and chromatin immunoprecipitation assays showing that Stat3 directly binds to and regulates the Survivin promoter. Furthermore, direct inhibition of Stat3 signaling blocked the expression of Survivin protein and induced apoptosis in breast cancer cells. Direct inhibition of Survivin expression also induced apoptosis. Increased Survivin protein expression correlates significantly (P = 0.001) with elevated Stat3 activity in primary breast tumor specimens from high-risk patients who were resistant to chemotherapy treatment. Conclusions: We identify Survivin as a direct downstream target gene of Stat3 in human breast cancer cells that is critical for their survival in culture. Our findings suggest that activated Stat3 signaling contributes to breast cancer progression and resistance to chemotherapy by, at least in part, inducing expression of the antiapoptotic protein, Survivin.

Experimentally Derived Metastasis Gene Expression Profile Predicts Recurrence and Death in Patients With Colon Cancer

BACKGROUND & AIMS Staging inadequately predicts metastatic risk in patients with colon cancer. We used a gene expression profile derived from invasive, murine colon cancer cells that were highly metastatic in an immunocompetent mouse model to identify patients with colon cancer at risk of recurrence. METHODS This phase 1, exploratory biomarker study used 55 patients with colorectal cancer from Vanderbilt Medical Center (VMC) as the training dataset and 177 patients from the Moffitt Cancer Center as the independent dataset. The metastasis-associated gene expression profile developed from the mouse model was refined with comparative functional genomics in the VMC gene expression profiles to identify a 34-gene classifier associated with high risk of metastasis and death from colon cancer. A metastasis score derived from the biologically based classifier was tested in the Moffitt dataset. RESULTS A high score was significantly associated with increased risk of metastasis and death from colon cancer across all pathologic stages and specifically in stage II and stage III patients. The metastasis score was shown to independently predict risk of cancer recurrence and death in univariate and multivariate models. For example, among stage III patients, a high score translated to increased relative risk of cancer recurrence (hazard ratio, 4.7; 95% confidence interval, 1.566-14.05). Furthermore, the metastasis score identified patients with stage III disease whose 5-year recurrence-free survival was >88% and for whom adjuvant chemotherapy did not increase survival time. CONCLUSION A gene expression profile identified from an experimental model of colon cancer metastasis predicted cancer recurrence and death, independently of conventional measures, in patients with colon cancer.

Metastasis-Associated Gene Expression Changes Predict Poor Outcomes in Patients with Dukes Stage B and C Colorectal Cancer.

Purpose: Colorectal cancer prognosis is currently predicted from pathologic staging, providing limited discrimination for Dukes stage B and C disease. Additional markers for outcome are required to help guide therapy selection for individual patients. Experimental Design: A multisite single-platform microarray study was done on 553 colorectal cancers. Gene expression changes were identified between stage A and D tumors (three training sets) and assessed as a prognosis signature in stage B and C tumors (independent test and external validation sets). Results: One hundred twenty-eight genes showed reproducible expression changes between three sets of stage A and D cancers. Using consistent genes, stage B and C cancers clustered into two groups resembling early-stage and metastatic tumors. A Prediction Analysis of Microarray algorithm was developed to classify individual intermediate-stage cancers into stage A–like/good prognosis or stage D–like/poor prognosis types. For stage B patients, the treatment adjusted hazard ratio for 6-year recurrence in individuals with stage D–like cancers was 10.3 (95% confidence interval, 1.3-80.0; P = 0.011). For stage C patients, the adjusted hazard ratio was 2.9 (95% confidence interval, 1.1-7.6; P = 0.016). Similar results were obtained for an external set of stage B and C patients. The prognosis signature was enriched for downregulated immune response genes and upregulated cell signaling and extracellular matrix genes. Accordingly, sparse tumor infiltration with mononuclear chronic inflammatory cells was associated with poor outcome in independent patients. Conclusions: Metastasis-associated gene expression changes can be used to refine traditional outcome prediction, providing a rational approach for tailoring treatments to subsets of patients. (Clin Cancer Res 2009;15(24):7642–51)

Transcriptional recapitulation and subversion of embryonic colon development by mouse colon tumor models and human colon cancer

BackgroundThe expression of carcino-embryonic antigen by colorectal cancer is an example of oncogenic activation of embryonic gene expression. Hypothesizing that oncogenesis-recapitulating-ontogenesis may represent a broad programmatic commitment, we compared gene expression patterns of human colorectal cancers (CRCs) and mouse colon tumor models to those of mouse colon development embryonic days 13.5-18.5.ResultsWe report here that 39 colon tumors from four independent mouse models and 100 human CRCs encompassing all clinical stages shared a striking recapitulation of embryonic colon gene expression. Compared to normal adult colon, all mouse and human tumors over-expressed a large cluster of genes highly enriched for functional association to the control of cell cycle progression, proliferation, and migration, including those encoding MYC, AKT2, PLK1 and SPARC. Mouse tumors positive for nuclear β-catenin shifted the shared embryonic pattern to that of early development. Human and mouse tumors differed from normal embryonic colon by their loss of expression modules enriched for tumor suppressors (EDNRB, HSPE, KIT and LSP1). Human CRC adenocarcinomas lost an additional suppressor module (IGFBP4, MAP4K1, PDGFRA, STAB1 and WNT4). Many human tumor samples also gained expression of a coordinately regulated module associated with advanced malignancy (ABCC1, FOXO3A, LIF, PIK3R1, PRNP, TNC, TIMP3 and VEGF).ConclusionCross-species, developmental, and multi-model gene expression patterning comparisons provide an integrated and versatile framework for definition of transcriptional programs associated with oncogenesis. This approach also provides a general method for identifying pattern-specific biomarkers and therapeutic targets. This delineation and categorization of developmental and non-developmental activator and suppressor gene modules can thus facilitate the formulation of sophisticated hypotheses to evaluate potential synergistic effects of targeting within- and between-modules for next-generation combinatorial therapeutics and improved mouse models.

Multi-Platform, Multi-Site, Microarray-Based Human Tumor Classification

The introduction of gene expression profiling has resulted in the production of rich human data sets with potential for deciphering tumor diagnosis, prognosis, and therapy. Here we demonstrate how artificial neural networks (ANNs) can be applied to two completely different microarray platforms (cDNA and oligonucleotide), or a combination of both, to build tumor classifiers capable of deciphering the identity of most human cancers. First, 78 tumors representing eight different types of histologically similar adenocarcinoma, were evaluated with a 32k cDNA microarray and correctly classified by a cDNA-based ANN, using independent training and test sets, with a mean accuracy of 83%. To expand our approach, oligonucleotide data derived from six independent performance sites, representing 463 tumors and 21 tumor types, were assembled, normalized, and scaled. An oligonucleotide-based ANN, trained on a random fraction of the tumors (n = 343), was 88% accurate in predicting known pathological origin of the remaining fraction of tumors (n = 120) not exposed to the training algorithm. Finally, a mixed-platform classifier using a combination of both cDNA and oligonucleotide microarray data from seven performance sites, normalized and scaled from a large and diverse tumor set (n = 539), produced similar results (85% accuracy) on independent test sets. Further validation of our classifiers was achieved by accurately (84%) predicting the known primary site of origin for an independent set of metastatic lesions (n = 50), resected from brain, lung, and liver, potentially addressing the vexing classification problems imposed by unknown primary cancers. These cDNA- and oligonucleotide-based classifiers provide a first proof of principle that data derived from multiple platforms and performance sites can be exploited to build multi-tissue tumor classifiers.

A chemical and phosphoproteomic characterization of dasatinib action in lung cancer

We describe a strategy to comprehend signaling pathways active in lung cancer cells and targeted by dasatinib employing chemical proteomics to identify direct interacting proteins combined with immunoaffinity purification of tyrosine phosphorylated peptides corresponding to activated tyrosine kinases. We identified nearly 40 different kinase targets of dasatinib. These include SFK members (LYN, SRC, FYN, LCK, YES), non-receptor tyrosine kinases (FRK, BRK, ACK), and receptor tyrosine kinases (Ephrin receptors, DDR1, EGFR). Using quantitative phosphoproteomics we identified peptides corresponding to autophosphorylation sites of these tyrosine kinases that are inhibited in a concentration-dependent manner by dasatinib. Using drug resistant gatekeeper mutants, we show that SFK kinases, particularly SRC and FYN, as well as EGFR are relevant targets for dasatinib action. The combined mass spectrometry based approach described here provides a system-level view of dasatinib action in cancer cells and suggests both functional targets and rationale combinatorial therapeutic strategies.

12-Chemokine Gene Signature Identifies Lymph Node-like Structures in Melanoma: Potential for Patient Selection for Immunotherapy?

We have interrogated a 12-chemokine gene expression signature (GES) on genomic arrays of 14,492 distinct solid tumors and show broad distribution across different histologies. We hypothesized that this 12-chemokine GES might accurately predict a unique intratumoral immune reaction in stage IV (non-locoregional) melanoma metastases. The 12-chemokine GES predicted the presence of unique, lymph node-like structures, containing CD20+ B cell follicles with prominent areas of CD3+ T cells (both CD4+ and CD8+ subsets). CD86+, but not FoxP3+, cells were present within these unique structures as well. The direct correlation between the 12-chemokine GES score and the presence of unique, lymph nodal structures was also associated with better overall survival of the subset of melanoma patients. The use of this novel 12-chemokine GES may reveal basic information on in situ mechanisms of the anti-tumor immune response, potentially leading to improvements in the identification and selection of melanoma patients most suitable for immunotherapy.

A Gene Expression Model of Intrinsic Tumor Radiosensitivity: Prediction of Response and Prognosis After Chemoradiation

PURPOSE Development of a radiosensitivity predictive assay is a central goal of radiation oncology. We reasoned a gene expression model could be developed to predict intrinsic radiosensitivity and treatment response in patients. METHODS AND MATERIALS Radiosensitivity (determined by survival fraction at 2 Gy) was modeled as a function of gene expression, tissue of origin, ras status (mut/wt), and p53 status (mut/wt) in 48 human cancer cell lines. Ten genes were identified and used to build a rank-based linear regression algorithm to predict an intrinsic radiosensitivity index (RSI, high index = radioresistance). This model was applied to three independent cohorts treated with concurrent chemoradiation: head-and-neck cancer (HNC, n = 92); rectal cancer (n = 14); and esophageal cancer (n = 12). RESULTS Predicted RSI was significantly different in responders (R) vs. nonresponders (NR) in the rectal (RSI R vs. NR 0.32 vs. 0.46, p = 0.03), esophageal (RSI R vs. NR 0.37 vs. 0.50, p = 0.05) and combined rectal/esophageal (RSI R vs. NR 0.34 vs. 0.48, p = 0.001511) cohorts. Using a threshold RSI of 0.46, the model has a sensitivity of 80%, specificity of 82%, and positive predictive value of 86%. Finally, we evaluated the model as a prognostic marker in HNC. There was an improved 2-year locoregional control (LRC) in the predicted radiosensitive group (2-year LRC 86% vs. 61%, p = 0.05). CONCLUSIONS We validate a robust multigene expression model of intrinsic tumor radiosensitivity in three independent cohorts totaling 118 patients. To our knowledge, this is the first time that a systems biology-based radiosensitivity model is validated in multiple independent clinical datasets.

Prediction of Radiation Sensitivity Using a Gene Expression Classifier

The development of a successful radiation sensitivity predictive assay has been a major goal of radiation biology for several decades. We have developed a radiation classifier that predicts the inherent radiosensitivity of tumor cell lines as measured by survival fraction at 2 Gy (SF2), based on gene expression profiles obtained from the literature. Our classifier correctly predicts the SF2 value in 22 of 35 cell lines from the National Cancer Institute panel of 60, a result significantly different from chance (P = 0.0002). In our approach, we treat radiation sensitivity as a continuous variable, significance analysis of microarrays is used for gene selection, and a multivariate linear regression model is used for radiosensitivity prediction. The gene selection step identified three novel genes (RbAp48, RGS19, and R5PIA) of which expression values are correlated with radiation sensitivity. Gene expression was confirmed by quantitative real-time PCR. To biologically validate our classifier, we transfected RbAp48 into three cancer cell lines (HS-578T, MALME-3M, and MDA-MB-231). RbAp48 overexpression induced radiosensitization (1.5- to 2-fold) when compared with mock-transfected cell lines. Furthermore, we show that HS-578T-RbAp48 overexpressors have a higher proportion of cells in G2-M (27% versus 5%), the radiosensitive phase of the cell cycle. Finally, RbAp48 overexpression is correlated with dephosphorylation of Akt, suggesting that RbAp48 may be exerting its effect by antagonizing the Ras pathway. The implications of our findings are significant. We establish that radiation sensitivity can be predicted based on gene expression profiles and we introduce a genomic approach to the identification of novel molecular markers of radiation sensitivity.