Transcriptomic analyses of ovarian clear-cell carcinoma with concurrent endometriosis.

Introduction: Endometriosis, a benign inflammatory disease whereby endometrial-like tissue grows outside the uterus, is a risk factor for endometriosis-associated ovarian cancers. In particular, ovarian endometriomas, cystic lesions of deeply invasive endometriosis, are considered the precursor lesion for ovarian clear-cell carcinoma (OCCC). Methods: To explore this transcriptomic landscape, OCCC from women with pathology-proven concurrent endometriosis (n = 4) were compared to benign endometriomas (n = 4) by bulk RNA and small-RNA sequencing. Results: Analysis of protein-coding genes identified 2449 upregulated and 3131 downregulated protein-coding genes (DESeq2, P< 0.05, log2 fold-change > |1|) in OCCC with concurrent endometriosis compared to endometriomas. Gene set enrichment analysis showed upregulation of pathways involved in cell cycle regulation and DNA replication and downregulation of pathways involved in cytokine receptor signaling and matrisome. Comparison of pathway activation scores between the clinical samples and publicly-available datasets for OCCC cell lines revealed significant molecular similarities between OCCC with concurrent endometriosis and OVTOKO, OVISE, RMG1, OVMANA, TOV21G, IGROV1, and JHOC5 cell lines. Analysis of miRNAs revealed 64 upregulated and 61 downregulated mature miRNA molecules (DESeq2, P< 0.05, log2 fold-change > |1|). MiR-10a-5p represented over 21% of the miRNA molecules in OCCC with endometriosis and was significantly upregulated (NGS: log2fold change = 4.37, P = 2.43e-18; QPCR: 8.1-fold change, P< 0.05). Correlation between miR-10a expression level in OCCC cell lines and IC50 (50% inhibitory concentration) of carboplatin in vitro revealed a positive correlation (R2 = 0.93). MiR-10a overexpression in vitro resulted in a significant decrease in proliferation (n = 6; P< 0.05) compared to transfection with a non-targeting control miRNA. Similarly, the cell-cycle analysis revealed a significant shift in cells from S and G2 to G1 (n = 6; P< 0.0001). Bioinformatic analysis predicted that miR-10a-5p target genes that were downregulated in OCCC with endometriosis were involved in receptor signaling pathways, proliferation, and cell cycle progression. MiR-10a overexpression in vitro was correlated with decreased expression of predicted miR-10a target genes critical for proliferation, cell-cycle regulation, and cell survival including [SERPINE1 (3-fold downregulated; P< 0.05), CDK6 (2.4-fold downregulated; P< 0.05), and RAP2A (2-3-fold downregulated; P< 0.05)]. Discussion: These studies in OCCC suggest that miR-10a-5p is an impactful, potentially oncogenic molecule, which warrants further studies.

ZEB2 regulates endocrine therapy sensitivity and metastasis in luminal a breast cancer cells through a non-canonical mechanism.

PURPOSE: The transcription factors ZEB1 and ZEB2 mediate epithelial-to-mesenchymal transition (EMT) and metastatic progression in numerous malignancies including breast cancer. ZEB1 and ZEB2 drive EMT through transcriptional repression of cell-cell junction proteins and members of the tumor suppressive miR200 family. However, in estrogen receptor positive (ER +) breast cancer, the role of ZEB2 as an independent driver of metastasis has not been fully investigated. METHODS: In the current study, we induced exogenous expression of ZEB2 in ER + MCF-7 and ZR-75-1 breast cancer cell lines and examined EMT gene expression and metastasis using dose-response qRT-PCR, transwell migration assays, proliferation assays with immunofluorescence of Ki-67 staining. We used RNA sequencing to identify pathways and genes affected by ZEB2 overexpression. Finally, we treated ZEB2-overexpressing cells with 17ß-estradiol (E2) or ICI 182,780 to evaluate how ZEB2 affects estrogen response. RESULTS: Contrary to expectation, we found that ZEB2 did not increase canonical epithelial nor decrease mesenchymal gene expressions. Furthermore, ZEB2 overexpression did not promote a mesenchymal cell morphology. However, ZEB1 and ZEB2 protein expression induced significant migration of MCF-7 and ZR-75-1 breast cancer cells in vitro and MCF-7 xenograft metastasis in vivo. Transcriptomic (RNA sequencing) pathway analysis revealed alterations in estrogen signaling regulators and pathways, suggesting a role for ZEB2 in endocrine sensitivity in luminal A breast cancer. Expression of ZEB2 was negatively correlated with estrogen receptor complex genes in luminal A patient tumors. Furthermore, treatment with 17ß-estradiol (E2) or the estrogen receptor antagonist ICI 182,780 had no effect on growth of ZEB2-overexpressing cells. CONCLUSION: ZEB2 is a multi-functional regulator of drug sensitivity, cell migration, and metastasis in ER + breast cancer and functions through non-canonical mechanisms.

