Nuclear export of circular RNA.

Circular RNAs (circRNAs), which are increasingly being implicated in a variety of functions in normal and cancerous cells1-5, are formed by back-splicing of precursor mRNAs in the nucleus6-10. circRNAs are predominantly localized in the cytoplasm, indicating that they must be exported from the nucleus. Here we identify a pathway that is specific for the nuclear export of circular RNA. This pathway requires Ran-GTP, exportin-2 and IGF2BP1. Enhancing the nuclear Ran-GTP gradient by depletion or chemical inhibition of the major protein exporter CRM1 selectively increases the nuclear export of circRNAs, while reducing the nuclear Ran-GTP gradient selectively blocks circRNA export. Depletion or knockout of exportin-2 specifically inhibits nuclear export of circRNA. Analysis of nuclear circRNA-binding proteins reveals that interaction between IGF2BP1 and circRNA is enhanced by Ran-GTP. The formation of circRNA export complexes in the nucleus is promoted by Ran-GTP through its interactions with exportin-2, circRNA and IGF2BP1. Our findings demonstrate that adaptors such as IGF2BP1 that bind directly to circular RNAs recruit Ran-GTP and exportin-2 to export circRNAs in a mechanism that is analogous to protein export, rather than mRNA export.

ESRP1 controls biogenesis and function of a large abundant multiexon circRNA.

While the majority of circRNAs are formed from infrequent back-splicing of exons from protein coding genes, some can be produced at quite high level and in a regulated manner. We describe the regulation, biogenesis and function of circDOCK1(2-27), a large, abundant circular RNA that is highly regulated during epithelial-mesenchymal transition (EMT) and whose formation depends on the epithelial splicing regulator ESRP1. CircDOCK1(2-27) synthesis in epithelial cells represses cell motility both by diverting transcripts from DOCK1 mRNA production to circRNA formation and by direct inhibition of migration by the circRNA. HITS-CLIP analysis and CRISPR-mediated deletions indicate ESRP1 controls circDOCK1(2-27) biosynthesis by binding a GGU-containing repeat region in intron 1 and detaining its splicing until Pol II completes its 157 kb journey to exon 27. Proximity-dependent biotinylation (BioID) assay suggests ESRP1 may modify the RNP landscape of intron 1 in a way that disfavours communication of exon 1 with exon 2, rather than physically bridging exon 2 to exon 27. The X-ray crystal structure of RNA-bound ESRP1 qRRM2 domain reveals it binds to GGU motifs, with the guanines embedded in clamp-like aromatic pockets in the protein.

Ceramide-induced integrated stress response overcomes Bcl-2 inhibitor resistance in acute myeloid leukemia.

Inducing cell death by the sphingolipid ceramide is a potential anticancer strategy, but the underlying mechanisms remain poorly defined. In this study, triggering an accumulation of ceramide in acute myeloid leukemia (AML) cells by inhibition of sphingosine kinase induced an apoptotic integrated stress response (ISR) through protein kinase R-mediated activation of the master transcription factor ATF4. This effect led to transcription of the BH3-only protein Noxa and degradation of the prosurvival Mcl-1 protein on which AML cells are highly dependent for survival. Targeting this novel ISR pathway, in combination with the Bcl-2 inhibitor venetoclax, synergistically killed primary AML blasts, including those with venetoclax-resistant mutations, as well as immunophenotypic leukemic stem cells, and reduced leukemic engraftment in patient-derived AML xenografts. Collectively, these findings provide mechanistic insight into the anticancer effects of ceramide and preclinical evidence for new approaches to augment Bcl-2 inhibition in the therapy of AML and other cancers with high Mcl-1 dependency.

The landscape of alternative polyadenylation during EMT and its regulation by the RNA-binding protein Quaking.

