Large-scale molecular annotation of epithelial ovarian cancer (EOC) indicates remarkable heterogeneity in the etiology of that disease. This diversity presents a significant obstacle against intervention target discovery. However, inactivation of miRNA biogenesis is commonly associated with advanced disease. Thus, restoration of miRNA activity may represent a common vulnerability among diverse EOC oncogenotypes. To test this, we employed genome-scale, gain-of-function, miRNA mimic toxicity screens in a large, diverse spectrum of EOC cell lines. We found that all cell lines responded to at least some miRNA mimics, but that the nature of the miRNA mimics provoking a response was highly selective within the panel. These selective toxicity profiles were leveraged to define modes of action and molecular response indicators for miRNA mimics with tumor-suppressive characteristics in vivo. A mechanistic principle emerging from this analysis was sensitivity of EOC to miRNA-mediated release of cell fate specification programs, loss of which may be a prerequisite for development of this disease.

Rhabdoid histology in clear-cell renal cell carcinoma is associated with a poor prognosis. The prognosis of patients with clear-cell renal cell carcinoma may also be influenced by molecular alterations. The aim of this study was to evaluate the association between histologic features and salient molecular changes in rhabdoid clear-cell renal cell carcinoma. We macrodissected the rhabdoid and clear-cell epithelioid components from 12 cases of rhabdoid clear-cell renal cell carcinoma. We assessed cancer-related mutations from eight cases using a clinical next-generation exome-sequencing platform. The transcriptome of rhabdoid clear-cell renal cell carcinoma (n=8) and non-rhabdoid clear-cell renal cell carcinoma (n=37) was assessed by RNA-seq and gene expression microarray. VHL (63%) showed identical mutations in all regions from the same tumor. BAP1 (38%) and PBRM1 (13%) mutations were identified in the rhabdoid but not in the epithelioid component and were mutually exclusive in 3/3 cases and 1 case, respectively. SETD2 (63%) mutations were discordant between different histologic regions in 2/5 cases, with mutations called only in the epithelioid and rhabdoid components, respectively. The transcriptome of rhabdoid clear-cell renal cell carcinoma was distinct from advanced-stage and high-grade clear-cell renal cell carcinoma. The diverse histologic components of rhabdoid clear-cell renal cell carcinoma, however, showed a similar transcriptomic program, including a similar prognostic gene expression signature. Rhabdoid clear-cell renal cell carcinoma is transcriptomically distinct and shows a high rate of SETD2 and BAP1 mutations and a low rate of PBRM1 mutations. Driver mutations in clear-cell renal cell carcinoma are often discordant across different morphologic regions, whereas the gene expression program is relatively stable. Molecular profiling of clear-cell renal cell carcinoma may improve by assessing for gene expression and sampling tumor foci from different histologic regions.

Because of its high rate of metastasis, inflammatory breast cancer (IBC) has a poor prognosis compared with non-inflammatory types of breast cancer (non-IBC). In a recent study, Lehmann and colleagues identified seven subtypes of triple-negative breast cancer (TNBC). We hypothesized that the distribution of TNBC subtypes differs between TN-IBC and TN-non-IBC. We determined the subtypes and compared clinical outcomes by subtype in TN-IBC and TN-non-IBC patients.

Cancer cells require both migratory and tumorigenic property to establish metastatic tumors outside the primary microenvironment. Identifying the characteristic features of migratory cancer stem cells with tumorigenic property is important to predict patient prognosis and combat metastasis. Here we established one epithelial and two mesenchymal cell lines from ascites of a bladder cancer patient (i.e. cells already migrated outside primary tumor). Analyses of these cell lines demonstrated that the epithelial cells with surface expression of PD-L1, E-cadherin, CD24, and VEGFR2 rapidly formed tumors outside the primary tumor microenvironment in nude mice, exhibited signatures of immune evasion, increased stemness, increased calcium signaling, transformation, and novel E-cadherin-RalBP1 interaction. The mesenchymal cells on the other hand, exhibited constitutive TGF-β signaling and were less tumorigenic. Hence, targeting epithelial cancer stem cells with rapid tumorigenesis signatures in future might help to combat metastasis.

