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.

Docker4Circ: A Framework for the Reproducible Characterization of circRNAs from RNA-Seq Data.

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.

The expression of LINE1-MET chimeric transcript identifies a subgroup of aggressive breast cancers.

Demethylation of the long interspersed nuclear element (LINE-1; L1) antisense promoter can result in transcription of neighboring sequences as for the L1-MET transcript produced by the L1 placed in the second intron of MET. To define the role of L1-MET, we investigated the sequence and the transcription of L1-MET in vitro models and heterogeneous breast cancers, previously reported to show other L1-derived transcripts. L1-MET expressing cell lines were initially identified in silico and investigated for L1-MET promoter methylation, cDNA sequence and cell fraction mRNA. The transcriptional level of L1-MET and MET were then evaluated in breast specimens, including 9 cancer cell lines, 41 carcinomas of different subtypes, and 11 normal tissues. In addition to a L1-MET transcript ending at MET exon 21, six novel L1-MET splice variants were identified. Normal breast tissues were negative for the L1-MET expression, whereas the triple-negative breast cancer (TNBC) and the high-grade carcinomas were enriched with the L1-MET mRNA (p = 0.005 and p = 0.018, respectively). In cancer cells and tissues the L1-MET expression was associated with its promoter hypomethylation (ρ = -0.8 and -0.9, respectively). No correlation was found between L1-MET and MET mRNA although L1-MET expressing tumors with higher L1-MET/MET ratio were negative for the MET protein expression (p = 0.006). Besides providing the first identification and detailed description of L1-MET in breast cancer, we clearly demonstrate that higher levels of this transcript specifically recognize a subset of more aggressive carcinomas, mainly TNBC. We suggest the possible evaluation of L1-MET in the challenging diagnosis of early TNBCs.

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.

Genome-wide activity of unliganded estrogen receptor-α in breast cancer cells.

Estrogen receptor-α (ERα) has central role in hormone-dependent breast cancer and its ligand-induced functions have been extensively characterized. However, evidence exists that ERα has functions that are independent of ligands. In the present work, we investigated the binding of ERα to chromatin in the absence of ligands and its functions on gene regulation. We demonstrated that in MCF7 breast cancer cells unliganded ERα binds to more than 4,000 chromatin sites. Unexpectedly, although almost entirely comprised in the larger group of estrogen-induced binding sites, we found that unliganded-ERα binding is specifically linked to genes with developmental functions, compared with estrogen-induced binding. Moreover, we found that siRNA-mediated down-regulation of ERα in absence of estrogen is accompanied by changes in the expression levels of hundreds of coding and noncoding RNAs. Down-regulated mRNAs showed enrichment in genes related to epithelial cell growth and development. Stable ERα down-regulation using shRNA, which caused cell growth arrest, was accompanied by increased H3K27me3 at ERα binding sites. Finally, we found that FOXA1 and AP2γ binding to several sites is decreased upon ERα silencing, suggesting that unliganded ERα participates, together with other factors, in the maintenance of the luminal-specific cistrome in breast cancer cells.

Genomic lens on neuroglobin transcription.

Neuroglobin is a brain globin with neuroprotective effects against ischemia and related pathological processes, acting as O2 sensor and antiapoptotic pathway transducer. Here, we survey data on neuroanatomical coexpression of transcription factors, epigenetic signature and predictive transcription factor binding sites at the neuroglobin gene locus to find hints of pathways to neuroglobin transcriptional regulation. These data provide a glimpse of how neuroglobin expression may translate into neuronal diversity and function, as well as disease.

miR148b is a major coordinator of breast cancer progression in a relapse-associated microRNA signature by targeting ITGA5, ROCK1, PIK3CA, NRAS, and CSF1.

