The formation of the proper number of functional nephrons requires a delicate balance between renal progenitor cell self-renewal and differentiation. The molecular factors that regulate the dramatic expansion of the progenitor cell pool and differentiation of these cells into nephron precursor structures (renal vesicles) are not well understood. Here we show that Sall1, a nuclear transcription factor, is required to maintain the stemness of nephron progenitor cells. Transcriptional profiling of Sall1 mutant cells revealed a striking pattern, marked by the reduction of progenitor genes and amplified expression of renal vesicle differentiation genes. These global changes in gene expression were accompanied by ectopic differentiation at E12.5 and depletion of Six2+Cited1+ cap mesenchyme progenitor cells. These findings highlight a novel role for Sall1 in maintaining the stemness of the progenitor cell pool by restraining their differentiation into renal vesicles.
During kidney development, the formation of the proper number of functional nephrons is crucial for the future health of the adult. The kidney has a limited capacity to repair existing nephrons upon injury or disease, but lacks the ability to generate new nephrons. A reduction in nephron number can cause hypoplastic disease, a leading cause of childhood kidney failure, or promote hypertension and chronic kidney disease later in life (Keller et al., 2003). Because nephrogenesis ceases shortly after birth, therapies to promote increased nephron number during this crucial window could reduce future kidney-related mortality. However, our understanding of factors that regulate nephron number has been limited.
One of the most important mechanisms to ensure a complete renal endowment is the proper maintenance of the renal progenitor cell pool. This pool gives rise to all the tubular segments of the nephron, except collecting ducts, and has been identified using a number of transgenic tools (Boyle et al., 2008; Kobayashi et al., 2008; Mugford et al., 2008). These studies have shown that the progenitor cells are induced within the metanephric mesenchyme (MM) by the outgrowth of the ureteric bud (UB), express Six2 and Cited1, and form an organized cap of cells around the tips of the UB (cap mesenchyme, CM). The CM receives signals from the UB and is triggered to adopt one of two independent fates: to undergo self-renewal to maintain the progenitor pool or to initiate a mesenchymal-to-epithelial transition to differentiate into the segmented nephron. Differentiation triggers sequential morphological and gene expression changes. First, a clustered pre-tubular aggregate (PTA) undergoes epithelialization into a circular renal vesicle (RV). Next, polarization of the RV forms the comma and S-shaped body, followed by segmentation into specialized podocyte (glomerular), proximal and distal structures of the nephron (Little and McMahon, 2012). The formation of new nephrons occurs in the nephrogenic zone, located in the periphery of the kidney, and is synchronized with the iterative branching of the UB to form a new nephron at the tip of each UB.
A delicate balance between the expansion of the progenitor cell pool and epithelialization to form RVs must be maintained throughout nephrogenesis to ensure the complete complement of 300,000 to 1 million human nephrons (Nyengaard and Bendtsen, 1992). Depletion of the progenitor cell pool by unrestrained differentiation can result in early arrest of nephrogenesis, causing renal hypoplasia or agenesis. Important studies have identified several factors responsible for the maintenance of the renal progenitor cell population. Wnt9b, secreted from the UB, acts as the initial inductive signal for the adjacent CM and simultaneously promotes the expansion of the progenitor pool (Carroll et al., 2005; Karner et al., 2011). Interpretation of the Wnt9b signal seems to be dependent on the level of expression of Six2 within the niche. Six2, expressed in the CM, maintains the nephron progenitors in the undifferentiated state. Deletion of Six2 causes the CM to prematurely differentiate into RV structures, exhausting the progenitor cell pool and halting kidney formation (Self et al., 2006). Fgf9 and Fgf20 are also important in the CM to maintain the multipotent state of progenitor cells (Barak et al., 2012). Simultaneous deletion of both factors causes premature differentiation, similar to the phenotype observed in the Six2 mutant. Together, these studies represent key advances in understanding what factors maintain the renal progenitor cell niche. However, apart from Six2, there is limited information about the transcriptional regulators that determine whether a cell in the CM will self-renew or undergo differentiation to form the RV.
Sall1 encodes a multi-zinc-finger transcription factor highly expressed in multipotent renal progenitor cells (Osafune et al., 2006) and CM-derived differentiating structures (PTA, RV, comma and S-shaped bodies) (Takasato et al., 2004). Mutations in human SALL1 cause Townes Brocks syndrome (TBS; Online Mendelian Inheritance in Man, no. 107408), an autosomal dominant disorder associated with a variety of multi-organ defects, including renal hypoplasia and renal agenesis (Kohlhase, 2000; Kohlhase et al., 1998). Sall1 plays an important role for the initial outgrowth of the UB into the surrounding MM early during metanephric development (Nishinakamura et al., 2001). However, in a majority of Sall1 mutants, the UB invades the mesenchyme and undergoes several rounds of branching before arresting. Using a colony-formation assay, Osafune and colleagues showed that isolated Sall1-expressing cells were able to form colonies from a single cell and were also competent to differentiate into all tubular structures of the nephron. The same population of isolated multipotent cells deficient in Sall1 also formed colonies and were competent to differentiate, although the colony size was significantly smaller (Osafune et al., 2006), indicating that Sall1 is important for the expansion of the progenitor population, but is not required for its differentiation.
We sought to identify the molecular mechanism by which Sall1 functions during kidney development. In addition to the early function of Sall1 for UB invasion (Nishinakamura et al., 2001), we show here that Sall1 is required to expand the renal progenitor cell pool by regulating the differentiation of these cells into RVs. In the absence of Sall1, the progenitor cell pool is exhausted as a result of unrestrained differentiation, leading to the early arrest of kidney development.
A majority of Sall1 mutants form small kidneys with reduced branching
In mice of the ICR background, the UB invaded the MM, branched and formed a small kidney in 88% of Sall1GFP homozygous mutant embryos (Fig. 1A,G). We took advantage of the mutant metanephroi that formed small kidneys to investigate the molecular mechanism of Sall1 in kidney development. Branching was significantly reduced in the mutant from embryonic day (E)12 to 16 compared with heterozygous controls (Fig. 1C,D). The expression of the UB tip markers Ret and Wnt11 were not altered in the mutant. In the wild type, Wnt9b was expressed in the stalks of the ureter and downregulated in the tips of the UB. In the mutant, Wnt9b expression was expanded from the stalks to the tips of the ureter (Fig. 1E).
We characterized the distribution of the UB branching events in E11.5 and E12.5 Sall1GFP heterozygous control and mutant kidneys (Fig. 1F,G). We found that a majority of mutant kidneys branched by E12.5, but were delayed. At E11.5, 66% of the heterozygous kidneys were at T-stage and 34% had branched to form six to eight UB tips. By contrast, we did not observe any mutant kidneys with six to eight UB tips at E11.5. Instead, 35% of mutant kidneys had ureters that were unbranched, with the remainder of the kidneys at the T-stage. The percentage of kidneys with unbranched ureters in the mutant dropped to 12% by E12.5. Half of the mutant kidneys at E12.5 were at the T-stage, 32% had six to eight UB tips, and 7% of kidneys had 10 or more UB tips. A majority of the heterozygous kidneys had 10 or more UB tips at E12.5. These data indicate that whereas UB invasion of MM is delayed in the Sall1 mutant, in most cases the UB invades the mutant MM and undergoes several rounds of branching to form a small kidney.
Sall1GFP mutants have ectopic renal vesicle formation
Histological analysis of E11.5-E17 kidneys revealed a number of RV or tubular-like structures (hereafter referred to as RV/tubule structures) throughout the mutant kidney, including some in ectopic locations toward the periphery of the kidney (Fig. 2A). Interestingly, at the T-stage, the mutant had ectopically located condensed clusters of cells with early lumen formation on the dorsal side of the ureter. At E13.5 we observed few UB tips, numerous differentiated RV-like structures, and very few differentiated comma and S-shaped bodies. At E17.5 we observed numerous differentiated structures ectopically located toward the periphery of the kidney, including structures that resembled glomeruli. To investigate this observation further we performed immunohistochemistry for Sall1-GFP (expressed in CM, stroma, RV and differentiating tubule structures) and cytokeratin (expressed in UB; Fig. 2B). We were able to use an N-terminal antibody to Sall1 because the Sall1GFP mouse expresses 52 N-terminal amino acids fused to GFP. Similar to the histological analysis, in the mutant we observed many differentiated RV/tubule structures that were Sall1-GFP positive (Sall1-GFP+) and cytokeratin negative. These structures were also laminin positive and cytokeratin negative at E17.5 (Fig. 2B). There were very few UB tips in the mutant. We observed Sall1+, laminin+ RV/tubule structures in ectopic locations, and a lack of Sall1-GFP+ CM.
Investigating earlier development, we observed numerous differentiating structures (PTAs or RVs) that formed on the dorsal side of the ureter (ectopic location) in the E12.5 Sall1GFP mutants that expressed Sall1, Lef1, NCAM1/2, R-cadherin (Cdh4), laminin and Jag1 (Fig. 3A). Control E11.5 heterozygous kidneys matched for the number of UB tips (stage matched) had a small number of condensing cells forming PTAs. Quantification of Lef1+NCAM+ structures with a morphologically distinct lumen in control and mutant stage-matched kidney (T-stage to six UB tips) revealed that the mutant had significantly more differentiating structures (0.15±0.06 versus 1.42±0.22, P=1.1×10-8). We analyzed the RV/tubule structures by immunofluorescence at E13 using markers of RVs (NCAM, Lef1, Jag1, Lhx1 and laminin; Fig. 3B). E13 wild-type kidneys had a highly structured pattern of ureteric bud tips surrounded by organized CM cells with RVs forming in the ‘arm pit’ of the ureter. These RVs were NCAM+, Lef1+, Lhx1+, Jag1+ and laminin+ and formed on the ventral side of the ureter, towards the interior of the kidney. In the mutant, some renal vesicles were located ectopically on the dorsal side of the ureter, towards the periphery of the forming kidney. Whereas the peripheral nephrogenic zone of the wild type was populated by CM cells, the mutant kidney lacked Pax2+ CM cells in the nephrogenic zone. RV formation was abnormal, with no apparent UB/RV pattern, and RVs formed without apparent contact with the UBs. The RV structures in the mutant expressed NCAM, Lef1, Lhx1, Jag1 and laminin, and at E13 the mutant lacked more differentiated comma and S-shaped bodies. The vesicle structures also expressed Wt1, indicative of a podocyte lineage; however, the mutant RV morphology was abnormal with no organized podocyte layer.
We performed whole-mount in situ hybridization for renal vesicle markers Pax8, Lef1, Lhx1, Fgf8, Ccnd1 and Wnt4 (Fig. 3C). We found that gene expression in renal vesicles was irregular in the mutant, and the expression of some genes (Pax8 and Lef1) was higher in larger foci compared with the controls. Quantitative real-time PCR analysis (qRT-PCR) using stage-matched kidneys revealed that similar to the findings with in situ analysis, Wnt4 expression was not altered in the mutant, but other genes expressed in RVs such as Fgf8, Jag1, Mecom, Ccnd1, Pax8, Lama4 and Npy (Barak et al., 2012; Brunskill et al., 2008; Grieshammer et al., 2005; Karner et al., 2011; Thiagarajan et al., 2011; Yu et al., 2012), were increased in the mutant (Fig. 3D), supporting the evidence for unrestrained differentiation in the Sall1 mutant.
We investigated whether these ectopic RVs could differentiate normally in the Sall1GFP mutant (supplementary material Fig. S1). At E16.5 the mutant had very few mature differentiated structures; however, we detected some Cdh6+lotus tetragonolobus lectin (LTL)+ proximal tubules in ectopic locations toward the periphery of the kidney in the nephrogenic zone. We also detected Wt1+ structures that were probably developing glomeruli, although these were also peripherally located and the cells making up the presumptive podocyte layer were poorly organized.
We also observed unrestrained differentiation in the Sall1GFP mutant kidney when cultured ex vivo in organ culture. At the T-stage, mutant kidneys displayed GFP+ CM cells around each UB tip, similar to stage-matched control kidneys (Fig. 4A). After 96 hours, the heterozygous kidney had branched many times with GFP+ CM cells surrounding the UB tips. Mutant rudiments exhibited fewer branches than heterozygous controls. However, we consistently observed that the GFP+ cells in the mutant formed condensed foci, presumably induced CM cells, that were larger and more numerous than those observed in the heterozygote (Fig. 4A). To investigate whether the condensed GFP+ cells were in fact differentiated RVs, we immunostained organ-culture kidneys for Lhx1. After 72 hours in culture, the heterozygous kidneys had few Lhx1+GFP+ structures, whereas the mutant had multiple differentiating structures (Fig. 4Ba-f). Many RVs in the mutant were ectopically located toward the periphery of the kidney and there were multiple RV structures surrounding a single UB tip (Fig. 4Bd). Quantification of Lhx1+GFP+ RV formations per UB tip revealed significantly more RV formations per UB tip in the mutant compared with the control (Fig. 4Bg). These organ-culture studies indicated that the CM is formed around UB tips at E12.5, but it is progressively lost through differentiation in the Sall1GFP mutant.
Renal progenitor cells are depleted in Sall1GFP mutants
Taken together, the lack of Pax2+ CM cells and overabundance of RV structures led us to hypothesize that the progenitor cell pool in the Sall1 mutant kidney was depleted because of the onset of rapid nephrogenesis. Progenitor cell markers Six2, Cited1, Bmp7 and Osr1 were significantly reduced in the mutant, indicating that the progenitor cells were depleted by E13 (Fig. 5A). This was confirmed by immunofluorescence, by which we observed markedly reduced expression of Six2, and Cited1 expression was undetectable (Fig. 5B). qRT-PCR analysis revealed that Cited1, Osr1, Bbx, Fgf20 and Six2 were significantly reduced (≥50%) in the mutant, and Gdnf, Eya1, Bmp7 and Myc were moderately reduced in the stage-matched mutant in comparison with the heterozygous control (Fig. 5C).
To determine whether reduced proliferation could account for progenitor cell depletion, we performed cell cycle analysis on sorted GFP+ cells (Fig. 5D). In the Sall1-GFP+ cells, the cell cycle was not significantly altered at E12.5 between control and mutant. However, at E13.5 the mutant had a higher percentage of cells in the G0-G1 stage (45.5±3.3 versus 70.7±1.4, P=0.02), and a decrease in the percentage of cells in the S and G2-M stages. Similar cell cycle alterations were observed at E14.5 in the mutant. We confirmed this result using pHH3 immunofluorescence to determine the percentage of cells in the G2-M phase of the cell cycle (Fig. 5E). We observed a significant decrease in the percentage of cells in the G2-M phase of the cell cycle only at E13.5, confirming our cell cycle results using the GFP-sorted cells. TUNEL staining revealed that there was no significant apoptosis in either the heterozygous control kidney or the mutant kidney at E12.5. Significant apoptosis was observed at E13.5 and E14.5 in the mutant, although a majority of the cells undergoing apoptosis in the mutant were located in RVs (Fig. 5F).
At E12.5, when we observed ectopic differentiation in the mutant, we did not observe an alteration in the cell cycle or an increase in apoptosis. Changes in cell cycle and apoptosis only became apparent at E13.5 and E14.5. At these stages the majority of the CM in the mutant was depleted, leaving stroma and differentiating structures (Figs 2, 3; Fig. 5E,F). These data, together with the findings of ectopic RV formation (Figs 3, 4), support the conclusion that Sall1 mutant nephron progenitors are biased towards differentiation, leading to their depletion.
To investigate the early (E11.5-12.5) mechanism of kidney development in the absence of Sall1, we performed whole genome transcriptional profiling using GFP+ cells from Sall1GFP heterozygous and mutant E11.5 and E12.5 kidney. We precisely matched the stage of development between the heterozygote and the mutant to compare gene expression based on developmental stage of the kidney and embryonic day. Sall1GFP heterozygous mice are indistinguishable from wild type, thus could serve as a suitable control, and enabled sorting of GFP+ cells. Samples collected were as follows: E11.5 heterozygous kidneys at the T-stage to 6 UB tip stage; E12.5 mutant kidneys at the T-stage to 6 UB tip stage; and E12.5 heterozygous kidneys with 12-18 UB tips (Fig. 1F).
Using the FACS strategy previously described to isolate Sall1GFP Hi multipotent progenitor cells (Osafune et al., 2006), ∼100,000 GFP+ cells from more than 350 embryos for each sample type were collected (Fig. 6A). We performed qRT-PCR to confirm that the sorted GFP+ pool reflected that of Sall1-expressing cells. The GFP+ population was enriched with the Sall1GFP transcript, and CM progenitor cell markers Six2 and Cited1. Stromal markers Raldh2 (Aldh1a2 - Mouse Genome Informatics) and Foxd1 were present at very low levels in the GFP+ population. The UB marker Ret and endothelial cell markers Flk1 (Kdr), VE-Cadherin (Cdh5) and CD45 (Ptprc) were enriched in the GFP-negative population (Fig. 6B). The precise morphological and developmental staging of this assay allowed us to analyze gene expression profiles for the normal developmental progression between the initial invasion of the UB at E11.5 and the beginning of differentiation at E12.5, as well as those genes that were altered in the mutant.
When we compared expression profiles of the E11.5 and E12.5 heterozygous samples, we found very similar global distributions, confirming that the overall composition of the CM differs very little at these stages (Fig. 6C). However, when we compared expression profiles of the E11.5 heterozygote and E12.5 mutant that were matched for developmental stage (Fig. 6D), we observed a significant deviation in the distribution of many genes. This was intriguing because although the control and mutant samples were specifically matched for stage, there were significant gene expression changes occurring in the mutant.
The stage-matched comparison between control and mutant GFP+ cells revealed 248 upregulated genes and 46 downregulated genes, by at least twofold. This finding is consistent with the predominant role of Sall1 as a transcriptional repressor through its association with the NuRD complex (Kiefer et al., 2002; Lauberth et al., 2007; Lauberth and Rauchman, 2006). Among the statistically significant changes were genes previously reported to be altered in Sall1 mutants (Itga8, Gdnf, Kif26b) (Nishinakamura et al., 2001; Uchiyama et al., 2010), supporting the validity of our data. We confirmed the gene expression changes observed in the microarray with qRT-PCR on a number of genes (supplementary material Fig. S2A). Table 1 lists the top 25 genes and their fold change. Many genes identified with altered expression in the mutant are associated with renal syndromes (Table 2), indicating that the data may reveal new genes important for kidney development and birth defects. Gene ontology analysis using DAVID (http://david.abcc.ncifcrf.gov/) (Table 3) revealed that many of the genes with altered expression in the mutant are associated with the extracellular region, indicating that Sall1 may affect genes expressed in the MM that signal to the ureter, surrounding CM or stromal cells. Biological adhesion, cellular component organization, locomotion and growth were among the biological processes affected with high statistical significance in the mutant. Pathway analysis revealed highly statistically significant values for cadherin signaling, axon guidance mediated by Slit-Robo, and Wnt signaling, which are all pathways important in the development of the kidney (Carroll et al., 2005; Grieshammer et al., 2004; Karner et al., 2011; Kispert et al., 1998).
When we analyzed the fold changes in our samples and focused on genes that were altered ≥1.5-fold in the mutant, we uncovered a significant trend (Fig. 6E). The control sample exhibited a subset of genes that had increased or decreased expression between E11.5 and E12.5. The same genes in the mutant had magnified gene expression changes. Genes that were upregulated in expression between E11.5 and E12.5, as the kidney developed, in the control progenitor pool had higher expression in the mutant, whereas genes that decreased in expression in the control progenitor pool as the kidney developed had significantly lower expression in the mutant (Table 1). This suggested that the normal developmental changes that occurred at the initiation of kidney development were amplified in Sall1-deleted cells, confirming our previous observations in vivo of a sensitized differentiation program.
Many of the genes that had decreased expression in the microarray analysis in the mutant are expressed in the CM (Cited1, Osr1, Meox1, Meox2, Bbx, Six2, Eya1, Bmp7, c-Myc, Fgf20, Fgf9, Crym, Hoxa11 and Gdnf). Importantly, this includes genes that have a defined functional role in maintaining or expanding self-renewing progenitors (Barak et al., 2012; Blank et al., 2009; Couillard and Trudel, 2009; Sajithlal et al., 2005; Self et al., 2006; Tena et al., 2007). This result indicated that the genetic program that defines self-renewing nephron progenitors is generally downregulated early in kidney development in the mutant. We also compared our genome-wide data to the list of genes that change in RVs relative to CM, identified by microarray analysis using laser capture microdissection of these cell populations (Brunskill et al., 2008). We found that there were 63 genes that displayed a statistically significant change (at least 1.5-fold) in expression, similar to induction of RV genes during normal development. We confirmed the bioinformatic analysis with qRT-PCR of whole embryonic kidney for a subset of these genes (supplementary material Fig. S2B; Fig. 3D). Strikingly, similar to the results of our array analysis, the change in gene expression (CM versus RV) was in the same direction as in the mutant; however, the fold change was greater, suggesting that Sall1 may be acting to control the dynamic range of gene expression responsible for the conversion of CM to RV (supplementary material Fig. S2B; Table 1).
Initial studies of Sall1 mutants focused on the impaired outgrowth and initial branching of the UB, the earliest step where the gene is required in the developing kidney (Nishinakamura et al., 2001). However, in a majority of Sall1 mutants, the UB invades the MM and initiates nephrogenesis. This observation, together with the finding that Sall1+ cells are multi-potent progenitors (Osafune et al., 2006), led us to investigate whether Sall1 has a direct role in regulating the fate of renal progenitor cells. Here, we present data revealing a crucial role for Sall1 in the progenitor cell niche, regulating the transition between the self-renewing CM and the differentiation of these cells into RVs, the first epithelial structure in the developing nephron. Genetic deficiency of Sall1 leads to a reduction of CM gene expression, whereas the differentiation program is accentuated. These gene expression changes are accompanied by depletion of Cited1+Six2+ progenitor cells and evidence of ectopic renal vesicles in the peripheral nephrogenic zone.
Recent studies have identified Fgf9, Fgf20 and Bmp7 as important secreted factors, which together, promote survival, proliferation and competence of CM to respond to inductive signals (Barak et al., 2012; Blank et al., 2009; Oxburgh et al., 2004). In Bmp7-null mice, progenitor cells are mainly lost through apoptosis soon after nephrogenesis initiates, but premature differentiation is not observed (Blank et al., 2009; Dudley et al., 1995). At E12.5 in Sall1 mutant kidney, Bmp7 expression was modestly (<20%) reduced (Fig. 5C), and we did not observe cell cycle alterations or apoptosis at this stage (Fig. 5D-F). Yet, at this stage, there was evidence of premature differentiation in Sall1 mutants, suggesting that reduced Bmp7 signaling is probably not contributing significantly to loss of CM. By contrast, Fgf9/Fgf20 compound mouse mutants display a phenotype which is reminiscent of Sall1-/- (Barak et al., 2012). Our microarray data showed a ∼60% reduction in Fgf20 expression at E12.5 in Sall1-/- mice, raising the possibility that once the MM is specified and invaded by the UB, Sall1 may regulate Fgf20 expression to maintain self-renewing progenitor cells in a differentiation-competent state.
Our transcriptional profiling in Sall1-/- renal progenitors showed many more genes were upregulated than downregulated. We have previously demonstrated that Sall1 associates with the nucleosome remodeling and deacetylase (NuRD) chromatin remodeling complex to repress gene transcription (Lauberth et al., 2007; Lauberth and Rauchman, 2006), suggesting a molecular mechanism for Sall1-mediated inhibition of gene expression in the kidney. Consistent with this finding, we recently showed that conditional deletion of the NuRD-specific component Mi2b (Chd4) in Six2+ cells leads to the depletion of CM (Denner and Rauchman, 2013). Moreover, Sall1 and Mi2b display a strong genetic interaction in the developing kidney, supporting a cooperative role for Sall1 and NuRD in maintaining renal progenitor cells. Investigation of lymphocyte differentiation suggests that proper genomic targeting of the Mi2-NuRD complex requires interaction with lineage-specific DNA-binding factors, such as Ikaros, to balance growth and differentiation of progenitors (Zhang et al., 2011). In embryonic stem cells (ESCs), Sall4 and NuRD are found in a macromolecular complex that has been shown to regulate pluripotency and lineage determination (Liang et al., 2008). A striking finding of our genome-wide analysis is that normal developmental gene expression is magnified in the Sall1-/- mutant (Fig. 6E; Table 1). This is reminiscent of recent studies in which NuRD was shown to modulate transcriptional heterogeneity and the dynamic range of expression of pluripotency genes in ESCs (Reynolds et al., 2012). These findings raise the intriguing possibility that Sall1 and NuRD might function together to balance growth and differentiation in the kidney by similarly setting the threshold of expression of key CM genes, such as Six2 and Fgf20 (Karantzali et al., 2011; Liang et al., 2008; Reynolds et al., 2012), that in turn, determine the decision of a cell to self-renew or differentiate.
Six2 deficiency leads to a phenotype that is similar to that observed in Sall1 mutants, in which there is ectopic RV formation and loss of progenitor cells (Self et al., 2006). In the Six2 mutant, Sall1 expression is maintained in the CM (Self et al., 2006); however, we observed reduction in Six2 in the Sall1 mutant (Fig. 5). The fact that Six2 expression is significantly decreased, but not undetectable in Sall1 mutants, may explain why premature differentiation is apparent at an earlier stage (E11.5 versus E12.5) in Six2-/- compared with Sall1-/- kidneys. We hypothesize that Sall1 is upstream of Six2 and might be needed to maintain an appropriate level of expression in the CM so that at each UB tip only a subset of cells are committed to differentiation, with the remainder undergoing self-renewal. Although Sall1 has been defined as a transcriptional repressor, it is likely to have co-activator properties through its association with NuRD (Zhang et al., 2011) or other co-factors. Future studies to define Sall1 genomic binding sites will provide insight into whether Six2 is a direct transcriptional target of Sall1.
Wnt9b signaling through the canonical β-catenin pathway induces CM to differentiate, but it is also required to promote proliferation of self-renewing progenitor cells. It has been proposed that interpretation of the Wnt9b signal in CM is dependent on the level of Six2 and Cited1 (Brown et al., 2013; Carroll et al., 2005; Karner et al., 2011) to define the cellular response to β-catenin transcriptional input. A current model suggests that relatively high levels of nuclear β-catenin together with reduced Six2 in Cited1-negative nephron progenitor cells favors differentiation (Brown et al., 2013; Park et al., 2012). By contrast, increased levels of Six2 in a transcriptional complex with Tcf/Lef, prevents differentiation induced by β-catenin. In the Sall1 mutant, Wnt9b expression is upregulated and expands to the tip of the Sall1 mutant UB (Kiefer et al., 2010) (Fig. 1E); therefore, it may serve as the signal to the CM that promotes unrestrained differentiation. Consistent with this notion, constitutively active β-catenin expressed in CM results in ectopic RV formation, indicating that at high levels this signal is sufficient for nephron induction (Kiefer et al., 2010; Park et al., 2012). Our previous studies showed that increased levels of Wnt9b, similar to those observed in Sall1 mutants (Kiefer et al., 2010) (Fig. 1E), was sufficient to increase expression of the RV genes Wnt4, Fgf8 and Pax8; yet this level of Wnt9b expression did not induce ectopic RVs (Kiefer et al., 2012). This suggests that in the absence of a very strong canonical Wnt signal (as seen with the stabilized β-catenin gain-of-function allele), the renal progenitor cell cannot fully overcome barriers that restrain differentiation. In the Sall1 mutant, the combination of reduced Six2 expression and elevated Wnt9b levels extending to the UB tip may synergize to cause ectopic RV formation. Recent studies suggest that stromal signals cooperate with Wnt9b to inhibit nephron progenitor expansion and promote expression of differentiation genes (Das et al., 2013). Future studies will determine whether Sall1-dependent signals from the stroma affect nephron progenitor cell fate. Alternatively, Sall1-deficient cells could bypass the requirement for a Wnt inductive signal as shown for constitutively active Notch (Boyle et al., 2011).
Together, our data show that Sall1 is a nuclear factor that is required to maintain the stemness of renal progenitor cells by restricting their differentiation. Loss of Sall1 biases the cellular response to Wnt9b toward differentiation rather than self-renewal. We propose a model whereby Sall1 functions to coordinate UB branching with nephron induction (Fig. 7). In the CM, there is an intrinsic requirement for Sall1 to influence the responsiveness of the nephron progenitors to a β-catenin signal to determine whether the output is self-renewal or differentiation, possibly through the modulation of Six2 expression levels. However, it is also possible that the Sall1-NuRD complex influences the β-catenin transcriptional output at the level of chromatin to determine the cell fate decision in CM, as NuRD has been shown to inhibit canonical Wnt signaling (Major et al., 2008). Sall1-dependent factors from the MM also regulate the expression of Wnt9b to control the strength and domain of the inductive signal from the UB. Together, these two activities of Sall1 ensure that as the UB branches and induces the CM to differentiate, the population of progenitor cells is not depleted and remains adequate to produce a normal complement of nephrons.
MATERIALS AND METHODS
Sall1GFP mice have been previously described (Takasato et al., 2004) and bred onto the ICR background for at least ten generations. Female mice were checked for the presence of a mucosal plug, and noon on the day of detection was considered to be E0.5. These studies were performed according to the St. Louis University and VA St. Louis Health Care System animal care guidelines.
Whole-mount in situ hybridization, histology and immunohistochemistry
In situ hybridization was performed as previously described (Kiefer et al., 2010). Primer sequences used to produce the probes are listed in supplementary material Table S2. Embryonic kidneys were embedded in paraffin, sectioned at 4 μm, and stained using Harris’ Hematoxylin and Eosin Y (Sigma-Aldrich) for histology. For immunofluorescence, sections were boiled in 10 mM citric acid, pH 6, for antigen retrieval, and incubated with the following primary antibodies: rabbit anti-Sall1 (1:1000, Rauchman Lab), mouse anti-NCAM (1:500, Sigma-Aldrich), mouse anti-cytokeratin (1:1000, Sigma-Aldrich), mouse anti-Lhx1 (1:50, Developmental Studies Hybridoma Bank, University of Iowa), rabbit anti-Wt1 (1:1000, Santa Cruz), rabbit anti-Pax2 (1:1000, Covance), rabbit anti-CDH6 (1:1000, Sigma-Aldrich), rabbit anti-Lef1 (1:250, Cell Signaling), rabbit anti-CITED1 (1:800, Thermo Fisher Scientific), rabbit anti-SIX2 (1:200, Lifespan Biosciences), mouse anti-pHH3 (1:100, Millipore), rabbit anti-Jag1 (1:250, Santa Cruz), rabbit anti-laminin (1:500, Sigma-Aldrich), rat anti-R-cadherin (1:5, Developmental Studies Hybridoma Bank, University of Iowa), and biotin anti-LTL (1:100, Vector Laboratories). Antibody reactivity was detected using Cy3-labeled anti-mouse IgG (1:2000, Jackson ImmunoResearch), Alexa-Fluor-488-labeled anti-rabbit IgG (1:400, Life Technologies), Cy3-labeled anti-rat IgG (1:2000, Jackson ImmunoResearch), and fluorescein streptavidin (1:100, Vector Laboratories). Antigen retrieval was not used for sections co-stained with the R-cadherin and Pax2 primary antibodies (Fig. 3A). Without antigen retrieval Pax2 staining is strong in the ureter and very weak in the CM and differentiating structures. Sections were imaged using an Olympus 41BX upright epifluorescence microscope.
Quantification of UB branching
Embryonic kidneys were immunostained and the number of cytokeratin+ UB tips was counted on at least five non-sequential sections (20× magnification) from at least two independent embryos for each stage and genotype. Results are reported as the average number of tips per section ± s.e.m. Statistical analysis using standard t-tests was performed.
Quantitative real-time PCR (qRT-PCR)
Total RNA was isolated from sorted cells using an RNeasy Micro Kit or from embryonic kidney tissue using an RNeasy Mini Kit according to the manufacturer’s protocol (Qiagen). cDNA was prepared using the High Capacity RNA-to-cDNA kit (Life Technologies). Primer sequences are listed in supplementary material Table S1. qRT-PCR was performed using an ABI 7300 Real-Time PCR Thermocycler and SYBR Green PCR Master Mix (Life Technologies) as described previously (Kiefer et al., 2012). Real-time reactions were performed in triplicate and relative expression was calculated using the delta CT method and normalized to Rpl19 or Hprt1 control transcripts (Kiefer et al., 2012).
Cell cycle analysis
Embryonic kidneys were disaggregated using 50 μl 0.05% Trypsin-EDTA for 5-10 minutes on ice. Cells were fixed for 5 minutes at room temperature with 0.25% PFA, pelleted, washed with PBS, and incubated with 75% ethanol for 10 minutes. Cells were then stained in propidium iodide (PI; 50 μg/ml PI, 0.1 mg/ml RNase A, 0.1% Triton X-100 and 0.1 mM EDTA) for 1 hour at room temperature, and the cell cycle was analyzed using a FACS Calibur and ModFitLT V3.0 (Win32) software.
Mitotic and apoptotic index
Mitotic index was determined by staining embryonic kidneys for pHH3 and Sall1. Nuclei were stained using DAPI. The ratio of (pHH3+Sall1+)/(total Sall1+) cells was calculated. At least 800 nuclei for each stage and genotype were counted on at least five non-sequential sections (20× magnification) from at least two independent embryos. The apoptotic index was determined by performing TUNEL analysis using the ApopTag Red In Situ Apoptosis Detection Kit (Millipore) according to the manufacturer’s instructions. The ratio of (TUNEL+Sall1+)/(total Sall1+) cells was calculated. At least 900 nuclei for each stage and genotype were counted on at least five non-sequential sections (20× magnification) from at least two independent embryos. Results were reported as the mitotic/apoptotic ratio ± s.e.m. Statistical analysis using standard t-tests was performed.
Embryonic kidney organ culture
Embryonic kidney culture was performed as previously described (Barak and Boyle, 2011) with the following exceptions: kidneys were cultured on 0.1 μm filters (Millipore) supported by a metal grid, in DMEM and 10% FBS. The kidneys were imaged daily using a Leica M165 FC fluorescent stereo microscope to monitor development. Explants were stained overnight at 4°C with mouse-anti Lhx1 antibody (1:10, Developmental Studies Hybridoma Bank, University of Iowa) and for 2 hours at room temperature with Cy3-labeled anti-mouse IgG (1:800, Jackson ImmunoResearch) secondary antibody. Kidneys were imaged with an Olympus FV1000MPE confocal microscope.
Sall1GFP cell sorting and transcriptional profiling
E11.5 and E12.5 kidneys were disaggregated and sorted on a BD FACSARIA III. GFP+ cells were collected by using wild-type kidneys as a GFP-negative control. Using this method we achieved 93% purity of the GFP+ population of cells by post-sorting analysis. Approximately 100,000 sorted cells were collected into three independent pools for each sample and RNA extracted using Qiagen RNeasy Micro Kit. RNA integrity and concentration were tested using the Agilent Bioanalyzer on an RNA Pico chip. Total RNA (5-10 ng) was used for probe synthesis using NuGen WT-Ovation Pico + Exon Assay and the probe hybridized to Affymetrix GeneChip Mouse Gene 1.0 ST arrays. Probe production and hybridization were performed at the Cincinnati Children’s Hospital Gene Expression Microarray Core. Array normalization, fold change, and statistical analysis were performed using R (http://www.r-project.org/) and Bioconductor (http://www.bioconductor.org/). Array data have been submitted to ArrayExpress with accession number E-MTAB-2159.
The authors thank the Flow Cytometry Core Facility and the Microscopy Core at St. Louis University for their services. We appreciate the assistance of Steven Potter and Shawn Smith at the Microarray Core Facility, Cincinnati Children’s Hospital, Cincinnati, OH. Ryuichi Nishinakamura generously provided Sall1GFP mice. We also thank Darcy Denner and Jake Basta for thoughtful reading of the manuscript and Akbar Siddiqui for valuable assistance with the bioinformatic analysis of the microarray data.
J.M.B., S.M.K. and M.R. developed concepts, performed data analysis, and prepared and revised the manuscript. J.M.B., S.M.K. and L.R. performed experiments. D.D. assisted with microarray data analysis.
This work was funded by the American Heart Association and the National Institute of Diabetes and Digestive and Kidney Diseases. Deposited in PMC for release after 12 months.
The authors declare no competing financial interests.