Microarrays and RNA-Seq identify molecular mechanisms driving the end of nephron production
© Brunskill et al; licensee BioMed Central Ltd. 2011
Received: 8 October 2010
Accepted: 12 March 2011
Published: 12 March 2011
Skip to main content
© Brunskill et al; licensee BioMed Central Ltd. 2011
Received: 8 October 2010
Accepted: 12 March 2011
Published: 12 March 2011
The production of nephrons suddenly ends in mice shortly after birth when the remaining cells of the multi-potent progenitor mesenchyme begin to differentiate into nephrons. We exploited this terminal wave of nephron production using both microarrays and RNA-Seq to serially evaluate gene transcript levels in the progenitors. This strategy allowed us to define the changing gene expression states following induction and the onset of differentiation after birth.
Microarray and RNA-Seq studies of the progenitors detected a change in the expression profiles of several classes of genes early after birth. One functional class, a class of genes associated with cellular proliferation, was activated. Analysis of proliferation with a nucleotide analog demonstrated in vivo that entry into the S-phase of the cell cycle preceded increases in transcript levels of genetic markers of differentiation. Microarrays and RNA-Seq also detected the onset of expression of markers of differentiation within the population of progenitors prior to detectable Six2 repression. Validation by in situ hybridization demonstrated that the markers were expressed in a subset of Six2 expressing progenitors. Finally, the studies identified a third set of genes that provide indirect evidence of an altered cellular microenvironment of the multi-potential progenitors after birth.
These results demonstrate that Six2 expression is not sufficient to suppress activation of genes associated with growth and differentiation of nephrons. They also better define the sequence of events after induction and suggest mechanisms contributing to the rapid end of nephron production after birth in mice.
In humans, the final number of nephrons that are produced during the formation of the kidney is extremely variable, ranging from 230,000 to 1,800,000 . The number has clinical relevance because a decrease in the number has been associated with hypertension . Although our knowledge about the molecular control of nephron formation has grown substantially in recent years, little is known about the mechanisms controlling the final number.
During renal development, multi-potential progenitors surround the branch tips of the ureteric bud. These self-renewing progenitors are maintained in an undifferentiated state at least in part by Six2 , and are induced to differentiate into nephrons by Wnt9b . In mice and humans the lifespan of this population of progenitors is limited providing a means to regulate the final endowment of nephrons. Even though little is known about the mechanism controlling the population's lifespan, it must be able to account for some important differences between humans and mice at the end of nephron production. In humans, for example, branching morphogenesis of the ureteric bud ends early at 22 weeks of gestation , whereas in mice it ends within the three days after birth . In humans, nephron production extends twelve weeks beyond the period of ureteric bud branching to 34 weeks gestation, whereas in mice production ends around the same time as branching morphogenesis ends, about three days after birth . Thus, the relationship of the completion of these two processes, branching morphogenesis and nephron production, markedly differ in humans and mice, suggesting the possibility of some small difference in the mechanism controlling the completion of nephron production.
To understand the molecular mechanisms that drive the completion of nephron production in mice, we defined the gene expression programs of the progenitors of nephrons during the first four days after birth. At this time, all remaining progenitors progress from a primarily un-induced to an induced state, and then form renal vesicles. This final wave of production offers a unique opportunity to evaluate the progenitors because of the near synchronous change in their behavior. Our results indicate that induction, defined by altered gene expression, occurs before significant decreases in Six2. We also show that proliferation increases prior to detectable increases in transcripts of genes associated with differentiation. Further, we show that expression of some genes, previously defined as markers of later developmental stages, is present in capping mesenchyme cells co-expressing Six2. Finally, we observed altered expression levels of genes encoding proteins in the glycolytic pathway, consistent with a change in the microenvironment of the population of progenitors after birth. The microarray results were independently validated and expanded by using a next generation deep sequencing RNA-Seq approach. The resulting profiles better define the order of events and the genes involved after induction of progenitors, as nephron production comes to an end.
The difference in the level of expression of GFP between capping mesenchyme and the renal vesicles allowed us to collect cells of the cap at birth, P1, and P2, and the developing renal vesicles at P3 and P4. We have shown that the multi-potential progenitors, which form the cap and express Cited1, disappear by three days of age in the mouse, and that nephron production ceases at that time . Thus, the cells that are captured at the later times are no longer considered self-renewing progenitors. By gating the FAC sorting to collect only the most highly expressing GFP-positive cells at each time (See additional file 3: Image of the FACS plot of cells collected for RNA measurements:), we could exploit this phase of development to determine the changing gene expression states during induction and early differentiation. The sudden loss of the high-expressing, un-induced mesenchyme after birth and the continued lower-level expression of GFP in the early developing nephron made it possible to collect a series of samples beginning with the un-induced progenitors, then induced progenitors, and finally cells in the early stages of differentiation. This period after birth, which contrasts the embryonic period when the dominant population of GFP-positive cells is un-induced, permits an enriched population of recently induced cells to be isolated.
Taking advantage of the series of arrays, we next examined a subset of developmental times to detect genes showing very early changes in transcript levels, representing the initial response to either induction or to the extra-uterine environment. By comparing P0 samples with P2 using unpaired t-test, we identified 70 genes (p < 0.05) with at least a 1.5-fold change (Figure 2a). Fifty genes were more highly expressed in the GFP(+) cells at birth than at P2 and included Tgfbi, Dock5, Robo2, and Mfap4. The expression patterns of these four genes in the cap have historic validation in the embryonic kidney by in situ hybridization . Expression was down regulated earlier than a significant down-regulation in expression of Cited1 (Figure 2b.), a gene that has been used to mark un-induced cells. The expression of this marker decreased by P2, but it decreased an additional four- to five-fold in the subsequent two days. The genes we identified therefore represent a very early response after birth.
Twenty genes were more highly expressed in the cap at P2 than at P0. They included Cdh6, Bmp2, Ccnd1, Pax8, Hey1, Lhx1, Npy, and Jag1. They too have historically validated embryonic expression, but in the early stage nephrons, rather than in the mesenchyme. Arrays also identified developmental control genes, such as Rspo3, a secreted regulator of beta-catenin signaling , and Fat3 (see discussion), two genes that have not been studied during renal development. It is interesting to note at P2, when capping mesenchyme is still abundant, that many of the genes showing higher expression compared to P0 are reported in the literature to be associated with more advanced stages of differentiation and not with the mesenchyme. These results could represent either contamination with differentiated cells or low-level expression of vesicle genes in the induced cap prior to vesicle formation.
We also used RNA-Seq to independently measure differences in gene expression between progenitors at P0 and P2. Deep sequencing of cDNA libraries lacks the bias related to microarrays, and it can be used to digitally quantify gene expression, even at the exon level, over a wide dynamic range. We aligned the sequences to the mm9 mouse sequence database. This is a subset of Ref-seq, which is an integrated, non-redundant database of linked nucleotide and protein sequences. The total number of single-end sequence reads from P0 and P2 was 3.1 and 3.3 million reads, respectively. Alignment with mitochondrial sequences gave an additional 0.65 and 0.4 million reads, respectively. Reads aligned to 19,168 mouse genes and counts ranged from 1 to almost 150,000. Interestingly, and providing a measure of the quality of the analysis, we found that the percent decrease (37%) in expression of GFP at P2 was similar to the decrease (49%) in expression of Crym, the gene controlling GFP expression.
RNA-Seq reads at P0 and P2 for genes expressed during renal development
P < 10-93
P < 10-200
Of interest, in our analysis of the microarray data we found that eight genes encoding enzymes in the glycolytic pathway had a significantly reduced level of expression (P < 0.05) at P2 compared to P0. Five of these were down regulated at least 1.5 fold (Pfkl, Aldoa, Tpi1, Eno3 and Pdk1). Three of these genes (Pfkl, Aldoa, and Tpi1) clustered together (Figure 2), indicating similar expression patterns after birth. The pattern of change was evident by 24 hours after birth. The change preceded the increase in Wnt4 levels, a marker of renal vesicles, and preceded the decrease in Cited1 levels (Figure 2 and See additional file 4). RNA-Seq validated this decrease in level of expression of seven genes (Gpi1, Pfkl, Aldoa, Tpi1, Eno3, Pkm2, and Pdk1). In contrast, the level of expression of genes encoding enzymes of the Kreb's cycle (Aco2, Idh2, Ogdh, Scla2, Sdha, Fh1, and Mdh1) remained constant in the interval between P0 and P2 by both microarray and RNA-Seq. We analyzed microarray data from the GUDMAP consortium  to identify changes in expression of glycolysis genes during the process of differentiation of cap mesenchyme in the embryo. When we compared the cap mesenchyme expression of E15.5 embryos to renal vesicles of E12.5 embryos, we could not detect similar changes in expression of the glycolysis genes.
The genome-wide evaluation of the multi-potent progenitor cell mesenchyme and early stage nephron transcriptome after birth showed sequential activation or inactivation of many genes. The use of tissue after birth to obtain expression data eliminates potential artifacts introduced during in vitro studies of sequential activation after induction. The period after birth allows serial measurements of transcripts early after induction because there is such an abrupt change in the behavior of the progenitors. By P3 all remaining capping mesenchyme has begun conversion into renal vesicles. Because the inductive period and the time to show signs of differentiation after induction each take 12-24 hours (reviewed in Saxen ), the final wave of induction to form nephrons had begun by P2. This strongly suggests, therefore, that the Cited1(+) cells at P2, which still constitute most of the GFP(+) capping mesenchyme (See additional file 1), were likely to have been induced and to be different from the capping mesenchyme on P0 and P1 when the bulk of the high GFP expressing cells was un-induced. There is clearly a progression of molecular events post-induction, with the down-regulation of Cited1 corresponding to a more advanced state of induction or possibly a state of commitment.
Our results confirmed differences in expression level between capping mesenchyme and renal vesicle found among a series of markers previously reported by Mugford et.al. . Seven of the nine markers, which distinguished between progenitor cell mesenchyme and vesicles by in situ (Cited1, Bbx, Eya1, Osr1, Six2, Dpf3, and Meox1), showed greater than a 2-fold change in level of expression between P0 and P4. The trend was correct for the remaining two markers (Hoxc5 and Brpf1), but did not reach a 2-fold change in level. Our data extend this list significantly by adding more than 400 genes that change in level between P0, when mesenchyme is abundant, and P4, when GFP primarily marks renal vesicles. This provides a significant resource that can be used when monitoring induction of the mesenchyme.
Discovery of genes that have a spatially restricted pattern of expression during the process of differentiation helps to identify pathways that may be needed for nephron formation. We found that Fat3 was expressed in the cap at P0 and was down regulated by P2. The Fat genes encode cadherin-type proteins with cell-cell adhesion properties. With mutation of Fat genes in Drosophila, there is overgrowth of tissues . It is interesting in Drosophila that Fat is part of a pathway involved in the suppression of wingless, a homolog of a murine gene (Wnt4) that is expressed after mesenchymal induction and is required for nephron development . This suggests that the down-regulation of Fat3 may be needed before activation of Wnt4 in the induced mesenchyme.
We also identified early up-regulated expression of genes in several different pathways and examined them by in situ hybridization. One gene, Bmp2, was activated early in the induced capping mesenchyme next to one side of the ureteric bud branch tips. The secreted protein encoded by Bmp2 is known to regulate both the branching of the ureteric bud tips and the proliferation of their cells . With the localized expression, Bmp2 may regulate the regional growth of cells within the tips. Its expression is also of interest because along with Clu and Lama4 it is activated in the induced cap before morphologic changes are evident. It seems reasonable that capping mesenchymal cells initiating renal vesicle formation would begin expressing genes associated with renal vesicles, although to our knowledge this has not been previously demonstrated.
It is interesting that these induced genes, which are associated with differentiation, are co-expressed with Six2. Six2 is necessary to maintain the population of multi-potential progenitors; however, it has not been shown to be sufficient. The co-expression of genes involved in differentiation in a subset of Six2 (+) cells suggests that Six2 is not sufficient to block transcription of some genes associated with differentiation in the presence of an inductive signal.
We also observed an interesting surge in cell proliferation that preceded the expression of markers of differentiation. The sequence of events suggests the possibility that proliferation promotes the reprogramming of renal progenitors. This type of mechanism has been described during reprogramming of somatic cells into pluri-potent stem cells . An increase by P2 in proliferation, evidence of a change in the behavior of progenitors, is also consistent with induction of the cells by P2.
Also of note, the gene expression profiles of early forming nephrons differed before and after birth. The genes included those encoding enzymes of the glycolytic pathway. The shift in transcription is compatible with a response by the population of progenitors to a change in the microenvironment, such as a post-natal increase in oxygen levels. Prior to birth the kidney is fed by deoxygenated blood, blood with the same oxygen content as that returning to the placenta. In addition, oxygen delivery to the nephrogenic region is further limited because the tissue is relatively avascular. Relative physiologic hypoxia, such as this, is known to cause an increase in transcription of genes encoding enzymes of the glycolytic pathway . Changes in levels of expression of other genes that are regulated in an oxygen-dependent manner, such as P4ha1, Bnip3L, and Txnip, provide further supportive evidence of an increase in oxygenation in the progenitors after birth.
Cellular fates of placental cytotrophoblasts , hematopoietic progenitors , human neural stem cells , bone marrow stromal cells , and human embryonic stem cells  have been shown to be altered by oxygen. The fate of murine embryonic stem cells also appears to be coupled to metabolism . It seems reasonable, therefore, to speculate that a change in the microenvironment, such as a change in oxygenation within the physiologic range, might also lead to a change in behavior of the multi-potential progenitors in the kidney in vivo.
Lastly, the final nephron endowment is clearly a result of regulation of the balance between the rate of progenitor renewal and the rate of differentiation. Simple geometry might play an important role in the balance in both mice and humans. In the early kidney, the capping mesenchyme layer is relatively thick. As branching morphogenesis proceeds, the number of branch tips will expand geometrically, subdividing the capping mesenchyme while also inducing it. Unless renewal of the progenitors similarly expands, the cap around each tip will thin. And, at some point it will no longer be able to promote further branching. After branching ends, if usage of cells to make nephrons exceeds the renewal rate, nephron production will consume the remaining cap mesenchyme. In humans this simple model fits nicely with the completion of nephron production because there is a prolonged period of nephron production without branching. In the mouse, however, there is an abrupt end to both nephron production and branching morphogenesis. The end coincides with birth and with a change in metabolism that is compatible with an increase in oxygenation of the progenitors. Coupled with the known effects of oxygen on cellular fate, the events suggest a possible trigger in mice at birth that shifts the balance between renewal and differentiation of progenitors favoring differentiation. We speculate that the trigger then limits the lifespan of the population of progenitors and causes the production of nephrons in mice to end abruptly.
In this study we used microarrays, RNA-Seq, in situ hybridization, and EdU nucleotide incorporation to examine the synchronous wave of nephron formation that occurs in mouse following birth. The results provide a global definition of the changing gene expression program that drives the transition from un-induced capping mesenchyme, to induced capping mesenchyme, and to renal vesicle. Several genes were found to change in expression before Cited1 was down regulated, suggesting a further molecular subdivision of induced capping mesenchyme. We also observed that cell proliferation preceded differentiation. Further, we observed the expression of renal vesicle associated genes, including Bmp2, within the Cited1(-), but Six2(+), region of the capping mesenchyme. This demonstrates that Six2 at the normal level of expression alone is not sufficient to block differentiation. Finally, genes essential for glycolysis were down regulated post-birth, compatible with an increase in levels of oxygen within the capping mesenchyme domain. Given the known influence of oxygen on the fate of cells this finding suggests a possible trigger that promotes the final burst of nephron production.
Tg(Crym-EGFP)82Gsat/Mmcd mice with GFP expression in nephron progenitors were used. To obtain developmentally uniform tissue, mice were housed in standard light-dark cycles, matings were started at midnight, and mice checked for vaginal plugs eight hours later. At noon on day 18 (e18.5) embryos were delivered by Caesarian section, resuscitated, and placed with a foster mother. Each litter of pups was used to obtain kidneys for multiple post-natal samples. Kidneys were dissected at birth (P0) and at 24-hour intervals afterwards (P1-P4). We used CD-1 mice for in situ hybridization and nucleotide analog incorporation. The use of experimental animals described in this study complies with the guidelines of the Institutional Animal Care and Use Committee at Cincinnati Children's Hospital.
Two pairs of kidneys were isolated for each sample and three samples isolated for each time. The kidneys were digested in 0.05% Trypsin - EDTA for 15 minutes at 37°C, treated with 2% FBS, and mechanically disaggregated at 4°C. The suspension of cells was filtered through a 70-micron filter and separated using a BD FACS Aria II Cell Sorter to collect GFP (+) cells with the highest fluorescence intensity (See additional file 3: Image of the FACS plot of cells collected for RNA measurements).
Fifteen samples, three for each age, were isolated by cell sorting and collected in RLT (Qiagen, CA). 5 ng of total RNA was used with the NuGen WT-Pico, Exon Module, and Fl-Ovation v.2 kits to generate target. 2.5 μg of target was hybridized to Affymetrix Mouse Gene 1.0 ST arrays. The data was deposited in the Gene Expression Omnibus database (GSM429020-GSM429034). Analysis was performed using GeneSpring GX and RMA normalization. COMBAT software was used to correct for batch effect . ToppGene Suite was used to analyze gene lists to identify functional groups of genes . Cytoscape was used to create Figure 6 and 7. To compare post-natal changes in expression of glycolysis genes during differentiation of progenitors to changes in the embryo we analyzed microarray data from the GUDMAP consortium which had been deposited in the Gene Expression Omnibus (GEO, GSE12588 and GSE6290 datasets).
A dual-labeled fluorescent whole mount in situ hybridization procedure was used with confocal microscopy to validate microarray results. Digoxigenin- or fluorescein-labeled riboprobes were made for each of the test genes (Clu, Fat3, Lama4, Tgfbi, Bmp2) and for both of the reference genes (Wnt11 and Six2). Perkin Elmer tyramide signal amplification and the Zeiss LSM 510 microscope were used to detect the hybridization signal. A series of 2-micron thick optical sections were obtained beginning at the surface and extending into the kidney at 5 micron intervals.
We collected GFP-positive progenitors by fluorescence activated cell sorting, purified the RNA, and synthesized cDNA from 150 ng of pooled RNA from each time, P0 and P2, using the Clontech SMART cDNA synthesis kit. cDNA was then amplified using Stratagene Herculase polymerase during a 24-cycle PCR reaction. The PCR products were purified and treated with the BAL-31 nuclease to remove the SMART cDNA primer end, and sequenced using the Illumina Genome Analyzer II standard 36-cycle paired-end protocol. Sequences were aligned to the mus musculus (mm9 sequence database) subset of RefSeq  in an unpaired alignment using the ELAND . The alignment could contain only two mismatches of a 32 bp read to be included in the datasets. GeneSpring was used to graphically present the raw RNA-Seq data and Partek Genomic Suites for statistical analysis.
We used the Invitrogen Click-iT EdU Alexa Fluor 488 to label cells in S-phase of the cell cycle. Pups were injected with EdU, a nucleotide analog, at 24 or at 48 hours of age. The kidneys from pups were excised 60-90 min later, fixed, and immunostained with antibodies to E-cadherin as described before . After blocking in 3% BSA in PBS, and washing in 0.1 M sodium phosphate buffer pH 8, they were incubated in complete Click-iT reaction buffer for 90 min at room temperature, washed and imaged by confocal microscopy.
Support was provided by grants from the Kidney Foundation of Greater Cincinnati to HL, and DK061916 to SS.
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.