RUNX1 marks a luminal castration-resistant lineage established at the onset of prostate development
Figures

RUNX1 marks a subpopulation of mouse proximal prostate luminal cells (PLCs).
(A) Co-immunostaining of RUNX1, NKX3.1, CDH1 in the mouse Anterior Prostate (AP). Higher magnification images of (i) proximal AP and (ii) distal AP are shown. Arrows indicate RUNX1+ NKX3.1- cells, arrowheads show RUNX1- NKX3.1+ cells. Scale bars: 500 µm (yellow) and 50 µm (white). (B) Quantification of RUNX1 and NKX3.1 nuclear intensity (log10) in CDH1+ epithelial cells by QBIC in proximal and distal AP. n = 6–8 mice. (C, D) Flow-cytometry analysis of P2-Runx1:RFP mice, and corresponding quantification of the percentages of RFP+ and RFP cells in the basal and luminal fractions of the proximal and distal AP. n = 7 mice. (E) Experimental strategy to grow organoids from sorted RFP+ and RFP cells from the basal (CD49fhigh) and luminal (CD24high) lineages of P2-Runx1:RFP mouse reporters. (F) Organoid Forming Capacity (OFC) of RFP+ and RFP- basal and luminal sorted cells after 7 days in culture. n = 4 mice. (G) Whole-mount immunostaining of unipotent K5+, unipotent K8+ or multipotent K5+ K8+ organoids. Scale bar: 50 µm. (H) Quantification of the type of organoids characterized by whole-mount immunostaining, as in G. Numbers of organoids quantified are shown above the graph. p Value is indicated for the proportion of K8+ organoids between Proximal AP Luminal RFP+ versus RFP- derived subset. Other comparisons were not statistically significant. n = 2 mice per group. Source files are available in Figure 1—source data 1.
-
Figure 1—source data 1
Source data files for Figure 1.
- https://cdn.elifesciences.org/articles/60225/elife-60225-fig1-data1-v2.xlsx

RUNX1 is enriched in the mouse prostate epithelium.
(A) Co-immunostaining of RUNX1, NKX3.1, CDH1 in the mouse DLP (top) and VP (bottom). Closed arrows indicate RUNX1+ NKX3.1- cells, chevron arrows indicate rare RUNX1+ NKX3.1+ cells, arrowheads show RUNX1- NKX3.1+ cells. Scale bars: 50 µm. (B) Quantification of RUNX1 and NKX3.1 nuclear intensity (log10) in CDH1+ epithelial cells by QBIC in proximal DLP (n = 3), distal DLP (n = 8), proximal VP (n = 3) and distal VP (n = 8). (C) H and E staining of the mouse AP. Higher magnification images of (i) proximal AP and (ii) distal regions are shown. Arrows show typical proximal luminal cells with condensed nuclei and cytoplasm, arrowheads show typical distal secretory luminal cells with enlarged cytoplasm and strong pink eosinophilic staining. Scale bars: 200 µm (blue) and 50 µm (black). Source files are available in Figure 1—figure supplement 1—source data 1.
-
Figure 1—figure supplement 1—source data 1
Source data files for Figure 1—figure supplement 1.
- https://cdn.elifesciences.org/articles/60225/elife-60225-fig1-figsupp1-data1-v2.xlsx

Runx1 expression is mediated by the P2 promoter in the mouse prostate epithelium.
(A) Schematic diagrams of the P2-Runx1:RFP reporter allele (top), Runx1 WT allele (middle), and P1-Runx1:GFP reporter allele (bottom). In P2-Runx1:RFP or P1-Runx1:GFP mice, expression of RFP or GFP is directed by Runx1 proximal P2 or distal P1 promoter respectively. WT Runx1 expression is directed by the remaining WT allele. (B) Flow-cytometry analysis of the epithelial EPCAM+ fraction of WT, P2-Runx1:RFP and P1-Runx1:GFP mice. Percentages are indicated in each quadrant gates. (C) Confocal images indicating P2-Runx1:RFP signal in the glandular epithelium (top), while P1-Runx1:GFP (bottom) activity is detected in the surrounding stroma. Scale bar: 50 µm. (D) Flow-cytometry analysis of P2-Runx1:RFP mice, and corresponding quantification of the percentages of RFP+ and RFP- cells in the basal and luminal fractions of DLP (n = 4 mice) and VP (n = 4 mice) lobes. Source files are available in Figure 1—figure supplement 2—source data 1.
-
Figure 1—figure supplement 2—source data 1
Source data files for Figure 1—figure supplement 2.
- https://cdn.elifesciences.org/articles/60225/elife-60225-fig1-figsupp2-data1-v2.xlsx

