Carboplatin-induced gene expression changes in vitro are prognostic of survival in epithelial ovarian cancer

<p>Abstract</p> <p>Background</p> <p>We performed a time-course microarray experiment to define the transcriptional response to carboplatin <it>in vitro</it>, and to correlate this with clinical outcome in epithelial ovarian cancer (EOC). RNA was isolated from carboplatin and control-treated 36M2 ovarian cancer cells at several time points, followed by oligonucleotide microarray hybridization. Carboplatin induced changes in gene expression were assessed at the single gene as well as at the pathway level. Clinical validation was performed in publicly available microarray datasets using disease free and overall survival endpoints.</p> <p>Results</p> <p>Time-course and pathway analyses identified 317 genes and 40 pathways (designated time-course and pathway signatures) deregulated following carboplatin exposure. Both types of signatures were validated in two separate platinum-treated ovarian and NSCLC cell lines using published microarray data. Expression of time-course and pathway signature genes distinguished between patients with unfavorable and favorable survival in two independent ovarian cancer datasets. Among the pathways most highly induced by carboplatin <it>in vitro</it>, the NRF2, NF-kB, and cytokine and inflammatory response pathways were also found to be upregulated prior to chemotherapy exposure in poor prognosis tumors.</p> <p>Conclusion</p> <p>Dynamic assessment of gene expression following carboplatin exposure <it>in vitro </it>can identify both genes and pathways that are correlated with clinical outcome. The functional relevance of this observation for better understanding the mechanisms of drug resistance in EOC will require further evaluation.</p

Analysis of Multiple Sarcoma Expression Datasets: Implications for Classification, Oncogenic Pathway Activation and Chemotherapy Resistance

Background: Diagnosis of soft tissue sarcomas (STS) is challenging. Many remain unclassified (not-otherwise-specified, NOS) or grouped in controversial categories such as malignant fibrous histiocytoma (MFH), with unclear therapeutic value. We analyzed several independent microarray datasets, to identify a predictor, use it to classify unclassifiable sarcomas, and assess oncogenic pathway activation and chemotherapy response. Methodology/Principal Findings: We analyzed 5 independent datasets (325 tumor arrays). We developed and validated a predictor, which was used to reclassify MFH and NOS sarcomas. The molecular “match” between MFH and their predicted subtypes was assessed using genome-wide hierarchical clustering and Subclass-Mapping. Findings were validated in 15 paraffin samples profiled on the DASL platform. Bayesian models of oncogenic pathway activation and chemotherapy response were applied to individual STS samples. A 170-gene predictor was developed and independently validated (80-85% accuracy in all datasets). Most MFH and NOS tumors were reclassified as leiomyosarcomas, liposarcomas and fibrosarcomas. “Molecular match” between MFH and their predicted STS subtypes was confirmed both within and across datasets. This classification revealed previously unrecognized tissue differentiation lines (adipocyte, fibroblastic, smooth-muscle) and was reproduced in paraffin specimens. Different sarcoma subtypes demonstrated distinct oncogenic pathway activation patterns, and reclassified MFH tumors shared oncogenic pathway activation patterns with their predicted subtypes. These patterns were associated with predicted resistance to chemotherapeutic agents commonly used in sarcomas. Conclusions/Significance: STS profiling can aid in diagnosis through a predictor tracking distinct tissue differentiation in unclassified tumors, and in therapeutic management via oncogenic pathway activation and chemotherapy response assessment

Proteomic Identification of Interleukin-2 Therapy Response in Metastatic Renal Cell Cancer

