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 specifc GapmeRs (or a control with scrambled GapmeR) was followed by diferentially expressed genes (DEGs) determination with Agilent microarrays and with Signifcance Analysis of Microarrays statistical test. DEGs were tested as a tumor risk classifer 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 signifcant expression upregulation/downregulation of hundreds of genes. Interestingly, PVT1 knockdown caused upregulation of 160 genes that were repressed by androgen, including a signifcantly 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 classifcation 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 signifcant 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 frst evidence that PVT1 was involved in signaling a genome-wide androgendependent transcriptional repressive program of tumor suppressor protein-coding genes in prostate cancer cells. Identifcation of transcriptional inhibition of tumor suppressor genes by PVT1 highlights the pathway to the investiga‑tion of mechanisms that lie behind the oncogenic role of PVT1 in cancer.

High-Quality Draft Genome Sequence Resources of Eight Xylella fastidiosa Strains Isolated from Citrus, Coffee, Plum, and Hibiscus in South America.

Xylella fastidiosa subsp. pauca, once confined to South America and infecting mainly citrus and coffee plants, has been found to be associated with other hosts and in other geographic regions. We present high-quality draft genome sequences of X. fastidiosa subsp. pauca strains J1a12, B111, U24D, and XRB isolated from citrus plants in Brazil, strain Fb7 isolated from a citrus plant in Argentina and strains 3124, Pr8x, and Hib4 isolated, respectively, from coffee, plum, and hibiscus plants in Brazil. Sequencing was performed using Roche 454-GS FLX, MiSeq-Illumina or Pacific Biosciences platforms. These high-quality genome assemblies will be useful for further studies about the genomic diversity, evolution, and biology of X. fastidiosa.

High-quality draft genome sequence resources of eight Xylella fastidiosa strains isolated from citrus, coffee, plum and hibiscus in South America.

Xylella fastidiosa subsp. pauca , once confined to South America and infecting mainly citrus and coffee plants, has been found to be associated with other hosts and in other geographic regions. We present high-quality draft genome sequences of X. fastidiosa subsp. pauca strains J1a12, B111, U24D and XRB isolated from citrus plants in Brazil, strain Fb7 isolated from a citrus plant in Argentina and strains 3124, Pr8x and Hib4 isolated, respectively, from coffee, plum and hibiscus plants in Brazil. Sequencing was performed using Roche 454-GS FLX, MiSeq-Illumina or Pacific Biosciences platforms. These high-quality genome assemblies will be useful for further studies about the genomic diversity, evolution, and biology of X. fastidiosa.

High-quality draft genome sequence resources of eight Xylella fastidiosa strains isolated from citrus, coffee, plum and hibiscus in South America.

Xylella fastidiosa subsp. pauca , once confined to South America and infecting mainly citrus and coffee plants, has been found to be associated with other hosts and in other geographic regions. We present high-quality draft genome sequences of X. fastidiosa subsp. pauca strains J1a12, B111, U24D and XRB isolated from citrus plants in Brazil, strain Fb7 isolated from a citrus plant in Argentina and strains 3124, Pr8x and Hib4 isolated, respectively, from coffee, plum and hibiscus plants in Brazil. Sequencing was performed using Roche 454-GS FLX, MiSeq-Illumina or Pacific Biosciences platforms. These high-quality genome assemblies will be useful for further studies about the genomic diversity, evolution, and biology of X. fastidiosa.

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 similar to 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.

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 similar to 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.

Exploring the Schistosoma mansoni adult male transcriptome using RNA-seq.

Schistosoma mansoni is one of the agents of schistosomiasis, a chronic and debilitating disease. Here we present a transcriptome-wide characterization of adult S. mansoni males by high-throughput RNA-sequencing. We obtained 1,620,432 high-quality ESTs from a directional strand-specific cDNA library, resulting in a 26% higher coverage of genome bases than that of the public ESTs available at NCBI. With a 15×-deep coverage of transcribed genomic regions, our data were able to (i) confirm for the first time 990 predictions without previous evidence of transcription; (ii) correct gene predictions; (iii) discover 989 and 1196 RNA-seq contigs that map to intergenic and intronic genomic regions, respectively, where no gene had been predicted before. These contigs could represent new protein-coding genes or non-coding RNAs (ncRNAs). Interestingly, we identified 11 novel Micro-exon genes (MEGs). These data reveal new features of the S. mansoni transcriptional landscape and significantly advance our understanding of the parasite transcriptome.