Characterization of P2-Runx1:RFP derived mouse prostate organoids.
(A) Organoid Forming Capacity (OFC) of RFP+ and RFP- basal and luminal sorted cells from DLP (left) and VP (right) cultured for 7 days. n = 4 mice. (B) Quantification of organoid morphologies determined by brightfield assessment. p Values are indicated only for statistically significant (p<0.05) comparisons between RFP+ versus RFP- derived subsets. n = 4 mice. (C) Quantification of the type of organoids (pink: unipotent K5+, green: unipotent K8+, orange: multipotent K5+ K8+ organoids) characterized by whole-mount immunostaining. p Value is indicated for the proportion of K5+K8+ organoids between VP Luminal RFP+ versus RFP- derived subset. Other comparisons were not statistically significant. n = 2 mice. (D) Quantification of organoid morphologies according to the type of organoids (unipotent K5+, unipotent K8+, multipotent K5+ K8+). Results are shown for luminal proximal and distal AP derived organoids. n = 2 mice. Numbers of organoids quantified are shown above the graph. Two-way ANOVA using Tukey’s multiple comparisons test was used for statistical analysis. (E) Quantification of the organoid area estimated from maximum projection of stacked images from basal (left) and luminal (right) fractions of the proximal AP. Numbers of organoids quantified are shown above the graph. Two-tailed Mann-Whitney test was used for statistical analysis. n = 2 mice. Source files are available in Figure 1—figure supplement 3—source data 1.
-
Figure 1—figure supplement 3—source data 1
Source data files for Figure 1—figure supplement 3.
- https://cdn.elifesciences.org/articles/60225/elife-60225-fig1-figsupp3-data1-v2.xlsx

RUNX1-expressing cells are enriched in the castrated prostate epithelium.
(A) P2-Runx1:RFP reporter mice were surgically castrated between 6 and 12 weeks of age and analyzed at least 4 weeks post-castration. (B, C) Flow-cytometry analysis and corresponding quantification of the proportion of RFP+ and RFP- cells in the EPCAM+ fraction of intact and castrated prostates of P2-Runx1:RFP mice. n = 3 mice per group. Int: Intact, Cas: Castrated. (D) Co-immunostaining of RUNX1, NKX3.1, CDH1 in the castrated wild-type mouse prostate. Higher magnification images of (i) proximal, (ii) intermediate, and (iii) distal AP are shown. Arrows indicate RUNX1- NKX3.1+ cells, arrowheads show a luminal cell co-stained for RUNX1 and NKX3.1. Amp: ampullary gland. Scale bars: 500 µm (yellow) and 50 µm (white). Int: Intact, Cas: Castrated. (E) Co-immunostaining of RUNX1 and TROP2 showing colocalization of the two markers in both proximal (bottom) and distal (top) castrated AP. Scale bars: 50 µm (white). Source files are available in Figure 2—source data 1.
-
Figure 2—source data 1
Source data files for Figure 2.
- https://cdn.elifesciences.org/articles/60225/elife-60225-fig2-data1-v2.xlsx

Characterization of RUNX1 expression in the castrated mouse prostate.
(A) Flow cytometry analysis and corresponding quantifications of the basal and luminal EPCAM+ fraction from AP, DLP, and VP lobes of castrated P2-Runx1:RFP mice (n = 3). (B, C) Co-immunostaining of RUNX1, NKX3.1, CDH1 in castrated DLP (B) and VP (C) lobes. Scale bars: 50 µm (white). Source files are available in Figure 2—figure supplement 1—source data 1.
-
Figure 2—figure supplement 1—source data 1
Source data files for Figure 2—figure supplement 1.
- https://cdn.elifesciences.org/articles/60225/elife-60225-fig2-figsupp1-data1-v2.xlsx

