Intersection of FOXO- and RUNX1-Mediated Gene PNAS PLUS Expression Programs in Single Breast Epithelial Cells During Morphogenesis and Tumor Progression
Total Page:16
File Type:pdf, Size:1020Kb
Intersection of FOXO- and RUNX1-mediated gene PNAS PLUS expression programs in single breast epithelial cells during morphogenesis and tumor progression Lixin Wanga,1, Joan S. Bruggeb, and Kevin A. Janesa,1,2 aDepartment of Biomedical Engineering, University of Virginia, Charlottesville, VA 22908; and bDepartment of Cell Biology, Harvard Medical School, Boston, MA 02115 Edited by Yoshiaki Ito, Institute of Molecular and Cell Biology, Singapore, and accepted by the Editorial Board July 28, 2011 (received for review March 3, 2011) Gene expression networks are complicated by the assortment of be obtained with conventional oligonucleotide microarrays after regulatory factors that bind DNA and modulate transcription PCR-based amplification of cDNA from 10 microdissected cells combinatorially. Single-cell measurements can reveal biological (2, 22, 23). By surveying hundreds of cells overall, stochastic mechanisms hidden by population averages, but their value has profiling captures the most reproducible cell to cell expression not been fully explored in the context of mRNA regulation. Here, heterogeneities in a cell population. Individual transcripts that we adapted a single-cell expression profiling technique to exam- are nonuniformly expressed are then organized by the pattern ine the gene expression program downstream of Forkhead box O of their expression fluctuations to reveal coordinated single- (FOXO) transcription factors during 3D breast epithelial acinar mor- cell programs. phogenesis. By analyzing patterns of mRNA fluctuations among In this work, we report a detailed analysis centering around individual matrix-attached epithelial cells, we found that a subset one class of transcription factors—the FOXOs—whose expres- of FOXO target genes was jointly regulated by the transcription sion was found to be strongly nonuniform in our initial study factor Runt-related transcription factor 1 (RUNX1). Knockdown of (2). FOXO proteins are a subgroup of the Forkhead family of RUNX1 causes hyperproliferation and abnormal morphogenesis, transcriptional regulators that play important roles in cell cycle both of which require normal FOXO function. Down-regulating arrest, stress responses, and cell death (24). All FOXO isoforms SYSTEMS BIOLOGY RUNX1 and FOXOs simultaneously causes widespread oxidative recognize a common (A/G)TAAA(T/C)A DNA consensus, which stress, which arrests proliferation and restores normal acinar mor- is frequently observed in the extended promoters of many genes phology. In hormone-negative breast cancers lacking human epi- (SI Appendix, Table S1) (25, 26). However, FOXOs are often HER2 fi fi dermal growth factor receptor 2 ( ) ampli cation, we nd that not redundant, and isoform-specific functions have been widely RUNX1 down-regulation is strongly associated with up-regulation documented (27–29). FOXO transcriptional activity is regulated of FOXO1, which may be required to support growth of RUNX1- both positively and negatively by phosphorylation and ubiquiti- negative tumors. The coordinate function of these two tumor sup- nation (30–34). Although FOXOs predominantly function as tran- pressors may provide a failsafe mechanism that inhibits cancer- scriptional activators, they may also act as repressors in certain progression. contexts (35). Last, FOXOs can interact with at least a dozen other transcription factors, resulting in altered binding specificity heterogeneity | triple-negative | stochastic | systems biology and transcriptional activity (SI Appendix, Table S2) (36). Thus, in many ways, FOXO proteins are an archetype for the complexity enetically identical mammalian cells often display patterns that lies between signal transduction and gene expression. Gof gene protein expression that are profoundly different (1– Here, we asked whether a focused examination of endogenous 3). Cell to cell heterogeneity has been recognized in cancer for FOXOs by stochastic profiling would yield insights into their nearly half a century (4, 5). However, only lately have heteroge- function at the network level. Instead of using constitutively neous cell populations been exploited as a means for uncover- active FOXO alleles to homogenize signaling (SI Appendix, ing new mechanisms of biological regulation (6–8). With recently Table S3), we quantified fluctuations in FOXO expression and developed techniques that can interrogate single-cell states, we activity that occur naturally during 3D organotypic culture of are now poised to embrace heterogeneity rather than average it breast epithelial cells (2, 37). By mapping these fluctuations onto out (9–11). a panel of FOXO target genes, we discovered that the measured Nonuniformities emerge at the earliest steps of gene expression FOXO expression signature divides into two groups with distinct (12–14). Therefore, measurements of mRNA expression hetero- single-cell patterns. One of these groups receives a key second geneity may help to uncover new biology if combined with systems input from another transcription factor, RUNX1, which is var- approaches for analysis (11). This possibility offers a particularly iably active. RUNX1 is required for the proper timing of pro- rich opportunity for unraveling transcriptional networks, where liferative suppression in 3D spheroids, and its stable knockdown complex combinations of factors work together to mediate pro- results in abnormal hyperproliferative acinar structures. In- grams of gene expression (15). Early pioneering studies initially terestingly, this phenotype requires native single-cell regulation focused on describing the fundamental kinetics and statistics of transcription in single cells (1, 12–14, 16–19). Now, an open question is how methods that provide single-cell information should be applied to help understand the coordination of tran- Author contributions: L.W., J.S.B., and K.A.J. designed research; L.W. and K.A.J. performed scriptional programs and their role in cell phenotype (20, 21). research; L.W., J.S.B., and K.A.J. analyzed data; and J.S.B. and K.A.J. wrote the paper. For this purpose, we recently developed a technique called The authors declare no conflict of interest. stochastic profiling (2). Stochastic profiling does not attempt to This article is a PNAS Direct Submission. Y.I. is a guest editor invited by the Editorial Board. quantify mRNA levels in single cells but instead, gleans single- 1L.W. and K.A.J. contributed equally to this work. cell information by analyzing the statistical fluctuations in 15–20 2To whom correspondence should be addressed. E-mail: [email protected]. repeated measurements of 10 cells. The increased amount of This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10. starting material allows high-quality, transcriptome-wide data to 1073/pnas.1103423108/-/DCSupplemental. www.pnas.org/cgi/doi/10.1073/pnas.1103423108 PNAS Early Edition | 1of10 Downloaded by guest on September 24, 2021 of FOXO function to withstand the increased oxidative stress expressed was the transcription factor FOXO1. The pattern of caused by RUNX1 knockdown. Combined inhibition of RUNX1 single-cell FOXO1 expression correlated with many transcripts and FOXO leads to an acute state of oxidative stress that induces that had been previously linked to oxidative stress and proliferative proliferation arrest and restores normal 3D morphology. In a suppression, suggesting a coordinated cellular response. Because large study of estrogen receptor (ER)-, progesterone receptor FOXO family transcriptional activity is itself linked to oxidative (PR)-, and HER2-negative breast cancers, we find that reduced stress and proliferation arrest (38, 39), we chose to investigate the RUNX1 expression correlates with FOXO1 up-regulation, which interplay between the expression of FOXOs and FOXO-dependent presumably enables tumor progression. Our results illustrate target genes. how careful single-cell analysis of gene expression can reveal The MCF10A clone 5E (MCF10A-5E) used for 3D culture functional interactions within transcriptional networks that are expresses not only FOXO1 but also FOXO3 (SI Appendix, Fig. important for cancer-relevant cell phenotypes. S1), and we found that both FOXO proteins were heteroge- neously expressed and localized in individual acini during mor- Results phogenesis (Fig. 1A). Together with the two FOXOs, we mon- Dissecting Single-Cell FOXO Coregulation by Stochastic Sampling. We itored time-dependent changes in the levels of 18 transcripts previously cataloged the transcriptional heterogeneities that during morphogenesis: eight genes were validated or reported emerge among individual matrix-attached breast epithelial cells FOXO targets (BTG1, CAV1, CDKN1A, FBXO32, SEMA3C, during acinar morphogenesis in a 3D culture model (2, 37). One SESN1, SOD2, and SOX4)(SI Appendix, Table S3), five genes gene that was predicted with high confidence to be nonuniformly were constitutively expressed (GAPDH, HINT1, PRDX6, S100A6, AB) FOXO1 FOXO3 BTG1 CDKN1A FBXO32 SOD2 SESN1 CAV1 SEMA3C SOX4 10 FOXO1 FOXO3 1 0 -1 Validated FOXO target d induction (log Reported FOXO target Fol 4 6 8 1012 4 6 8 1012 4 6 8 1012 4 6 8 1012 4 6 8 1012 4 6 8 1012 4 6 8 1012 4 6 8 1012 4 6 8 1012 4 6 8 1012 E-cadherin DAPI ) PRDX6 HDPAG HINT1 S100A6 UBC CDKN1C KRT10 1BNCC BCL2L13 CCNI 10 1 0 -1 Loading controls d induction (log Other Fol 4 6 8 1012 4 6 8 1012 4 6 8 1012 4 6 8 1012 4 6 8 1012 4 6 8 1012 4 6 8 1012 4 6 8 1012 4 6 8 1012 4 6 8 1012 Day of morphogenesis C D p E Log10 (homogeneously expressed) 0 -4 -8 -12 -16 CCNB1 CCNB1 3 KRT10 3 SOD2 Candidate SOD2 BCL2L13 heterogeneities KRT10 2 2 CCNI FBXO32 SOX4 SESN1 FOXO1 1 1 CDKN1C FOXO3 BTG1 CDKN1A CCNI 0 0 FOXO3 BCL2L13 FOXO1 CDKN1A -1 CAV1 -1 FBXO32 CAV1 UBC change from geometric mean SOX4 2 -2 GAPDH -2 SEMA3C BTG1 Log False-discovery rate PRDX6 CDKN1C -3 HINT1 Standard deviations from geometric mean -3 S100A6 SESN1 61 11 9 3412 14 7813 5618 71 12 1015 123456789101112131415161718 10-cell stochastic samplings 10-cell stochastic samplings Fig. 1. Focused stochastic sampling of a FOXO expression dichotomy in 3D breast epithelial cultures. (A) Heterogeneous expression and localization of FOXO1 and FOXO3 proteins at day 10 of morphogenesis.