NEK5 activity regulates the mesenchymal and migratory phenotype in breast cancer cells.

PURPOSE: Breast cancer remains a prominent global disease affecting women worldwide despite the emergence of novel therapeutic regimens. Metastasis is responsible for most cancer-related deaths, and acquisition of a mesenchymal and migratory cancer cell phenotypes contributes to this devastating disease. The utilization of kinase targets in drug discovery have revolutionized the field of cancer research but despite impressive advancements in kinase-targeting drugs, a large portion of the human kinome remains understudied in cancer. NEK5, a member of the Never-in-mitosis kinase family, is an example of such an understudied kinase. Here, we characterized the function of NEK5 in breast cancer. METHODS: Stably overexpressing NEK5 cell lines (MCF7) and shRNA knockdown cell lines (MDA-MB-231, TU-BcX-4IC) were utilized. Cell morphology changes were evaluated using immunofluorescence and quantification of cytoskeletal components. Cell proliferation was assessed by Ki-67 staining and transwell migration assays tested cell migration capabilities. In vivo experiments with murine models were necessary to demonstrate NEK5 function in breast cancer tumor growth and metastasis. RESULTS: NEK5 activation altered breast cancer cell morphology and promoted cell migration independent of effects on cell proliferation. NEK5 overexpression or knockdown does not alter tumor growth kinetics but promotes or suppresses metastatic potential in a cell type-specific manner, respectively. CONCLUSION: While NEK5 activity modulated cytoskeletal changes and cell motility, NEK5 activity affected cell seeding capabilities but not metastatic colonization or proliferation in vivo. Here we characterized NEK5 function in breast cancer systems and we implicate NEK5 in regulating specific steps of metastatic progression.

Constitutive activation of MEK5 promotes a mesenchymal and migratory cell phenotype in triple negative breast cancer.

Triple negative breast cancer (TNBC) is an aggressive subtype of breast cancer with limited targeted therapeutic options. A defining feature of TNBC is the propensity to metastasize and acquire resistance to cytotoxic agents. Mitogen activated protein kinase (MAPK) and extracellular regulated kinase (ERK) signaling pathways have integral roles in cancer development and progression. While MEK5/ERK5 signaling drives mesenchymal and migratory cell phenotypes in breast cancer, the specific mechanisms underlying these actions remain under-characterized. To elucidate the mechanisms through which MEK5 regulates the mesenchymal and migratory phenotype, we generated stably transfected constitutively active MEK5 (MEK5-ca) TNBC cells. Downstream signaling pathways and candidate targets of MEK5-ca cells were based on RNA sequencing and confirmed using qPCR and Western blot analyses. MEK5 activation drove a mesenchymal cell phenotype independent of cell proliferation effects. Transwell migration assays demonstrated MEK5 activation significantly increased breast cancer cell migration. In this study, we provide supporting evidence that MEK5 functions through FRA-1 to regulate the mesenchymal and migratory phenotype in TNBC.

Dual inhibition of MEK1/2 and MEK5 suppresses the EMT/migration axis in triple-negative breast cancer through FRA-1 regulation.