scRNA-seq profiling of intact and castrated Runx1+ cells reveals transcriptomic similarity between proximal luminal cells and castration-resistant cells.
(A) Experimental strategy for scRNA-seq on RFP+ and RFP- cells individually dissected lobes of intact and castrated prostates isolated from P2-Runx1:RFP reporter mice. (B) UMAP visualization (left) and graph-abstracted representation (PAGA, right) of prostate epithelial cells (n = 3,825 cells from three independent experiments). Colors represent different clusters. In PAGA, clusters are linked by weighted edges that represent a statistical measure of connectivity. (C) Dot plot showing the expression of selected marker genes associated with each cluster. (D–H) UMAP visualization of prostate epithelial cells. Cells in D and G are colored by a gradient of log-normalized expression levels for each gene indicated. Cell colors in E represent the treatment of origin (intact, castrated), in F individual lobes of origin (AP, DLP, VP), and in H RFP FACS gate of origin (RFP+, RFP-).

Pre-processing of the scRNA-seq dataset of adult intact and castrated mouse prostates.
(A) UMAP visualization in MULTI-seq barcodes space for each independent experiment (Run 1: left, Run 2: middle, Run 3: right) before (top) and after (bottom) quality control and barcode filtration (see Materials and methods for details). (B) UMAP visualization of the integrated batch-corrected dataset (n = 4499 cells from three independent experiments). Colors represent different clusters. (C) UMAP visualization of the integrated batch-corrected dataset split and colored by individual experiment (Run 1: left, Run 2: middle, Run 3: right). (D, E) Example key marker genes used for cell type characterization. (D) UMAP visualization of the integrated batch-corrected dataset. Cells are colored by a gradient of log-normalized expression levels for each gene indicated. (E) Dot plot showing the expression of selected marker genes associated with each cluster. Cells corresponding to the prostatic clusters were kept for downstream analyses.

Characterization of the scRNA-seq prostate epithelial subset.
(A–D) Number of cells present in the prostate epithelial subset. In A colors correspond to phenotypic sorted populations, in B treatment of origin, in C RFP FACS gate of origin, in D individual lobes. (E, G) UMAP visualization of prostate epithelial cells. Cells in E are colored by clusters. Clusters 0, 1, and 5 were merged together. In F and G colors represent phenotypic sorted populations. (H–J) Cells are colored by a gradient of log-normalized expression levels for each gene. Canonical epithelial markers are shown in H, luminal lineage markers in I, and basal lineage markers in J.

Characterization of the scRNA-seq prostate epithelial dataset.
(A) Heatmap of the 10 most differentially upregulated genes per cluster. (B, C) Sankey diagrams summarizing the proportion of cells and between phenotypic groups. (D, E) Cells are colored by a gradient of log-normalized expression levels for each gene. D shows the expression of Ly6 family members, and in E RUNX family members are shown.

Gene Ontology and differential expression analysis within the scRNA-seq prostate epithelial dataset.
Bar plots of the eight most significantly (g:Profiler adjusted p-value<0.05) enriched gene ontology (GO) terms (GO:BP, Biological Processes) using the list of differentially upregulated genes specific to (A) Lum-A, (B) Lum-B, (C) Lum-C, (D) Lum-D, (E) Lum-E, (F) Lum-F, and (G) Basal clusters. (H, I) Fifteen most enriched GO terms (g:Profiler adjusted p-value<0.05) for genes upregulated in the Lum-D cluster against combined Lum-A, Lum-B, and Lum-C clusters in H. Results of the opposite comparison are shown in I. (J) Violin plots showing the expression levels of specific genes within the Lum-D cluster. Data are split and colored by treatment for each lobe. The expression of Krt4, Tacstd2/Trop2 and Runx1 does not vary upon treatment, while Androgen Receptor (AR) regulated genes such as Psca and Tspan1 are downregulated after castration.

Lineage tracing of Runx1-expressing cells establishes the intrinsic castration-resistant properties of the proximal luminal lineage.
(A) Schematic summary of the genetic lineage-tracing system employed. (B) Experimental strategy for lineage-tracing experiments. (C) Co-immunostaining of RFP, RUNX1, CDH1 in the proximal AP. Arrows indicate RFP labeled RUNX1+ cells. Scale bar: 50 µm. (D) Quantification of the percentage of epithelial RFP+ cells in proximal and distal regions of the prostate in intact (n = 5), castrated (n = 4) and regenerated (n = 4) mice. (E) Quantification of the percentage of epithelial RFP+ cells in intact (n = 5), castrated (n = 5) and regenerated (n = 4) mice. (F) Co-immunostaining of RFP, K5, K8 in the proximal AP, distal AP, and DLVP (DLP + VP). Arrowheads indicate RFP labeled basal cells (K5+) found in distal AP, the white arrow indicates a luminal (K8+) RFP+ clone made of two cells. Scale bar: 50 µm. (G) Quantification of the percentage of epithelial RFP+ clones comprising between two and four cells in intact (n = 5), castrated (n = 5) and regenerated (n = 4) mice. (H, I) Quantification of the percentage of RFP+ cells being K5+ or K8+ in H, or NKX3.1+ or NKX3.1- in I, in intact (n = 5), castrated (n = 5) and regenerated (n = 4) mice. Int: Intact, Cas: Castrated, Rgn: Regenerated. Source files are available in Figure 4—source data 1.
-
Figure 4—source data 1
Source data files for Figure 4.
- https://cdn.elifesciences.org/articles/60225/elife-60225-fig4-data1-v2.xlsx

