SOX11 promotes epithelial/mesenchymal hybrid state and alters tropism of invasive breast cancer cells
Abstract
SOX11 is an embryonic mammary epithelial marker that is normally silenced prior to birth. High SOX11 levels in breast tumours are significantly associated with distant metastasis and poor outcome in breast cancer patients. Here, we show that SOX11 confers distinct features to ER-negative DCIS.com breast cancer cells, leading to populations enriched with highly plastic hybrid epithelial/mesenchymal cells, which display invasive features and alterations in metastatic tropism when xenografted into mice. We found that SOX11+DCIS tumour cells metastasize to brain and bone at greater frequency and to lungs at lower frequency compared to cells with lower SOX11 levels. High levels of SOX11 leads to the expression of markers associated with mesenchymal state and embryonic cellular phenotypes. Our results suggest that SOX11 may be a potential biomarker for breast tumours with elevated risk of developing metastases and may require more aggressive therapies.
Introduction
SOX11 is an embryonic mammary factor that is not expressed in normal breast epithelial cells after birth (Wansbury et al., 2011). However, SOX11 is expressed in many triple negative and HER2+ invasive breast cancers (Wansbury et al., 2011). SOX11 expression in invasive breast cancer is associated with increased distant metastasis formation (Oliemuller et al., 2017). Inhibition of SOX11 by siRNA suppressed growth and proliferation of ER- breast cancer cell lines, but had no significant effect on growth and proliferation of ER+ breast cancer cell lines (Shepherd et al., 2016). SOX11 repression using siRNA reduced both cell migration and invasion in basal-like breast cancer (BLBC) cell lines, supporting a role for SOX11 in promoting breast cancer progression. In addition, SOX11 inhibition in MDA-MB-468, a BLBC line, resulted in reduced expression of FOXC1, CCNE1, KRT14, MIA and SFRP1, all of which are PAM50 genes that are highly expressed in BLBC. This finding suggests SOX11 may modulate key features of basal-like cancer cells, including Keratin 14 expression, a marker expressed by basal cells and some luminal cells within the terminal ductal lobular unit, where many breast cancers arise (Gusterson and Eaves, 2018).
Our previous studies showed that when constitutively expressed by the human breast epithelial cell line MCF10A, SOX11 increased the number of basal/myoepithelial clones formed in standard colony-formation assays and led to increased mammosphere formation, suggesting that SOX11 can modulate mammary progenitor features and stem cell activity when expressed in mature postnatal mammary epithelial cells (Oliemuller et al., 2017). Furthermore, we also showed that SOX11 promotes invasive transition of DCIS.com cells that were injected into the mammary duct to mimic formation of DCIS-like lesions within the mammary duct prior to progression to invasive state and tumour formation (Oliemuller et al., 2017). Expression of SOX11 in preinvasive breast lesions and potential sites of microinvasion in samples from DCIS cases supports a role for SOX11 in promoting in situ to invasive breast carcinoma transition. These findings led us to hypothesise that reactivation of embryonic mammary developmental programmes mediated by SOX11 in postnatal breast epithelial cells would promote invasive progression of pre-invasive DCIS breast lesions and acquisition of features associated with poor patient outcome, including formation of distant metastasis.
Therefore, to further explore the role of SOX11 in regulating cancer stem cell (CSC) states and to determine whether SOX11 promotes metastatic dissemination of invasive breast cancer cells, we have created an inducible model to study the roles of SOX11 in progenitor/stem cell regulation and breast cancer progression in vitro and in vivo. In the model described here, SOX11 is expressed at higher levels in DCIS.com cells driven from the EF1A promoter after induction with doxycycline (DOX), when compared to another model we have previously used to study DCIS progression, in which low levels of SOX11 expression is constitutively driven by the CMV promoter (Oliemuller et al., 2017). In the inducible model presented here, SOX11 expression confers epithelial/mesenchymal hybrid state, features typical of embryonic mammary cell phenotypes, and influences organ-specific metastatic tropism to breast cancer cells. Finally, we also identify several novel SOX11 targets, many of which are highly correlated with SOX11 expression in primary breast cancers and breast cancer metastases.
Results
Inducible expression of SOX11 leads to changes in stem cell profiles of DCIS.com cells
To investigate the role of SOX11 in breast cancer progression, we used the pINDUCER21 system to stably transduce DCIS.com cells, an invasive cell line from the MCF10A breast cancer progression series, so that SOX11 was expressed only when induced with Doxycycline (DOX) (referred to as iSOX11 cells) (Figure 1A–B). The results show a significantly higher, sustained expression of SOX11 levels compared with the previous constitutive model we have used to study DCIS progression which lost SOX11 expression over time (Figure 1—figure supplement 1; Oliemuller et al., 2017). As expected, SOX11 localised mostly to the nuclei in iSOX11 cells, similar to that observed in SOX11+ DCIS case samples (Figure 1A–C and Figure 1—figure supplement 1). SOX11 is also detected in the cytoplasm of iSOX11 cells using western blotting (Figure 1A), a location that was not observed in the DCIS-SOX11 cells (data not shown), showing that some differences exist when SOX11 is expressed at different levels in the two models.
To study the role of SOX11 in regulating stem cell state, we used the inducible iSOX11 model and assessed CSC profiles. CSCs are subpopulations of cancer cells sharing similar characteristics as normal stem or progenitor cells such as self-renewal ability and multi-lineage differentiation to drive tumour growth and heterogeneity. ALDH1 and CD24 are widely used CSC markers in breast cancer (Liu et al., 2014). A higher proportion of CD24+ and ALDH+ cells are detected in iSOX11 cells compared to control iEV cells (DCIS.com cells stably transduced with pINDUCER21 empty vector and induced with DOX) by flow cytometry (Figure 1D–E and Figure 1—figure supplements 1–2). Moreover, CD24 expression levels were significantly increased in a DOX dose-dependent manner in iSOX11 cells (Figure 1—figure supplement 1). CD24 is predominantly located in the membrane and CD24 was detected in the cytoplasmic fraction of both iEV and iSOX11 cells, and also detected in the nuclear fraction of iSOX11 cells, which has been previously reported in breast cancer (Figure 1F; Duex et al., 2017).
When cells were grown on collagen, we detected loss of membranous E-Cadherin localisation and acquisition of N-Cadherin expression in iSOX11 cells and formation of bud-like structures (Figure 1G–H). Nuclear expression of N-Cadherin has been reported previously, primarily in poorly differentiated breast tumours (Rezaei et al., 2012). We also detected a significant increase in iSOX11 cells expressing Vimentin by western blotting; increased Vimentin was detected in both the cytoplasmic and nuclear fraction. Vimentin was not observed in the nucleus of cells grown in the same conditions when stained by immunofluorescence, suggesting that the nuclear fraction contains proteins attached to the external part of the nuclear envelope since both Lamin B and Vimentin are known to associate with the nuclear envelope (Figure 1G–H and Figure 1—figure supplement 3; Georgatos and Blobel, 1987a; Georgatos and Blobel, 1987b). We assessed the expression of markers associated with mammary epithelial and mesenchymal states in iEV and iSOX11 as well as non-induced (denoted ni) niEV and niSOX11 DCIS cells. K5, K8, and VIM levels are elevated in iSOX11 cells compared to iEV cells (Figure 1—figure supplement 3). These results are consistent with a role for SOX11 in promoting a mesenchymal state and influencing mammary epithelial phenotypes when expressed in epithelial cells.
To further analyse the effects of SOX11 on mammary epithelial phenotypes, we assessed the expression of markers associated with the two major mammary lineages in DCIS cells grown as spheroids that formed from iSOX11 and to iEV, niEV, and niSOX11 control cells. We detected significant increases in expression of both luminal (K8) and basal (K14, SMA) lineage markers in SOX11-expressing spheroids as well as a significant increase in cells co-expressing markers of both lineages, which is suggestive of embryonic mammary phenotypes (Figure 1I and Figure 1—figure supplements 4–5).
DCIS cells expressing SOX11 grow more slowly and form more invasive spheroids
Reduced cell growth was detected in iSOX11 cells compared to control iEV cells grown in both 2D or as spheroids (Figure 2A–B). Control iEV cells were more clonogenic than iSOX11 cells whether plated in colony-forming assays or when grown from single cells (Figure 2—figure supplement 1). Although iEV cells have greater colony-forming capacity than iSOX11 when large colonies (greater than 50 microns) were counted, a greater number of small colonies (less than 50 microns) formed (Figure 2—figure supplement 1). We found DCIS-LacZ control cells also exhibited both greater colony and sphere-forming capacity than DCIS-SOX11 cells (Figure 2—figure supplement 1). As with iSOX11 cells, a greater number of small colonies form in colony-formation assays with DCIS-SOX11 cells. Spheroids formed from PKH-labelled DCIS-SOX11 cells are smaller and retain PKH dye more than labelled DCIS-LacZ cells (Figure 2—figure supplement 1). These results show that SOX11 expression leads to a reduction of larger spheroids formed from DCIS.com cells and increased formation of smaller spheroids that retain label. Multiple attempts to transduce primary human breast epithelial cells to express SOX11 constitutively with a CMV-driven vector were unsuccessful but are consistent with a possible role for SOX11 in regulating a quiescent state (data not shown).
SOX11 expression is enriched in ER- and HER2+ invasive breast cancers (Figure 2—figure supplement 2). To explore whether sustained reactivation of SOX11 promotes tumour progression, we assessed features of iSOX11 spheroids. Spheroids formed from iSOX11 DCIS cells were smaller and formed a higher number of peripheral microspheres than iEV cells (Figure 2C). When overlaid with Collagen I, spheroids formed from iSOX11 cells are more invasive compared to iEV, niEV, and niSOX11 control cells (Figure 2D). Tumours that formed after mammary fat pad xenografts of iSOX11 cells were smaller than iEV tumours when mice were fed chow with moderate levels (200 or 625 ppm) of DOX (Figure 2E–F). Despite their smaller volume, tumours originated from cells with high SOX11 levels displayed greater bioluminescence than control tumours, suggesting iSOX11 tumours contained more viable cells (Figure 2G). IHC staining of mammary tumours was performed to detect Cleaved Caspase 3 (CC3). Larger necrotic and CC3+ (apoptotic) areas were observed in the EV tumour tumours and niSOX11 tumours when compared to the iSOX11 tumours, which showed little central necrosis and fewer CC3+ cells (Figure 2H). We also observed that extremely small and less luminescent tumours were formed in mice xenografted with iSOX11 cells, compared to control cells that were fed higher levels (1250 or 2000 ppm) of DOX-chow (Figure 2—figure supplement 3). Higher levels of CC3 were detected in these iEV tumours, which displayed central necrosis, whilst almost no CC3 or necrosis was detected in iSOX11 tumours (Figure 2—figure supplement 3). Tumours formed from iSOX11 DCIS cells grew out quickly when (1250 or 2000 ppm) DOX-chow was replaced with normal chow, suggesting that high levels of SOX11 could keep tumours in a non-proliferative state and that upon DOX withdrawal, proliferation resumed and the ostensible quiescent state is a reversible condition (Figure 2—figure supplement 3).
When injected directly into the mammary duct, iSOX11 DCIS cells formed slightly larger tumours, with similar bioluminescence levels as iEV tumours, which indicates the microenvironment highly influences behaviour of SOX11+ tumour cells (Figure 2I–K). A greater proportion of mammary tumours formed from iSOX11 cells expressed moderate to high levels of ALDH1 (6/6) compared to control tumours (1/8). In addition, CD24+ cells were observed at greater frequency in iSOX11 tumours compared to control tumours (Figure 2—figure supplement 4). Of note, we observed nuclear CD24 staining in iSOX11 tumour cells, but not in iEV, niEV, or niSOX11 tumours, in line with the observed in vitro result (Figure 2—figure supplement 4). It has been suggested that cells designated CD24- using flow cytometry maybe expressing CD24 in the nucleus and this promotes aggressive tumour properties, since it has been shown that nuclear CD24 is able to drive tumour growth (Duex et al., 2017).
SOX11 expression promotes expression of developmental pathways frequently activated in cancer
To identify potential SOX11 targets that could confer features of embryonic, epithelial/mesenchymal hybrid state, and metastasis-promoting features to mature breast cells, RNA from iEV, niEV, or niSOX11 control and iSOX11 cells grown in 2D and in 3D (from spheroids at two time-points: 2 days and 5 days after spheroid formation) were sequenced (Figure 3A, Supplementary file 1 and Figure 3—figure supplement 1). Uninduced niSOX11 control cells express very low levels of SOX11, and suggest the vector is slightly leaky. iSOX11 DCIS cells cultured in 2D showed enrichment of genes regulating actin filament sequestration, phospholipid catabolism, ERBB signaling, chemotaxis, and epithelial differentiation (Figure 3B and Supplementary file 2). iSOX11 DCIS cells grown in 3D showed enrichment of genes regulating ECM disassembly, collagen biosynthesis, glycosaminoglycan metabolism, and platelet degranulation (Figure 3—figure supplement 1 and Supplementary file 3). Platelet activation is one of the first steps of tissue repair as part of the wound healing process and Sox4 and Sox11 have recently been shown to reactivate an embryonic epidermal programme during wound repair in mice (Miao et al., 2019). We found substantial overlap of the embryonic wound signature that was shown to be directly regulated by Sox11 and Sox4 (Miao et al., 2019) with iSOX11 spheroids (Table 1). In particular, we detected upregulation of embryonic wound signature components with links to actin polymerisation and cell adhesion, and one known regulator of embryonic stem cell pluripotency, RCOR2, which can function with other transcription factors to induce pluripotent stem cells (Figure 3—figure supplement 1; Yang et al., 2011).
RNA sequencing analysis detected CDH2 (encoding N-Cadherin) expressed at significantly higher levels in iSOX11 cells grown in 2D and as spheroids, whilst other EMT markers, including CDH1 (encoding E-cadherin) and VIM were not significantly changed compared to control cells (Supplementary file 1). Other notable downstream targets of SOX11 included MEX3A, which encodes an RNA-binding protein, that marks slowly proliferating multipotent stem cells in mouse intestine (Barriga et al., 2017) and totipotent cells in C. elegans (Pereira et al., 2013); MMP11, ST6GALNAC5, and TUBB3, which are highly expressed in breast cancers that metastasize to brain (Bos et al., 2009; Kim et al., 2015; Lee et al., 2016). We confirmed that a number of putative SOX11 targets of interest (Figure 3C and Figure 3—figure supplement 1, Supplementary file 2), including MEX3A and TUBB3 were upregulated in iSOX11 cells when measured by qPCR (Figure 3D). We also detected upregulation of MEX3A and TUBB3 protein in iSOX11 cells (Figure 3E–F). We found many SOX11+ breast cancer cell lines express high levels of MEX3A or TUBB3 compared to DCIS.com cell line (Figure 3G). MEX3A and TUBB3 levels are correlated with SOX11 expression in breast cancers in the TCGA dataset (Figure 3H) and in breast cancer cell lines in the Broad dataset (Ghandi et al., 2019; Figure 3—figure supplement 1). Notably, with greater increases of SOX11 levels, a higher percentage of samples with increased MEX3A are observed. These findings support a role for SOX11 in mediating developmental signals during breast cancer progression.
DCIS cells expressing SOX11 show alterations in metastatic tropism
To explore whether sustained reactivation of SOX11 promotes tumour progression, we injected luciferase-tagged iEV and iSOX11 cells into the mammary fat pad. Four weeks after orthotopic xenografting with iSOX11 cells, brain micrometastases were detected in two out of six mice by IVIS imaging, whilst none were observed in six mice xenografted with control cells (Figure 4A). Liver and lung micrometastases were detected in both cohorts when assessed by IVIS imaging.
After injection into the tail vein, IVIS imaging detected tumour cells in lungs in seven out of the eight mice, bone in three out of eight mice, and in brain in one of eight mice xenografted with iSOX11 cells, whilst no bone or brain metastases were observed in eight mice xenografted with control cells (Figure 4B). It was noted that the frequencies of mice with lung metastasis were similar, but a significant reduction of iSOX11 DCIS cells accumulated in the lungs was observed when compared to mice engrafted with iEV cells when quantified by IVIS (Figure 4C–D). Macroscopic examination of metastatic lesions confirmed the reduction in tumour burden (data not shown).
SOX11 expression in primary breast cancer is associated with increased metastasis formation at distant sites (Figure 4E). SOX11 is highly expressed in brain metastases and is also detected in bone metastasis from breast cancer patients (Figure 4—figure supplement 1). SOX11 amplification and overexpression has been detected in breast cancer brain metastasis in patients with ER-, ER+, and BRCA1-/- tumours in a recent study that comprehensively profiled a small number of cases (Figure 4—figure supplement 1; Saunus et al., 2015). In another dataset in which samples from 21 breast cancer brain metastasis (BCBM) patients were transcriptionally profiled by RNA sequencing, elevated levels of SOX11 are detected in the brain metastasis in one third of cases compared to the primary tumour (Figure 4—figure supplement 2; Varešlija et al., 2019).
SOX11+ DCIS.com cells from brain metastasis display a colonisation and growth advantage after intracranial xenografting iSOX11 cells were isolated from mouse brain after orthotopic injections and expanded (designated iSOX11Br). Induction of SOX11 with DOX led to increased expression of SOX11, MEX3A and TUBB3 compared to untreated cells (Figure 5A). Other cells isolated and expanded from bone and lung metastasis after orthotopic injections of iSOX11 did not express SOX11 after induction with DOX.
After xenografting iSOX11Br cells into the tail vein, mice that had been fed normal chow had a greater metastatic burden in the lungs compared to mice fed DOX-containing chow to induce SOX11 expression, similar to results observed using the parental iSOX11 DCIS cell line (Figure 5B–D). After intracranial injections of iSOX11Br cells, higher levels of bioluminescence were detected in the brains of mice fed DOX chow (Figure 5E). Induction of SOX11 expression led to larger tumour burden in the brain and reduced survival (Figure 5F). These results indicate that iSOX11Br cells have a colonisation and growth advantage in the brain compared to that niSOX11Br cells lacking SOX11 expression.
SOX11 regulates proliferative state of ER- breast cancer cells
Next, we examined the effect of reducing SOX11 levels in CAL-148 ER- breast cancer cell line that expresses very high SOX11 levels. Using DOX-inducible lentiviral vectors containing shRNAs to SOX11 and control non-targeting shRNAs, reductions in SOX11 levels were obtained, as well as decreases of TUBB3, MEX3A, GPC2, MPK4, OLFM2, ST6GALNAC5, and NCAD levels and an increase in SERPINA3 in CAL-148 cells after SOX11 knockdown when compared to control (Figure 6A). Reduced levels of SOX11, MEX3A and TUBB3 protein were detected after SOX11 knockdown (Figure 6B). Cell viability assays detected greater cell numbers after SOX11 knockdown in CAL-148 cells, compared to control cells (Figure 6C). Colony-formation assays detected an increase in clonogenetic potential of CAL-148 cells with reduced SOX11 levels (Figure 6D). These results suggest SOX11 could regulate proliferative state of stem cells in ER- breast cancer cells.
TUBB3, an established SOX11 target, regulates proliferation and invasion of ER- breast cancer cells
As with SOX11, both distant metastasis-free survival and overall survival of breast cancer patients are reduced when high levels of TUBB3, an established SOX11 target in neural cells, are expressed in primary tumours (Figure 7—figure supplement 1). We next examined the effects of TUBB3 on triple negative breast cancer (TNBC) growth and invasion. Using siRNA-mediated knockdown, we found that reducing TUBB3 levels in BT-20 cells resulted in reduced growth in both 2D culture and spheroid culture (Figure 7A–C). More cells were arrested in G2/M phase of the cell cycle when TUBB3 levels were reduced (Figure 7D and Figure 7—figure supplement 2), consistent with a known role of TUBB3 in regulating cell cycle progression of tumour cells. As a result of cell cycle arrest, an increase of multinucleated cells in the >G2/M phase and a higher proportion of dead cells in sub-G1 phase were detected. Spheroid invasion assay s detected less invasion when TUBB3 levels were reduced (Figure 7E–F). These results suggest SOX11 can regulate proliferation and invasive growth through TUBB3 in ER- breast cancer cells.
MEX3A, a novel potential SOX11 downstream effector, regulates cell growth and E/M state of ER- breast cancer cells
Both distant metastasis-free survival and overall survival of breast cancer patients are reduced when high levels of MEX3A are expressed in primary tumours, as would be expected for a SOX11 target (Figure 7—figure supplement 1). Due to MEX3A’s established links with regulation of both EMT and proliferation of intestinal and various types of cancer cells, we knocked down MEX3A in ER- BT-20, CAL-148 and HCC1187 breast cancer cells to determine if MEX3A regulates similar processes. In HCC1187 and BT-20 cells, a slight but significant reduction of invasive growth was detected when MEX3A levels were reduced with some, but not all siRNAs (Figure 8A–B and Figure 8—figure supplement 1). An increase in cell numbers were observed after knockdown of MEX3A in both CAL-148 and HCC1187 cells grown in 2D or 3D (Figure 8C–D). Cell cycle analysis detected a reduction in S phase after MEX3A knockdown in both cell lines (Figure 8E, Figure 8—figure supplement 1). Several candidate cell cycle regulators from a consensus stem cell quiescence signature (Cheung and Rando, 2013) were downregulated in iSOX11 cell signatures, including RRM2 and SURVIVIN (Figure 8—figure supplement 1 and Supplementary file 3). Western blotting found that both RRM2 and SURVIVIN were upregulated when MEX3A levels were reduced in CAL-148 and BT-20 cells (Figure 8F and Figure 8—figure supplement 1). High expression levels of MEX3A co-occur with increased levels of E2F3, CCNE1 and CDKN2A, and decreases in RB1 levels in the TCGA dataset (Figure 8—figure supplement 1); MEX3A levels show strong correlation with the levels of E2F3 and CCNE1 (Supplementary file 4). The addition of EGF to CAL-148 cells growing in either serum-free or low-serum media led to reduction of MEX3A levels (Figure 8—figure supplement 2). CAL-148 cells normally form aggregates when grown in 2D culture conditions and after MEX3A knockdown, CAL-148 cells displayed a reduced ability to form aggregates and an acquired ability to adhere to plastic (Figure 8G). After MEX3A knockdown, CAL-148 cells displayed increased expression of E-CADHERIN and EPCAM (Figure 8H). BT-20 and HCC1187 cells lack expression of E-CADHERIN but showed an increase in EPCAM levels after MEX3A knockdown (Figure 8—figure supplement 2). Together, these results are consistent with roles for MEX3A in regulation of cell growth and EMT in SOX11+ ER- breast cancer cells.
Discussion
We previously showed that SOX11 promotes invasive transition of DCIS cells using both in vitro and in vivo models where SOX11 is expressed in DCIS.com cells under the control of a constitutive CMV-driven promoter (Oliemuller et al., 2017). Here, we investigated the role of SOX11 in breast cancer progression after tumour formation using a doxycycline-inducible EF1A promoter to express SOX11 in DCIS.com cells at a level comparable to that observed in clinical DCIS and breast cancer samples. DCIS cells induced to express SOX11 prior to spheroid formation form smaller spheroids, which display more invasion compared to control spheroids, recapitulating the phenotypes we observed with constitutive expression of SOX11 (Oliemuller et al., 2017). In addition, we observed unique phenotypes using these experimental conditions, including substantial morphological changes (cell detachment and formation of multiple satellite spheroids indicative of a highly invasive phenotype without hydrogel), suggesting breast lesions expressing very high levels of SOX11 possess an inherent high potential to form metastatic lesions.
Stem cells and progenitors may serve as a cell of origin for breast cancers. Many SOX (SRY-related HMG-box) transcription factors are expressed in the postnatal breast and some have been shown to control normal and/or CSCs (Domenici et al., 2019; Kogata et al., 2018; Mehta et al., 2019). SOX11 is a unique transcription factor, since it does not appear to be expressed in normal mammary epithelial cells after birth in either mouse or human and therefore is well-poised to reactivate developmental pathways when expressed in breast cancers (Oliemuller et al., 2017; Tsang et al., 2020; Wansbury et al., 2011; Zvelebil et al., 2013).
CD24 and ALDH1 are widely used CSC markers in breast cancer (Liu et al., 2014). It is now widely accepted that stem cell states fluctuate and are not fixed (Liu et al., 2014) and cells transition between the two states (Liu et al., 2014; O'Brien-Ball and Biddle, 2017). Our results suggest a dose-dependent effect of SOX11 on CD24 levels. It is possible that CMV-driven versus EF1A-driven SOX11 has distinct cell context effects. Lower levels of CMV-driven SOX11 expression results in expansion of the ALDH+/CD24- CSC population in DCIS.com cells (Oliemuller et al., 2017). Meanwhile, EF1A-driven inducible SOX11 expression, presented here, leads to expansion of the CD24+ CSC population and higher levels of SOX11 are associated with an increased size of the ALDH+/CD24+ CSC population. Triple negative breast cancer (TNBC) cells exhibit robust expression of CD24, suggesting that the inducible model described here mimics what is observed in human TNBC and may present a potential therapeutic opportunity for some SOX11+ breast cancers (Barkal et al., 2019).
Sox11 is highly expressed in prenatal mammary epithelial cells from the time the mammary organ initially forms (Zvelebil et al., 2013). scRNA-sequencing also detects Sox11 expression by the majority of Lgr5+ embryonic mammary epithelial cells isolated at E14.5, a stage when embryonic mammary epithelial cells are multipotent and most cells express markers associated with the two major mammary lineages (the basal/myoepithelial and luminal lineages) (Lilja et al., 2018; Wuidart et al., 2018). Prenatal human breast tissues also express markers associated with both basal and luminal mammary lineages (Jolicoeur et al., 2003). Co-expression of basal and luminal markers are observed in iSOX11 cells and spheroids formed from them. Our findings are consistent with SOX11 marking populations enriched with embryonic mammary cell phenotypes (Bland and Howard, 2018; Wuidart et al., 2018).
Epithelial cells with mesenchymal features have been detected at high frequency during organogenesis in the developing mouse embryo using single cell RNA-seq analysis. Sox11 positively regulates expression of mesenchymal markers, including N-cadherin (Cdh2) and Fibronectin1 (Fn1) in several developing organs, including intestine, liver, lung, and skin (Dong et al., 2018; Halbleib and Nelson, 2006). A significant association of SOX11 and N-CADHERIN (CDH2) expression is also observed in genomic analyses across six cancer types, suggesting that mesenchymal pathways may be activated by SOX11 in both normal and cancer cells (Vervoort et al., 2013). SOX11 has previously been implicated in participating in the regulation of epithelial-mesenchymal transition (EMT) (Venkov et al., 2011). EMT, stemness and plasticity are intertwined (Nieto et al., 2016; Wahl and Spike, 2017) and SOX11 is well poised to function as a key regulator of epithelial/mesenchymal cell states during development and in cancers.
The regulation of N-CADHERIN (CDH2) by SOX11 is significant since CDH2 promotes motility in human breast cancer cells regardless of their E-CADHERIN expression (Nieman et al., 1999). CDH2 is a predictive biomarker for distant metastasis in early-stage breast cancer (Aleskandarany et al., 2015), that is commonly detected in breast cancer cells and provides a mechanism for transendothelial migration (Qi et al., 2005).
Micrometastasis to brain is detected in ~30% of mice within 1 month of xenografting DCIS cells induced to express high levels of SOX11. SOX11 is amplified and overexpressed in ~30% of brain metastases in a recent study of a small cohort of breast cancer patients that were profiled using an integrated genomic and transcriptomic analysis of fresh frozen tumour samples (Saunus et al., 2017; Saunus et al., 2015) and highly expressed by ~30% of BCBM in another independent study (Varešlija et al., 2019). Our finding that iSOX11 mammary tumours spontaneously metastasize to the brain is clinically relevant.
We have identified a number of genes that are regulated by SOX11 in our mouse models of breast cancer that are also associated with SOX11 expression in both primary breast cancer and breast cancer brain metastasis. A number of these, such as TUBB3, encode targetable molecules that are being used in novel pharmacogenetic approaches in combination with other factors (Karki et al., 2013). Others, such as MEX3A, may provide novel biomarkers or therapeutic targets (Bufalieri et al., 2020; Wang et al., 2020; Yang et al., 2020). MEX3A controls the polarity and stemness and affects the cell cycle of intestinal epithelial cells through the downregulation of the mRNA encoding the CDX2 transcription factor (Pereira et al., 2013) and, in addition, exhibits a transforming activity when overexpressed in gastric epithelial cells (Jiang et al., 2012). The intestine has a reservoir of quiescent stem cells that are resistant to chemotherapy that are marked by Mex3a and are multipotent so have the capacity to produce any kind of intestinal cell and contribute to tumour heterogeneity (Barriga et al., 2017). MEX3A may mark a rare mammary stem cell in the human breast that could escape traditional chemotherapy treatments, but this remains to be demonstrated.
The analysis of co-occurrence of SOX11 with MEX3A and TUBB3 in breast cancers profiled in TCGA suggests that MEX3A is expressed predominantly in samples with higher levels of SOX11. The differences in cell growth observed both in vitro and in vivo with increasing levels of SOX11 suggest a dose-dependent regulation of SOX11 downstream targets and/or a possible post-translational modification. Also, the discrepancies between the tumour size observed after xenografting iSOX11 cells via the mammary fat pad model and the MIND model, as well as the finding that EGF can decrease MEX3A levels, suggest an influence of microenvironmental factors on SOX11 tumour cell behaviour. This is in agreement with differential expression of TUBB3 and MEX3A, depending if DCIS-SOX11 cells were xenografted in the mammary fat pad or via MIND (Oliemuller et al., 2017). Together, this data suggests it may be necessary to classify SOX11+ tumours, depending on SOX11 expression level, as well as the expression of its effectors in order to stratify breast cancer patients. The subcellular localisation of SOX11 can be informative for cancer classification since high SOX11 mRNA levels and detection of the nuclear protein are reliable markers of mantle cell lymphoma (Mozos et al., 2009). During neurogenesis, SOX11 is detected in both the nucleus and cytoplasm (Balta et al., 2018). Phosphorylation of specific serine residues can modulate SOX11 subcellular localisation and prevent its nuclear localisation (Balta et al., 2018). We have not determined whether SOX11 localises to the cytoplasm, as well as the nucleus, in breast cancer case samples, and if SOX11 localisation could be useful for classification of breast cancers.
SOX11+ breast cancer cells express markers indicative of phenotypic plasticity and have a high tendency to undergo metastasis (O'Brien-Ball and Biddle, 2017). Together, these data suggest that patients whose DCIS and primary breast cancers express high levels of SOX11 are among a high-risk metastasis subgroup that should be considered for aggressive therapies in neo-adjuvant settings.
Materials and methods
Cell culture
Request a detailed protocolDCIS.com‐Luc cells were generated by transducing cells with lentiviral expression particles for firefly luciferase 2 (LVP325; Amsbio, Abingdon, UK). DCIS.com-Luc-mCherry were created by lentiviral transduction of the mCherry sequence in the DCIS.com-Luc cells. Supplementary material (Supplementary file 6) provides details and sources of cell lines and media used. All cell lines were tested and were mycoplasma-free. The DCIS.com have been extensively profiled in Maguire et al., 2016 and we confirmed PIK3CA mutation status by PCR and Western blotting. Cell lines were authenticated by STR profiling (Eurofins).
Expression vectors
Request a detailed protocolThe SOX11 coding sequence (GENEID: 6664) from clone HsCD00295480 19 in the pENTR223.1 plasmid (DNASU) 20 was subcloned into pInducer21 (ORF-EG) plasmid gift from Stephen Elledge and Thomas Westbrook (Addgene plasmid # 46948; http://n2t.net/addgene: 46948; RRID:Addgene_46948) (Meerbrey et al., 2011). pInducer13-shRNA and pInducer13-NS shRNA were made by subcloning the specific shRNA174 for SOX11 and the NS shRNA from pGIPZ plasmids from Horizon into pInducer13.
RNA isolation
Request a detailed protocolRNA from cells grown in 2D for 48 hr in presence or absence of 1 μM doxycycline (DOX) and from spheroids treated for 2 or 5 days (DOX 2 days or DOX 5 days) or not treated with DOX (DOX 0 Days) (n = 3 for each time point) was isolated with an RNAeasyPlus Micro kit (74034; Qiagen, Manchester, UK) and DNase treatment. RNAClean and concentrator-5 (Zymo Research, Irvine, CA) were used. RNA concentration and purity were determined with a Qubit fluorometer (Invitrogen, Carlsbad, C) and a nanodrop spectrophotometer. RNA integrity number was measured with a bioanalyzer and an Agilent RNA Pico kit (Agilent Technologies, Cheshire, UK).
cDNA synthesis and qPCR
Request a detailed protocolOne microgram of each RNA sample was reverse transcribed with QuantiTect Reverse Transcription kit (Qiagen, Manchester, UK) in a final volume of 20 μl. cDNA was diluted ten times for subsequent quantitative polymerase chain reaction (qPCR) analysis, as described previously (Oliemuller et al., 2017), with the probes and methods listed in supplementary material (Supplementary file 7).
Western blotting
Request a detailed protocolWestern blotting was performed as previously described (Zvelebil et al., 2013). Details of the antibodies used are provided in supplementary material (Supplementary file 8).
Immunohistochemistry (IHC)
Request a detailed protocolIHC was performed on formalin-fixed paraffin embedded samples. Samples were stained with antibodies as described previously (Oliemuller et al., 2017).
Immunofluorescence
Request a detailed protocolAntibodies and staining protocols are detailed in supplementary material (Supplementary file 9). EVOS fluorescence microscope was used for imaging. Confocal images were captured with a Leica Microsystems (Cambridge, UK) TCS‐SP2 confocal microscope.
Spheroid formation
Request a detailed protocolFive thousand cells, untreated or treated with 1 μM doxycycline (DOX) for 48 hr, were plated per well in 96-well ultra-low-attachment plates (Corning 7007, Corning, NY, USA) in media containing DOX or not. After 24 hr, when the spheroids were formed, new media containing 1 μM doxycycline was added for 48 hr to spheres formed in absence of DOX for 2 days (DOX 2 days) or spheres that were formed in presence of DOX to a total of 5 days (DOX 5 days). Media without doxycycline was added to the control spheroids (No DOX). Images were obtained with a Celigo cytometer (Nexcelom, Manchester, UK).
Colony-formation assays
Request a detailed protocolDCIS.com cells were plated at 250 per well in six‐well (Falcon F3046, Corning, NY) plates. After 7 days, plates were stained with 0.2% crystal violet dissolved in 20% methanol in PBS. Area was measured, and the percentage relative to number of cells plated was calculated.
For single-cell colony assays and single-cell mammosphere assays, cells were FACS sorted and a single living cell plated per well in a 96-well or ultra-low attachment plate respectively. After 14 days, plates were analysed in a Celigo cytometer (Nexcelom, Manchester, UK).
For mammosphere assays, 5000 DCIS cells/ml were plated in low‐attachment six‐well plates (Corning 3471) and incubated in medium supplemented with 2% NeuroCult SM1 without vitamin A (StemCell Technologies) and 0.65% methylcellulose (R and D Systems, Abingdon, UK). After 14 days, wells were scanned with a Celigo cytometer. Mammosphere‐forming efficiency was calculated by dividing the number of mammospheres by the number of cells plated per well.
Cell viability assays
Request a detailed protocolThree thousand cells per well were plated in 96‐well plates (655098; Greiner Bio‐one, Stonehouse, UK) or in ultra-low attachment plates for 24 hr before starting the experiments. CellTiter‐Glo (Promega, Southampton, UK) was used according to the manufacturer's protocol. Luminescence was measured with a Victor X5 58 plate reader (Perkin‐Elmer, Seer Green, UK). In DCIS.com cells, CellTiter-Glo assays were performed at the time of adding the doxycycline and after 1, 2 and 3 days. When cells were transfected with siRNAs, the transfection was done overnight and next day (day 1), full supplemented media was added. CellTiter‐Glo assays were performed on days 1 and 5 to allow time for the siRNAs to knock down the genes of interest. CAL-148 cells stably transduced with shRNAs were measured at 0, 3 and 7 days, since this cell line has a long doubling time.
Invasion assays
Request a detailed protocolThree days after plating 5000 DCIS.com cells, spheroids were embedded in collagen I (354249; Corning) at 2.2 mg/ml diluted in medium. Complete medium was added on top after 1 hr. Images were acquired at this time and after 48 hr with a Celigo cytometer. The total area of matrix invaded by cells was calculated with ImageJ after marking of the area manually. For BT-20 and HCC1187 cells, the spheres used in the invasion assays were formed at the same time that they were transfected with the TUBB3 or MEX3A specific siRNAs in Opti-MEM media supplemented with 10% FBS. After 24 hr, this media was changed for the normal media of this cells. The cells were invaded in collagen 72 hr after plating.
Tranfections and siRNA
Request a detailed protocolBT-20, HCC1187 and CAL-148 cells were transfected with 50 pmol of each TUBB3 or MEX3A siRNA (siGENOME SMARTpool and four individual siRNAs), control nontargeting siRNA (Thermo Scientific, Waltham, MA) by using Lipofectamine RNAiMAX (Invitrogen, Life Technologies Corporation, Carlsbad, CA) in Opti-MEM media supplemented with 10% FBS (Gibco, Life Technologies Corporation, Carlsbad, CA) media according to the manufacturer's instructions for 16 hr in a six-well or 96-well plate, depending on the experiment, and then incubated with complete media. For sphere transfections, 5000 cells per well were incubated with 2.5 pmol of each specific siRNA in 150 µl of normal media per well in an ultra-low attachment plate for 16 hr. The next day, media was replaced with fresh media.
Flow cytometry analyses and fluorescence‐activated cell sorting (FACS) Aldehyde dehydrogenase (ALDH) activity was measured with the Aldefluor assay (StemCell Technologies, Cambridge, UK) as described before (Oliemuller et al., 2017). Cells were also co‐stained with Aldefluor and anti‐CD24–PE–Cy7 (561646) (1:100) and anti‐CD44–APC (559942) (1:20) (BD Biosciences, Oxford, UK). A BD FACS LSRII flow cytometer was used and samples were analysed with BD FACS Diva software (BD Biosciences). mCherry+ DCIS cells from the tumours in the tail vein xenograft experiments were sorted by a FACSAriaIII (BD). FACSAriaIII was used for clonogenic and mammosphere assays that required plating single cells per well. Living and dead cells were distinguished with DAPI 1:5000. Cell cycle analysis was described previously (Zvelebil et al., 2013).
PKH staining
Request a detailed protocoltwo × 10−6 DCIS-LacZ or DCIS-SOX11 cells were resuspended in Diluent C and stained with PKH26 Dye Solution to a final concentration of 5 mM following the manufacturer’s protocol (MINI26, Sigma). After confirming by flow cytometry that the 100% of the population stained, spheres with 5000 cells were formed in ultra-low attachment 96-well plates.
Animal experiments
Request a detailed protocolAll animal work was carried out under UK Home Office project and personal licenses following local ethical approval from The Institute of Cancer Research Ethics Committee and in accordance with local and national guidelines. For xenograft tumour assays, DCIS.com labelled with Luc2-mCherry and stably transduced with empty vector (EV) control or SOX11 (in pInducer21 backbone) were resuspended in PBS for implantation into female NSG-Foxn1null mice. 1.0 × 106 cells/site in both sides were injected in each mammary fat pad number 4 and 2.0 × 104 cells/site into the mammary duct via the nipple of mammary gland number four as previously described (Oliemuller et al., 2017).
Intracranial injections were performed using a stereotaxic frame by injecting 1 × 10−5 cells (Stoelting, Wood Lane, IL, USA) into the striatum (2 mm right from the midline, 2 mm anterior from bregma, 3 mm deep). Six NSG-Foxn1null mice were used per for each condition.
For tail vein injections, the mice were placed into a hot box set at 38°C for up to 5 min, to dilate the tail veins. The mice were then placed into a restrainer, the lateral tail vein identified, and 2.5 × 105 cells were slowly infused (over 30 s) through the tail vein using a venoflux 25 g butterfly needle. Before treatment, mice used were randomised in groups based on their weight.
For the doxycycline induction experiments using iSOX11 or iEV, a week before the injections, the animals injected with transduced cells were separated into two cohorts and maintained with or without chow containing doxycycline (0.2 g/kg, 0.625, 1.250 or 2.0 g/kg from Envigo) for the duration of the experiment.
For in vivo imaging, mice were injected with 200 μl of luciferin (XenoLight D-Luciferin Potassium Salt, Perkin Elmer, 30 mg/mL). After 5 min, the mice were anaesthetised with isoflurane and the bioluminescence was measured at least at three different time points on an IVIS Lumina imaging systems (Perkin-Elmer). For ex vivo imaging, the organs were resected and imaged for 3 min to detect any signal. In the analysis, identical square regions of interest (ROI) were drawn around tumours to measure total and average bioluminescence signal.
Tumour dissociation iSOX11 or iEV cells from brain, bone and lung metastasis were isolated with a cell dissociation kit following the manufacturer instructions (Miltenyi Biotec, Bergisch Gladbach, Germany) and using the gentleMACS Octo Dissociator with Heaters and its 37°C_h_TDK_3 protocol. To select human cells and discard murine cells, the result of this dissociation was FACS sorted in a FACSAriaIII and mCherry+ cells were selected and grown in normal DCIS media.
RNA sequencing cDNA library preparation was carried out at Oxford Genomics Centre, The Wellcome Trust Centre for Human Genetics using PolyA+ RNA enrichment method for total RNA from cultured cells. mRNA fraction was selected from the total RNA before conversion to cDNA. Second strand cDNA synthesis incorporated dUTP. The cDNA was end-repaired, A-tailed and adapter-ligated. Prior to amplification, samples underwent uridine digestion. The prepared libraries were size selected, multiplexed and quality checked before paired end sequencing over four lanes of a flow cell.
Sequence files were trimmed by the use of trim_galore (http://www.bioinformatics.babraham.ac.uk/projects/trim_galore/) with default settings. Trimmed data were separately mapped to the GRCh38 and GRCm38 genome assemblies by the use of hisat2 (v2.0.5) with options ‐‐sp 1000,1000 ‐‐ omixed‐‐no‐discordant, and were filtered to remove non‐primary alignments. Species‐specific read sets were generated by removing any read that produced a valid alignment in both human and mouse from the results for both species. The remaining data were imported into SeqMonk (http://www.bioinformatics.babraham.ac.uk/projects/seqmonk/) with a filter of mapping quality (MAPQ) score ≥20. Reads were quantified over the transcript set from Ensembl v78 with annotated mis‐spliced, pseudogene and unannotated transcripts removed. Initial quantification was raw read counts from the opposing strand to the transcript, with all exons for each gene being collated into a single measure. This allowed gene‐level differential expression to be assessed by the use of DESeq2 (https://bioconductor.org/packages/release/bioc/html/DESeq2.html), with a cutoff of a false discovery rate of <0.05. Subsequent visualisation was performed by requantifying expression as log2 fragments per million reads of library. RNA sequencing files were submitted to ArrayExpress.
Survival analysis
Request a detailed protocolThe prognostic importance of SOX11 mRNA expression was assessed by the use of survival data using the cBioPortal for Cancer Genomics (http://cbioportal.org) (Cerami et al., 2012; Gao et al., 2013). Data obtained from The Cancer Genome Atlas (https://www.cancer.gov/tcga) was examined by use of the Kaplan–Meier Plotter survival analysis tool (http://kmplot.com) and METABRIC (METABRIC Group et al., 2012), and statistical significance was determined with the Wald test.
Statistical analysis
Request a detailed protocolThe data in the graphs are presented as mean and standard deviation, unless specified otherwise. Experiments were analysed with a two‐tailed Student's t‐test with a confidence interval of 95% when the number of groups equalled 2, or with a parametric ANOVA and post hoc test when the number of groups was >2, unless otherwise specified. p-value≤0.0001 is considered as extremely significant (****), p≤0.001 as highly significant (***), p≤0.01 as very significant (**), p≤0.05 as significant (*), and p>0.05 as not significant (ns), respectively.
Data availability
Sequencing data have been deposited in ArrayExpress as accession E-MTAB-9108. All data generated or analysed during this study are included in the manuscript and supporting files.
-
ArrayExpressID EBI. RNA-seq of DCIS-pInducer21-SOX11 cells grown in 2D and 3D.
References
-
Phosphorylation modulates the subcellular localization of SOX11Frontiers in Molecular Neuroscience 11:211.https://doi.org/10.3389/fnmol.2018.00211
-
Mammary lineage restriction in developmentNature Cell Biology 20:637–639.https://doi.org/10.1038/s41556-018-0111-6
-
Molecular regulation of stem cell quiescenceNature Reviews Molecular Cell Biology 14:329–340.https://doi.org/10.1038/nrm3591
-
Lamin B constitutes an intermediate filament attachment site at the nuclear envelopeThe Journal of Cell Biology 105:117–125.https://doi.org/10.1083/jcb.105.1.117
-
Basal-like breast cancers: from pathology to biology and back againStem Cell Reports 10:1676–1686.https://doi.org/10.1016/j.stemcr.2018.04.023
-
Cadherins in development: cell adhesion, sorting, and tissue morphogenesisGenes & Development 20:3199–3214.https://doi.org/10.1101/gad.1486806
-
Basal cells of second trimester fetal breasts: immunohistochemical study of myoepithelial precursorsPediatric and Developmental Pathology 6:398–413.https://doi.org/10.1007/s10024-003-1125-y
-
βIII-Tubulin: biomarker of taxane resistance or drug target?Expert Opinion on Therapeutic Targets 17:461–472.https://doi.org/10.1517/14728222.2013.766170
-
Gene expression profiling of breast Cancer brain metastasisScientific Reports 6:28623.https://doi.org/10.1038/srep28623
-
Emerging role of SOX proteins in breast Cancer development and maintenanceJournal of Mammary Gland Biology and Neoplasia 24:213–230.https://doi.org/10.1007/s10911-019-09430-6
-
N-cadherin promotes motility in human breast Cancer cells regardless of their E-cadherin expressionJournal of Cell Biology 147:631–644.https://doi.org/10.1083/jcb.147.3.631
-
Reprogramming to developmental plasticity in Cancer stem cellsDevelopmental Biology 430:266–274.https://doi.org/10.1016/j.ydbio.2017.07.025
-
SOX11 promotes invasive growth and ductal carcinoma in situ progressionThe Journal of Pathology 243:193–207.https://doi.org/10.1002/path.4939
-
CDX2 regulation by the RNA-binding protein MEX3A: impact on intestinal differentiation and stemnessNucleic Acids Research 41:3986–3999.https://doi.org/10.1093/nar/gkt087
-
Breast Cancer brain metastases: clonal evolution in clinical contextInternational Journal of Molecular Sciences 18:152.https://doi.org/10.3390/ijms18010152
-
Regulatory roles for SOX11 in development, stem cells and CancerSeminars in Cancer Biology 20:15.https://doi.org/10.1016/j.semcancer.2020.06.015
-
Transcriptome characterization of matched primary breast and brain metastatic tumors to detect novel actionable targetsJNCI: Journal of the National Cancer Institute 111:388–398.https://doi.org/10.1093/jnci/djy110
-
MEX3A knockdown inhibits the development of pancreatic ductal adenocarcinomaCancer Cell International 20:63.https://doi.org/10.1186/s12935-020-1146-x
-
Early lineage segregation of multipotent embryonic mammary gland progenitorsNature Cell Biology 20:666–676.https://doi.org/10.1038/s41556-018-0095-2
Article and author information
Author details
Funding
Breast Cancer Now (Programme Funding to the Breast Cancer Now Toby Robins Research Centre)
- Erik Oliemuller
- Richard Newman
- Siu Man Tsang
- Gareth Muirhead
- Farzana Noor
- Syed Haider
- Beatrice A Howard
Cancer Research UK (CRUK: A21855)
- Shane Foo
Spanish Ministry of Education and Science (SAF2017-84934-R)
- Maria dM Vivanco
Government of the Autonomous Community of the Basque Country, the Department of Industry, Tourism and Trade (KK-2018/00054)
- Iskander Aurrekoetxea-Rodríguez
- Maria dM Vivanco
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Breast Cancer Now for funding this work as part of Programme Funding to the Breast Cancer Now Toby Robins Research Centre. This work was supported by a grant from CRUK (CRUK: A21855). This work was supported by grants from the Spanish Ministry of Education and Science (SAF2017-84934-R, to MV) and the Government of the Autonomous Community of the Basque Country, the Department of Industry, Tourism and Trade (Elkartek: KK-2018/00054, IA-R and MV). We acknowledge the role of the Breast Cancer Now Tissue Bank in collecting and making available the normal breast samples used in the generation of this publication. We thank the Oxford Genomics Centre at the Wellcome Centre for Human Genetics (funded by Wellcome Trust grant reference 203141/Z/16/Z) for the generation and initial processing of the sequencing data.
Ethics
Animal experimentation: All animal work was carried out under UK Home Office PPL number: PB0FA698C (BAH Project licence holder) and personal licenses following local ethical approval of all protocols from The Institute of Cancer Research Ethics Committee and in accordance with local and national guidelines. All surgery was performed under isoflurane anesthesia with appropriate analgesia, and every effort was made to minimise suffering.
Copyright
© 2020, Oliemuller et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,215
- views
-
- 287
- downloads
-
- 30
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Cancer Biology
- Computational and Systems Biology
Effects from aging in single cells are heterogenous, whereas at the organ- and tissue-levels aging phenotypes tend to appear as stereotypical changes. The mammary epithelium is a bilayer of two major phenotypically and functionally distinct cell lineages: luminal epithelial and myoepithelial cells. Mammary luminal epithelia exhibit substantial stereotypical changes with age that merit attention because these cells are the putative cells-of-origin for breast cancers. We hypothesize that effects from aging that impinge upon maintenance of lineage fidelity increase susceptibility to cancer initiation. We generated and analyzed transcriptomes from primary luminal epithelial and myoepithelial cells from younger <30 (y)ears old and older >55y women. In addition to age-dependent directional changes in gene expression, we observed increased transcriptional variance with age that contributed to genome-wide loss of lineage fidelity. Age-dependent variant responses were common to both lineages, whereas directional changes were almost exclusively detected in luminal epithelia and involved altered regulation of chromatin and genome organizers such as SATB1. Epithelial expression of gap junction protein GJB6 increased with age, and modulation of GJB6 expression in heterochronous co-cultures revealed that it provided a communication conduit from myoepithelial cells that drove directional change in luminal cells. Age-dependent luminal transcriptomes comprised a prominent signal that could be detected in bulk tissue during aging and transition into cancers. A machine learning classifier based on luminal-specific aging distinguished normal from cancer tissue and was highly predictive of breast cancer subtype. We speculate that luminal epithelia are the ultimate site of integration of the variant responses to aging in their surrounding tissue, and that their emergent phenotype both endows cells with the ability to become cancer-cells-of-origin and represents a biosensor that presages cancer susceptibility.
-
- Cancer Biology
Metastasis is the leading cause of cancer-related mortality. Paneth cells provide stem cell niche factors in homeostatic conditions, but the underlying mechanisms of cancer stem cell niche development are unclear. Here, we report that Dickkopf-2 (DKK2) is essential for the generation of cancer cells with Paneth cell properties during colon cancer metastasis. Splenic injection of Dkk2 knockout (KO) cancer organoids into C57BL/6 mice resulted in a significant reduction of liver metastases. Transcriptome analysis showed reduction of Paneth cell markers such as lysozymes in KO organoids. Single-cell RNA sequencing analyses of murine metastasized colon cancer cells and patient samples identified the presence of lysozyme positive cells with Paneth cell properties including enhanced glycolysis. Further analyses of transcriptome and chromatin accessibility suggested hepatocyte nuclear factor 4 alpha (HNF4A) as a downstream target of DKK2. Chromatin immunoprecipitation followed by sequencing analysis revealed that HNF4A binds to the promoter region of Sox9, a well-known transcription factor for Paneth cell differentiation. In the liver metastatic foci, DKK2 knockout rescued HNF4A protein levels followed by reduction of lysozyme positive cancer cells. Taken together, DKK2-mediated reduction of HNF4A protein promotes the generation of lysozyme positive cancer cells with Paneth cell properties in the metastasized colon cancers.