Introduction—To detect a predictive protein profile that distinguishes between IL-2 therapy responders and non-responders among metastatic RCC patients we used surface-enhanced laser desorption/ionization time-of-flight mass spectrometry (SELDI TOF-MS). Materials and Methods—Protein extracts of 56 metastatic clear cell RCC patients obtained from radical nephrectomy specimens and prior to IL-2 therapy were applied to protein chip arrays of different chromatographic properties and analyzed using SELDI TOF-MS. A class prediction algorithm was applied to identify a subset of protein peaks whose expression values were associated with IL-2 response status. Multivariate analysis was performed to assess the association between the proteomic profile and the IL-2 response status controlling for the effect of lymphadenopathy. Results—From a total of 513 protein peaks we discovered a predictor set of 11 peaks that performed optimally for predicting IL-2 response status (86 % accuracy, Fisher’s p\u3c0.004, permutation p\u3c0.01). The results were validated on an independent data set with an overall accuracy of 72% (p \u3c 0.05, permutation p\u3c0.01). On multivariate analysis the proteomic profile was significantly associated with IL-2 response when corrected for lymph node status (p\u3c 0.04). Conclusions—We have identified and validated a proteomic pattern that is an independent predictor of IL-2 response. The ability to predict the probability of IL-2 response could permit targeted selection of patients most likely to respond to IL-2, while avoiding unwanted toxicities in patients less likely to respond. This proteomic predictor has the potential to significantly aid clinicians in the decision making of appropriate therapy for metastatic RCC patients

Integrated Analysis of Multiple Microarray Datasets Identifies a Reproducible Survival Predictor in Ovarian Cancer

BACKGROUND: Public data integration may help overcome challenges in clinical implementation of microarray profiles. We integrated several ovarian cancer datasets to identify a reproducible predictor of survival. METHODOLOGY/PRINCIPAL FINDINGS: Four microarray datasets from different institutions comprising 265 advanced stage tumors were uniformly reprocessed into a single training dataset, also adjusting for inter-laboratory variation ("batch-effect"). Supervised principal component survival analysis was employed to identify prognostic models. Models were independently validated in a 61-patient cohort using a custom array genechip and a publicly available 229-array dataset. Molecular correspondence of high- and low-risk outcome groups between training and validation datasets was demonstrated using Subclass Mapping. Previously established molecular phenotypes in the 2(nd) validation set were correlated with high and low-risk outcome groups. Functional representational and pathway analysis was used to explore gene networks associated with high and low risk phenotypes. A 19-gene model showed optimal performance in the training set (median OS 31 and 78 months, p < 0.01), 1(st) validation set (median OS 32 months versus not-yet-reached, p = 0.026) and 2(nd) validation set (median OS 43 versus 61 months, p = 0.013) maintaining independent prognostic power in multivariate analysis. There was strong molecular correspondence of the respective high- and low-risk tumors between training and 1(st) validation set. Low and high-risk tumors were enriched for favorable and unfavorable molecular subtypes and pathways, previously defined in the public 2(nd) validation set. CONCLUSIONS/SIGNIFICANCE: Integration of previously generated cancer microarray datasets may lead to robust and widely applicable survival predictors. These predictors are not simply a compilation of prognostic genes but appear to track true molecular phenotypes of good- and poor-outcome

MicroRNA paraffin-based studies in osteosarcoma reveal reproducible independent prognostic profiles at 14q32

Background: Although microRNAs (miRNAs) are implicated in osteosarcoma biology and chemoresponse, miRNA prognostic models are still needed, particularly because prognosis is imperfectly correlated with chemoresponse. Formalin-fixed, paraffin-embedded tissue is a necessary resource for biomarker studies in this malignancy with limited frozen tissue availability. Methods: We performed miRNA and mRNA microarray formalin-fixed, paraffin-embedded assays in 65 osteosarcoma biopsy and 26 paired post-chemotherapy resection specimens and used the only publicly available miRNA dataset, generated independently by another group, to externally validate our strongest findings (n = 29). We used supervised principal components analysis and logistic regression for survival and chemoresponse, and miRNA activity and target gene set analysis to study miRNA regulatory activity. Results: Several miRNA-based models with as few as five miRNAs were prognostic independently of pathologically assessed chemoresponse (median recurrence-free survival: 59 months versus not-yet-reached; adjusted hazards ratio = 2.90; P = 0.036). The independent dataset supported the reproducibility of recurrence and survival findings. The prognostic value of the profile was independent of confounding by known prognostic variables, including chemoresponse, tumor location and metastasis at diagnosis. Model performance improved when chemoresponse was added as a covariate (median recurrence-free survival: 59 months versus not-yet-reached; hazard ratio = 3.91; P = 0.002). Most prognostic miRNAs were located at 14q32 - a locus already linked to osteosarcoma - and their gene targets display deregulation patterns associated with outcome. We also identified miRNA profiles predictive of chemoresponse (75% to 80% accuracy), which did not overlap with prognostic profiles. Conclusions: Formalin-fixed, paraffin-embedded tissue-derived miRNA patterns are a powerful prognostic tool for risk-stratified osteosarcoma management strategies. Combined miRNA and mRNA analysis supports a possible role of the 14q32 locus in osteosarcoma progression and outcome. Our study creates a paradigm for formalin-fixed, paraffin-embedded-based miRNA biomarker studies in cancer