Epithelial-mesenchymal transition (EMT) plays important roles in tumour progression and is orchestrated by dynamic changes in gene expression. While it is well established that post-transcriptional regulation plays a significant role in EMT, the extent of alternative polyadenylation (APA) during EMT has not yet been explored. Using 3' end anchored RNA sequencing, we mapped the alternative polyadenylation (APA) landscape following Transforming Growth Factor (TGF)-ß-mediated induction of EMT in human mammary epithelial cells and found APA generally causes 3'UTR lengthening during this cell state transition. Investigation of potential mediators of APA indicated the RNA-binding protein Quaking (QKI), a splicing factor induced during EMT, regulates a subset of events including the length of its own transcript. Analysis of QKI crosslinked immunoprecipitation (CLIP)-sequencing data identified the binding of QKI within 3' untranslated regions (UTRs) was enriched near cleavage and polyadenylation sites. Following QKI knockdown, APA of many transcripts is altered to produce predominantly shorter 3'UTRs associated with reduced gene expression. These findings reveal the changes in APA that occur during EMT and identify a potential role for QKI in this process.

miR-200/375 control epithelial plasticity-associated alternative splicing by repressing the RNA-binding protein Quaking.

Members of the miR-200 family are critical gatekeepers of the epithelial state, restraining expression of pro-mesenchymal genes that drive epithelial-mesenchymal transition (EMT) and contribute to metastatic cancer progression. Here, we show that miR-200c and another epithelial-enriched miRNA, miR-375, exert widespread control of alternative splicing in cancer cells by suppressing the RNA-binding protein Quaking (QKI). During EMT, QKI-5 directly binds to and regulates hundreds of alternative splicing targets and exerts pleiotropic effects, such as increasing cell migration and invasion and restraining tumour growth, without appreciably affecting mRNA levels. QKI-5 is both necessary and sufficient to direct EMT-associated alternative splicing changes, and this splicing signature is broadly conserved across many epithelial-derived cancer types. Importantly, several actin cytoskeleton-associated genes are directly targeted by both QKI and miR-200c, revealing coordinated control of alternative splicing and mRNA abundance during EMT These findings demonstrate the existence of a miR-200/miR-375/QKI axis that impacts cancer-associated epithelial cell plasticity through widespread control of alternative splicing.

Extensive transcriptional responses are co-ordinated by microRNAs as revealed by Exon-Intron Split Analysis (EISA).

Epithelial-mesenchymal transition (EMT) has been a subject of intense scrutiny as it facilitates metastasis and alters drug sensitivity. Although EMT-regulatory roles for numerous miRNAs and transcription factors are known, their functions can be difficult to disentangle, in part due to the difficulty in identifying direct miRNA targets from complex datasets and in deciding how to incorporate 'indirect' miRNA effects that may, or may not, represent biologically relevant information. To better understand how miRNAs exert effects throughout the transcriptome during EMT, we employed Exon-Intron Split Analysis (EISA), a bioinformatic technique that separates transcriptional and post-transcriptional effects through the separate analysis of RNA-Seq reads mapping to exons and introns. We find that in response to the manipulation of miRNAs, a major effect on gene expression is transcriptional. We also find extensive co-ordination of transcriptional and post-transcriptional regulatory mechanisms during both EMT and mesenchymal to epithelial transition (MET) in response to TGF-ß or miR-200c respectively. The prominent transcriptional influence of miRNAs was also observed in other datasets where miRNA levels were perturbed. This work cautions against a narrow approach that is limited to the analysis of direct targets, and demonstrates the utility of EISA to examine complex regulatory networks involving both transcriptional and post-transcriptional mechanisms.

Genome-wide identification of miR-200 targets reveals a regulatory network controlling cell invasion.

