PVT1 signals an androgen-dependent transcriptional repression program in prostate cancer cells and a set of the repressed genes predicts high-risk tumors.

BACKGROUND: Androgen receptor (AR) and polycomb repressive complex 2 (PRC2) are known to co-occupy the loci of genes that are downregulated by androgen-stimulus. Long intergenic non-coding RNA (lincRNA) PVT1 is an overexpressed oncogene that is associated with AR in LNCaP prostate cancer cells, and with PRC2 in HeLa and many other types of cancer cells. The possible involvement of PVT1 in mediating androgen-induced gene expression downregulation in prostate cancer has not been explored. METHODS: LNCaP cell line was used. Native RNA-binding-protein immunoprecipitation with anti-AR or anti-EZH2 was followed by RT-qPCR with primers for PVT1. Knockdown of PVT1 with specific GapmeRs (or a control with scrambled GapmeR) was followed by differentially expressed genes (DEGs) determination with Agilent microarrays and with Significance Analysis of Microarrays statistical test. DEGs were tested as a tumor risk classifier with a machine learning Random Forest algorithm run with gene expression data from all TCGA-PRAD (prostate adenocarcinoma) tumors as input. ChIP-qPCR was performed for histone marks at the promoter of one DEG. RESULTS: We show that PVT1 knockdown in androgen-stimulated LNCaP cells caused statistically significant expression upregulation/downregulation of hundreds of genes. Interestingly, PVT1 knockdown caused upregulation of 160 genes that were repressed by androgen, including a significantly enriched set of tumor suppressor genes, and among them FAS, NOV/CCN3, BMF, HRK, IFIT2, AJUBA, DRAIC and TNFRSF21. A 121-gene-set (out of the 160) was able to correctly predict the classification of all 293 intermediate- and high-risk TCGA-PRAD tumors, with a mean ROC area under the curve AUC = 0.89 ± 0.04, pointing to the relevance of these genes in cancer aggressiveness. Native RIP-qPCR in LNCaP showed that PVT1 was associated with EZH2, a component of PRC2. PVT1 knockdown followed by ChIP-qPCR showed significant epigenetic remodeling at the enhancer and promoter regions of tumor suppressor gene NOV, one of the androgen-repressed genes that were upregulated upon PVT1 silencing. CONCLUSIONS: Overall, we provide first evidence that PVT1 was involved in signaling a genome-wide androgen-dependent transcriptional repressive program of tumor suppressor protein-coding genes in prostate cancer cells. Identification of transcriptional inhibition of tumor suppressor genes by PVT1 highlights the pathway to the investigation of mechanisms that lie behind the oncogenic role of PVT1 in cancer. Video Abstract.

Atlas of Schistosoma mansoni long non-coding RNAs and their expression correlation to protein-coding genes.

Long non-coding RNAs (lncRNAs) have been widely discovered in several organisms with the help of high-throughput RNA sequencing. LncRNAs are over 200 nt-long transcripts that do not have protein-coding (PC) potential, having been reported in model organisms to act mainly on the overall control of PC gene expression. Little is known about the functionality of lncRNAs in evolutionarily ancient non-model metazoan organisms, like Schistosoma mansoni, the parasite that causes schistosomiasis, one of the most prevalent infectious-parasitic diseases worldwide. In a recent transcriptomics effort, we identified thousands of S. mansoni lncRNAs predicted to be functional along the course of parasite development. Here, we present an online catalog of each of the S. mansoni lncRNAs whose expression is correlated to PC genes along the parasite life-cycle, which can be conveniently browsed and downloaded through a new web resource http://verjolab.usp.br. We also provide access now to navigation on the co-expression networks disclosed in our previous publication, where we correlated mRNAs and lncRNAs transcriptional patterns across five life-cycle stages/forms, pinpointing biological processes where lncRNAs might act upon.Database URL: http://verjolab.usp.br.

Chromatin Landscape Distinguishes the Genomic Loci of Hundreds of Androgen-Receptor-Associated LincRNAs From the Loci of Non-associated LincRNAs.

