STAT1-dependent expression of energy metabolic pathways links tumour growth and radioresistance to the Warburg effect

<p>Abstract</p> <p>Background</p> <p>The Signal Transducer and Activator of Transcription 1 (STAT1) has traditionally been regarded as a transmitter of interferon signaling and a pro-apoptotic tumour suppressor. Recent data have identified new functions of STAT1 associated with tumourigenesis and resistance to genotoxic stress, including ionizing radiation (IR) and chemotherapy. To investigate the mechanisms contributing to the tumourigenic functions of STAT1, we performed a combined transcriptomic-proteomic expressional analysis and found that STAT1 is associated with regulation of energy metabolism with potential implication in the Warburg effect.</p> <p>Methods</p> <p>We generated a stable knockdown of STAT1 in the SCC61 human squamous cell carcinoma cell line, established tumour xenografts in athymic mice, and compared transcriptomic and proteomic profiles of STAT1 wild-type (WT) and knockdown (KD) untreated or irradiated (IR) tumours. Transcriptional profiling was based on Affymetrix Human GeneChip^®Gene 1.0 ST microarrays. Proteomes were determined from the tandem mass spectrometry (MS/MS) data by searching against the human subset of the UniProt database. Data were analysed using Significance Analysis of Microarrays for ribonucleic acid and Visualize software for proteins. Functional analysis was performed with Ingenuity Pathway Analysis with statistical significance measured by Fisher's exact test.</p> <p>Results</p> <p>Knockdown of STAT1 led to significant growth suppression in untreated tumours and radio sensitization of irradiated tumours. These changes were accompanied by alterations in the expression of genes and proteins of glycolysis/gluconeogenesis (GG), the citrate cycle (CC) and oxidative phosphorylation (OP). Of these pathways, GG had the most concordant changes in gene and protein expression and demonstrated a STAT1-dependent expression of genes and proteins consistent with tumour-specific glycolysis. In addition, IR drastically suppressed the GG pathway in STAT1 KD tumours without significant change in STAT1 WT tumours.</p> <p>Conclusion</p> <p>Our results identify a previously uncharacterized function of STAT1 in tumours: expressional regulation of genes encoding proteins involved in glycolysis, the citrate cycle and mitochondrial oxidative phosphorylation, with predominant regulation of glycolytic genes. STAT1-dependent expressional regulation of glycolysis suggests a potential role for STAT1 as a transcriptional modulator of genes responsible for the Warburg effect.</p

Expression Signature of IFN/STAT1 Signaling Genes Predicts Poor Survival Outcome in Glioblastoma Multiforme in a Subtype-Specific Manner

Previous reports have implicated an induction of genes in IFN/STAT1 (Interferon/STAT1) signaling in radiation resistant and prosurvival tumor phenotypes in a number of cancer cell lines, and we have hypothesized that upregulation of these genes may be predictive of poor survival outcome and/or treatment response in Glioblastoma Multiforme (GBM) patients. We have developed a list of 8 genes related to IFN/STAT1 that we hypothesize to be predictive of poor survival in GBM patients. Our working hypothesis that over-expression of this gene signature predicts poor survival outcome in GBM patients was confirmed, and in addition, it was demonstrated that the survival model was highly subtype-dependent, with strong dependence in the Proneural subtype and no detected dependence in the Classical and Mesenchymal subtypes. We developed a specific multi-gene survival model for the Proneural subtype in the TCGA (the Cancer Genome Atlas) discovery set which we have validated in the TCGA validation set. In addition, we have performed network analysis in the form of Bayesian Network discovery and Ingenuity Pathway Analysis to further dissect the underlying biology of this gene signature in the etiology of GBM. We theorize that the strong predictive value of the IFN/STAT1 gene signature in the Proneural subtype may be due to chemotherapy and/or radiation resistance induced through prolonged constitutive signaling of these genes during the course of the illness. The results of this study have implications both for better prediction models for survival outcome in GBM and for improved understanding of the underlying subtype-specific molecular mechanisms for GBM tumor progression and treatment response

MUC1-associated proliferation signature predicts outcomes in lung adenocarcinoma patients

Background: MUC1 protein is highly expressed in lung cancer. The cytoplasmic domain of MUC1 (MUC1-CD) induces tumorigenesis and resistance to DNA-damaging agents. We characterized MUC1-CD-induced transcriptional changes and examined their significance in lung cancer patients. Methods: Using DNA microarrays, we identified 254 genes that were differentially expressed in cell lines transformed by MUC1-CD compared to control cell lines. We then examined expression of these genes in 441 lung adenocarcinomas from a publicly available database. We employed statistical analyses independent of clinical outcomes, including hierarchical clustering, Student's t-tests and receiver operating characteristic (ROC) analysis, to select a seven-gene MUC1-associated proliferation signature (MAPS). We demonstrated the prognostic value of MAPS in this database using Kaplan-Meier survival analysis, log-rank tests and Cox models. The MAPS was further validated for prognostic significance in 84 lung adenocarcinoma patients from an independent database. Results: MAPS genes were found to be associated with proliferation and cell cycle regulation and included CCNB1, CDC2, CDC20, CDKN3, MAD2L1, PRC1 and RRM2. MAPS expressors (MAPS+) had inferior survival compared to non-expressors (MAPS-). In the initial data set, 5-year survival was 65% (MAPS-) vs. 45% (MAPS+, p < 0.0001). Similarly, in the validation data set, 5-year survival was 57% (MAPS-) vs. 28% (MAPS+, p = 0.005). Conclusions: The MAPS signature, comprised of MUC1-CD-dependent genes involved in the control of cell cycle and proliferation, is associated with poor outcomes in patients with adenocarcinoma of the lung. These data provide potential new prognostic biomarkers and treatment targets for lung adenocarcinoma

