By integrating the expression profiles of long noncoding RNAs (lncRNAs) with clinical outcome and somatic copy-number alteration, the authors identified new lncRNAs that are associated with certain cancer subtypes and clinical prognoses. Experimental validation of the prostate cancer cell growth dependence of two new lncRNAs demonstrates the power of this approach for discovering disease-related lncRNAs. Despite growing appreciation of the importance of long noncoding RNAs (lncRNAs) in normal physiology and disease, our knowledge of cancer-related lncRNAs remains limited. By repurposing microarray probes, we constructed expression profiles of 10,207 lncRNA genes in approximately 1,300 tumors over four different cancer types. Through integrative analysis of the lncRNA expression profiles with clinical outcome and somatic copy-number alterations, we identified lncRNAs that are associated with cancer subtypes and clinical prognosis and predicted those that are potential drivers of cancer progression. We validated our predictions by experimentally confirming prostate cancer cell growth dependence on two newly identified lncRNAs. Our analysis provides a resource of clinically relevant lncRNAs for the development of lncRNA biomarkers and the identification of lncRNA therapeutic targets. It also demonstrates the power of integrating publically available genomic data sets and clinical information for discovering disease-associated lncRNAs.

The androgen receptor (AR) is a key factor that regulates the behavior and fate of prostate cancer cells. The AR-regulated network is activated when AR binds enhancer elements and modulates specific enhancer–promoter looping. Kallikrein-related peptidase 3 (KLK3), which codes for prostate-specific antigen (PSA), is a well-known AR-regulated gene and its upstream enhancers produce bidirectional enhancer RNAs (eRNAs), termed KLK3e. Here, we demonstrate that KLK3e facilitates the spatial interaction of the KLK3 enhancer and the KLK2 promoter and enhances long-distance KLK2 transcriptional activation. KLK3e carries the core enhancer element derived from the androgen response element III (ARE III), which is required for the interaction of AR and Mediator 1 (Med1). Furthermore, we show that KLK3e processes RNA-dependent enhancer activity depending on the integrity of core enhancer elements. The transcription of KLK3e was detectable and its expression is significantly correlated with KLK3 (R² = 0.6213, P < 5 × 10⁻¹¹) and KLK2 (R² = 0.5893, P < 5 × 10⁻¹⁰) in human prostate tissues. Interestingly, RNAi silencing of KLK3e resulted in a modest negative effect on prostate cancer cell proliferation. Accordingly, we report that an androgen-induced eRNA scaffolds the AR-associated protein complex that modulates chromosomal architecture and selectively enhances AR-dependent gene expression. Significance We report that enhancer RNAs (eRNAs), a class of long noncoding RNAs, participate in the androgen receptor (AR)-dependent looping complex that enhances spatial communication of distal enhancers and target promoters, leading to transcriptional activation events. Furthermore, our data show that KLK3 eRNA (KLK3e) selectively enhances the gene expression of AR-regulated genes, and provide evidence for a positive regulatory loop in which AR-dependent transcription is modulated by an intermediate eRNA. These findings may translate into improved RNA-based therapy (eRNA suppression) to enhance the durability of androgen deprivation therapy (ADT) and prediction of the efficacy of ADT by measuring the enhancer-derived activity (eRNA expression) in prostate tumors.

Significance: Alternative RNA splicing and the spliceosome machinery have been implicated in cancer progression. A genome-wide CRISPR screen identified the RNA processing factor heterogeneous nuclear ribonucleoprotein L (HNRNPL) as required for prostate cancer growth by regulating alternative RNA splicing and circular RNA formation. HNRNPL and its RNA clients are overexpressed during prostate cancer progression, supporting their potential role as therapeutic targets.

Mounting evidence suggests that long noncoding RNAs (lncRNAs) can function as microRNA sponges and compete for microRNA binding to protein-coding transcripts. However, the prevalence, functional significance and targets of lncRNA-mediated sponge regulation of cancer are mostly unknown. Here we identify a lncRNA-mediated sponge regulatory network that affects the expression of many protein-coding prostate cancer driver genes, by integrating analysis of sequence features and gene expression profiles of both lncRNAs and protein-coding genes in tumours. We confirm the tumour-suppressive function of two lncRNAs (TUG1 and CTB-89H12.4) and their regulation of PTEN expression in prostate cancer. Surprisingly, one of the two lncRNAs, TUG1, was previously known for its function in polycomb repressive complex 2 (PRC2)-mediated transcriptional regulation, suggesting its sub-cellular localization-dependent function. Our findings not only suggest an important role of lncRNA-mediated sponge regulation in cancer, but also underscore the critical influence of cytoplasmic localization on the efficacy of a sponge lncRNA.

Detailed analysis of DNase-seq protocols reveals the importance of choosing the right enzyme concentration and fragment length and cautions that many transcription factor footprints may represent cutting bias. Sequencing of DNase I hypersensitive sites (DNase-seq) is a powerful technique for identifying cis-regulatory elements across the genome. We studied the key experimental parameters to optimize performance of DNase-seq. Sequencing short fragments of 50–100 base pairs (bp) that accumulate in long internucleosome linker regions was more efficient for identifying transcription factor binding sites compared to sequencing longer fragments. We also assessed the potential of DNase-seq to predict transcription factor occupancy via generation of nucleotide-resolution transcription factor footprints. In modeling the sequence-specific DNase I cutting bias, we found a strong effect that varied over more than two orders of magnitude. This indicates that the nucleotide-resolution cleavage patterns at many transcription factor binding sites are derived from intrinsic DNase I cleavage bias rather than from specific protein-DNA interactions. In contrast, quantitative comparison of DNase I hypersensitivity between states can predict transcription factor occupancy associated with particular biological perturbations.

Smad7 plays an essential role in the negative-feedback regulation of transforming growth factor β (TGF-β) signaling by inhibiting TGF-β signaling at the receptor level. It can interfere with binding to type I receptors and thus activation of receptor-regulated Smads or recruit the E3 ubiquitin ligase Smurf to receptors and thus target them for degradation. Here, we report that Smad7 is predominantly localized in the nucleus of Hep3B cells. The targeted expression of Smad7 in the nucleus conferred superior inhibitory activity on TGF-β signaling, as determined by reporter assay in mammalian cells and by its effect on zebrafish embryogenesis. Furthermore, Smad7 repressed Smad3/4-, Smad2/4-, and Smad1/4-enhanced reporter gene expression, indicating that Smad7 can function independently of type I receptors. An oligonucleotide precipitation assay revealed that Smad7 can specifically bind to the Smad-responsive element via its MH2 domain, and DNA-binding activity was further confirmed in vivo with the promoter of PAI-1, a TGF-β target gene, by chromatin immunoprecipitation. Finally, we provide evidence that Smad7 disrupts the formation of the TGF-β-induced functional Smad-DNA complex. Our findings suggest that Smad7 inhibits TGF-β signaling in the nucleus by a novel mechanism.