Lineage tracing of RUNX1-expressing cells labeled in intact mice.
(A) Scheme of the tamoxifen labeling strategy to evaluate the labeling efficiency of Runx1-expressing cells. (B) Flow-cytometry analysis of intact (n = 2), castrated (n = 2), and regenerated (n = 2) Runx1CreER Rosa26LSL-RFP mice, and corresponding quantification of the percentage of RFP+ cells in the epithelial EPCAM+ fraction. (C) Estimated absolute number of epithelial RFP+ cells in intact (n = 5) and castrated (n = 5) and regenerated (n = 4) prostates based on IHC quantifications. Int: Intact, Cas: Castrated, Rgn: Regenerated. (D) Quantification of the percentage of epithelial RFP+ cells being K4+ or K4 in intact (n = 3), castrated (n = 4) and regenerated (n = 5) mice. (F) Co-immunostaining showing the enrichment of K4+ RFP labeled cells in the proximal region of the AP (top) and DLVP (bottom) prostate lobes. Scale bar: 50 µm. Source files are available in Figure 4—figure supplement 1—source data 1.
-
Figure 4—figure supplement 1—source data 1
Source data files for Figure 4—figure supplement 1.
- https://cdn.elifesciences.org/articles/60225/elife-60225-fig4-figsupp1-data1-v2.xlsx

RUNX1 marks proximal cells during embryonic prostate development.
(A–C) Co-immunostainings of the mouse urogenital sinus at E18.5 for RUNX1, K5, K8 in A, RUNX1, p63, CDH1 in B, RUNX1, NKX3.1, CDH1 in C. Scale bar: 50 µm. (D) Co-immunostainings of RUNX1, NKX3.1, CDH1 at postnatal (P) day 14. Higher magnification images of (p) proximal, (i) intermediate, and (d) distal regions are shown. Scale bars: 200 µm (yellow) and 50 µm (white). Amp: ampullary gland; Sv: seminal vesicles; Ur: urethra; p: proximal; i: intermediate; d: distal. (E) Scheme of the protocol to culture ex vivo explants of mouse UGS harvested at E15.5. (F) Representative images of UGS explants at E15.5 (day 0), day 3, day 5, and day 7 of culture showing the formation of premature prostate buds. (G) Co-immunostaining of RUNX1, NKX3.1, CDH1 in UGS explants harvested at day 0, day 1, day 2, and day 7. Higher magnification images of each square (left) are shown for each time point. Chevron arrows show RUNX1+ NKX3.1+ cells, closed arrows indicate RUNX1- NKX3.1+ cells, arrowheads show RUNX1+ NKX3.1- cells. Scale bars: 200 µm (yellow) and 50 µm (white). (H) Quantification of RUNX1 and NKX3.1 nuclear intensity (log10) in CDH1+ epithelial cells of UGS explants by QBIC. Quantification was performed within the boundaries delimited in G by dotted lines, at day 0 (n = 3 explants), day 1 (n = 7 explants), day 2 (n = 6 explants), and day 7 (n = 6 explants). Source files are available in Figure 5—source data 1.
-
Figure 5—source data 1
Source data files for Figure 5.
- https://cdn.elifesciences.org/articles/60225/elife-60225-fig5-data1-v2.xlsx

