1. Stem Cells and Regenerative Medicine
Download icon

Distinct skeletal stem cell types orchestrate long bone skeletogenesis

  1. Thomas H Ambrosi
  2. Rahul Sinha
  3. Holly M Steininger
  4. Malachia Y Hoover
  5. Matthew P Murphy
  6. Lauren S Koepke
  7. Yuting Wang
  8. Wan-Jin Lu
  9. Maurizio Morri
  10. Norma F Neff
  11. Irving L Weissman
  12. Michael T Longaker
  13. Charles KF Chan  Is a corresponding author
  1. Institute for Stem Cell Biology and Regenerative Medicine, Stanford University School of Medicine, United States
  2. Chan Zuckerberg BioHub, United States
  3. Ludwig Center for Cancer Stem Cell Biology and Medicine at Stanford University, United States
  4. Department of Surgery, Stanford University School of Medicine, United States
  5. Hagey Laboratory for Pediatric Regenerative Medicine, Stanford University School of Medicine, Stanford University, United States
Research Article
Cite this article as: eLife 2021;10:e66063 doi: 10.7554/eLife.66063
5 figures, 1 table and 1 additional file

Figures

Figure 1 with 1 supplement
Two cell populations with skeletal stem cell characteristics in postnatal long bones.

(A) Diagram showing two previously described skeletal stem cell (SSC) populations and the downstream populations they generate, which were defined by the specific expression patterns of cell surface proteins. Top: SSC lineage tree of the osteochondral SSC (ocSSC) that gives rise to bone, cartilage, and stromal populations. Bottom: the lineage tree of the perivascular SSC (pvSSC) able to give rise to bone, cartilage, adipose tissue, and stromal populations. BCSP: bone cartilage stroma progenitor; CP: cartilage progenitor; APC: adipogenic progenitor cell; Pre-Ad.: pre-adipocyte; OPC: osteochondrogenic progenitor cell. (B) Representative flow cytometric gating strategy for the isolation of ocSSCs (top) and pvSSCs (bottom). (C) Representative confocal microscopy images of in vivo derived single-color clonal colonies of renal capsule-transplanted purified ocSSCs (left) and pvSSCs (right) derived from Actin-CreERt Rainbow mice. Three independent transplants per cell type under renal capsules were performed. (D) Renal capsule transplant-derived ossicles of purified GFP-labeled ocSSCs (top) and pvSSCs (bottom). Images show the photograph of the kidney with transplant (top left) and GFP signal of graft tissue (bottom left) as well as Movat pentachrome cross-section staining (right). (E) Quantification of ocSSC (top) and pvSSC (bottom) graft composition. Results of three separate experiments with n = 3 per SSC type. All data are shown as mean ± SEM. (F) Representative images of staining of clonally derived cultures that underwent tri-lineage differentiation assays in vitro. Alcian blue (chondrogenesis), Alizarin red S (osteogenesis), and oil red O (adipogenesis) stainings are shown along with the number of clones that stained positive for each differentiation type (ocSSC n = 7; pvSSC n = 8 clones). (G) Representative immunohistochemistry images for GFP (green) and perilipin (red) of tissue derived from intratibially transplanted purified GFP-labeled ocSSCs (top) and pvSSCs (bottom) 2 weeks after injection. White arrowheads: GFP+Perilipin+ cells. TB: trabecular bone. Three separate experiments with ocSSC n = 3 and pvSSC n = 4. Scale bars, 30 µm.

Figure 1—figure supplement 1
Perivascular SSCs, but not ocSSCs, are a source of bone marrow adipose tissue.