A deeper mechanistic understanding of tumour angiogenesis regulation is needed to improve current anti-angiogenic therapies. Here we present evidence from systems-based miRNA analyses of large-scale patient data sets along with in vitro and in vivo experiments that miR-192 is a key regulator of angiogenesis. The potent anti-angiogenic effect of miR-192 stems from its ability to globally downregulate angiogenic pathways in cancer cells through regulation of EGR1 and HOXB9. Low miR-192 expression in human tumours is predictive of poor clinical outcome in several cancer types. Using 1,2-dioleoyl-sn-glycero-3-phosphatidylcholine (DOPC) nanoliposomes, we show that miR-192 delivery leads to inhibition of tumour angiogenesis in multiple ovarian and renal tumour models, resulting in tumour regression and growth inhibition. This anti-angiogenic and anti-tumour effect is more robust than that observed with an anti-VEGF antibody. Collectively, these data identify miR-192 as a central node in tumour angiogenesis and support the use of miR-192 in an anti-angiogenesis therapy.

Two major identifiable sources of variation in data derived from the Serial Analysis of Gene Expression (SAGE) are within-library sampling variability and between-library heterogeneity within a group. Most published methods for identifying differential expression focus on just the sampling variability. In recent work, the problem of assessing differential expression between two groups of SAGE libraries has been addressed by introducing a beta-binomial hierarchical model that explicitly deals with both of the above sources of variation. This model leads to a test statistic analogous to a weighted two-sample t-test. When the number of groups involved is more than two, however, a more general approach is needed.

Anti-angiogenesis therapy has shown clinical benefit in patients with high-grade serous ovarian cancer (HGSC), but adaptive resistance rapidly emerges. Thus, approaches to overcome such resistance are needed. We developed the setting of adaptive resistance to anti-VEGF therapy, and performed a series of in vivo experiments in both immune competent and nude mouse models. Given the pro-angiogenic properties of tumor-associated macrophages (TAMs) and the dominant role of CSF1R in macrophage function, we added CSF1R inhibitors following emergence of adaptive resistance to anti-VEGF antibody. Mice treated with a CSF1R inhibitor (AC708) after anti-VEGF antibody resistance had little to no measurable tumor burden upon completion of the experiment while those that did not receive a CSF1R inhibitor still had abundant tumor. To mimic clinically used regimens, mice were also treated with anti-VEGF antibody and paclitaxel until resistance emerged, and then a CSF1R inhibitor was added. The addition of a CSF1R inhibitor restored response to anti-angiogenesis therapy, resulting in 83% lower tumor burden compared to treatment with anti-VEGF antibody and paclitaxel alone. Collectively, our data demonstrate that the addition of a CSF1R inhibitor to anti-VEGF therapy and taxane chemotherapy results in robust anti-tumor effects.

Investigations into the function of nonpromoter DNA methylation have yielded new insights into epigenetic regulation of gene expression. Previous studies have highlighted the importance of distinguishing between DNA methylation in discrete functional regions; however, integrated nonpromoter DNA methylation and gene expression analyses across a wide number of tumor types and corresponding normal tissues have not been performed. Through integrated analysis of gene expression and DNA methylation profiles, we examined 32 tumor types and identified 57 tumor suppressors and oncogenes out of 260 genes exhibiting a correlation of > 0.5 between gene body methylation and gene expression in at least one tumor type. The lymphocyte-specific gene CARD11 exhibits robust association between gene body methylation and expression across 19 of 32 tumor types examined. It is significantly overexpressed in kidney renal cell carcinoma (KIRC) and lung adenocarcinoma (LUAD) tumor tissues in comparison with respective control samples; and is significantly associated with lower overall survival in KIRC. Contrary to its canonical function in lymphocyte NFκB activation, CARD11 activates the mTOR pathway in KIRC and LUAD, resulting in suppressed autophagy. Furthermore, demethylation of a CpG island within the gene body of CARD11 decreases gene expression. Collectively, our study highlights how DNA methylation outside the promoter region can impact tumor progression. IMPLICATIONS: Our study describes a novel regulatory role of gene body DNA methylation-dependent CARD11 expression on mTOR signaling and its impact on tumor progression.