Breast cancer is often fatal during its metastatic dissemination. To unravel the role of microRNAs (miRs) during malignancy, we analyzed miR expression in 77 primary breast carcinomas and identified 16 relapse-associated miRs that correlate with survival and/or distinguish tumor subtypes in different datasets. Among them, miR-148b, down-regulated in aggressive breast tumors, was found to be a major coordinator of malignancy. In fact, it is able to oppose various steps of tumor progression when overexpressed in cell lines by influencing invasion, survival to anoikis, extravasation, lung metastasis formation, and chemotherapy response. miR-148b controls malignancy by coordinating a novel pathway involving over 130 genes and, in particular, it directly targets players of the integrin signaling, such as ITGA5, ROCK1, PIK3CA/p110α, and NRAS, as well as CSF1, a growth factor for stroma cells. Our findings reveal the importance of the identified 16 miRs for disease outcome predictions and suggest a critical role for miR-148b in the control of breast cancer progression.

The role of Transposable Elements in shaping the combinatorial interaction of Transcription Factors.

BACKGROUND: In the last few years several studies have shown that Transposable Elements (TEs) in the human genome are significantly associated with Transcription Factor Binding Sites (TFBSs) and that in several cases their expansion within the genome led to a substantial rewiring of the regulatory network. Another important feature of the regulatory network which has been thoroughly studied is the combinatorial organization of transcriptional regulation. In this paper we combine these two observations and suggest that TEs, besides rewiring the network, also played a central role in the evolution of particular patterns of combinatorial gene regulation. RESULTS: To address this issue we searched for TEs overlapping Estrogen Receptor α (ERα) binding peaks in two publicly available ChIP-seq datasets from the MCF7 cell line corresponding to different modalities of exposure to estrogen. We found a remarkable enrichment of a few specific classes of Transposons. Among these a prominent role was played by MIR (Mammalian Interspersed Repeats) transposons. These TEs underwent a dramatic expansion at the beginning of the mammalian radiation and then stabilized. We conjecture that the special affinity of ERα for the MIR class of TEs could be at the origin of the important role assumed by ERα in Mammalians. We then searched for TFBSs within the TEs overlapping ChIP-seq peaks. We found a strong enrichment of a few precise combinations of TFBS. In several cases the corresponding Transcription Factors (TFs) were known cofactors of ERα, thus supporting the idea of a co-regulatory role of TFBS within the same TE. Moreover, most of these correlations turned out to be strictly associated to specific classes of TEs thus suggesting the presence of a well-defined "transposon code" within the regulatory network. CONCLUSIONS: In this work we tried to shed light into the role of Transposable Elements (TEs) in shaping the regulatory network of higher eukaryotes. To test this idea we focused on a particular transcription factor: the Estrogen Receptor α (ERα) and we found that ERα preferentially targets a well defined set of TEs and that these TEs host combinations of transcriptional regulators involving several of known co-regulators of ERα. Moreover, a significant number of these TEs turned out to be conserved between human and mouse and located in the vicinity (and thus candidate to be regulators) of important estrogen-related genes.

Glucocorticoid receptor activity discriminates between progesterone and medroxyprogesterone acetate effects in breast cells.

The purpose of this article is to determine the tumorigenic potential of estradiol treatment (E2) when combined with either progesterone (P4) or medroxyprogesterone acetate (MPA) in normal luminal human breast cells (HBE) and in human breast cancer cells (T47-D, MCF-7). Proliferation profiles were evaluated, along with the gene transactivation activity between the progesterone and glucocorticoid receptors (PR, GR) in HBE, T47-D, and MCF-7 cells treated by E2 + P4 or E2 + MPA. High throughput transcriptome analysis was performed on RNA from HBE cells treated by E2, E2 + MPA and E2 + P4. GR content was analyzed in normal breast cells as well. In HBE cells, E2 + P4 treatment was antiproliferative and promoted cellular differentiation. In contrast, E2 + MPA displayed mitogenic, antiapoptotic effects in HBE cells and did not influence cellular differentiation. The effect of P4 and MPA on cell proliferation was, however, variable in breast cancer cells. In cells containing GR or/and PR, MPA decreased proliferation whereas P4 antiproliferative effect needed the presence of PR. In HBE cells, the regulation of genes by E2 + P4, and E2 + MPA was significantly different, particularly in cell proliferation and cell death gene families. Further analysis revealed a modulation of the glucocorticoid receptor gene expression pathway by E2 + MPA. Predominant MPA glucocorticoid activity in normal and breast cancer cells was demonstrated using a glucocorticoid antagonist and the down-regulation of the GR by RNA interference. In normal luminal breast cells and in breast cancer cells, P4 and MPA combined with E2 treatment have opposing mitogenic effects due to GR. The consequences of MPA glucocorticoid potencies as well as the importance of GR in breast tissue merit a reappraisal.