(A) Clonally derived osteochondrogenic (OC) and adipogenic (Ad) cell types of transplanted purified osteochondrogenic skeletal stem cells (ocSSCs) (top) and perivascular SSCs (pvSSCs) (bottom) isolated from ‘Rainbow’ mice detected in derived renal grafts. Left, Movat pentachrome staining of graft tissue and right, adjacent section with immunofluorescence of a green and a red clone. (B) Flow cytometric profile of GFP-labeled ocSSC (top) or pvSSC (bottom)-derived graft tissue of renal transplants for expression of CD45 and GFP. (C) Representative flow cytometric profiles of in vitro (left) derived primary cultures of purified ocSSCs (top) and pvSSCs (bottom) as well as in vivo generated grafts (right). (D) Fibroblast colony-forming unit (CFU-F) assay results of purified ocSSCs and pvSSCs from long bones of E17.5 embryos and 8-week-old mice (n = 3 per age group). Representative colonies derived from ocSSCs and pvSSCs stained for crystal violet are shown above respective data columns. (E) Tri-lineage differentiation outcome of bulk-sorted ocSSCs (left) and pvSSCs (right) from E17.5- or 8-week-old-mouse-derived long bones. Alcian blue (chondrogenesis), Alizarin red S (osteogenesis), and oil red O (adipogenesis) stainings are shown. APCs (adipogenic progenitor cells) that underwent chondrogenesis are shown as control. (F) Oil red O stainings of ocSSC and pvSSC derived from long bones of 24-month-old mice (24-mo) that underwent adipogenic differentiation are shown. (G) Immunohistochemistry fluorescence image of GFP-labeled ocSSC (top)- and pvSSC (bottom)-derived cells 2 weeks after intratibial injection. GP = growth plate. (H) Co-staining of ocSSC (top)- and pvSSC (bottom)-derived cells (green) with the osteogenic label osteocalcin (OCN) and (I) the chondrogenic label COL2. Scale bars, 30 µm. All data are shown as mean + SEM. Significance between groups assessed by unpaired, two-tailed Student’s t-test and corrected with Welch’s test for unequal distribution if needed.

Figure 1—figure supplement 1—source data 1

Functional characterization of ocSSCs and pvSSCs.

https://cdn.elifesciences.org/articles/66063/elife-66063-fig1-figsupp1-data1-v1.xlsx
Figure 2 with 1 supplement
Osteochondral and perivascular SSCs are anatomically and developmentally distinct.

(A) Flow cytometric quantification of micro-dissected long bone regions of 8-week-old mice for the prevalence of osteochondrogenic skeletal stem cells (ocSSCs), bone cartilage and stromal progenitors (BCSPs), perivascular SSCs (pvSSCs), and adipogenic progenitor cells (APCs). EP: epiphysis, MP: metaphysis, DP: diaphysis, PE: periosteum (n = 3 mice). (B) Frequency of ocSSCs, BCSPs (top) and pvSSCs, APCs (bottom) in the long bones of mouse embryos developing into postnatal life assessed by flow cytometry (n = 3 mice per age group). (C) Dissected limb bones of GFP-expressing embryos between E12-E13.5 (≤E13.5) and postnatal day 1–10 (≥P1) were transplanted under the renal capsule of NSG mice. 4 weeks after transplantation, bones were dissected out and sectioned for hematoxylin and eosin staining (left) and immunohistological staining for GFP (green) and perilipin (red) quantified as the percentage of GFP+Perilipin+ (right) cells (n = 4 bones per age group) below the distal growth plate. (D) Representative flow cytometric plots showing Sca1 expression in cells gated for CD45-Ter119-Tie2-CD51+6C3-Thy1- (ocSSC/BCSP) in long bones at postnatal day 1 and at 8 weeks of age (uninjured) as well as during regeneration at days post fracture (dpf) 5 and 15. (E) Flow cytometric analysis of lineage output in ocSSC (n = 4)- and pvSSC (n = 5)-derived renal grafts displayed as the percentage of all donor-derived cells from at least two independent experiments. All data are shown as mean + SEM. Scale bars, 30 µm.

Figure 2—source data 1

Anatomical and developmental assessment of ocSSCs and pvSSCs.

https://cdn.elifesciences.org/articles/66063/elife-66063-fig2-data1-v1.xlsx
Figure 2—figure supplement 1
No overlap between bone-resident SSC subtypes.

(A) Representative flow cytometric plots of single-cell suspensions from 8-week-old long bones dissected into metaphysis/GP (growth plate) and flushed bone marrow stained for osteochondrogenic skeletal stem cells (ocSSCs) (top) and perivascular SSCs (pvSSCs) (bottom). (B) Image of GFP-labeled tibia and thoracic spine transplanted under the renal capsule (left). Image of GFP-labeled postnatal day 1 (P1) and E13.5 tibia bones dissected from the renal capsule of NSG mice 4 weeks after transplantation (right). (C) Flow cytometric analysis of transplanted and dissected out tibia bones for GFP+ocSSCs and bone cartilage and stromal progenitors (BCSPs). (D) Flow cytometric analysis of transplanted and dissected out tibia bones for GFP+Sca1+ cells. (E) Representative flow cytometric plot showing gating for ocSSC lineage tree from CD45-Ter119-Sca1+ cells. (F) Quantification of flow cytometric data of ocSSC lineage tree cell populations as percentage of CD45-Ter119-Sca1+ cells (n = 3 mice). (G) Renal transplants of GFP-labeled E13.5 long bone-derived ocSSCs showing flow cytometric plot (left) and quantification of derived phenotypic pvSSC/APC (right) (n = 3). All data shown as mean + SEM.