Characterization of RUNX1 expression during prostate development in vivo.
(A) Co-immunostainings of the mouse urogenital sinus at E18.5 for RUNX1, NKX3.1, K4. p: proximal, d: distal. Scale bar: 50 µm. Dotted lines indicate the urogenital epithelium. (B) Co-immunostainings of the mouse urogenital sinus at E18.5 for RUNX1, NKX3.1, LY6D. p: proximal, d: distal. Scale bar: 50 µm. (C) Co-immunostaining of RUNX1, NKX3.1, K4 at postnatal (P) day 14. Higher magnification images of (p) proximal, (i) intermediate, and (d) distal regions are shown. Scale bars: 200 µm (yellow) and 50 µm (white). (D, E) Co-immunostainings at P14 for RUNX1, K5, K8 in B RUNX1, p63, CDH1 in C. p: proximal, i: intermediate, d: distal. Scale bar: 50 µm.

Characterization of RUNX1 expression during prostate development in UGS explant cultures.
(A) Co-immunostaining of RUNX1, K5, K8 in UGS explants harvested at day 2 and day 7. p: proximal, d: distal. Scale bar: 50 µm. (E) Quantification of the percentage of K5+, K8+, K5+ K8+ cells during UGS explant cultures at day 0 (n = 2 explants), day 1 (n = 3 explants), day 2 (n = 3 explants), day 7 (n = 4 explants), showing a progressive reduction in the proportion of K5+ K8+ double positive cells. (C) Co-immunostaining of RUNX1, NKX3.1, LY6D in UGS explants harvested at day 7. p: proximal, d: distal. Scale bar: 50 µm. (D) Co-immunostaining of RUNX1, NKX3.1, K4 in UGS explants harvested at day 1 (top), 3 (middle), and 7 (bottom). Higher magnification images of (p) proximal, (i) intermediate, and (d) distal regions are shown for day 7. Scale bars: 200 µm (yellow) and 50 µm (white). Dotted lines indicate the urogenital epithelium. (E) Co-immunostaining of RUNX1, Ki67, CDH1 in UGS explants harvested at day 2 and day 7. p: proximal, d: distal. Scale bar: 50 µm. Source files are available in Figure 5—figure supplement 2—source data 1.
-
Figure 5—figure supplement 2—source data 1
Source data files for Figure 5—figure supplement 2.
- https://cdn.elifesciences.org/articles/60225/elife-60225-fig5-figsupp2-data1-v2.xlsx

scRNA-seq of UGS explant cultures shows specification of the proximal luminal lineage during embryonic prostate development.
(A) Experimental strategy for scRNA-seq of UGS explant cultures at day 0, day 1, day 3 and day 6. (B, C) Force directed visualization of the developing prostatic epithelium in UGS explant cultures. In B cells are colored by experimental time points, and in C cells are colored by clusters. (D) PAGA representation of the clusters as in C. Weighted edges between cluster nodes represent a statistical measure of connectivity. (E) Fraction of cells per cluster at each experimental time point, displaying a progressive cellular diversification. (F–H) PAGA representations with cluster nodes colored by a gradient representing the mean log-normalized expression levels of each gene. (I) Force directed visualization of the developing prostatic epithelium in UGS explant cultures. Color gradient represents AUC scores per cell. Per-cell AUC scores were calculated using the ‘AUCell’ package. Gene signatures for ‘Lum-D’ (left) and ‘Basal’ (right) were generated using the list of differentially upregulated genes previously obtained from our adult mouse prostate clusters.

Pre-processing of the scRNA-seq dataset of UGS explant cultures.
(A) UMAP visualization in MULTI-seq barcodes space before (top) and after (bottom) quality control and barcode filtration (see Materials and methods for details). (B) UMAP visualization of the dataset colored by estimated cell cycle phase (top). Due to the strong underlying impact on clustering, the cell cycle effect was regressed out (bottom). (C, D) UMAP visualization of the dataset colored in C by experimental time point and in D by labeled populations. (E) Dot plot showing the expression of selected marker genes associated with each labeled population. (F) Bar plots of the 10 most significantly (g:Profiler adjusted p-value<0.05) enriched gene ontology (GO) terms (GO:BP, Biological Processes) on the differentially upregulated genes specific to the ‘Hypoxic/Stressed’ cluster. (G–I) UMAP visualization of the prostatic subset (from D), labeled by time points in G, clusters in H, estimated cell cycle phase in I. (J) Percentage of cells in the different estimated cell cycle phases in each clusters. (K) Diffusion maps representation of the prostatic subset using the first two diffusion components with (top) and without C9 (bottom). C9 was excluded for downstream analyses due to its outlying profile.

