Abstract
Abnormal patterns of head and brain growth are a replicated finding in a subset of individuals with autism spectrum disorder (ASD). It is not known whether risk factors associated with ASD and abnormal brain growth (both overgrowth and undergrowth) converge on common biological pathways and cellular mechanisms in the developing brain. Heterozygous mutations in PTEN (PTEN+/−), which encodes a negative regulator of the PI3K-Akt-mTOR pathway, are a risk factor for ASD and macrocephaly. Here we use the developing cerebral cortex of Pten+/− mice to investigate the trajectory of brain overgrowth and underlying cellular mechanisms. We find that overgrowth is detectable from birth to adulthood, is driven by hyperplasia, and coincides with excess neurons at birth and excess glia in adulthood. β-Catenin signaling is elevated in the developing Pten+/− cortex, and a heterozygous mutation in Ctnnb1 (encoding β-catenin), itself a candidate gene for ASD and microcephaly, can suppress Pten+/− cortical overgrowth. Thus, a balance of Pten and β-catenin signaling regulates normal brain growth trajectory by controlling cell number, and imbalance in this relationship can result in abnormal brain growth.
SIGNIFICANCE STATEMENT We report that Pten haploinsufficiency leads to a dynamic trajectory of brain overgrowth during development and altered scaling of neuronal and glial cell populations. β-catenin signaling is elevated in the developing cerebral cortex of Pten haploinsufficient mice, and a heterozygous mutation in β-catenin, itself a candidate gene for ASD and microcephaly, suppresses Pten+/− cortical overgrowth. This leads to the new insight that Pten and β-catenin signaling act in a common pathway to regulate normal brain growth trajectory by controlling cell number, and disruption of this pathway can result in abnormal brain growth.
Introduction
Reports that some children with autism spectrum disorder (ASD) have relatively large heads date back to the first description of autism (Kanner, 1943). In intervening years, the finding of early overgrowth of the head (macrocephaly) and brain in ASD has been reported in some, but not all, studies (for review, see Courchesne et al., 2007; Raznahan et al., 2013). This discrepancy likely reflects the fact that head and brain overgrowth occurs in a subset of the remarkably heterogeneous population of individuals with ASD, which includes individuals with microcephaly. Because of rapidly expanding knowledge of contributing genetic factors, extreme variations in head and brain size have been proposed as a candidate biomarker to help stratify the ASD population into biologically and clinically meaningful subgroups (Klein et al., 2013; Stessman et al., 2014). However, several fundamental gaps exist in our knowledge regarding abnormal brain growth and ASD, including the precise trajectory of atypical growth, the cell types and brain areas impacted, and the relationship between altered growth of the brain and the symptoms of ASD. Addressing these gaps in a clinical population presents challenges given the genetic and phenotypic heterogeneity among ASD individuals and the paucity of postmortem samples across developmental stages. Thus, a complementary approach of investigating well-defined genetic risk factors associated with ASD and extreme variations in head and brain size can help us understand whether these risk factors may influence brain development via common biological pathways and cellular mechanisms.
The gene PTEN (Phosphatase and tensin homolog) provides an entry point into this problem: PTEN encodes a phosphatidylinositol-3,4,5-trisphosphate 3-phosphatase that acts as a canonical negative regulator of the PI3K-Akt-mTOR pathway. Mutations in PTEN are reported in ∼20% of individuals with ASD and macrocephaly (Butler et al., 2005; Buxbaum et al., 2007; Varga et al., 2009; McBride et al., 2010). Extreme variations in head circumference are associated with de novo mutations in PTEN (macrocephaly) and components of the Wnt-β-catenin pathway (microcephaly and macrocephaly) in individuals with ASD (O'Roak et al., 2012; Stessman et al., 2014). Importantly, tissue growth and cell number are influenced by both Pten and β-catenin via Akt and GSK-3β signaling (Cross et al., 1995; Aberle et al., 1997; Sun et al., 1999). Furthermore, Pten and β-catenin have independently been shown to play roles in both neurogenesis and gliogenesis (Groszer et al., 2001; Chenn and Walsh, 2002; Fraser et al., 2004), and mice with elevated β-catenin signaling display brain overgrowth due to excess cell number (Chenn and Walsh, 2002). We have previously reported that adult mice with germline haploinsufficiency for Pten (Pten+/−) demonstrate both brain overgrowth and social behavioral deficits that are relevant to ASD (Page et al., 2009; Clipperton-Allen and Page, 2014).
Based on these observations, we hypothesize that Pten haploinsufficiency alters brain growth trajectory and scaling of neuronal and glial populations through elevated β-catenin signaling. Here we make use of Pten haploinsufficient mice to test this hypothesis and address three related questions: (1) What is the trajectory of brain overgrowth caused by Pten haploinsufficiency? (2) What drives brain overgrowth caused by Pten haploinsufficiency? (3) Does Pten act in a common pathway with other risk genes associated with ASD and extreme variation in head and brain size to regulate brain growth?
Materials and Methods
Mice.
All animal experiments were conducted in accordance with National Institutes of Health and Association for Assessment and Accreditation of Laboratory Animal Care guidelines and were approved by The Scripps Research Institute's Institutional Animal Care and Use Committee.
Mice of the B6.129-Ptentm1Rps line (Podsypanina et al., 1999) were obtained from the repository at the National Cancer Institute at Frederick, where they were already backcrossed onto a congenic C57BL/6J background by the donating investigator. The line has been maintained by backcrossing to C57BL/6J mice for >10 generations. Mice used in this study were generated by crossing Ptentm1Rps/+ (Pten+/−) mice with WT (Pten+/+) mice. All other mouse lines have been described previously and were obtained from the Jackson Laboratory; these consisted of Gt(ROSA)26Sortm34.1(CAG-Syp/tdTomato)Hze (Synaptophysin-tdTomato, stock #012570), Tg(TCF/Lef1-HIST1H2BB/EGFP)61Hadj/J (TCF/Lef:H2B-GFP, stock #013752) (Ferrer-Vaquer et al., 2010), Rptortm1Lex (Rptor+/−, stock #013191) (Guertin et al., 2006), Ptentm1Hwu (PtenloxP/loxP, stock #006440) (Lesche et al., 2002), Ctnnb1tm2Kem (Ctnnb1loxP/loxP, stock #004152) (Brault et al., 2001), Mtortm1.2Koz (MtorloxP/loxP, stock #011009) (Risson et al., 2009), and Emx1tm1(cre)Krj (Emx1-Cre+; stock #005628) (Gorski et al., 2002). To obtain Emx1-Cre+; PtenloxP/+; Ctnnb1loxP/+ mice, PtenloxP/+; Ctnnb1loxP/+ male mice were mated with Emx1-Cre+; PtenloxP/+; Ctnnb1loxP/+ female mice. WT, Emx1-Cre−, and Emx1-Cre+; Pten+/+; Ctnnb1+/+ mice were used as controls. To obtain Emx1-Cre+; PtenloxP/+; MtorloxP/+ mice, PtenloxP/+; MtorloxP/+ male mice were mated with Emx1-Cre+; PtenloxP/+; MtorloxP/+ female mice. WT, Emx1-Cre−, and Emx1-Cre+; Pten+/+; Mtor+/+ mice were used as controls. Mouse genotypes were confirmed by PCR using genomic DNA isolated from tail or ear samples.
Isotropic fractionator.
Isotropic fractionator was used to determine total cell number and cell type ratio in the cortex as described previously (Herculano-Houzel and Lent, 2005). Briefly, following perfusion with 4% PFA, mouse brains were fixed in 4% PFA for 1 week. Brains at P0 were then dissected into the cerebral cortex and remainder, whereas adult brains were dissected into the cerebral cortex, cerebellum, olfactory bulbs, and remainder. After dissection, the cerebral cortex was mechanically dissociated using a 7 ml glass tissue homogenizor (Kontes Glass) in standard solution (1% Triton X-100 in 40 mm sodium citrate). Nuclei in the homogenate were then stained with DAPI (Invitrogen, catalog #D3571), and the total nuclei number was counted using a hemocytometer (Fisher Scientific, catalog #0267110). For each brain, 10 μl nuclei homogenate samples was loaded into a hemocytometer, and nuclei from 0.1 μl of each sample were counted. This process was done four times for each animal, and the average was used to calculate the total cortical nuclei number of each mouse. To obtain the ratio of neuronal to non-neuronal cells, nuclei were stained with both anti-NeuN antibody (1:500, Abcam, catalog #104225) and DAPI; at least 500 nuclei stained with DAPI were counted, and the ratio of neurons was calculated by dividing the number of NeuN+ nuclei by the total nuclei number. To obtain the ratio of different glial types, microglia marker anti-Iba1 (1:500, Wako, catalog #019-19741), astrocyte marker anti-S100β (1:500, Sigma, catalog #S2532), and oligodendrocyte marker anti-Olig2 (1:1000, Millipore, catalog #AB9610) were used to stain nuclei, and the number of each glial cell type was counted. Anti-Ki67 (1:1000, Abcam, catalog #AB15580) and anti-NeuN (1:500, Abcam, catalog #AB104225) were used to determine the number of proliferating cells and neurons within the cerebral cortex at E14.5. Anti-Ctip2 (1:1000, Abcam, catalog #AB18465) and anti-Cux1 (1:500, Santa Cruz Biotechnology, catalog #sc-13024) were used to count populations of deep- and upper-layer neurons, respectively, within the cerebral cortex.
To estimate the rate of “false positives” we might encounter using this approach, we double-stained nuclei with anti-NeuN and anti-Iba1 antibodies and counted nuclei positive for both markers. We found that ∼2% of nuclei are positive for both anti-NeuN and anti-Iba1 antibodies. Because NeuN+ (neurons) and Iba1+ (microglia) cells are two distinct populations, we think this provides us with an estimate of error that occurs during the quantification process, likely due to immunoreactive cellular debris overlying nuclei.
To confirm that mice carrying only Emx1-Cre or the loxP alleles used in this study show no effects on brain mass or cell number and were thus valid to use as controls, Emx1-Cre−; loxP− (WT), Emx1-Cre−; loxP+, and Emx1-Cre+; loxP− mice were compared. No significant differences in brain mass, brain anatomical region mass, cortex cell number, or cell density were observed among these groups (data not shown); thus, data from these mice were pooled and defined as controls.
Immunohistochemistry.
Mouse brains were fixed in 4% PFA, incubated in a 20% sucrose/PBS solution at 4°C for 2 d, and embedded in Tissue-Tek OCT compound (Sakura). Coronal or sagittal sections were collected on Superfrost/Plus slides and immunostained with the following antibodies. Primary antibodies used in this study include anti-NeuN (1:2000), anti-Iba1 (1:500), anti-GFAP (1:2000, Invitrogen, catalog #130300), anti-Olig2 (1:1000), anti-Ki67 (1:5000), anti-BrdU (1:5000, Abcam, catalog #AB6326), anti-cleaved caspase 3 (1:5000, Cell Signaling Technology Technology, catalog #9661S), anti-Ctip2 (1:2000), and anti-Cux1 (1:500). AlexaFluor-488, -594, and -647-conjugated secondary antibodies (1:2000) were purchased from Invitrogen. DAPI was used for nuclear labeling (Invitrogen, catalog #D3571). Images were obtained with an Olympus VS120 microscope and processed using the VS-DESKTOP software (Olympus).
Laminar plots of anti-NeuN-immunostained sections were generated as follows: after image acquisition, ROIs were cropped from the cortex across the radial thickness using Photoshop (Adobe). Cortex thickness was normalized, and images were then auto-adjusted for each channel and saved. Plot profiles for each channel of each image were obtained using ImageJ (http://imagej.nih.gov/ij/). Three brains from each genotype, and 2 images and 2–3 ROIs (a total of 5–6 ROIs from each brain), were analyzed.
To estimate neuronal and total cell density within the histological sections of the somatosensory cerebral cortex at P0 and adulthood, 250-μm-wide columns that span the radial thickness of the somatosensory cortex were defined as ROIs. These columns were divided into 10 bins of equal size in the radial axis. For each animal, two adjacent ROIs from each of three plane-matched sagittal sections were analyzed. In each bin, neuron number (NeuN+) and total cell number (DAPI+) were counted, and density for each marker was calculated by dividing the counts by the area within each bin. Three brains from each genotype were analyzed.
To assay neuronal apoptosis in the developing cerebral cortex, sagittal sections of P4 cortex were immunostained using antibodies recognizing cleaved caspase 3 (CC3) and NeuN. The whole of the cortex, or the upper and deep layers, were defined as ROIs; and for each animal, counts of cells positive for both markers were obtained from three plane-matched sections. Three animals from each genotype were analyzed.
Analysis of cell cycle exit.
Dams carrying E14.5 embryos were injected intraperitoneally with 100 μg BrdU/g body weight and killed 24 h after injection. Immunohistochemical staining with anti-BrdU and anti-Ki67 was performed following the above-mentioned procedure. Two hundred fifty-micrometer-wide columns that span the entire radial thickness of the medial cortex were defined as ROIs. For each animal, two adjacent ROIs from each of three plane-matched sagittal sections were analyzed. The ratio of cells exiting the cell cycle was calculated by dividing the number of BrdU+, Ki67− cells by the total number of BrdU+ cells (BrdU+, Ki67−/BrdU+). Three animals per genotype were analyzed.
Western blot analysis.
P0 pups from each genotype were decapitated and cortical tissue was dissected, rapidly frozen, and stored at −80°C until used. Tissue was homogenized in RIPA buffer containing Protease Inhibitor Mixture (Roche) and phosphatase inhibitors. Total protein concentrations were measured using Pierce BCA Protein Assay Kit (ThermoScientific), and samples were diluted for equal amounts of protein in each lane. Proteins were electrophoresed onto NuPAGE 4%–12% Bis-Tris Gel (Novex, Invitrogen) and were subsequently transferred to PVDF membranes (Immobilon, Millipore). Immunoblots were blocked with 5% BSA (Sigma) in PBS (Invitrogen) for 1 h, and incubated with primary antibodies in 1% BSA at 4°C overnight. All antibodies were purchased from Cell Signaling Technology: phospho-Akt (S473; 1:1000, #4060), Akt (1:1000, #2920), phospho-S6 (S235/236; 1:2000, #4858), S6 (1:1000, #2217), phospho-GSK-3α/β (1:1000, #9331), GSK-3β (1:1000, #9315), and β-actin (1:2000, #4970). Membranes were then incubated with peroxidase-conjugated anti-mouse (Akt) or anti-rabbit (all others) IgG secondary antibodies (1:20,000, Jackson ImmunoResearch Laboratories). Proteins were visualized with enhanced chemiluminescence using WesternBright Quantum kit (Advansta) and quantified using ImageJ. For each protein of interest, a ratio of phospho-protein/total protein was calculated. Normalized values were expressed as a ratio of the mean expression level of controls. Graphs represent an average of at least two different technical replicates from at least three different animals.
Analysis of β-catenin reporter mice.
Images were acquired using identical acquisitions parameters on an Olympus VS120 microscope equipped with a 20× objective. For each genotype, two sections, matched for sectioning plane, from each of three animals were analyzed. An experimenter blind to the identity of images opened them in Photoshop and applied auto-adjustment of the image levels. The cerebral cortex was defined as an ROI, and gray value (mean) and integrated density were measured. This same ROI was then exported to ImageJ and analyzed as follows: Image>Adjust>Threshold>Auto, and then Analyze>Analyze Particles. Particle size was selected so that only individual GFP+ nuclei were counted.
Statistical analysis.
Independent-sample t tests were used to compare the Pten+/− and Pten+/+ mice, and adult Emx1-Cre+/−; PtenloxP/+ and control mice. For the P0 Emx1-Cre; PtenloxP line, one-way ANOVAs were used to assess genotype effects, and planned comparison independent-sample t tests compared heterozygous and homozygous conditional Pten mutations, as well as each genotype to controls. Genotype effects were assessed in Emx1-Cre; PtenloxP; MtorloxP and Emx1-Cre; PtenloxP; Ctnnb1loxP conditional mouse lines, as well as the Pten+/−; Rptor+/− line, by one-way ANOVA, with Tukey's post hoc tests as appropriate. All statistics were performed using Excel or PASW 18 (IBM), with significance set at p < 0.05. Statistics for nonsignificant results are not shown.
Results
Trajectory of brain growth in germline Pten+/− mice
Adult germline Pten+/− mice of both sexes show elevated brain mass without a corresponding change in body mass (Page et al., 2009; Clipperton-Allen and Page, 2014), and thus are a model of absolute and relative brain overgrowth in adulthood. To obtain a trajectory of altered brain growth associated with germline Pten haploinsufficiency, we measured brain mass across developmental stages. Increased mass was observed in male Pten+/− mice at all postnatal ages (Fig. 1A), including birth (postnatal day 0, P0; t(48) = 5.91, p < 0.001), P4 (t(30) = 4.32, p < 0.001), P7 (t(17) = 2.81, p = 0.012), P14 (t(12) = 4.62, p < 0.001), P30 (t(7) = 6.35, p < 0.001), and adulthood (P56 or older; t(23) = 15.30, p < 0.001). The magnitude of increase in brain mass relative to WT littermates varied with age, with the smallest difference between genotypes at P7 and the greatest difference between genotypes at P0 and adulthood (Fig. 1B). Body mass, however, was similar between Pten+/− mice and WT controls at all ages except birth, when Pten+/− mice weighed more than WT (t(48) = 4.64, p < 0.001; Fig. 1C). Pten+/− led to overgrowth across brain regions at P0 (cerebral cortex: t(14) = 2.68, p = 0.018; remainder: t(14) = 2.48, p = 0.026; Fig. 1D,E) and adulthood (cerebral cortex: t(12) = 9.30, p < 0.001; cerebellum: t(12) = 3.03, p = 0.010; olfactory bulb: t(12) = 2.64, p = 0.022; remainder: t(12) = 7.92, p < 0.001; Fig. 1H,I). A similar magnitude of brain overgrowth was also observed in female Pten+/− mice at P0 (brain mass: 93.73 ± 1.45 mg in WT controls, and 103.19 ± 1.32 mg in Pten+/− mice; t(14) = 4.82, p < 0.001) and adulthood (brain mass: 419.34 ± 6.35 mg in WT controls, and 497.93 ± 5.14 mg in Pten+/− mice; t(14) = 9.61, p < 0.001), suggesting that the trajectory of brain growth is grossly similar in male and female Pten+/− mice.
Among the potential underlying causes for brain overgrowth in Pten+/− mice are cellular hypertrophy (increased cell size) or hyperplasia (increased cell number). Mutations in Pten are known to influence both of these processes (for review, see Stiles et al., 2004). In the following experiments, we dissociated the relative contributions of altered cell number (hyperplasia) and density (a reflection of cell size and spacing) to Pten+/− brain overgrowth.
Cell number is increased in the cerebral cortex of newborn and adult germline Pten+/− mice
Abnormal structure and function of the cerebral cortex have been reported across studies of individuals with ASD, including delayed frontal cortex maturation (Zilbovicius et al., 1995), increased cortex size (Piven et al., 1996; Courchesne, 2002; Carper and Courchesne, 2005), increased neuronal number in the prefrontal cortex (Courchesne et al., 2011), and altered molecular networks (Boylan et al., 2007; Voineagu et al., 2011). Our previous and current results also indicate that the size and mass of Pten+/− cortex are increased, so we focused our study on the cerebral cortex. To ascertain whether cell number is changed in Pten+/− cortex, we applied the isotropic fractionator technique (Herculano-Houzel and Lent, 2005) to the cortices of WT and Pten+/− mice at P0 and adulthood. Total cell number was increased in Pten+/− cortex at both stages (P0: t(14) = 3.63, p = 0.003; adult: t(12) = 5.20, p < 0.001; Fig. 1F,J). No significant genotype difference in cortical cell density was observed (Fig. 1G,K). These results show that hyperplasia coincides with Pten haploinsufficiency-induced cortex overgrowth, whereas a change in gross cell density was not detected.
Excess neurons in the cerebral cortex of newborn, but not adult, germline Pten+/− mice
The dynamic changes in brain mass in developing Pten+/− mice compared with WT (Fig. 1B) could be due to differential allocation of neuronal and glial cell types in Pten+/− cortex, reflecting the sequential production of neurons and glia during corticogenesis. To test this, we measured neuron number by staining nuclei with anti-NeuN antibody (which labels neuronal nuclei) together with DAPI (which labels all nuclei). We found that the NeuN+/total nuclei (DAPI+) ratio was significantly higher in Pten+/− cortex at birth (t(12) = 3.24, p = 0.007; Fig. 2A), with a corresponding increase in the total number of NeuN+ cells (t(12) = 4.45, p < 0.001; Fig. 2B), while the number of NeuN negative (NeuN−) cells was similar to WT (Fig. 2C). Thus, hyperplasia in the newborn Pten+/− cortex is caused by an expanded neuronal population.
We next assessed neuronal number in adult Pten+/− cortex. Strikingly, and in contrast to our findings in the newborn cortex, at adulthood we found that the NeuN+/DAPI+ ratio was significantly lower in Pten+/− than WT cortex (t(12) = 4.21, p = 0.001; Fig. 2I). Although the total number of NeuN+ cells did not differ between genotypes (Fig. 2J), the number of NeuN− cells was significantly increased in Pten+/− cortex (t(12) = 6.44, p < 0.001; Fig. 2K).
Immunohistochemistry and laminar plot analysis of NeuN signal suggested that the gross organization of the cortex was intact in Pten+/− animals at P0 and in adulthood (Fig. 2E,F,M,N). Interestingly, overall cortex thickness was increased in Pten+/− mice at P0 (Fig. 2D), but not in adulthood (Fig. 2L). Counting cell density across the radial thickness of the cortex by bins, we found that at P0, total cell (DAPI+) density was similar between WT control and Pten+/− cortices, with the exception of a slight difference in the deepest bins (Fig. 2G). In contrast, neuronal (NeuN+) density was significantly increased in a bin that includes the upper layers of the cortical plate in Pten+/− mice (Fig. 2H). In adulthood, total cell density was similar between WT control and Pten+/− cortices, with the exception of a significant decrease in bins corresponding to upper layers in Pten+/− mice (Fig. 2O). Neuronal density was significantly decreased in bins corresponding to upper layers, as well as deep layers, in adult Pten+/− cortex (Fig. 2P). These results are consistent with the isotropic fractionator results, which indicate that the total neuron number is not changed in adulthood. In the context of an enlarged cortex with grossly normal total cell density, this implies a decrease in neuronal density.
Markers of proliferation, lamination, and apoptosis in the developing cerebral cortex of germline Pten+/− mice
One possible explanation for the increased neuron number at birth is increased neural progenitor cell proliferation during neurogenesis. To test this possibility, we quantified the number of proliferating (Ki67+) cells at embryonic day 14.5 (E14.5) using isotropic fractionator. No differences were found between WT and Pten+/− mice for brain mass (Fig. 1A,B), total cell number, cell density, and NeuN+/DAPI+ ratio (Fig. 3A–D) at E14.5. In contrast, the Ki67+/DAPI+ ratio and the total number of Ki67+ cells were significantly increased in Pten+/− cortex (Fig. 3E,F), indicating that Pten+/− mice have more neural progenitors undergoing proliferation at E14.5. In addition, when we examined exit from the cell cycle by BrdU pulse treatment at E14.5, and immunostaining for BrdU and Ki67 at E15.5, we found that the ratio of BrdU+, Ki67−/BrdU+ (ratio of progenitors exiting cell cycle) was significantly decreased in Pten+/− cortex (p < 0.001; Fig. 3G,H). Thus, Pten+/− mice have more proliferating cells and a decreased ratio of cells exiting the cell cycle.
These observations led us to hypothesize that Pten haploinsufficiency results in a relative increase in upper-layer neurons, the majority of which are born after E14.5. To test this hypothesis, we counted Cux1+ (upper-layer) and Ctip2+ (deep-layer) nuclei at P0. We found an increase in the Cux1+/DAPI+ ratio and the total number of Cux1+ nuclei in Pten+/− cortex at birth. The Ctip2+/DAPI+ ratio was not significantly changed, but the total number of Ctip2+ neurons displayed a trend toward an increase (p = 0.07) in Pten+/− cortex relative to WT (Fig. 3I–K). Because total neuron number is not significantly different between WT and Pten+/− cortices in adulthood, we predicted that the excess Cux1+ neurons are eliminated during postnatal development. Consistent with this possibility, in adult animals, we found a decrease in the Cux1+/DAPI ratio, whereas the total Cux1+ cell number was similar between WT and Pten+/− cortices (Fig. 3L–N).
Because apoptosis is a mechanism by which excess neurons generated during embryonic neurogenesis are removed during postnatal development to arrive at a mature neuron population size (Naruse and Keino, 1995), we examined apoptosis marker cleaved caspase 3 (CC3) at postnatal day 4. We found a significant increase in CC3+ neurons within Pten+/− cortex (Fig. 3O–Q), even though CC3+ neuronal density is not significantly changed. This result indicates that apoptosis is one mechanism by which excess neurons within Pten+/− cortex are removed. In addition, we also found a trend toward an increased ratio of upper-layer (II-IV)/deep-layer (V, VI, and subplate) CC3+ neurons in Pten+/− cortex (p = 0.06; Fig. 3R), suggesting that, along with a general increase in the number of apoptotic neurons throughout the Pten+/− cortex, there may also be a relative increase in upper-layer neurons undergoing apoptosis in these animals.
Excess glia in the cerebral cortex of adult germline Pten+/− mice
To determine the identity of excess non-neuronal cells in adulthood, we adapted the isotropic fractionator technique to allow for estimates of cells positive for the following glial cell type markers: Olig2 for oligodendrocytes (Yokoo et al., 2004), S100β for astrocytes (Bonaguidi et al., 2011), and Iba1 for microglia (Imai et al., 1996; Morgan et al., 2010). Importantly, all of these proteins have subcellular distributions that include the nucleus (Migheli et al., 1999; Ohsawa et al., 2000; Arnett et al., 2004). We found that both the ratio (S100β: t(6) = 2.85, p = 0.029; Olig2: t(6) = 2.51, p = 0.045; Iba1: t(6) = 3.82, p = 0.009; Fig. 4A) and the total number of nuclei positive for each of these glial markers (S100β: t(6) = 3.00, p = 0.024; Olig2: t(6) = 4.61, p = 0.004; Iba1: t(6) = 7.46, p < 0.001; Fig. 4B) were significantly increased in adult Pten+/− cortex.
To confirm the presence of excess glia, adult brain sections were immunostained with the three markers listed above. We found that the density of all three glial markers (Fig. 4C), as well as GFAP, an additional marker for astrocytes (data not shown), was increased in Pten+/− cortex. We speculated that this increase in glia might result in increased size of white matter tracts in Pten+/− brains, and we found that the thickness of the corpus callosum was significantly increased in adult Pten+/− brains (t(4) = 10.15, p < 0.001; Fig. 4D,E). This is consistent with reports of alterations in the white matter and corpus callosum of individuals with heterozygous PTEN mutations (Frazier et al., 2014).
Markers of elevated β-catenin activity in the developing cerebral cortex of newborn germline Pten+/− mice
Pten functions as a negative regulator of cell growth through the PI3K-Akt-mTOR-S6 pathway (Fig. 5A). Pten also regulates proliferation through GSK-3β, β-catenin, and cyclin D (Radu et al., 2003; Huang et al., 2007). To understand whether these pathways are dysregulated in Pten+/− cortex, we examined key readouts of these pathways. Western blot analysis of P0 cerebral cortex revealed a significant increase in the level of phospho-Akt (p-Akt; t(8) = 2.53, p < 0.05) and phospho-GSK-3β (p-GSK-3β; t(8) = 2.95, p < 0.05) in Pten+/− mice, but we did not detect a significant change in levels of phospho-S6 (p-S6) at this stage (Fig. 5B).
Because we observed elevated p-GSK-3β in newborn Pten+/− cortex, and given that GSK-3β regulates β-catenin activity (Aberle et al., 1997; Clevers and Nusse, 2012), we hypothesized that Pten haploinsufficiency leads to increased β-catenin activity in the developing brain. To test this hypothesis, we crossed β-catenin reporter TCF/Lef:H2B-GFP (Ferrer-Vaquer et al., 2010), which GFP labels the nuclei of cells that have undergone active β-catenin signaling, into a germline Pten+/− background. We then compared GFP signal between WT and Pten+/− cortices at P0 (Fig. 5C). We found that there were significant increases in GFP mean signal (t(5) = −5.11, p = 0.004; Fig. 5D) and total signal (t(5) = −5.77, p = 0.002; Fig. 5E) in Pten+/− cortex compared with WT cortex. In addition, the total number (t(5) = −3.07, p = 0.028; Fig. 5F) and density of GFP+ cells (t(5) = −2.79, p = 0.038; Fig. 5G) were also increased in Pten+/− cortex at birth. Given the lag between β-catenin activity and reporter expression (Ferrer-Vaquer et al., 2010), this likely reflects abnormally high levels of β-catenin signaling in the prenatal Pten+/− brain.
Pten conditional haploinsufficiency recapitulates cerebral cortex overgrowth
Our results thus far demonstrate that the cerebral cortex of germline Pten+/− mice captures the overgrowth seen throughout the brains of these animals and that this phenotype coincides with increased β-catenin signaling in the developing cortex. To determine whether modulation of mTOR or β-catenin activity modifies Pten+/−-induced cortical overgrowth, we wanted to isolate the effects of Pten mutations to the developing cerebral cortex. This would give greater control over cell type-specific effects of Pten mutations and would allow us to circumvent the potential problem of embryonic lethality. Homozygous germline loss-of-function mutations in Pten are lethal in early embryogenesis (Podsypanina et al., 1999); we predict that a second site genetic modifier that acts as a strong enhancer of germline Pten+/− phenotypes may also be embryonic lethal.
To test whether conditional haploinsufficiency of Pten in the cerebral cortex is sufficient to induce cortex overgrowth, we recombined one copy of a floxed allele of Pten using Emx1-Cre. This driver is expressed in cortical neural progenitor cells and descendants, including pyramidal neurons, astrocytes, and oligodendrocytes, as well as in a subset of cells in the olfactory bulb (Gorski et al., 2002). Whole-brain mass of adult male Emx1-Cre+; PtenloxP/+ mice was significantly higher than controls (t(10) = 3.37, p = 0.007), and upon dissection, cerebral cortex (t(10) = 6.15, p < 0.001) showed a significant increase in mass, whereas the mass of cerebellum, olfactory bulb, and remainder was similar to control mice (Fig. 6A). Emx1-Cre+; PtenloxP/+ female mice also exhibited similar overgrowth of the cortex (data not shown). Adult Emx1-Cre+; PtenloxP/+ cortex also showed an increase in total cell number (t(10) = 5.17, p < 0.001; Fig. 6B), no significant change in cell density (Fig. 6C,E), and a decrease in NeuN+/DAPI+ ratio (t(10) = 3.92, p = 0.003; Fig. 6D,E), confirming that Emx1-Cre+; PtenloxP/+ mice phenocopy germline Pten+/− mice and thus are a useful model to isolate the effects of Pten haploinsufficiency to the cerebral cortex.
Genetic reduction of Mtor does not suppress overgrowth in conditional Pten+/− cerebral cortex
Although our Western blot analysis of newborn Pten+/− cortex did not indicate a significant change in p-S6 levels (a readout for mTOR activity; Fig. 5B), there are likely spatial or temporal alterations in mTOR activity that may not be detected using this approach. Thus, we also wanted to test whether a genetic reduction of mTOR itself may modify Pten+/−-induced cortical overgrowth. We found that Emx1-Cre+; PtenloxP/+; MtorloxP/+ whole brain and cortex mass were not significantly different from that of Emx1-Cre+; PtenloxP/+ mice, based on one-way ANOVA and post hoc tests (Fig. 7A). We also found that neither total number nor density of cortical cells significantly differed between these two groups (Fig. 7B–D). Given the ample evidence that suppressing mTOR complex 1 (mTORC1) activity can rescue brain overgrowth in conditional homozygous Pten mutants (for review, see Zhou and Parada, 2012), we also crossed a germline loss-of-function allele for Rptor (Guertin et al., 2006), which encodes the mTORC1 effector Raptor, into a germline Pten+/− background. Consistent with our previous results, we found that genetic reduction of Raptor did not modify cortical overgrowth caused by Pten haploinsufficiency (Fig. 7E–G).
Genetic reduction of β-catenin suppresses overgrowth in conditional Pten+/− cerebral cortex
Because we found evidence for elevated β-catenin activity in Pten+/− cortex (Fig. 5), we hypothesized that genetically reducing levels of β-catenin (encoded by the gene Ctnnb1) could suppress Pten+/−-induced brain overgrowth. Thus, we generated Emx1-Cre+; PtenloxP/+; Ctnnb1loxP/+ mice and examined the cerebral cortex of these animals in adulthood. We found a significant effect of genotype on cortex mass (F(3,17) = 45.33, p < 0.001), with post hoc tests indicating that the cortex of Emx1-Cre+; PtenloxP/+; Ctnnb1loxP/+ mice weighed significantly less than that of Emx1-Cre+; PtenloxP/+ mice (p = 0.001) and more than that of control mice (p < 0.001; Fig. 8A).
To determine whether conditional haploinsufficiency of Ctnnb1 suppresses cortex hyperplasia in adulthood, we measured the total number of cells in Emx1-Cre+; PtenloxP/+; Ctnnb1loxP/+ cortex. There was a significant effect of genotype on cortical cell number (F(3,17) = 18.82, p < 0.001), with post hoc tests indicating that Emx1-Cre+; PtenloxP/+; Ctnnb1loxP/+ cortex had significantly fewer cells than Emx1-Cre+; PtenloxP/+ cortex (p = 0.018) and significantly more cells than control cortex (p = 0.013; Fig. 8B). Cortical cell density was unchanged across genotypes (Fig. 8C). Finally, we investigated whether Ctnnb1 conditional haploinsufficiency could reverse the decreased NeuN+/DAPI+ ratio seen in adult Pten+/− cortex. We found significant effects of genotype on NeuN+/DAPI+ ratio (F(3,17) = 21.83, p < 0.001) and NeuN− cells (F(3,17) = 32.07, p < 0.001), with post hoc tests indicating that the decreased NeuN+/DAPI+ ratio seen in Emx1-Cre+/−; PtenloxP/+ cortex was suppressed in Emx1-Cre+; PtenloxP/+; Ctnnb1loxP/+ cortex (p = 0.036; Fig. 8D), due to a significant decrease in NeuN− cells (p = 0.001; Fig. 8F). The number of cortical NeuN+ cells was similar among all four genotypes (Fig. 8E). Representative cortex images with stained with NeuN and DAPI are shown in Figure 8G. Our results indicate that genetic reduction of β-catenin acts as a suppressor of the cortex overgrowth, hyperplasia, and altered allocation of neuronal/non-neuronal cells caused by Pten haploinsufficiency.
Comparison of heterozygous versus homozygous conditional Pten mutations on cell number and density in the newborn cerebral cortex
Since previous reports have demonstrated that homozygous Pten mutations lead to cellular hypertrophy (Goberdhan et al., 1999; Backman et al., 2001; Groszer et al., 2001; Kwon et al., 2001; Fraser et al., 2004; Kazdoba et al., 2012), we hypothesized that heterozygous and homozygous Pten mutations may have differential effects on cell density in the developing brain. We first examined brain mass in Emx1-Cre+; PtenloxP/+ and Emx1-Cre+; PtenloxP/loxP mice at P0. We did not examine adult animals because Emx1-Cre+; PtenloxP/loxP mice were not viable after P10. One-way ANOVA revealed significant effects of genotype on brain mass (F(2,14) = 29.55, p < 0.001) and cortex mass (F(2,14) = 65.35, p < 0.001). Planned comparisons showed that brain mass of both Emx1-Cre+; PtenloxP/+ and Emx1-Cre+; PtenloxP/loxP mice was significantly increased relative to control mice (Emx1-Cre+; PtenloxP/+: t(10) = 4.03, p = 0.002; Emx1-Cre+; PtenloxP/loxP: t(10) = 6.87, p < 0.001; Fig. 9A). This increase was restricted to the cortex (Emx1-Cre+; PtenloxP/+: t(10) = 4.26, p = 0.002; Emx1-Cre+; PtenloxP/loxP: t(10) = 9.40, p < 0.001; Fig. 9A). The Emx1-Cre+; PtenloxP/loxP cortex was noticeably larger than Emx1-Cre+; PtenloxP/+ and control cortices, with Emx1-Cre+; PtenloxP/loxP having an average cortex mass almost double that of control mice. Representative sagittal brain sections show a massive increase in cortex size in Emx1-Cre+; PtenloxP/loxP mice compared with Emx1-Cre+; PtenloxP/+ and control brains (Fig. 9E).
These results indicate that cortex size does not increase linearly with decreased Pten dosage. To determine the relative contributions of altered cell number and density to the dramatic increase in cortical growth caused by homozygous Pten mutations, and to compare this to the effects of heterozygous Pten mutations, we applied isotropic fractionator to the cerebral cortex of P0 Emx1-Cre+; PtenloxP/loxP, Emx1-Cre+; PtenloxP/+, and control animals. One-way ANOVA indicated significant effects of genotype on total cell number (F(2,14) = 11.21, p = 0.001), cell density (F(2,14) = 35.68, p < 0.001), and NeuN+/DAPI+ ratio (F(2,14) = 13.54, p = 0.001; Fig. 9B–D). Planned comparisons revealed that the total number of cells in Emx1-Cre+; PtenloxP/+ cortex was significantly higher than that of control cortex (t(10) = 3.62, p = 0.005; Fig. 9B). Interestingly, although we also found increases in total cell number in Emx1-Cre+; PtenloxP/loxP cortex compared with control (t(10) = 4.69, p < 0.001; Fig. 9B), there was no significant difference in total cell number between Emx1-Cre+; PtenloxP/+ and Emx1-Cre+; PtenloxP/loxP cortices. Consistent with the large increase in mass and relatively small increase in cell number, the cell density was significantly lower in the Emx1-Cre+; PtenloxP/loxP cortex than the control (t(10) = 6.69, p < 0.001) and Emx1-Cre+; PtenloxP/+ cortices (t(8) = 8.76, p < 0.001; Fig. 9C). We also found that the NeuN+/DAPI+ ratio was significantly decreased in Emx1-Cre+; PtenloxP/loxP cortex compared with that of controls (t(10) = 2.62, p = 0.025; Fig. 9D), in contrast to the significantly increased NeuN+/DAPI+ ratio in Emx1-Cre+; PtenloxP/+ cortex (t(10) = 3.99, p = 0.003; Fig. 9D). These results suggest that Pten homozygous versus heterozygous mutations differentially impact cell density. This also gave us confidence that measurement of cell density via isotropic fractionator is capable of detecting the cellular hypertrophy that is well associated with homozygous mutations in Pten.
This nonlinear change in cortex mass, cell number, and neuronal/non-neuronal cell type allocation as a result of homozygous versus heterozygous Pten mutations led us to hypothesize that downstream molecular pathways may also be differentially impacted. We therefore examined p-Akt, p-S6, and p-GSK-3β levels via Western blot. One-way ANOVA revealed a significant effect of genotype on levels of p-Akt (F(2,12) = 21.58, p < 0.001), p-S6 (F(2,15) = 24.82, p < 0.001), and p-GSK-3β (F(2,12) = 9.69, p = 0.003). Planned comparisons showed that levels of p-Akt (t(8) = 4.76, p = 0.008), p-S6 (t(5) = 5.02, p = 0.002), and p-GSK-3β (t(5) = 3.54, p = 0.004) were significantly increased in Emx1-Cre+; PtenloxP/loxP cortex relative to control, and that p-Akt (t(8) = 2.32, p = 0.016) and p-GSK-3β (t(8) = 2.34, p = 0.046), but not p-S6, were increased in Emx1-Cre+; PtenloxP/+ cortex relative to control (Fig. 9F). Together, these data suggest that Pten gene dosage has nonlinear effects on downstream molecular pathways (Fig. 9G), with homozygous Pten mutations resulting in a robust hyperactivation of mTOR-S6 signaling in the newborn cerebral cortex.
Discussion
The brains of germline Pten+/− mice are overgrown from birth, with the magnitude of overgrowth showing a dynamic pattern through postnatal development
Heterozygous mutations in PTEN cause macrocephaly and increase risk for ASD (Butler et al., 2005; Buxbaum et al., 2007; Varga et al., 2009; McBride et al., 2010). MRI suggests that there may be a developmental trajectory to the effects of PTEN mutations on brain structure (Frazier et al., 2014). The relative rarity of PTEN mutations within the population of individuals with ASD and scarcity of postmortem brain samples, together with heterogeneity in genetic background among individuals with PTEN mutations, mean that animal models can be a valuable complement to these studies. Our use of isogenic Pten haploinsufficient mice in the present study allowed us to assess the trajectory of brain overgrowth caused by heterozygous Pten mutations and elucidate underlying cellular and molecular mechanisms, as discussed below.
Overgrowth of the Pten+/− cortex coincides with excess neurons at birth, excess glia in adulthood
Studies of individuals with ASD and macrocephaly have revealed overgrowth of the cerebral cortex in early childhood (Carper et al., 2002; Carper and Courchesne, 2005). Examination of postmortem samples from a group of children with ASD and macrocephaly indicates that hyperplasia may be a cause for brain overgrowth, with increased neuron number in prefrontal cortex and no change in neuron size or glia number (Courchesne et al., 2011). Genes regulating the cell cycle have also been reported as enriched in postmortem brains of individuals with ASD (Chow et al., 2012). MRI has revealed alterations in gray and white matter volumes in individuals with ASD (Hazlett et al., 2005; Wassink et al., 2007; Schumann et al., 2010), and those with PTEN mutations show significantly increased white matter compared with individuals with idiopathic ASD or healthy controls (Frazier et al., 2014).
Our data show that overgrowth of the Pten+/− brain coincides with cortical hyperplasia, with excess neurons at birth and excess cortical glia and thickening of the corpus callosum in adulthood. We also find that the magnitude of brain overgrowth in Pten+/− animals is dynamic over the course of development (Fig. 1B). Thus, it is worth considering the normal trajectory of neurogenesis and gliogenesis in the context of these observations. During cortical development, the timing of neurogenesis is largely prenatal, whereas gliogenesis is largely postnatal (Kriegstein and Alvarez-Buylla, 2009). The scale of the mature cortical neuronal population is influenced by a process of overproduction of neurons during embryonic neurogenesis, followed by an elimination of excess neurons by apoptosis during postnatal development (Naruse and Keino, 1995).
Pten and downstream PI3K-Akt signaling are associated with cell cycle regulation (Diao and Chen, 2007) and apoptosis (Chung and Eng, 2005), which are two cellular functions relevant for the control of cell number in the developing brain. Mutations in Pten have been reported to influence cell number in a variety of contexts, particularly in tumorigenesis (Li et al., 1997; Teng et al., 1997; Duerr et al., 1998). In the developing nervous system, conditional knock-out of Pten in neural progenitor cells leads to increased proliferation, as well as altered apoptosis (Groszer et al., 2001). Likewise, Pten conditional deletion in astrocyte precursors leads to increased proliferation (Fraser et al., 2004) and cultured Pten+/− SVZ precursor cells proliferate at a higher rate than WT cells (Li et al., 2002).
Consistent with these cellular functions of Pten, we find more proliferating cells during embryonic neurogenesis, excess neurons at birth but not adulthood, more cells undergoing apoptosis during early postnatal development, and excess glia in adulthood in the Pten+/− cortex. These observations raise the intriguing possibility that Pten haploinsufficiency leads to sequential overproduction of neurons and glia, with a correction of excess neurons (e.g., via apoptosis), but not glia, occurring during postnatal development (Fig. 10). This may help explain the relatively small difference in brain mass between Pten+/− and WT at P4 and P7 (Fig. 1B).
It remains to be determined whether the excess neurons and glia seen in mutant animals reflect a cell-autonomous role for Pten in both neurogenesis and gliogenesis, or whether excess neurogenesis and subsequent apoptosis caused by Pten mutations may trigger an overproduction of glia (e.g., to remove cellular debris via phagocytosis). Future work will also investigate how the development of neural circuitry underlying ASD-relevant behaviors is impacted by sequential overproduction and trimming of neurons during corticogenesis, and by overproduction of glia during postnatal development. Relevant to this, it has been reported that overproduction of upper-layer cortical neurons, generated by pharmacological modulation of tankyrase/Axin signaling during embryonic neurogenesis, is associated with altered excitatory inhibitory balance and ASD-relevant behavioral phenotypes (Fang et al., 2014). Elevated markers of glia, including astrocytes and microglia, and neuroinflammation have been reported in the brains of ASD patients (Vargas et al., 2005; Suzuki et al., 2013; Edmonson et al., 2014). Astrocytes and microglia are involved in brain development processes potentially relevant to ASD pathogenesis, including synaptic pruning and plasticity (Zeidán-Chuliá et al., 2014).
Pten+/− cortex shows elevated β-catenin signaling at birth, and genetic reduction of β-catenin suppresses cortical overgrowth
There is evidence that Pten regulates β-catenin activity though Akt and GSK-3β in the context of cancer biology (Stambolic, 2002). However, a connection between Pten and β-catenin has not been established in the developing nervous system. In our study, we provide the first evidence of a genetic interaction between Pten and Ctnnb1, two factors linked with ASD and extreme variations in head circumference, in regulating growth in the developing brain via control of cell number, thus placing these in a common pathway.
Our findings that β-catenin signaling is elevated in the newborn Pten+/− cortex and that genetic reduction of β-catenin can suppress Pten+/−-induced cortical overgrowth have implications for understanding the genetic landscape of autism. Extreme variations in head circumference in individuals with ASD are associated with de novo mutations in PTEN (macrocephaly), as well as in components of the Wnt-β-catenin pathway including CTNNB1 (microcephaly) and CHD8 (macrocephaly) (O'Roak et al., 2012; Stessman et al., 2014). In addition, mutations in CTNNB1 have been found to cause microcephaly and intellectual disability (de Ligt et al., 2012; Dubruc et al., 2014). Our present results inform these clinical genetic findings and indicate that Pten and β-catenin signaling act in a balance to regulate normal brain growth trajectory via control of cell number, and that disruption of this balance can lead to abnormal patterns of brain growth. Considering that CHD8 binds β-catenin and negatively regulates transcription of β-catenin target genes (Thompson et al., 2008), it would appear logical to predict that CHD8 also participates in this signaling network. It is possible that a similar mechanism might also be acting in other neurodevelopmental disorders that feature altered brain growth in subgroups of individuals, such as intellectual disability.
The focus of this study is on the trajectory and mechanisms of abnormal brain growth, with behavioral testing being beyond the scope of the present study. Future experiments will test whether ASD-relevant behaviors are modified in Pten haploinsufficient mice where cortical overgrowth and hyperplasia have been suppressed by reducing β-catenin activity. This will allow us to explore whether there is a causal relationship between hyperplasia and behavioral symptoms of ASD, and whether interventions targeting the Wnt/β-catenin signaling system might be a complement to therapies that target mTOR signaling in individuals with ASD and macrocephaly.
Homozygous versus heterozygous Pten mutations differentially impact cell density and mTOR-S6 activity in the newborn cortex
Similar to our findings in germline Pten haploinsufficient mice, we find that conditional haploinsufficiency for Pten in the developing cortex leads to an increase in total cell number but does not change cell density. Although an effect on cell size would appear likely given the function of Pten, and indeed, should be examined in future experiments, these observations indicate that hyperplasia is the primary driver of increased cortical mass caused by Pten haploinsufficiency. We find this result rather surprising given that increased neuronal soma size (hypertrophy) is well linked with brain overgrowth caused by homozygous mutations in Pten (Goberdhan et al., 1999; Backman et al., 2001; Groszer et al., 2001; Kwon et al., 2001; Fraser et al., 2004; Kazdoba et al., 2012). Consistent with these reports, we find a profound decrease in cellular density in the cerebral cortex of mice with conditional homozygous Pten mutation, indicating hypertrophy. Strikingly, there is no significant difference in total cortical cell number between conditional heterozygous and homozygous Pten mutations, despite the large increase in cortical mass between these two genotypes. We interpret this as follows: mutations that result in Pten levels or function ≥50% of WT (e.g., heterozygous loss-of-function alleles) drive changes in overall brain size primarily via hyperplasia, whereas mutations that result in Pten levels or function <50% of WT (e.g., homozygous loss-of-function alleles) drive changes in overall brain size via combined hyperplasia and hypertrophy. In agreement with this possibility, haploinsufficiency for Ctnnb1, but not Mtor or Rptor, is capable of suppressing cortical overgrowth in a Pten+/− background. Furthermore, mTORC1 activity, as indicated by p-S6 levels, is strongly potentiated in the newborn cerebral cortex of mice with conditional homozygous, but not heterozygous, Pten mutations. Our findings indicate that homozygous and heterozygous Pten mutations result in categorically different effects as regards cell density and mTOR-p-S6 signaling in the developing brain.
Footnotes
This work was supported by gift funds from Ms. Nancy Lurie Marks, startup funds from The Scripps Research Institute, and Simons Foundation Grant Award 360712 to D.T.P. Y.C. was supported by an O'Keeffe Neuroscience Scholars Award. W.-C.H. was supported by the Fraternal Order of Eagles. We thank Ms. Trina Kemp for administrative support.
The authors declare no competing financial interests.
- Correspondence should be addressed to Dr. Damon T. Page, Department of Neuroscience, The Scripps Research Institute, 130 Scripps Way, Jupiter, FL 33458. paged{at}scripps.edu