Altered temporal sequence of transcriptional regulators in the generation of human cerebellar granule cells
Abstract
Brain development is regulated by conserved transcriptional programs across species, but little is known about the divergent mechanisms that create species-specific characteristics. Among brain regions, human cerebellar histogenesis differs in complexity compared with nonhuman primates and rodents, making it important to develop methods to generate human cerebellar neurons that closely resemble those in the developing human cerebellum. We report a rapid protocol for the derivation of the human ATOH1 lineage, the precursor of excitatory cerebellar neurons, from human pluripotent stem cells (hPSCs). Upon transplantation into juvenile mice, hPSC-derived cerebellar granule cells migrated along glial fibers and integrated into the cerebellar cortex. By Translational Ribosome Affinity Purification-seq, we identified an unexpected temporal shift in the expression of RBFOX3 (NeuN) and NEUROD1, which are classically associated with differentiated neurons, in the human outer external granule layer. This molecular divergence may enable the protracted development of the human cerebellum compared to mice.
Editor's evaluation
Your paper beautifully addresses the divergent mechanisms that create species-specific features of brain development, focusing on transcriptional programs and their timing, that generate human cerebellar cells compared to those in mouse. You describe a rapid protocol for the derivation of the human ATOH1 lineage that generates excitatory cerebellar neurons from human embryonic stem cells (hESCs), and study them in vitro and in vivo. You observed transcription factors classically associated with mouse differentiated neurons expressed in the human outer external granule layer where granule cell precursors reside. These results argue that the prolonged development of the cerebellum in the human is linked to its increased size in evolution.
https://doi.org/10.7554/eLife.67074.sa0Introduction
Understanding the development of the human brain is an emerging area of neuroscience. The human cerebellum is now recognized to contribute to cognitive functions (Allen et al., 1997; Carta et al., 2019; Fiez, 1996; Schmahmann et al., 2019; Stoodley et al., 2017; Wagner et al., 2017), in addition to a critical role in motor control and motor learning. As one of the most ancient cortical regions, the cerebellum also appears central to human cognitive evolution, having rapidly expanded in absolute size, and relative to the neocortex in apes and humans (Barton and Venditti, 2014). Recent studies provide evidence of the complexity of human cerebellar histogenesis compared with nonhuman primates and rodents (Haldipur et al., 2019), making it important to develop methods to generate human cerebellar neurons to model human development and disease. The cerebellar cortex develops from rhombomere 1 (Wingate and Hatten, 1999) with a primary germinal zone that produces the Purkinje neurons and interneurons, and a secondary germinal zone, marked by the ATOH1 transcription factor, that emerges from the rhombic lip and generates cerebellar granule cells (GCs) (Hatten and Heintz, 1995). Importantly, prior to the specification of granule cell progenitors (GCPs), the ATOH1 lineage gives rise to a subset of hindbrain nuclei and cerebellar nuclei in mice (Machold and Fishell, 2005; Wang et al., 2005). While previous studies have mostly focused on generating human Purkinje cells (Buchholz et al., 2020; Erceg et al., 2010; Muguruma et al., 2015; Nayler et al., 2017; Silva et al., 2020; Wang et al., 2015; Watson et al., 2018) from human pluripotent stem cells (hPSCs), little attention has been given to defining the molecular pathways that generate human GCs and other rhombic lip derivatives. The importance of understanding the human ATOH1 lineage is underscored by the fact that GCPs are a known cell of origin for medulloblastoma, the most common metastatic childhood brain tumor (Behesti and Marino, 2009; Marino et al., 2000), and GCs are implicated in neurodevelopmental disorders including autism (Bauman, 1991; Kloth et al., 2015; Menashe et al., 2013).
Here, we report a rapid and simple protocol for the directed derivation of the human ATOH1 lineage, the precursor of excitatory cerebellar neurons, by the sequential addition of six factors, in a chemically defined culture medium. Using transgenic reporter lines and Translational Ribosome Affinity Purification (TRAP)-seq adapted to hPSCs, we tracked developmental gene expression in culture in a lineage-specific manner. Compared to previously reported studies, this method accelerates the production of ATOH1+ cells (day in vitro [DIV] 16 versus DIV35) in previous studies (Muguruma et al., 2015), with a dramatic increase in yield (80% vs. 17%). Strategies to overcome culture variability in terms of gene expression, a common limitation of hPSC-derived models, included the addition of BMP7 to stabilize ATOH1 gene expression, patterning of cells from a single-cell stage, and growth on transwell membranes where cells have access to medium in two dimensions (above and below). The translational profile of the hPSC-derived ATOH1 lineage most closely resembled human cerebellar tissue in the second trimester compared to other brain regions. Finally, we report the discovery of a shift in the expression of transcriptional regulators (RBFOX3 [NeuN] and NEUROD1) in the progenitor zone of the human external granule cell layer (EGL). NeuN and to a large extent NeuroD1 are expressed in postmitotic neurons in vertebrates (Miyata et al., 1999; Mullen et al., 1992). This molecular divergence may provide the mechanism whereby the GCP pool persists into year 2 post birth in humans, but only for 2 weeks in mice.
Results
Directed derivation of the human ATOH1 lineage from hPSCs
Addition of dual SMAD inhibitors for neuralization (Chambers et al., 2009), followed by the addition of fibroblast growth factor (FGF) and a small-molecule agonist of WNT signaling (CHIR99021, ‘CHIR’ from hereon) for posteriorization, in a chemically defined serum-free medium, induced the expression of anterior hindbrain markers EN2, MEIS2, GBX2 and repressed midbrain (OTX2), and spinal cord-level (HOXA2) markers (Figure 1—figure supplement 1A) by DIV11. FGF + CHIR treatment was superior at inducing mid/hindbrain markers (EN2, GBX2) and reducing mid/forebrain markers (OTX2/PAX6) compared to a previously reported combination of insulin + FGF2 (Muguruma et al., 2015; Figure 1—figure supplement 1A). Empirical testing of the timing, duration, and concentrations of CHIR + FGF, necessary for a ‘cerebellar territory’ expression profile of EN2+;GBX2+;OTX2-, revealed that the addition of 2.5 µM CHIR from DIV1 until at least DIV9–11 (Figure 1A, Figure 1—figure supplement 1B) was necessary and that dual SMAD inhibition until DIV7 was sufficient (Figure 1—figure supplement 1B; data not shown). While both FGF8b (Chi et al., 2003; Guo et al., 2010; Martinez et al., 1999) and FGF2 (Muguruma et al., 2015) treatments induced cerebellar territory, FGF2 induced higher EN2 levels (Figure 1—figure supplement 1B) and improved cell survival (Figure 1—figure supplement 1C).
Although this method consistently yielded cerebellar territory in two hPSC lines (Figure 1—figure supplement 1D), gene expression levels were variable between experiments. To reduce variability, a common limitation of 2D and 3D stem cell differentiation methods, we tested different culture surfaces (material and area), and the addition of signaling molecules to override stochastic gene expression. Sparse single-cell plating of hPSCs on transwell membranes allowed the formation of similarly sized and shaped colonies, resembling the growing neural plate, as patterning proceeded (Figure 1—figure supplement 1E). This improved the consistency of the geometric arrangement of the patterned cells as early colonies, reduced gene expression variability (Figure 1C), and increased cell survival, likely due to cells receiving medium in 2D (above and below) compared to 1D in regular culture dishes.
Previous work has shown that roof plate-derived BMP7 induces Atoh1 expression during GCP specification in the mouse hindbrain (Alder et al., 1999). Without knowledge of the effective concentration of secreted BMP7 in the developing human brain, we tested a range of concentrations and importantly found that low BMP7 concentrations increased ATOH1 expression, while higher levels, previously used for mouse ESC-derived GCs (Salero and Hatten, 2007), induced low levels of ATOH1 expression and caused cell death (Figure 1B; data not shown). BMP7 treatment further reduced ATOH1 expression variability (Figure 1C). Therefore, the addition of BMP7 appears to override the stochasticity of ATOH1 expression in cultures. Finally, BDNF was added to improve GC survival (Lindholm et al., 1993). Together, these experiments defined optimal conditions for deriving a human cerebellar territory and ATOH1 lineage specification as first steps towards derivation of the ATOH1 lineage and GC differentiation (Figure 1D).
To monitor the dynamics of ATOH1 lineage specification and differentiation in culture, we derived a clonal ATOH1-EGFP hPSC line, expressing nuclear EGFP under a human ATOH1 enhancer (Figure 1—figure supplement 2A). Fluorescence-activated cell (FAC)-sorting of ATOH1-EGFP+ cells at DIV16 followed by RT-qPCR revealed coexpression of genes associated with the cerebellar territory in mice (Allan Brain Atlas, Morales and Hatten, 2006) including ATOH1, PAX6, EN2, ID4, LHX2, LHX9, and MEIS2 (Figure 1—figure supplement 1F) in ATOH1-EGFP+ cells. Time-series analysis of ATOH1-EGFP by imaging and flow cytometry revealed that by DIV16, 80 ± 9% of the cells are ATOH1-EGFP+ (Figure 1E). This is a fivefold increase in the efficiency of ATOH1+ cell derivation compared to a previously reported 3D protocol (80 vs. 17%; Muguruma et al., 2015). By DIV23, the percentage of ATOH1-EGFP+ cells decreased, while ATOH1-EGFP- cells increased, indicating the onset of differentiation (Figure 1E). To investigate when the ATOH1-EGFP+ cells coexpress known GCP markers, gene expression was analyzed in FAC-sorted ATOH1-EGFP cells at three time points between DIV16-23. In the mouse, the selected genes are dynamically expressed within the Atoh1+ domain, changing between embryonic day (E) 11.5, when the Atoh1+ population gives rise to hindbrain and cerebellar nuclei, and E15.5, when the Atoh1+ population forms the EGL (Figure 1—figure supplement 1F, right; Machold and Fishell, 2005; Wang et al., 2005; Allan Brain Atlas). Consistent with the reduction in ATOH1-EGFP+ cells by DIV23, gene expression shifted between DIV19 and 23, with markers expressed in the mouse Atoh1 domain at E11.5, prior to EGL establishment (ID4, LHX2, LHX9), decreasing in expression, while EGL markers increased (PAX6, ATOH1) (Figure 1—figure supplement 1F). Thus, at DIV23 the hPSC-ATOH1 lineage initiates GCP production but likely contains a mixture of progenitors. Indeed, prior to DIV23, at DIV16, a small subset of cells expressed Calretinin by immunohistochemistry (Figure 1—figure supplement 3A), a marker of excitatory cerebellar nuclei produced by the Atoh1 lineage prior to GC production. By immunohistochemistry, at DIV28, the ATOH1+ cells coexpressed PAX6 (mean ~70%) and EN2 (mean ~50%), and ~40% coexpressed all three markers within the EGFP+ population. The percentage of coexpression of these three markers within the EN2+ (marker of mid/hindbrain) population was ~70% (Figure 1F). Thus, a great majority of the EN2+ cells are GCPs. Midbrain progenitors that express EN2 do not express ATOH1/PAX6 (Akazawa et al., 1995; Allan Brain Atlas), and dorsal posterior hindbrain/spinal cord Atoh1+ never express EN2 (Davidson et al., 1988).
To examine the identities of cells that are negative for ATOH1-EGP, we performed double immunolabeling with several other markers at DIV28. These analyses revealed the presence of ATOH1-EGFP-;PAX6+,NeuN+ cells (Figure 1—figure supplement 3B, asterisk), indicative of differentiated GCs. A small number of Calretinin+ cells, indicative of cerebellar nuclei or unipolar brush cells (Abbott and Jacobowitz, 1995; Fujita et al., 2020; Wizeman et al., 2019), and ATOH1-EGFP-;SOX2+ cells (28 ± 16% [mean ± 1 SD of the SOX2+ cells, n = 2]), indicative of cells deriving from the ventricular zone neuroepithelium in the developing cerebellum (GABAergic lineage, Figure 1—figure supplement 3C). Interestingly, a great majority of the SOX2+ cells were ATOH1-EGFP+ (72 ± 16.7% [mean ±1 SD, n = 2]; Figure 1—figure supplement 3C and D). In the developing mouse cerebellum, SOX2 is abundantly expressed throughout the ventricular zone neuroepithelium, but is a rare find in the EGL (Sutter et al., 2010).
A defining feature of GCPs is their extensive proliferative capacity in response to sonic hedgehog (SHH) in the postnatal mouse cerebellum (Dahmane and Ruiz i Altaba, 1999; Lewis et al., 2004; Wallace, 1999; Wechsler-Reya and Scott, 1999). EdU uptake, a measure of cell proliferation, was significantly increased in the ATOH1-EGFP+ cells (Figure 1G) after treatment with SAG (an agonist of SHH signaling) compared to control. In conclusion, the hPSC-ATOH1+ lineage displayed similar dynamics in gene expression over time as seen in the mouse embryo (albeit extended in period) and produced GCPs, which respond to SAG, by DIV23–28. Compared to previously reported methods, both the speed (35 days versus 16 days) and the yield (from 17% to 80%) of GCP production were increased.
Differentiation of hPSCs into cerebellar granule cells
In the mouse, early postmitotic GCs switch off Atoh1 and transiently express TAG1 in the inner EGL. In our cultures, TAG1 was expressed already at DIV18, after which its expression increased (Figure 2A). Magnetic-activated cell sorting (MACS) using an antibody against TAG1, which is expressed on the cell surface, yielded 13 ± 7% (n = 7) TAG1+ cells at DIV28 (Figure 2B and C). TAG1+ cells were then co-cultured with either mixed cerebellar mouse neurons or glia to provide a permissive differentiation environment. By DIV3 in co-culture, the majority of the TAG1+ cells were PAX6+ (Figure 2D, top panel). By contrast, the TAG1- fraction contained mostly cells with larger nuclei that were PAX6- (Figure 2D, bottom panel). By DIV20, co-culture (DIV48 total days) of TAG1+ cells with mouse neurons or glia resulted in a majority of the human cells displaying characteristic GC morphology, namely a small (<10 um) round nucleus and a limited number of extended processes (3–4) including bifurcated axons (Figure 2E, 136/230 examined, 60%). Mouse GCs isolated at P0 and grown concurrently in the same dish as the human cells displayed similar features (Figure 2E). The human cells expressed the GC markers NEUROD1 (42/56 examined, 75%, Figure 2E) and PAX6 (174/218 examined, 80%). In double-labeling experiments, all examined cells that expressed PAX6 (18/22, 81%) also expressed NEUROD1. Interestingly, presynaptic specializations were apparent in neurons grown in co-culture with mouse glia only, indicating that cross-species neuron-glia interactions are sufficient to trigger synapse initiation (Figure 2F). In addition, the expression of vesicular glutamate receptor 1 (VGLUT1), a presynaptic marker that in the cerebellum is localized to GC parallel fiber-Purkinje cell excitatory synapses, was detected (Figure 2G).
While we carried out the most comprehensive characterization of the MACsorted cells at DIV28 and their differentiation in co-culture with mouse cells, beyond DIV28, it should be noted that the human cultures continued to express ATOH1-EGFP. TAG1 came on transiently in subsets of cells also at later stages (~13% at DIV35; data not shown). At DIV28 + 20, the TAG1- fraction contained cells with several morphologies (Figure 2—figure supplement 1), including cells with small round nuclei that expressed NEUROD1 (~47%) and NeuN (~70%). A small fraction was Calretinin+ (~5%). The majority of the Calretinin+ cells at this stage (70 ± 36% of the 5%) were double positive for NeuN.
To assess whether hPSC-GCs would integrate into the mouse cerebellar cortex, especially whether they would undergo the classic glial-guided migration of GCs, we implanted human cells into the juvenile mouse cerebellum at the time when mouse GCs are undergoing migration. Importantly, TAG1-sorted cells integrated into the neonatal mouse cerebellar cortex and migrated along glial fibers with stereotypical morphology of migrating neurons passaging through the EGL to settle in the internal granule cell layer (IGL) (Figure 3A, Figure 3—figure supplement 1). Moreover, we observed the stereotypical elongated morphology of migrating neurons on glia when TAG1-sorted cells were co-cultured with glia isolated from the mouse cerebellum (Figure 3B) for 36 hr, similar to previous observations with mouse GCs (Hatten, 1985). By contrast, cells that were not attached to glia displayed more rounded morphologies (Figure 3B). Together, these experiments demonstrate that the DIV28 hPSC-GCs can undergo glial-guided migration both in vitro and upon transplantation in vivo.
hPSC-ATOH1+ cells match the molecular profile of the human cerebellum in the second trimester
To transcriptionally profile the hPSC-ATOH1+ progenitors, we adapted the TRAP methodology, first developed in transgenic mice for cell-type-specific translational profiling (Doyle et al., 2008; Heiman et al., 2008), for hPSCs. An EGFP-tagged L10a ribosomal subunit was driven by the human ATOH1 enhancer (Figure 1—figure supplement 2B), enabling GFP-mediated immunoprecipitation (IP) of ribosomally attached mRNAs in the ATOH1 lineage specifically. Importantly, this method bypasses the need for cell dissociation and provides information about transcripts, including the 3′UTR, that are present in both the cell soma and processes. TRAP followed by RNA sequencing (TRAP-seq) was performed at DIV28 when the ATOH1 lineage coexpressed GCP markers and displayed increased proliferation in response to SAG (Figure 1). The ATOH1 transcript was indeed enriched in IPs compared to the input (Figure 1—figure supplement 2B, Figure 4—figure supplement 1A). Several other GCP markers were similarly enriched, while markers of differentiated GCs were depleted (Figure 4A). Heatmaps of the major signaling pathways in development (WNT, BMP, SHH, FGF, HIPPO) highlighted enrichment of the WNT pathway in particular (Figure 4—figure supplement 1B). Interestingly, gene ontology analysis revealed axon guidance, the WNT pathway, neuronal migration, and cell division among the top significantly enriched developmental processes (Supplementary file 1).
To investigate how the hPSC-ATOH1+ cells compare to the molecular signature of the developing human cerebellum, we compared our dataset to RNA-seq data from the PsychEncode study (Li et al., 2018), which samples different brain regions during a developmental timeline covering 8 post coitus week (PCW) until after birth in the human. Gene Set Enrichment Analysis (GSEA) revealed that the DIV28 ATOH1 lineage most closely matched the profile of the 13–17 PCW human cerebellum (Figure 4B and C, Figure 4—source data 2). Moreover, comparisons to published single-cell RNA-seq data from the developing mouse cerebellum revealed most resemblance to various cerebellar glutamatergic lineages (Figure 4—figure supplement 2, Figure 4—figure supplement 2—source data 1). There was no resemblance to glia or endothelial cells. Together, these analyses provide the first translational dataset for the hPSC-derived ATOH1 lineage. While the data closely match the developing human cerebellum, development appears accelerated in culture.
Molecular profiling reveals a transcriptional shift in the human versus mouse cerebellum
In a surprising finding, TRAP-seq analysis of the ATOH1 lineage revealed the expression of several genes that are classically associated with postmitotic GCs in mice including RBFOX3 (encodes for the NeuN antigen) and NEUROD1 (Figure 5A). The coexpression of NEUROD1 and ATOH1 was confirmed by RT-PCR in IPs (Figure 5A) and in FAC-sorted EGFP+ cells derived from the ATOH1-EGFP line (data not shown). By immunohistochemistry, ATOH1-EGFP+ cells expressed NeuN and NEUROD1 (Figure 1—figure supplement 3B). To investigate if the expression of these factors extends into the proliferative zone of the human EGL in vivo or if this is an in vitro phenomenon, we performed immunohistochemistry in the human cerebellum at 17 PCW as our data most closely matched the human 13–17 PCW. We performed concurrent immunohistochemical analyses on P0 mouse cerebella because the 17 PCW human cerebellum most closely resembles P0 in mice, based on cerebellar foliation depth and pattern (Biran et al., 2012; Haldipur et al., 2019). In contrast to the mouse, many outer EGL cells in the human, where the ATOH1+ progenitors reside (Haldipur et al., 2019), expressed both NEUROD1 (Figure 5B and D, Figure 5—figure supplement 1A) and NeuN (Figure 5C and D). Moreover, Ki67, a marker of proliferating cells, was not as prevalent at this stage in the human EGL as it is in the P0 mouse cerebellum (Figure 5B, high magnification, and C). In the human EGL, Ki67 expression was sparser, Ki67 and NEUROD1 coexpression was clearly evident, and there was no clear separation of cells into two layers as seen in the mouse. Extensive examination of multiple stages of development in the mouse (from E15.5 until P6) and multiple regions of the developing cerebellum revealed very few NEUROD1+ cells at the pial surface of the EGL (Figure 5—figure supplement 1B, C). Only at P6, a relatively late stage of cerebellar development in the mouse, when GCP proliferation peaks (Ki67 throughout the EGL) and the EGL is at its thickest, did we detect an increase in NEUROD1+ cells in the oEGL (Figure 5—figure supplement 1B and C). Punctate NeuN expression was detected in the human oEGL but not at any stages or regions examined in the mouse (Figure 5D, Figure 5—figure supplement 1C). Finally, since a large number of ATOH1-EGFP+ cells were SOX2+ in our cultures (Figure 1—figure supplement 3C and D) we examined SOX2 expression in the 17 PCW human cerebellum and found that SOX2 is extensively expressed in the human EGL (Figure 5—figure supplement 1D), in contrast to the mouse, where SOX2+ cells in the EGL are a rare find (Sutter et al., 2010). Together, these data reveal marked species differences and extensive coexpression of transcriptional regulators that have classically until now been ascribed to either progenitors or differentiating neurons based on mouse studies.
Discussion
We report a method for the scalable derivation of the human ATOH1 neuronal lineage from hPSCs that yields cerebellar progenitors by day 16 and GCs within 48 days in chemically defined medium. In contrast to previous methods for derivation of cerebellar neurons with a focus on Purkinje cells (Buchholz et al., 2020; Erceg et al., 2010; Muguruma et al., 2015; Nayler et al., 2017; Silva et al., 2020; Wang et al., 2015; Watson et al., 2018), we provide a timeline of developmental progression of the human ATOH1 lineage by gene expression and an unbiased comparison of the transcriptional profile to the developing human brain. While our focus was derivation of GCs, the expression of markers of the cerebellar nuclei (LHX9, LHX2) and the presence of Calretinin+ cells indicate the production of additional ATOH1 derivatives. Future identification of cell surface markers for each class of glutamatergic neuron known to derive from the ATOH1 lineage will allow isolation and further characterization of the various populations yielded by our method. MAC-sorted early postmitotic GCs (TAG1+) integrated into the mouse cerebellum and migrated through the EGL into the IGL in close apposition to glial fibers, providing an in vivo system to study human GC migration and migration defects. This is the first successful demonstration of migration followed by integration of hPSC-GCs into a stratified cortex, suggesting that molecular pathways required for glial-guided migration have developed in the human cells.
Adaptation of TRAP-seq for use in hPSCs allowed lineage-specific molecular profiling of the ATOH1 lineage, providing a superior depth of read compared to scRNA-seq methods. Comparison of our data to RNA-seq data from the developing human brain matched the DIV28 hPSC-ATOH1 lineage to the human cerebellum at 13–17 PCW. Analysis of key developmental pathways revealed that the WNT pathway is particularly enriched. Indeed, mouse genetic studies show that WNT signaling is critical for both early cerebellar development and later circuit establishment (Lucas and Salinas, 1997; McMahon and Bradley, 1990).
Interestingly, using TRAP-seq methodology we identified a molecular divergence in the temporal expression of transcriptional regulators in progenitors in the human EGL, which in the mouse are expressed in cells that are further along the differentiation path. Genetic differences in developmental timing can cause birth defects or give rise to a novel morphology that could confer an evolutionary advantage. The divergent expression of NEUROD1 and RBFOX3 in the human oEGL suggests an expansion of an ‘intermediate’ cellular state that may serve to allow the progenitor pool to persist much longer, to enable the protracted period of cerebellar development in humans compared to mouse. Both NEUROD1 and RBFOX3 are expressed in early postmitotic as well as fully differentiated neurons, rarely in neuroblasts, across species and brain regions (D’Amico et al., 2013; Lee et al., 1995; Miyata et al., 1999; Zhang et al., 2016). Overexpression of either gene in progenitors induces neuronal differentiation (Boutin et al., 2010; Butts et al., 2014; Lee et al., 1995; Pataskar et al., 2016; Zhang et al., 2016). Interestingly, in the Xenopus embryo, the coexpression of Atoh1 and NeuroD1 in the EGL has been reported, and it was hypothesized that both the timing of expression and the gene regulatory function of NEUROD1 have adapted to support development of a nonproliferative EGL (Butts et al., 2014). In the human, we propose that the coexpression of these factors in GCPs, together with the more extensive expression of SOX2 in the human EGL, may enable an extended immature (stem cell) or quiescent state to preserve the progenitor pool. We speculate that different levels of each factor, in combination, may be important for proliferation versus quiescent states in the GCPs. Indeed, Ki67 labeling showed fewer positive cells in the human EGL at 17 PCW than the comparable P0 mouse EGL, hinting at differences in proliferative regulation. Our findings, together with findings in the Xenopus and the differential proliferative capacity of the EGL in several vertebrate species (Iulianella et al., 2019), support the existence of evolutionary adaptations in EGL development across species. In future work, analyses of cell cycle length and comparison of molecular features in human versus mouse GCPs should provide further insights into the molecular basis of human cerebellar expansion and human-specific molecular mechanisms that will be critical for understanding medulloblastoma pathogenesis as well as cerebellar-mediated human cognitive evolution.
Materials and methods
Human tissue collection, fixation, and embedding
Request a detailed protocolHuman prenatal brain tissue (two cerebella; 17 PCW) were acquired from the Human Developmental Biology Resource (http://www.hdbr.org/) following institutional policies. Tissues were fixed in 4% PFA for 7–10 days and washed multiple times in PBS. Samples were then cryoprotected in increasing concentrations of 5, 15, and 30% sucrose (in PBS) at 4°C. Samples were embedded in Tissue-Tek O.C.T. compound (VWR, 25608-930) and stored at –80°C until use.
Mice
All procedures with mice were performed according to guidelines approved by the Rockefeller University Institutional Animal Care and Use Committee. C57Bl/6J mice (Jackson Laboratory) were maintained on a 12 hr light/dark cycle and a regular diet. Timed matings generated a mixture of male and female pups that were used for all described studies.
hPSC culture
Request a detailed protocolHuman embryonic/pluripotent stem cells were used under approved institutional ESCRO committee protocols (The Rockefeller University). The RUES2 line (NIH #0013) was created, authenticated, and provided by the lab of Dr. Ali Brivanlou at the Rockefeller University. The line was further karyotyped to confirm the sex (female) and karyotype (normal). The W9 line (NIH #0062) was obtained from Dr. Brivanlou’s lab where it was routinely used in stem cell differentiation experiments along the three main embryonic lineages. hPSCs were maintained in growth media (HUESM medium conditioned with mouse embryonic fibroblasts and supplemented with 20 ng/ml bFGF [Invitrogen]) (Deglincerti et al., 2016). The growth medium was exchanged daily. For transgenic lines, Puromycin (Gibco) was added to the medium (1 μg/ml) during maintenance culture. Cells were grown as colonies on tissue culture dishes coated with hESC qualified Matrigel solution (Corning) in a 37°C humidified incubator with 5% CO2.
hPSC differentiation
Request a detailed protocolhPSCs maintained as colonies were dissociated from plates with Trypsin-EDTA (0.25%, Gibco) for 4 min at 37°C in a humidified incubator. Cells were washed once with growth media (see previous section) and then resuspended in growth media with 10 μM ROCK-inhibitor Y-27632 (Abcam). Single cells were plated at 900 cells/ml on Transwell 6-well plates with permeable 24 mm polyester membrane inserts (Corning, 3450, matrigel-coated) or regular tissue culture-treated plates (matrigel-coated). On Transwell dishes, cells were plated on top of the membrane with 1 ml growth medium plus ROCK-inhibitor added below and above the membrane (2 ml total). The next day, 1 ml growth medium plus ROCK-inhibitor was added to the top part of the Transwell (3 ml total). On day 2, the medium was switched to differentiation medium: DMEM/F12 with sodium bicarbonate (Invitrogen, 11320-033), 0.5% BSA, 0.1 mM β-mercaptoethanol, 2 mM glutamate, 10 μM NEAA, 1× N2 supplement, 1× B27 without retinoic acid (all from Gibco) and Primocin (0.1 mg/ml, InvivoGen). This marked day 0 of differentiation. For differentiation experiments with transgenic lines, Puromycin (Gibco) was added to the medium (0.5 μg/ml) until DIV28. The following small molecules and growth factors were added to the differentiation medium on days indicated in Figure 1D as follows: 10 μM SB431542 (SB, Tocris) and 100 nM LDN-193189 (Stemgent) at days 0–7, 2.5 μM CHIR99021 (Stemgent) at days 1–11, 20 ng/ml bFGF (Invitrogen) at days 1–28, 25 ng/ml BDNF at days 11–28, and 50 ng/ml recombinant mouse BMP7 (R&D Systems) at days 7–15. Medium was exchanged for fresh differentiation medium plus appropriate factors every other day until day 28. The following conditions were also tested (not all data are included in the article) to arrive to the optimized protocol for the derivation of the ATOH1 lineage: CHIR99021 (Stemgent) at days 1–11 (range tested: 1–3 μM), 7 μg/ml human recombinant Insulin (Tocris), 100 ng/ml recombinant human/mouse FGF8b (Tocris) between days 1–11 (range tested: 50–500 ng/ml), recombinant mouse BMP7 (R&D Systems) at days 7–15 (range tested: 20–1000 ng/ml), 20 ng/ml BMP6 (R&D Systems), BMP4 (range tested: 4–160 ng/ml, R&D Systems), and 100 ng/ml GDF7 (R&D Systems). Depending on culture conditions used, such as culture surface, the concentration of CHIR99021 may need to be adjusted to achieve optimal anterior hindbrain patterning.
Generation of transgenic hESC lines
Request a detailed protocolFor derivation of the ATOH1-EGFP line, a human ATOH1 enhancer sequence (GenBank accession number AF218259.1; Helms et al., 2000) was amplified from human genomic DNA using two in-fusion primers (forward primer: 5′-TTCAAAATTTTATCGATaaggttcttCTATGGAGTTTGCA-3′; reverse primer: 5′- AATAGGGCCCTCTAGAGAATTCCTGAACAACCCCAC-3′). The amplicon was cloned into a modified version of the self-inactivating lentiviral vector pPS-EF1α-GFP-RFP (System Biosciences, LV603PA-1) where GFP and RFP had been removed and replaced by an hPGK-Puromycin cassette (pSIN-EF1a promoter-BGH polyA-hPGK-Puromycin). The vector was digested with ClaI and XbaI enzymes to remove the EF1a promoter, and the ATOH1 enhancer sequence was cloned using in-fusion HD cloning (Clontech). Subsequently, a sequence containing the human beta globin minimal promoter (hβ−globin) followed by a nuclear localization signal and EGFP (nls-EGFP) was amplified from the J2XnGFP plasmid DNA (a gift from Dr. Jane Johnson, UT Southwestern) and cloned downstream of the human ATOH1 enhancer by in-fusion HD cloning with the following primers (nls-EFGFP forward primer: 5′-GGTTGTTCAGGAATTCGATGGGCTGGGCATAAAAGT-3′; nls-EFGFP reverse primer: 5′- GCCCTCTAGAGAATTCAACTAGAGGCACAGTCGAGGC-3′) to obtain the following lentiviral construct: pSIN-hATOH1 enhancer-hβ-globin-nlsEGFP-bGH polyA-hPGK-Puromycin. For derivation of the ATOH1-EGFP-L10a line, the nlsEGFP was replaced with EGFP-L10a. Briefly, the pSIN-hATOH1 enhancer- hβ-globin-nlsEGFP-bGH polyA-hPGK-Puromycin construct was digested with SbfI/BsrGI enzymes to remove the nlsEGFP. An EGFP-L10a fusion fragment was amplified from a template plasmid (mPCP2-A box-s296, a gift from Dr. Nathaniel Heintz, The Rockefeller University) with the following in-fusion PCR primers (EGFP-L10a-forward primer: 5′- CATTTGCTTCTAGCCTGCAGGTCGCCACCATGGTGAG-3′; and EGFP-L10a reverse primer: 5′- CCGCTTTACTTGTACATTATCTAGATCCGGTGGATCC-3′) and cloned by in-fusion HD cloning to obtain pSIN-hATOH1 enhancer-EGFP-L10a-bGH polyA-hPGK-Puromycin. Lentiviral-mediated delivery was used according to established protocols to introduce the vectors into the RUES2 hESC line to generate transgenic lines. Two clonal lines with normal karyotype (Molecular Cytogenetics Core facility, MSKCC) and pluripotency characteristics were selected and expanded per transgenic line. Upon differentiation, EGFP expression was examined by microscopy. In the ATOH1-EGFP line, EGFP was broadly expressed in the nuclei of a subset of cells, while in the ATOH1-EGFP-L10a line, EGFP puncta were localized to the nucleoli, the site of ribosomal biogenesis (Figure 1—figure supplement 2A and B). A majority of the labeled cells coexpressed Ki67 (proliferation marker, data not shown), and the ATOH1 transcript was enriched in FAC-sorted EGFP+ cells from the ATOH1-EGFP line compared to EGFP- cells, and upon IP in the ATOH1-EGFP-L10a line compared to input. By contrast, housekeeping genes were not enriched (Figure 1—figure supplement 2A and B).
Fluorescence-activated cell sorting (FACS)
Request a detailed protocolCells at days 16, 19, 23, and 28 of differentiation were FAC-sorted on a BD FACSAriaII with BD FACSDiva 8.0.1 software (BD Biosciences) using a 100 μm nozzle and a 488 nm laser according to standard procedure. Briefly, differentiation cultures were washed once in Ca2+/Mg2+-free PBS and the cells were dissociated by incubation in Accutase (Millipore, SCR005) for 5 min at 37°C. Dissociated cells were resuspended in MACS buffer (see MACS sections) and put through a cell strainer (BD Falcon, 352235). Gating was performed on EGFP-positive and -negative control cells, and propidium iodide (Sigma), at an appropriate dilution, was used for dead cell exclusion. EGFP+ cells were collected in MACS buffer for further downstream analyses.
MACS of TAG1+ cells
Request a detailed protocolOn days 28–32 of differentiation, TAG1+ cells were isolated by MACS (Miltenyi Biotec) according to the manufacturer’s protocol. Briefly, cells were washed 1× in MACS buffer (0.5% BSA, 0.9% glucose in Ca2+/Mg2+-free PBS) and then gently scraped off from Transwell membranes in MACS buffer using a cell scraper (USA Scientific, CC7600-0220). Cells were collected by centrifugation (300 × g, 10 min). 1 ml Trypsin (1 g/ml)-DNase (100 mg/ml) solution (Worthington Biochemical, 3703 and 2139) was added to the pellet for 1.5 min at 37°C without disturbing the pellet. The Trypsin-DNase was then exchanged for 1 ml DNase (100 mg/ml), and the cell pellet was triturated using a fine-bore pulled glass pipette until a uniform cell suspension was obtained. The suspension was put through a 40 μm cell strainer (BD Biosciences, 352340) to remove remaining clumps. Serum containing ‘cerebellum’ media (see subsequent section, plus 10% horse serum, Invitrogen 26050-088) was added to inactivate the Trypsin-DNase, and the single-cell suspension was washed and spun (300 × g, 10 min) in 50 ml of CMF-PBS (PBS Thermo Fisher, 14190-250; 0.2% w/v glucose Millipore, G8769; 0.004% v/v NAHCO3 Millipore, S8761; 0.00025% Phenol Red Millipore, P0290). The cells were resuspended in 1 ml serum containing ‘cerebellum’ medium and spun in a tabletop centrifuge (300 × g, 5 min) in an Eppendorf tube. This step helped reduce dead cells and debris, which accumulated in the supernatant. The cell pellet was then resuspended in fresh ‘cerebellum’ medium plus serum and incubated in a bacterial dish in a 35°C, 5% CO2 incubator for 1 hr. This step allowed surface antigens to reappear after the enzymatic dissociation of cells. The number of live cells was counted, and the cells were incubated in TAG1 antibody in cerebellum medium (1:2) for 20 min at room temperature (RT), followed by 1× wash in MACS buffer, and then incubated in anti-mouse IgM microbeads (Miltenyi Biotec, 130-047-302) for 15 min at RT, and TAG1+ cells were sorted through MS columns (Miltenyi Biotec, 130-042-201) following the manufacturer’s description.
Purification of cerebellar neurons and glia and co-culture with human cells
Request a detailed protocolMixed cerebellar cultures or glial fractions were prepared from P0-1 pups and cultured in serum-free ‘cerebellum’ media (BME, Gibco; 2 mM L-glutamate, Gibco; 1% v/v BSA, Sigma; ITS liquid media supplement, Sigma; 0.9% v/v glucose, Sigma; 0.1 mg/ml Primocin, InvivoGen) as previously described (Baptista et al., 1994; Hatten, 1985). Mixed cerebellar cultures (no separation step) were plated at 2.8 × 106 cells/ml on poly-D-lysine (Millipore)-coated glass coverslips (Fisher, 12-545-81) placed in 24-well plates. Glia-only fractions (separated by Percoll gradient) were plated at 0.8 × 106 cells/ml. TAG1+ human cells isolated at DIV28 by MACS were then plated on top of mixed cerebellar cultures the next day (at 0.3–0.5 × 106 cells/ml). Glial cultures were allowed to form a monolayer (5–7 days) upon which DIV28 TAG1+ cells were then plated. As previously reported for their mouse counterparts (Hatten, 1985), the ratio of glia to human neurons in co-culture was crucial. At 1:4 (glia:neuron), the neurons induced detachment of glia from the culture dish and co-aggregated into attached spheres while at a 1:2 (glia:neuron) ratio the glia remained as a bed upon which neurons attached and extended processes. Half of the medium (serum-free cerebellum media) was replaced with fresh medium every 4 days. For the in vitro migration assays, glia were isolated from P0-1 pups as described above and plated on poly-D-lysine-coated glass coverslips. The next day, DIV28 TAG1-sorted human cells were plated on top in a 1:2 ratio (glia:neuron) in the above described medium. Cells were fixed at 36 hr in 4% PFA for 15 min at RT and processed for immunolabeling.
Gene expression analysis by RT-PCR and RT-qPCR
Request a detailed protocolmRNAs were extracted using the RNeasy Plus mini kit (QIAGEN) with on column genomic DNA elimination, and cDNAs transcribed with the Transcription First Strand cDNA Synthesis Kit (Roche) according to the manufacturer’s description. Reverse transcription and PCR were carried out according to the manufacturer’s descriptions using the HotStarTaq PLUS DNA Polymerase kit (QIAGEN) on a PTC-200 Peltier Thermal Cycler (MS Research). To catch cDNA amplification in the exponential phase, experiments were run for 30 cycles only as follows: initial heat activation at 95°C for 5 min, denaturation at 94°C for 30 s, annealing at primer-specific temperatures listed in Table 1 for 40 s, extension at 72°C for 1 min, final extension at 72°C for 7 min. RT-qPCR was performed using the SYBR Green method according to the manufacturers’ descriptions (Roche) using the default SYBR Green program on a Roche LightCycler 480 (Roche). All experiments were performed on at least three independent biological replicates, and each sample was run in triplicate for RT-qPCR. Data were normalized to housekeeping genes for comparisons and fold change calculated by the 2ΔΔCT method. Figure 1B and Figure 1—figure supplement 2A show technical replicates. As ATOH1 was not detected by the end of the RT-qPCR at cycle number 45 (but the housekeeping gene was) in the EGFP- samples, 45 was assigned as the ATOH1 cycle number value (CT) to signify a nondetected signal and calculate the fold change in the EGFP+ samples compared to EGFP- samples. For primer sequences and annealing temperatures, see Table 1. Primers were designed using Primer3 (https://bioinfo.ut.ee/primer3-0.4.0/) and validated to give single amplicons in a concentration-dependent manner at temperatures used.
Immunocyto/histochemistry
Request a detailed protocolCells grown on 1.5 thickness cover glass (Fisher) were fixed for 15 min at RT with 4% PFA and washed 3× PBS. Cells were then blocked in PBS containing 1% normal horse serum (Gibco) and 0.1% Triton for 1 hr and then incubated with primary antibodies in blocking solution for 1 hr at RT to overnight at 4°C. Cells were washed 3× PBS for 10 min, and Alexa Fluor conjugated secondary antibody incubations were performed in blocking solution for 1 hr at RT. Dapi was sometimes added as nuclear counterstain (1 μg/ml, Molecular Probes). For immunohistochemistry on frozen brain sections, human cerebella were fixed and embedded as described earlier. Mouse brains at E15.5, E17.5, P0, and P6 were fixed in 4% PFA overnight at 4°C, then cryoprotected in 20% sucrose overnight at 4°C and embedded in OCT. 14-μm-thick sagittal sections were prepared for all brains on a Leica CM 3050S cryostat. Frozen sections were thawed, postfixed for 10 min in 4% PFA at RT, and immunohistochemistry was carried out as above, except that the blocking solution contained 10% normal horse serum (Gibco) and 0.2% Triton in PBS. For analysis of transplantation experiments, brains were fixed in 4% PFA overnight at 4°C, washed in PBS, and embedded in 3% agarose. 50-μm-thick vibratome sections were postfixed with 4% PFA for 15 min at RT followed by blocking in 10% normal horse serum (Gibco) and 0.2% Triton in PBS overnight at 4°C. Primary antibody incubations were carried out in blocking solution for two nights at 4°C followed by extensive washes (4 × 15 min each) in PBS containing 0.1% Triton, and sections were then incubated in secondary antibodies overnight at 4°C. Sections/cells were mounted with ProLong Gold anti-fade mounting media (Invitrogen) and 1.5 thickness Fisherbrand cover glass. For antibody sources and dilutions, see Table 2.
Proliferation assay
Request a detailed protocolCell proliferation was measured by EdU incorporation using the Click-iT EdU Cell Proliferation Kit for Imaging (Invitrogen, C10338) according to the manufacturer’s description. On day 11 of differentiation, cells from the ATOH1-EGFP line were gently scraped off Transwell membranes and plated on poly-D-lysine (Millipore, A-003-E) and Laminin (Invitrogen, 23017-015)-coated glass coverslips, and differentiation was resumed as described previously (Figure 1D). On day 28, EdU was added to the culture medium according to the manufacturer’s description and cells were treated with either SAG (0.5 μM, Cayman Chemicals 11914) or DMSO. Cells were fixed with 4% PFA (15 min at RT) on day 30 and processed for Click-iT EdU detection and immunocytochemistry. See below for quantification and statistical analysis.
Transplantation of hPSC-derived cerebellar granule cells in the mouse cerebellum
Request a detailed protocolAll procedures were approved by the Rockefeller University Institutional Animal Care and Use Committee. TAG1+ cells were isolated by MACS on DIV28–32 as described earlier. Cells were counted and resuspended in ‘transplantation medium’ containing: BME, Gibco; 10% v/v horse serum, Invitrogen 26050-088; 0.9% v/v glucose, Sigma; and 0.5% w/v Fast Green for visualization of injection volumes. Cells were kept on ice while mouse pups were prepared for injections (up to 2 hr). Neonatal mouse pups (P1-4) were cryo-anesthetized, the heads were cleaned with alcohol, and a 1 μl single-cell suspension (5 × 106 cells/ml) was manually injected directly into the left cerebellar hemisphere using a glass microcapillary (Eppendorf, 5195 000.079) controlled with an Eppendorf CellTram Vario manual microinjector. The procedure was performed under a Zeiss OMPI-1FC surgical microscope (ENT) with Zeiss eyepieces. The capillary directly pierced through both skin and skull. The positioning of the capillary was guided by the left earlobe and Lambda, and the tip of the capillary was placed just under the skull to target the cerebellar surface. The capillary was held in place for 1 min after the completion of injection to minimize backflow after which it was gently pulled out and the pups were warmed up on a heating pad (Sunbeam) before being returned to their mothers. Injected animals were analyzed 48 hr post injection. Brains were dissected out and fixed overnight in 4% PFA. The entire cerebellum of each animal (N = 5 analyzed) was sectioned coronally at 50 μm thickness on a vibratome (Leica VT 1000S), and all sections were processed for immunohistochemistry using a human nuclear antigen antibody (HuNu) to detect human cells plus additional antibodies as described in the article.
TRAP and RNA sequencing
Request a detailed protocolTRAP was performed as previously reported (Heiman et al., 2014) on three independent differentiation experiments at DIV28. Briefly, polysomes were stabilized by adding 100 µg/ml cycloheximide to cell culture media for 10 min prior to homogenization of cells with polysome extraction buffer. Following clearing by centrifugation, supernatants were incubated at 4°C with end-over-end rotation for 16–18 hr with biotinylated Streptavidin T1 Dynabeads (Thermo Fisher, 65601) previously conjugated with GFP antibodies (Sloan Kettering Institute Antibody Core, HtzGFP-19C8 and HtzGFP-19F7). The beads were collected on a magnetic rack, washed, and resuspended in lysis buffer with β-mercaptoethanol (Agilent, 400753) to extract bound RNA from polysomes. RNA was purified using the RNeasy micro kit (QIAGEN, 74004). RNA quantity and quality were measured using an Agilent 2100 Bioanalyzer with the 6000 Pico Kit (Agilent, 5067-1513). Full-length cDNA was prepared using Clontech’s SMART-Seq v4 Ultra Low Input RNA Kit (634888) from 0.5 ng RNA with an RIN ≥ 8.3. 1 ng cDNA was then used to prepare libraries using the Illumina Nextera XT DNA sample preparation kit (FC-131-1024). Libraries with unique barcodes were pooled at equal molar ratios and sequenced on Illumina NextSeq 500 sequencer to generate 75 bp single reads, following the manufacturer’s protocol.
RNA sequencing analysis
Request a detailed protocolSequence and transcript coordinates for human hg19 UCSC genome and gene models were retrieved from the Bioconductor Bsgenome.Hsapiens.UCSC.hg19 (version 1.4.0) and TxDb.Hsapiens.UCSC.hg19.knownGene (version 3.2.2) Bioconductor libraries, respectively. Transcript expressions were calculated using the Salmon quantification software (Patro et al., 2017) (version 0.8.2) from raw FastQ files. Gene expression levels as TPMs and counts were retrieved using Tximport (Love et al., 2016) (version 1.8.0). Normalization and rlog transformation of raw read counts in genes were performed using DESeq2 (Love et al., 2018) (version 1.20.0). For visualization in genome browsers, RNA-seq reads were aligned to the genome using Rsubread’s subjunc method (version 1.30.6) (Liao et al., 2013) and exported as bigWigs normalized to reads per million using the rtracklayer package (version 1.40.6). Genes significantly enriched or depleted in IP over input were identified using DESeq2 with a Benjamini–Hochberg adjusted p-value cutoff of 0.05 and absolute log fold change cut-offs of both 0 and 2. The PsychEncode’s 'Human mRNA-seq processed data' as counts was retrieved from the PsychEncode’s portal (http://development.psychencode.org). GSEA of significantly enriched or depleted gene sets (absolute logFC > 2, padj<0.05) were performed using GSVA (version 1.34.0) (Hänzelmann et al., 2013) against DESeq2 normalized PsychEncode midbrain and cerebellum RNA-seq counts. Statistical significance of gene set enrichment within samples was determined using Limma’s geneSetTest with normalized median scaled expression values. Visualization of genes and gene sets as heatmaps was performed using the Pheatmap R package (version 1.0.10) (Subramanian et al., 2005). GO term enrichment was obtained for all genes differentially expressed between IP and input (absolute logFC > 0, adjusted p-value<0.05) using the Fisher test in the topGO Bioconductor package and ranked using the elim algorithm and functional annotation from the org.Hs.eg.db Bioconductor package (version 3.10). For comparison to mouse scRNA-seq data, cell-type marker gene sets were retrieved from Wizeman et al. and mouse symbols translated to human ortholog symbols. Genes significantly upregulated in IP over input (logFC > 0, adjusted p-value<0.05) were tested for enrichment of cell type markers using the GOseq Bioconductor package (Young et al., 2010).
Imaging
Request a detailed protocolSingle z-plane images (512 × 512 pixels) were acquired using an inverted Zeiss LSM 880 NLO laser scanning confocal microscope operated with ZEN imaging software (Zeiss) and fitted with a Plan-Apochromat 40×/1.4 NA objective oil immersion lens, a Nomarski prism, and HeNe and Argon lasers for excitation at 405, 488, 561, and 633 nm. Phase-contrast images in Figure 1—figure supplement 1D were acquired with a 20× objective on a Leica DMIL LED microscope fitted with a Leica MC120 HD camera.
Quantification and statistical analyses
Request a detailed protocolSample size estimations were based on prior pilot experiments. Independent experiments, performed on different days (cell culture experiments) or independent samples (e.g., individual mice), were considered biological replicates. Repeat measurements on the same sample were considered technical replicates. Samples were randomly allocated to control and treatment groups. For immunocyto-/histochemistry experiments, cells were manually counted in ImageJ (version 2.1.0/1.53c) on single z-plane confocal images (512 × 512 pixels) from three independent (in vitro) experiments (unless stated otherwise). Data were checked for normal distribution by Shapiro–Wilk’s test of normality and analyzed by parametric tests as described below using SPSS software (IBM). In Figure 1E, the mean percentages of ATOH1-EGFP+ cells were compared by ANOVA. The percentages of PAX6, EN2, and EGFP-positive cells were calculated per PAX6, EN2, and EGFP populations as indicated in Figure 1F. Multiple representative images per coverslip were quantified, and a total of 1787 EN2+, 3487 GFP+, and 2845 PAX6+ cells were counted. For SAG treatment experiments, the percentage of EdU/EGFP double-positive cells was calculated as a subset of the Dapi population. A total of 84 images (6209 EdU+ cells) were analyzed with five outliers removed. Outliers were defined as datapoints that were 1.5× outside of the interquartile range of box plots. Data were log transformed for normality, and the SAG and control treatment groups were compared by ANCOVA with the number of Dapi+ cells as a covariate. The number of TAG1+ cells was counted after MACS using a hemocytometer under a Leica microscope (Leica DMIL LED) and expressed as a percentage of the cells at the start of the sort (input). All other markers (NEUROD1, NeuN, Calretinin, PAX6, SOX2) were quantified on multiple confocal z-plane images in ImageJ from three independent experiments (except SOX2, n = 2) as described earlier. For the Tag1+ fraction, the number of cells quantified for each marker is indicated in the main text. For the SOX2 quantifications, a total of 2376 cells were counted. For the TAG1- fraction (Figure 2—figure supplement 1), a total of 2714 cells were counted. The percentage of NEUROD1-positive cells at the pial surface of the cerebellum was counted and normalized to the number of Dapi+ nuclei along the length of the pial surface. Multiple sagittal sections from two brains per time point and species were analyzed with the total numbers of cells/sample indicated in the description of Figure 5—figure supplement 1.
Data availability
Sequencing data have been deposited in GEO under accession code: GSE163710.
-
NCBI Gene Expression OmnibusID GSE163710. TRAP seq of the human pluripotent stem cell derived ATOH1 lineage.
-
Human mRNA seq processed data: Gene expression in countsID psychencode. Integrative functional genomic analysis of human brain development and neuropsychiatric risks.
-
elifesciencesSpecification of diverse cell types during early neurogenesis of the mouse cerebellum.https://doi.org/10.7554/eLife.42388.018
References
-
Rapid evolution of the cerebellum in humans and other great apesCurrent Biology 24:2440–2444.https://doi.org/10.1016/j.cub.2014.08.056
-
Cerebellar granule cells: insights into proliferation, differentiation, and role in medulloblastoma pathogenesisThe International Journal of Biochemistry & Cell Biology 41:435–445.https://doi.org/10.1016/j.biocel.2008.06.017
-
Perinatal cerebellar injury in human and animal modelsNeurology Research International 2012:858929.https://doi.org/10.1155/2012/858929
-
A gene with sequence similarity to Drosophila engrailed is expressed during the development of the neural tube and vertebrae in the mouseDevelopment 104:305–316.
-
Self-organization of human embryonic stem cells on micropatternsNature Protocols 11:2223–2232.https://doi.org/10.1038/nprot.2016.131
-
Efficient differentiation of human embryonic stem cells into functional cerebellar-like cellsStem Cells and Development 19:1745–1756.https://doi.org/10.1089/scd.2009.0498
-
Neuronal regulation of astroglial morphology and proliferation in vitroThe Journal of Cell Biology 100:384–396.https://doi.org/10.1083/jcb.100.2.384
-
Mechanisms of neural patterning and specification in the developing cerebellumAnnual Review of Neuroscience 18:385–408.https://doi.org/10.1146/annurev.ne.18.030195.002125
-
The generation of granule cells during the development and evolution of the cerebellumDevelopmental Dynamics 248:506–513.https://doi.org/10.1002/dvdy.64
-
The Subread aligner: fast, accurate and scalable read mapping by seed-and-voteNucleic Acids Research 41:e108.https://doi.org/10.1093/nar/gkt214
-
WNT-7a induces axonal remodeling and increases synapsin I levels in cerebellar neuronsDevelopmental Biology 192:31–44.https://doi.org/10.1006/dbio.1997.8734
-
Induction of medulloblastomas in p53-null mutant mice by somatic inactivation of Rb in the external granular layer cells of the cerebellumGenes & Development 14:994–1004.
-
Co-expression profiling of autism genes in the mouse brainPLOS Computational Biology 9:e1003128.https://doi.org/10.1371/journal.pcbi.1003128
-
Molecular markers of neuronal progenitors in the embryonic cerebellar anlageThe Journal of Neuroscience 26:12226–12236.https://doi.org/10.1523/JNEUROSCI.3493-06.2006
-
NeuN, a neuronal specific nuclear protein in vertebratesDevelopment 116:201–211.
-
Human iPSC-Derived Cerebellar Neurons from a Patient with Ataxia-Telangiectasia Reveal Disrupted Gene Regulatory NetworksFrontiers in Cellular Neuroscience 11:321.https://doi.org/10.3389/fncel.2017.00321
-
The Theory and Neuroscience of Cerebellar CognitionAnnual Review of Neuroscience 42:337–364.https://doi.org/10.1146/annurev-neuro-070918-050258
-
Maturation of Human Pluripotent Stem Cell-Derived Cerebellar Neurons in the Absence of Co-cultureFrontiers in Bioengineering and Biotechnology 8:70.https://doi.org/10.3389/fbioe.2020.00070
-
The role of the rhombic lip in avian cerebellum developmentDevelopment 126:4395–4404.
Article and author information
Author details
Funding
National Institute of Neurological Disorders and Stroke (1R21NS093540-01)
- Mary E Hatten
Rockefeller University (Pilot award)
- Hourinaz Behesti
- Mary E Hatten
Starr Foundation (Tri-Institutional Stem Cell Initiative Grant)
- Mary E Hatten
US Army Medical Research Acquisition Activity (W81XWH1510189)
- Mary E Hatten
Rockefeller University
- Mary E Hatten
Renate, Hans, and Maria Hofmann Trust
- Mary E Hatten
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Dr. Ali Brivanlou and Dr. Zeeshan Ozair (The Rockefeller University) for the provision of vital reagents and critical discussions throughout the project, Dr. Nathaniel Heintz and Jie Xing (The Rockefeller University) for the kind gift of the EGFP-L10a containing construct, helpful discussions regarding the adaptation of bacTRAP methodology for use in hPSCs and technical support. We also thank Dr. Bobak Mosadegh (Cornell University) for discussions on the use of transwells and minimization of culture variability, Dr. Jane Johnson (UT Southwestern) for the kind gift of the J2XnGFP plasmid, and the late Dr. Tom Jessell (Columbia University) for the provision of the TAG1 antibody. We are also grateful to staff at the Rockefeller Core facilities including imaging, flow cytometry, and high-throughput sequencing. Finally, we thank members of the Hatten lab for helpful discussions. This work was supported by NIH 1R21NS093540-01 (MEH), a Rockefeller University Center for Clinical and Translational Science Pilot grant (MEH, HB), a Tri-Institutional Stem Cell Initiative grant from the Starr Foundation (MEH), and a gift from the Renate Hans and Maria Hofmann Trust (MEH).
Ethics
Human subjects: Fixed de-identified human tissue were acquired from the Human Developmental Biology Resource (http://www.hdbr.org/) following institutional policies.
This study was performed in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health. All of the animals were handled according to approved institutional animal care and use committee (IACUC) protocol (#14746-H) of the Rockefeller University. All surgery was performed under hypothermia, and every effort was made to minimize suffering.
Copyright
© 2021, Behesti et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,943
- views
-
- 282
- downloads
-
- 17
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Developmental Biology
- Neuroscience
Exposure to ketone bodies in early development can reduce neurological impairments in a strain of the nematode C. elegans with PTEN defects.
-
- Developmental Biology
- Neuroscience
A finely tuned balance between excitation and inhibition (E/I) is essential for proper brain function. Disruptions in the GABAergic system, which alter this equilibrium, are a common feature in various types of neurological disorders, including autism spectrum disorders (ASDs). Mutations in Phosphatase and Tensin Homolog (PTEN), the main negative regulator of the phosphatidylinositol 3-phosphate kinase/Akt pathway, are strongly associated with ASD. However, it is unclear whether PTEN deficiencies can differentially affect inhibitory and excitatory signaling. Using the Caenorhabditis elegans neuromuscular system, where both excitatory (cholinergic) and inhibitory (GABAergic) inputs regulate muscle activity, we found that daf-18/PTEN mutations impact GABAergic (but not cholinergic) neurodevelopment and function. This selective impact results in a deficiency in inhibitory signaling. The defects observed in the GABAergic system in daf-18/PTEN mutants are due to reduced activity of DAF-16/FOXO during development. Ketogenic diets (KGDs) have proven effective for disorders associated with E/I imbalances. However, the mechanisms underlying their action remain largely elusive. We found that a diet enriched with the ketone body β-hydroxybutyrate during early development induces DAF-16/FOXO activity, therefore improving GABAergic neurodevelopment and function in daf-18/PTEN mutants. Our study provides valuable insights into the link between PTEN mutations and neurodevelopmental defects and delves into the mechanisms underlying the potential therapeutic effects of KGDs.