Figure 3 with 1 supplement
Molecular differences of SSC types infer specific niche functions.

(A) Single-cell RNA-sequencing analysis results of 143 osteochondrogenic skeletal stem cells (ocSSCs) and 169 perivascular SSCs (pvSSCs) shown as clustering by Uniform Manifold Approximation and Projection (UMAP). (B) Top Gene Ontology (GO) Biological Processes by all differentially expressed genes between ocSSCs and pvSSCs as determined through EnrichR. (C) Track plots of top 10 differentially expressed genes in ocSSCs and pvSSCs showing individual peaks per cell as the degree of their expression. (D) Dot plots showing expression of selected genes in ocSSCs and pvSSCs previously reported to characterize specific cell types. (E) Flow cytometric quantification of ocSSCs and pvSSCs 2 weeks after whole-body sublethal irradiation of 8-week-old mice (n = 6 control; n = 7 5Gy, from two independent experiments). (F) Flow cytometric quantification of ocSSCs and pvSSCs at various days after stabilized bi-cortical femoral fractures of 8-week-old mice shown as fold change of uninjured (n = 3 per timepoint from two independent experiments). (G) Flow cytometric quantification of ocSSCs, bone cartilage and stromal progenitors (BCSPs), pvSSCs, and adipogenic progenitor cells (APCs) in long bones of young (8 weeks; n = 3) and aged (30 months; n = 3) mice shown as fold change of young. (H) Renal capsule-derived grafts of co-transplants of equal numbers of ocSSCs and pvSSCs showing Movat pentachrome-stained cross-section (left) and the corresponding immunohistological staining (right) for GFP (ocSSC-derived cells) and red fluorescent protein (RFP) (pvSSC-derived cells). Light blue arrowhead: cartilage (Ca); yellow arrowhead: bone (Bo); brown arrowhead: marrow (Ma) lining cells. All data are shown as mean + SEM. Significance between groups was assessed by unpaired, two-tailed Student’s t-test and corrected with Welch’s test for unequal distribution if needed. Scale bars, 30 µm.

Figure 3—source data 1

Changes in SSC abundance in response to injury and aging.

https://cdn.elifesciences.org/articles/66063/elife-66063-fig3-data1-v1.xlsx
Figure 3—figure supplement 1
Two genetically and functionally distinct SSC lineages of long bones.

(A) Total read counts (top) and uniquely mapped reads (bottom) of quality-filtered osteochondrogenic skeletal stem cells (ocSSCs) (143 cells) and perivascular SSCs (pvSSCs) (169 cells). (B) ERCC (External RNA Controls Consortium) fraction of all read counts (top) and gene count (bottom) of quality-filtered ocSSCs (143 cells) and pvSSCs (169 cells). (C) Plots showing ribosomal (top) and mitochondrial (bottom) genes as percentage of total gene expression in filtered datasets. (D) Expression of Ly6a/Sca1 in single-cell RNA-sequencing data from ocSSCs and pvSSCs. (E) Flow cytometric quantification of bone cartilage and stromal progenitors (BCSPs) and adipogenic progenitor cells (APCs) 2 weeks after whole-body sublethal irradiation of 8-week-old mice (n = 6 control; n = 7 5Gy, from two independent experiments). (F) Flow cytometric quantification of BCSPs and APCs at various days after stabilized bi-cortical femoral fractures of 8-week-old mice shown as fold change of uninjured (n = 3 per timepoint from two independent experiments). Statistical testing was performed by unpaired, two-tailed Student’s t-test between groups. All data are shown as mean + SEM.

Figure 3—figure supplement 1—source data 1

Quality control of single cell RNA-sequencing data.

https://cdn.elifesciences.org/articles/66063/elife-66063-fig3-figsupp1-data1-v1.xlsx
Variable expression of commonly used reporter genes in skeletal stem cell types.