Triple-negative breast cancer (TNBC) presents a clinical challenge due to the aggressive nature of the disease and a lack of targeted therapies. Constitutive activation of the mitogen-activated protein kinase (MAPK)/extracellular signal-regulated kinase (ERK) pathway has been linked to chemoresistance and metastatic progression through distinct mechanisms, including activation of epithelial-to-mesenchymal transition (EMT) when cells adopt a motile and invasive phenotype through loss of epithelial markers (CDH1), and acquisition of mesenchymal markers (VIM, CDH2). Although MAPK/ERK1/2 kinase inhibitors (MEKi) are useful antitumor agents in a clinical setting, including the Food and Drug Administration (FDA)-approved MEK1,2 dual inhibitors cobimetinib and trametinib, there are limitations to their clinical utility, primarily adaptation of the BRAF pathway and ocular toxicities. The MEK5 (HGNC: MAP2K5) pathway has important roles in metastatic progression of various cancer types, including those of the prostate, colon, bone and breast, and elevated levels of ERK5 expression in breast carcinomas are linked to a worse prognoses in TNBC patients. The purpose of this study is to explore MEK5 regulation of the EMT axis and to evaluate a novel pan-MEK inhibitor on clinically aggressive TNBC cells. Our results show a distinction between the MEK1/2 and MEK5 cascades in maintenance of the mesenchymal phenotype, suggesting that the MEK5 pathway may be necessary and sufficient in EMT regulation while MEK1/2 signaling further sustains the mesenchymal state of TNBC cells. Furthermore, additive effects on MET induction are evident through the inhibition of both MEK1/2 and MEK5. Taken together, these data demonstrate the need for a better understanding of the individual roles of MEK1/2 and MEK5 signaling in breast cancer and provide a rationale for the combined targeting of these pathways to circumvent compensatory signaling and subsequent therapeutic resistance.

Regulation of cellular sterol homeostasis by the oxygen responsive noncoding RNA lincNORS.

We hereby provide the initial portrait of lincNORS, a spliced lincRNA generated by the MIR193BHG locus, entirely distinct from the previously described miR-193b-365a tandem. While inducible by low O2 in a variety of cells and associated with hypoxia in vivo, our studies show that lincNORS is subject to multiple regulatory inputs, including estrogen signals. Biochemically, this lincRNA fine-tunes cellular sterol/steroid biosynthesis by repressing the expression of multiple pathway components. Mechanistically, the function of lincNORS requires the presence of RALY, an RNA-binding protein recently found to be implicated in cholesterol homeostasis. We also noticed the proximity between this locus and naturally occurring genetic variations highly significant for sterol/steroid-related phenotypes, in particular the age of sexual maturation. An integrative analysis of these variants provided a more formal link between these phenotypes and lincNORS, further strengthening the case for its biological relevance.

A novel screening approach comparing kinase activity of small molecule inhibitors with similar molecular structures and distinct biologic effects in triple-negative breast cancer to identify targetable signaling pathways.

Breast cancer affects women globally; the majority of breast cancer-related mortalities are due to metastasis. Acquisition of a mesenchymal phenotype has been implicated in the progression of breast cancer cells to an invasive, metastatic state. Triple-negative breast cancer (TNBC) subtypes have high rates of metastases, recurrence, and have poorer prognoses compared to other breast cancer types, partially due to lack of commonly targeted receptors. Kinases have diverse and pivotal functions in metastasis in TNBC, and discovery of new kinase targets for TNBC is warranted. We previously used a screening approach to identify intermediate-synthesis nonpotent, nonselective small-molecule inhibitors from the Published Kinase Inhibitor Set that reversed the mesenchymal phenotype in TNBC cells. Two of these inhibitors (GSK346294A and GSK448459A) are structurally similar, but have unique kinase activity profiles and exhibited differential biologic effects on TNBC cells, specifically on epithelial-to-mesenchymal transition (EMT). Here, we further interrogate these effects and compare activity of these inhibitors on transwell migration, gene (qRT-PCR) and protein (western blot) expressions, and cancer stem cell-like behavior. We incorporated translational patient-derived xenograft models in these studies, and we focused on the lead inhibitor hit, GSK346294A, to demonstrate the utility of our comparative analysis as a screening modality to identify novel kinase targets and signaling pathways to pursue in TNBC. This study introduces a new method for discovering novel kinase targets that reverse the EMT phenotype; this screening approach can be applied to all cancer types and is not limited to breast cancer.

Pten and Dicer1 loss in the mouse uterus causes poorly differentiated endometrial adenocarcinoma.

