Dissecting the genomic activity of a transcriptional regulator by the integrative analysis of omics data

Novel estrogen-responsive genes (ERGs) for the evaluation of estrogenic activity

Estrogen action is mediated by various genes, including estrogen-responsive genes (ERGs). ERGs have been used as reporter-genes and markers for gene expression. Gene expression profiling using a set of ERGs has been used to examine statistically reliable transcriptomic assays such as DNA microarray assays and RNA sequencing (RNA-seq). However, the quality of ERGs has not been extensively examined. Here, we obtained a set of 300 ERGs that were newly identified by six sets of RNA-seq data from estrogen-treated and control human breast cancer MCF-7 cells. The ERGs exhibited statistical stability, which was based on the coefficient of variation (CV) analysis, correlation analysis, and examination of the functional association with estrogen action using database searches. A set of the top 30 genes based on CV ranking were further evaluated quantitatively by RT-PCR and qualitatively by a functional analysis using the GO and KEGG databases and by a mechanistic analysis to classify ERα/β-dependent or ER-independent types of transcriptional regulation. The 30 ERGs were characterized according to (1) the enzymes, such as metabolic enzymes, proteases, and protein kinases, (2) the genes with specific cell functions, such as cell-signaling mediators, tumor-suppressors, and the roles in breast cancer, (3) the association with transcriptional regulation, and (4) estrogen-responsiveness. Therefore, the ERGs identified here represent various cell functions and cell signaling pathways, including estrogen signaling, and thus, may be useful to evaluate estrogenic activity.

A Regulatory Axis between Epithelial Splicing Regulatory Proteins and Estrogen Receptor α Modulates the Alternative Transcriptome of Luminal Breast Cancer

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.

The Estrogen Receptor α Signaling Pathway Controls Alternative Splicing in the Absence of Ligands in Breast Cancer Cells

Background: The transcriptional activity of estrogen receptor α (ERα) in breast cancer (BC) is extensively characterized. Our group has previously shown that ERα controls the expression of a number of genes in its unliganded form (apoERα), among which a large group of RNA-binding proteins (RBPs) encode genes, suggesting its role in the control of co- and post-transcriptional events. Methods: apoERα-mediated RNA processing events were characterized by the analysis of transcript usage and alternative splicing changes in an RNA-sequencing dataset from MCF-7 cells after siRNA-induced ERα downregulation. Results: ApoERα depletion induced an expression change of 681 RBPs, including 84 splicing factors involved in translation, ribonucleoprotein complex assembly, and 3′end processing. ApoERα depletion results in 758 isoform switching events with effects on 3′end length and the splicing of alternative cassette exons. The functional enrichment of these events shows that post-transcriptional regulation is part of the mechanisms by which apoERα controls epithelial-to-mesenchymal transition and BC cell proliferation. In primary BCs, the inclusion levels of the experimentally identified alternatively spliced exons are associated with overall and disease-free survival. Conclusion: Our data supports the role of apoERα in maintaining the luminal phenotype of BC cells by extensively regulating gene expression at the alternative splicing level.

Steroid Receptor Signallings as Targets for Resveratrol Actions in Breast and Prostate Cancer

Extensive research over the past 25 years in hormone-dependent cancers, such as breast cancer and prostate cancer, has identified the molecular mechanisms driven by steroid receptors, elucidating the interplay between genomic and non-genomic steroid receptors mechanism of action. Altogether, these mechanisms create the specific gene expression programs that contribute to endocrine therapy resistance and cancer progression. These findings, on the bidirectional molecular crosstalk between steroid and growth factor receptors pathways in endocrine resistance, suggest the use of multi-target inhibitors together with endocrine therapies, for treating resistant disease. In this review we will discuss the novel understanding on the chemopreventive and anti-cancer activities of Resveratrol (3,5,4′-trihydroxy-stilbene) (RSV), a phytoalexin found in grapes acting on a plethora of targets. We will highlight Resveratrol effect on steroid receptors signalling and its potential use in the treatment of hormone-dependent cancer. Understanding the molecular mechanisms by which the bioactive compound influences cancer cell behaviour, by interfering with steroid receptors functional activity, will help to advance the design of combination strategies to increase the rate of complete and durable clinical response in patients.

Luminal breast cancer-specific circular RNAs uncovered by a novel tool for data analysis

Circular RNAs are highly stable molecules present in all eukaryotes generated by distinct transcript processing. We have exploited poly(A-) RNA-Seq data generated in our lab in MCF-7 breast cancer cells to define a compilation of exonic circRNAs more comprehensive than previously existing lists. Development of a novel computational tool, named CircHunter, allowed us to more accurately characterize circRNAs and to quantitatively evaluate their expression in publicly available RNA-Seq data from breast cancer cell lines and tumor tissues. We observed and confirmed, by ChIP analysis, that exons involved in circularization events display significantly higher levels of the histone post-transcriptional modification H3K36me3 than non-circularizing exons. This result has potential impact on circRNA biogenesis since H3K36me3 has been involved in alternative splicing mechanisms. By analyzing an Ago-HITS-CLIP dataset we also found that circularizing exons overlapped with an unexpectedly higher number of Ago binding sites than non-circularizing exons. Finally, we observed that a subset of MCF-7 circRNAs are specific to tumor versus normal tissue, while others can distinguish Luminal from other tumor subtypes, thus suggesting that circRNAs can be exploited as novel biomarkers and drug targets for breast cancer.

Luminal lncRNAs Regulation by ERα-Controlled Enhancers in a Ligand-Independent Manner in Breast Cancer Cells

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.