Single-Molecule Sequencing Reveals Estrogen-Regulated Clinically Relevant lncRNAs in Breast Cancer

Estrogen receptor (ER)α-positive tumors are commonly treated with ERα antagonists or inhibitors of estrogen synthesis, but most tumors develop resistance, and we need to better understand the pathways that underlie the proliferative and tumorigenic role of this estrogen-activated transcription factor. We here present the first single-molecule sequencing of the estradiol-induced ERα transcriptome in the luminal A-type human breast cancer cell lines MCF7 and T47D. Sequencing libraries were prepared from the polyadenylated RNA fraction after 8 hours of estrogen or vehicle treatment. Single-molecule sequencing was carried out in biological and technical replicates and differentially expressed genes were defined and analyzed for enriched processes. Correlation analysis with clinical expression and survival were performed, and follow-up experiments carried out using time series, chromatin immunoprecipitation and quantitative real-time PCR. We uncovered that ERα in addition to regulating approximately 2000 protein-coding genes, also regulated up to 1000 long noncoding RNAs (lncRNAs). Most of these were up-regulated, and 178 lncRNAs were regulated in both cell lines. We demonstrate that Long Intergenic Non-protein Coding RNA 1016 (LINC01016) and LINC00160 are direct transcriptional targets of ERα, correlate with ERα expression in clinical samples, and show prognostic significance in relation to breast cancer survival. We show that silencing of LINC00160 results in reduced proliferation, demonstrating that lncRNA expression have functional consequences. Our findings suggest that ERα regulation of lncRNAs is clinically relevant and that their functions and potential use as biomarkers for endocrine response are important to explore.