The microRNAs of the miR-200 family maintain the central characteristics of epithelia and inhibit tumor cell motility and invasiveness. Using the Ago-HITS-CLIP technology for transcriptome-wide identification of direct microRNA targets in living cells, along with extensive validation to verify the reliability of the approach, we have identified hundreds of miR-200a and miR-200b targets, providing insights into general features of miRNA target site selection. Gene ontology analysis revealed a predominant effect of miR-200 targets in widespread coordinate control of actin cytoskeleton dynamics. Functional characterization of the miR-200 targets indicates that they constitute subnetworks that underlie the ability of cancer cells to migrate and invade, including coordinate effects on Rho-ROCK signaling, invadopodia formation, MMP activity, and focal adhesions. Thus, the miR-200 family maintains the central characteristics of the epithelial phenotype by acting on numerous targets at multiple levels, encompassing both cytoskeletal effectors that control actin filament organization and dynamics, and upstream signals that locally regulate the cytoskeleton to maintain cell morphology and prevent cell migration.

Specificity protein 1 (Sp1) maintains basal epithelial expression of the miR-200 family: implications for epithelial-mesenchymal transition.

Epithelial-mesenchymal transition (EMT) is required for the specification of tissues during embryonic development and is recapitulated during the metastatic progression of tumors. The miR-200 family plays a critical role in enforcing the epithelial state with their expression lost in cells undergoing EMT. EMT can be mediated by activation of the ZEB1 and ZEB2 (ZEB) transcription factors, which repress miR-200 expression via a self-reinforcing double negative feedback loop to promote the mesenchymal state. However, it remains unclear what factors drive and maintain epithelial-specific expression of miR-200 in the absence of EMT-inducing factors. Here, we show that the transcription factor Specificity Protein 1 (Sp1) binds to the miR-200b∼200a∼429 proximal promoter and activates miR-200 expression in epithelial cells. In mesenchymal cells, Sp1 expression is maintained, but its ability to activate the miR-200 promoter is perturbed by ZEB-mediated repression. Reduction of Sp1 expression caused changes in EMT-associated markers in epithelial cells. Furthermore, we observed co-expression of Sp1 and miR-200 during mouse embryonic development wherein miR-200 expression was only lost in regions with high ZEB expression. Together, these findings indicate that miR-200 family members require Sp1 to drive basal expression and to maintain an epithelial state.

Epigenetic modulation of the miR-200 family is associated with transition to a breast cancer stem-cell-like state.

The miR-200 family is a key regulator of the epithelial-mesenchymal transition, however, its role in controlling the transition between cancer stem-cell-like and non-stem-cell-like phenotypes is not well understood. We utilized immortalized human mammary epithelial (HMLE) cells to investigate the regulation of the miR-200 family during their conversion to a stem-like phenotype. HMLE cells were found to be capable of spontaneous conversion from a non-stem to a stem-like phenotype and this conversion was accompanied by the loss of miR-200 expression. Stem-like cell fractions isolated from metastatic breast cancers also displayed loss of miR-200 indicating similar molecular changes may occur during breast cancer progression. The phenotypic change observed in HMLE cells was directly controlled by miR-200 because restoration of its expression decreased stem-like properties while promoting a transition to an epithelial phenotype. Investigation of the mechanisms controlling miR-200 expression revealed both DNA methylation and histone modifications were significantly altered in the stem-like and non-stem phenotypes. In particular, in the stem-like phenotype, the miR-200b-200a-429 cluster was silenced primarily through polycomb group-mediated histone modifications whereas the miR-200c-141 cluster was repressed by DNA methylation. These results indicate that the miR-200 family plays a crucial role in the transition between stem-like and non-stem phenotypes and that distinct epigenetic-based mechanisms regulate each miR-200 gene in this process. Therapy targeted against miR-200 family members and epigenetic modifications might therefore be applicable to breast cancer.

The inducible tissue-specific expression of the human IL-3/GM-CSF locus is controlled by a complex array of developmentally regulated enhancers.