Characterization of the developing prostatic epithelium in the scRNA-seq dataset of UGS explant culture.
(A) Heatmap of the 10 most differentially upregulated genes per cluster. (B–E) PAGA representations with cluster nodes colored by a gradient representing the mean log-normalized expression levels of each genes indicated. B shows the expression of common prostate lineage markers, C shows known regulators of prostate development, D shows basal markers, E shows proliferation markers. (F) Box plots of per-cell AUC signature scores calculated using the ‘AUCell’ package for individual UGS explant clusters. Gene signatures were generated using the list of differentially upregulated genes previously obtained from our adult mouse prostate clusters. (G) Force-directed visualization of the developing prostatic epithelium in UGS explant cultures. Color gradient represents AUC scores per cell.

RUNX1+ cells contribute to the establishment of the proximal luminal lineage during embryonic prostate development.
(A) Strategy for lineage tracing of RUNX1+ cells in UGS explant cultures. Tamoxifen was applied on day 0 and day 1 and washed out on day 2. (B, C) Co-immunostaining of RFP, RUNX1, CDH1 in UGS explants harvested at day 2 (B) and day 7 (C). Higher magnification images of proximal (i) and (ii) distal regions are shown for day 7. Arrows show RFP+ RUNX1-low cells, arrowheads show RFP+ RUNX1+ cells. Scale bars: 200 µm (yellow) and 50 µm (white). (C) Quantification of the percentage of epithelial RUNX1+ cells in the RFP subset at day 2 (n = 7) and day 7 (n = 3) of UGS explant cultures. Quantification was performed within the boundaries delimited in B by dotted lines. (E, F) Co-immunostaining of RFP, NKX3.1, CDH1 in UGS explants harvested at day 2 (E) and day 7 (F). Higher magnification images of (i) proximal and (ii) distal regions are shown for day 7. Arrows show RFP+ NKX3.1+ cells, arrowheads show RFP+ NKX3.1- cells. Scale bars: 200 µm (yellow) and 50 µm (white). (G) Quantification of the percentage of epithelial NKX3.1+ cells in the RFP subset at day 2 (n = 6) and day 7 (n = 4) of UGS explant cultures. Quantification was performed within the boundaries delimited in F by dotted lines. (H, I) Co-immunostaining of RFP, K4, CDH1 in UGS explants harvested at day 2 (H) and day 7 (I). Higher magnification images of (i) proximal and (ii) distal regions are shown for day 7. Arrows show RFP+ K4 cells, arrowheads show RFP+ K4+ cells. Scale bars: 200 µm (yellow) and 50 µm (white). (J) Quantification of the percentage of epithelial K4+ cells in the RFP subset at day 2 (n = 3) and day 7 (n = 3) of UGS explant cultures. Quantification was performed within the boundaries delimited in I by dotted lines. Source files are available in Figure 7—source data 1.
-
Figure 7—source data 1
Source data files for Figure 7.
- https://cdn.elifesciences.org/articles/60225/elife-60225-fig7-data1-v2.xlsx