Despite having similar histologic features, patients with high-grade serous ovarian carcinoma (HGSC) often experience highly variable outcomes. The underlying determinants for long-term survival (LTS, ≥10 years) versus short-term survival (STS, <3 years) are largely unknown. The present study sought to identify molecular predictors of LTS for women with HGSC. A cohort of 24 frozen HGSC samples was collected (12 LTS and 12 STS) and analyzed at DNA, RNA, and protein levels. OVCAR5 and OVCAR8 cell lines were used for in vitro validation studies. For in vivo studies, we injected OVCAR8 cells into the peritoneal cavity of female athymic nude mice. From RNAseq analysis, 11 genes were found to be differentially expressed between the STS and LTS groups (fold change > 2; false discovery rate < 0.01). In the subsequent validation cohort, transmembrane protein 62 (TMEM62) was found to be related to LTS. CIBERSORT analysis showed that T cells (follicular helper) were found at higher levels in tumors from LTS than STS groups. In vitro data using OVCAR5 and OVCAR8 cells showed decreased proliferation with TMEM62 overexpression and positive correlation with a longevity-regulating pathway (KEGG HSA04213) at the RNA level. In vivo analysis using the OVCAR8-TMEM62-TetON model showed decreased tumor burden in mice with high- vs. low-expressing TMEM62 tumors. Our results demonstrate that restoring TMEM62 may be a novel approach for treatment of HGSC. These findings may have implications for biomarker and intervention strategies to help improve patient outcomes

Sarcomatoid transformation, wherein an epithelioid carcinomatous tumour component coexists with a sarcomatoid histology, is a predictor of poor prognosis in clear cell renal cell carcinoma. Our understanding of sarcomatoid change has been hindered by the lack of molecular examination. Thus, we sought to characterize molecularly the biphasic epithelioid and sarcomatoid components of sarcomatoid clear cell renal cell carcinoma and compare them to non-sarcomatoid clear cell renal cell carcinoma. We examined the transcriptome of the epithelioid and sarcomatoid components of advanced stage sarcomatoid clear cell renal cell carcinoma (n=43) and non-sarcomatoid clear cell renal cell carcinoma (n=37) from independent discovery and validation cohorts using the cDNA microarray and RNA-seq platforms. We analyzed DNA copy number profiles, generated using SNP arrays, from patients with sarcomatoid clear cell renal cell carcinoma (n=10) and advanced non-sarcomatoid clear cell renal cell carcinoma (n=155). The epithelioid and sarcomatoid components of sarcomatoid clear cell renal cell carcinoma had similar gene expression and DNA copy number signatures that were, however, distinct from those of high-grade, high-stage non-sarcomatoid clear cell renal cell carcinoma. Prognostic clear cell renal cell carcinoma gene expression profiles were shared by the biphasic components of sarcomatoid clear cell renal cell carcinoma and the sarcomatoid component showed a partial epithelial-to-mesenchymal transition signature. Our genome-scale microarray-based transcript data were validated in an independent set of sarcomatoid and non-sarcomatoid clear cell renal cell carcinomas using RNA-seq. Sarcomatoid clear cell renal cell carcinoma is molecularly distinct from non-sarcomatoid clear cell renal cell carcinoma, with its genetic programming largely shared by its biphasic morphological components. These data explain why a low percentage of sarcomatoid histology augurs a poor prognosis; suggest the need to modify the pathological grading system and introduce the potential for candidate biomarkers to detect sarcomatoid change preoperatively without specifically sampling the histological sarcomatoid component.

MicroRNAs (miRNAs) simultaneously modulate different oncogenic networks, establishing a dynamic system of gene expression and pathway regulation. In this study, we analyzed global miRNA and messenger RNA (mRNA) expression profiles of 17 cell lines representing different molecular breast cancer subtypes. Spearman's rank correlation test was used to evaluate the correlation between miRNA and mRNA expression. Hierarchical clustering and pathway analysis were also performed. Publicly available gene expression profiles (n = 699) and tumor tissues (n = 80) were analyzed to assess the relevance of key miRNA-regulated pathways in human breast cancer. We identified 39 significantly deregulated miRNAs, and the integration between miRNA and mRNA data revealed the importance of immune-related pathways, particularly the Oncostatin M (OSM) signaling, associated with mesenchymal-like breast cancer cells. OSM levels correlated with genes involved in the inflammatory response, epithelial-to-mesenchymal transition (EMT), and epidermal growth factor (EGF) signaling in human estrogen receptor (ER)-negative/human epidermal growth factor receptor 2 (HER2)-negative breast cancer. Our results suggest that the deregulation of specific miRNAs may cooperatively impair immune and EMT pathways. The identification of the OSM inflammatory pathway as an important mediator of EMT in triple-negative breast cancer (TNBC) may provide a novel potential opportunity to improve therapeutic strategies.