(A) Microarray data of freshly purified bulk cell populations of the previously defined osteochondrogenic skeletal stem cell (ocSSC) and perivascular SSC (pvSSC) lineage trees from 8-week-old mice showing gene expression of markers commonly known to trace and/or label cell populations enriched for SSCs. Expression is shown as normalized activity as processed by GEXC (Gene Expression Commons). Each cell-type represents a pool of cells derived from three to four mice. (B) Single-cell heatmaps of gene expression of markers commonly used to trace and/or label cell populations enriched for SSCs as observed in single-cell RNA-sequencing (scRNAseq) results of ocSSCs and pvSSCs. (C) Flow cytometric analysis of antibody labeling for leptin receptor (LepR) in different populations of the ocSSC and pvSSC lineage trees in differently aged mice (n = 3). (D) Fibroblast colony-forming unit (CFU-F) assay of freshly isolated ocSSCs and pvSSCs separated by their expression of LepR (ocSSC n = 9; pvSSC n = 6, from two independent experiments). (E) Expression of osteocalcin (OCN, green) in in vitro osteogenically differentiated ocSSCs and pvSSCs in LepR-positive and -negative fractions. (F) Quantification of percentage of cells expressing OCN as determined by antibody labeling at 2 weeks of osteogenic differentiation (n = 5). Significance between groups was assessed by unpaired, two-tailed Student’s t-test. All data are shown as mean + SEM. Scale bars, 30 µm.

Figure 4—source data 1

Differences of LepR expression in SSC subtypes.

https://cdn.elifesciences.org/articles/66063/elife-66063-fig4-data1-v1.xlsx
Figure 5 with 1 supplement
Stem cell crosstalk by distinct SSC subpopulations facilitates niche interactions.

(A) Single-cell RNA-sequencing analysis results of osteochondrogenic skeletal stem cells (ocSSCs) and perivascular SSCs (pvSSCs) shown as Leiden clustering to reveal heterogeneity within cell populations (top). Composition of clusters by SSC type (bottom). Cluster 1: 155 cells; cluster 2: 57 cells; cluster 3: 70 cells; cluster 4: 30 cells. (B) Heatmap of top 40 differentially expressed genes for each of the four Leiden clusters. (C) Violin plots of selected marker genes for each of the four Leiden clusters. (D) Expression of selected ligands and their receptors in Leiden clusters of SSC populations. Connected ligand-receptor genes in pairs have scaled z-score expression >0.5. (E) Expression of cluster 2-specific marker Wif1 in UMAP plot. (F) Expression of cluster 1-specific marker Cdh13 in UMAP plot. (G) Representative in situ RNAscope images showing detection of the Wif1 RNA transcripts in the growth plate (left) and their absence in diaphyseal bone marrow (right). (H) Representative in situ RNAscope images showing detection of the Cdh13 RNA transcripts in diaphyseal bone marrow (right) and their absence in the growth plate (left). Red arrowheads: RNA transcript-expressing cells. RZ: resting zone; PZ: proliferative zone; HZ: hypertrophic zone. Scale bars, 10 µm.

Figure 5—figure supplement 1
Single-cell heterogeneity and crosstalk of SSC subtypes.

(A) Expression of genes previously associated with stem cell-like populations and their expression in the four Leiden clusters. (B) Pathway enrichment analysis of top 100 genes expressed in Leiden clusters as determined by BioPlanet 2019 through EnrichR. (C) Cell cycle analysis of single-cell RNA-sequencing (scRNAseq) data with cell cycle status distribution in different Leiden clusters. (D) CytoTrace scores for cells of each cluster. Data are shown as box plots with min to max distribution (cluster 1: 155 cells; cluster 2: 57 cells; cluster 3: 70 cells; cluster 4: 30 cells). (E) Gene expression of receptors reported to be involved in skeletal regulation in the four Leiden clusters. (F) Representative in situ RNAscope images showing detection of the Wif1 (top) and Cdh13 (bottom) RNA transcripts in periosteal cells. (G) Schematic summary of the presented work. Single-cell suspension of mechanically and enzymatically digested bone tissue with subsequent selection by plastic adherence has long been used as a gold standard to enrich for skeletal stem and progenitor cells. Going beyond single-gene label selection for skeletal stem cell (SSC) populations, this work prospectively isolated SSC subtypes using a combination of surface markers. Functional in vitro and in vivo studies combined with SmartSeq2 scRNAseq revealed the existence of two distinct bona fide SSC types, an osteochondral and a perivascular SSC, with differences with respect to their anatomical localization, developmental appearance, differentiation capacity, and function. Scale bars, 10 µm.