ERalpha as ligand-independent activator of CDH-1 regulates determination and maintenance of epithelial morphology in breast cancer cells.

Estrogen receptor alpha (ERalpha) and E-cadherin are primary markers of luminal epithelial breast cancer cells with E-cadherin being a main caretaker of the epithelial phenotype. E-cadherin repression is needed for cancer cells to acquire motile and invasive properties, and it is known that in ER-positive breast cancer cells, estrogen down-regulate E-cadherin gene transcription. We report here that ERalpha is bound to the E-cadherin promoter in both the presence and the complete absence of estrogen, suggesting an unexpected role for unliganded ERalpha in E-cadherin transcription. Indeed, our data reveal that activation by unliganded ERalpha and repression by estrogen-activated ERalpha require direct binding to a half-estrogen response element within the E-cadherin promoter and exchange from associated coactivators to corepressors. Therefore, these results suggest a pivotal role for unliganded ERalpha in controlling a fundamental caretaker of the epithelial phenotype in breast cancer cells. Here, we show that ERalpha-positive breast cancer T47D cells transduced with the sfRON kinase undergo a full epithelial-mesenchymal conversion and lose E-cadherin and ERalpha expression. Our data show that, although the E-cadherin gene becomes hypermethylated and heterochromatic, kinase inhibitors can restore E-cadherin expression, together with an epithelial morphology in an ERalpha-dependent fashion. Similarly, transfection of ERalpha, in the absence of ligands, was sufficient to restore E-cadherin transcription in both sfRON-T47D and other ERalpha-, E-cadherin-negative cells. Therefore, our results suggest a novel role for the ERalpha that plays the dual role of ligand-independent activator and ligand-dependent repressor of E-cadherin in breast cancer cells.

Estrogen receptor alpha controls a gene network in luminal-like breast cancer cells comprising multiple transcription factors and microRNAs.

Luminal-like breast tumor cells express estrogen receptor alpha (ERalpha), a member of the nuclear receptor family of ligand-activated transcription factors that controls their proliferation, survival, and functional status. To identify the molecular determinants of this hormone-responsive tumor phenotype, a comprehensive genome-wide analysis was performed in estrogen stimulated MCF-7 and ZR-75.1 cells by integrating time-course mRNA expression profiling with global mapping of genomic ERalpha binding sites by chromatin immunoprecipitation coupled to massively parallel sequencing, microRNA expression profiling, and in silico analysis of transcription units and receptor binding regions identified. All 1270 genes that were found to respond to 17beta-estradiol in both cell lines cluster in 33 highly concordant groups, each of which showed defined kinetics of RNA changes. This hormone-responsive gene set includes several direct targets of ERalpha and is organized in a gene regulation cascade, stemming from ligand-activated receptor and reaching a large number of downstream targets via AP-2gamma, B-cell activating transcription factor, E2F1 and 2, E74-like factor 3, GTF2IRD1, hairy and enhancer of split homologue-1, MYB, SMAD3, RARalpha, and RXRalpha transcription factors. MicroRNAs are also integral components of this gene regulation network because miR-107, miR-424, miR-570, miR-618, and miR-760 are regulated by 17beta-estradiol along with other microRNAs that can target a significant number of transcripts belonging to one or more estrogen-responsive gene clusters.

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.

Computational Analysis of circRNA Expression Data.

Analysis of circular RNA (circRNA) expression from RNA-Seq data can be performed with different algorithms and analysis pipelines, tools allowing the extraction of heterogeneous information on the expression of this novel class of RNAs. Computational pipelines were developed to facilitate the analysis of circRNA expression by leveraging different public tools in easy-to-use pipelines. This chapter describes the complete workflow for a computationally reproducible analysis of circRNA expression starting for a public RNA-Seq experiment. The main steps of circRNA prediction, annotation, classification, sequence reconstruction, quantification, and differential expression are illustrated.

CircuitsDB: a database of mixed microRNA/transcription factor feed-forward regulatory circuits in human and mouse.