The closely linked human IL-3 and GM-CSF genes are tightly regulated and are expressed in activated T cells and mast cells. In this study, we used transgenic mice to study the developmental regulation of this locus and to identify DNA elements required for its correct activity in vivo. Because these two genes are separated by a CTCF-dependent insulator, and the GM-CSF gene is regulated primarily by its own upstream enhancer, the main objective in this study was to identify regions of the locus required for correct IL-3 gene expression. We initially found that the previously identified proximal upstream IL-3 enhancers were insufficient to account for the in vivo activity of the IL-3 gene. However, an extended analysis of DNase I-hypersensitive sites (DHSs) spanning the entire upstream IL-3 intergenic region revealed the existence of a complex cluster of both constitutive and inducible DHSs spanning the -34- to -40-kb region. The tissue specificity of these DHSs mirrored the activity of the IL-3 gene, and included a highly inducible cyclosporin A-sensitive enhancer at -37 kb that increased IL-3 promoter activity 40-fold. Significantly, inclusion of this region enabled correct in vivo regulation of IL-3 gene expression in T cells, mast cells, and myeloid progenitor cells.

Post-transcriptional Gene Regulation by MicroRNA-194 Promotes Neuroendocrine Transdifferentiation in Prostate Cancer.

Potent therapeutic inhibition of the androgen receptor (AR) in prostate adenocarcinoma can lead to the emergence of neuroendocrine prostate cancer (NEPC), a phenomenon associated with enhanced cell plasticity. Here, we show that microRNA-194 (miR-194) is a regulator of epithelial-neuroendocrine transdifferentiation. In clinical prostate cancer samples, miR-194 expression and activity were elevated in NEPC and inversely correlated with AR signaling. miR-194 facilitated the emergence of neuroendocrine features in prostate cancer cells, a process mediated by its ability to directly target a suite of genes involved in cell plasticity. One such target was FOXA1, which encodes a transcription factor with a vital role in maintaining the prostate epithelial lineage. Importantly, a miR-194 inhibitor blocked epithelial-neuroendocrine transdifferentiation and inhibited the growth of cell lines and patient-derived organoids possessing neuroendocrine features. Overall, our study reveals a post-transcriptional mechanism regulating the plasticity of prostate cancer cells and provides a rationale for targeting miR-194 in NEPC.

A modular enhancer is differentially regulated by GATA and NFAT elements that direct different tissue-specific patterns of nucleosome positioning and inducible chromatin remodeling.

We investigated alternate mechanisms employed by enhancers to position and remodel nucleosomes and activate tissue-specific genes in divergent cell types. We demonstrated that the granulocyte-macrophage colony-stimulating factor (GM-CSF) gene enhancer is modular and recruits different sets of transcription factors in T cells and myeloid cells. The enhancer recruited distinct inducible tissue-specific enhanceosome-like complexes and directed nucleosomes to different positions in these cell types. In undifferentiated T cells, the enhancer was activated by inducible binding of two NFAT/AP-1 complexes which disrupted two specifically positioned nucleosomes (N1 and N2). In myeloid cells, the enhancer was remodeled by GATA factors which constitutively displaced an upstream nucleosome (N0) and cooperated with inducible AP-1 elements to activate transcription. In mast cells, which express both GATA-2 and NFAT, these two pathways combined to activate the enhancer and generate high-level gene expression. At least 5 kb of the GM-CSF locus was organized as an array of nucleosomes with fixed positions, but the enhancer adopted different nucleosome positions in T cells and mast cells. Furthermore, nucleosomes located between the enhancer and promoter were mobilized upon activation in an enhancer-dependent manner. These studies reveal that distinct tissue-specific mechanisms can be used either alternately or in combination to activate the same enhancer.

Insufficiently complex unique-molecular identifiers (UMIs) distort small RNA sequencing.

