Abstract
The cerebral cortex is built during embryonic neurogenesis, a period when excitatory neurons are generated from progenitors. Defects in neurogenesis can cause acute neurodevelopmental disorders, such as microcephaly (reduced brain size). Altered dosage of the 1q21.1 locus has been implicated in the etiology of neurodevelopmental phenotypes; however, the role of 1q21.1 genes in neurogenesis has remained elusive. Here, we show that haploinsufficiency for Rbm8a, an exon junction complex (EJC) component within 1q21.1, causes severe microcephaly and defective neurogenesis in the mouse. At the onset of neurogenesis, Rbm8a regulates radial glia proliferation and prevents premature neuronal differentiation. Reduced Rbm8a levels result in subsequent apoptosis of neurons, and to a lesser extent, radial glia. Hence, compared to control, Rbm8a-haploinsufficient brains have fewer progenitors and neurons, resulting in defective cortical lamination. To determine whether reciprocal dosage change of Rbm8a alters embryonic neurogenesis, we overexpressed human RBM8A in two animal models. Using in utero electroporation of mouse neocortices as well as zebrafish models, we find RBM8A overexpression does not significantly perturb progenitor number or head size. Our findings demonstrate that Rbm8a is an essential neurogenesis regulator, and add to a growing literature highlighting roles for EJC components in cortical development and neurodevelopmental pathology. Our results indicate that disruption of RBM8A may contribute to neurodevelopmental phenotypes associated with proximal 1q21.1 microdeletions.
Introduction
During embryonic cortical development, progenitor populations undergo neurogenesis, producing excitatory neurons that ultimately contribute to the six-layered mature neocortex (Franco and Müller, 2013; Greig et al., 2013). In mice, neurogenesis begins at ∼E10.5 and continues through the end of embryogenesis. The earliest progenitors are neuroepithelial cells, which differentiate into radial glia by E11.5. Radial glial cells undergo self-renewal divisions, and as development proceeds, switch to neurogenic divisions to generate new radial glia and either a neuron or a neuron-producing intermediate progenitor (IP; Noctor et al., 2001; Englund et al., 2005; Kowalczyk et al., 2009; Vasistha et al., 2014). Defective progenitor proliferation or excessive apoptosis can skew neuron number and ultimately alter mature brain size and function (Silver et al., 2010; Lancaster et al., 2013; Marthiens et al., 2013; Insolera et al., 2014). Thus, aberrant neurogenesis can cause severe neurodevelopmental disorders including microcephaly (small brain size accompanied by cognitive impairment; Hu et al., 2014). Understanding the genes regulating neurogenesis is a high priority to elucidate the etiology of neurodevelopmental diseases.
Recent studies from our lab and others have demonstrated critical roles for the exon junction complex (EJC) in neurogenesis and neurodevelopmental pathology. The core EJC, composed of MAGOH, EIF4A3, and RBM8A, binds mRNA and regulates RNA metabolism (Kataoka et al., 2001; Palacios et al., 2004; Shibuya et al., 2004; Ashton-Beaucage et al., 2010; Roignant and Treisman, 2010). We showed previously that Magoh haploinsufficiency in mice causes microcephaly due to IP depletion and neuronal apoptosis (Silver et al., 2010; McMahon et al., 2014). Subsequent human genetic studies revealed copy number variations (CNVs) of EJC components, EIF4A3 and RBM8A, associated with intellectual disability and brain malformations (Nguyen et al., 2013). RBM8A is one of 16 genes in the 1q21.1 proximal region [breakpoint (BP) 2-BP3], located adjacent to the 1q21.1 distal region (BP3-BP4). Microdeletions and duplications of either region are associated with neurodevelopmental phenotypes including brain size abnormalities and autism spectrum disorders; however, the causative gene(s) remain unknown (Brunetti-Pierri et al., 2008; International Schizophrenia Consortium, 2008; Mefford et al., 2008; Rosenfeld et al., 2012). Additionally, compound RBM8A mutations (null plus a noncoding mutation) cause thrombocytopenia-absent radius (TAR) syndrome, a blood skeletal disorder in which 7% of patients also exhibit cognitive impairments and brain malformations (Albers et al., 2012; Rosenfeld et al., 2012; Nguyen et al., 2014). Although there is strong evidence that RBMA8A causes TAR syndrome, it remains undetermined whether RBMA8A also contributes to neurodevelopmental phenotypes.
Here we examined the requirement of Rbm8a in embryonic brain development. Using a novel conditional mouse allele, we show that Rbm8a depletion in the dorsal telencephalon causes severe microcephaly due to neurogenesis defects affecting progenitor proliferation, progenitor and neuron number, and apoptosis. Additionally, using zebrafish models and in utero electroporation of mouse embryos, we demonstrate RBM8A overexpression does not significantly perturb progenitor number or head size. Our findings show Rbm8a is critical for proper cortical development, revealing that RBM8A haploinsufficiency may be a major driver for microcephaly in individuals with proximal 1q21.1 deletions.
Materials and Methods
Mouse husbandry and generation of conditional Rbm8a allele.
All experiments were performed in agreement with the guidelines from the Division of Laboratory Animal Resources from Duke University School of Medicine and Institutional Animal Care and Use Committee. Plug dates were defined as E0.5 on the morning the plug was identified. The conditional targeting vector for ES cell targeting was designed and generated by the Transgenic Facility at Duke University Cancer center. Positive ES clones were selected by performing long-range PCR of both arms. For long-range PCR of 5′ arms, the following conditions were used: 94°C × 2 min (1×); 98°C × 15 s, 60°C × 15 s, 68°C × 7.5 min (40×); 72°C × 7 min. 5′ F1: ATGCCTCCCTTCTAAGACAGGCTG; 5′ R1: AAGGGTTATTGAATATGATCGGAATTGG. For long-range PCR of 3′ arms, the following conditions were used: 94°C × 2 min (1×); 98°C × 15 s, 60°C × 15 s, 68°C × 2.5 min (40×); 72°C × 10 min. 3′ F1:CATTCGCCTTCTTGACGAGTTCTTC; 3′ R1:GTCTGCTCTTCCAGCTCACAACTG. Positive clones were electroporated into C57BL/6J blastocysts, and the chimeras were mated to C57BL/6J females to obtain germline transmission. All experiments to analyze this line were done on an inbred C57BL/6J background. For genotyping Rbm8aloxp mice, the following conditions were used: 94°C × 2 min (1×); 94°C × 15 s, 58°C × 30 s, 72°C × 30 s (30×); 72°C × 7 min (1×). LoxF1: CGGACGTGCTGGATCTTCAC; LoxR1: GCACACAGACTCCCCATAGG. The following strain was acquired from The Jackson Laboratory: Emx1-Cre (B6.129S2-Emx1tm1(cre)Krj/J).
Western blot and qRT-PCR analysis.
E11.5 embryonic cortices were collected from Emx1-Cre and Emx1-Cre;Rbm8aloxp/+ mice and lysed in RIPA lysis buffer with protease inhibitors (Pierce). Cortical lysates were run on 4–20% precasted SDS-polyacrylamide gels (Bio-Rad). Gels were transferred onto nitrocellulose membranes and blotted using the following primary antibodies: rabbit anti-Rbm8a (1:200; Bethyl Laboratories), mouse anti-Rbm8a (1:200; Novus), and mouse anti-α-tubulin (1:10,000; Sigma). Blots were developed using ECL reagent (Pierce). Densitometry was performed using ImageJ. Final values were quantified by normalizing Rbm8a levels to loading controls (tubulin) and analyzed for significance using a Student's t test. E10.5, E11.5, E12.5, and E14.5 embryonic cortices were collected from C57BL/6J (wild-type), Emx1-Cre, and Emx1-Cre;Rbm8aloxp/+ embryos and RNA was extracted using the RNeasy kit (Qiagen). At E10.5 the entire cortex was used, whereas at all other ages, only dorsal neocortex was included. DNA was prepared according to the iScript kit (Bio-Rad). qPCR was performed in triplicates using TaqMan probes (Life Technologies): Rbm8a (Mm04214345_s1), Pax6 (Mm00443081_m1), and Gapdh (4352339E). For wild-type samples at different developmental stages, absolute qRT-PCR was performed. A standard curve was generated with a five serial 10-fold dilution of cDNA from an independent E14.5 wild-type embryo. For E10.5 and E12.5 Emx1-Cre and Emx1-Cre;Rbm8aloxp/+ samples comparative qRT-PCR was performed. Values were normalized to Gapdh control. For each genotype three embryos were examined.
Immunohistochemistry and quantification.
Brains were fixed overnight in 4% PFA at 4°C, followed by submersion in 30% sucrose until sinking, as previously described (Silver et al., 2010). Brain cryostat sections (20 μm) were prepared and stored at −80°C until use. Sections were permeabilized with 0.25% Triton X-100 for 10 min and blocked with MOM block reagent (Vector Laboratories) for 1 h at room temperature (RT). Sections were incubated with primary antibodies for 2 h at RT or overnight at 4°C. Sections were then incubated in species-appropriate secondary antibodies and Hoechst for 15 min at room temperature. The following primary antibodies were used: rabbit anti-CC3 (diluted 1:200; Cell Signaling Technology), rabbit anti-Tbr2 (1:1000; Abcam), rabbit anti-Pax6 (1:1000; Millipore), rabbit anti-Rbm8a (1:200; Bethyl Laboratories), mouse anti-Pax6 (1:50; Developmental Studies Hybridoma Bank), rabbit anti-Tbr1 (1:200; Abcam), rabbit anti-Foxp1 (1:200; Abcam), rabbit anti-Cux1 (1:100; Santa Cruz Biotechnology), mouse anti-TuJ1 (1:400; Covance), mouse anti-Flag (1:200; Sigma), and chicken anti-GFP (1:2000; Abcam). The following secondary antibodies were used: Alexa Fluor 488, 568, 594, 647 (1:200; Invitrogen). TUNEL staining was done in accordance with the manufacturer's directions using an ApopTag fluorescein in situ apoptosis kit (Millipore). High-magnification images were captured using a Zeiss Axio Observer Z.1 microscope coupled with an apotome. Cortical thickness was measured with Zen software. Cell quantification was performed with ImageJ. Three sections from anatomically comparable regions per embryo and three biological replicates from control (Emx1-Cre) and mutant (Emx1-Cre;Rbm8aloxp/+) were measured/quantified.
Primary progenitor analysis.
E11.5–E14.5 embryonic dorsal cortices were dissected and dissociated with 0.25% Trypsin-EDTA (Life Technologies) at 37°C for 5 min and triturated with flame-polished Pasteur glass pipette. Dissociated cells were washed once with RT PBS and plated on poly-d-lysine-coated coverslips in 6-well culture dishes. Cells were incubated in DMEM (1×) culture media (Life Technologies) supplemented with B27, N2, bFGF, and N-acetyl-l-cysteine (Sun et al., 2005) at 37°C with 5% CO2 for 3 h. Cells were then fixed with 4% PFA for 10 min before staining with any of the following: rabbit anti-Rbm8a (1:200; Bethyl Laboratories), mouse anti-Pax6 (1:50; Developmental Studies Hybridoma Bank), rabbit anti-Pax6 (1:1000; Millipore), mouse anti-TuJ1 (1:1000; Covance), and Hoechst. For EdU analysis, cortical cells were dissociated and allowed to settle for an hour in culture media. Cells were pulsed with 10 μm EdU for 30 min according to manufacturer's protocol (Life Technologies). Cells were incubated for 24 h then EdU detection was performed using Click-iT EdU Alexa Fluor 647 Imaging Kit (Life Technologies) and costained with rabbit anti-Ki67 (1:200; Vector Laboratories), mouse anti-TuJ1 (1:1000; Covance), and Hoechst. For all primary cell analysis, high-magnification images were captured using a Zeiss Axio Observer Z.1 microscope coupled with an apotome. Cell quantification was performed with ImageJ.
In utero electroporation.
Gene transfer into embryonic wild-type mice was done by in utero electroporation as described with the Harvard apparatus BTX electroporator (Saito, 2006). E13.5 embryos were electroporated with 2.7 μg of DNA in a total volume of 1 μl (2 μg of either pCAGGS-EX or pCAGGS-Flag-RBM8A accompanied by 0.7 μg of pCAGGS-EGFP). E16.5 embryonic brains were prepared as described above. Flag-RBM8A was subcloned into a pCAGGS vector from a pCL-Flag-RBM8A vector, a gift from Dr. Niels Gehring. Graphs include analyses of the following: pCAGGS-EX plus pCAGGS-EGFP (n = 9–11 sections, three to four brains) and pCAGGS-Flag-RBM8A plus pCAGGS-EGFP (n = 9 sections, three brains).
Zebrafish analysis.
The full-length cDNA of RBM8A was obtained from Life Technologies (Clone ID 10383, ultimate ORF collection) and was cloned into pCS2 vector. Capped mRNA was synthesized in vitro using the mMessage mMachine SP6 kit (Life Technologies; AM1340). We injected 100 ng of RBM8A mRNA into wild-type zebrafish embryos at the one- to two-cell stage. Injected embryos were either fixed at 2 d post fertilization (dpf) for immunostaining or fixed at 4.5 dpf for head-size measurement; the distance across the convex tips of the eyecups was measured and compared with an age-matched control group from the same clutch. Whole-mount immunostaining with HuC/D (postmitotic neurons) was performed for investigating neuronal development and head-size regulation at a cellular level. Embryos were fixed in 4% PFA overnight and stored in 100% methanol at −20°C. After rehydration in PBS, PFA-fixed embryos were washed in immunofluorescence (IF) buffer (0.1% Tween 20 and 1% BSA in PBS 1×) for 10 min at room temperature. The embryos were incubated in the blocking buffer (10% FBS and 1% BSA in PBS 1×) for 1 h at room temperature. After two washes in IF buffer for 10 min each, embryos were incubated in the primary antibody solution, 1:1000 anti-HuC/D (A21271; Life Technologies), in blocking solution overnight at 4°C. After two washes in IF buffer for 10 min each, embryos were incubated in the secondary antibody solution, 1:1000 Alexa Fluor goat anti-mouse IgG (A11001; Life Technologies), in blocking solution for 1 h at room temperature.
Results
Rbm8a is expressed in the developing mouse neocortex
In this study we asked whether the EJC protein, RBM8A, is required for cortical development. To determine the role for Rbm8a in corticogenesis, we first examined its expression in the developing mouse cerebral cortex. At E10.5, the onset of neurogenesis, we detected Rbm8a mRNA in E10.5 embryonic neocortices using qRT-PCR (Fig. 1A). Relative to E10.5, at both E12.5 and E14.5, Rbm8a levels increased significantly (p = 0.001 and p = 0.050). mRNA in situ hybridization of E14.5 mouse sagittal sections detected robust Rbm8a expression throughout the embryo, with an enrichment in the ventricular zone (VZ) and subventricular zone (SVZ) of the developing neocortex, where radial glia and IPs reside (Fig. 1B–D; Visel et al., 2004; Englund et al., 2005; Fietz and Huttner, 2011; genepaint.org). The relatively higher Rbm8a RNA levels in the VZ/SVZ layers compared with intermediate zone (IZ)/cortical plate (CP) layers is consistent with Rbm8a transcriptome data of E14.5 neocortical layers (Ayoub et al., 2011). Given the abundance of Rbm8a in the VZ/SVZ, we assessed RBM8A protein expression at E10.5, E12.5, and E14.5. Immunofluorescence of coronal neocortical sections revealed RBM8A expression throughout the neuroepithelium and neocortex at all three stages (Fig. 1E–H). We then used both immunofluorescence of E14.5 tissue sections and primary cell culture from E14.5 dorsal cortices to determine whether RBM8A is expressed in progenitors and/or neurons. This experiment revealed Rbm8a expression in both radial glia (Pax6 positive) and neurons (Tuj1 positive; Fig. 1H–P). In E10.5 and E12.5 tissue sections and in dissociated E14.5 cells, RBM8A localization was primarily nuclear (Fig. 1I–P). This is consistent with RBM8A's localization in immortalized cells and known function in regulating nucleocytoplasmic mRNA populations, supporting the reliability of this antibody (Bono and Gehring, 2011; Daguenet et al., 2012). These results demonstrate that RBM8A is highly expressed in the developing neocortex, consistent with a potential role in early stages of neurogenesis.
Rbm8a haploinsufficiency causes microcephaly
To determine the requirement of Rbm8a in cortical development, we generated a mouse carrying a loxp allele of Rbm8a (Rbm8aloxp/+; Fig. 2A; see Materials and Methods for details). We crossed Rbm8aloxp/+ mice with Emx1-Cre, which drives Cre expression in neural progenitors of the dorsal telencephalon beginning at E9.5 (Gorski et al., 2002; Chou et al., 2009; Sahara and O'Leary, 2009; cre.jax.org). Genotyping of genomic DNA confirmed two bands of the predicted size in heterozygous Emx1-Cre;Rbm8aloxp/+ mice but only one band in Emx1-Cre (referred to as control hereafter; Fig. 2B). Following Cre recombination, the loxp transcript retains exons 1 and 6 but is predicted to undergo nonsense-mediated decay (NMD). qRT-PCR analysis of E10.5 and E12.5 neocortices revealed a 50% reduction of Rbm8a mRNA levels in Emx1-Cre;Rbm8aloxp/+ brains at both stages (p = 0.05, and p < 0.0001, respectively; Fig. 2C). The reduced Rbm8a mRNA levels are consistent with haploinsufficiency and suggest the recombined Rbm8a loxp transcript is degraded by NMD. Western blot analysis of E11.5 Emx1-Cre;Rbm8aloxp/+ neocortices revealed a 70% reduction of Rbm8a protein, as evidenced with two independent polyclonal and monoclonal antibodies (p < 0.0001; Fig. 2D,E; data not shown). No truncated protein was evident at 3.8 kDa, the expected size if a protein product were produced by exons 1 and 6. Given RBM8A's function in translation (Shibuya et al., 2004), reduced RBM8A in Emx1-Cre;Rbm8aloxp/+ mice may interfere with translation of remaining Rbm8a transcripts, resulting in a 70% loss of Rbm8a protein. Consistent with the Western blotting and qRT-PCR results, RBM8A signal was also reduced by immunofluorescence of primary E12.5 cell cultures of Emx1-Cre;Rbm8aloxp/+ dorsal neocortices (Fig. 2F–I). These results indicate that conditional deletion of a single Rbm8a allele in the developing mouse neocortex induces a 70% reduction in RBM8A protein levels by E11.5.
Emx1-Cre;Rbm8aloxp/+ mice were viable, with normal Mendelian ratios at weaning (χ2 analysis, p = 0.564, n = 22 mice Emx1-Cre, n = 26 mice Emx1-Cre;Rbm8aloxp/+). We examined the cortical size of Emx1-Cre;Rbm8aloxp/+ mice at P12, by which time embryonic neurogenic events in the cortex are complete. In agreement with an essential role for Rbm8a in corticogenesis, P12 Emx1-Cre;Rbm8aloxp/+ brains were significantly microcephalic compared with their littermate controls (Fig. 2J,K). Importantly, Rbm8aloxp/+ mouse (control) without Cre exhibited no phenotype and were comparable to C57BL/6J mice (data not shown). Measurement of the dorsal surface area of the cortex showed Emx1-Cre;Rbm8aloxp/+ brains were ∼70% smaller than control (p < 0.0001; Fig. 2L). These results demonstrate that haploinsufficiency for Rbm8a causes severe microcephaly.
To determine the onset of microcephaly in Rbm8a haploinsufficient mice, we next examined the cortical thickness of embryonic brains using coronal sections. At E11.5, Emx1-Cre;Rbm8aloxp/+ and control brains exhibited comparable cortical thickness (p = 0.801; Fig. 3A,B,I). However, by E12.5, Emx1-Cre;Rbm8aloxp/+ cortices were 21% thinner than control (p = 0.005; Fig. 3C,D,I). By E13.5, Emx1-Cre;Rbm8aloxp/+ dorsal cortices were 50% thinner than control, and this was evident in both dorsal and more lateral cortical regions (p < 0.0001, both; Fig. 3E–I). By E14.5, Emx1-Cre;Rbm8aloxp/+ brains were markedly microcephalic, and coronal sections revealed the ventricle and dorsal telencephalon were severely reduced (Fig. 3J–O). Staining for radial glia, IPs, and neurons within the lateral region adjacent to the ganglionic eminence revealed the presence of these cell populations in the remaining neocortex (Fig. 3P–U). These analyses demonstrate that in Emx1-Cre;Rbm8aloxp/+ mice, the microcephaly arises during embryogenesis between E11.5 and E12.5, a period corresponding to progenitor proliferation and the onset of neuron production (Greig et al., 2013).
Rbm8a is required for proper number of radial glia and intermediate progenitors
Given the timing of microcephaly onset, the Rbm8a expression pattern, and the correlation between progenitor dysfunction and microcephaly, we next assessed progenitors in more detail in Emx1-Cre;Rbm8aloxp/+ cortices. We quantified radial glia number (Pax6-positive cells) in mutant and control brains. At E11.5, we observed similar densities of radial glia in control and Emx1-Cre;Rbm8aloxp/+ brains (p = 0.435; Fig. 4A,B,I), which is consistent with the normal dorsal cortical thickness at E11.5. Similarly, Pax6 mRNA levels in E10.5 Emx1-Cre;Rbm8aloxp/+ and control brains were also similar (100 ± 27.69% vs 114 ± 5.05%, p = 0.618). However, by E12.5 we noted a significant 54% reduction in the density of radial glia in Emx1-Cre;Rbm8aloxp/+ brains compared with control (p = 0.002; Fig. 4C,D,I). At E13.5, the remaining radial glia in the neocortex appeared disorganized (Fig. 4E,F). These results indicate Pax6-positive radial glial cells are initially produced in Rbm8a-haploinsufficient E11.5 brains, but are decreased significantly in number within a day of development. Hence the reduced cortical thickness of E12.5 mutants may be due in part to fewer radial glia. As radial glia produce IPs, we assessed the impact of radial glia depletion upon IP number. At E13.5 we observed an 87% reduction in the density of Tbr2-positive IPs in Emx1-Cre;Rbm8aloxp/+ cortices (p = 0.0001; Fig. 4G,H,J). These analyses indicate that Rbm8a regulates the number of both radial glia and intermediate progenitors in the developing neocortex. As both populations produce excitatory neurons, ultimately this would be predicted to have a severe impact on neuronal output.
Rbm8a haploinsufficient embryonic neocortices contain ectopic neurons
Having established Rbm8a is required for the proper number of radial glia and IPs, we next examined the impact of Rbm8a depletion upon neurons. Loss of radial glia can sometimes be explained by precocious production of neurons at the expense of progenitors, as seen in many microcephaly mutants, including Magoh+/− (Yingling et al., 2008; Silver et al., 2010; Xie et al., 2013; Insolera et al., 2014). To evaluate neurons, we first performed immunofluorescence of coronal sections using the neuronal marker Tuj1, which labels immature neurons of the cortical plate. At E11.5 we observed no obvious difference in the relative distribution of Tuj1-positive neurons of Emx1-Cre;Rbm8aloxp/+ and control brains (Fig. 5A,B). However, by E12.5 the layer of Tuj1-positive neurons was strikingly thicker in Emx1-Cre;Rbm8aloxp/+ cortices compared with control (Fig. 5C,D). This phenotype persisted at E13.5, at which stage Tuj1-positive cells were detectable throughout the thickness of the dorsal Emx1-Cre;Rbm8aloxp/+ neocortex (Fig. 5E,F). These results demonstrate that loss of Rbm8a leads to aberrant distribution of postmitotic neurons.
Altered neuron distribution could be caused by defective neuronal migration and/or increased neuron number. To further ascertain the cause, we quantified deep-layer neurons in E12.5 cortices using Tbr1 immunostaining (Fig. 5G–I). This analysis revealed a significant increase in Tbr1-positive neurons in the Emx1-Cre;Rbm8aloxp/+ neocortex compared with control (p = 0.001). We also used E12.5 dissociated primary cultures to quantify neurons and progenitors (Fig. 5J–L). The overall distribution of cell types was different between control and Emx1-Cre;Rbm8aloxp/+ (p < 0.0001). Compared with control, Rbm8a haploinsufficiency induced a 2.7-fold increase in neurons (p = 0.007; Tuj1+ only; Fig. 5L). Concomitantly, Pax6-positive radial glia were reduced by 1.5-fold in Emx1-Cre;Rbm8aloxp/+ cortices (p = 0.002; Pax6+ only; Fig. 5L). These differences were sustained at E13.5, when there was a 2.5-fold increase in neurons and a twofold decrease in radial glia (p = 0.003 and p = 0.036, respectively; data not shown). These results are consistent with the 50% loss of radial glia quantified in E12.5 Emx1-Cre;Rbm8aloxp/+ brain sections. Moreover, these fold changes are consistent with the presence of more neurons in Rbm8a-deficient cortices.
Rbm8a haploinsufficiency impacts the balance between proliferation and differentiation of neural progenitors
As Rbm8a haploinsufficiency altered the number of neurons and progenitors, we posited this could be due to premature neuronal differentiation of radial glia. At E11.5 radial glial cells undergo self-renewing divisions to produce new progenitors, and they also begin to produce neurons via neurogenic divisions. If loss of Rbm8a triggers premature neuronal differentiation, we expected that in Rbm8a mutant brains, a higher proportion of progenitors would exit the cell cycle. To test this possibility, we quantified cell-cycle exit of E11.5 progenitors using dissociated cells. Progenitors were pulsed with EdU for 30 min to label S-phase and analyzed 24 h later. We quantified both cycling progenitors (EdU+Ki67+) and cells that had exited the cell cycle (EdU+Ki67−; Fig. 6A–F). Compared with control, Emx1-Cre;Rbm8aloxp/+ showed a significantly lower fraction of EdU+Ki67+ cycling progenitors (p = 0.015; Fig. 6G). This indicates that Rbm8a deficiency causes a higher fraction of progenitors to exit the cell cycle, which is consistent with fewer radial glia observed in E12.5 cultures and sections. We also observed a significantly higher fraction of EdU+Tuj1+ cells in Emx1-Cre;Rbm8aloxp/+ progenitors (p = 0.002; Fig. 6H). This result indicates that the increased cell-cycle exit is associated with more neurons, consistent with neuronal analyses of cultures and tissue sections. These data show that Rbm8a functions in radial glia to prevent premature exit from the cell cycle.
Given the expression pattern of Rbm8a and its function in regulating radial glia number, we next asked whether proliferation was impacted in Rbm8a haploinsufficient brains. We previously showed that RBM8A knockdown in human cell culture induces an abnormally high fraction of G2/M cells (Silver et al., 2010); therefore, we assessed whether Rbm8a haploinsufficiency affects G2/M phases of E11.5 cortical progenitors by performing phospho-histone3 (PH3) staining. Consistent with prior published cell-culture data (Silver et al., 2010), we observed a significant increase in PH3 staining in E11.5 Emx1-Cre;Rbm8aloxp/+ brains, compared with control (p = 0.042; Fig. 6I–K). Increased PH3 staining has also been seen in other microcephaly mutants (Feng and Walsh, 2004; Yingling et al., 2008; Marthiens et al., 2013; Chen et al., 2014). These data indicate that in E11.5 Rbm8a haploinsufficient brains, progenitor proliferation is altered before any significant change in radial glia number or cortical thickness.
Rbm8a haploinsufficiency causes widespread apoptosis
Rbm8a could influence neuron number in two ways: by producing neurons at the expense of progenitors or by preferential apoptosis of progenitors leading to disproportionately higher neuronal fractions. Based on the cell-cycle exit analysis we favored the former mechanism; however, to address the second possibility we assessed apoptosis in Rbm8a haploinsufficient brains. To examine apoptosis, we performed immunofluorescence for cleaved-caspase3 (CC3) and TUNEL, mid- and late-stage apoptosis markers, respectively. At E11.5, there was a slight increase in apoptosis of Emx1-Cre;Rbm8aloxp/+ compared with control (Fig. 7A,B). This mild phenotype was consistent with normal cortical thickness and progenitor number observed at this age. Beginning at E12.5 much more extensive apoptosis was observed only in Emx1-Cre;Rbm8aloxp/+ dorsal cortices (Fig. 7C–G). Apoptosis was evident throughout the thickness of the neocortex, and was especially high in basal layers (CP and subplate), where neurons reside following migration.
Having established that extensive apoptosis occurs at E12.5 coincident with progenitor loss, we next examined which populations of cells undergo apoptosis. Based on the distribution of apoptotic cells, we posited that neurons undergo apoptosis. Indeed, colocalization of CC3 and Tuj1 revealed neuronal apoptosis in Emx1-Cre;Rbm8aloxp/+ mutants (Fig. 7H–K). We also detected CC3+Pax6+ cells in E12.5 brains, demonstrating some radial glial cells are eliminated by apoptosis (Fig. 7L–O). If higher neuronal fractions were explained by preferential progenitor apoptosis, we predicted to observe more dying radial glial cells than dying neurons. However, radial glia apoptosis was less extensive than neuronal apoptosis, with ∼7% of all Pax6+ radial glia undergoing apoptosis, compared with 21% of all Tuj1+ neurons (p = 0.005). One technical explanation for this result is that Pax6 is a nuclear marker and its epitope could be lost in apoptotic cells. However, the high abundance of apoptotic cells in basal layers of the neocortex, together with this colocalization analysis, indicates apoptosis is not skewed toward radial glia. Our developmental analyses of Rbma8a haploinsufficient embryos indicate proliferation and cell-cycle exit defects precede both altered cell populations and apoptosis. This suggests Rbm8a impacts radial glia and neuron number by first influencing progenitor proliferation. In the mutant, aberrant populations are then eliminated by apoptosis.
Rbm8a haploinsufficiency impacts neuronal layers
Having observed both proliferation and apoptosis defects early in neurogenesis, we next examined the impact of Rbm8a depletion upon neuronal layers at the end of neurogenesis. Specifically we asked whether there were differences in upper versus lower layers of P0 Rbm8a haploinsufficient brains. Analysis of P0 brains revealed a marked microcephaly in Emx1-Cre;Rbm8aloxp/+ mutants, consistent with analysis at other stages (Fig. 8A,B). Remarkably, most of the neocortex was missing in Rbm8a haploinsufficient brains, as evidenced by coronal tissue sections (Fig. 8C,D). To determine which neuronal populations were present in remaining neocortical tissue, we performed immunofluorescence of medial sections for Cux1 (layers II/III), Foxp1 (layers III–V), and Tbr1 (layer VI; Fig. 8E–J). Strikingly, Rbm8a-deficient brains had a significant dearth of both Cux1 and Foxp1-positive neurons, compared with control (p < 0.0001 each; Fig. 8E–H,K). In contrast, the density of Tbr1-positive layer VI neurons was not significantly reduced in Emx1-Cre;Rbm8aloxp/+ brains (Fig. 8I–K). However, although Tbr1 neuronal density was not affected, Tbr1-positive neurons were inappropriately distributed throughout the cortical tissue (p < 0.0001; Fig. 8I,J,L). The dramatic reduction in outer layer neurons is consistent with depleted progenitor populations seen at E12.5 and E13.5. Our analyses demonstrate an essential requirement for Rbm8a in neurogenesis and brain size, with microcephaly explained by aberrant progenitor proliferation and severe apoptosis.
RBM8A overexpression does not dramatically perturb neurogenesis or increase head size
The data presented herein show that Rbm8a is essential for proper neurogenesis and brain size. RBM8A is located in the proximal 1q21.1 CNV, in which both microdeletions and microduplications are associated with brain malformations including microcephaly and macrocephaly (Brunetti-Pierri et al., 2008; Mefford et al., 2008; Rosenfeld et al., 2012). Given this association and the report of micro/macrocephaly in some patients with TAR syndrome, we next asked whether Rbm8a overexpression could have a reciprocal impact upon neurogenesis to that observed with Rbm8a haploinsufficiency.
To first test this hypothesis, we performed in utero electroporation to express Flag-tagged human RBM8A in E13.5 dorsal cortices and subsequently analyzed neurogenesis at E16.5. Between E13.5 and E16.5, radial glial progenitors mainly produce neurons, which migrate to the cortical plate. Therefore we reasoned that if RBM8A overexpression impacts neurogenesis, at these stages dramatic differences in progeny distribution and/or number should be evident. Aberrant distribution can be indicative of defective neuronal migration or progenitor proliferation. We coelectroporated EGFP to mark electroporated cells with either RBM8A or empty vector, as a control (Fig. 9A–D). Using coronal sections of E16.5 electroporated brains, we quantified the distribution of EGFP-positive cells, which marks both electroporated progenitors and their progeny. Comparison of control versus RBM8A electroporated brains revealed no significant difference in the distribution of EGFP-positive cells in the CP, IZ, or SVZ/VZ layers (p = 0.999; Fig. 9D). Expression of Flag-tagged RBM8A was detectable in EGFP-positive cells as evidenced with either anti-Flag staining or anti-RBM8A staining (Fig. 9E–Q). This indicates RBM8A overexpression in E13.5 progenitors does not overtly impact distribution of progenitors and progeny.
Because Rbm8a haploinsufficiency caused depletion of both radial glia and IPs, we asked if RBM8A overexpression affects these populations. We quantified the proportion of EGFP-positive cells that colocalize with either Pax6 or Tbr2 (Fig. 9R–W,BB). Neither radial glial progenitor nor IP number were significantly impacted by RBM8A overexpression (p = 0.261 and 0.893, respectively; Fig. 9BB). Although we quantified no difference in distribution of transfected cells, it is possible there could be defects in either proliferation and/or apoptosis associated with RBM8A overexpression. For example, decreased apoptosis is associated with macrocephaly (Haydar et al., 1999). However, analysis of CC3 staining revealed no evidence of decreased apoptosis upon RBM8A overexpression (data not shown). We also examined whether RBM8A could impact proliferation by quantifying the fraction of EGFP-positive cells that colocalize with PH3 (Fig. 9X–BB). We observed no significant difference between control and RBM8A overexpression brains (p = 0.525). Due to the timing of electroporation, we cannot exclude the possibility that RBM8A overexpression impacts earlier stages of neurogenesis. However, the lack of any defects in EGFP distribution, progenitor fractions, or proliferation suggests overexpression of RBM8A in E13.5 progenitors does not grossly impact the number of radial glia or IPs.
To independently assess whether RBM8A overexpression leads to head size and neuronal defects, we used zebrafish to overexpress the human RBM8A mRNA at the one- to two-cell stage. We have shown previously that head size with concomitant neurogenesis defects in zebrafish embryos represents a useful surrogate for the evaluation of loci associated with neurocognitive traits in humans (Golzio et al., 2012; Beunders et al., 2013; Sugathan et al., 2014). For example, overexpression and knockdown of KCTD13, a locus within the 16p11.2 CNV, induces microcephaly and macrocephaly in zebrafish, respectively, and these data have been since corroborated independently by association studies in large cohorts (Golzio et al., 2012; Fromer et al., 2014). We overexpressed human RBM8A mRNA in fertilized eggs and evaluated gross morphometric and cell-specific characteristics of overexpressants at 2 and 4.5 dpf. Masked quantitative scoring of RBM8A RNA-injected embryos at 4.5 dpf showed no significant difference of head size between RBM8A RNA-injected embryos and controls (p = 0.46; Fig. 10A–C). Because RBM8A overexpression could perturb neurogenesis without affecting head size, we stained embryos at 2 dpf with anti-HuC/D, to mark postmitotic neurons (Grandel et al., 2006; Golzio et al., 2012; Beunders et al., 2013). We observed a normal bilateral distribution of HuC/D in the forebrain in RBM8A overexpressants that was indistinguishable from that of control embryos from the same clutch (p = 0.75; Fig. 10D–F). This is consistent with lack of any defect in early stages of neuronal differentiation or migration. Finally, to exclude the possibility that the absence of phenotype was due to degradation of the RBM8A RNA, we extracted total RNA from injected embryos and performed RT-PCR with human RBM8A primers. We could detect the RBM8A-injected RNA at 2 dpf and up to 4.5 dpf (Fig. 10G; data not shown). These data indicate that overexpression of RBM8A does not perturb neuronal differentiation at early developmental stages and is unlikely responsible for head size defects in patients with a duplication of RBM8A.
Discussion
In this study we demonstrate that the exon junction complex component Rbm8a is essential for embryonic neurogenesis and proper brain size in mice. Rbm8a haploinsufficiency causes microcephaly, due to progenitor depletion and massive neuronal apoptosis. Hence Rbm8a is a critical regulator for neural progenitor maintenance, differentiation, and brain size. Our study highlights RBM8A as a critical gene that may contribute to neurodevelopmental phenotypes associated with proximal 1q21.1 microdeletions.
Mechanisms of Rbm8a function in neurogenesis and microcephaly
We show that Rbm8a haploinsufficiency causes microcephaly due to depletion of both progenitors and neurons, ultimately leading to severe cortical lamination defects including loss of virtually all upper layer neurons. Our study elucidates two mechanisms to explain how these populations are perturbed. First we discovered that Rbm8a regulates radial glia proliferation, preventing premature cell-cycle exit and neuronal differentiation. As a result Rbm8a depletion causes accelerated neurogenesis by radial glia, resulting in a striking reduction in progenitor pools along with significantly more neurons. A second mechanism by which Rbm8a regulates brain size is via apoptosis, which in the mutant occurs largely in ectopic neurons, but is also evident to a lesser extent within radial glia. Apoptosis is a common microcephaly phenotype and is associated frequently with ectopic neuron production (Yingling et al., 2008; Kim et al., 2010; Silver et al., 2010; Marthiens et al., 2013; Insolera et al., 2014). Our data indicate that a combination of altered progenitor division and apoptosis of progeny contribute to the severe loss of dorsal cortex and ultimately microcephaly.
How do defects in radial glia proliferation initially impact progenitor and neuron number? We speculate that Rbm8a influences the balance of proliferative and neurogenic divisions of radial glia, leading to production of neurons at the expense of progenitors. As shown previously by our group, Rbm8a is required for proper cell division (Silver et al., 2010). More G2/M progenitors could indicate a faster overall cell cycle with a similar G2/M length or a prolonged G2/M with identical cell-cycle duration. As cell cycle has been linked to cell fate in the developing brain (Lange et al., 2009; Pilaz et al., 2009; Arai et al., 2011), it is plausible that cell-cycle defects could influence production of neurons and progenitors. Interestingly, while PH3+ progenitors were increased, the proportion of cycling Ki67+ progenitors decreased. This reduction in total cycling progenitors may be attributed to any of the following possibilities, which are not mutually exclusive: massive progenitor apoptosis, increased cell-cycle exit, and depletion of progenitor production. Future studies will help determine the detailed cellular and molecular mechanisms by which Rbm8a influences progenitor number as well as apoptosis. Of note, RBM8A joins a growing repertoire of RNA binding proteins implicated in these aspects of corticogenesis, as discussed in a recent review (DeBoer et al., 2013).
RBM8A forms a tight heterodimer with MAGOH as part of the EJC (Lau et al., 2003; Tange et al., 2005; Bono et al., 2006; Gehring et al., 2009); therefore, it is reasonable to predict that haploinsufficiency phenotypes would be similar. In support of this hypothesis, both Rbm8a and Magoh have similar expression patterns in the developing cortex. In addition, the temporal onset of Rbm8a and Magoh haploinsufficiency phenotypes, including apoptosis, reduced intermediate progenitors, and ectopic neurons are also similar (Silver et al., 2010). However, in contrast to Magoh haploinsufficiency models (Silver et al., 2010; McMahon et al., 2014), Rbm8a haploinsufficient mice display a more severe microcephaly. A careful comparison of neurogenesis phenotypes indicates that the major difference is Rbm8a mutants have more extensive apoptosis of both neurons and radial glia. But why might Rbm8a and Magoh haploinsufficiency induce varying severity of apoptosis? One possibility is protein levels; Emx1-Cre;Magohloxp/+ and germline Magoh+/− brains have 50% reduced MAGOH protein levels (Silver et al., 2010; McMahon et al., 2014), whereas Emx1-Cre;Rbm8aloxp/+ cortices showed a 70% RBM8A reduction. Clearly the dosage of these genes is important, thus a higher fold reduction of either Magoh or Rbm8a could be more detrimental for neurogenesis. Another explanation for the difference is redundancy; mice have only one Rbm8a homolog, yet they have two Magoh homologs, which are virtually identical in sequence (Silver et al., 2010, 2013; Singh et al., 2013). Although Magohb is expressed at lower levels than Magoh in the embryonic brain, it could compensate for reduced Magoh levels. A third possibility is that either Rbm8a and/or Magoh have functions independent of the EJC or each other. Indeed there are some indications of noncanonical functions for both genes, in regulating signaling (Muromoto et al., 2009; Togi et al., 2013). Future genetic studies perturbing additional EJC components will be valuable to determine whether RBM8A impacts neurogenesis via the EJC or not and if so which aspects of EJC function are dysregulated during cortical development.
Rbm8a dosage in neurodevelopmental disease
We demonstrate that Rbm8a haploinsufficiency perturbs brain size, which is consistent with human genetic studies showing a strong association between reduced RBM8A copy number and neurodevelopmental disorders. Copy number losses of RBM8A are significantly enriched in patients with neurodevelopmental disorders, including microcephaly (Albers et al., 2012; Nguyen et al., 2013). Moreover, microdeletions of the proximal 1q21.1 region, where RBM8A resides, are associated with microcephaly, intellectual disability, behavioral abnormalities, and seizures (Sharp et al., 2006; Brunetti-Pierri et al., 2008; International Schizophrenia Consortium, 2008; Mefford et al., 2008; Rosenfeld et al., 2012; Nguyen et al., 2013). Patients with TAR syndrome, a blood and skeletal disorder attributed to compound mutations of RBM8A, present with and without these brain malformations (Nguyen et al., 2013, 2014). We speculate that different genetic backgrounds and potential effects (trans- and cis-) from genes expressed within the larger 1q21.1 region (BP3-BP4) could contribute to the co-occurrence of blood/skeletal defects and neurodevelopmental phenotypes in TAR patients. For example, there is precedence for genetic interactions between genes located within CNVs and for trans-effects of deletion/duplications upon neighboring gene expression (Blumenthal et al., 2014; Carvalho et al., 2014). Regardless, our study highlights RBM8A as a strong causal candidate gene for brain-related phenotypes associated with the proximal 1q21.1 microdeletions. We obtained 70% reduction of Rbm8a expression in the mouse haploinsufficiency model, simulating a heterozygous deletion of RBM8A in humans. This new mouse model will be a valuable in vivo tool to determine how RBM8A loss impacts the etiology of neurodevelopmental phenotypes, including the cellular and behavioral basis.
Remarkably, dosage gains of RBM8A are also associated with neurodevelopmental phenotypes. CNV gains in RBM8A, as well as microduplications of the proximal 1q21.1 region (BP2-BP3), are associated with intellectual disability, autism, and brain malformations such as macrocephaly (Sharp et al., 2006; Brunetti-Pierri et al., 2008; International Schizophrenia Consortium, 2008; Mefford et al., 2008; Rosenfeld et al., 2012; Nguyen et al., 2013). We found that although loss of one Rbm8a allele led to severe neurogenesis defects, neurogenesis or head size were obviously affected by Rbm8a overexpression. Our data from two model organisms suggest that Rbm8a overexpression might not be a major contributor to macrocephaly in human patients with microduplications of the proximal region of 1q21.1. Instead, given that Rbm8a is also expressed in postmitotic neurons, we speculate that copy number gains might impact mature neuronal functions such as synapse development (Alachkar et al., 2013).
EJC components have collectively emerged as critical for neurodevelopmental disorders. In addition to RBM8A, altered gene dosage of additional EJC components, including EIF4A3, UPF3A, and UPF3B, are seen in patients presenting with intellectual disability and brain malformations (Tarpey et al., 2007; Addington et al., 2011; Nguyen et al., 2013). Noncoding EIF4A3 mutations cause Richieri-Costa-Pereira syndrome, a craniofacial disorder associated with neurological deficits (Favaro et al., 2014). In addition, coding mutations in UPF3B cause intellectual disability, schizophrenia, and autism (Tarpey et al., 2007; Addington et al., 2011). However, mechanisms underlying neurodevelopmental pathologies for these disorders still remain elusive. By uncovering a critical requirement for the EJC component RBM8A in cortical development, our study helps provide important insights into how this essential molecular complex contributes to the etiology of neurodevelopmental disorders.
Footnotes
This work was supported by National Institutes of Health (NIH) grants to D.L.S. (R01NS083897 and R00NS064197), AHA 13PRE17070078 Predoctoral Fellowship and Broad Graduate student Award to H.M., SFARI 239983 and NIH P50MH094268 to N.K., and National Alliance for Research on Schizophrenia and Depression Young Investigator Award to C.G., N.K. is a Distinguished Brumley Professor. We thank Autumn Rorrer for assistance with mouse husbandry, Chad Russel and Jason Willer for technical assistance for zebrafish experiments, members of the Silver lab for useful insights and helpful discussions, and Cagla Eroglu for critical reading of the manuscript.
The authors declare no competing financial interests.
- Correspondence should be addressed to Debra L. Silver, 224 CARL Building, Box 3175, 213 Research Drive, Durham, NC 27710. debra.silver{at}duke.edu