Endometrial cancer remains the most common gynecological malignancy in the United States. While the loss of the tumor suppressor, PTEN (phosphatase and tensin homolog), is well studied in endometrial cancer, recent studies suggest that DICER1, the endoribonuclease responsible for miRNA genesis, also plays a significant role in endometrial adenocarcinoma. Conditional uterine deletion of Dicer1 and Pten in mice resulted in poorly differentiated endometrial adenocarcinomas, which expressed Napsin A and HNF1B (hepatocyte nuclear factor 1 homeobox B), markers of clear-cell adenocarcinoma. Adenocarcinomas were hormone-independent. Treatment with progesterone did not mitigate poorly differentiated adenocarcinoma, nor did it affect adnexal metastasis. Transcriptomic analyses of DICER1 deleted uteri or Ishikawa cells revealed unique transcriptomic profiles and global miRNA downregulation. Computational integration of miRNA with mRNA targets revealed deregulated let-7 and miR-16 target genes, similar to published human DICER1-mutant endometrial cancers from TCGA (The Cancer Genome Atlas). Similar to human endometrial cancers, tumors exhibited dysregulation of ephrin-receptor signaling and transforming growth factor-beta signaling pathways. LIM kinase 2 (LIMK2), an essential molecule in p21 signal transduction, was significantly upregulated and represents a novel mechanism for hormone-independent pathogenesis of endometrial adenocarcinoma. This preclinical mouse model represents the first genetically engineered mouse model of poorly differentiated endometrial adenocarcinoma.

Changes in mRNA/protein expression and signaling pathways in in vivo passaged mouse ovarian cancer cells.

The cure rate for late stage epithelial ovarian cancer (EOC) has not significantly improved over several decades. New and more effective targets and treatment modalities are urgently needed. RNA-seq analyses of a syngeneic EOC cell pair, representing more and less aggressive tumor cells in vivo were conducted. Bioinformatics analyses of the RNA-seq data and biological signaling and function studies have identified new targets, such as ZIP4 in EOC. Many up-regulated tumor promoting signaling pathways have been identified which are mainly grouped into three cellular activities: 1) cell proliferation and apoptosis resistance; 2) cell skeleton and adhesion changes; and 3) carbohydrate metabolic reprograming. Unexpectedly, lipid metabolism has been the major down-regulated signaling pathway in the more aggressive EOC cells. In addition, we found that hypoxic responsive genes were at the center stage of regulation and detected functional changes were related to cancer stem cell-like activities. Moreover, our genetic, cellular, biochemical, and lipidomic analyses indicated that cells grown in 2D vs. 3D, or attached vs. suspended had dramatic changes. The important clinical implications of peritoneal cavity floating tumor cells are supported by the data proved in this work. Overall, the RNA-seq data provide a landscape of gene expression alterations during tumor progression.

Epigenetic Targeting of Adipocytes Inhibits High-Grade Serous Ovarian Cancer Cell Migration and Invasion.

Ovarian cancer (OC) cells frequently metastasize to the omentum, and adipocytes play a significant role in ovarian tumor progression. Therapeutic interventions targeting aberrant DNA methylation in ovarian tumors have shown promise in the clinic, but the effects of epigenetic therapy on the tumor microenvironment are understudied. Here, we examined the effect of adipocytes on OC cell behavior in culture and impact of targeting DNA methylation in adipocytes on OC metastasis. The presence of adipocytes increased OC cell migration and invasion, and proximal and direct coculture of adipocytes increased OC proliferation alone or after treatment with carboplatin. Treatment of adipocytes with hypomethylating agent guadecitabine decreased migration and invasion of OC cells toward adipocytes. Subcellular protein fractionation of adipocytes treated with guadecitabine revealed decreased DNA methyltransferase 1 (DNMT1) levels even in the presence of DNA synthesis inhibitor, aphidicolin. Methyl-Capture- and RNA-sequencing analysis of guadecitabine-treated adipocytes revealed derepression of tumor-suppressor genes and epithelial-mesenchymal transition inhibitors. SUSD2, a secreted tumor suppressor downregulated by promoter CpG island methylation in adipocytes, was upregulated after guadecitabine treatment, and recombinant SUSD2 decreased OC cell migration and invasion. Integrated analysis of the methylomic and transcriptomic data identified pathways associated with inhibition of matrix metalloproteases and fatty acid α-oxidation, suggesting a possible mechanism of how epigenetic therapy of adipocytes decreases metastasis. In conclusion, the effect of DNMT inhibitor on fully differentiated adipocytes suggests that hypomethylating agents may affect the tumor microenvironment to decrease cancer cell metastasis.Implications: Epigenetic targeting of tumor microenvironment can affect metastatic behavior of ovarian cancer cells. Mol Cancer Res; 16(8); 1226-40. ©2018 AACR.