The attachment of unique molecular identifiers (UMIs) to RNA molecules prior to PCR amplification and sequencing, makes it possible to amplify libraries to a level that is sufficient to identify rare molecules, whilst simultaneously eliminating PCR bias through the identification of duplicated reads. Accurate de-duplication is dependent upon a sufficiently complex pool of UMIs to allow unique labelling. In applications dealing with complex libraries, such as total RNA-seq, only a limited variety of UMIs are required as the variation in molecules to be sequenced is enormous. However, when sequencing a less complex library, such as small RNAs for which there is a more limited range of possible sequences, we find increased variation in UMIs are required, even beyond that provided in a commercial kit specifically designed for the preparation of small RNA libraries for sequencing. We show that a pool of UMIs randomly varying across eight nucleotides is not of sufficient depth to uniquely tag the microRNAs to be sequenced. This results in over de-duplication of reads and the marked under-estimation of expression of the more abundant microRNAs. Whilst still arguing for the utility of UMIs, this work demonstrates the importance of their considered design to avoid errors in the estimation of gene expression in libraries derived from select regions of the transcriptome or small genomes.

The microRNA-200 family regulates epithelial to mesenchymal transition.

Epithelial-mesenchymal transition (EMT) is implicated in metastasis initiation and has recently been shown to be regulated by the miRNA-200 family and miR-205. Expression of these miRNAs was lost in invasive breast cancer cell lines displaying mesenchymal-like morphology suggesting these microRNAs may play a role in cancer metastasis.

Granulocyte-macrophage colony-stimulating factor enhancer activation requires cooperation between NFAT and AP-1 elements and is associated with extensive nucleosome reorganization.

The human granulocyte-macrophage colony-stimulating factor (GM-CSF) gene is activated by an NFAT-dependent enhancer forming an inducible DNase I hypersensitive (DH) site. The enhancer core comprising the DH site contains the GM330 and GM420 elements that bind NFAT and AP-1 cooperatively. Here we demonstrate that both elements are essential for enhancer activity and that Sp1 and AML1 sites in the enhancer become occupied in vivo only after activation. Chromatin structure analysis revealed that the GM-CSF enhancer core elements are divided between two adjacent nucleosomes that become destabilized and highly accessible after activation. Inducible chromatin reorganization was not restricted to the enhancer core but extended across a 3-kb domain of mobilized nucleosomes, within which the nucleosome repeat length was compressed from approximately 185 to 150 bp. The GM420 element is a high-affinity site that binds NFAT independently of AP-1 but depends on the linked AP-1 site for enhancer function. Nevertheless, just the NFAT motif from the GM420 element was sufficient to form a DH site within chromatin even in the absence of the AP-1 site. Hence, NFAT has the potential to cooperate with other transcription factors by promoting chromatin remodelling and increasing accessibility at inducible regulatory elements.

Phosphorylation of cold shock domain/Y-box proteins by ERK2 and GSK3beta and repression of the human VEGF promoter.

The hypoxia responsive region (HRR) of the VEGF promoter plays a key role in regulating VEGF expression. We found that the cold shock domain (Y-box) repressor proteins, dbpA and dbpB/YB-1, bind distinct strands of the human VEGF HRR. We find both dbpA and dbpB are phosphorylated by ERK2 and GSK3beta in vitro, and the binding of dbpB to single-strand VEGF HRR DNA is regulated by this phosphorylation. These findings suggest the ERK/MAPK and PI3K pathways may regulate VEGF expression in part through regulating the action of these repressor proteins.

The tyrosine phosphatase PTPN14 (Pez) inhibits metastasis by altering protein trafficking.