The bone marrow microenvironment is known to provide a survival advantage to residual acute myeloid leukemia cells, possibly contributing to disease recurrence. The mechanisms by which stroma in the microenvironment regulates leukemia survival remain largely unknown. Using reverse-phase protein array technology, we profiled 53 key protein molecules in 11 signaling pathways in 20 primary acute myeloid leukemia samples and two cell lines, aiming to understand stroma-mediated signaling modulation in response to the targeted agents temsirolimus (MTOR), ABT737 (BCL2/BCL-XL), and Nutlin-3a (MDM2), and to identify the effective combination therapy targeting acute myeloid leukemia in the context of the leukemia microenvironment. Stroma reprogrammed signaling networks and modified the sensitivity of acute myeloid leukemia samples to all three targeted inhibitors. Stroma activated AKT at Ser473 in the majority of samples treated with single-agent ABT737 or Nutlin-3a. This survival mechanism was partially abrogated by concomitant treatment with temsirolimus plus ABT737 or Nutlin-3a. Mapping the signaling networks revealed that combinations of two inhibitors increased the number of affected proteins in the targeted pathways and in multiple parallel signaling, translating into facilitated cell death. These results demonstrated that a mechanism-based selection of combined inhibitors can be used to guide clinical drug selection and tailor treatment regimens to eliminate microenvironment-mediated resistance in acute myeloid leukemia.

The prevailing pathomechanistic paradigm for myotonic dystrophy (DM) is that aberrant expression of embryonic/fetal mRNA/protein isoforms accounts for most aspects of the pleiotropic phenotype. To identify aberrant isoforms in skeletal muscle of DM1 and DM2 patients, we performed exon-array profiling and RT-PCR validation on the largest DM sample set to date, including Duchenne, Becker and tibial muscular dystrophy (NMD) patients as disease controls, and non-disease controls. Strikingly, most expression and splicing changes in DM patients were shared with NMD controls. Comparison between DM and NMD identified almost no significant differences. We conclude that DM1 and DM2 are essentially identical for dysregulation of gene expression, and DM expression changes represent a subset of broader spectrum dystrophic changes. We found no evidence for qualitative splicing differences between DM1 and DM2. While some DM-specific splicing differences exist, most of the DM splicing differences were also seen in NMD controls. SSBP3 exon 6 missplicing was observed in all diseased muscle and led to reduced protein. We conclude there is no widespread DM-specific spliceopathy in skeletal muscle and suggest that missplicing in DM (and NMD) may not be the driving mechanism for the muscle pathology, since the same pathways show expression changes unrelated to splicing.

The IDAT file format is used to store BeadArray data from the myriad of genomewide profiling platforms on offer from Illumina Inc. This proprietary format is output directly from the scanner and stores summary intensities for each probe-type on an array in a compact manner. A lack of open source tools to process IDAT files has hampered their uptake by the research community beyond the standard step of using the vendor's software to extract the data they contain in a human readable text format. To fill this void, we have developed the illuminaio package that parses IDAT files from any BeadArray platform, including the decryption of files from Illumina's gene expression arrays. illuminaio provides the first open-source package for this task, and will promote wider uptake of the IDAT format as a standard for sharing Illumina BeadArray data in public databases, in the same way that the CEL file serves as the standard for the Affymetrix platform.

Loss of genome-wide methylation is a common feature of cancer, and the degree of hypomethylation has been correlated with genomic instability. Global methylation of repetitive elements possibly arose as a defense mechanism against parasitic DNA elements, including retrotransposons and viral pathogens. Given the alterations of global methylation in both viral infection and cancer, we examined genome-wide methylation levels in head and neck squamous cell carcinoma (HNSCC), a cancer causally associated with human papilloma virus (HPV). We assayed global hypomethylation levels in 26 HNSCC samples, compared with their matched normal adjacent tissue, using Pyrosequencing-based methylation assays for LINE repeats. In addition, we examined cell lines derived from a variety of solid tumors for LINE and SINE (Alu) repeats. The degree of LINE and Alu hypomethylation varied among different cancer cell lines. There was only moderate correlation between LINE and Alu methylation levels, with the range of variation in methylation levels being greater for the LINE elements. LINE hypomethylation was more pronounced in HPV-negative than in HPV-positive tumors. Moreover, genomic instability, as measured by genome-wide loss-of-heterozygosity (LOH) single nucleotide polymorphism (SNP) analysis, was greater in HNSCC samples with more pronounced LINE hypomethylation. Global hypomethylation was variable in HNSCC. Its correlation with both HPV status and degree of LOH as a surrogate for genomic instability may reflect alternative oncogenic pathways in HPV-positive versus HPV-negative tumors.