Cell signaling events triggered by androgen hormone in prostate cells is dependent on activation of the androgen receptor (AR) transcription factor. Androgen hormone binding to AR promotes its displacement from the cytoplasm to the nucleus and AR binding to DNA motifs, thus inducing activatory and inhibitory transcriptional programs through a complex regulatory mechanism not yet fully understood. In this work, we performed RNA-seq deep-sequencing of LNCaP prostate cancer cells and found over 7000 expressed long intergenic non-coding RNAs (lincRNAs), of which ∼4000 are novel lincRNAs, and 258 lincRNAs have their expression activated by androgen. Immunoprecipitation of AR, followed by large-scale sequencing of co-immunoprecipitated RNAs (RIP-Seq) has identified in the LNCaP cell line a total of 619 lincRNAs that were significantly enriched (FDR < 10%, DESeq2) in the anti-Androgen Receptor (antiAR) fraction in relation to the control fraction (non-specific IgG), and we named them Androgen-Receptor-Associated lincRNAs (ARA-lincRNAs). A genome-wide analysis showed that protein-coding gene neighbors to ARA-lincRNAs had a significantly higher androgen-induced change in expression than protein-coding genes neighboring lincRNAs not associated to AR. To find relevant epigenetic signatures enriched at the ARA-lincRNAs' transcription start sites (TSSs) we used a machine learning approach and identified that the ARA-lincRNA genomic loci in LNCaP cells are significantly enriched with epigenetic marks that are characteristic of in cis enhancer RNA regulators, and that the H3K27ac mark of active enhancers is conspicuously enriched at the TSS of ARA-lincRNAs adjacent to androgen-activated protein-coding genes. In addition, LNCaP topologically associating domains (TADs) that comprise chromatin regions with ARA-lincRNAs exhibit transcription factor contents, epigenetic marks and gene transcriptional activities that are significantly different from TADs not containing ARA-lincRNAs. This work highlights the possible involvement of hundreds of lincRNAs working in synergy with the AR on the genome-wide androgen-induced gene regulatory program in prostate cells.

Inhibition of histone methyltransferase EZH2 in Schistosoma mansoni in vitro by GSK343 reduces egg laying and decreases the expression of genes implicated in DNA replication and noncoding RNA metabolism.

BACKGROUND: The possibility of emergence of praziquantel-resistant Schistosoma parasites and the lack of other effective drugs demand the discovery of new schistosomicidal agents. In this context the study of compounds that target histone-modifying enzymes is extremely promising. Our aim was to investigate the effect of inhibition of EZH2, a histone methyltransferase that is involved in chromatin remodeling processes and gene expression control; we tested different developmental forms of Schistosoma mansoni using GKS343, a selective inhibitor of EZH2 in human cells. METHODOLOGY/PRINCIPAL FINDINGS: Adult male and female worms and schistosomula were treated with different concentrations of GSK343 for up to two days in vitro. Western blotting showed a decrease in the H3K27me3 histone mark in all three developmental forms. Motility, mortality, pairing and egg laying were employed as schistosomicidal parameters for adult worms. Schistosomula viability was evaluated with propidium iodide staining and ATP quantification. Adult worms showed decreased motility when exposed to GSK343. Also, an approximate 40% reduction of egg laying by GSK343-treated females was observed when compared with controls (0.1% DMSO). Scanning electron microscopy showed the formation of bulges and bubbles throughout the dorsal region of GSK343-treated adult worms. In schistosomula the body was extremely contracted with the presence of numerous folds, and growth was markedly slowed. RNA-seq was applied to identify the metabolic pathways affected by GSK343 sublethal doses. GSK343-treated adult worms showed significantly altered expression of genes related to transmembrane transport, cellular homeostasis and egg development. In females, genes related to DNA replication and noncoding RNA metabolism processes were downregulated. Schistosomula showed altered expression of genes related to cell adhesion and membrane synthesis pathways. CONCLUSIONS/SIGNIFICANCE: The results indicated that GSK343 presents in vitro activities against S. mansoni, and the characterization of EZH2 as a new potential molecular target establishes EZH2 inhibitors as part of a promising new group of compounds that could be used for the development of schistosomicidal agents.

The Schistosoma mansoni genome encodes thousands of long non-coding RNAs predicted to be functional at different parasite life-cycle stages.

Next Generation Sequencing (NGS) strategies, like RNA-Seq, have revealed the transcription of a wide variety of long non-coding RNAs (lncRNAs) in the genomes of several organisms. In the present work we assessed the lncRNAs complement of Schistosoma mansoni, the blood fluke that causes schistosomiasis, ranked among the most prevalent parasitic diseases worldwide. We focused on the long intergenic/intervening ncRNAs (lincRNAs), hidden within the large amount of information obtained through RNA-Seq in S. mansoni (88 libraries). Our computational pipeline identified 7029 canonically-spliced putative lincRNA genes on 2596 genomic loci (at an average 2.7 isoforms per lincRNA locus), as well as 402 spliced lncRNAs that are antisense to protein-coding (PC) genes. Hundreds of lincRNAs showed traits for being functional, such as the presence of epigenetic marks at their transcription start sites, evolutionary conservation among other schistosome species and differential expression across five different life-cycle stages of the parasite. Real-time qPCR has confirmed the differential life-cycle stage expression of a set of selected lincRNAs. We have built PC gene and lincRNA co-expression networks, unraveling key biological processes where lincRNAs might be involved during parasite development. This is the first report of a large-scale identification and structural annotation of lncRNAs in the S. mansoni genome.