Abstract
Valproic acid (VPA), a widely used antiepileptic drug, is an inhibitor of histone deacetylases, which epigenetically modify cell proliferation/differentiation in developing tissues. A series of recent clinical studies in humans reported that VPA exposure in utero impaired histogenesis and the development of the central nervous system, leading to increased risks of congenital malformation and the impairment of higher brain functions in children. In the present study conducted in mice, we report that VPA exposure in utero (1) increases the amount of acetylated histone proteins, (2) alters the expression of G1-phase regulatory proteins, (3) inhibits the cell cycle exit of neural progenitor cells during the early stage of neocortical histogenesis, and (4) increases the production of projection neurons distributed in the superficial neocortical layers in embryonic brains. Together, our findings show that VPA exposure in utero alters proliferation/differentiation characteristics of neural progenitor cells and hence leads to the neocortical dysgenesis.
SIGNIFICANCE STATEMENT This study provides new insight into the mechanisms of how an altered in utero environment, such as drug exposure, affects the generation of neurons prenatally. The antiepileptic drug valproic acid (VPA) is a good target molecule as in utero exposure to VPA has been repeatedly reported to increase the risk of nervous system malformations and to impair higher brain functions in children. We show that VPA decreases the probability of differentiation of the neural progenitor cells (NPCs) in mice, resulting in an abnormally increased number of projection neurons in the superficial layers of the neocortex. Further, we suggest that histone deacetylase inhibition by VPA may be involved in the dysregulation of proliferation/differentiation characteristics of NPCs.
Introduction
Valproic acid (VPA) is a widely used antiepileptic drug, which when used during pregnancy has been reported to increase the risks of central nervous system anomalies (Jentink et al., 2010), of autism spectrum disorders (Bromley et al., 2013; Christensen et al., 2013), and of lowered intelligence quotient scores (Meador et al., 2013) in postnatal children. Thus, the International League Against Epilepsy recommends avoiding the administration of VPA to women of childbearing age (Tomson et al., 2015).
In addition to its antiepileptic activity, VPA inhibits histone deacetylases (HDACs; Göttlicher et al., 2001; Phiel et al., 2001). HDACs compress chromatin structures through the deacetylation of lysine residues in the N termini of core histone proteins, thereby decreasing the transcriptive activities of genes (Strahl and Allis, 2000; Wade, 2001). Especially in developing mammalian tissues, VPA exposure has been reported to alter the expression patterns of cell cycle-related genes in an epigenetic manner and to dysregulate the cell proliferation/differentiation characteristics of stem/progenitor cells (Tou et al., 2004).
We hypothesized that VPA exposure in utero may induce neocortical dysgenesis by altering the cell cycle kinetics of neural progenitor cells (NPCs) via its inhibitory activity on HDACs. In this study, we analyzed the proliferation/differentiation characteristics of NPCs under VPA exposure in utero, and conducted a quantitative analysis of the neocortical architecture in postnatal mice based on a mathematical model of neocortical histogenesis (Takahashi et al., 1995a, 1996b).
Materials and Methods
Animal maintenance and VPA administration.
Pregnant CD1 mice (Sankyo Labo Service) were maintained under a 12 h light/dark schedule. Since the half-life of VPA is reported to be <1 h in mice (Nau, 1985), we administered 0.4% VPA (Sigma-Aldrich) or distilled water as a control in drinking water from embryonic day 1 (E1) until birth. We limited the number of pups reared by a single mother to eight to standardize the nutritive conditions of the pups born to these prolific CD-1 dams. Embryos and postnatal mice of either sex were used. All the experimental procedures were in full compliance with the institutional guidelines of the Laboratory Animal Center, Keio University School of Medicine.
Tissue processing for histological analysis.
Postnatal mice were perfused transcardially with 4% paraformaldehyde (Sigma-Aldrich), unless specified. The analysis was conducted using serial coronal sections with 4 μm thickness, which included the dorsomedial primary somatosensory area of neocortices (field 1; Caviness, 1975).
As for the analysis in embryos, whole heads were immersion fixed with 4% paraformaldehyde overnight at 4°C. Serial coronal sections of forebrains/telencephala were cut at a 4 μm thickness, including the dorsomedial regions of the cerebral walls that correspond to the future field 1 (Takahashi et al., 1995a).
Measurement of the surface length of the telencephala.
The neocortical sections of postnatal day 4 (P4) and P21 mice were stained with 0.1% cresyl violet (Sigma-Aldrich), and were photographed by a digital camera (D700, Nikon) attached to an upright fluorescence microscope (BX51, Olympus). The surface length from the longitudinal cerebral fissure to the rhinal fissure was measured by ImageJ 1.47 software (National Institutes of Health) as an index of the surface area, which approximately represents the tangential growth of the neocortices. Eight or more brains from ≥4 litters were analyzed for the VPA-exposed mice and controls. For each brain, four nonconsecutive sections were measured.
Measurement of the number of neurons in the neocortices.
The numbers of (1) GABAergic neurons, (2) non-GABAergic neurons, and (3) glial cells in the neocortices of P21 mice were analyzed, as previously described (Goto et al., 2004; Mitsuhashi et al., 2010). Briefly, the P21 mice were perfused transcardially with 4% paraformaldehyde and 0.5% glutaralaldehyde (Kanto Chemical) in 0.1 m PBS. The sections were immunohistochemically stained using anti-GABA antibody (1:1000; A2052, Sigma-Aldrich, RRID: AB_477652; 4°C, overnight) as a primary antibody, Alexa Fluor 555 anti-rabbit IgG antibody (1:600; A-31572, Thermo Fischer Scientific, RRID: AB_10562716; room temperature, 1.5 h) as a secondary antibody, and bisbenzimide trihydrochloride (1:300 of 1% solution; H33342, Sigma-Aldrich; room temperature, 1.5 h) as a counterstaining reagent. The immunologically stained slides were photographed by digital camera, and images of the sections were examined by Zen lite 2011 software (Carl Zeiss Microscopy). Six or more brains from ≥3 litters were analyzed for the VPA-exposed mice and controls. For each brain, four nonconsecutive sections were analyzed. Additionally, GABA/GAD67 double staining and Cux1/CTIP2 double staining were conducted using anti-GAD67 antibody (1:250; MAB5406, Millipore Bioscience Research Reagents, RRID: AB_2278725; 4°C, overnight), anti-Cux1 antibody (1:50; sc-13024, Santa Cruz Biotechnology, RRID: AB_2261231; 4°C, overnight), anti-CTIP2 antibody (1:250; ab18465, Abcam, RRID: AB_2064130; 4°C, overnight), and the aforementioned anti-GABA antibody as primary antibodies, and Alexa Fluor 488 anti-mouse IgG antibody (1:600; A-11029, Thermo Fischer Scientific, RRID: AB_2534088; room temperature, 1.5 h) and Alexa Fluor 488 anti-rat IgG antibody (1:600; A-11006, Thermo Fischer Scientific, RRID: AB_2534074; room temperature, 1.5 h) in addition to the aforementioned secondary antibody in the sections of P21 neocortices fixed with 4% paraformaldehyde.
As for the analysis on P4 rats, we counted the number of neurons in sections stained only with 0.1% cresyl violet. This was because the expression level of GABA on P4 rats had been reported to be ∼50% compared with that in adult rats, and the sensitivity of immunohistochemical staining against GABA was considered to be insufficient for the analysis (Coyle and Enna, 1976). Eight brains from ≥4 litters were analyzed for the VPA-exposed mice and controls. For each brain, four nonconsecutive sections were analyzed.
Golgi staining.
Brains of P21 mice were removed and were immersed in an impregnation solution using a commercially available kit (FD Rapid GolgiStain kit, FD Neurotechnologies) based on the Golgi's silver staining technique. The brains were cut into coronal sections of 200 μm thickness and were stained according to the manufacturer's protocol. Images were captured by a microscope (BZ-9000, Keyence). Eight brains from four litters were analyzed for the VPA-exposed mice and controls.
Measurement of the size of embryos and forebrains/telencephala and of the thickness of cerebral walls.
E11 embryos were photographed using a digital camera attached to a dissecting microscope (SZX9, Olympus). The crown–rump length and the maximum rostrocaudal length of the forebrains were measured by ImageJ 1.47 software. Fifteen embryos from five litters were measured for the VPA-exposed embryos and controls.
Coronal sections, including the mediolateral cerebral walls of E10, E11, E12, E14, E16, E17, and E18 brains, were stained with 0.1% basic fuchsin (Sigma-Aldrich), and were photographed with a slide scanner (Nano Zoomer-XR, Hamamatsu Photonics). The maximum mediolateral length (width) of the forebrains/telencephala and the thickness of the cerebral walls were measured by NDP.view2 software (Hamamatsu Photonics). Eight or more forebrains/telencephala from ≥5 litters were measured for the VPA-exposed embryos and controls on each embryonic day. Four nonconsecutive sections were measured for each forebrain/telencephalon.
Estimation of the cell cycle length of NPCs in the ventricular zone.
Cumulative 5-bromo-2′deoxyuridine (BrdU) labeling analysis was conducted for the estimation of the cell cycle length of NPCs in the ventricular zone (VZ) on E10, E11, E12, E14, and E16, as previously described (Takahashi et al., 1995a). After exposure to BrdU (Sigma-Aldrich), coronal sections of the dorsomedial cerebral walls were immunohistochemically stained using anti-BrdU antibody (1:50; MCA2060, AbD Serotec, RRID: AB_323427; 4°C, overnight) as a primary antibody, and Alexa Fluor 555 anti-rat IgG antibody (1:600; A-21434, Thermo Fischer Scientific, RRID: AB_2535855; room temperature, 1.5 h) as a secondary antibody. Images of the sections were captured by a laser confocal microscope (LSM 700, Carl Zeiss Microscopy) and were examined by ZEN lite 2011 software. BrdU labeling indices (LIs; ratios of BrdU-positive nuclei to total number of nuclei) were calculated within the cerebral wall (Takahashi et al., 1995a). The LIs of the VZ were plotted against the experimental interval, and regression lines were calculated. The maximum LI corresponded to the growth fraction of the VZ. We calculated the total cell cycle length (TC) and the length of the S-phase (TS) from the regression lines of the LIs. TC − TS corresponds to the period required for the LI to reach its maximum. TC/TS corresponds to the y-intercept of the regression line. The combined length of the G2 phase and the M phase (TG2 + M) corresponds to the time required to label all mitotic cells at the ventricular surface. The length of G1 phase (TG1) was calculated as TC − (TS + TG2 + M). A total of ≥20 forebrains/telencephalons from ≥10 litters were analyzed for the VPA-exposed embryos and controls on each embryonic day. Four nonconsecutive sections were measured for each forebrain/telencephalon.
Identification of the end of neuronogenetic period in the VZ.
To identify the end of the neuronogenetic period, the disappearance of the VZ was confirmed by (1) the disappearance of the S-phase zone and (2) the attenuation of immunostaining for Pax6, a marker for progenitors in the VZ. Coronal sections of dorsomedial cerebral walls were prepared 2 h after BrdU injection on E17 and E18, and then immunohistochemically double-stained for Pax6/BrdU using an anti-Pax6 antibody (1:500; PD022, Medical and Biological Laboratories, RRID: AB_1520876; 4°C, overnight) and Alexa Fluor 647 anti-rabbit IgG antibody (1:600; A-21245, Thermo Fischer Scientific, RRID: AB_2535813; room temperature, 1.5 h) in addition to the aforementioned anti-BrdU antibody. Three brains from three litters were analyzed for the VPA-exposed embryos and controls on each embryonic day.
Analysis of the probability of cell cycle exit.
Analysis of the probability of cell cycle exit (quiescent or Q fraction) was conducted on E10, E11, E12, E14, and E16 as previously described (Takahashi et al., 1996b; Tarui et al., 2005; Mitsuhashi et al., 2010). To define a 2 h cohort of NPCs that exit the S phase, 5-iodo-2′ deoxyuridine (IdU; Sigma-Aldrich) and BrdU were administered to pregnant mice by intraperitoneal injections at 7:00 A.M. and 9:00 A.M., respectively. This time difference was chosen because BrdU is not incorporated into NPCs beyond 3.5 h following a single injection (Takahashi et al., 1992), leading to a decline in the BrdU LIs. Then, the mice were divided into two groups: (1) one in which BrdU injection was administered every 3 h (Q experiment), and (2) a second group in which BrdU was not administered again (P + Q experiment; P: proportion of daughter cells that remain proliferative, P = 1 − Q). In both groups, mice were sacrificed after more than TC − TS hour had passed after the initial BrdU injection. Coronal sections of the dorsomedial cerebral walls were immunohistochemically double stained using anti-IdU/BrdU antibody (1:25; 347580, Becton Dickinson, RRID: AB_400326; 4°C, overnight) in addition to the aforementioned anti-BrdU antibody. A total of ≥8 forebrains/telencephala from ≥3 litters were analyzed each for the VPA-exposed embryos and controls on each embryonic day. For each forebrain/telencephalon, 4–6 nonconsecutive sections were measured.
Measurement of Pax6-positive and/or Tbr2-positive nuclei.
Coronal sections of E16 dorsomedial cerebral walls were immunohistochemically double stained for Pax6 and Tbr2 using an anti-Pax6 antibody (1:100; AB528427, Developmental Studies Hybridoma Bank, RRID: AB_528427; 4°C, overnight) and an anti-Tbr2 antibody (1:200; ab23345, Abcam, RRID: AB_778267; 4°C, overnight). Images of the sections were captured by the laser confocal microscope and the numbers of Pax6-positive and/or Tbr2-positive nuclei were counted with ZEN lite 2011 software using its range indicator. Tbr2 expression was considered positive when its fluorescence level was clearly above background level; Pax6 expression was considered positive when its fluorescence level was clearly above background level and was detectable over 50% of the nucleus (Englund et al., 2005; Kowalczyk et al., 2009). Six telencephala from three litters were analyzed for the VPA-exposed embryos and controls. For each telencephalon, four nonconsecutive sections were evaluated.
Number and laminar fate of the Q cells in the postnatal neocortices.
IdU and BrdU were administered to the pregnant mice in the same manner as in the Q experiment on E16, and the coronal sections of the field 1 of neocortices on P21 were cut (Takahashi et al., 1999). In addition to the antibodies used in the Q fraction analysis, anti-calbindin antibody (1:500; AB1778, Millipore Bioscience Research Reagents, RRID: AB_2068336; 4°C, overnight) was applied to serve as a marker of layer II in the neocortices. The number of IdU-positive/BrdU-negative neurons, i.e., the “E16-born” Q cells, was counted. Four brains from three litters were examined for the VPA-exposed mice and controls. For each brain, four nonconsecutive sections were analyzed. Additionally, the E16-born Q cells were costained for Cux1 using the aforementioned anti-Cux1 antibody, to assess the association between the E16-born Q cells and the Cux1-positive superficial layer neurons.
Estimation of the cell cycle length and the population volume of secondary proliferative population.
Percent labeled mitosis method was conducted to estimate the TC of the secondary proliferative population (SPP) on E16, as previously described (Takahashi et al., 1995b; Miyama et al., 1997). A total of ≥29 telencephala from 12 litters were analyzed for the VPA-exposed embryos and controls. For each telencephalon, four nonconsecutive sections were analyzed.
Based on the fact that the SPP does not undergo interkinetic nuclear movement, separation of the VZ and the SPP population was conducted by the 1 h cohort method as previously described, provided that IdU/BrdU double-staining method was used for detection of the 1 h cohort cells (Takahashi et al., 1995b; Miyama et al., 1997). Additionally, the aforementioned anti-Tbr2 antibody was used to assess the association between the SPP cells and the Tbr2-positive basal progenitors (BPs; Englund et al., 2005). The number of proliferative cells in the VZ and the SPP was calculated by multiplication of the number of 1 h cohort cells by their estimated TCs. Nine or more telencephala from ≥4 litters were analyzed for the VPA-exposed embryos and controls. For each telencephalon, 4–6 nonconsecutive sections were analyzed.
Furthermore, the number of BrdU-positive nuclei was counted in the sections that were prepared 2 h after BrdU injection on E18. Eight telencephala from ≥3 litters were analyzed for the VPA-exposed embryos and controls. For each telencephalon, four nonconsecutive sections were analyzed.
Terminal deoxynucleotidyl transferase-mediated dUTP nick-end labeling staining.
Postnatal mice were perfused with 10% neutral buffered formalin (Nacalai Tesque) on P4. Terminal deoxynucleotidyl transferase-mediated dUTP nick-end labeling (TUNEL)-positive nuclei were detected by a commercially available kit (ApopTag Plus Peroxydase In Situ Apoptosis detection kit, Millipore Bioscience Research Reagents). Fifteen brains from five litters were analyzed for the VPA-exposed mice and controls. For each brain, four nonconsecutive sections were analyzed.
Immunoblot analysis.
The cerebral walls of E12 embryos were homogenized manually in 10 volumes of an ice-cold low-osmolality lysis buffer (0.5 m HEPES-KOH, pH 7.8, 1 m KCl, 0.5 m EDTA, 1 mm DTT, proteinase inhibitor mixture) for separation of nuclear fraction, followed by application of a high osmolality lysis buffer (0.5 m HEPES-KOH, pH 7.8, 1 m KCl, 1 m MgCl2, 68% glycerol, 1 mm DTT, proteinase inhibitor mixture). After fragmentation of genomic DNA by sonication, SDS-polyacylamidegel electrophoresis was performed using 10 μg of the nuclear protein. Precision Plus Protein WesternC standards (#161-0376, Bio-Rad Laboratories) were applied to the gels as molecular weight markers. The following antibodies were used for immunoblot analysis: anti-cyclinD1 antibody (1:250; RM-9104, Thermo Fischer Scientific, RRID: AB_720758; 4°C, overnight), anti-cyclin-dependent kinase (cdk) 2 antibody (1:250; sc-163, Santa Cruz Biotechnology, RRID: AB_631215; 4°C, overnight), anti-cdk4 antibody (1:250; sc-260, Santa Cruz Biotechnology, RRID: AB_631219; 4°C, overnight), and anti-p27Kip1 antibody (1:500; 610241, Becton Dickinson, RRID: AB_397636; 4°C, overnight). Signal intensities were standardized by comparison with those of β-actin protein. For each test, the cerebral walls of 12–15 embryos were collected from a single mother for each of the VPA-exposed embryos and controls. Three independent tests were conducted for analysis of each protein.
Measurement of the amount of total acetylated histone H3 protein.
The cerebral walls of E12 embryos were homogenized in cold PBS, and the number of cells was counted for standardization. The cells were homogenized in 10 volumes of acid lysis buffer (10 mm HEPES-KOH, pH 7.8, 1.5 mm MgCl2, 10 mm KCl, 0.5 mm DTT, 0.2 m HCl, proteinase inhibitor mixture) for extraction of histone proteins, followed by application of ice-cold acetone (Sigma-Aldrich) for protein purification. The amount of total acetylated histone H3 protein was measured by a commercially available kit based on an enzyme-linked immunosorbent assay (EpiQuik total histone H3 acetylation detection fast kit, Epigentek). For each test, the cerebral walls of 12–15 embryos were collected from a single control and VPA-exposed mother. Four independent measurements were conducted.
Statistical analysis.
Statistical significance was evaluated by two-sided Student's t tests in which p values of <0.05 were regarded as statistically significant. Additionally, analysis using a linear mixed-effects model was conducted when indicated.
Results
VPA exposure in utero increased the thickness of superficial neocortical layers on P21
The plasma concentration of VPA in the dams was 23.5 ± 5.7 μg/ml on gestation day 7 (mean ± SEM, n = 9), 18.3 ± 5.0 μg/ml on day 10 (n = 8), 37.9 ± 7.7 μg/ml on day 13 (n = 9), and 22.3 ± 5.1 μg/ml on day 16 (n = 9; Fig. 1). These concentrations were lower than the level that reportedly increased the seizure threshold by 50% in mice (120–150 μg/ml; Löscher, 1999), though marked variance was observed due to the short half-life of VPA in mice (<1 h; Nau, 1985; Fig. 1). VPA exposure did not alter the body weight of the pregnant mothers on gestation day 12 (45.7 ± 1.1 vs 47.4 ± 3.0 g in controls; p = 0.622, n = 4) or the number of mice born from a single mother (11.2 ± 0.7 vs 12.8 ± 0.6 mice in controls; p = 0.108, n = 20), compared with those in controls. The body weight of the VPA-exposed postnatal mice was decreased by 6.8%, compared with that in controls, on P17 (11.8 ± 0.2 vs 12.7 ± 0.3 g in controls; p = 0.024, n = 14). However, the body weights were not altered on P21 compared with controls (16.3 ± 0.3 vs 16.8 ± 0.2 g in controls; p = 0.111, n = 32), though greater variance was observed in the VPA-exposed mice.
We did not detect any external malformations, such as macroscopic spina bifida, in the in utero VPA-exposed mice. The dorsomedial primary somatosensory area of the neocortices (field 1; Fig. 2A; Caviness, 1975) in the VPA-exposed mice comprised six layers, similar to that in controls (Fig. 2B–D). No apparent differences were observed in the morphology of pyramidal neurons between the VPA-exposed mice and controls (Fig. 2E–G).
We measured the surface length of the telencephala as an index of the surface area, which approximately represents the tangential growth of the neocortices. The surface length of the telencephala was increased in the VPA-exposed mice by 4.6%, compared with that in controls on P21 (7.69 ± 0.10 vs 7.35 ± 0.07 mm in controls; p = 0.012, n = 8). The full thickness of the neocortices (layers I–VI) and the thickness of the superficial layers (layers II–IV) were increased by 10.0% (781.8 ± 11.4 vs 710.5 ± 13.8 μm in controls; p = 0.004, n = 6) and 14.9% (331.8 ± 6.4 vs 288.7 ± 6.5 μm in controls; p = 0.001, n = 6), respectively, in the VPA-exposed mice, compared with those in controls on P21 (Fig. 3A). In contrast, the thickness of the deeper layers (layers V and VI) was not significantly different between the two groups (383.3 ± 14.0 vs 362.5 ± 10.8 μm in controls; p = 0.290, n = 6; Fig. 3A). These results were consistent with the results of immunohistochemical staining for Cux1, a marker for the superficial layer neurons, and CTIP2, a marker for the deeper layer neurons (Fig. 2D).
VPA exposure in utero increased the number of non-GABAergic projection neurons in the superficial layers of the P21 neocortices
VPA exposure in utero increased the number of non-GABAergic projection neurons by 13.9% (431.6 ± 7.0 vs 378.9 ± 8.0 cells per counting area in controls; 250 μm in width and the total neocortical thickness in height; p < 0.001, n = 6) in field 1 of the neocortices on P21, compared with that in controls (Fig. 3B). This increase was caused solely by a 21.1% increase in the number of non-GABAergic projection neurons in the superficial layers (219.3 ± 5.9 vs 181.1 ± 7.4 cells in controls; p = 0.004, n = 6); no significant difference was observed in the deeper layers (211.1 ± 5.4 vs 196.8 ± 4.3 cells in controls; p = 0.074, n = 6; Fig. 3B). Further, analysis using a linear mixed-effects model showed a significant interaction between the increased number of non-GABAergic neurons and the superficial layers in the VPA-exposed mice (p = 0.004, n = 6). The numbers of GABAergic interneurons (28.2 ± 1.1 vs 26.4 ± 1.4 cells in controls; p = 0.389, n = 6) and glial cells (44.9 ± 3.1 vs 43.1 ± 2.0 cells in controls; p = 0.638, n = 6), or their distribution patterns were not different between the two groups (Fig. 3C–E).
VPA exposure in utero decreased the size of the forebrains/telencephala during the early development of the neocortices
When measured on E11, VPA exposure decreased (1) the crown–rump length by 7.2% (5.67 ± 0.06 vs 6.11 ± 0.09 mm in controls; p < 0.001, n = 15) and (2) the rostrocaudal length of the forebrains by 10.0% (1.31 ± 0.02 vs 1.46 ± 0.03 mm in controls; p < 0.001, n = 15), compared with those in controls (Fig. 4A). The maximum mediolateral length (width) of the forebrains/telencephala was decreased by 3.7–11.7% in the VPA-exposed embryos on E10, E11, E12, E14, and E16, compared with that in controls, whereas the width was equivalent between the two groups on E17 and E18 (Fig. 4B; Table 1).
The thickness of the cerebral walls was decreased by 6.9–14.2% in the VPA-exposed embryos, compared with those in controls, on E10, E11, E12, E14, and E16, while the thickness was equivalent between the two groups on E17 and E18 (Fig. 4C; Table 2). The layer structure of the cerebral walls was not different between the two groups (Fig. 4D).
VPA exposure in utero did not modify the length of neuronogenetic period or the total number of cell divisions of NPCs in the VZ
Newly produced neurons, identified as mitotically quiescent (Q) cells, were identified only on and after E11, but not on E10, both in the VPA-exposed embryos and in controls (Fig. 5A; Takahashi et al., 1996b). Further, the ventricular zone (VZ), defined as a homogeneous pseudostratified ventricular epithelium in which the nuclei show interkinetic nuclear movement (Bystron et al., 2008), disappeared and was displaced by the subventricular zone (SVZ) during E17 in both groups (Fig. 5B; Takahashi et al., 1995b; Englund et al., 2005), indicating the termination of the neuronogenetic period within the VZ. Namely, the neuronogenetic period in the VZ lasted for 6 d both in the VPA-exposed embryos and in controls. The TC of the NPCs in the VZ was not different between the two groups on E10, E11, E12, E14, and E16 (Fig. 5C,D; Table 3). Thus, the total number of cell divisions of the NPCs was ∼11 cycles in both groups.
VPA exposure in utero increased the maximum number of NPCs by decreasing the probability of cell cycle exit (Q fraction) during the early-middle neuronogenetic period
To investigate the effect of VPA exposure in utero on the proliferation/differentiation profiles of the NPCs in the VZ, we analyzed the probability of cell cycle exit (Q fraction) for a quantitative analysis labeling 2 h cohort of the NPCs (Takahashi et al., 1996b). The distribution pattern of the proliferative NPCs (P cells) and Q cells within the cerebral walls were almost identical between the VPA-exposed embryos and controls (Fig. 6A). However, the Q fraction, calculated as the number of Q cells divided by the number of P + Q cells, was significantly decreased in the VPA-exposed embryos on E11, E12, and E14, while it was identical between the two groups on E16 (Fig. 6B; Table 4; Takahashi et al., 1996b).
The number of NPC divisions before the Q fraction reached 0.5, the critical parameter for the maximum size of the NPC population, increased in the VPA-exposed embryos (7.6 vs 7.1 cycles in controls; Fig. 6C). This indicates an augmented expansion of the NPC population during the early and middle phase of neuronogenesis (when Q is <0.5) and increased neuronal production destined for the superficial layers during the latter phase (when Q is >0.5; Takahashi et al., 1996b, 1999). Indeed, the total number of nuclei that were positive for Pax6, a transcription factor expressed in the NPCs of the VZ, was increased on E16 in the VPA-exposed embryos by 15.0% compared with controls (130.5 ± 3.6 vs 113.5 ± 3.3 nuclei per counting area in controls; 100 μm in width and the total thickness of the VZ, SVZ, and intermediate zone (IZ) in height; p = 0.006, n = 6; Fig. 6D–F). However, the total number of nuclei that were positive for Tbr2, a transcription factor expressed in the basal progenitor cells (BPs) of the SVZ, was not different between the two groups on E16 (95.8 ± 3.1 vs 88.0 ± 3.6 nuclei per counting area in controls; 100 μm in width and the total thickness of the VZ, SVZ, and IZ in height; p = 0.135, n = 6; Fig. 6D–F).
To investigate the effects of the altered pattern of Q fraction ascent by VPA exposure in utero on neuronal production during the latter phase of neuronogenesis, we analyzed the distribution and the number of Q cells “born” on E16 in neocortical field 1 at P21. The E16-born Q cells distributed normally in layer II of the P21 neocortices (Fig. 7A,B) and were positive for Cux1 in both groups (Fig. 7C). However, the number of these E16-born Q cells was increased by 134% by VPA exposure (11.3 ± 0.9 vs 4.8 ± 0.8 cells per counting area in controls; 250 μm in width and the total thickness of neocortices in height; p = 0.003, n = 4; Fig. 7B,D).
VPA exposure in utero increased the number of neurons in the superficial layers of the neocortices on P4
We investigated neocortical architecture on P4 to confirm that the enhanced neocortical thickness in VPA-exposed mice was a result of increased neuronal production rather than overgrowth of the neuropil. Rapid growth of the neuropil takes place after P4 by glial proliferation and synapse formation, and thus the neocortical architecture of P4 directly reflects the production and distribution of projection neurons (Takahashi et al., 1996a).
VPA exposure in utero increased the total neocortical thickness by 10.4% on P4 (516.2 ± 4.3 vs 467.8 ± 7.9 μm in controls; p < 0.001, n = 8). The number of neurons in the VPA-exposed mice was increased by 20.9% in the superficial layers (future layers II–IV; 243.3 ± 3.6 vs 201.2 ± 7.7 cells per counting area in controls; 125 μm in width and the total neocortical thickness in height; p < 0.001, n = 8), while it was not altered in the deeper layers (future layers V and VI; 221.8 ± 4.0 vs 219.8 ± 6.9 cells in controls; p = 0.798, n = 8) on P4 (Fig. 8A,B). Analysis using a linear mixed-effects model showed a significant interaction between the increased number of neurons and the superficial layers in the VPA-exposed mice (p < 0.001, n = 8). There were no differences in the number of glial cells (46.3 ± 1.2 vs 46.3 ± 1.9 cells per counting area in controls; p = 1.000, n = 8; Fig. 8C), in the overall densities of neurons, or in the surface length of the telencephala (4.64 ± 0.04 vs 4.68 ± 0.05 mm in controls; p = 0.553, n = 8) between the two groups on P4.
VPA exposure in utero leads neither to the altered cell cycle kinetics of secondary proliferative population nor to the increased occurrence of apoptosis
We analyzed the cell cycle kinetics of the secondary proliferative population (SPP), since the SPP might be involved in the production of the superficial neocortical neurons (Boulder Committee, 1970; Takahashi et al., 1995b; Miyama et al., 1997; Götz and Huttner, 2005; Bystron et al., 2008; Kowalczyk et al., 2009). VPA exposure in utero did not alter the following indices of the SPP on E16, compared with those in controls: (1) the TC (13.9 vs 14.0 h in controls), (2) the distribution pattern in the cerebral wall (Fig. 8D,E), and (3) the size of population (52.4 ± 6.7 vs 59.7 ± 9.5 cells per counting area in controls; 100 μm in width and the total thickness of cerebral walls in height; p = 0.544, n = 9). Additionally, the number of BrdU-positive nuclei was not different between the two groups after a 2 h exposure to BrdU on E18 (34.4 ± 1.2 vs 31.5 ± 0.8 nuclei per counting area in controls; 100 μm in width and the total thickness of the SVZ and IZ in height; p = 0.065, n = 8), indicating that the population size of the SPP after the disappearance of the VZ was not altered by VPA exposure in utero.
The number of TUNEL-positive neurons was not different between the two groups in the neocortices on P4 (1.40 ± 0.18 vs 1.48 ± 0.16 cells per counting area in controls; 125 μm in width and the total neocortical thickness in height; p = 0.757, n = 15), at which the apoptotic activity had reached its maximum in the murine neocortices (Verney et al., 2000). Additionally, no difference in the number of pyknotic nuclei was detected in the embryonic cerebral walls between the two groups.
VPA exposure in utero increased the amount of cell cycle regulatory proteins and total acetylated histone H3 protein in NPCs
VPA exposure in utero increased the amount of cell cycle regulatory proteins cyclinD1, cyclin-dependent kinase (cdk)2, cdk4, and p27Kip1 within the nuclear fractions of the NPCs by 30%–50% on E12 (n = 3; Fig. 9A,B), and increased the amount of total acetylated histone H3 protein extracted from E12 embryonic cerebral walls by 23 ± 9% compared with controls (p = 0.044, n = 4; Fig. 9C).
Discussion
VPA exposure in utero more than doubled the number of E16-born neurons distributed in the superficial layers of the neocortices (Fig. 7D). Additionally, VPA exposure in utero increased the amounts of acetylated histone proteins and G1-phase regulatory proteins in the embryonic cerebral walls (Fig. 9A–C) and inhibited NPC differentiation during the early to middle phases of neuronogenesis (Fig. 6B,C). Our findings are compatible with previous reports finding that (1) VPA exposure in utero increases the neocortical thickness both in humans (Wood et al., 2014) and rats (Sabers et al., 2014), and (2) VPA exposure to cultured mice embryonic stem cells increased the production of superficial neocortical neurons (Juliandi et al., 2012). The administration regimen used in our study, i.e., low-dose VPA administration throughout the pregnancy, resulted in plasma concentrations lower than those needed to increase the seizure threshold in mice (Löscher, 1999); our administration regimen resembles that for pregnant women with epilepsy and better represents the pathology in human embryos than previous reports.
Mechanisms of neocortical dysgenesis by VPA exposure in utero
VPA exposure in utero modified the ascending pattern of the Q fraction (Fig. 6B,C) without altering the onset or the end of neuronogenesis or the total number of NPC divisions in the VZ. The Q alteration increased the number of NPC division cycles before Q reached 0.5 (∼E11–E14; Fig. 6C) and thus increased the peak population size of the NPCs and subsequent neuronal production, leading to neocortical thickening (Fig. 10A,B; Takahashi et al., 1996b, 1999).
NPC proliferation during the early to middle phases of neuronogenesis leads to tangential expansion of the VZ, which increases the number of ontogenic columns according to the radial unit hypothesis (Rakic, 1988; Jones and Rakic, 2010). The increased neocortical surface area in the VPA-exposed mice suggests a greater number of ontogenic columns due to augmented NPC proliferation. Contrary to a previous in vitro study, the number of deeper layer neurons was not decreased by VPA exposure in utero in the present study, suggesting that the abnormally increased NPC population was able to produce deeper layer neurons even with the inhibition of NPC differentiation, during the early to middle phases of neuronogenesis (Figs. 2B,D, 3B; Juliandi et al., 2012).
Some of the NPCs continue to divide after exiting the VZ, giving rise to projection neurons and glial cells (Boulder Committee, 1970; Takahashi et al., 1995b; Miyama et al., 1997; Götz and Huttner, 2005; Bystron et al., 2008; Kowalczyk et al., 2009). Whereas neuronal production from the SPP of the outer SVZ plays an important role in gyrus formation, particularly in primates (Lui et al., 2011; Florio and Huttner, 2014; Dehay et al., 2015), its role in lissencephalic animals is a subject of further investigation (Takahashi et al., 1995b; Mitsuhashi and Takahashi, 2016). In this context, we have shown that VPA exposure in utero did not alter the population size of SPP (determined here as the S-phase marker-positive cells in the SVZ/VZ interface or as the Tbr2-positive BPs) on E16 (Fig. 6F) or on E18 (determined as the S-phase marker-positive cells in the SVZ). Thus, we have concluded that Q alteration of the VZ NPCs is the predominant mechanism of the increase of projection neurons.
The neocortices contain not only projection neurons but interneurons produced in the ganglionic eminences as well. There was no alteration in the number of interneurons by VPA exposure in utero in the present study (Fig. 3C), although our method may not be sensitive enough to detect subtle changes in this relatively small population. However, since the great majority of neocortical neurons (≤80%) are projection neurons (Coyle and Enna, 1976), we concluded that the neocortical dysgenesis induced by VPA exposure in utero arises mainly, if not exclusively, from the increase of projection neurons produced from the VZ NPCs.
Increased expression of cell cycle regulatory proteins in NPCs by VPA exposure in utero
The G1-phase regulatory proteins control the cell proliferation/differentiation characteristics of the NPCs as evidenced by previous findings that (1) the cdk4/cyclinD1 complex prevents G1 lengthening and NPC differentiation (Lange et al., 2009), and (2) cdk inhibitor p27Kip1 inhibits cell cycle progression and induces NPC differentiation (Mitsuhashi et al., 2001; Tarui et al., 2005). A quantitative and/or qualitative balance between G1-phase regulatory proteins with opposing functions may regulate the progression of the cell cycle and the induction of differentiation. Our results showed that VPA exposure in utero induced the expression of both differentiation-inhibitory proteins (cyclinD1, cdk2, and cdk4) and a differentiation-inducing protein (p27Kip1; Fig. 9A,B). This nonspecific increase in G1-phase regulatory proteins may explain the discrepancy between our results on the SPP and those of a previous study reporting increased generation of the BPs by specific induction of cdk4/cyclinD1 (Lange et al., 2009). We speculate that the differentiation-inhibitory proteins were predominant in amount or effect compared with the differentiation-inducing proteins, during the early to middle phases of neuronogenesis when VPA exposure in utero markedly inhibited NPC differentiation (Fig. 6B,C). Additionally, previous studies have reported that (1) nonsense mutation in HDAC1 encoding genes and exposure to trichostatin A, another HDAC inhibitor, both inhibited NPC differentiation in the developing zebrafish retina (Yamaguchi et al., 2005) and (2) inactivation of histone acetyltransferase by Trrap deletion downregulates the expression of cell cycle regulatory genes and promotes NPC differentiation in mice (Tapias et al., 2014). Both results indicate regulation of the proliferation/differentiation characteristics of the NPCs by acetylation of histone proteins. Although the underlying mechanisms are still unclear, we suggest that the HDAC inhibitory activity of VPA may be involved in the inhibition of NPC differentiation.
Growth retardation during the earliest stages of embryogenesis by VPA exposure in utero
A decrease in the number of NPCs “before” the onset of neuronogenesis in the VPA-exposed embryos was indicated by the reduced rostrocaudal length of the forebrains on E11 (Fig. 4A). Furthermore, growth retardation was not restricted to the central nervous system in the VPA-exposed embryos, as we also observed decreased crown–rump length on E11 (Fig. 4A). In previous reports, VPA exposure increased caspase-3 expression in murine embryonic head tissues on E9 (Tung and Winn, 2011) and p21Cip1/WAF1 expression in a human cancer cell line (Zupkovitz et al., 2010), which might have occurred in the VPA-exposed embryos before the onset of neuronogenesis in the present study. However, VPA exposure after E12, even when these early effects of VPA were excluded, reportedly increased the neocortical thickness and neuronal number (Sabers et al., 2014).
Furthermore, as estimated from the regression curves of the Q fraction in the VPA-exposed embryos (Fig. 6C), Q alteration was capable of increasing both the maximum number of NPCs (at Q = 0.5) and neuronal production during the latter phase of neuronogenesis in the VZ, even if the number of NPCs at the onset of neuronogenesis was lower than that in controls (Fig. 10A,B). Thus, the altered pattern of the Q fraction ascent during neuronogenesis affected the neuron production sufficiently to compensate for the reduction in the NPC population size at the outset of neocortical histogenesis.
Possible mechanism for impairment of higher brain functions by VPA exposure in utero
As previously mentioned, VPA exposure in utero increases the risk of impairments in higher brain functions, including autism spectrum disorders in humans (Bromley et al., 2013; Christensen et al., 2013). Autistic features are also reported in rodents exposed to high doses of VPA over a short term during the early phase of neuronogenesis (Schneider and Przewłocki, 2005). A relative increase in excitatory to inhibitory synaptic function is suggested as a pathogenetic mechanism of autism both in humans and in mice exposed to VPA in utero (Rubenstein and Merzenich, 2003; Gogolla et al., 2009). Indeed, a previous study has reported an increased neuronal number in the postmortem brains of children with autism (Courchesne et al., 2011).
Interestingly, VPA exposure in utero increased the number of projection neurons without altering the number of interneurons in the present study, though this result could be a matter of methodological sensitivity (Fig. 3B,C). In other words, VPA exposure in utero may impair the excitatory/inhibitory ratio of neuronal function by specifically increasing the number of excitatory projection neurons (and excitatory synapses). Analyzing the higher brain functions of mice exposed to low doses of VPA throughout their embryonic period will be an important topic of future investigation.
This report showed that in utero exposure to an HDAC inhibitory agent, VPA, increased the number of excitatory projection neurons in the postnatal neocortices by inhibiting NPC differentiation in embryos. Such abnormalities may not only cause neocortical dysgenesis but also impairments in higher brain functions. We expect that further studies on the epigenetic modification mechanisms of NPC proliferation/differentiation will yield insights into the pathogenesis of “environment-induced” neuronal dysgenesis, possibly leading to novel treatments.
Footnotes
This work was supported by Grants-in-Aid for Scientific Research (B: 20390299, 23390276, and 26293248 for T.T.; C: 25461560 for T.M.) and a Grant-in-Aid for Young Scientists (B: 22791001 for T.M. and 22791038 for S. Shimozato) of the Japan Society for the Promotion of Science, and by the Mother and Child Health Foundation, the Japan Epilepsy Research Foundation, Keio Gijuku Fukuzawa Memorial Fund for the Advancement of Education and Research, Kawano Masanori Memorial Public Interest Incorporated Foundation for Promotion of Pediatrics, and the Japan Foundation for Pediatric Research Grant (for T.M.). We thank Dr. Takayuki Abe for suggestions regarding statistical analyses, Dr. Satoshi Narumi, and Ms. Namiko Saito for supporting manuscript preparation.
The authors declare no competing financial interests.
- Correspondence should be addressed to Takao Takahashi, MD, PhD, Department of Pediatrics, Keio University School of Medicine, 35 Shinanomachi, Shinjuku-ku, Tokyo, 160-8582, Japan. ttakahashi{at}keio.jp