An integrative transcriptomics approach identifies miR-503 as a candidate master regulator of the estrogen response in MCF-7 breast cancer cells

Estrogen receptor α (ERα) is an important biomarker of breast cancer severity and a common therapeutic target. In response to estrogen, ERα stimulates a dynamic transcriptional program including both coding and noncoding RNAs. We generate a fine-scale map of expression dynamics by performing a temporal profiling of both messenger RNAs (mRNAs) and microRNAs (miRNAs) in MCF-7 cells (an ER+ model cell line for breast cancer) in response to estrogen stimulation. We identified three primary expression trends—transient, induced, and repressed—that were each enriched for genes with distinct cellular functions. Integrative analysis of mRNA and miRNA temporal expression profiles identified miR-503 as the strongest candidate master regulator of the estrogen response, in part through suppression of ZNF217—an oncogene that is frequently amplified in cancer. We confirmed experimentally that miR-503 directly targets ZNF217 and that overexpression of miR-503 suppresses MCF-7 cell proliferation. Moreover, the levels of ZNF217 and miR-503 are associated with opposite outcomes in breast cancer patient cohorts, with high expression of ZNF217 associated with poor survival and high expression of miR-503 associated with improved survival. Overall, these data indicate that miR-503 acts as a potent estrogen-induced candidate tumor suppressor miRNA that opposes cellular proliferation and has promise as a novel therapeutic for breast cancer. More generally, our work provides a systems-level framework for identifying functional interactions that shape the temporal dynamics of gene expression.