EWSR1-ATF1 dependent 3D connectivity regulates oncogenic and differentiation programs in Clear Cell Sarcoma.

Oncogenic fusion proteins generated by chromosomal translocations play major roles in cancer. Among them, fusions between EWSR1 and transcription factors generate oncogenes with powerful chromatin regulatory activities, capable of establishing complex gene expression programs in permissive precursor cells. Here we define the epigenetic and 3D connectivity landscape of Clear Cell Sarcoma, an aggressive cancer driven by the EWSR1-ATF1 fusion gene. We find that EWSR1-ATF1 displays a distinct DNA binding pattern that requires the EWSR1 domain and promotes ATF1 retargeting to new distal sites, leading to chromatin activation and the establishment of a 3D network that controls oncogenic and differentiation signatures observed in primary CCS tumors. Conversely, EWSR1-ATF1 depletion results in a marked reconfiguration of 3D connectivity, including the emergence of regulatory circuits that promote neural crest-related developmental programs. Taken together, our study elucidates the epigenetic mechanisms utilized by EWSR1-ATF1 to establish regulatory networks in CCS, and points to precursor cells in the neural crest lineage as candidate cells of origin for these tumors

Identifying microRNA/mRNA dysregulations in ovarian cancer

Abstract Background MicroRNAs are a class of noncoding RNA molecules that co-regulate the expression of multiple genes via mRNA transcript degradation or translation inhibition. Since they often target entire pathways, they may be better drug targets than genes or proteins. MicroRNAs are known to be dysregulated in many tumours and associated with aggressive or poor prognosis phenotypes. Since they regulate mRNA in a tissue specific manner, their functional mRNA targets are poorly understood. In previous work, we developed a method to identify direct mRNA targets of microRNA using patient matched microRNA/mRNA expression data using an anti-correlation signature. This method, applied to clear cell Renal Cell Carcinoma (ccRCC), revealed many new regulatory pathways compromised in ccRCC. In the present paper, we apply this method to identify dysregulated microRNA/mRNA mechanisms in ovarian cancer using data from The Cancer Genome Atlas (TCGA). Methods TCGA Microarray data was normalized and samples whose class labels (tumour or normal) were ambiguous with respect to consensus ensemble K-Means clustering were removed. Significantly anti-correlated and correlated genes/microRNA differentially expressed between tumour and normal samples were identified. TargetScan was used to identify gene targets of microRNA. Results We identified novel microRNA/mRNA mechanisms in ovarian cancer. For example, the expression level of RAD51AP1 was found to be strongly anti-correlated with the expression of hsa-miR-140-3p, which was significantly down-regulated in the tumour samples. The anti-correlation signature was present separately in the tumour and normal samples, suggesting a direct causal dysregulation of RAD51AP1 by hsa-miR-140-3p in the ovary. Other pairs of potentially biological relevance include: hsa-miR-145/E2F3, hsa-miR-139-5p/TOP2A, and hsa-miR-133a/GCLC. We also identified sets of positively correlated microRNA/mRNA pairs that are most likely result from indirect regulatory mechanisms. Conclusions Our findings identify novel microRNA/mRNA relationships that can be verified experimentally. We identify both generic microRNA/mRNA regulation mechanisms in the ovary as well as specific microRNA/mRNA controls which are turned on or off in ovarian tumours. Our results suggest that the disease process uses specific mechanisms which may be significant for their utility as early detection biomarkers or in the development of microRNA therapies in treating ovarian cancers. The positively correlated microRNA/mRNA pairs suggest the existence of novel regulatory mechanisms that proceed via intermediate states (indirect regulation) in ovarian tumorigenesis.</p

C-Fos expression is a molecular predictor of progression and survival in epithelial ovarian carcinoma