A Comprehensive and Universal Method for Assessing the Performance of Differential Gene Expression Analyses

The number of methods for pre-processing and analysis of gene expression data continues to increase, often making it difficult to select the most appropriate approach. We present a simple procedure for comparative estimation of a variety of methods for microarray data pre-processing and analysis. Our approach is based on the use of real microarray data in which controlled fold changes are introduced into 20% of the data to provide a metric for comparison with the unmodified data. The data modifications can be easily applied to raw data measured with any technological platform and retains all the complex structures and statistical characteristics of the real-world data. The power of the method is illustrated by its application to the quantitative comparison of different methods of normalization and analysis of microarray data. Our results demonstrate that the method of controlled modifications of real experimental data provides a simple tool for assessing the performance of data preprocessing and analysis methods

Image-Based Assessment of Growth and Signaling Changes in Cancer Cells Mediated by Direct Cell-Cell Contact

Many important biological processes are controlled through cell-cell interactions, including the colonization of metastatic tumor cells and the control of differentiation of stem cells within their niche. Despite the crucial importance of the cellular environment in regulating cellular signaling, in vitro methods for the study of such interactions are difficult and/or indirect.We report on the development of an image-based method for distinguishing two cell types grown in coculture. Furthermore, cells of one type that are in direct contact with cells of a second type (adjacent cells) can be analyzed separately from cells that are not within a single well. Changes are evaluated using population statistics, which are useful in detecting subtle changes across two populations. We have used this system to characterize changes in the LNCaP prostate carcinoma cell line when grown in contact with human vascular endothelial cells (HUVECs). We find that the expression and phosphorylation of WWOX is reduced in LNCaP cells when grown in direct contact with HUVECs. Reduced WWOX signaling has been associated with reduced activation or expression of JNK and p73. We find that p73 levels are also reduced in LNCaP cells grown in contact with HUVECs, but we did not observe such a change in JNK levels.We find that the method described is statistically robust and can be adapted to a wide variety of studies where cell function or signaling are affected by heterotypic cell-cell contact. Ironically, a potential challenge to the method is its high level of sensitivity is capable of classifying events as statistically significant (due to the high number cells evaluated individually), when the biological effect may be less clear. The methodology would be best used in conjunction with additional methods to evaluate the biological role of potentially subtle differences between populations. However, many important events, such as the establishment of a metastatic tumor, occur through rare but important changes, and methods such as we describe here can be used to identify and characterize the contribution of the environment to these changes

Identification of molecular mechanisms for cellular drug resistance by combining drug activity and gene expression profiles

Acquired drug resistance is a major problem in cancer treatment. To explore the genes involved in chemosensitivity and resistance, 10 human tumour cell lines, including parental cells and resistant subtypes selected for resistance against doxorubicin, melphalan, teniposide and vincristine, were profiled for mRNA expression of 7400 genes using cDNA microarray technology. The drug activity of 66 cancer agents was evaluated on the cell lines, and correlations between drug activity and gene expression were calculated and ranked. Hierarchical clustering of drugs based on their drug–gene correlations yielded clusters of drugs with similar mechanism of action. Genes correlated with drug sensitivity and resistance were imported into the PathwayAssist software to identify putative molecular pathways involved. A substantial number of both proapoptotic and antiapoptotic genes such as signal transducer and activator of transcription 1, mitogen-activated protein kinase 1 and focal adhesion kinase were found to be associated to drug resistance, whereas genes linked to cell cycle control and proliferation, such as cell division cycle 25A and signal transducer of activator of transcription 5A, were associated to general drug sensitivity. The results indicate that combined information from drug activity and gene expression in a resistance-based cell line panel may provide new knowledge of the genes involved in anticancer drug resistance and become a useful tool in drug development

Down-Regulation of miR-101 in Endothelial Cells Promotes Blood Vessel Formation through Reduced Repression of EZH2