Recrutamento do complexo repressivo polycomb 2 pelo RNA não codificador longo antissenso ANRASSF1 modula a expressão do gene RASSF1A e a proliferação celular / Recruitment of polycomb repressive complex 2 by intronic long noncoding RNA ANRASSF1 modulates RASSF1A expression and cell proliferation

O gene supressor tumoral RASSF1A tem sido associado com redução da proliferação celular em diversos tumores. Sua expressão é regulada por eventos epigenéticos que envolvem o complexo repressivo polycomb (PRC2), no entanto os mecanismos moleculares da modulação do recrutamento deste modificador epigenético para este locus ainda são desconhecidos. Neste trabalho identificamos e caracterizamos ANRASSF1, um RNA não codificador longo (lncRNA) intrônico unspliced, que é transcrito na fita oposta do gene RASSF1A, em várias linhagem celulares e tecidos, e se liga a PRC2. ANRASSF1 é transcrito pela RNAPII, possui cap-5´ e cauda poli-A, além de localizar-se no núcleo e possuir uma meia-vida em média quatro vezes menor comparada com outros lncRNAs ligados à PRC2. A super-expressão ectópica de ANRASSF1 reduziu os níveis de RASSF1A e aumentou a taxa de proliferação em células HeLa, enquanto seu silenciamento provocou efeito oposto. Essas mudanças nos níveis de ANRASSF1 não afetaram a abundância da isoforma RASSF1C em nenhuma das condições. A super-expressão de ANRASSF1 provocou um grande aumento tanto da ocupação de PRC2 como da marca de histona repressiva H3K27me3 especificamente na região promotora RASSF1A. Nenhum efeito da super-expressão de ANRASSF1 foi detectado na ocupação de PRC2 e na histona H3K27me3 nas regiões promotoras de RASSF1C e de outros quatro genes vizinhos, incluindo dois genes supressores tumorais bem caracterizados. Além disso, foi demonstrado que ANRASSF1 forma um híbrido de RNA/DNA e recruta SUZ12, um componente do PRC2, para o promotor de RASSF1A. Notavelmente, foi detectado pelo ensaio de RNase-ChIP que a degradação de ANRASSF1 diminui a ocupação de PRC2 neste promotor. Esses resultados demonstram um novo mecanismo de repressão epigenética do supressor tumoral RASSF1A, envolvendo um lncRNA unspliced antissenso, onde ANRASSF1 reprime seletivamente a expressão da isoforma de RASSF1 que sobrepõe o transcrito antissenso de modo local e específico...

Tumor-suppressor RASSF1A gene down-regulation has been implicated in increasing cell proliferation in several tumors. Its expression is regulated by epigenetic events involving polycomb repressive complex 2 (PRC2), however the molecular mechanisms modulating recruitment of this epigenetic modifier to the locus remain largely unknown. Here, we identify and characterize ANRASSF1, an endogenous unspliced long noncoding RNA (lncRNA) that is transcribed from the opposite strand of RASSF1 gene in several cell lines and tissues, and binds to PRC2. ANRASSF1 is transcribed by RNA Polymerase II, 5'-capped, polyadenylated, displays nuclear localization, and has on average a four-fold shorter half-life compared to other lncRNAs that bind PRC2. ANRASSF1 ectopic overexpression decreases RASSF1A abundance and increases the proliferation rate of HeLa cells, whereas its silencing causes opposite effects. These changes in NRASSF1 levels do not affect RASSF1C isoform abundance. ANRASSF1 overexpression causes a marked increase both in PRC2 occupancy and in histone H3K27me3 repressive mark specifically at the RASSF1A promoter region. No effect of ANRASSF1 overexpression is detected on PRC2 occupancy and on histone H3K27me3 at the promoter regions of RASSF1C and of four other neighbor genes, including two well-characterized tumor suppressor genes. Additionally, we demonstrate that ANRASSF1 forms an RNA/DNA hybrid, and recruits SUZ12, a PRC2 component, to the RASSF1A promoter. Notably, depletion of ANRASSF1 disrupts SUZ12 occupancy on RASSF1A promoter as measured by RNAse-ChIP assay. Together, these results show a new mechanism of epigenetic repression of RASSF1A tumor suppressor gene involving an antisense unspliced lncRNA, in which ANRASSF1 selectively represses expression of the RASSF1 isoform overlapping the antisense transcript in a location-specific manner. In a broader perspective, our findings suggest that other non-characterized unspliced intronic lncRNAs transcribed in the...