BACKGROUND: Transcription Factors (TFs) and microRNAs (miRNAs) are key players for gene expression regulation in higher eukaryotes. In the last years, a large amount of bioinformatic studies were devoted to the elucidation of transcriptional and post-transcriptional (mostly miRNA-mediated) regulatory interactions, but little is known about the interplay between them. DESCRIPTION: Here we describe a dynamic web-accessible database, CircuitsDB, supporting a genome-wide transcriptional and post-transcriptional regulatory network integration, for the human and mouse genomes, based on a bioinformatic sequence-analysis approach. In particular, CircuitsDB is currently focused on the study of mixed miRNA/TF Feed-Forward regulatory Loops (FFLs), i.e. elementary circuits in which a master TF regulates an miRNA and together with it a set of Joint Target protein-coding genes. The database was constructed using an ab-initio oligo analysis procedure for the identification of the transcriptional and post-transcriptional interactions. Several external sources of information were then pooled together to obtain the functional annotation of the proposed interactions. Results for human and mouse genomes are presented in an integrated web tool, that allows users to explore the circuits, investigate their sequence and functional properties and thus suggest possible biological experiments. CONCLUSIONS: We present CircuitsDB, a web-server devoted to the study of human and mouse mixed miRNA/TF Feed-Forward regulatory circuits, freely available at: http://biocluster.di.unito.it/circuits/

AP-2alpha regulates migration of GN-11 neurons via a specific genetic programme involving the Axl receptor tyrosine kinase.

BACKGROUND: Neuronal migration is a crucial process that allows neurons to reach their correct target location to allow the nervous system to function properly. AP-2alpha is a transcription factor essential for neural crest cell migration and its mutation results in apoptosis within this cell population, as demonstrated by genetic models. RESULTS: We down-modulated AP-2alpha expression in GN-11 neurons by RNA interference and observe reduced neuron migration following the activation of a specific genetic programme including the Adhesion Related Kinase (Axl) gene. We prove that Axl is able to coordinate migration per se and by ChIP and promoter analysis we observe that its transcription is directly driven by AP-2alpha via the binding to one or more functional AP-2alpha binding sites present in its regulatory region. Analysis of migration in AP-2alpha null mouse embryo fibroblasts also reveals an essential role for AP-2alpha in cell movement via the activation of a distinct genetic programme. CONCLUSION: We show that AP-2alpha plays an essential role in cell movement via the activation of cell-specific genetic programmes. Moreover, we demonstrate that the AP-2alpha regulated gene Axl is an essential player in GN-11 neuron migration.

DSCAM-AS1-Driven Proliferation of Breast Cancer Cells Involves Regulation of Alternative Exon Splicing and 3'-End Usage.

DSCAM-AS1 is a cancer-related long noncoding RNA with higher expression levels in Luminal A, B, and HER2-positive Breast Carcinoma (BC), where its expression is strongly dependent on Estrogen Receptor Alpha (ERα). DSCAM-AS1 expression is analyzed in 30 public datasets and, additionally, by qRT-PCR in tumors from 93 BC patients, to uncover correlations with clinical data. Moreover, the effect of DSCAM-AS1 knockdown on gene expression and alternative splicing is studied by RNA-Seq in MCF-7 cells. We confirm DSCAM-AS1 overexpression in high grade Luminal A, B, and HER2+ BCs and find a significant correlation with disease relapse. In total, 908 genes are regulated by DSCAM-AS1-silencing, primarily involved in the cell cycle and inflammatory response. Noteworthily, the analysis of alternative splicing and isoform regulation reveals 2085 splicing events regulated by DSCAM-AS1, enriched in alternative polyadenylation sites, 3'UTR (untranslated region) shortening and exon skipping events. Finally, the DSCAM-AS1-interacting splicing factor heterogeneous nuclear ribonucleoprotein L (hnRNPL) is predicted as the most enriched RBP for exon skipping and 3'UTR events. The relevance of DSCAM-AS1 overexpression in BC is confirmed by clinical data and further enhanced by its possible involvement in the regulation of RNA processing, which is emerging as one of the most important dysfunctions in cancer.

AP-2alpha and AP-2gamma regulate tumor progression via specific genetic programs.