Lineage tracing of RUNX1-expressing cells in UGS explants.
(A, B) Co-immunostaining of RFP, p63, CDH1 in UGS explants harvested at day 2 (A) and day 7 (B). Higher magnification images of (i) proximal and (ii) distal regions are shown for day 7. Arrows show RFP+ p63+ cells, arrowheads show RFP+ p63- cells. Scale bars: 200 µm (yellow) and 50 µm (white). (C) Quantification of the percentage of epithelial RFP+ p63+ cells at day 2 (n = 7) and day 7 (n = 4) of UGS explant cultures. Quantification was performed within the boundaries delimited in B by dotted lines. (D) Co-immunostaining of RFP, K5, K8 in UGS explants harvested at day 2 (top) and day 7 (bottom). Higher magnification images of proximal (i) and (ii) distal regions are shown for day 7. Arrows show RFP+ K5+ K8+ cells, chevron arrows show RFP+ K5- K8+, arrowheads show RFP+ K5+ K8- cells. Scale bars: 200 µm (yellow) and 50 µm (white). (E) Quantification of the percentage of epithelial K5+ K8+, K5- K8+ cells and K5+ K8- cells in the RFP subset at day 2 (n = 4) and day 7 (n = 3) of UGS explant cultures. Quantifications were performed within the boundaries delimited in D by dotted lines. (F, G) Co-immunostaining of RFP, Nupr1 (mRNA), CDH1 in UGS explants harvested at day 2 (F) and day 7 (G). Higher magnification images of (i) proximal and (ii) distal regions are shown for day 7. Arrows show RFP+ Nupr1+ cells, arrowheads show RFP+ Nupr1- cells. Scale bars: 200 µm (yellow) and 50 µm (white). (J) Quantification of the percentage of epithelial K4+ cells in the RFP subset at day 2 (n = 3) and day 7 (n = 4) of UGS explant cultures. Quantification was performed within the boundaries delimited in G by dotted lines. Source files are available in Figure 7—figure supplement 1—source data 1.
-
Figure 7—figure supplement 1—source data 1
Source data files for Figure 7—figure supplement 1.
- https://cdn.elifesciences.org/articles/60225/elife-60225-fig7-figsupp1-data1-v2.xlsx
Tables
Reagent type (species) or resource | Designation | Source or reference | Identifiers | Additional information |
---|---|---|---|---|
Strain, strain background (Mus musculus, male) | ICR (CD-1) wild-type | Envigo | Hsd:ICR (CD-1) | 7–15 week old males |
Strain, strain background (Mus musculus, male) | P1-Runx1:GFP | Georges Lacaud lab | 7–15 week old males | |
Strain, strain background (Mus musculus, male) | P2-Runx1:RFP | Georges Lacaud lab | 7–15 week old males | |
Strain, strain background (Mus musculus, male) | Runx1mER-CRE-mER Runx1CreER Rosa26LSL-RFP | RIKEN (Japan) Samokhvalov et al., 2007 | Runx1-MER-Cre-MER | C57Bl/6J background 7–15 week old males |
Strain, strain background (Mus musculus, male) | Rosa26Lox-Stop-Lox-tdRFP Runx1CreER Rosa26LSL-RFP | European Mouse Mutant Archive Luche et al., 2007 | B6.Cg-Thy1 Gt(ROSA)26Sortm1Hjf | C57Bl/6J background 7–15 week old males |
Antibody | Anti-RUNX1 (rabbit monoclonal) | Cell Signaling | Cat: 8529 RRID:AB_10950225 | IHC/IF (1:100) |
Antibody | Anti-NKX3.1 (rabbit polyclonal) | Athenaes | Cat: AES-0314 | IHC/IF (1:200) |
Antibody | Anti-CDH1 (goat polyclonal) | R and D Systems | Cat: AF748 AB_355568 | IHC/IF (1:400) |
Antibody | Anti-p63 (rabbit monoclonal) | Cell Signaling | Cat: 39692 RRID:AB_2799159 | IHC/IF (1:800) |
Antibody | Anti-K5 (rabbit monoclonal) | Abcam | Cat: ab52635 RRID:AB_869890 | IHC/IF (1:400) |
Antibody | Anti-K8 (rabbit monoclonal) | Abcam | Cat: ab53280 RRID:AB_869901 | IHC/IF (1:400) |
Antibody | Anti-K4 (mouse monoclonal) | Abcam | Cat: Ab9004 RRID:AB_306932 | IHC/IF (1:100) |
Antibody | Anti-LY6D (rabbit polyclonal) | Proteintech | Cat: 17361–1-AP | IHC/IF (1:100) |
Antibody | Anti-TROP-2 (goat polyclonal) | R and D Systems | Cat: AF1122 RRID:AB_2205662 | IHC/IF (1:200) |
Antibody | Anti-BrdU (rat monoclonal) | Abcam | Cat: ab6326 RRID:AB_305426 | IHC/IF (1:400) |
Antibody | Anti-Ki67 (rabbit monoclonal) | Abcam | Cat: ab15580 RRID:AB_443209 | IHC/IF (1:800) |
Antibody | Anti-RFP (rabbit polyclonal) | Rockland | Cat: 600-402-379 RRID:AB_828391 | IHC/IF (1:400) |
Antibody | Anti-RFP (rabbit monoclonal) | MBL | Cat: PM005 RRID:AB_591279 | IF (1:200) |
Antibody | Anti-GFP (rabbit polyclonal) | MBL | Cat: 598 RRID:AB_591816 | IF (1:200) |
Antibody | EnVision+/HRP Anti-Rabbit | Dako (Agilent) | Cat: K4003 RRID:AB_2630375 | IHC/IF Ready to use |
Antibody | EnVision+/HRP Anti-Rabbit | Dako (Agilent) | Cat: K4001 RRID:AB_2827819 | IHC/IF Ready to use |
Antibody | ImmPRESS HRP Anti-Goat | Vector Laboratories | Cat: MP-7405 RRID:AB_2336526 | IHC/IF Ready to use |
Antibody | ImmPRESS HRP Anti-Rat | Vector Laboratories | Cat: MP-7444 RRID:AB_2336530 | IHC/IF Ready to use |
Antibody | Donkey anti-Goat IgG 647 | ThermoFischer Scientific | Cat: A-21447 RRID:AB_141844 | IF (1:400) |
Antibody | Anti-CD16/32 Fc block | Biolegend | Cat: 101301 Clone: 93 RRID:AB_312800 | FACS (1:200) |
Antibody | Anti-CD45 SB436 | ThermoFischer Scientific | Cat: 62-0451-82 Clone: 30-F11 RRID:AB_2744774 | FACS (1:200) |
Antibody | Anti-EPCAM BV421 | Biolegend | Cat: 118225 Clone: G8.8 RRID:AB_2563983 | FACS (1:200) |
Antibody | Anti-EPCAM APC | Biolegend | Cat: 118214 Clone: G8.8 RRID:AB_1134102 | FACS (1:200) |
Antibody | Anti-CD49f FITC | Biolegend | Cat: 313606 Clone: GoH3 RRID:AB_345300 | FACS (1:200) |
Antibody | Anti-CD49f APC | Biolegend | Cat: 313616 Clone: GoH3 RRID:AB_1575047 | FACS (1:200) |
Antibody | Anti-CD24 BV786 | BD Biosciences | Cat: 744470 Clone: M1/69 RRID:AB_2742258 | FACS (1:200) |
Sequence-based reagent | MULTI-seq reagents | Zev Gartner lab McGinnis et al., 2019b | ||
Software, algorithm | R v3.6.3 | CRAN R Project | SCR_001905 | https://cran.r-project.org |
Software, algorithm | deMULTIplex | McGinnis et al., 2019b | https://github.com/chris-mcginnis-ucsf/MULTI-seq | |
Software, algorithm | DoubletFinder | McGinnis et al., 2019a | SCR_018771 | https://github.com/chris-mcginnis-ucsf/DoubletFinder |
Software, algorithm | Seurat v3.1.5 | Satija et al., 2015; Rahul Satija lab | SCR_016341 | https://github.com/satijalab/seurat |
Software, algorithm | Scanpy v1.4.6 PAGA | Wolf et al., 2019 | SCR_018139 | https://scanpy.readthedocs.io/en/stable/ |
Software, algorithm | AUCell v1.8.0 | Aibar et al., 2017 | https://github.com/aertslab/AUCell | |
Software, algorithm | scater v1.14.6 | Bioconductor | SCR_015954 | https://bioconductor.org/packages/release/bioc/html/scater.html |
Software, algorithm | QuPath v0.2 | Bankhead et al., 2017 | SCR_018257 | https://qupath.github.io/ |
Software, algorithm | Cellranger v3.1.0 | 10x Genomics | SCR_017344 | |
Software, algorithm | FlowJo v10 | BD Life Sciences | SCR_008520 | |
Software, algorithm | Harmony | PerkinElmer | SCR_018809 | |
Software, algorithm | Graphpad Prism v8.4.2 | Graphpad | SCR_002798 |
Additional files
-
Supplementary file 1
Quality control metrics and metadata of scRNAseq experiments.
- https://cdn.elifesciences.org/articles/60225/elife-60225-supp1-v2.xlsx
-
Supplementary file 2
Genes differentially expressed between adult clusters.
- https://cdn.elifesciences.org/articles/60225/elife-60225-supp2-v2.xlsx
-
Supplementary file 3
Genes differentially expressed in intact versus castrated Lum-D cells.
- https://cdn.elifesciences.org/articles/60225/elife-60225-supp3-v2.xlsx
-
Supplementary file 4
Genes differentially expressed between UGS explants clusters.
- https://cdn.elifesciences.org/articles/60225/elife-60225-supp4-v2.xlsx
-
Supplementary file 5
Gene lists used for the analysis of gene set activity.
- https://cdn.elifesciences.org/articles/60225/elife-60225-supp5-v2.xlsx
-
Transparent reporting form
- https://cdn.elifesciences.org/articles/60225/elife-60225-transrepform-v2.docx