Figure 5—figure supplement 1—source data 1

Cytotrace analysis of leiden clusters.

https://cdn.elifesciences.org/articles/66063/elife-66063-fig5-figsupp1-data1-v1.xlsx

Tables

Key resources table
Reagent type
(species)
or resource
DesignationSource
or reference
IdentifiersAdditional
information
StrainMouse: B6 (C57BL/Ka-Thy1.1-CD45.1)The Jackson LaboratoryJAX: 000406RRID:IMSR_JAX:000406
StrainMouse: NSG (NOD.Cg-Prkdcscid Il2rgtm1Wjl/SzJ)The Jackson LaboratoryJAX: 005557RRID:IMSR_JAX:005557
StrainMouse: GFP (C57BL/6-Tg(CAG-EGFP)1Osb/J)The Jackson LaboratoryJAX: 003291RRID:IMSR_JAX:003291
StrainMouse: Actin-CreERt Rosa26-Rainbow (homozygous)In-houseN/A
AntibodyAnti-Mouse Ly-6A/E (Sca-1) APC (clone: D7)ThermoFisherCat#: 17–5981FACS 1:200 (RRID:AB_469487)
AntibodyAnti-Mouse Ly-6A/E (Sca-1) Alexa Fluor 700 (clone: D7)ThermoFisherCat#: 56–5981FACS 1:200 (RRID:AB_657837)
AntibodyAnti-Mouse Ly-6A/E (Sca-1) PE/Cy7 (clone: D7)ThermoFisherCat#: 25–5981FACS 1:200 (RRID:AB_469669)
AntibodyAnti-Mouse CD45 FITC (clone: 30-F11)ThermoFisherCat#: 11–0451FACS 1:200 (RRID:AB_465050)
AntibodyAnti-Mouse CD45 APC (clone: 30-F11)ThermoFisherCat#: 17–0451FACS 1:200 (RRID:AB_469393)
AntibodyAnti-Mouse CD45 PE/Cy5 (clone: 30-F11)ThermoFisherCat#: 15–0451FACS 1:200 (RRID:AB_468751)
AntibodyAnti-Mouse CD31 (PECAM-1) FITC (clone: 390)ThermoFisherCat#: 11–0311FACS 1:200 (RRID:AB_465011)
AntibodyAnti-Mouse CD31 (PECAM-1) APC (clone: 390)ThermoFisherCat#: 17–0311FACS 1:200 (RRID:AB_657735)
AntibodyAnti-Mouse TER-119 PE/Cy5 (clone: TER-119)ThermoFisherCat#: 15–5921FACS 1:200 (RRID:AB_468811)
AntibodyAnti-Mouse CD140a (PDGFRA) APC (clone: APA5)ThermoFisherCat#: 17–1401FACS 1:100 (RRID:AB_529482)
AntibodyAnti-Mouse CD24 APC-eFluor 780 (clone: M1/69)ThermoFisherCat#: 47–0242FACS 1:200 (RRID:AB_10853172)
AntibodyAnti-Mouse CD51 PE (clone: RMV7)ThermoFisherCat#: 12–0512FACS 1:100 (RRID:AB_465703)
AntibodyAnti-Mouse CD90.1 APC-eFluor 780 (clone: HIS51)ThermoFisherCat#: 47–0900FACS 1:100 (RRID:AB_1272256)
AntibodyAnti-Mouse CD90.2 APC-eFluor 780 (clone: 53–2.1)ThermoFisherCat#: 47–0902FACS 1:100 (RRID:AB_1272187)
AntibodyAnti-Mouse BP1 APC (clone: 6C3)ThermoFisherCat#: 17–5891FACS 1:100 (RRID:AB_2762697)
AntibodyAnti-Mouse CD105 (Endoglin) Biotin (clone: MJ7/18)ThermoFisherCat#: 13–1051FACS 1:100 (RRID:AB_466555)
AntibodyAnti-Mouse Tie2 (clone: Tek4)ThermoFisherCat#: 14–5987FACS 1:20 (RRID:AB_467792)
AntibodyAnti-Mouse Leptin R Biotinylated Antibody (goat polyclonal)R and D SystemsCat#: BAF497FACS 1:50 (RRID:AB_2296953)
AntibodyGoat anti-GFP (polyclonal)NovusBiologicalsCat#: NB100-1770IF 1:200 (RRID:AB_10128178)
AntibodyRabbit anti-Perilipin (polyclonal)ThermoFisherCat#: PA5-72921IF 1:200 (RRID:AB_2718775)
AntibodyRabbit anti-Osteocalcin (polyclonal)AbcamCat#: ab93876IF 1:200 (RRID:AB_10675660)
AntibodyRabbit anti-Collagen II (polyclonal)AbcamCat#: ab34712IF 1:200 (RRID:AB_731688)
AntibodyRat anti-Endomucin (Clone: V.7C7)ThermoFisherCar#: 14-5851-82IF 1:400 (RRID:AB_891527)
AntibodyAlexa Fluor 488 donkey anti-goat (polyclonal)ThermoFisherCat#: A32814IF 1:500 (RRID:AB_2762838)
AntibodyAlexa Fluor 594 donkey anti-rabbit (polyclonal)ThermoFisherCat#: A21207IF 1:500 (RRID:AB_141637)
AntibodyAlexa Fluor 700 goat anti-rabbit (polyclonal)ThermoFisherCat#: A21038IF 1:500 (RRID:AB_2535709)
AntibodyAlexa Flour 488 goat anti-rat (polyclonal)AbcamCat#: ab150157IF 1:500 (RRID:AB_2722511)
Sequence-based reagentOligo-dT30VNIDTN/A(5’-AAGCAGTGGTATCAACGCA
GAGTACT30VN-3’)
Sequence-based reagentTemplate-switching oligonucleotide (TSO)ExiqonN/A(5’-AAGCAGTGGTATCAACGCAGAGTACATrGrG+G-3’)
Sequence-based reagentISPCR primersIDTN/A(5’-AAGCAGTGGTATCAACGCAGAGT-3’)
Sequence-based reagentdNTP Set (100 mM)ThermoFisherCat#: 10297-018A
Sequence-based reagentERCC (External RNA Controls Consortium) ExFold
RNA Spike-In Mixes
ThermoFisherCat#: 4456740
Peptide, recombinant proteinSMARTScribe reverse transcriptaseClontechCat#: 639538
Peptide, recombinant proteinStreptavidin PE-Cyanine7 ConjugateThermoFisherCat#: 25-4317-82
Peptide, recombinant proteinEpidermal growth factorPeproTechCat#: 315–09
Peptide, recombinant proteinPlatelet-derived growth factor BBPeproTechCat#: 315–18
Peptide, recombinant proteinBasic fibroblast growth factorSigma-AldrichCat#: F0291
Peptide, recombinant proteinTransforming growth factor β1PeproTechCat#: 100–21
Commercial assay or kitKAPA HiFi HotStart ReadyMixKapa BiosystemsCat#: KK2602
Commercial assay or kitHS NGS Fragment Kit (1–6000 bp), 1000AgilentCat#: DNF-474–1000
Commercial assay or kitNextera XT DNA Library Preparation kitIlluminaCat#: FC-131–1096
Commercial assay or kitRNeasy Micro KitQiagenCat#: 74004
Commercial assay or kitRNAScope Multiplex Fluorescent V2 AssayAcdbioCat#: 323100
Chemical compoundEDTAThermoFisherCat#: 15573–038
Chemical compoundFluoromount-GThermoFisherCat#: 00-4958-02
Chemical compoundParaformaldehyde (PFA)Electron Microscopy
Sciences
Cat#: 15710
Chemical compoundMatrigelCorningCat#: CB40234A
Chemical compoundBovine serumalbumine (BSA)Sigma-AldrichCat#: A9647
Chemical compoundTriton X-100 (10%)ThermoFisherCat#: 85111
Chemical compoundRecombinant RNase inhibitorClontechCat#: 2313B
Chemical compoundSaffronSigma-AldrichCat#: S8381-5G
Chemical compoundAcid FuchsinSigma-AldrichCat#: F8129-25G
Chemical compoundAcid Red 73Sigma-AldrichCat#: 49823–25 MG
Chemical compoundPhosphotungstic AcidSigma-AldrichCat#: P4006
Chemical compoundHematoxylinSigma-AldrichCat#: MHS32-1L
Chemical compoundShandon Eosin YThermoFisherCat#: 6766009
Chemical compoundType II collagenaseSigma-AldrichCat#: C6885
Chemical compound100 U/ml DNase IWorthingtonCat#: NC9199796
Chemical compoundFetal bovine serum (FBS)ThermoFisherCat#: 16000–069
Chemical compoundOptimal Cutting Temperature compound (OCT)ThermoFisherCat#: 23-730-571
Chemical compoundPenicillin-Streptomycin SolutionThermoFisherCat#: 15140–122
Chemical compoundMedia 199 (M199)Sigma-AldrichCat#: C6885
Chemical compoundAgencourt AMPure XP beadsBeckman CoulterCat#: A63882
Chemical compoundUltraPure DNase/RNase-Free Distilled WaterThermoFisherCat#: 10977023
Chemical compoundTRIzol LSThermoFisherCat#: 10296028
Chemical compoundOil Red OSigma-AldrichCat#: O0625
Chemical compoundAlizarin Red SCarl RothCat#: A5533-25G
Chemical compoundAlcian Blue 8GXSigma-AldrichCat#: A3157
Chemical compoundCrystal VioletSigma-AldrichCat#: C0775
Chemical compoundMCDB201 MediaSigma-AldrichCat#: M6770
Chemical compoundDexamethasoneSigma-AldrichCat#: D-4902
Chemical compoundL-Ascorbic acid 2-phosphateSigma-AldrichCat#: A8960
Chemical compoundInsulin-transferrin-selenium (ITS) mixSigma-AldrichCat#: I3146
Chemical compoundLinoleic acid-AlbuminSigma-AldrichCat#: L9530
Chemical compoundIndomethacinSigma-AldrichCat#: I7378
Peptide, recombinant proteinRecombinant Human InsulinRocheCat#: 11376497001
Chemical compoundIsobutylmethylxanthineSigma-AldrichCat#: I5879
Chemical compound3,3′,5-triiodo-L-thyronine (T3)Sigma-AldrichCat#: T6397
Chemical compoundβ-glycerophosphateSigma-AldrichCat#: G9891
Chemical compoundL-thyroxineSigma-AldrichCat#: T0397
Chemical compoundBetaineSigma-AldrichCat#: B0300
Chemical compoundDTT (DL-dithiothreitol) 100 mMPromegaCat#: P1171
SoftwareImageJNIHhttp://wsr.imagej.net/distros/osx/ij152-osx-java8.zipRRID:SCR_003070
SoftwareFlowJoFLOWJ LLChttps://www.flowjo.com/RRID:SCR_008520
SoftwareBD FACSAria IIBD Bioscienceshttp://www.bdbiosciences.com/cn/home
SoftwareGene Expression Commons (GEXC) databaseSeita et al., 2012https://gexc.riken.jp
Softwarebcl2fastq2 2.18Illuminahttps://support.illumina.com/sequencing/sequencing_software/bcl2fastq-conversion-software.htmlRRID:SCR_015058
SoftwareSkewerJiang et al., 2014https://sourceforge.net/projects/skewerRRID:SCR_001151
SoftwareSTAR 2.4Dobin et al., 2013https://github.com/alexdobin/STAR
SoftwareRSEM 1.2.21Li and Dewey, 2011https://deweylab.github.io/RSEM/RRID:SCR_013027
SoftwareScanpy 1.8.Wolf et al., 2018https://github.com/theislab/scanpyRRID:SCR_018139
SoftwareGraphPad Prism 9.02GraphPad Softwarehttp://www.graphpad.com/scientificsoftware/prismRRID:SCR_002798
OtherRNAscope Probe-Mm-Cdh13-C3AcdbioCat#: 443251RNAscope 1:1500
OtherRNAscope Probe-Mm-Wif1-C3AcdbioCat#: 412361-C3RNAscope 1:1500
Other4′ , 6-diamidino-2-phenylindole (DAPI)BiolegendCat#: 4228011 ug/ml
OtherPropidium iodide (PI)BiolegendCat#: 4213011 ug/ml

Additional files

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Download citations (links to download the citations from this article in formats compatible with various reference manager tools)

Open citations (links to open the citations from this article in various online reference manager services)