Search

Nanopore direct RNA Sequencing data contain information about the presence of RNA modifications, but their detection poses substantial challenges. Here the authors introduce Nanocompore, a new methodology for modification detection from Nanopore data.

Epithelial splicing regulatory proteins 1 and 2 (ESRP1/2) control the splicing pattern during epithelial to mesenchymal transition (EMT) in a physiological context and in cancer, including breast cancer (BC). Here, we report that ESRP1, but not ESRP2, is overexpressed in luminal BCs of patients with poor prognosis and correlates with estrogen receptor α (ERα) levels. Analysis of ERα genome-binding profiles in cell lines and primary breast tumors showed its binding in the proximity of ESRP1 and ESRP2 genes, whose expression is strongly decreased by ERα silencing in hormone-deprived conditions. The combined knock-down of ESRP1/2 in MCF-7 cells followed by RNA-Seq, revealed the dysregulation of 754 genes, with a widespread alteration of alternative splicing events (ASEs) of genes involved in cell signaling, metabolism, cell growth, and EMT. Functional network analysis of ASEs correlated with ESRP1/2 expression in ERα+ BCs showed RAC1 as the hub node in the protein–protein interactions altered by ESRP1/2 silencing. The comparison of ERα- and ESRP-modulated ASEs revealed 63 commonly regulated events, including 27 detected in primary BCs and endocrine-resistant cell lines. Our data support a functional implication of the ERα-ESRP1/2 axis in the onset and progression of BC by controlling the splicing patterns of related genes.

Abstract In the study of genomic regulation, strategies to integrate the data produced by Next Generation Sequencing (NGS)-based technologies in a meaningful ensemble are eagerly awaited and must continuously evolve. Here, we describe an integrative strategy for the analysis of data generated by chromatin immunoprecipitation followed by NGS which combines algorithms for data overlap, normalization and epigenetic state analysis. The performance of our strategy is illustrated by presenting the analysis of data relative to the transcriptional regulator Estrogen Receptor alpha (ERα) in MCF-7 breast cancer cells and of Glucocorticoid Receptor (GR) in A549 lung cancer cells. We went through the definition of reference cistromes for different experimental contexts, the integration of data relative to co-regulators and the overlay of chromatin states as defined by epigenetic marks in MCF-7 cells. With our strategy, we identified novel features of estrogen-independent ERα activity, including FoxM1 interaction, eRNAs transcription and a peculiar ontology of connected genes.

Recent improvements in cost-effectiveness of high-throughput technologies has allowed RNA sequencing of total transcriptomes suitable for evaluating the expression and regulation of circRNAs, a relatively novel class of transcript isoforms with suggested roles in transcriptional and post-transcriptional gene expression regulation, as well as their possible use as biomarkers, due to their deregulation in various human diseases. A limited number of integrated workflows exists for prediction, characterization, and differential expression analysis of circRNAs, none of them complying with computational reproducibility requirements. We developed Docker4Circ for the complete analysis of circRNAs from RNA-Seq data. Docker4Circ runs a comprehensive analysis of circRNAs in human and model organisms, including: circRNAs prediction; classification and annotation using six public databases; back-splice sequence reconstruction; internal alternative splicing of circularizing exons; alignment-free circRNAs quantification from RNA-Seq reads; and differential expression analysis. Docker4Circ makes circRNAs analysis easier and more accessible thanks to: (i) its R interface; (ii) encapsulation of computational tasks into docker images; (iii) user-friendly Java GUI Interface availability; and (iv) no need of advanced bash scripting skills for correct use. Furthermore, Docker4Circ ensures a reproducible analysis since all its tasks are embedded into a docker image following the guidelines provided by Reproducible Bioinformatics Project.

Tab2, originally described as a component of the inflammatory pathway, has been implicated in phenomena of gene de-repression in several contexts, due to its ability to interact with the NCoR corepressor. Tab2 interacts also with steroid receptors and dismisses NCoR from antagonist-bound Estrogen and Androgen Receptors on gene regulatory regions, thus modifying their transcriptional activity and leading to pharmacological resistance in breast and prostate cancer cells. We demonstrated previously that either Tab2 knock-down, or a peptide mimicking the Estrogen Receptor alpha domain interacting with Tab2, restore the antiproliferative response to Tamoxifen in Tamoxifen-resistant breast cancer cells. In this work, we map the domain of Tab2 responsible of Estrogen Receptor alpha interaction. First, using both co-immunoprecipitation and pull-down with recombinant proteins, we found that the central part of Tab2 is primarily responsible for this interaction, and that this region also interacts with Androgen Receptor. Then, we narrowed down the essential interaction region by means of competition assays using recombinant protein pull-down. The interaction motif was finally identified as a small region adjacent to, but not overlapping, the Tab2 MEKK1 phosphorylation sites. A synthetic peptide mimicking this motif efficiently displaced Tab2 from interacting with recombinant Estrogen Receptor alpha in vitro, prompting us to test its efficacy using derivatives of the MCF7 breast carcinoma cell lines that are spontaneously resistant to Tamoxifen. Indeed, we observed that this mimic peptide, made cell-permeable by addition of the TAT minimal carrier domain, reduced the growth of Tamoxifen-resistant MCF7 cells in the presence of Tamoxifen. These data indicate a novel functional domain of the Tab2 protein with potential application in drug design.

Estrogen receptor-α (ERα) is a ligand-inducible protein which mediates estrogenic hormones signaling and defines the luminal BC phenotype. Recently, we demonstrated that even in absence of ligands ERα (apoERα) binds chromatin sites where it regulates transcription of several protein-coding and lncRNA genes. Noteworthy, apoERα-regulated lncRNAs marginally overlap estrogen-induced transcripts, thus representing a new signature of luminal BC genes. By the analysis of H3K27ac enrichment in hormone-deprived MCF-7 cells, we defined a set of Super Enhancers (SEs) occupied by apoERα, including one mapped in proximity of the DSCAM-AS1 lncRNA gene. This represents a paradigm of apoERα activity since its expression is largely unaffected by estrogenic treatment, despite the fact that E2 increases ERα binding on DSCAM-AS1 promoter. We validated the enrichment of apoERα, p300, GATA3, FoxM1 and CTCF at both DSCAM-AS1 TSS and at its associated SE by ChIP-qPCR. Furthermore, by analyzing MCF-7 ChIA-PET data and by 3C assays, we confirmed long range chromatin interaction between the SE and the DSCAM-AS1 TSS. Interestingly, CTCF and p300 binding showed an enrichment in hormone-depleted medium and in the presence of ERα, elucidating the dynamics of the estrogen-independent regulation of DSCAM-AS1 expression. The analysis of this lncRNA provides a paradigm of transcriptional regulation of a luminal specific apoERα regulated lncRNA.