Complete Transcriptome RNA-Seq.

RNA-Seq is the leading technology for analyzing gene expression on a global scale across a broad spectrum of sample types. However, due to chemical modifications by fixation or degradation due to collection methods, samples often contain an abundance of RNA that is no longer intact, and the capability of current RNA-Seq protocols to accurately quantify such samples is often limited. We have developed an RNA-Seq protocol to address these key issues as well as quantify gene expression from the whole transcriptome. Furthermore, for compatibility with improved sequencing platforms, we use restructured adapter sequences to generate libraries for Illumina HiSeq, MiSeq, and NextSeq platforms. Our protocol utilizes duplex-specific nuclease (DSN) to remove abundant ribosomal RNA sequences while retaining other types of RNA for superior transcriptome profiling from low quantity input. We employ the Illumina sequencing platform, but this method is described in sufficient detail to adapt to other platforms.

MicroRNA-335-5p and -3p synergize to inhibit estrogen receptor alpha expression and promote tamoxifen resistance.

microRNAs (miRNAs) are small noncoding RNA molecules involved in the regulation of gene expression and play critical roles in human malignancies. Next-generation sequencing analysis of the MCF-7 breast cancer cell line overexpressing miR-335-5p and miR-335-3p demonstrated that the miRNA duplex repressed genes involved in the ERα signaling pathway, and enhanced resistance of MCF-7 cells to the growth inhibitory effects of tamoxifen. These data suggest that despite its conventional role in tumor suppression, the miR-335 transcript can also play an oncogenic role in promoting agonistic estrogen signaling in a cancerous setting.

Discrete Responses to Limitation for Iron and Manganese in Agrobacterium tumefaciens: Influence on Attachment and Biofilm Formation.

UNLABELLED: Transition metals such as iron and manganese are crucial trace nutrients for the growth of most bacteria, functioning as catalytic cofactors for many essential enzymes. Dedicated uptake and regulatory systems have evolved to ensure their acquisition for growth, while preventing toxicity. Transcriptomic analysis of the iron- and manganese-responsive regulons of Agrobacterium tumefaciens revealed that there are discrete regulatory networks that respond to changes in iron and manganese levels. Complementing earlier studies, the iron-responsive gene network is quite large and includes many aspects of iron-dependent metabolism and the iron-sparing response. In contrast, the manganese-responsive network is restricted to a limited number of genes, many of which can be linked to transport and utilization of the transition metal. Several of the target genes predicted to drive manganese uptake are required for growth under manganese-limited conditions, and an A. tumefaciens mutant with a manganese transport deficiency is attenuated for plant virulence. Iron and manganese limitation independently inhibit biofilm formation by A. tumefaciens, and several candidate genes that could impact biofilm formation were identified in each regulon. The biofilm-inhibitory effects of iron and manganese do not rely on recognized metal-responsive transcriptional regulators, suggesting alternate mechanisms influencing biofilm formation. However, under low-manganese conditions the dcpA operon is upregulated, encoding a system that controls levels of the cyclic di-GMP second messenger. Mutation of this regulatory pathway dampens the effect of manganese limitation. IMPORTANCE: Responses to changes in transition metal levels, such as those of manganese and iron, are important for normal metabolism and growth in bacteria. Our study used global gene expression profiling to understand the response of the plant pathogen Agrobacterium tumefaciens to changes of transition metal availability. Among the properties that are affected by both iron and manganese levels are those required for normal surface attachment and biofilm formation, but the requirement for each of these transition metals is mechanistically independent from the other.

Dual regulation by microRNA-200b-3p and microRNA-200b-5p in the inhibition of epithelial-to-mesenchymal transition in triple-negative breast cancer.