The events occurring during tumor formation and progression display similarities to some of the steps in embryonic morphogenesis. The family of AP-2 proteins consists of five different transcription factors (alpha, beta, gamma, delta, and epsilon) that play relevant roles in embryonic development, as demonstrated by the phenotypes of the corresponding knockout mice. Here, we show that AP-2alpha and AP-2gamma proteins play an essential role in tumorigenesis. Down-modulation of AP-2 expression in tumor cells by RNA interference (RNAi) led to enhanced tumor growth and reduced chemotherapy-induced cell death, as well as migration and invasion. Most of these biological modulations were rescued by AP-2 overexpression. We observed that increased xenotransplant growth was mostly due to highly enhanced proliferation of the tumor cells together with reduced innate immune cell recruitment. Moreover, we showed that migration impairment was mediated, at least in part, by secreted factors. To identify the genetic programs involved in tumorigenesis, we performed whole genome microarray analysis of AP-2alpha knockdown cells and observed that AP-2alpha regulates specific genes involved in cell cycle, cell death, adhesion, and migration. In particular, we showed that ESDN, EREG, and CXCL2 play a major role in AP-2 controlled migration, as ablation of any of these genes severely altered migration.

Quantitative expression profiling of highly degraded RNA from formalin-fixed, paraffin-embedded breast tumor biopsies by oligonucleotide microarrays.

Microarray-based gene expression profiling is well suited for parallel quantitative analysis of large numbers of RNAs, but its application to cancer biopsies, particularly formalin-fixed, paraffin-embedded (FFPE) archived tissues, is limited by the poor quality of the RNA recovered. This represents a serious drawback, as FFPE tumor tissue banks are available with clinical and prognostic annotations, which could be exploited for molecular profiling studies, provided that reliable analytical technologies are found. We applied and evaluated here a microarray-based cDNA-mediated annealing, selection, extension and ligation (DASL) assay for analysis of 502 mRNAs in highly degraded total RNA extracted from cultured cells or FFPE breast cancer (MT) biopsies. The study included quantitative and qualitative comparison of data obtained by analysis of the same RNAs with genome-wide oligonucleotide microarrays vs DASL arrays and, by DASL, before and after extensive in vitro RNA fragmentation. The DASL-based expression profiling assay applied to RNA extracted from MCF-7 cells, before or after 24 h stimulation with a mitogenic dose of 17beta-estradiol, consistently allowed to detect hormone-induced gene expression changes following extensive RNA degradation in vitro. Comparable results where obtained with tumor RNA extracted from FFPE MT biopsies (6 to 19 years old). The method proved itself sensitive, reproducible and accurate, when compared to results obtained by microarray analysis of RNA extracted from snap-frozen tissue of the same tumor.

Identification of new genes associated with breast cancer progression by gene expression analysis of predefined sets of neoplastic tissues.

Gene expression profiles were studied by microarray analysis in 2 sets of archival breast cancer tissues from patients with distinct clinical outcome. Seventy-seven differentially expressed genes were identified when comparing 30 cases with relapse and 30 cases without relapse within 72 months from surgery. These genes had a specific ontological distribution and some of them have been linked to breast cancer in previous studies: AIB1, the two keratin genes KRT5 and KRT15, RAF1, WIF1 and MSH6. Seven out of 77 differentially expressed genes were selected and analyzed by qRT-PCR in 127 cases of breast cancer. The expression levels of 6 upregulated genes (CKMT1B, DDX21, PRKDC, PTPN1, SLPI, YWHAE) showed a significant association to both disease-free and overall survival. Multivariate analysis using the significant factors (i.e., estrogen receptor and lymph node status) as covariates confirmed the association with survival. There was no correlation between the expression level of these genes and other clinical parameters. In contrast, SERPINA3, the only downregulated gene examined, was not associated with survival, but correlated with steroid receptor status. An indirect validation of our genes was provided by calculating their association with survival in 3 publicly available microarray datasets. CKMT1B expression was an independent prognostic marker in all 3 datasets, whereas other genes confirmed their association with disease-free survival in at least 1 dataset. This work provides a novel set of genes that could be used as independent prognostic markers and potential drug targets for breast cancer.

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.