Osterix-Cre marks distinct subsets of CD45- and CD45+ stromal populations in extra-skeletal tumors with pro-tumorigenic characteristics
Abstract
Cancer-associated fibroblasts (CAFs) are a heterogeneous population of mesenchymal cells supporting tumor progression, whose origin remains to be fully elucidated. Osterix (Osx) is a marker of osteogenic differentiation, expressed in skeletal progenitor stem cells and bone-forming osteoblasts. We report Osx expression in CAFs and by using Osx-cre;TdTomato reporter mice we confirm the presence and pro-tumorigenic function of TdTOSX+ cells in extra-skeletal tumors. Surprisingly, only a minority of TdTOSX+ cells expresses fibroblast and osteogenic markers. The majority of TdTOSX+ cells express the hematopoietic marker CD45, have a genetic and phenotypic profile resembling that of tumor infiltrating myeloid and lymphoid populations, but with higher expression of lymphocytic immune suppressive genes. We find Osx transcript and Osx protein expression early during hematopoiesis, in subsets of hematopoietic stem cells and multipotent progenitor populations. Our results indicate that Osx marks distinct tumor promoting CD45- and CD45+ populations and challenge the dogma that Osx is expressed exclusively in cells of mesenchymal origin.
Introduction
In the past couple of decades the common beliefs of how an incipient tumor grows and progresses to a metastatic stage have drastically changed due to the increasing findings of a tight crosstalk between tumor cells and the surrounding stroma, known as the tumor microenvironment (TME) (Hanahan and Weinberg, 2011). Tumor stroma comprises a variety of different cells from the mesenchymal and hematopoietic compartments, with pro- and anti-tumor functions. Cancer-associated fibroblasts (CAFs) are cells of the mesenchymal lineage involved in supporting various stages of tumorigenesis from growth to metastatic dissemination (Kalluri and Zeisberg, 2006; Chen and Song, 2019). Their abundance and phenotypic markers differ among the various types of cancers. CAFs support tumor growth by stimulating tumor cell proliferation, inhibiting apoptotic signals and anti-tumor immune responses, altering the extra-cellular matrix (ECM) to favor invasiveness, and by offering protection from chemotherapeutic approaches. Their heterogeneity probably accounts for the different tumor promoting effects and reflects the fact that CAFs can originate from multiple cellular sources (Chen and Song, 2019). Several studies indicate that CAFs derive from bone marrow mesenchymal stem and progenitor cells recruited at the tumor site. Additional evidence indicates that they can also derive from myofibroblasts, as well as from the trans-differentiation of local pericytes, endothelial and epithelial cells. In breast cancer, subsets of bone marrow-derived CAFs can exhibit a unique inflammatory profile depending on the location to which they are recruited, thus being functionally distinct from the resident CAFs and having better tumor-promoting functions (Raz et al., 2018). Adding more complexity to the origin of CAFs, an additional population expressing several markers commonly found in fibroblasts and CAFs, is the fibrocyte (Abe et al., 2001; McDonald and LaRue, 2012; van Deventer et al., 2013). Fibrocytes originate from the bone marrow, where they have been described to contribute to bone marrow fibrosis under pathological conditions (Ohishi et al., 2012). These cells have also been implicated in supporting tumor progression, but they differ from the bone marrow-derived CAF populations in that they also express the hematopoietic marker CD45, along with markers associated with the monocyte/macrophage lineage populations (Abe et al., 2001). Fibrocytes support tumor growth by making collagen, albeit at lower levels than the CAFs, releasing growth factors, commonly produced by immune suppressive myeloid cells, and by modulating resistance to anti-angiogenic therapy (Goto and Nishioka, 2017). Recent single-cell RNA sequencing (scRNAseq) studies further confirmed the complexity of the CAF populations, suggesting tumor specificity as well as functional differences among the various subsets (Bartoschek et al., 2018; Costa et al., 2018; Elyada et al., 2019). However, these studies did not directly address the origin of the various subsets, which could account for their specific cellular phenotype.
Because of their ability to produce collagen and other extracellular matrix proteins, we sought to determine whether a subset of bone marrow-derived CAFs could share markers of committed osteolineage cells. Osteolineage cells derive from bone marrow mesenchymal stem cells (MSCs) and their differentiation program is driven by sequential activation of two specific transcription factors, Runx2 and Sp7(Osx gene), and subsequent acquisition of a phenotype characterized by the ability to produce bone matrix, primarily type I collagen, and to mineralize (Nakashima et al., 2002). Although Osx is largely thought as a marker of differentiated osteoblasts, emerging data demonstrate that in the embryonic and perinatal bone marrow Sp7 is expressed in definitive MSCs that give rise to the marrow stroma, including osteoblasts and adipocytes (Liu et al., 2013; Mizoguchi et al., 2014). Furthermore, during embryogenesis, Osx is present in extra-skeletal tissues, including the olfactory bulb, the intestine and the kidney (Chen et al., 2014; Jia et al., 2015).
Based on the above observations, we hypothesized that a subset of CAFs, derived from Osx+ cells in the bone marrow, contributes to ECM (i.e. collagen) production at the tumor site, thereby creating a tumor supporting stroma. Using a cell tracking system, we found the presence of cells targeted by the Osx promoter within the TME; these TdTOSX+ cells favor tumor growth when co-injected with tumor cells in mice. Surprisingly, only a minority of tumor-resident cells derived from Osx+ cells expresses fibroblast markers, extracellular matrix and matrix remodeling genes. The majority of these newly identified TdTOSX+ tumor infiltrating cells are also positive for CD45, a marker of hematopoietic lineage, and share markers expressed by tumor-infiltrating immune cells. Importantly, we confirmed Sp7 transcripts and Osx protein in a subset of hematopoietic stem cells (HSC), giving rise to TdTOSX+;CD45+ tumor infiltrating immune populations. This study further identifies new populations of TME cells targeted by Osx and challenges the use of Osx-cre driven lineage tracing mouse models to exclusively study mesenchymal lineage cell fate.
Results
Embryonic and adult-derived osteolineage Osx+ cells are present in extra-skeletal tumors
To determine whether osteolineage cells may be present in the TME, we crossed the established tetracycline-dependent Tg(Sp7-tTA,tetO-EGFP/cre)1Amc/J (Osx-cre) to the B6.Cg-Gt(ROSA)26Sortm9(CAG-tdTomato)Hze/J (TdT) to generate the Osx-cre;TdT reporter mouse model (Rodda and McMahon, 2006). When constitutively activated, TdT marks the entire osteolineage, including bone surface osteoblasts, osteocytes and bone marrow cells with mesenchymal stem and osteoprogenitor cell features; while delaying Osx-cre expression until postnatally restricts TdT targeting to committed osteoblasts and osteocytes (Mizoguchi et al., 2014; Fontana et al., 2017). Therefore, Osx-cre;TdT mice and control animals carrying only the TdT transgene (WT;TdT) were kept on standard chow to allow constitutive embryonic transgene activation, or fed a doxycycline (doxy)-containing diet until weaning to suppress transgene activation until one month of age (Figure 1A). We previously reported that doxy-fed mice display less than 1% of spontaneous recombination in the bone residing osteoblasts at weaning, but full transgene activation 1 month thereafter (Fontana et al., 2017).
Osx-cre;TdT reporter mice and WT;TdT fed a normal diet (no doxy) were inoculated with either 105 B16-F10 melanoma cells subcutaneously or 105 PyMT breast cancer cells in the mammary fat pad (MFP). To account for the possible growth delay sometimes observed in no-doxy fed Osx-cre mice (Davey et al., 2012), tumor cells were injected at 12 weeks of age, when the growth delay is fully recovered. Importantly, no differences in tumor growth were observed in age and sex matched Osx-cre;TdT and WT;TdT mice, excluding any potential confounding effect of the Osx-cre (Figure 1—figure supplement 1). To assess the presence of osteolineage Osx+ derived cells within the tumor stroma of no doxy-fed mice, we performed flow cytometry analysis on B16-F10 and PyMT tumors isolated 2 weeks post-inoculation, and found that about 10–18% of the total cells were TdTOSX+ (Figure 1B and D). Histological analysis of frozen sections from the same tumors confirmed the presence of TdTOSX+ cells. These cells had either round or elongated shape, indicating morphologic heterogeneity within the TdTOSX+ population (Figure 1C and E). These TdTOSX+ cells within the tumor stroma of no-doxy mice may derive from resident Osx-expressing embryonic progenitors or from mobilization of bone marrow Osx+ cells.
Remarkably, TdTOSX+ cells were also present in the B16-F10 and PyMT TME of doxy-fed animals 4 weeks post weaning, although their numbers were lower than in tumors isolated from no doxy-treated mice (Figure 1F and H). Fluorescence micrographs also confirmed the presence of TdTOSX+ cells in the stroma in both tumor models (Figure 1G and I), with the same pleiotropic morphology as seen in no doxy-treated mice. Thus, embryonic and adult-derived osteolineage Osx+ cells are present in the stroma of extra-skeletal tumors. Based on these findings, the doxy-fed Osx-cre;TdT mice were used for all the subsequent studies.
TdTOSX+ cells from primary tumors express mesenchymal markers
To determine the phenotype of TdTOSX+ cells more in depth, we isolated TdTOSX+ cells by FACS sorting 14 days after inoculation of B16-F10 tumor cells into doxy-fed Osx-cre;TdT reporter mice (Figure 2A). Due to the limited number of TdTOSX+, tumors from 2 to 4 mice were pooled for FACS sorting and reported as one experimental replicate (single data point). Messenger RNA from an immortalized murine CAF cell line (i-CAFs) and cortical long bone devoid of bone marrow cells were used as controls. Real-Time PCR confirmed Sp7 expression in TdTOSX+ cells but not in the TdTOSX– fraction, which includes the remaining stroma and the tumor cells (Figure 2B). Interestingly, Sp7 expression was also detected in i-CAFs (Figure 2B), and as expected in bone extracts in higher abundance. Expression of osteoblast-specific markers such as Runx2, Bglap (Ocn) (Osteocalcin), Ibsp (Bsp, bone sialoprotein) and Alpl (Tnap, tissue non-specific alkaline phosphatase) were negligible in TdTOSX+ cells and undetectable or barely detectable in the i-CAFs (Figure 2C–F). However, S100a4(Fsp1) (Fibroblast-specific protein 1) and Acta2 (α-SMA) (α-smooth muscle actin), two CAFs specific markers, were detected in TdTOSX+ cells (Figure 2G and H). TdTOSX+ also expressed Col1a1 (Collagen Type Ia1) and Col1a2 (Collagen Type Ia2) similar to i-CAFs and bone extracts (Figure 2I and J). Expression of matrix metalloproteinases, Mmp2 and Mmp9, involved in matrix remodeling and highly expressed in CAFs and osteoblasts respectively, was also detected in the TdTOSX+ cells (Figure 2K and L).
To validate expression of Sp7 in CAFs from primary tumors, we orthotopically injected B16-F10 or PyMT cells in 8 week old WT mice and the CAFs were isolated as the adherent fraction of a single cell suspension of the tumor mass. We detected Sp7 expression only in CAFs but not in the tumor cell lines nor in the non-adherent fraction of the single cell suspension from the tumor mass (Figure 2M and N). Importantly, corroborating our findings in mice, using Oncomine (Rhodes et al., 2004) we found that microarray analysis of human breast stromal cells collected from patients with invasive breast carcinoma from the Finak gene set (Finak et al., 2008) showed significantly higher Sp7 expression in the tumor associated stroma compared to the adjacent normal breast tissue (Figure 2O). Thus, a subset of cells in the TME of mice and humans express Sp7, an osteolineage cell marker, but also markers associated with the CAFs. TdTOSX+ cells may represent a subpopulation of CAFs with specific features.
Tumor resident but not bone marrow resident TdTOSX+ cells increase tumor growth
To determine whether tumor resident TdTOSX+ cells can affect tumor growth as do CAFs, TdTOSX+ cells were sorted from primary B16-F10 tumors inoculated in 8 week old doxy-fed Osx-cre;TdT mice and re-injected along with B16-F10 tumor cells at 5:1 ratio into age-matched WT recipient animals (Figure 3A). Remarkably, co-injecting TdTOSX+ with tumor cells resulted in larger tumors (Figure 3B) compared to injection of B16-F10 cells alone.
Because TdTOSX+ cells normally reside in the bone microenvironment (Fontana et al., 2017; Strecker et al., 2013), we next examined bone sections from 8 week old doxy-fed Osx-cre;TdT reporter mice, naïve or bearing soft tissue B16-F10 or PyMT tumors. We confirmed presence of TdTOSX+ cells on the surface of cortical and trabecular bone, TdTOSX+ osteocytes, and also abundant staining in the bone marrow. TdT was not detected WT;TdT mice (Figure 3C and Figure 3—figure supplement 1A). We further analyzed the percentage of TdTOSX+ cells in the bone marrow by FACS and found that TdTOSX+ cells accounted for about 10% of total marrow cells in naïve mice and their number increased in mice injected with B16-F10 or PyMT tumor cells (Figure 3D and Figure 3—figure supplement 2B). Similar findings were observed in the Tamoxifen (TAM) inducible Sp7-creERT2; Rosa26<fs-TdTomato> (Osx-creERT2;TdT) mouse model. TAM was administered to 7 week old mice via IP injections starting 3 days prior to B16-F10 tumor inoculation, and 1, 6, and 9 days post tumor cell injection (Figure 4—figure supplement 1A). No tumor bearing mice were used as control. Similar to the Tet-OFF Osx-cre model, we found that bone marrow TdTOSX+ cells significantly increased in presence of a tumor (Figure 3E).
To test their capacity to promote tumor growth, bone marrow-derived TdTOSX+ cells were isolated from tumor bearing mice and re-injected with B16-F10 tumor cells into WT mice (5:1 ratio). In contrast to what observed with tumor-derived TdTOSX+ cells, bone marrow-derived TdTOSX+ cells did not enhance tumor growth compared to mice injected with tumor cells alone; if anything, there was a trend towards a smaller tumor size (Figure 3F). This result further prompted us to isolate bone marrow-derived TdTOSX+ cells and analyze the expression levels of the same mesenchymal genes analyzed in the tumor-derived TdTOSX+ cells (Figure 2). Interestingly, TdTOSX+ cells from the TME expressed higher levels of Sp7, CAF markers S100a4 and Acta2, and matrix proteins Col1a1 and Ibsp compared to TdTOSX+ from the bone marrow (Figure 3—figure supplement 2). Bone marrow TdTOSX+ cells instead expressed significantly higher levels of the early osteoblast transcription factor Runx2 and of the protease Mmp9, while no differences in other osteoblast markers such as Bglap and Alpl were detected between the tumor and bone marrow sorted TdTOSX+ cells (Figure 3—figure supplement 2).
These results suggest that Osx+ cells infiltrating a tumor are functionally and phenotypically distinct from bone marrow resident Osx+ cells.
The majority of TdTOSX+ cells express CD45
To gain insights on the possible origin of tumor-derived TdTOSX+ cells, we performed FACS analysis of circulating cells in naïve mice and two weeks post B16-F10 or PyMT tumor inoculation. Interestingly, we found 5–8% of TdTOSX+ cells in the blood of tumor-free doxy-fed Osx-cre;TdT mice, and this percentage was 3- to 4-fold higher in tumor-bearing animals (Figure 4A and B). The unexpected high number of circulating TdTOSX+ cells prompted us to determine whether they may express the immune marker CD45, which is also present in fibrocytes. Surprisingly, 95% of TdTOSX+ cells in circulation of tumor bearing mice were also positive for CD45, and they represented about 13–18% of total blood cells (Figure 4C and D). Similarly, about 92–95% of total TdTOSX+ cells in the bone marrow were CD45+ (Figure 4E and F), representing 12–20% of total marrow cells. Such finding explained the low expression levels of mesenchymal markers in the bone marrow-derived TdTOSX+ cells (Figure 3—figure supplement 2). Intriguingly, also the majority of TdTOSX+ cells in the TME expressed CD45 (Figure 4G and H), data further confirmed by RT-PCR in TdTOSX+ fraction sorted from B16-F10 tumors (not shown).
To exclude a possible off-target effect of the doxy-dependent Tet-OFF system, we turned to the TAM inducible Osx-creERT2;TdT model. FACS analysis confirmed presence of both CD45- and CD45+ TdTOSX+ populations in the bone marrow of B16-F10 bearing Osx-creERT2;TdT mice, but not in the cre negative mice (Figure 4—figure supplement 1B).
Next, to confirm that expression of CD45 in the TdTOSX+ cells was not due a tight interaction between a CD45+ immune cell and TdTOSX+ mesenchymal cell, we performed CD45 immunostaining on paraffin-embedded B16-F10 tumor sections and on bone marrow single cell suspension from Osx-cre;TdT reporter mice. We found co-localization of Osx-driven TdTomato and CD45 in the same cell, confirming that CD45 is indeed expressed in subsets of TdTOSX+ cells (Figure 4I and Figure 4—figure supplement 2). Thus, Osx not only marks mesenchymal cells, but also cells expressing the hematopoietic marker CD45, suggesting that they could represent bone marrow fibrocytes or an immune cell subset.
Presence of functionally distinct populations of TdTOSX+ cells in the tumor microenvironment
To further characterize the phenotype of the CD45- and CD45+ TdTOSX+ populations, we next performed RNAseq analysis. We subcutaneously injected the GFP-labeled PyMT-BO1-tumor line in Osx-cre;TdT reporter mice. Two weeks post tumor inoculation, the stromal cells were separated from the GFP+ tumor cells and sorted based on the expression, or lack of thereof, of TdT and CD45 markers. We sorted 4 groups of cells: double negative (TdTOSX-;CD45-), representing the non-immune tumor stroma, CD45 single positive (TdTOSX-;CD45+), representing the tumor infiltrating immune populations, TdTOSX single positive (TdTOSX+;CD45-), and double positive (TdTOSX+;CD45+).
Principal components analysis of RNA-seq expression patterns across all four groups revealed clustering that was uniquely dependent on CD45 expression and suggested that there were only two very distinct cell populations regardless of TdTOSX (Figure 5A). Inspection of each cell population in a one versus all other approach for only statistically significant up-regulated genes (FDR >= 0.05, log2 fold-change >= 2) to identify robustly expressed biomarkers, showed that all four groups shared many genes in common based on CD45 expression alone. The CD45 negative populations shared 1248 genes in common with 536 in the TdTOSX single positive and 477 in the double negative cells uniquely up-regulated (Figure 5B). Likewise, the CD45 single and double positive cells shared 495 significantly up-regulated genes with an additional 77 and 316 uniquely expressed, respectively.
Subsequent Gene Ontology (GO) analysis of the log2 fold-changes for each cell population revealed that the CD45 negative subsets showed upregulation of GO pathways related to cellular responses to growth factors, extracellular structure and matrix organization and skeletal system development, confirming their mesenchymal nature (Figure 5C). Among the top upregulated genes in the TdTOSX single positive cells versus the double negative (Table 1) we found several markers expressed by osteolinage cells, such as fibromodulin (Fmod), a binding protein regulating bone mineralization (Gori et al., 2001) and also involved in TGFβ signaling during cancer pathogenesis (Pourhanifeh et al., 2019), Collagen 24a1 (Col24a1), known to modulate collagen chain trimerization (Koch et al., 2003; Matsuo et al., 2006), Frem1, a protein involved in the formation and organization of basement membranes (Vissers et al., 2011), and Msx1, a transcription factor important in craniofacial bone development (Orestes-Cardoso et al., 2002). TdTOSX single positive cells also expressed high levels of fibulin 7 (Fbln7), a cell adhesion molecule overexpressed in glioblastoma by pericytes and involved in neovascularization (de Vega et al., 2019; Ikeuchi et al., 2018). Such result confirmed the osteolineage nature of the TdTOSX+ cells and indicated their tumor supporting role.
Conversely, both CD45+ subsets were significantly up-regulated for GO pathways related to regulation of lymphocyte activation and proliferation, and innate immune cell responses (Figure 5D). Interestingly, among the top genes upregulated in the TdTOSX+ CD45+ subsets versus the CD45 single positive populations (Table 2), we found many genes expressed by tumor promoting immune cells, such as regulatory T cells (i.e. Foxp3) and gamma-delta T cells (i.e. Tcrg-C1). The most upregulated gene in the TdTOSX+;CD45+ subset (>80 fold versus CD45 single positive) was Cd163l1, a marker of gamma-delta IL-17 producing cells (Tan et al., 2019) and M2 macrophages (González-Domínguez et al., 2015), two tumor-promoting immune populations. We confirmed Cd163l1 higher expression by RT-PCR together with other two genes, Wdr78 and Spock2, also upregulated in the double positive cells (Figure 5—figure supplement 1). The second most expressed gene was Lymphotoxin beta (Ltb), a cytokine produced by lymphocytes and NK cells and associated with carcinogenesis (Wolf et al., 2010). The third one, Klrg1, negatively regulates cytotoxic lymphocytes and is associated with lymphocyte senescence and dysfunction (Henson and Akbar, 2009). Thus, these analyses revealed that two main transcriptional programs predominate based on CD45 expression, and suggest that TdTOSX+;CD45- comprise a subset of CAFs with some characteristics of skeletal cells, while TdTOSX+;CD45+ cells represent a heterogeneous subset of tumor-promoting immune cells.
Double positive TdTOSX+;CD45+ cells are a heterogeneous immune population enriched in lymphoid cells
Next we performed flow cytometry analysis to better characterize the TdTOSX+;CD45+ population and validate some of the RNAseq findings using the B16-F10 melanoma model. We confirmed that the double positive TdTOSX+;CD45+ represented about 20% of the total CD45+ cells in both the primary tumor and the bone marrow of doxy-fed Osx-cre;TdT mice (Figure 6A and B). Next, we analyzed the percentage of TdTOSX+ expressing the common myeloid and lymphoid markers (gate strategy in Figure 6—figure supplement 1), such as CD11b (monocytes), F4/80 (macrophages), Gr1 (granulocytes/neutrophils), CD3 (T cells, further divided into CD4+ and CD8+), and NK1.1 (Natural Killer cells). To standardize the comparison between the TdTOSX+;CD45+ and TdTOSX-;CD45+ (single positive) populations within each sample, data were represented as percentage of either total TdTOSX+;CD45+ or total CD45 single positive cells (Figure 6C and D). TdTOSX+;CD45+ subset expressed both myeloid and lymphoid markers, with a similar pattern as the CD45 single positive cells. In both subsets, the majority of the tumor immune infiltrate included myeloid populations, which represented about 60% of the TdTOSX+;CD45+ and over 70% of the total CD45 single positive cells. Lymphoid cells appeared to be more abundant among the TdTOSX+;CD45+ (32.3 ± 6.5%), compared to the TdTOSX-;CD45+ cells (23.3 ± 4.3%). Based on this observation, we calculated the ratio of lymphoid over myeloid populations and observed that the tumor-infiltrating TdTOSX+;CD45+ double positive population was enriched in lymphoid cells (Figure 6E). Flow cytometric analysis of the bone marrow from tumor bearing Osx-cre;TdT mice revealed similar pattern of distribution between the TdTOSX+;CD45+ and the CD45 single positive immune populations (Figure 6F and G), although the frequencies of each immune subset in the bone marrow differed compared to the cells at tumor site. These data suggest that Osx marks multiple immune cell types but that TdTOSX+-derived immune cells infiltrating a tumor are skewed towards lymphoid populations.
Since adoptive transfer of tumor-derived TdTOSX+ cells increases tumor growth (Figure 3B), we asked whether TdTOSX+;CD45+ immune populations contribute to tumor progression. We isolated by FACS sorting TdTOSX+;CD45+ and the TdTOSX-;CD45+ from B16-F10 subcutaneous tumors from doxy-fed Osx-Cre;TdT mice. Sorted cells were then re-injected along with B16-F10 tumor cells (ratio 5:1) into new WT recipient mice. The co-injection of CD45 single positive cells decreased tumor growth relative to B16-F10 injected alone at day 15, while TdTOSX+;CD45+ did not show anti-tumor effects (Figure 6H).
Collectively, these data indicate that Sp7 marks cells of both mesenchymal and hematopoietic origin to support the development of tumor promoting populations.
Double positive TdTOSX+;CD45+ cells derive from TdTOSX+ HSCs
To determine whether Sp7 may be activated at early stages of hematopoiesis as a potential explanation for the heterogeneity of TdTOSX+;CD45+ cells, we performed flow cytometric analysis for hematopoietic stem cell (HSC) markers in the bone marrow of doxy-fed Osx-cre;TdT and TAM-induced Osx-creERT2;TdT mice bearing B16-F10 tumors subcutaneously. We analyzed the lineage negative (Lin-) LSK and LK populations (Figure 7A) gated based on the expression of Sca1 and c-kit and found 10.9–11.45% TdTOSX+ LSK and 7.08–8.94% TdTOSX+ LK in Osx-cre;TdT and Osx-creERT2;TdT mice, respectively (Figure 7B and D). The LSK population was further divided into HSC and multipotent progenitor MPP(1-4) subsets (Figure 7A and gate strategy in Figure 7—figure supplement 1). Interestingly, TdTOSX+ cells represented about 30% of MPP1, and between 7–11% of the MPP2-MPP3 and MPP4 subsets in Osx-cre;TdT mice (Figure 7C). Similar results were also obtained in Osx-creERT2;TdT mice (Figure 7E). The slight difference in frequency between the two mouse models could depend on the different timing of Cre activation, but both models succeeded to find TdTOSX+ cells in HSCs.
Importantly, RT-PCR of isolated HSCs from the bone marrow of 7–8 week old wild-type naïve mice showed Sp7 expression compared to the TdTOSX- fraction isolated from the tumor mass, used as negative control (Figure 8A). Similarly, Sp7 transcripts were found in the MPP1-4 subsets, and in the common lymphoid progenitors CLP, but not in the common myeloid progenitors CMP (Figure 8A). In contrast, RT-PCR analysis of mature immune populations from naïve mice (macrophages, dendritic cells and T cells) failed to detect Sp7 expression (Figure 8B). Finally, we performed immunostaining for Osx to confirm protein expression on a single cell level. To enrich for Osx+ cells, we FACS sorted TdTOSX+ HSCs from the bone marrow of Osx-cre;TdT reporter mice after enrichment for c-kit+ cells. Sorted cells were plated on serum-coated coverslips for two days before immunostaining for Osx and with DAPI for nuclear localization. We confirmed Osx localization in the nuclei of the majority of TdTOSX+ HSCs (Figure 8D). As positive control for Osx staining, we used bone marrow stromal cells isolated from WT mice cultured for 4 days in osteogenic medium (Figure 8C). Whole bone marrow was used as additional control confirming that only a small number of cells expressed Osx while the majority was negative (Figure 8E). Thus, for the first time we show that subsets of HSCs express the transcription factor Sp7, giving rise to tumor infiltrating immune populations.
Discussion
Numerous studies have established the importance of the TME for primary tumor growth and metastatic dissemination (Goubran et al., 2014). In skeletal metastases, crosstalk between bone residing cells, osteoblasts and osteoclasts, and tumor cells drives tumor growth (Weilbaecher et al., 2011). Recent studies also suggest that osteoblasts and their secretory products stimulate the expansion of bone marrow-derived myeloid populations, which in turn escape the bone marrow and reach distant sites to support tumor growth by inhibiting anti-tumor immune responses (D'Amico et al., 2016; Engblom et al., 2017). Our study demonstrates that a transcription factor required for osteoblast differentiation, Osterix, marks pro-tumorigenic stromal cells infiltrating extra-skeletal tumors. A possible role of Osx in tumorigenesis is emerging from human studies as well showing that presence of Sp7 in the tumor is associated with poor patient survival (Yao et al., 2019). Our own analysis of published gene microarray data of breast cancer stroma (Finak et al., 2008), reveals higher Sp7 expression in tumor stroma relative to healthy tissue from the same subject, corroborating the notion that Osx+ stromal cells regulate tumor growth. Therefore, Osx function may well extend beyond its known role in bone development and homeostasis.
Osterix has been primarily studied in the context of osteoblast differentiation and regulation of bone mass. Sp7 deficient mice die within 1 hr of birth with a complete absence of intramembranous and endochondral bone formation (Nakashima et al., 2002; Baek et al., 2009). During embryogenesis and perinatally, Sp7 is expressed in MSCs in the bone marrow, and it is necessary for the full osteogenic program. In adult mice, Sp7 is mainly confined to committed osteoblasts and osteocytes. There is increasing evidence of Sp7 expression in cells residing outside the skeleton, including synovial fibroblasts (Miura et al., 2019), dental pulp (Monterubbianesi et al., 2019), olfactory glomerular cells, gastric and intestinal epithelium (Chen et al., 2014), and kidney (Strecker et al., 2013). Here we demonstrate that Sp7 is also present in a CAF subset, and in hematopoietic precursors, thus challenging the dogma that Sp7 is exclusively expressed by mesenchymal cells. Our cell tracking studies using the repressible Tet-OFF Osx-cre;TdT reporter and the tamoxifen-inducible Osx-creERT2;TdT model demonstrate the presence of TdTOSX+;CD45+ cells in the bone marrow of adult mice. At least two other independent studies have also shown abundant and persistent Osx+ cells in the bone marrow when the reporter gene is activated either embryonically or neonatally (Liu et al., 2013; Mizoguchi et al., 2014), and at 3 weeks of age (Strecker et al., 2013), but no signal in the bone marrow when the Osx-cre is activated after 8 weeks of age (Mizoguchi et al., 2014). Notably, most of the studies with Osx-cre reporter mice have excluded the hematopoietic marker CD45 in the analysis of the Osx+ populations, and no report to date has shown Sp7 expression in the hematopoietic compartment. Importantly, we also report that bone marrow residing TdTOSX+ populations expand during tumor progression. Such result is highly unlikely linked to systemic factors released from extra-skeletal tumors enhancing the Tet-transactivator in the Tet-OFF;Osx-cre;TdT mice. In fact, we observed increased TdTOSX+ cells (CD45- and CD45+ subsets) also in the bone marrow of tumor bearing Osx-creERT2;TdT animals. Our findings are instead consistent with previous observations reporting increased numbers of bone marrow residing MSCs and osteoblast precursors in mice and patients with lung adenocarcinomas (Engblom et al., 2017), and with the reprogramming of hematopoiesis towards increased myelopoiesis during tumor progression (Capietto et al., 2013; Meyer et al., 2018).
Here, we report for the first time that a subset of hematopoietic lineage cells in the bone marrow and at tumor site derives from an Osx+ progenitor, which is present embryonically and persists until 7–8 weeks of age. Importantly, we detected Sp7 transcripts and Osx nuclear localization in a subset of HSCs. These results support the observation that TdTOSX+;CD45+ cells are very heterogeneous and express monocyte, macrophage, granulocyte, T cell and NK cell markers. The heterogeneity of the TdTOSX+;CD45+ populations is most likely derived from expression of Sp7 in HSC and MPP subsets. These precursors can give rise to the lineage committed common lymphoid (CLP) and myeloid (CMP) progenitors that differentiate into all the mature immune cells (Pietras et al., 2015). Sp7 transcripts are also detected in CLP but not in CMP or mature immune populations, indicating that Osx is activated early during hematopoiesis, with higher expression in cells committed to the lymphoid lineage. Despite expressing both myeloid and lymphoid markers, we noted increased in lymphoid over myeloid ratio in the TdTOSX+;CD45+ tumor infiltrating cells relative to the TdTOSX-;CD45+ cells. TdTOSX+;CD45+ subset is also enriched in genes expressed by immunosuppressive or exhausted lymphocytes. Such result is consistent with the functional tumor co-injection studies showing that TdTOSX-;CD45+ cells reduce tumor growth while the TdTOSX+;CD45+ subset does not.
Co-injection of tumor cells with cells sorted from the primary tumor has been used to validate the pro-tumorigenic roles of CAFs and certain immune suppressive CD45+ populations. Here, we show that tumor-derived TdTOSX+ cells, comprising both CD45 positive and negative subsets, increase tumor growth compared to tumor cells injected alone. By contrast, co-injection of bone marrow-derived TdTOSX+ cells with the tumor cells does not exert any pro-tumor effect suggesting that bone marrow-derived TdTOSX+ cells are functionally distinct from tumor-derived TdTOSX+ cells. Hence, if tumor-derived TdTOSX+ cells originate from the bone marrow in response to an incipient tumor, as it would be anticipated by their increase in the circulation of tumor-bearing mice, they must undergo a conditioning process within the TME to acquire a pro-tumorigenic function. Because Osx regulates expression of proteins required for the generation of the bone tissue, it is likely that a main role of Osx+ mesenchymal cells (TdTOSX+;CD45-) in the TME is to produce and remodel extracellular matrix. Fibroblasts with an osteoblast signature involved in matrix production have been recently identified by single cell RNA sequencing in rheumatoid arthritis (Croft et al., 2019). ECM accumulation is quite frequent within the TME, causing in more severe cases an intense fibrotic response, or desmoplasia, and tumor stiffening (Gkretsi and Stylianopoulos, 2018). Stiffening is not only required for a primary tumor to displace the host tissue and grow in size, but also contributes to cell–ECM interactions and can promote cancer cell invasion to surrounding tissues (Gkretsi and Stylianopoulos, 2018). As noted earlier, upregulation of Sp7 in breast cancer cells is associated with increased invasion and bone metastasis, and this may occur by upregulation of ECM modifying metalloproteinases, MMP9 and MMP13, and other factors that increase vascularization, and affect bone cell function (Yao et al., 2019).
Since Sp7 is not expressed in mature immune populations, it is difficult to envision a role for Osx in modulating anti-tumor immune responses. Co-injection experiments indicate that TdTOSX+;CD45+ cells have better pro-tumorigenic function than TdTOSX–;CD45+ cells. However, only the co-injection of the bulk TdTOSX+ cells increases tumor growth over the B16-F10 tumor cells alone, suggesting that the Osx+ mesenchymal population (TdTOSX+;CD45-) is responsible for enhancing tumor growth. This finding is intriguing since the TdTOSX+;CD45- mesenchymal cells represent only 5% of the total tumor infiltrating TdTOSX+ populations. One could speculate that TdTOSX+;CD45- cells represent a subset of highly pro-tumorigenic CAFs, which can support tumor progression even when present in very limited numbers. Another possibility is that the TdTOSX+;CD45- mesenchymal cells might be required to create and maintain an immune suppressive environment where the CD45+ populations fail to exert anti-tumor effects.
In conclusion, we demonstrate that Sp7 is expressed in a subset of tumor infiltrating mesenchymal cells with CAF and osteogenic cell features. Surprisingly, Sp7 expression is also found in hematopoietic precursors, and marks tumor infiltrating immune populations enriched in immune suppressive markers. Considering the emerging data that Sp7 expression in the tumor cells is linked to tumor progression, our results emphasize the importance of Osx in the TME and the need to evaluate the prognostic value of stromal Sp7 expression in the patients.
Materials and methods
Tumor cell lines, CAF isolation and generation of primary cell cultures
Request a detailed protocolB16-F10 (C57BL/6 mouse melanoma cell line, ATCC # CRL-6475-LUC2), PyMT and PyMT-BO1-GFP (Su et al., 2016) (C57BL/6 mouse breast cancer cell lines), kindly provided by Prof. Katherine N Weilbaecher (Washington University in St. Louis), were cultured at 37°C in complete media (DMEM supplemented with 2 mM l-glutamine, 100 μg/ml streptomycin, 100 IU/ml penicillin, and 1 mM sodium pyruvate) containing 10% FBS. Immortalized CAF cell line was kindly provided by Prof. Gregory D Longmore (Washington University in St. Louis). All the cell lines used have been tested for mycoplasma and were mycoplasma free. CAFs were isolated from primary tumors as described in Corsa et al., 2016. Tumors were minced, digested and plated as single cell suspension on a tissue plastic dish for 30 min to separate the adherent fraction, composed of CAFs, from the cells in suspension, comprising of tumor cells, tumor infiltrating immune populations and endothelial cells, among few others. CAFs were then left 24 hr in culture and RNA was extracted afterwards.
Bone marrow stromal cells (BMSCs) were cultured in complete alpha-MEM (without ascorbic acid) containing 10% FBS and differentiated for 4 days in osteogenic medium (complete alpha-MEM supplemented with 50 μg/ml ascorbic acid and 10 μM beta-glycerophosphate) to obtain pre-osteoblasts.
Sorted hematopoietic stem cells (HSCs) were plated overnight on a glass coverslip coated with serum in the presence of StemPRO (Gibco) media and then fixed and used for immunostaining.
Mouse strains and tumor models
Request a detailed protocolAnimals were housed in a pathogen-free animal facility at Washington University (St. Louis, MO, USA). Age and sex-matched animals were used in all experiments according to IACUC guidelines. B6.Cg-Tg(Sp7-tTA,tetO-EGFP/cre)1Amc/J (catalog #006361; The Jackson laboratory, ME USA) mice, which carry a tetracycline-responsive Osx promoter driving Cre (Osx-cre) were mated with reporter mice B6.Cg-Gt(ROSA)26SorTm9(CAG-tdTomato)Hze/J (TdT) (catalog #007909; The Jackson Laboratory), to generate Osx-cre;TdT mice. To suppress Cre expression, 200ppm doxycyline (doxy) was added to the chow (Test Diet #1816332–203, Purina, MO USA) and fed to some groups of mice until weaning (P28); pups were switched to standard rodent chow at weaning. Sp7-creERT2/Rosa26<fs-TdTomato> mice were kindly provided by Prof. Matthew D Silva (Washington University in St. Louis). Tamoxifen (#T5648, Millipore Sigma, MO USA) was dissolved in corn oil and injected at 100 mg/kg intra-peritoneally. Wild-Type (WT) C57BL/6 mice were purchased from The Jackson Laboratory.
Subcutaneous (sq) injections were performed using 105 B16-F10 tumor cells suspended in 100 µl of sterile PBS and Matrigel Matrix (Corning #354234), while 105 PyMT were suspended in 50 µl of sterile PBS and Matrigel Matrix and injected into the MFP. Adoptive transfer and co-injection with tumor cells were performed using a 5:1 ratio of TdTOSX+:B16-F10 tumor cells. The number of cells injected in each experiment varied depending upon the number of cells obtained from sorting in each experiment (hence, the differences in tumor size among the different experiments), but within (2.5×105:5×104) and (5×105:1×105) cells. B16-F10 cells only were injected as control using the same number of tumor cells co-injected with TdTOSX+ cells. Tumors were monitored by caliper measurements and mice were sacrificed between day 14 and 16 post tumor inoculation. Tumor volume was calculated according to the formula: 0.5236*length*(width [Kalluri and Zeisberg, 2006]).
Flow cytometry and sorting
Request a detailed protocolSingle cell suspensions were prepared from fresh bone marrow and tumors upon sacrifice. Bone marrow cells were obtained from femurs and tibias by centrifugation, while tumors were minced and digested with 3.0 mg/ml collagenase A (Roche, Basel Switzerland) and 50 U/ml DNase I (Millipore Sigma) in serum-free media for 45 min at 37°C. Cells were filtered through 70 μm strainers, red blood cells (RBC) were then removed with RBC lysis buffer (Millipore Sigma), washed twice in PBS and stained in FACS buffer (0.5% BSA, 2 mM EDTA, 0.01%NaN3) with the following antibodies: CD45 (clone 30-F11, eBioscience), CD11b (clone M1/70, eBioscience), Ly6G (Gr-1) (clone RB6-8C5, eBioscience), F4/80 (BM8, Biolegend), CD3 (145–2 C11, eBioscience), CD4 (RM4-5, BD Biosciences), CD8 (clone 53–6.7, BD Biosciences), NK1.1 (clone PK136, Biolegend). Once stained, cells were sorted for the adoptive transfer, RT-PCR and RNAseq studies or were fixed for flow cytometry (gate strategy in Figure 7—figure supplement 1) in BD Cytofix fixation buffer (BD Biosciences, 554655). Samples were acquired using the BD LSR-Fortessa cytometer and analyzed with FlowJo version 9.3.2.
Hematopoietic stem cell subsets were analyzed from fresh bone marrow isolated from femurs, tibias and hip bones crushed with mortar and pestle in 3 ml of PBS, filtered in a 70μm-nylon strainer to remove bone fragments and lysed of red blood cells (RBC). 10 million cells were stained for each sample using the following antibodies: CD3, Gr1, B220 (clone RA3-6B2, eBioscience), Ter119 (clone TER119, e Bioscience), CD41 (clone eBioMWReg30, eBioscience) for lineage negative gating, Sca-1 (clone D7, BD Biosciences), CD117/c-kit (clone ACK2, eBioscience), CD16/32 (clone 2.4G2, BD Biosciences), CD34 (clone RAM34, eBioscience), CD150/SLAM (clone TC15-12F12.2, Biolegend), CD48 (clone HM48-1, eBioscience), CD135/Flt3 (clone A2F10, eBioscience), gate strategy in Figure 7—figure supplement 1. Once stained, cells were fixed in BD Cytofix fixation buffer, acquired using the Bio-Rad ZE5 (YETI) cytometer and analyzed with FlowJo version 10.4.1.
Isolation of hematopoietic progenitor subsets for Real Time analysis was obtained by FACS sorting using the following markers: HSCs (Lineage- c-Kit+Sca-1+CD48-CD150+), MPPs (Lineage- c-Kit+Sca-1+CD48-CD150-), CLPs (Lineage- c-Kit+Sca-1+IL7rα+Flt3+Ly6D-) and CMPs (Lineage- c-Kit+Sca-1-CD34+CD16/32). For IF, TdTomato was also used to sort Osx+ HSCs.
Isolation of hematopoietic stem cells for Osx immunostaining was performed as following: whole bone marrow was isolated from femurs, tibias and iliac crests and incubated with mouse CD117-conjugated microbeads (Miltenyi Biotec #130-091-224). CD117-enrichment was then carried out using the AutoMACS Pro Seperator (Miltenyi Biotec). Post CD117 enrichment, the positive cell fraction was stained to sort phenotypically defined TdTomato+ HSCs (Lineage- c-Kit+Sca-1+CD48-CD150+) by FACS.
Histology
Request a detailed protocolTumors and bones (femur and tibia) were fixed in 4% PFA over-night (ON) at 4°C. The fixed bones were partially decalcified in 14% EDTA for 3 days, washed in a sucrose gradient (1 hr in 10% sucrose, 1 hr in 20% sucrose, ON in 30% sucrose) before snap-freezing them in OCT embedding medium. Tumors were directly put in sucrose gradient and then embedded in OCT. Frozen sections were cut at 5 μm thickness and kept at −20°C until analysis.
Thawed sections were washed in PBS and mounted with VECTASHIELD Mounting Medium containing DAPI (Vector Laboratories, CA USA).
Images were taken using the Leica DMi8 Confocal Microscopy at the Musculoskeletal Research Center (Washington University in St Louis), magnification 200X.
Immunohistochemical staining
Request a detailed protocolTissues were fixed in 10% neutral-buffered formalin for 18 hr, embedded in paraffin after graded-ethanol dehydration, and sectioned into 6 μm sections using a microtome. Automated staining was carried out on the BondRxm (Leica Biosystems). Following dewaxing and citrate-buffered antigen retrieval, sections were stained with the primary antibody for 1 hr at RT. Sequential chromogenic detection was performed with the Bond IntenseR (Rat primary) and Bond Polymer Refine Red (Rabbit primary) detection kits (Leica Biosystems). The primary antibodies used were rat anti-CD45 (Invitrogen clone 30-F11, #14-0451-82) followed by rabbit anti-RFP/TdT (Rockland, # 600-401-379). Stained sections were dehydrated in graduated ethanol and xylene washes then mounted with xylene-based Cytoseal (Thermo Fisher).
Stained slides were imaged at 20X magnification with the Zeiss Axioscan slide scanner. For image analysis, Halo v3 (Indica Labs) was used to deconvolve dual stained images into single channels and pseudo colored with blue representing nuclei (hematoxylin), green representing CD45+ cells (DAB stained), and red representing Osx+ cells (RFP stained).
Immunohistochemestry and immunofluorescence staining
For colocalization of TdT+ cells and CD45 in B16-F10 soft tissue tumors, subcutaneous tumors were dissected and fixed in 10% neutral-buffered formalin for 18 hr, embedded in paraffin after graded-ethanol dehydration, and sectioned into 6 μm sections using a microtome. Sections were dewaxed in xylene and then hydrated through graded ethanol washes. Endogenous peroxidases were quenched by incubating in hydrogen peroxide (1% in PBS) for 10 min at RT. Antigen retrieval was then performed by microwave treatment in a citrate buffer for 15 min. Sections were washed in TBS-T (1X TBS with 0.05%Tween-20) followed by blocking for 30 min at RT in blocking buffer (5% goat serum, 2.5% BSA in TBS). Slides were then blocked using the Avidin/Biotin Blocking Kit (Vector Labs). For staining, slides were incubated overnight in a humidified chamber at 4°C with 1:200 rat anti-CD45 (Invitrogen clone 30-F11, #14-0451-82) and 1:400 rabbit anti-RFP (Rockland, #600-401-379) antibodies diluted in 50% blocking buffer. After primary staining, slides were washed in TBS-T and then incubated with HRP-conjugated anti-rat IgG secondary (1:500 in 50% blocking buffer) for 30 min at RT. For fluorescent detection of CD45, slides were washed and then incubated in 1:50 FITC Tyramide reagent (PerkinElmer, # NEL701A001KT) for 8 min at RT. For fluorescent detection of RFP, slides were washed and incubated in AlexaFluor 594-conjugated anti-rabbit IgG secondary (1:500 in 50% blocking buffer) for 30 min at RT. Stained slides were washed and mounted using ProLong Gold with DAPI (Invitrogen).
Osx staining
Request a detailed protocolHSCs, BMSC or whole bone marrow cells were plated overnight on coverslips coated with FBS for 30 min at RT then cells. Cells were fixed in 4% PFA for 15 min, washed and permeabilized with PBS/0.3%TritonX-100 for 3 min. Blocking was performed with 5% normal serum in permeabilization buffer for 1 hr at RT. Anti-Osx antibody (Abcam, #ab227820) was dissolved 1:500 in a buffer containing PBS/0.3%TritonX-100/1%BSA and incubated ON at 4°C. For fluorescent detection, coverslips were washed and incubated in AlexaFluor 488-conjugated anti-rabbit secondary antibody (1:1000) in PBS/0.3%TritonX-100/1%BSA for 1 hr at RT. Stained slides were washed and mounted using VECTASHIELD Mounting Medium containing DAPI (Vector Laboratories, CA USA). Fluorescent signals were captured by using a Nikon Eclipse 80i microscope and a Nikon DS-Qi1MC camera (Nikon, CO, USA).
Real Time PCR
Request a detailed protocolTotal RNA was extracted with TRIzol (Invitrogen, CA USA) and quantified on a ND-1000 spectrophotometer (NanoDrop Technologies). The cDNA was synthesized with 1 μg RNA using High Capacity cDNA Reverse Transcription Kit (#4368814, Applied Biosystems, CA USA).
For purified hematopoietic precursors (HSC, MPP, CLP and CMP) total RNA was extracted using the NucleoSpin RNA XS kit (Macherey-Nagel #740902) and reverse transcribed with the SuperScript VILO kit (Invitrogen #11754–050).
The amount of each gene was determined using Power SYBR Green mix on 7300 Real-Time PCR System (Applied Biosystems). Cyclophilin mRNA was used as housekeeping control. Specific primers for mice were as follows: Cyclophilin, 5’-AGC ATA CAG GTC CTG GCA TC-3’ and 5’-TTC ACC TTC CCA AAG ACC AC-3’; Osterix, 5’-AAG GGT GGG TAG TCA TTT GCA-3’ and 5’-CCC TTC TCA AGC ACC AAT GG-3’; S100a4 (Fsp1), 5’-TGA GCA ACT TGG ACA GCA ACA-3’ and 5’-TTC CGG GGT TCC TTA TCT GGG-3’; Acta2 (a-SMA), 5’-GTC CCA GAC ATC AGG GAG TAA-3’ and 5’-TCG GAT ACT TCA GCG TCA GGA-3’; Col1a1, 5’-GGC CTT GGA GGA AAC TTT GC-3’ and 5’-GGG ACC CAT TGG ACC TGA AC-3’; Bglap (Osteocalcin), 5’-GGA CTG AGG CTC TGT GAG GT-3’ and 5’-CAG ACA CCA TGA GGA CCA TC-3’; Runx2, 5’-GTT ATG AAA AAC CAA GTA GCC AGG-3’ and 5’-GTA ATC TGA CTC TGT CCT TGT GGA-3’; Ibsp (Bsp), 5’-AGG ACT AGG GGT CAA ACA C-3’ and 5’-AGT AGC GTG GCC GGT ACT TA-3’; Alpl (Tnap), 5’-GGG GAC ATG CAG TAT GAG TT-3’ and 5’-GGC CTG GTA GTT GTT GTG AG-3’. The relative quantification in gene expression was determined using the 2-ΔΔCt method.
RNA sequencing and analysis
Request a detailed protocolSingle cell suspensions were prepared from fresh tumor upon sacrifice. Tumors were minced and digested as described above. Tumor stromal cells were separated from the GFP+ tumor cells and sorted based on the expression, or lack of thereof, of TdT and CD45 markers.
Library preparation was performed with 10 ng of total RNA and integrity was determined using an Agilent bioanalyzer. ds-cDNA was prepared using the SMARTer Ultra Low RNA kit for Illumina Sequencing (Clontech) per manufacturer’s protocol. The cDNA was fragmented using a Covaris E220 sonicator using peak incident power 18, duty factor 20%, cycles/burst 50, for a 120 s. cDNA was blunt ended, had an A base added to the 3’ ends, and then had Illumina sequencing adapters ligated to the ends. Ligated fragments were then amplified for 12 cycles using primers incorporating unique index tags. Fragments were sequenced on an Illumina HiSeq 3000 using single end reads extending 50 bases. Basecalls and demultiplexing were performed with Illumina’s bcl2fastq software and a custom python demultiplexing program with a maximum of one mismatch in the indexing read. RNA-seq reads were then aligned to the Mus musculus Ensembl release 76 top-level assembly with STAR version 2.0.4b (Dobin et al., 2013). Gene counts were derived from the number of uniquely aligned unambiguous reads by Subread:featureCount version 1.4.5 (Liao et al., 2014). Isoform expression of known Ensembl transcripts were estimated with Sailfish version 0.6.13 (Patro et al., 2017). Sequencing performance was assessed with RSeQC version 2.3 (Wang et al., 2012).
All gene counts were then imported into the R/Bioconductor package EdgeR (Robinson et al., 2010) and TMM normalization size factors were calculated to adjust for samples for differences in library size. Ribosomal genes and genes not expressed in at least two samples greater than one count-per-million were excluded from further analysis. The TMM size factors and the matrix of counts were then imported into the R/Bioconductor package Limma (Ritchie et al., 2015) and weighted likelihoods based on the observed mean-variance relationship of every gene and sample were then calculated with Limma’s voomWithQualityWeights (Liu et al., 2015). The data were then fitted to a Limma generalized linear model to test for changes in a single cell population relative to the mean of all other populations to find only genes that were uniquely up-regulated with Benjamini-Hochberg false-discovery rate adjusted p-values less than or equal to 0.05. For each group of cells, global perturbations in known Gene Ontology (GO) terms were then measured using the R/Bioconductor package GAGE (Luo et al., 2009) to quantify the mean log two fold-changes of ll genes in a given term versus the background log two fold-changes of all genes found outside the respective term. Only globally up-regulated GO terms with Benjamini-Hochberg false-discovery rate adjusted p-values less than or equal to 0.05 were considered for comparison across populations.
Data availability
Sequencing data are deposited in GEO under accession code GSE143586 https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE143586.
-
NCBI Gene Expression OmnibusID GSE143586. GSE143586.
References
-
Peripheral blood fibrocytes: differentiation pathway and migration to wound sitesThe Journal of Immunology 166:7556–7562.https://doi.org/10.4049/jimmunol.166.12.7556
-
Positive regulation of adult bone formation by osteoblast-specific transcription factor osterixJournal of Bone and Mineral Research 24:1055–1065.https://doi.org/10.1359/jbmr.081248
-
Down-regulation of PLCγ2-β-catenin pathway promotes activation and expansion of myeloid-derived suppressor cells in CancerThe Journal of Experimental Medicine 210:2257–2271.https://doi.org/10.1084/jem.20130281
-
Turning foes to friends: targeting cancer-associated fibroblastsNature Reviews Drug Discovery 18:99–115.https://doi.org/10.1038/s41573-018-0004-1
-
Dickkopf-related protein 1 (Dkk1) regulates the accumulation and function of myeloid derived suppressor cells in CancerJournal of Experimental Medicine 213:827–840.https://doi.org/10.1084/jem.20150950
-
Fibulin-7 is overexpressed in glioblastomas and modulates glioblastoma neovascularization through interaction with angiopoietin-1International Journal of Cancer 145:2157–2169.https://doi.org/10.1002/ijc.32306
-
STAR: ultrafast universal RNA-seq alignerBioinformatics 29:15–21.https://doi.org/10.1093/bioinformatics/bts635
-
Stromal gene expression predicts clinical outcome in breast CancerNature Medicine 14:518–527.https://doi.org/10.1038/nm1764
-
N-cadherin regulation of bone growth and homeostasis is osteolineage Stage-SpecificJournal of Bone and Mineral Research 32:1332–1342.https://doi.org/10.1002/jbmr.3112
-
CD163L1 and CLEC5A discriminate subsets of human resident and inflammatory macrophages in vivoJournal of Leukocyte Biology 98:453–466.https://doi.org/10.1189/jlb.3HI1114-531R
-
Fibromodulin is expressed by both chondrocytes and osteoblasts during fetal bone developmentJournal of Cellular Biochemistry 82:46–57.https://doi.org/10.1002/jcb.1115
-
Fibrocytes: a novel stromal cells to regulate resistance to Anti-Angiogenic therapy and Cancer progressionInternational Journal of Molecular Sciences 19:98.https://doi.org/10.3390/ijms19010098
-
Regulation of tumor growth and metastasis: the role of tumor microenvironmentCancer Growth and Metastasis 7:9–18.https://doi.org/10.4137/CGM.S11285
-
Collagen XXIV, a vertebrate fibrillar collagen with structural features of invertebrate collagens: selective expression in developing cornea and boneThe Journal of Biological Chemistry 278:43236–43244.https://doi.org/10.1074/jbc.M302112200
-
CREB-AP1 protein complexes regulate transcription of the collagen XXIV gene (Col24a1) in osteoblastsJournal of Biological Chemistry 281:5445–5452.https://doi.org/10.1074/jbc.M509923200
-
Hematopoietic stem cell derived carcinoma-associated fibroblasts: a novel originInternational Journal of Clinical and Experimental Pathology 5:863–873.https://doi.org/10.1016/j.neo.2015.04.004
-
The role of fibromodulin in Cancer pathogenesis: implications for diagnosis and therapyCancer Cell International 19:157.https://doi.org/10.1186/s12935-019-0870-6
-
Bone marrow-derived fibroblasts are a functionally distinct stromal cell population in breast CancerJournal of Experimental Medicine 215:3075–3093.https://doi.org/10.1084/jem.20180818
-
Antagonizing integrin β3 increases immunosuppression in CancerCancer Research 76:3484–3495.https://doi.org/10.1158/0008-5472.CAN-15-2663
-
Circulating fibrocytes prepare the lung for Cancer metastasis by recruiting Ly-6C+ monocytes via CCL2The Journal of Immunology 190:4861–4867.https://doi.org/10.4049/jimmunol.1202857
-
RSeQC: quality control of RNA-seq experimentsBioinformatics 28:2184–2185.https://doi.org/10.1093/bioinformatics/bts356
-
Cancer to bone: a fatal attractionNature Reviews Cancer 11:411–425.https://doi.org/10.1038/nrc3055
Article and author information
Author details
Funding
National Institutes of Health (R01 CA235096)
- Roberta Faccio
National Institutes of Health (R01 AR066551)
- Roberta Faccio
Shriners Hospital (P19-07408CR)
- Roberta Faccio
Siteman Cancer Center
- Roberto Civitelli
- Roberta Faccio
JIT Award by Washinghton University Institute for Clinical and Translational Sciences (TR000448)
- Francesca Fontana
- Roberto Civitelli
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank the Musculoskeletal Histology Core of Musculoskeletal Research Center supported by the National Institutes of Health P30 Grants AR057235 and P30 AR074992, and the Siteman Flow Cytometry Core at Washington University in St Louis. We also thank the Genome Technology Access Center in the Department of Genetics at Washington University School of Medicine for help with genomic analysis. The Center is partially supported by NCI Cancer Center Support Grant #P30 CA91842 to the Siteman Cancer Center and by ICTS/CTSA Grant# UL1TR002345 from the National Center for Research Resources (NCRR), a component of the National Institutes of Health (NIH), and NIH Roadmap for Medical Research. This publication is solely the responsibility of the authors and does not necessarily represent the official view of NCRR or NIH.
Data analysis was performed in part through the use of Washington University Center for Cellular Imaging (WUCCI) supported by Washington University School of Medicine, The Children’s Discovery Institute of Washington University and St. Louis Children’s Hospital (CDI-CORE-2015–505 and CDI-CORE-2019–813) and the Foundation for Barnes-Jewish Hospital (3770 and 4642).
We thank Danielle Ketterer, Sahil Mahajan PhD, Giulia Leanza PhD and Kristann Magee for the technical support in collecting in vivo samples.
This research was supported by grants from National Institutes of Health Grants R01 AR066551 (to RF), R01 CA235096 (to RF), grant from Shriners Hospital P19-07408 CR (to RF) and grants from Siteman Investment Program (to RF, RF and RC).
Partial support to this project was provided by a Just-in-Time Award by the Washington University Institute for Clinical and Translational Sciences (to RC and FF), Washington University (supported by NIH CTSA UL1 TR000448); and by funds by the Barnes-Jewish Hospital Foundation (to RC).
Ethics
Animal experimentation: This study was performed in accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health. All the animals were handled according to the Institutional Animal Care and Use Committee (IACUC) protocols (#2016-0228 and # 2019-0982 to RF, #2014-0279 and #2017-0095 to RC).
Copyright
© 2020, Ricci 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,139
- views
-
- 217
- downloads
-
- 13
- 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
- Evolutionary Biology
In growing cell populations such as tumours, mutations can serve as markers that allow tracking the past evolution from current samples. The genomic analyses of bulk samples and samples from multiple regions have shed light on the evolutionary forces acting on tumours. However, little is known empirically on the spatio-temporal dynamics of tumour evolution. Here, we leverage published data from resected hepatocellular carcinomas, each with several hundred samples taken in two and three dimensions. Using spatial metrics of evolution, we find that tumour cells grow predominantly uniformly within the tumour volume instead of at the surface. We determine how mutations and cells are dispersed throughout the tumour and how cell death contributes to the overall tumour growth. Our methods shed light on the early evolution of tumours in vivo and can be applied to high-resolution data in the emerging field of spatial biology.
-
- Cancer Biology
- Evolutionary Biology
In asexual populations that don’t undergo recombination, such as cancer, deleterious mutations are expected to accrue readily due to genome-wide linkage between mutations. Despite this mutational load of often thousands of deleterious mutations, many tumors thrive. How tumors survive the damaging consequences of this mutational load is not well understood. Here, we investigate the functional consequences of mutational load in 10,295 human tumors by quantifying their phenotypic response through changes in gene expression. Using a generalized linear mixed model (GLMM), we find that high mutational load tumors up-regulate proteostasis machinery related to the mitigation and prevention of protein misfolding. We replicate these expression responses in cancer cell lines and show that the viability in high mutational load cancer cells is strongly dependent on complexes that degrade and refold proteins. This indicates that the upregulation of proteostasis machinery is causally important for high mutational burden tumors and uncovers new therapeutic vulnerabilities.