Abstract
Interkinetic nuclear migration (INM) is a key feature of cortical neurogenesis. INM functions to maximize the output of the neuroepithelium, and more importantly, balance the self-renewal and differentiation of the progenitors. Although INM has been reported to be highly correlated with the cell cycle, little is known about the effects of cell cycle regulators on INM. In this study, by crossing Foxm1fl/fl mice with Emx1-Cre line, we report that a conditional disruption of forkhead transcription factor M1 (Foxm1) in dorsal telencephalon results in abnormal cell cycle progression, leading to impaired INM through the downregulation of Cyclin b1 and Cdc25b. The impairment of INM disturbs the synchronization of apical progenitors (APs) and promotes the transition from APs to basal progenitors (BPs) in a cell-autonomous fashion. Moreover, ablation of Foxm1 causes anxiety-related behaviors in adulthood. Thus, this study provides evidence of linkages among the cell cycle regulator Foxm1, INM, and adult behavior.
Introduction
The neuroepithelium gives rise to a great number of neurons and glial cells during mammalian neocortical development. Both neuroepithelial cells and radial glial cells (RGCs), which form later, are defined as apical progenitors (APs) because these progenitors are localized at the apical surface of the neuroepithelium. APs exhibit a bipolar shape and share an evolutionarily conserved feature: interkinetic nuclear migration (INM). During the cell cycle, the nuclei of the APs move toward the basal side of the ventricular zone (VZ) during G1 phase, undergo DNA synthesis (S phase), migrate rapidly to the apical surface during G2 phase, and then undergo mitosis (Götz and Huttner, 2005; Taverna and Huttner, 2010; Willardsen and Link, 2011). The nuclei of APs at various positions along the apical-basal axis in the VZ give rise to the remarkable pseudostratification of the mammalian neuroepithelium. This feature is thought to maximize the output of the neuroepithelium (Taverna and Huttner, 2010) and to be important for tissue-shape formation (Reiner et al., 2012).
Genetic and pharmacological studies have revealed that the cell cycle is required for INM, although INM is dispensable for normal cell-cycle progression. Furthermore, nuclear migration remains coordinated with the cell cycle when the cell cycle is modified (Murciano et al., 2002; Schenk et al., 2009; Kosodo et al., 2011). Subsequent studies demonstrated that INM is also highly correlated with neurogenesis (Fietz and Huttner, 2011; Willardsen and Link, 2011; Kosodo, 2012). Additionally, INM is important for regulating proliferation versus differentiation (Murciano et al., 2002; Buchman and Tsai, 2008; Kageyama et al., 2009; Latasa et al., 2009; Willardsen and Link, 2011). Forkhead transcription factor M1 (Foxm1), a key cell cycle regulator, is a novel candidate that might link proliferation and differentiation (Yokomine et al., 2010; Kalin et al., 2011). Although the functions of Foxm1 in cell-cycle progression have been well studied (Krupczak-Hollis et al., 2004; Costa, 2005; Laoukili et al., 2007; Schüller et al., 2007; Fu et al., 2008; Radhakrishnan and Gartel, 2008), the link between the cell-cycle regulator Foxm1 and INM has not yet been examined, and the possible effects of Foxm1 on behavior remain unknown. In the present study, we observed that the conditional disruption of Foxm1 in the telencephalon results in an impairment of INM and abnormal cell-cycle progression; these effects are accompanied by a decrease in the number of APs, which may be due to the downregulation of Cyclin b1 and Cdc25b. Furthermore, impaired basal-to-apical INM disturbs AP synchronization, resulting in an increased number of basal progenitors (BPs) in the developing cortex; this leads to a slight reduction in the number of cortical neurons, and the Foxm1 mutants exhibited decreased anxiety-related behavior. Our study demonstrates that the cell-cycle regulator Foxm1 is required for INM, and more importantly, deficiency in Foxm1 has effects on adult behavior.
Materials and Methods
Generation of mice with dorsal telencephalic-specific Foxm1 deletion
Foxm1fl/fl mice (Kalinichenko et al., 2002) were a kind gift from Dr Maureen Gannon (Vanderbilt University Medical Center) and Dr Pradip Raychaudhuri (University of Illinois at Chicago). Emx1-Cre mice were purchased from The Jackson Laboratory. Exons 4–7 of Foxm1 are excised upon exposure to Cre recombinase (Gorski et al., 2002), generating the recombinant Foxm1 allele. A dorsal telencephalon-specific disruption of Foxm1 was accomplished using Emx1-Cre, which drives recombination in both progenitor and projection neurons. Genotyping of the floxed allele was performed by genomic PCR using primer 1 (5′-TAGGAGATACACTGTTATAT-3′) and primer 2 (5′-TGTGGGAAAATGCTTACAAAAG-3′), which generated a 210 bp floxed band. Because Cre-mediated recombination also occurred in the olfactory bulb, the Foxm1 recombined allele was assayed using olfactory bulb DNA with primer 1 and primer 3 (5′-CTCATGTAGCATAGAGGGCTG-3′), which generated a 510 bp recombined band (Kalinichenko et al., 2002). Foxm1fl/fl; Emx1-Cre mice are referred to as Foxm1 cko/mutant mice. The day that a vaginal plug was detected was defined as embryonic day (E) 0.5. All animals were bred in the animal facility at Southeast University. All experiments were performed according to guidelines approved by Southeast University.
In situ hybridization and double labeling with immunohistochemistry
Brains from E10.5 to adult mice were dissected, fixed in 4% DEPC-PFA overnight, cryoprotected in 30% sucrose/DEPC-PBS at 4°C and embedded in OCT. Coronal sections (10 μm thick) were obtained using a Leica CM 3050S cryostat and stored at −70°C until use. For Nocht1, Delta-like 1, Hes5, Cyclin b1, and Cdc25b probe preparation, total RNA was extracted from E14.5 brains, and 400–600 bp fragments were subcloned into pBluescript SK vectors and linearized (details available upon request). The Foxm1 probe was a kind gift from Dr Ulrich Schüller (Ludwig Maximilian University; Schüller et al., 2007). In situ hybridization was performed as previously described (Tian et al., 2012). For double labeling with Tbr2 immunostaining, in situ hybridization was performed first; after the color reaction with NBT and BICP, the sections were refixed in 4% PFA for 30 min, and immunostaining was then performed per standard procedures. Images of in situ hybridization were captured by differential interference contrast microscopy (FV1000; Olympus), and pseudocoloration was performed using the FV10-ASW software (Olympus, v 02.01.02.05).
Immunohistochemistry
For immunohistochemistry, brains from E12.5 to adult mice were fixed in 4% PFA, cryoprotected in 30% sucrose at 4°C, and sectioned to 10 μm using a Leica CM 3050S cryostat. Sections were washed in PBS, blocked with 10% normal goat serum, permeabilized in PBS containing 0.1% Triton X-100 (PBT) for 2 h and incubated in primary antibody diluted in blocking solution overnight at 4°C. Sections were subsequently washed in PBT and incubated in secondary antibody for 2 h at 37°C. Finally, coverslips were applied. The primary antibodies used for immunostaining were mouse anti-Ki67 (Leica, 1:100), rabbit anti-Ki67 (Abcam, 1:500), rabbit anti-Tuj1 (Covance, 1:1000), rabbit anti-Tbr2 (Abcam, 1:500), rat anti-Ctip2 (Abcam, 1:1000), rabbit anti-Foxp2 (Abcam, 1:1000), rat anti-PH3 (Abcam, 1:1000), mouse anti-PH3 (Cell Signaling Technology, 1:500), rabbit anti-Pax6 (Covance, 1:200), mouse anti-Nestin (DSHB, 1:200), rabbit anti-β-Catenin (Sigma-Aldrich, 1:1000), mouse anti-N-Cadherin (Invitrogen, 1:1000), and mouse anti-ZO1 (Invitrogen 1:1000). AlexaFluor 555 donkey anti-mouse IgG (Invitrogen; A31570, 1:500), AlexaFluor 633 goat anti-mouse IgG (Invitrogen; A21050, 1:500), AlexaFluor 488 goat anti-rabbit IgG (Invitrogen, 1:500), AlexaFluor 555 donkey anti-rabbit IgG (Invitrogen, 1:500), and AlexaFluor 488 goat anti-rat IgG (Invitrogen, 1:500) were used as secondary antibodies. DAPI was purchased from Sigma-Aldrich.
BrdU labeling and S- phase calculation
BrdU (Sigma-Aldrich) and IdU (Sigma-Aldrich) were dissolved in physiological saline at concentrations of 10 and 5 mg/ml, respectively. Nuclei in the early S, late S, G2/M, and G1 phases were labeled by intraperitoneal injection of BrdU (50 mg/kg) at 0.5, 1.5, 5, or 14 h before harvest. Antigen retrieval was required for the immunostaining of BrdU and IdU with the rat anti-BrdU (Abcam, 1:1000) and mouse anti-IdU/BrdU (BD Bioscience, 1:150) primary antibodies. The protocol for S-phase compartment calculation was performed as previously described (Martynoga et al., 2005).
Cumulative BrdU labeling
Cumulative BrdU labeling in E14.5 embryos was performed by repeated intraperitoneal injections of BrdU into pregnant mice according to the following schedule. The first injection was given between 9:00 A.M. and 10:00 A.M. For the 0–8 h time points, there was a 2 h interval between each injection; for the 16 h time points, there were 4 h intervals, and for the 24 h time points, there were 4 h intervals. The mice were killed 2 h (0–8 h time points) or 4 h (16 and 24 h time points) after the last injection (Nowakowski et al., 1989; Calegari et al., 2005). Antigen retrieval was required for the immunostaining of BrdU and IdU with the rat anti-BrdU (Abcam, 1:1000) and mouse anti-IdU/BrdU (BD Bioscience, 1:150) primary antibodies. Nuclei were stained with DAPI.
Measurement of cell-cycle parameters and statistical analyses
For each time point, the BrdU labeling index was obtained from two independent litters, and three pairs of brains from these two litters were sectioned for immunohistochemistry; at least four nonconsecutive cryosections of the medial cortex were imaged with a confocal microscope (FV1000; Olympus) using a 40× objective for cell counting. Nuclei in the VZ were identified by DAPI staining, and the cells were counted in a 50-μm-wide field in each section. Mathematical analysis was performed as previously described (Nowakowski et al., 1989).
Quantitative real-time PCR
Total RNA was isolated from E14.5 dorsal telencephalons using the RNeasy Plus Mini Kit (Qiagen) according to the manufacturer's instructions. Each sample was reverse transcribed using MultiScribe reverse transcriptase (Fermentas). The quantitative PCRs were performed using SYBR Green fluorescent master mix (Roche) on a StepOne Plus Real-Time PCR System (Applied Biosystems). Primers were synthesized as shown in Table 1.
Primers for quantitative real-time PCR
Microscopic analysis
Sections were examined with an Olympus BX61 microscope equipped with the appropriate filter sets and a digital camera (DP71) or with a confocal microscope (FV1000; Olympus). Differential interference contrast (DIC) images were obtained from the confocal microscope (FV1000, Olympus) in our laboratory, and we changed the color of the DIC images into red fluorescence using the FV10-ASW software (Olympus, v 02.01.02.05).
Behavioral tests
Mice were 7 months of age when the tests began. The tests were performed in the following sequence: open-field test, elevated zero maze (O-maze) and fear conditioning.
Open-field test.
For the open-field test, each mouse was placed into a square chamber (40 × 40 cm). The movement of the mice was recorded using a Unibrain Fire-i digital camera. Data were collected for 15 min. The chamber was cleaned between each trial.
Elevated O-maze.
For this test, the equipment was composed of a 6-cm-wide ring with an outer diameter of 45 cm containing two equal walled (closed) sections and two unwalled (open) sections. The entire ring was elevated to a height of 100 cm. Mice were placed at the boundary of a walled section and an unwalled section, facing the unwalled section. The time spent in the open and closed arms and the total number of entries into the open arm were recorded for 6 min. The open- field test and elevated O-maze results were analyzed using Noldus EthoVision 8.5.
Fear conditioning.
Fear conditioning tests were performed following a standard training protocol (Ugo Basile). Mice were placed in a Plexiglas shock chamber and were allowed to explore the chamber for 150 s. A 2800 Hz tone was then sounded for 30 s as the conditioning stimulus; during the last 2 s of the tone, a 0.8 Ma footshock was delivered. All tests were performed 24 h after training. For contextual conditioning, mice were placed into the same chamber used for the training test, and freezing times were recorded for 4 min. For cued conditioning, mice were placed into a novel chamber for 3 min; the test lasted 3 min after the tone used in the training test was played. The freezing time during the pre-tone and post-tone periods was recorded.
Statistical analysis
Student's t tests were performed to compare region- and layer-specific changes in the developing cortex in absolute and relative percentages of cell numbers. For the BrdU labeling, four pairs of mice were obtained from three different litters; for the cumulative BrdU labeling, three pairs of mice were obtained from two different litters for each time point. In total, 24 pairs of mice from 16 litters were used. For the P8 and adult immunohistochemistry analyses, three pairs of mice were obtained from three different litters; for the in situ hybridization, five control and four FoxM1 cko mice were obtained from three different litters. Three pairs of mice were obtained from three different litters for the real-time PCR analysis. For the embryonic experiments, the gender of the mice was not determined. For the behavior tests, 19 control mice (10 males, 9 females) and 14 Foxm1 mutants (7 males, 7 females) were obtained from 12 litters. Statistical tests were performed with GraphPad 5.0. Differences were considered significant when p < 0.05.
Results
The conditional ablation of Foxm1 in the developing telencephalon results in a decrease in cell proliferation and an abnormal AP pool
To elucidate the function of Foxm1 in the developing telencephalon, we first examined Foxm1 expression. Foxm1 mRNA is expressed at high levels throughout the neuroepithelium (E10.5–E16.5; Fig. 1A–D). Interestingly, high levels of Foxm1 expression were detected in the area of overlap between the VZ and subventricular zone (SVZ; Fig. 1B–D, insets). We further performed in situ hybridization for Foxm1 combined with immunostaining for Tbr2, a marker of SVZ BPs, at E14.5. We detected Foxm1 mRNA in APs within the VZ but did not detect colocalization with Tbr2 (Fig. 1E). This expression pattern suggests possible roles for Foxm1 in cell proliferation and differentiation. Next, we conditionally disrupted Foxm1 in the developing dorsal telencephalon by crossing Foxm1fl/fl mice with an Emx1-Cre line. As expected, the level of Foxm1 mRNA was markedly reduced in the dorsal telencephalon (Fig. 1F,G); the expression of Foxm1 in the ventral telencephalon was unchanged (Fig. 1F).
Expression patterns of Foxm1 during embryonic development and conditional disruption of Foxm1 in the telencephalon. A–D, In situ hybridization of Foxm1 in coronal sections of E10.5–E16.5 telencephalons. Foxm1 mRNA was detectable at E10.5 (A) and covers the entire neuroepithelium (arrow). From E12.5 onward, high levels of Foxm1 mRNA were detected in the overlapping area between the VZ and SVZ (B–D, insets and arrows). E, Double staining for Foxm1 mRNA by in situ hybridization and for Tbr2 by immunostaining. No colocalization of Foxm1 mRNA and Tbr2 was observed in the basal side of the VZ. F, No Foxm1 mRNA was detected in the Foxm1 cko dorsal telencephalon (asterisk), whereas normal expression in the basal forebrain was not affected. G, Relative mRNA level of Foxm1. The removal of FoxM1 was demonstrated by quantitative real time-PCR (Foxm1fl/fl N = 5, Foxm1 cko N = 4, p < 0.001). H–J, Immunohistochemistry for 1.5 h BrdU labeling (H, I) showed that cell proliferation was decreased. The arrows indicate the range of VZ and SCZ. J, A reduction (∼16.7%) in BrdU+ proliferating cell density was observed in the dorsal cortex of the Foxm1 cko mice (Foxm1fl/fl N = 5, Foxm1 cko N = 5, p = 0.0146). Pax6 staining (K, L) revealed that the AP pool was affected in the Foxm1 cko mice. M, A reduction (12.9%) in the Pax6+ AP density was observed in the Foxm1 cko mice (Foxm1fl/fl N = 4, Foxm1 cko N = 4, p = 0.0005). The data are presented as the mean ± SEM. Scale bar, 50 μm.
Previous studies have shown that Foxm1 plays an important role in mitotic progression (Wang et al., 2011; Zhang et al., 2011), suggesting that Foxm1 might be involved in the proliferation of neuronal progenitor cells. To test this hypothesis, we quantified the proliferative capacity of the neuronal progenitors using 1.5 h BrdU labeling. We detected a decreased number of BrdU-labeled progenitors in Foxm1 mutants (Fig. 1H–J). Furthermore, immunostaining for Pax6, a marker of APs in the VZ, was performed on E14.5. Foxm1 mutants showed a 12.9% decrease in the number of Pax6+ progenitors compared with the control (Fig. 1K–M). Together, these results suggest that Foxm1 is required for the proliferation of neuronal progenitors and the maintenance of the AP pool.
Disruption of Foxm1 impairs INM in APs
Based on the strong Foxm1 expression observed in the VZ and its role in regulating the cell cycle, we suspected that Foxm1 might play a role in INM. Because the total length of the cell cycle at E14.5 is ∼14–16 h, including a 4 h S phase, 2 h G2/M phase, and 8–10 h G1 phase, we performed BrdU labeling at time points from 0.5–14 h to track the locations of the nuclei at different phases of the cell cycle. The labeled nuclei of APs that were observed 30 min after BrdU injection were located in the basal zone of the VZ in both the control and Foxm1 mutants (Fig. 2A–A″), consistent with the known position of AP nuclei during S phase (Latasa et al., 2009). No obvious differences were detected at this point, except that the nuclei of mutant APs appeared to be more condensed than those of control APs. Surprisingly, none of the 1.5 h labeled nuclei in the mutants reached the apical surface, although a large portion of the labeled nuclei in the control animals did (Fig. 2B–B″). To more carefully measure nuclear migration, the entire cortex was divided into eight bins, and the number of nuclei in each bin was calculated. We found that in control cortices, 17.6 ± 2.9% of the nuclei migrated into bin 1; however, the percentage of basal-to-apical migrating nuclei in bin 1 was markedly decreased in the mutant (5.3 ± 1.7%; Fig. 2B′,B″). Meanwhile, the percentage of BrdU-labeled nuclei remaining in bin 3 was increased in the mutants (41.8 ± 2.1%) compared with the controls (33.0 ± 1.4%; Fig. 2B,B″). These results indicate that basal-to-apical nuclear migration is disturbed by the deletion of Foxm1. We then performed 5 h BrdU labeling. In the control animals, a large proportion (40 ± 3.3%) of BrdU-labeled nuclei had accumulated in bin 1, the apical zone, and only a small proportion (8.2 ± 0.6%) remained in bin 4, the basal area of the VZ (Fig. 2C,C″). In mutant animals, however, the nuclei were dispersed throughout the VZ, with 31 ± 1.4% reaching bin 1 and 13 ± 1.2% remaining in bin 4. There was an obvious slowing of basal-to-apical nuclear migration (Fig. 2C′,C″). Similar defects were also observed for the 1.5 h and 5 h BrdU labeling at both E12.5 and E16.5 (Fig. 2E–H′). Moreover, a 14 h BrdU labeling was performed to further analyze whether differences in migration exist at later time points in the cell cycle. A similar basal enrichment of BrdU-labeled nuclei was observed in both the controls and mutants, but the control animals displayed a population of BrdU-labeled nuclei at the apical surface, representing the first group of APs that had completed one cell cycle and entered a second G2 phase. In the Foxm1 mutants, all of the BrdU-labeled G2 phase nuclei remained located in the upper area, and none of them had arrived at the apical surface (Fig. 2D–D″). Together, these results indicate that Foxm1 is required for INM, especially basal-to-apical migration.
Impairment in the basal-to-apical nuclear migration of neural progenitors. A–D, Immunostaining for BrdU in coronal sections of the dorsal telencephalon in E14.5 after 0.5 h (A, A′), 1.5 h (B, B′), 5 h (C, C′), and 14 h (D, D′) BrdU injection. The arrows indicate the most apical parts of the ventricular zone. The entire cortex was divided into 8 bins, as indicated, with bin 1 at the VZ and bin 8 at the marginal zone. The percentage of BrdU+ nuclei in each bin relative to the total number of BrdU+ cells was measured (A″, B″, C″, D″; Foxm1fl/fl N = 4, Foxm1 cko N = 4, p < 0.001). E, F, Coronal sections through the E12.5 dorsal telencephalon. Immunostaining of nuclei labeled at 1.5 h (E, E′) and 5 h (F, F′) after BrdU injections. G, H, Coronal sections through the E16.5 dorsal telencephalon. Immunostaining of nuclei labeled at 1.5 h (G, G′) and 5 h (H, H′) after BrdU injections. Similar abnormal basal-to-apical migration was observed at both E12.5 and E16.5. The arrows in E–H′ indicate BrdU labeled nuclei that have arrived at the apical surface. The two dashed lines in each figure indicate the apical and basal surfaces of the VZ. The data are presented as the mean ± SEM. Scale bar, 50 μm.
INM is impaired, but not arrested, and APs are able to pass through M phase in Foxm1 mutants
Previous studies have shown that a chemical blockade of INM disperses mitotic nuclei from their normal apical locations to random sites throughout the VZ during M phase (Murciano et al., 2002; Del Bene et al., 2008). To investigate whether Foxm1 gene deletion would have similar consequences due to its arrest of INM, we performed a time course analysis with double staining for BrdU and phospho-histone H3 (PH3), an M phase marker that labels APs undergoing mitosis at the apical surface. After a 1.5 h BrdU injection, PH3+ nuclei were found at the apical surface in both the controls and mutants, and no randomly dispersed mitotic nuclei were observed in the VZ in the mutants, which indicated that INM was not blocked (Fig. 3G,H). In control animals, ∼89.7 ± 3.7% of PH3+ nuclei were colabeled with BrdU (Fig. 3A–A‴,C), indicating that these progenitors were poised to complete DNA synthesis at the onset of BrdU injection, passed through G2 phase and underwent mitosis within 1.5 h. As expected, no PH3 and BrdU double-labeled nuclei were observed in the Foxm1 mutants, as a result of impaired basal-to-apical nuclear migration (Fig. 3B–B‴,C). After a 5 h BrdU labeling, a smaller subset of BrdU-labeled nuclei reached the apical surface in the mutants, indicating that the basal-to-apical migration was slower, but all of the PH3+ nuclei were found to be colabeled with BrdU in both the controls and mutants (Fig. 3D–F).
Ablation of Foxm1 impairs basal-to-apical migration but does not arrest the cell cycle at M phase. A, B, D, E, Sections stained for anti-BrdU and PH3 antibodies. A′–A‴, B′–B‴, D′–D‴, and E′–E‴ are enlarged figures of dashed outlines in A, B, D, and E, respectively. In the Foxm1 cko mice, after the 1.5 h BrdU labeling, none of the BrdU+ nuclei reached the apical surface or were colabeled with PH3 (A, B). By 5 h after BrdU labeling, the BrdU+ nuclei had arrived at the apical surface and were colabeled with PH3 in both the wild-type and Foxm1 cko mice (D, E). C, F, Percentage of PH3+ cells at the ventricular surface that are BrdU+PH3+ (Foxm1fl/fl N = 3, Foxm1 cko N = 3, p = 0.017). The arrows indicate BrdU+PH3+ nuclei, the arrowheads indicate BrdU−PH3+ nuclei. Immunostaining for PH3 and DAPI (G–J). The progenitors from both wild-type and Foxm1 cko mice exhibited normal M phase nuclei. The arrowheads indicate prophase nuclei. Asterisks indicate metaphase nuclei. The arrows indicate anaphase nuclei. The data are presented as the mean ± SEM. Scale bar, 50 μm.
Previous studies showed that Foxm1 is required for spindle assembly and centrosome replication (Costa, 2005; Wierstra and Alves, 2007; Xue et al., 2010). Foxm1-deficient hepatocytes (Wang et al., 2002) and myocardial cells (Ramakrishna et al., 2007) exhibit enlarged, polyploid nuclei as a result of mitotic arrest. To investigate whether mitotic progression was disrupted, we performed double staining for DAPI and PH3. Identical results were obtained from the controls and Foxm1 mutants. In prophase, the nuclei exhibited dispersed chromosomes, as shown by DAPI staining. The chromosomes were captured by spindle fibers, aligned at the equatorial plate and colabeled with PH3 during metaphase (Fig. 3I,J,I″,J″), and the chromatids were separated during anaphase (Fig. 3I‴,J‴). These data show that at least some Foxm1-deficient APs are competent to go through M phase, although there might be problems in M phase.
Cell-cycle parameters are abnormal, and the synchronization of APs during S and G2 phases is disturbed
The normal length of one cell cycle at E14.5 is ∼14–16 h. We performed a 14 h BrdU labeling to investigate whether Foxm1 disruption affects the total cell-cycle length. Double immunostaining for BrdU and Ki76, a marker of proliferating cells in all active phases of the cell cycle, showed that the fraction of BrdU+Ki76+ progenitors among all BrdU+ cells was 29.0 ± 1.3% in the control animals but 37.9 ± 2.1% in the mutants (Fig. 4A–C), indicating that a greater number of progen-itors remained in the cell cycle (Fig. 4B–B′) and that the cell cycle length might be longer in Foxm1 mutants than in controls. Because INM is dependent on cell cycle progression, we measured the S-phase duration using IdU and BrdU double labeling to precisely analyze the length of S phase, as previously described (Martynoga et al., 2005). The mice were first treated with IdU to label the progenitors in S phase. After 1.5 h, the mice were treated with BrdU. Thus, the progenitors that remained in S phase would be colabeled by BrdU and IdU, whereas those that had already left S phase and entered G2 would be BrdU−IdU+ (Fig. 4D). By calculating the length of S phase (Ts) using the formula
Abnormal cell cycle parameters and disrupted synchronization of APs in the Foxm1-ablated telencephalon. A, D, H, Schematic of BrdU and IdU injections. B, E, I, J–O, Coronal sections of E14.5 telencephalon. B, B′, Double staining of anti-BrdU and Ki67 at 14 h after BrdU injection. C, Percentage of BrdU+ cells that are BrdU+Ki67+. More BrdU+ and Ki67+ APs remained in the cell cycle in the Foxm1 cko mice than in the wild-type mice (Foxm1fl/fl N = 6, Foxm1 cko N = 5, p = 0.0002). E, E′, Double staining of BrdU and IdU at 2 h after IdU and 0.5 h after BrdU injections. The arrows indicate BrdU−IdU+ G2-phase cells. F, Foxm1 cko mice displayed a ∼2.7 h longer S phase than the wild-type mice (Foxm1fl/fl N = 3, Foxm1 cko N = 3, p = 0.0014). G, Number of BrdU+ APs that exited S phase within 1.5 h. In the Foxm1 cko mice, fewer S phase cells (∼25.6%) entered G2 phase during the 1.5 h interval (Foxm1fl/fl N = 3, Foxm1 cko N = 3, p = 0.0385). I, I′, Double staining of BrdU and IdU after the 5 h IdU and 1 h BrdU injection. The arrows indicate BrdU−IdU+ G2-phase cells. The INM-mediated synchronization of APs was disturbed. J–O, Double staining of BrdU and PH3 after the 1.5 h (J, K), 3 h (L, M) and 5 h (N, O) injection. The BrdU+ cells were not colabeled with PH3 until 5 h after BrdU injection. The data are presented as the mean±SEM. The arrows indicate the forefront nuclei. P, P′, Identification of VZ by DAPI staining. The solid line in P and P′ indicate the boundary between VZ and SVZ; the dashed line indicate apical surface. Q, BrdU labeling index and cumulative BrdU labeling experiments were performed at 0.5 h, 2 h, 4 h, 6 h, 8 h, 10 h, 20, and 28 h. The regression equations were y = 0.0715x + 0.2527, 2–10 h r2 = 0.974 (control) and y = 0.048x + 0.2345, 2–10 h r2 = 0.984 (Foxm1 mutants). The value of TS, TC and TC − TS were determined as previously described (Nowakowski et al., 1989). Scale bar, 50 μm.
The increased proportion of BrdU+/Ki67+, together with the delayed apical-to-basal INM at the 14 h BrdU injection, raises the possibility that not only S phase, but also other phases of the cell cycle might also be affected in the Foxm1 mutants. To precisely measure the cell-cycle parameters, cumulative BrdU labeling was performed (Fig. 4P,P′), and the proportion of BrdU+ cells among all VZ cells was counted at 0.5, 2, 4, 6, 8, 10, 20, and 28 h. As expected, the proliferative population in the VZ grew at a steady rate in both the control and Foxm1 mutant mice; however, the proportion of BrdU+ cells in the control was ∼10–18% higher than that of Foxm1 at each time point until 20 h (Fig. 4Q). Using linear regression, we were able to acquire a regression curve and calculate the cell-cycle parameters using the regression equation. As shown in Table 2, the length of S phase increased by 38.2% to 4.9 h compared with 3.5 h in the control, whereas the total length of the cell cycle (Tc) was 6.8 h longer in the Foxm1 mutants (20.8 h) than in the control (14 h). Moreover, the value of Tc − Ts, which represents the time that a proliferating, newly BrdU-injected AP takes to pass through G2, M, G1, and re-enter S phase, increased from 10.5 h in the control to 16 h in the Foxm1 mutants, which indicated that the lengths of both G2/M and G1 phase, as well as S phase, were prolonged. These results indicate that the entire cell cycle was affected after Foxm1 ablation. The difference in the S phase lengths determined in the BrdU/IdU double labeling and the cumulative BrdU labeling experiments may result from the different mathematical and statistical methods used; the length of S phase was directly calculated in the BrdU/IdU double-labeling experiments but was indirectly calculated from the regression equation in the cumulative labeling experiments following regression analysis (Nowakowski et al., 1989; Calegari et al., 2005).
Cell cycle parameters of the APs
Another interesting phenomenon was that, in the control animals, a large portion of the BrdU−IdU+ G2 phase nuclei had separated from the BrdU+IdU+ S phase nuclei and arrived at the apical surface; however, in the Foxm1 mutants, the BrdU−IdU+ G2 phase nuclei remained in the basal area in a mosaic pattern mixed with S phase nuclei. This result indicates that the synchronization of APs was disturbed (Fig. 4E–E′). To observe this asynchronization more precisely, progenitors were first labeled with IdU for 4 h and colabeled with BrdU for another 1 h (Fig. 4H). Consistent with the previous experiments, the control BrdU+IdU+ S phase nuclei remained in the basal half of the VZ after 4 h, whereas most of the BrdU−IdU+ G2 phase nuclei had arrived at the apical surface; thus, nuclei in the S and G2 phases were clearly separated. In the Foxm1 mutants, a similar mosaic distribution of S and G2 phase nuclei was observed (Fig. 4I–I′). Moreover, the first group of IdU-labeled nuclei that entered G2 phase and arrived at the apical surface exhibited a mottled shape; these are referred to as “forefront” nuclei (Fig. 4E,I). In our model, double labeling of PH3- and BrdU-labeled forefront nuclei was also observed as early as the 1.5 h BrdU labeling (Fig. 4J), whereas in the Foxm1 mutants, double labeling was not detected until the 5 h BrdU labeling (Fig. 4K,M,O), indicating an abnormally slow G2 progression. Together, these results suggested that G2 phase might also be affected in the Foxm1 mutants, causing the observed defects in AP synchronization during basal-to-apical INM.
Foxm1 is important for maintaining the balance between AP self-renewal and the transition to BPs
The progenitor cell decision to undergo self-renewal or differentiation is one of the key programs in neurogenesis (Calegari et al., 2005; Dehay and Kennedy, 2007; Fietz and Huttner, 2011), and both INM and the cell cycle are involved in this process. During cortical neurogenesis, APs generate BPs, also called intermediate progenitors (Götz and Huttner, 2005), through asymmetric cell division (Zhong and Chia, 2008; Shitamukai and Matsuzaki, 2012). Our results have shown that Foxm1 regulates the cell cycle, INM, and AP synchronization. Therefore, we investigated whether the impairments in the cell cycle and INM affect the BPs. As shown by immunostaining for Tbr2, the deletion of Foxm1 resulted in an increase in the total number of BPs in the entire cortex (Fig. 5A–C), and the number of BPs in the VZ also increased (Fig. 5D–F). Tbr2+ BPs are normally generated through either BP self-renewal within the SVZ or the transition from APs in the VZ. Double immunostaining with BrdU (14 h injection) and Tbr2 showed that the percentage of Tbr2+ BrdU+ APs in the VZ, which corresponded to the APs that were poised to leave G1 phase and differentiate into BPs, was significantly increased in the mutants (Fig. 5G–K), indicating that the transition from APs to BPs was enhanced in the VZ. Interestingly, we also observed a twofold increase in the number of nonsurface mitotic PH3+ proliferating cells in the SVZ of the mutants (Fig. 5N). As expected, the number of Tbr2+PH3+ nonsurface mitotic BPs showed similar increase in Foxm1 cko (Fig. 5O). The increased number of Tbr2+ PH3+ cells may indicate that the BPs in the SVZ that were generated from Foxm1-deficient APs also exhibited a similarly slow cell cycle progression, and the prolonged cell cycle allowed more BPs to remain in the cell cycle before they migrated away from the SVZ (Fig. 5L,M). Our results indicate that Foxm1 is important for maintaining the balance between AP self-renewal and the transition to BPs; the disruption of Foxm1 promotes embryonic neurogenesis by generating more BPs.
Increased transition from APs to BPs. A, B, D, E, Immunostaining for anti-Tbr2. The Foxm1 cko mice exhibited a larger number of Tbr2+ cells in both the SVZ (A, B) and VZ (D, E). C, The total of Tbr2+ cells was increased by ∼19.3% in the Foxm1 cko mice compared with the wild-type mice (Foxm1fl/fl N = 6, Foxm1 cko N = 5, p = 0.0001). F, The density of the Tbr2+ cells in the VZ was increased 25.2% in the Foxm1 cko mice compared with the wild-type mice (Foxm1fl/fl N = 5, Foxm1 cko N = 5, p = 0.0048). G–J, Double staining for BrdU and Tbr2. I and J are enlarged figure of dashed outlines in G and H, respectively. There were more APs transitioning into BPs in the Foxm1 cko mice. K, The density of the BrdU+Tbr2+ cells in the VZ of the Foxm1 cko mice was ∼31.8% greater than in the wild-type mice (Foxm1fl/fl N = 3, Foxm1 cko N = 3, p = 0.0259). L, M, Double staining for Tbr2 and PH3. N, The density of the PH3+ cells in the SVZ (Foxm1fl/fl N = 7, Foxm1 cko N = 7, p < 0.0001). O, The density of the Tbr2+PH3+ cells in the SVZ (Foxm1fl/fl N = 5, Foxm1 cko N = 4, p = 0.0092). The density of both the PH3+ and Tbr2+PH3+ proliferating BPs in the SVZ was increased by ∼2-fold in the Foxm1 cko mice compared with the wild-type mice. The data are presented as the mean ± SEM. Scale bars: A, B, G, H, L, M, 50 μm; D, E, I, J, 25 μm.
The G2 phase regulators Cyclin b1 and Cdc25b are downregulated in the Foxm1 cko dorsal telencephalon
To identify possible molecular mechanisms in Foxm1 regulation of INM and the defects in cell cycle progression, we analyzed the mRNA levels of cell cycle regulators important for different cell-cycle phases, such as the S phase regulators Cyclin a and Cdk2, G2/M phase regulators Cyclin b1 and Cdc25b, and G1 phase regulators Cyclin d1 and Aurora a and b (Laoukili et al., 2007, 2008; Wierstra and Alves, 2007; Kalin et al., 2011). Interestingly, there was a 70% downregulation of the mRNA level of the G2/M-phase regulator Cyclin b1 (Fig. 6A). Moreover, another important regulator, Cdc25b, which usually binds to Cyclin b1, was similarly reduced in the mutants (Fig. 6A). This finding was further confirmed by in situ hybridization. In control animals, both Cyclin b1 and Cdc25b showed a specific apical-high/basal-low expression pattern in the VZ, consistent with the fact that the expression of these two factors begins at S phase, reaching their maximum levels in G2/M phase and function in G2/M phase (Fig. 6C–F). In the mutants, very weak signals were detected only at the basal germinal zone, with a loss of high expression levels at the apical surface (Fig. 6C′, D′, E′, F′). Our data suggest that Foxm1, together with its downstream factors Cyclin b1 and Cdc25b, is required for normal cell-cycle progression. Loss of Foxm1 downregulated the expression of Cyclin b1 and Cdc25b in apical half of VZ and might be why basal-to-apical INM and the length of the cell cycle were disturbed in the Foxm1 mutants.
Downregulation of G2/M regulators in the dorsal telencephalon. A, Relative mRNA levels of cell cycle regulators were determined using quantitative real time-PCR. The mRNA levels of the G2/M phase regulators Cyclin b1 and Cdc25b mRNA were reduced in the Foxm1 cko mice (Foxm1fl/fl N = 3, Foxm1 cko N = 3; Cyclin b1, p < 0.05; Cdc25b, p < 0.001), whereas other cell cycle regulators were not affected. B, The relative mRNA levels of Notch1 and its direct targets were measured using quantitative real time-PCR (Foxm1fl/fl N = 3, Foxm1 cko N = 3); a similar level of Notch1 mRNA was observed in both the wild-type and Foxm1 cko mice. C–F′, A significant decrease in the expression levels of Cyclin b1 (C–D′) and Cdc25b (E–F′) were observed at E14.5 in the dorsal (arrow) but not basal forebrain (asterisk). G–J′, Normal expression of Notch1 (G, G′), Hes5 (H, H′), Delta-like 1 (I, I′), and Ngn2 (J, J′). The data are presented as the mean ± SEM. Scale bar, 50 μm.
Notch–Delta signaling is not affected in the Foxm1 mutants
Previous studies have shown that cell fate determination in the neuroepithelium is controlled by lateral inhibition through the selective activation of Notch-Delta signaling. AP synchronization is favored by lateral inhibition, and groups of APs are synchronized at specific mitotic stages through INM (Murciano et al., 2002; Buchman and Tsai, 2008; Kageyama et al., 2009; Latasa et al., 2009; Willardsen and Link, 2011). We observed that both Notch1 and Delta-like 1 were expressed in normal apical-high and basal-low gradients in the control and Foxm1 cko mice. Moreover, no differences in the expression levels of their downstream targets Ngn2 and Hes5 were detected (Fig. 6B, G–J′). These results indicate that the increased transition from APs to BPs in the Foxm1 mutants does not involve Notch signaling and that Foxm1 acts cell-autonomously in INM and fate determination.
Disruption of Foxm1 results in decreased anxiety-related behavior in adults
Our results suggest a predominantly early embryonic role for Foxm1 in INM, the cell cycle and the AP to BP transition. To investigate whether abnormal early development affects adult cortical morphology, we performed immunostaining for NeuN, a transcription factor expressed in mature neurons. We found a slight but significant (13.3%) decrease in the number of NeuN+ neurons in mutant cortices (Fig. 7A–C); a similar defect was observed by Nissl staining (data not shown). Immunostaining for distinct cortical layer-specific markers (CDP, Satb2, Ctip2, and FoxP2) in P8 brains revealed similar decreases in the numbers of neurons in each layer (Fig. 7D–K). Furthermore, neurons generated at E12.5 (Fig. 7L,M) and E15.5 (Fig. 7N,O) were labeled by BrdU injection and observed at P8. As expected, we found a 22.6% decrease in the number of E15.5-labeled BrdU+ cells in the Foxm1 cko mice (Fig. 7P). We then tested cell death at E14.5 and P8 by immunostaining of cleaved-caspase 3 and found no significant difference between the controls and Foxm1 mutants, which excluded the possibility that the decrease in NeuN+ neurons resulted from cell death (data not shown). We next asked whether the disruption of Foxm1 had behavioral consequences. We tested adult mice for open field activity for over 15 min. Foxm1 mutants were more active than the control mice, as measured by total movement and velocity (Fig. 8A,B). The traces of the open field trajectory for each genotype showed that mutants were more likely than controls to enter the center zone (Fig. 8I). Furthermore, Foxm1 mutants tended to spend more time in the center than the controls, although this difference was not statistically significant (Fig. 8C). Because anxiety plays an important role in thigmotaxis, which drives mice to move along the wall in the open-field test, and the behavior within the first 5 min of the test best reflects anxiety levels, we measured the duration in the center zone during the first 5 min. Foxm1 mutants spent 59.0% more time in the center zone compared with the controls (Fig. 8D,E). Identical results were observed for the distance to the center zone, which reflects movement away from the center zone (Fig. 8F–H). This result raised the possibility that Foxm1 mutants have defects in anxiety behavior. The elevated O-maze test was performed to further elucidate the role of Foxm1 in anxiety-related behavior. Significant differences in anxiety levels were detected between the mutant and control mice; the frequency of entering and time spent in the open arms was significantly greater for the mutants than for the controls (Fig. 8J,K), supporting the hypothesis that Foxm1 mutants are less anxious. Because Foxm1 was also ablated in the hippocampus in our Foxm1fl/fl; Emx1-Cre model, we examined the morphogenesis of the hippocampus and found no differences between the controls and the Foxm1 mutants (data not shown). We also tested hippocampus-related contextual and cued fear conditioning behaviors. Foxm1 mutants performed similarly to the controls (Fig. 8L,M). Thus, the loss of Foxm1 has no effects on fear memory but leads to reduced anxiety-related behavior.
Decreased numbers of neurons in the Foxm1 mutant adult cortex. A, B, Coronal sections of the adult motor-sensory cortex. Immunostaining for NeuN (A, B and insets). C, The number of NeuN+ neurons in an 800 × 400 μm region of the adult sensory cortex. The Foxm1 cko mice exhibited a ∼13% decrease in NeuN+ neurons compared with the controls (Foxm1fl/fl N = 3, Foxm1 cko N = 3, p = 0.0397). D–O, Coronal sections of the P8 sensory cortex. Immunostaining for the following layer-specific transcription factors was performed: CDP (layers II/III and IV; D, E), Satb2 (layers II-V; F, G), Ctip2 (layers V–VI; H, I) and FoxP2 (layers IV–VI; J, K). The neurons generated at E12.5 and E15.5 were labeled with BrdU and observed at P8 (E12.5, L, M; E15.5, N, O). P, Number of E15.5-labeled BrdU+ neurons in a 1680 × 1260 μm area of the sensory cortex. The Foxm1 cko tissue showed a ∼22.6% decrease in the number of BrdU+ neurons compared with the controls (Foxm1fl/fl N = 3, Foxm1 cko N = 3, p = 0.0335). The data are presented as the mean ± SEM. Scale bar, 50 μm.
Foxm1 cko mice display relatively high activity levels and decreased anxiety-related behavior. A, B, The general activity levels of the Foxm1 cko mice were higher than those of the wild-type mice. The total movement was 5858 ± 308.5 cm and 7062 ± 256.5 cm for the wild-type and Foxm1 cko mice, respectively (p = 0.0077), and the average velocity was 6.5 ± 0.3 cm/s and 7.9 ± 0.3 cm/s, respectively (p = 0.0078). C–E, The Foxm1 cko mice showed a 59% increase in the duration spent in the center zone within the first 5 min of an open-field test by straight comparison of the (time spent by the Foxm1 mutants in the center zone)/(time spent by the control in the center zone) (p = 0.009). F–H, The Foxm1 cko mice traveled a shorter distance to the center zone within the first 5 min (p = 0.0042) and 15 min (p = 0.0152) of an open-field test. I, The trajectories of the Foxm1fl/fl and Foxm1 cko mice over 5 min and 15 min. Solid line outlines indicate the center zone. J, K, In the elevated O-maze test, the Foxm1 cko mice spent significantly more time in the open arms than the wild-type mice (p = 0.0057) and made more entries into the open arms (p = 0.0018). L, M, There were no significant differences in the freezing time in either the contextual or cued conditioning tests. For the behavior tests, 19 Foxm1fl/fl mice (10 males, 9 females) and 14 Foxm1 cko mice (7 males, 7 females) were tested. The data are presented as the mean ± SEM. Scale bar, 50 μm.
Discussion
In the present study, we describe new roles for Foxm1 in controlling INM and cell-cycle progression. The inactivation of Foxm1 leads to impaired INM due to a prolonged cell cycle. Impaired INM disturbs AP synchronization and further enhances the transition from APs to BPs. Our study demonstrates that Foxm1 is important for maintaining the balance between self-renewal and differentiation by controlling INM and revealed a link between INM and the cell-cycle regulator Foxm1. More importantly, Foxm1 ablation resulted in reduced anxiety-related behavior.
It has been demonstrated that the cell cycle is required for INM (Huttner and Kosodo, 2005; Schenk et al., 2009; Kosodo et al., 2011; Reiner et al., 2012). The nuclei of progenitors arrested in S phase are distributed basally in the VZ rather than moving toward the apical side, whereas nuclei arrested in G2/M phase stop migrating toward the basal side and accumulate apically (Murciano et al., 2002; Del Bene, 2011). The inhibition of G1 phase results in a basal accumulation of APs (Ueno et al., 2006). Interestingly, the converse is not true; INM is not required for cell-cycle progression. Myosin II disruption prevents apical-to-basal migration in G1 phase but has no effects on cell cycle length. In our present study, the disruption of the cell cycle regulator Foxm1 in the dorsal telencephalon delayed INM. Further analysis of cell-cycle parameters by BrdU/Ki67 and BrdU/IdU double labeling and cumulative BrdU labeling demonstrated that both S phase and the entire cell cycle were longer in the Foxm1 mutants. Moreover, 1.5/2 h, 1/5 h BrdU/IdU, and BrdU/PH3 double labeling suggested that G2 phase may also be extended in the Foxm1 mutants. Our study has established a link between INM and the cell cycle regulator Foxm1.
Foxm1 is upregulated at the onset of tumorigenesis (Raychaudhuri and Park, 2011; Koo et al., 2012), consistent with its role in controlling cell cycle progression, including mitotic entry, spindle formation, centrosome replication and the G2/M checkpoint (Laoukili et al., 2007). By controlling the expression of CDKs at the checkpoints (Kalin et al., 2011), Foxm1 regulates the G1/S and G2/M transitions. The loss of Foxm1 leads to defects in spindle segregation and, consequently, polyploidy in hepatocytes and myocardial cells (Krupczak-Hollis et al., 2004; Ramakrishna et al., 2007). However, polyploidy was not observed in neuronal progenitors in the developing telencephalon of the Foxm1 mutants. Interestingly, we found that the mRNA levels of Cyclin b1 and Cdc25b, downstream targets of Foxm1, were dramatically decreased in the Foxm1 mutants. In situ hybridization revealed an apical half-specific decrease in Cyclin B1 and Cdc25b, which suggest that the proper expression of the Cyclin b1 and Cdc25b complex may be important for INM.
A key aspect of mammalian cortical neurogenesis is cell fate determination, accompanied by the balance between self-renewal and differentiation (Bertrand et al., 2002; Miyata et al., 2004; Calegari et al., 2005; Dehay and Kennedy, 2007). INM functions to synchronize APs and to organize APs into groups, ensuring that lateral inhibition occurs in the apical region of the VZ (Ueno et al., 2006; Kageyama et al., 2009; Taverna and Huttner, 2010). In the present study, we observed increased numbers of Tbr2+ BrdU+ progenitors in the VZ, which suggested that defects in INM and the cell cycle promote the transition from APs to BPs in the VZ and enhance neurogenesis. Furthermore, the deletion of Foxm1 not only impairs INM but also disturbs AP synchronization. APs expressing neurogenic genes are dispersed, which may cause lateral inhibition throughout the whole VZ and result in increased neurogenesis (Götz and Huttner, 2005; Molnár et al., 2011). However, we did not detect any change in the expression of the proneural gene Ngn2. These data indicate that the abnormal transition from APs to BPs in Foxm1 mutants was not caused by changes in Ngn2 or Notch/Delta signaling. Interestingly, the increased number of Tbr2+ PH3+ cells in the SVZ may reflect that the BPs in the SVZ also exhibited a slower cell-cycle progression, similar to the Foxm1 deficient APs. Further studies are required to reveal the precise underlying mechanism. Together, our data suggest that Foxm1 is required to maintain the balance of AP self-renewal and differentiation. Disrupting Foxm1 disturbs INM and the synchronization of APs and promotes the transition from APs to BPs.
The Foxm1 mutants exhibited a slightly decreased number of mature neurons, and more interestingly, a relatively high level of activity, as well as reduced anxiety-related behaviors. Although Foxm1 mutants exhibit no severe cortical abnormalities, it remains unknown whether synapse number, electrophysiological characteristics, or subtle features, such as axon or dendrite positioning, are affected. The prefrontal and cingulate cortices and the amygdala play important roles in processing anxiety-related behavior (Hall et al., 2001; Viosca et al., 2009), and these three regions contain excitatory neurons originating from Emx1-positive progenitors, in which Foxm1 is normally expressed. It is possible that subtle structures in these regions are altered. Together, our data reveal a new role for Foxm1 in anxiety-related behavior. Further studies are needed to reveal the detailed mechanism.
Footnotes
This work was supported by Grant 91232301 from the National Natural Science Foundation of China (C.Z.), Grant 2012AA022402 from the National High Technology Research and Development Program of China (C.Z. and X.G.), Grant 2013CB733801 from the National Basic Research Program of China (X.G.), Grants 31171040 and 81070904 from the National Natural Science Foundation of China (C.Z.), and Grant 31200823 from the National Natural Science Foundation of China (G.C.). We thank Dr Maureen Gannon (Vanderbilt University Medical Center) and Dr Pradip Raychaudhuri (University of Illinois at Chicago) for providing the Foxm1fl/fl mice, Dr Ulrich Schüller (Ludwig Maximilian University) for the Foxm1 750 bp fragment plasmid, Yiquan Wei, Li Liu, and Yang Yang for their assistance with laboratory and animal care, and other members of the laboratory for discussion, as well as Dr Xiaojin Yu (School of Public Healthy, Southeast University) for her assistance with the mathematical analysis.
The authors declare no competing financial interests.
- Correspondence should be addressed to Chunjie Zhao, Southeast University, 87 Dingjiaoqiao Road, Nanjing, Jiangsu, 210009, China. zhaocj{at}seu.edu.cn