Epithelial to mesenchymal transition (EMT) involves loss of an epithelial phenotype and activation of a mesenchymal one. Enhanced expression of genes associated with a mesenchymal transition includes ZEB1/2, TWIST, and FOXC1. miRNAs are known regulators of gene expression and altered miRNA expression is known to enhance EMT in breast cancer. Here we demonstrate that the tumor suppressive miRNA family, miR-200, is not expressed in triple negative breast cancer (TNBC) cell lines and that miR-200b-3p over-expression represses EMT, which is evident through decreased migration and increased CDH1 expression. Despite the loss of migratory capacity following re-expression of miR-200b-3p, no subsequent loss of the conventional miR-200 family targets and EMT markers ZEB1/2 was observed. Next generation RNA-sequencing analysis showed that enhanced expression of pri-miR-200b lead to ectopic expression of both miR-200b-3p and miR-200b-5p with multiple isomiRs expressed for each of these miRNAs. Furthermore, miR-200b-5p was expressed in the receptor positive, epithelial breast cancer cell lines but not in the TNBC (mesenchymal) cell lines. In addition, a compensatory mechanism for miR-200b-3p/200b-5p targeting, where both miRNAs target the RHOGDI pathway leading to non-canonical repression of EMT, was demonstrated. Collectively, these data are the first to demonstrate dual targeting by miR-200b-3p and miR-200b-5p and a previously undescribed role for microRNA processing and strand expression in EMT and TNBC, the most aggressive breast cancer subtype.

Effects of the endocrine-disrupting chemical DDT on self-renewal and differentiation of human mesenchymal stem cells.

BACKGROUND: Although the global use of the endocrine-disrupting chemical DDT has decreased, its persistence in the environment has resulted in continued human exposure. Accumulating evidence suggests that DDT exposure has long-term adverse effects on development, yet the impact on growth and differentiation of adult stem cells remains unclear. OBJECTIVES: Human mesenchymal stem cells (MSCs) exposed to DDT were used to evaluate the impact on stem cell biology. METHODS: We assessed DDT-treated MSCs for self-renewal, proliferation, and differentiation potential. Whole genome RNA sequencing was performed to assess gene expression in DDT-treated MSCs. RESULTS: MSCs exposed to DDT formed fewer colonies, suggesting a reduction in self-renewal potential. DDT enhanced both adipogenic and osteogenic differentiation, which was confirmed by increased mRNA expression of glucose transporter type 4 (GLUT4), lipoprotein lipase (LpL), peroxisome proliferator-activated receptor gamma (PPARγ), leptin, osteonectin, core binding factor 1 (CBFA1), and FBJ murine osteosarcoma viral oncogene homolog (c-Fos). Expression of factors in DDT-treated cells was similar to that in estrogen-treated MSCs, suggesting that DDT may function via the estrogen receptor (ER)-mediated pathway. The coadministration of ICI 182,780 blocked the effects of DDT. RNA sequencing revealed 121 genes and noncoding RNAs to be differentially expressed in DDT-treated MSCs compared with controls cells. CONCLUSION: Human MSCs provide a powerful biological system to investigate and identify the molecular mechanisms underlying the effects of environmental agents on stem cells and human health. MSCs exposed to DDT demonstrated profound alterations in self-renewal, proliferation, differentiation, and gene expression, which may partially explain the homeostatic imbalance and increased cancer incidence among those exposed to long-term EDCs.

A new method for stranded whole transcriptome RNA-seq.

This report describes an improved protocol to generate stranded, barcoded RNA-seq libraries to capture the whole transcriptome. By optimizing the use of duplex specific nuclease (DSN) to remove ribosomal RNA reads from stranded barcoded libraries, we demonstrate improved efficiency of multiplexed next generation sequencing (NGS). This approach detects expression profiles of all RNA types, including miRNA (microRNA), piRNA (Piwi-interacting RNA), snoRNA (small nucleolar RNA), lincRNA (long non-coding RNA), mtRNA (mitochondrial RNA) and mRNA (messenger RNA) without the use of gel electrophoresis. The improved protocol generates high quality data that can be used to identify differential expression in known and novel coding and non-coding transcripts, splice variants, mitochondrial genes and SNPs (single nucleotide polymorphisms).