Angiogenesis is a balanced process controlled by pro- and anti-angiogenic molecules of which the regulation is not fully understood. Besides classical gene regulation, miRNAs have emerged as post-transcriptional regulators of angiogenesis. Furthermore, epigenetic changes caused by histone-modifying enzymes were shown to modulate angiogenesis as well. However, a possible interplay between miRNAs and histone-modulating enzymes during angiogenesis has not been described. Here we show that VEGF-mediated down-regulation of miR-101 caused pro-angiogenic effects. We found that the pro-angiogenic effects are partly mediated through reduced repression by miR-101 of the histone-methyltransferase EZH2, a member of the Polycomb group family, thereby increasing methylation of histone H3 at lysine 27 and transcriptome alterations. In vitro, the sprouting and migratory properties of primary endothelial cell cultures were reduced by inhibiting EZH2 through up-regulation of miR-101, siRNA-mediated knockdown of EZH2, or treatment with 3-Deazaneplanocin-A (DZNep), a small molecule inhibitor of EZH2 methyltransferase activity. In addition, we found that systemic DZNep administration reduced the number of blood vessels in a subcutaneous glioblastoma mouse model, without showing adverse toxicities. Altogether, by identifying a pro-angiogenic VEGF/miR-101/EZH2 axis in endothelial cells we provide evidence for a functional link between growth factor-mediated signaling, post-transcriptional silencing, and histone-methylation in the angiogenesis process. Inhibition of EZH2 may prove therapeutic in diseases in which aberrant vascularization plays a role

Identification of novel small molecules that inhibit STAT3-dependent transcription and function

Activation of Signal Transducer and Activator of Transcription 3 (STAT3) has been linked to several processes that are critical for oncogenic transformation, cancer progression, cancer cell proliferation, survival, drug resistance and metastasis. Inhibition of STAT3 signaling has shown a striking ability to inhibit cancer cell growth and therefore, STAT3 has become a promising target for anti-cancer drug development. The aim of this study was to identify novel inhibitors of STAT-dependent gene transcription. A cellular reporter-based system for monitoring STAT3 transcriptional activity was developed which was suitable for high-throughput screening (Z’ = 0,8). This system was used to screen a library of 28,000 compounds (the ENAMINE Drug-Like Diversity Set). Following counter-screenings and toxicity studies, we identified four hit compounds that were subjected to detailed biological characterization. Of the four hits, KI16 stood out as the most promising compound, inhibiting STAT3 phosphorylation and transcriptional activity in response to IL6 stimulation. In silico docking studies showed that KI16 had favorable interactions with the STAT3 SH2 domain, however, no inhibitory activity could be observed in the STAT3 fluorescence polarization assay. KI16 inhibited cell viability preferentially in STAT3-dependent cell lines. Taken together, using a targeted, cell-based approach, novel inhibitors of STAT-driven transcriptional activity were discovered which are interesting leads to pursue further for the development of anti-cancer therapeutic agents

Co-Expression of α9β1 Integrin and VEGF-D Confers Lymphatic Metastatic Ability to a Human Breast Cancer Cell Line MDA-MB-468LN

INTRODUCTION AND OBJECTIVES: Lymphatic metastasis is a common occurrence in human breast cancer, mechanisms remaining poorly understood. MDA-MB-468LN (468LN), a variant of the MDA-MB-468GFP (468GFP) human breast cancer cell line, produces extensive lymphatic metastasis in nude mice. 468LN cells differentially express α9β1 integrin, a receptor for lymphangiogenic factors VEGF-C/-D. We explored whether (1) differential production of VEGF-C/-D by 468LN cells provides an autocrine stimulus for cellular motility by interacting with α9β1 and a paracrine stimulus for lymphangiogenesis in vitro as measured with capillary-like tube formation by human lymphatic endothelial cells (HMVEC-dLy); (2) differential expression of α9 also promotes cellular motility/invasiveness by interacting with macrophage derived factors; (3) stable knock-down of VEGF-D or α9 in 468LN cells abrogates lymphangiogenesis and lymphatic metastasis in vivo in nude mice. RESULTS: A comparison of expression of cyclo-oxygenase (COX)-2 (a VEGF-C/-D inducer), VEGF-C/-D and their receptors revealed little COX-2 expression by either cells. However, 468LN cells showed differential VEGF-D and α9β1 expression, VEGF-D secretion, proliferative, migratory/invasive capacities, latter functions being stimulated further with VEGF-D. The requirement of α9β1 for native and VEGF-D-stimulated proliferation, migration and Erk activation was demonstrated by treating with α9β1 blocking antibody or knock-down of α9. An autocrine role of VEGF-D in migration was shown by its impairment by silencing VEGF-D and restoration with VEGF-D. 468LN cells and their soluble products stimulated tube formation, migration/invasiveness of HMVEC-dLy cell in a VEGF-D dependent manner as indicated by the loss of stimulation by silencing VEGF-D in 468LN cells. Furthermore, 468LN cells showed α9-dependent stimulation of migration/invasiveness by macrophage products. Finally, capacity for intra-tumoral lymphangiogenesis and lymphatic metastasis in nude mice was completely abrogated by stable knock-down of either VEGF-D or α9 in 468LN cells. CONCLUSION: Differential capacity for VEGF-D production and α9β1 integrin expression by 468LN cells jointly contributed to their lymphatic metastatic phenotype