Members of the Fos protein family dimerise with Jun proteins to form the AP-1 transcription factor complex. They have a central function in proliferation and differentiation of normal tissue as well as in oncogenic transformation and tumour progression. We analysed the expression of c-Fos, FosB, Fra-1 and Fra-2 to investigate the function of Fos transcription factors in ovarian cancer. A total of 101 patients were included in the study. Expression of Fos proteins was determined by western blot analysis, quantified by densitometry and verified by immunohistochemistry. Reduced c-Fos expression was independently associated with unfavourable progression-free survival (20.6, 31.6 and 51.2 months for patients with low, moderate and high c-Fos expression; P=0.003) as well as overall survival (23.8, 46.0 and 55.5 months for low, moderate and high c-Fos levels; P=0.003). No correlations were observed for FosB, Fra-1 and Fra-2. We conclude that loss of c-Fos expression is associated with tumour progression in ovarian carcinoma and that c-Fos may be a prognostic factor. These results are in contrast to the classic concept of c-Fos as an oncogene, but are supported by the recently discovered tumour-suppressing and proapoptotic function of c-Fos in various cancer types

Expression profiling identifies genes involved in neoplastic transformation of serous ovarian cancer

Background: The malignant potential of serous ovarian tumors, the most common ovarian tumor subtype, varies from benign to low malignant potential (LMP) tumors to frankly invasive cancers. Given the uncertainty about the relationship between these different forms, we compared their patterns of gene expression. Methods: Expression profiling was carried out on samples of 7 benign, 7 LMP and 28 invasive (moderate and poorly differentiated) serous tumors and four whole normal ovaries using oligonucleotide microarrays representing over 21,000 genes. Results: We identified 311 transcripts that distinguished invasive from benign tumors, and 20 transcripts that were significantly differentially expressed between invasive and LMP tumors at p < 0.01 (with multiple testing correction). Five genes that were differentially expressed between invasive and either benign or normal tissues were validated by real time PCR in an independent panel of 46 serous tumors (4 benign, 7 LMP, 35 invasive). Overexpression of SLPI and WNT7A and down-regulation of C6orf31, PDGFRA and GLTSCR2 were measured in invasive and LMP compared with benign and normal tissues. Over-expression of WNT7A in an ovarian cancer cell line led to increased migration and invasive capacity. Conclusion: These results highlight several genes that may play an important role across the spectrum of serous ovarian tumorigenesis

Evaluation of Candidate Stromal Epithelial Cross-Talk Genes Identifies Association between Risk of Serous Ovarian Cancer and TERT, a Cancer Susceptibility “Hot-Spot”

We hypothesized that variants in genes expressed as a consequence of interactions between ovarian cancer cells and the host micro-environment could contribute to cancer susceptibility. We therefore used a two-stage approach to evaluate common single nucleotide polymorphisms (SNPs) in 173 genes involved in stromal epithelial interactions in the Ovarian Cancer Association Consortium (OCAC). In the discovery stage, cases with epithelial ovarian cancer (n = 675) and controls (n = 1,162) were genotyped at 1,536 SNPs using an Illumina GoldenGate assay. Based on Positive Predictive Value estimates, three SNPs—PODXL rs1013368, ITGA6 rs13027811, and MMP3 rs522616—were selected for replication using TaqMan genotyping in up to 3,059 serous invasive cases and 8,905 controls from 16 OCAC case-control studies. An additional 18 SNPs with Pper-allele<0.05 in the discovery stage were selected for replication in a subset of five OCAC studies (n = 1,233 serous invasive cases; n = 3,364 controls). The discovery stage associations in PODXL, ITGA6, and MMP3 were attenuated in the larger replication set (adj. Pper-allele≥0.5). However genotypes at TERT rs7726159 were associated with ovarian cancer risk in the smaller, five-study replication study (Pper-allele = 0.03). Combined analysis of the discovery and replication sets for this TERT SNP showed an increased risk of serous ovarian cancer among non-Hispanic whites [adj. ORper-allele 1.14 (1.04–1.24) p = 0.003]. Our study adds to the growing evidence that, like the 8q24 locus, the telomerase reverse transcriptase locus at 5p15.33, is a general cancer susceptibility locus