Factors secreted by tumor cells shape the local microenvironment to promote invasion and metastasis, as well as condition the premetastatic niche to enable secondary-site colonization and growth. In addition to this secretome, tumor cells have increased abundance of growth-promoting receptors at the cell surface. We found that the tyrosine phosphatase PTPN14 (also called Pez, which is mutated in various cancers) suppressed metastasis by reducing intracellular protein trafficking through the secretory pathway. Knocking down PTPN14 in tumor cells or injecting the peritoneum of mice with conditioned medium from PTPN14-deficient cell cultures promoted the growth and metastasis of breast cancer xenografts. Loss of catalytically functional PTPN14 increased the secretion of growth factors and cytokines, such as IL-8 (interleukin-8), and increased the abundance of EGFR (epidermal growth factor receptor) at the cell surface of breast cancer cells and of FLT4 (vascular endothelial growth factor receptor 3) at the cell surface of primary lymphatic endothelial cells. We identified RIN1 (Ras and Rab interactor 1) and PRKCD (protein kinase C-δ) as binding partners and substrates of PTPN14. Similar to cells overexpressing PTPN14, receptor trafficking to the cell surface was inhibited in cells that lacked PRKCD or RIN1 or expressed a nonphosphorylatable RIN1 mutant, and cytokine secretion was decreased in cells treated with PRKCD inhibitors. Invasive breast cancer tissue had decreased expression of PTPN14, and patient survival was worse when tumors had increased expression of the genes encoding RIN1 or PRKCD. Thus, PTPN14 prevents metastasis by restricting the trafficking of both soluble and membrane-bound proteins.

Down-regulation of the miRNA-200 family at the invasive front of colorectal cancers with degraded basement membrane indicates EMT is involved in cancer progression.

Cancer progression is a complex series of events thought to incorporate the reversible developmental process of epithelial-to-mesenchymal transition (EMT). In vitro, the microRNA-200 family maintains the epithelial phenotype by posttranscriptionally inhibiting the E-cadherin repressors, ZEB1 and ZEB2. Here, we used in situ hybridization and immunohistochemistry to assess expression of miR-200 and EMT biomarkers in formalin-fixed paraffin-embedded human colorectal adenocarcinomas. In addition, laser capture microdissection and quantitative real-time polymerase chain reaction were employed to quantify levels of miR-200 in the normal epithelium, tumor core, invasive front, and stroma. We find that miR-200 is downregulated at the invasive front of colorectal adenocarcinomas that have destroyed and invaded beyond the basement membrane. However, regional lymph node metastases and vascular carcinoma deposits show strong expression of miR-200, suggesting this family of miRNAs is involved in the recapitulation of the primary tumor phenotype at metastatic sites. In contrast, adenomas and adenocarcinomas with intact basement membranes showed uniform miR-200 expression from the tumor core to the tumor-host interface. Taken together, these data support the involvement of EMT and mesenchymal-to-epithelial transition (MET) in the metastasis cascade and show that miR-200 is downregulated in the initial stages of stromal invasion but is restored at metastatic sites.

Identification of an enhancer that increases miR-200b~200a~429 gene expression in breast cancer cells.

The miR-200b~200a~429 gene cluster is a key regulator of EMT and cancer metastasis, however the transcription-based mechanisms controlling its expression during this process are not well understood. We have analyzed the miR-200b~200a~429 locus for epigenetic modifications in breast epithelial and mesenchymal cell lines using chromatin immunoprecipitation assays and DNA methylation analysis. We discovered a novel enhancer located approximately 5.1kb upstream of the miR-200b~200a~429 transcriptional start site. This region was associated with the active enhancer chromatin signature comprising H3K4me1, H3K27ac, RNA polymerase II and CpG dinucleotide hypomethylation. Luciferase reporter assays revealed the upstream enhancer stimulated the transcription of the miR-200b~200a~429 minimal promoter region approximately 27-fold in breast epithelial cells. Furthermore, we found that a region of the enhancer was transcribed, producing a short, GC-rich, mainly nuclear, non-polyadenylated RNA transcript designated miR-200b eRNA. Over-expression of miR-200b eRNA had little effect on miR-200b~200a~429 promoter activity and its production did not correlate with miR-200b~200a~429 gene expression. While additional investigations of miR-200b eRNA function will be necessary, it is possible that miR-200b eRNA may be involved in the regulation of miR-200b~200a~429 gene expression and silencing. Taken together, these findings reveal the presence of a novel enhancer, which contributes to miR-200b~200a~429 transcriptional regulation in epithelial cells.