Abstract

Estrogen regulates diverse physiological effects and drives breast tumor progression by directly activating estrogen receptor (ER). However, due to the stochastic nature of gene transcription and the resulting heterogeneous cellular response, it is important to investigate estrogen-stimulated gene expression profiles at the single-cell level in order to fully understand how ER regulates transcription in breast cancer cells. In this study, we performed single-cell transcriptome analysis on ER-positive breast cancer cell lines following 17{beta}-estradiol stimulation. Overall, we observed robust gene expression diversity between individual cells. Moreover, we found over two thirds of the genes in breast cancer cells displayed a bimodal expression pattern, which caused averaging artifacts and masked the identification of potential estrogen-regulated genes. We overcame this issue by reconstructing a dynamic estrogen-responsive transcriptional network from discrete time points into a pseudotemporal continuum. Pathway analysis of the differentially expressed genes derived from the pseudotemporal analysis showed an estrogen-stimulated metabolic switch that favored biosynthesis and cell proliferation but reduced estrogen degradation. In addition, we identified folate-mediated one-carbon metabolism as a novel estrogen-regulated pathway in breast cancer cells. Notably, estrogen stimulation reprogramed this pathway through the mitochondrial folate pathway to coordinately fuel polyamine and de novo purine synthesis. Finally, we showed AZIN1 and PPAT, key regulators in the above pathways, are direct ER target genes and essential for breast cancer cell survival and growth. In summary, our single-cell study illustrated a dynamic transcriptional heterogeneity in ER-positive breast cancer cells in response to estrogen stimulation and uncovered a novel mechanism of an estrogen-mediated metabolic switch.