Tumor-predominant splice isoforms were identified during comparative in silico sequence analysis of EST clones, suggesting that global aberrant alternative pre-mRNA splicing may be an epigenetic phenomenon in cancer. We used an exon expression array to perform an objective, genome-wide survey of glioma-specific splicing in 24 GBM and 12 nontumor brain samples. Validation studies were performed using RT-PCR on glioma cell lines, patient tumor and nontumor brain samples.

The diversity and heterogeneity within high-grade serous ovarian cancer (HGSC), which is the most lethal gynecologic malignancy, is not well understood. Here, we perform comprehensive multi-platform omics analyses, including integrated analysis, and immune monitoring on primary and metastatic sites from highly clinically annotated HGSC samples based on a laparoscopic triage algorithm from patients who underwent complete gross resection (R0) or received neoadjuvant chemotherapy (NACT) with excellent or poor response. We identify significant distinct molecular abnormalities and cellular changes and immune cell repertoire alterations between the groups, including a higher rate of NF1 copy number loss, and reduced chromothripsis-like patterns, higher levels of strong-binding neoantigens, and a higher number of infiltrated T cells in the R0 versus the NACT groups.

Three dimensional (3D) growths of cancer cells in vitro are more reflective of in situ cancer cell growth than growth in monolayer (2D). The present study is designed to determine changes in protein and phosphoprotein that reflect adaptation of tumor cells to 3D as compared to 2D. Since relative hypoxia is a common feature of most solid tumors, the present study also aims to look at the impact of transition from normoxia to hypoxia in these two growth conditions.

The small molecule anticancer agent NSC-743380 modulates functions of multiple cancer-related pathways and is highly active in a subset of cancer cell lines in the NCI-60 cell line panel. It also has promising in vivo anticancer activity. However, the mechanisms underlying NSC-743380's selective anticancer activity remain uncharacterized. To determine biomarkers that may be used to identify responders to this novel anticancer agent, we performed correlation analysis on NSC-743380's anticancer activity and the gene expression levels in NCI-60 cell lines and characterized the functions of the top associated genes in NSC-743380-mediated anticancer activity. We found sulfotransferase SULT1A1 is causally associated with NSC-743380's anticancer activity. SULT1A1 was expressed in NSC-743380-sensitive cell lines but was undetectable in resistant cancer cells. Ectopic expression of SULT1A1 in NSC743380 resistant cancer cells dramatically sensitized the resistant cells to NSC-743380. Knockdown of the SULT1A1 in the NSC-743380 sensitive cancer cell line rendered it resistance to NSC-743380. The SULT1A1 protein levels in cell lysates from 18 leukemia cell lines reliably predicted the susceptibility of the cell lines to NSC-743380. Thus, expression of SULT1A1 in cancer cells is required for NSC-743380's anticancer activity and can be used as a biomarker for identification of NSC-743380 responders.

We analyzed molecular data on 2,579 tumors from The Cancer Genome Atlas (TCGA) of four gynecological types plus breast. Our aims were to identify shared and unique molecular features, clinically significant subtypes, and potential therapeutic targets. We found 61 somatic copy-number alterations (SCNAs) and 46 significantly mutated genes (SMGs). Eleven SCNAs and 11 SMGs had not been identified in previous TCGA studies of the individual tumor types. We found functionally significant estrogen receptor-regulated long non-coding RNAs (lncRNAs) and gene/lncRNA interaction networks. Pathway analysis identified subtypes with high leukocyte infiltration, raising potential implications for immunotherapy. Using 16 key molecular features, we identified five prognostic subtypes and developed a decision tree that classified patients into the subtypes based on just six features that are assessable in clinical laboratories.