# Translational derepression of Elavl4 isoforms at their alternative 5' UTRs determines neuronal development

Tatiana Popovitchenko, Yongkyu Park, Nicholas F Page, Xiaobing Luo, Zeljka Krsnik, Yuan Liu, Iva Salamon, Jessica D Stephenson, Matthew L Kraushar, Nicole L Volk, et al.

► To cite this version:

Tatiana Popovitchenko, Yongkyu Park, Nicholas F Page, Xiaobing Luo, Zeljka Krsnik, et al. Translational derepression of Elavl4 isoforms at their alternative 5' UTRs determines neuronal development. Nature Communications, 2020, 11 (1), pp.1674. 10.1038/s41467-020-15412-8. hal-02565573

HAL Id: hal-02565573

https://univ-rennes.hal.science/hal-02565573

Submitted on 6 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
ARTICLE

https://doi.org/10.1038/s41467-020-15412-8

OPEN

<figure><img src="image_2.png" /></figure>

Check for updates

# Translational derepression of *Elavl4* isoforms at their alternative 5' UTRs determines neuronal development

Tatiana Popovitchenko<fnref n="1" />,², Yongkyu Park<fnref n="1" />, Nicholas F. Page<fnref n="1" />,³, Xiaobing Luo<fnref n="1" />, Zeljka Krsnik⁴, Yuan Liu<fnref n="1" />,², Iva Salamon<fnref n="1" />,²,⁴, Jessica D. Stephenson<fnref n="1" />, Matthew L. Kraushar<fnref n="1" />,², Nicole L. Volk<fnref n="1" />, Sejal M. Patel<fnref n="1" />, H. R. Sagara Wijeratne<fnref n="1" />, Diana Li<fnref n="1" />, Kandarp S. Suthar<fnref n="1" />, Aaron Wach<fnref n="1" />, Miao Sun<fnref n="1" />, Sebastian J. Arnold⁵, Wado Akamatsu⁶, Hideyuki Okano⁶, Luc Paillard⁷, Huaye Zhang<fnref n="1" />, Steven Buyske⁹, Ivica Kostovic⁴, Silvia De Rubeis⁹,¹⁰,¹¹,¹², Ronald P. Hart³ & Mladen-Roko Rasin<fnref n="1" />✉

Neurodevelopment requires precise regulation of gene expression, including post-transcriptional regulatory events such as alternative splicing and mRNA translation. However, translational regulation of specific isoforms during neurodevelopment and the mechanisms behind it remain unknown. Using RNA-seq analysis of mouse neocortical polysomes, here we report translationally repressed and derepressed mRNA isoforms during neocortical neurogenesis whose orthologs include risk genes for neurodevelopmental disorders. We demonstrate that the translation of distinct mRNA isoforms of the RNA binding protein (RBP), *Elavl4*, in radial glia progenitors and early neurons depends on its alternative 5' UTRs. Furthermore, 5' UTR-driven *Elavl4* isoform-specific translation depends on upstream control by another RBP, Celf1. Celf1 regulation of *Elavl4* translation dictates development of glutamatergic neurons. Our findings reveal a dynamic interplay between distinct RBPs and alternative 5' UTRs in neuronal development and underscore the risk of post-transcriptional dysregulation in co-occurring neurodevelopmental disorders.

<footnote n="1">¹Department of Neuroscience and Cell Biology, Rutgers University, Robert Wood Johnson Medical School, Piscataway, NJ 08854, USA. ²Graduate Program in Neurosciences, Rutgers University, Piscataway, NJ 08854, USA. ³Department of Cell Biology and Neuroscience, Rutgers University, Piscataway, NJ 08854, USA. ⁴Croatian Institute for Brain Research, Center of Research Excellence for Basic, Clinical and Translational Neuroscience, University of Zagreb, School of Medicine, Zagreb 10000, Croatia. ⁵Institute of Experimental and Clinical Pharmacology and Toxicology, Faculty of Medicine, Signaling Research Centers BIOSS and CIBSS, University of Freiburg, Freiburg D-79104, Germany. ⁶Department of Physiology, Keio University School of Medicine, Tokyo 160-8582, Japan. ⁷Univ Rennes, CNRS, IGDR (Institut de génétique et développement de Rennes)-UMR 6290, F-35000 Rennes, France. ⁸Department of Statistics, Rutgers University, Piscataway, NJ 08854, USA. ⁹Department of Psychiatry, Icahn School of Medicine at Mount Sinai, 1 Gustave L. Levy Pl, New York, NY 10029, USA. ¹⁰Seaver Autism Center, Icahn School of Medicine at Mount Sinai, New York, NY 10029, USA. ¹¹Mindich Child Health and Development Institute, Icahn School of Medicine at Mount Sinai, New York, NY 10029, USA. ¹²Friedman Brain Institute, Icahn School of Medicine at Mount Sinai, New York, NY 10029, USA. ✉email: roko.rasin@rutgers.edu</footnote>
The mature neocortex generates complex behaviors including high-level cognition and voluntary motor functions. In the prenatal neocortex, neural stem cells called radial glia progenitors (RG) sequentially give rise to distinct subpopulations of glutamatergic neurons that are critical for normal circuits and functions¹⁻³. RG and glutamatergic neurons execute transcriptional and post-transcriptional programs that drive their development¹⁻⁷. Post-transcriptional regulatory programs include alternative splitting of mRNA isoforms that determine cell fate⁸ and the regulation of mRNA translation (protein synthesis)²,³,⁵,⁷,⁹⁻¹⁵. These two mechanisms are interconnected, as alternative mRNA sequences may yield distinct interactions with translational machinery.

Translational control is facilitated by RNA-binding proteins (RBPs), which often act on regulatory elements found in 5' and 3' untranslated regions (UTRs)³,⁵,⁷,¹⁶⁻¹⁸. UTRs are sites where multiple trans-acting factors may compete or cooperate to affect translation¹⁹,²⁰. Alternatively spliced 3' UTRs not only influence translation efficiency, but also modify the localization and stability of mRNAs in neurons¹⁸. Alternative 5' UTRs are upstream of the start codon and are in a powerful position to dictate translational efficiency of a given mRNA¹⁷,²¹. Yet, the effect of alternative 5' UTRs on translation and how multiple RBPs interact with them to dictate neuronal development are unknown.

Two RBPs in the Embryonic Lethal, Abnormal Vision, Droysophila, homolog-Like 4 family, Elavl4 (HuD) and Elavl1 (HuR), regulate neocortical development and mRNA translation⁹,¹⁰,¹³. Elavl4 deletion in mice results in seizures and repetitive behaviors⁹ reminiscent of phenotypes observed in individuals with comorbid epileptic encephalopathy (EE) and autism spectrum disorder (ASD). These behavioral anomalies are accompanied by anatomically disrupted glutamatergic neurons⁹. Another RBP, the CUGBP Elavl Family Member 1 (Celf1/Cugbp1), is highly expressed during neocortical development⁵,²². While Celf1 mediates multiple post-transcriptional steps in other systems²³⁻²⁵, its role in developing neocortical progenitors and neurons remains unknown. In addition, it is still poorly understood if and how distinct mRNA isoforms are translationally regulated in space and time during neuronal development.

Here, we provide a genome-wide screen of mRNA isoforms during prenatal neocortical neurogenesis and their translational patterns of repression and de-repression. Remarkably, human orthologues of many of these translationally-regulated mRNAs are encoded by known risk genes for comorbid neurodevelopmental disorders (e.g., EE, ASD, and intellectual disability). We demonstrate that uncontrolled translation of *Elavl4* isoforms in RG and early neurons leads to defective development of glutamatergic neurons. Moreover, selective translation of distinct *Elavl4* isoforms via their alternative 5' UTRs determines neuronal development. We also report that Celf1 is an upstream regulator of 5' UTR-driven translation and plays a role in neuronal development through *Elavl4*. Finally, we provide the first evidence in both mouse and human neocortex that CELF1 and ELAVL4 proteins are markers of early and late ventral RG, respectively. Our work represents an advancement in understanding the neurodevelopmental mechanisms of RBP-regulated translation of distinct isoforms via alternative 5' UTRs. Moreover, we provide evidence for the implications of this novel mechanism for neurodevelopment.

## Results

RBP isoforms are translationally regulated during development. Our previous data has suggested that there are dynamic changes in transcript translation between E13 and E16 and especially of mRNA translation regulators⁵,¹⁰,¹¹,¹³. Yet, the exact mRNAs and associated isoforms affected by these changes were

unknown. To identify translationally-regulated isoforms, we performed an unbiased RNA-seq screen by polysome profiling of dissected wild type (WT) E13 (early neurogenesis) and E16 (mid-to-late neurogenesis) neocortices. We compared isoform levels in total mRNA and in translationally distinct fractions: 40S-60S-80S (including monosome) and polysomal fractions.

Multidimensional scaling revealed clusters distinguished by embryonic stage and fraction, confirming the reproducibility of replicates and the unique profiles of free and polysome-associated isoforms (Fig. 1a). The volcano plot of these data contains groups of isoforms that are potentially in diverse states of transcriptional or translational regulation (Fig. 1b). One group includes isoforms that are differentially expressed in their steady state levels, suggesting transcriptional regulation or changes in mRNA stability between E13 (Fig. 1b, gold dots) and E16 (Fig. 1b, green dots). The vast majority of isoforms showed steady state expression between E13 and E16 (30,073 isoforms; Fig. 1b, black dots). Most of these had comparable association with either monosome or polysome fractions (20,643 isoforms; Fig. 1c). However, within the 2114 isoforms that have steady state levels and altered their association with polysomes, 1143 mRNAs decreased association with polysomes between E13 and E16. We define these mRNAs as our set of potentially "translationally repressed mRNA isoforms" for further exploration (Fig. 1d, top). The 971 mRNA isoforms that were unchanged in total steady state levels but exhibited increased association with polysomes between E13 and E16 are defined as our set of "translationally derepressed mRNAs" (Fig. 1d, bottom).

To better understand the function of mRNAs that are translationally repressed or derepressed during neocortical neurogenesis, the 2114 steady state isoforms that change in their polysome association were further analyzed by gene ontology (GO) and KEGG using g:Profiler²⁶. The terms enriched in repressed mRNAs were diverse in comparison to derepressed mRNAs (Supplementary Fig. 1). Enriched GO and KEGG groups included "RNA binding" and "Ribosome" in agreement with our previous report²⁷. The clear enrichment of translational control pathways among translationally derepressed mRNAs at E16 implicated these mRNAs as being particularly dynamic between E13 and E16.

We also discovered several mRNAs encoding RBPs in our dataset to be translationally dynamic between E13 and E16 (Fig. 1e). Taken together, these findings reinforce the importance of precise translation control during neocortical development and especially as it proceeds from early to late neurogenesis³,⁵,⁷,¹⁰,¹¹,¹³,²⁷.

### *Elavl4 is derepressed as neocortical neurogenesis proceeds.*

Among the translationally derepressed RBP mRNAs was *Elavl4*. Elavl4 is widely used as an early marker of post-mitotic neuronal differentiation and Elavl4-deficient embryonic neurospheres and adult subventricular zone progenitors demonstrate reduced capacity to differentiate into neuronal lineages²⁸,²⁹. However, the role of this RBP in prenatal neocortical development is unknown. We first established where Elavl4 protein was expressed at key time points during neurogenesis. Immunohistochemistry (IHC) using a pan-isoform antibody showed that Elavl4 is enriched in the cortical plate (CP) where post-mitotic neurons typically reside at E13 and E16 (Fig. 1f, left). Elavl4 was also enriched in the ventricular zone (VZ), where neocortical RG reside, at E16, but not at E13 (Fig. 1f, right). Using laser capture microdissection (LCM) of the VZ, we identified total *Elavl4* mRNA levels to be unchanged between E13 and E16 (Fig. 1g), supporting the idea that *Elavl4* is translationally derepressed between these two stages. Elavl4 colocalized with Nestin and N-Cadherin in RG end-feet in the E16 VZ (Fig. 1h, left, arrowheads). Elavl4 IHC was
a

<figure><img src="image_3.png" /></figure>

b

<figure><img src="image_4.png" /></figure>

c

<figure><img src="image_5.png" /></figure>

d

<figure><img src="image_6.png" /></figure>

e

<table><caption>Translation of RBPs at E16</caption><thead><tr><th>Derepressed</th><th>Repressed</th></tr></thead><tbody><tr><td>CNOT7</td><td>CNOT3</td></tr><tr><td>CPSF3</td><td>CNOT8</td></tr><tr><td>EIF3H</td><td>CSTF2</td></tr><tr><td>EIF4A1</td><td>DDX1</td></tr><tr><td>ELAVL4 (HUD)</td><td>FXR2</td></tr><tr><td>FAM98A</td><td>HNRNPA0</td></tr><tr><td>HNRNPA1</td><td>HNRNPAB</td></tr><tr><td>MBNL2</td><td>HNRNPD</td></tr><tr><td>NANOS3</td><td>HNRNPH3</td></tr><tr><td>PAPOLG</td><td>HNRNPK</td></tr><tr><td>STAU1</td><td>HNRNPR</td></tr><tr><td>STAU2</td><td>KHSRP</td></tr><tr><td>TNRC6A</td><td>MEX3C</td></tr><tr><td>U2AF1L4</td><td>MOV10</td></tr><tr><td>YTHDC2</td><td>PCBP3</td></tr><tr><td></td><td>PUF60</td></tr><tr><td></td><td>QKI</td></tr><tr><td></td><td>RALY</td></tr><tr><td></td><td>SRSF11</td></tr><tr><td></td><td>SYMPK</td></tr></tbody></table>

f

<figure><img src="image_7.png" /></figure>

g

<figure><img src="image_8.png" /></figure>

h

<figure><img src="image_9.png" /></figure>

i

<figure><img src="image_10.png" /></figure>

j

<figure><img src="image_11.png" /></figure>

k

<figure><img src="image_12.png" /></figure>

undetectable in Nestin and N-Cadherin-expressing RG in the VZ of *Elavl4* KO, demonstrating specificity of the antibody (Fig. 1h, right)⁹,²⁸.

Since the developing neocortex has a heterogeneous cell population, we performed fluorescence-activated cell sorting (FACS) on *Nestin-GFP* transgenic mice to selectively isolate RG that typically express *Nestin*³⁰,³¹ (Supplementary Fig. 2a). We also isolated intermediate progenitors and early neurons that typically express *Tbr2* using *Tbr2-GFP* transgenic mice³²,³³. The mRNA levels of *pan-Elavl4* mRNA were unchanged between those two stages in both FACS populations (Fig. 1i; Supplementary Fig. 2b), mirroring the LCM results obtained in the entire VZ (Fig. 1g). However, *Elavl4* protein significantly increased with time in *Nestin-GFP+* RG (Fig. 1j, k) with an upward trend in *Tbr2-GFP+* cells (Supplementary Fig. 2c). These data confirm that

*Elavl4* is translationally derepressed as neurogenesis proceeds, especially in RG.

**Distinct 5' UTRs dictate neocortical expression of *Elavl4*.** *Elavl4* is spliced into 4 different isoforms in mice that have alternative 5' UTRs and a common 3' UTR (Fig. 2a). We took the advantage of the isoform-specific 5' UTRs to differentiate between their expression in vivo using qRT-PCR (Fig. 2b). *Atxn1* mRNA is expressed in both progenitors and early neurons and was used as a positive control; its mRNA levels increased between E13 and E16 in *GFP+* cells obtained from *Nestin-GFP* neocortices. *Elavl4-v2* and *-v3* steady state levels were maintained between the two FACS populations at both stages, while *Elavl4-v1&4* and *Elavl4-v4* increased in *Nestin-GFP+* cells (Fig. 2b). We confirmed these findings with *Atxn1* and 5' UTR -specific fluorescent in situ
**Fig. 1** *Elav4* is translationally repressed during early neocortical neurogenesis. **a** Multidimensional scaling of neocortical tissue samples from E13 and E16 input (total RNA), 40S-60S-80S, and polysome fractions (n = 3 spins). Based on RNA-seq transcript abundances, M1 and M2 represent the first two dimensions. **b** Volcano plot of differentially expressed isoforms (gold and green dots) between E13 and E16 from input RNA-seq. Isoforms that are significantly higher at E13 are demarcated by gold circles, versus those higher at E16 are demarcated as green circles. Unchanged isoforms are demarcated as black circles. **c** Venn diagram of mRNAs that do not significantly change in their total expression levels (blue circle) compared to those that differ in association with 40S-60S-80S (purple circle), and/or polysomes (green circle). **d** Schematic: mRNA transcripts selected for further analysis were translationally repressed (top) or derepressed (bottom); these transcripts change from association with the 40S-60S-80S fraction (top) to the polysome fraction (bottom). **e** Translationally repressed and derepressed RNA-binding proteins (RBPs) at E16. **f** IHC for pan-Elav4 (green) in WT E13 and E16 neocortices (n = 12 animals per age). Arrowheads denote detectable Elav4 signal in VZ at E16. DAPI shown in blue. Scale bar: 100 µm. **g** pan-Elav4 mRNA levels determined by qRT-PCR at E13 and 16 from LCM dissected VZ (n = 6 animals). Data represent the mean and SD. Data normalized to Gapdh and then to E13. Statistics: Student's t-test. n.s. = not significant. **h** IHC for pan-Elav4 (green), N-Cadherin (red) and Nestin (blue) in WT and Elav4 KO E16 VZ (n = 3 animals). Arrowheads denote end feet. Scale bar: 50 µm. **i** Relative pan-Elav4 mRNA levels from FACS-sorted GFP+ cells determined by qRT-PCR (n = 3 FACS sorts per age and strain). Data represent the mean and SD. Data normalized to Gapdh and then to E13. Statistics: Student's t-test. n.s. = not significant. **j** Western blot analysis of Elav4-protein (top) from E13 and E16 Nestin-GFP+ FACS cells. Gapdh (bottom) was used as loading control. IB = immunoblotting. **k** Densitometry quantification of j. Data represent the mean and SEM. Data normalized to Gapdh and then to E13. Statistics: Welch's t-test, 2-tailed, unpaired. \*p < 0.05.

hybridization (FISH) probes in the E13 and E16 VZ (Fig. 2c). We further found increased association of distinct *Elav4* isoforms with heavy polysomes at E16 compared to monosome or light polysome fractions (Supplementary Fig. 2d, e); *Elav4-v3* could not be reliably measured given low expression. These findings corroborated E16 translational derepression of Elav4, inferred from RNAseq data, and the divergence between mRNA and protein levels (Fig. 1). Collectively, these results indicate that precise translational control of distinct *Elav4* isoforms occurs as neocortical neurogenesis proceeds.

We next aimed to understand Elav4 isoform-specific protein expression. As there are no isoform-specific Elav4 antibodies available, we co-transfected E12 and E16 RG with CAG-GFP and constructs bearing specific *Elav4* isoform 5' UTRs upstream of the *Renilla* coding sequence (Fig. 2d, f). In the VZ at E13, only sporadic *Elav4-v3* 5'UTR-Renilla was expressed and no *Elav4-v2* or *-v1&4* 5' UTR-Renilla expression was observed (Fig. 2d). In contrast, *Elav4-v1&4* 5' UTR-Renilla was expressed in post-mitotic neurons (Fig. 2e), which was expected since Elav4 is a known post-mitotic marker. *Elav4-v2* and *-v3* 5' UTR-Renilla were absent from post-mitotic neurons (Fig. 2e), suggesting that *Elav4-v1&4* are the isoforms detected by the Elav4 antibody in early neurons⁹. Indeed, we confirmed that the GFP in the *Elav4-GFP* transgenic line previously used by our lab⁹ is inserted in place of *-v4's* exon 1 (Fig. 2a) and is expressed in early post-mitotic neurons of developing neocortices (Supplementary Fig. 7d).

Remarkably, *Elav4-v2* and *-v3* 5' UTR-Renilla were abundantly expressed in the VZ at E17 (Fig. 2f), a stage when we detected Elav4 protein in the VZ (Fig. 1f). *Elav4-v1&4* 5' UTR-Renilla was only sporadically observed in VZ mitotic cells at this time (Fig. 2f). *Elav4-v3* 5' UTR-Renilla expression was low or absent in migrating neurons in the intermediate zone (IZ), while *Elav4-v2* and *-v1&4* 5' UTR-Renilla remained expressed in neurons (Fig. 2g). These data are in line with the polysome fractionation experiment as well as with the expression of Elav4 protein in the VZ, corroborating specificity of Elav4 IHC. Utilizing this approach, we revealed isoform and cell-type specific regulation of *Elav4* translation through its 5' UTRs in the developing neocortex.

**Elav4** isoforms are bound by Celf1 in developing neocortices. Based on these results, we hypothesized that a translational repressor acts on *Elav4* at its isoform-specific alternative 5' UTRs. We first tested whether the 5' UTRs of all translationally repressed or derepressed mRNAs between E13 and E16 bear specific binding motifs (Supplementary Fig. 3a, b). We found the binding motif of

RBP Celf1 to be present in 4.71% of the 5' UTRs of translationally derepressed mRNAs, including *Elav4s'* (Supplementary Fig. 3a). *Celf1* mRNA expression has previously been reported in developing neocortices⁵,²². We observed that Celf1 protein is high in Pax6+ RG at E13 (Fig. 3a, left), but decreases by E16 (Fig. 3a, right). Western blot on LCM-dissected VZ and CP corroborated these observations in the VZ and showed an increase in Celf1 in the CP from E13 to E16 (Fig. 3b, c, Supplementary Fig. 3c). There was no IHC signal in *Celf1* KO VZ at E13 (Supplementary Fig. 3d), confirming the specificity of the Celf1 antibody. Celf1 colocalized with Nestin-GFP in the VZ at E13 and decreased by E16 contrasting Elav4 IHC (Fig. 3d). The contrasting VZ expression patterns of Celf1 and Elav4 support the hypothesis that Celf1 translationally represses *Elav4* translation in vivo.

To assess whether Celf1 binds *Elav4* isoforms, RNA immunoprecipitation (RIP) was performed on E13 and E16 neocortices using RIP-certified Celf1 antibody and corresponding IgG control (Fig. 3e). The resulting RNAs were assessed for the expression levels of all *Elav4* isoforms along with negative controls (*β-actin* and *Gapdh*). Celf1 bound *Elav4* transcripts in E13/E16 neocortices (Fig. 3e) at different levels compared to the input. These findings suggest that Celf1 may regulate *Elav4* mRNA translation in developing neocortices on an isoform-specific level.

Therefore, our next goal was to understand how Celf1 acts on *Elav4* mRNA, as RBPs can affect all steps of post-transcriptional processing from splicing to translation (Supplementary Fig. 4a). Celf1 exists in at least two isoforms: “long” representing the full-length protein and a “short” isoform with a 27 amino acid truncation at the N-terminus (Supplementary Fig. 4b, c). Functional differences between the two isoforms have not been previously established. Structurally, the RNA recognition motif 1 (RRM1) near the N-terminus is truncated in Celf1 short (Celf1S), while RRM3 near the C-terminus is shared between both isoforms. Comparison of RRM1³⁴ and RRM3³⁵ structures suggests their binding modes have meaningful differences.

To investigate functional differences between the two isoforms, we turned to the Neuroblastoma (N2a) cell line for over-expression (OE) experiments. We first confirmed that Celf1 binds *Elav4* in N2a cells (Fig. 3f). Upon OE of *Celf1L* or *Celf1S* (Supplementary Fig. 5a), there were no changes in *Elav4* mRNA steady state levels (Fig. 4a, Supplementary Fig. 5b), splicing (Supplementary Fig. 5b), or decay (Supplementary Fig. 5d). However, *Celf1S* suppressed Elav4 protein levels while *Celf1L* did not (Fig. 4b), suggesting that the shorter N-terminus imbues *Celf1S* with its repressive function. Combined with our RIP data, these results suggest that Celf1 suppresses *Elav4* translation.
<figure><img src="image_13.png" /><figcaption>b</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>c</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>d</figcaption></figure>

IUE @ E12 → Analysis @ E13

<figure><img src="image_16.png" /><figcaption>f</figcaption></figure>

IUE @ E16 → Analysis @ E17

<figure><img src="image_17.png" /><figcaption>e</figcaption></figure>

<figure><img src="image_18.png" /><figcaption>g</figcaption></figure>

<figure><img src="image_19.png" /></figure>

To test whether *Celf1* acts on *Elavl4* mRNA through the 3' UTR or 5' UTR, we cloned the 3' UTR shared amongst all *Elavl4* isoforms downstream of *Luciferase* (*Luc*) and the 5' UTRs specific to *Elavl4-v2*, -v3, -v1&4 (the *Elavl4-v1* 5' UTR is indistinguishable from the -v4 5' UTR; Fig. 2a) upstream of *Renilla* (*Ren*). LUC or

REN relative light units (RLU) represent protein expression levels and were normalized to *Luc* or *Ren* mRNA levels, respectively. The mRNA to protein ratio then serves as a proxy for translation level¹¹,¹³. We found that *Celf1S* did not have an effect on translation mediated by the 3' UTR (Fig. 4c) nor on *Elavl4-v2*
**Fig. 2** *Elavl4* isoforms are translationally derepressed in the VZ by mid-neurogenesis. **a** Schematic of the four *Elavl4* isoforms in mice. All isoforms share a common 3' UTR sequence (black). Minor differences are in coding regions (exons in gray) and major differences in 5' UTRs (distinguished by color). NCBI Accession numbers for each variant is indicated underneath common name and the location of qRT-PCR probes used is shown below (red). **b** Relative mRNA levels of positive control *Atxn1* and the different *Elavl4* isoforms from FACS-sorted Tbr2-GFP+ and Nestin-GFP+ cells determined by qRT-PCR (n = 3 FACS sorts per age). Data represent the mean and SEM. Data normalized to *Gapdh*. Statistics: Two-way ANOVA with Tukey's. \*p < 0.05. **c** Fluorescent in situ hybridization of E13 and E16 VZ for *Nestin* (blue), *Atxn1* (green, left), *Elavl4*-v2 (red, left), *Elavl4*-v3 (red, right), and *Elavl4*-v4 (green, right). Scale bar: 20 µm. **d, f** E13 and E17 neocortices were in *utero* electroporated (IUE) at E12 (**d**) or E16 (**f**), respectively. VZ was transfected with CAG-GFP (green) and either *Elavl4*-v2-5' *UTR-Renilla*, *Elavl4*-v3-5' *UTR-Renilla* or *Elavl4*-v1&4-5' *UTR-Renilla* (red) (n = 3 or 4 animals, respectively). The merged image is on the far right. Arrowheads represent colocalization in VZ. **e** and **g** show zoomed-in regions above VZ as indicated. Arrows represent colocalization in cortical plate (CP), intermediate zone (IZ) and subventricular zone (SVZ). DAPI is in blue. **d, f** scale bar: 40 µm. e, g scale bar: 20 µm.

<figure><img src="image_20.png" /></figure>

**Fig. 3** *Celf1* binds *Elavl4* developing neocortices. **a** Celf1 (red) and Pax6 (green) VZ, IHC at E13 and E16 (n = 3 animals per age). DAPI shown in blue. **b** Western blot of Celf1 protein (upper panel) expression in VZ and CP from LCM neocortices at E13 and E16 (n = 3 animals per age). The Celf1 upper band is putatively Celf1 Long (Celf1L); the Celf1 lower band is putatively Celf1 Short (Celf1S), indicated at left. Gapdh was used as a loading control (bottom panel). **c** Densitometry quantification of **b**. Data represent the mean and SEM. Data normalized to *Gapdh* and then to E13. Statistics: Student's t-test. \*p < 0.05. **d** Celf1 (red, top) or Elavl4 (red, bottom) IHC on *Nestin*-GFP (green) neocortices at E13 and E16 (n = 4 animals per age). **e** RNA immunoprecipitation (RIP) from E13 and E16 neocortices. Relative mRNA levels were determined by qRT-PCR from Celf1 vs. IgG RIPs (n = 4 separate RIPs, 6 neocortices per age). *Pan-Elavl4* and specific *Elavl4* isoforms are shown in addition to *ActB* and *Gapdh* as non-binding controls. **f** RIP from N2a cells (n = 3 RIPs from separate transfections). Data represent the mean and SEM. Data normalized to *Gapdh*. Celf1 was used as positive control<sup>60</sup>. 18S was used as negative control. Statistics: Two-way ANOVA with Tukey's. \*\*\*\*p < 0.0001.

5' UTR based translation (Fig. 4d). However, *Celf1S* suppressed translation of *Ren* via the *Elavl4*-v3 and -v1&4 5' UTRs (Fig. 4e, f; Supplementary Fig. 5f). We then sought to narrow down the effective regulatory region(s) in *Elavl4*-v3 and -v1&4 5' UTRs. Using the same translational assay and a series of deletion constructs, we identified the region between nucleotides 13-24 of *Elavl4*-v3 5' UTR (Fig. 4e) and the region between nucleotides 273-284 of *Elavl4*-v1&4 (Fig. 4f) from the 5' UTR start site to be necessary for Celf1 regulation.

Celf1 binding sites in *Elavl4* 5' UTRs were verified in a RIP experiment. When Celf1 binding sites were either deleted or mutated, the Celf1 binding of 5' UTRs was significantly diminished

(Fig. 4g, h). These data indicate that Celf1 can suppress *Elavl4* mRNA translation in an isoform-specific manner and through specific regulatory elements in the 5' UTRs of *Elavl4*-v3 and -v1&4.

**Celf1 silencing derepresses Elavl4 translation.** Robust translation occurs when multiple ribosomes (a polysome complex) simultaneously translate a single mRNA. If *Celf1S* OE results in decreased translation of the *Elavl4* protein, then it is possible that some aspect of *Elavl4* polysome incorporation could be stymied in the presence of Celf1. Polysome fractionation under *Celf1S* OE in N2a cells did not yield significant shifts of *Elavl4*-v4
<figure><img src="image_21.png" /><figcaption>b</figcaption></figure>

<figure><img src="image_22.png" /><figcaption>c</figcaption></figure>

<figure><img src="image_23.png" /><figcaption>d</figcaption></figure>

<figure><img src="image_24.png" /></figure>

<figure><img src="image_25.png" /><figcaption>e</figcaption></figure>

<figure><img src="image_26.png" /></figure>

GUGAGAAGCGAC

<figure><img src="image_27.png" /><figcaption>f</figcaption></figure>

UGUCAGCUUUUA

<figure><img src="image_28.png" /><figcaption>g</figcaption></figure>

<figure><img src="image_29.png" /><figcaption>h</figcaption></figure>

<figure><img src="image_30.png" /></figure>

from heavier to lighter fractions (Supplementary Fig. 5e), while *Elavl4-v3* was too low in polysome fractions to yield meaningful results. This suggests that excess *Celf1* in our experiment is not sufficient to disturb *Elavl4* mRNA in polysomes.

Meanwhile, in the *Celf1* knockdown condition in N2a cells (Supplementary Figs. 4d, 6a), *Elavl4-v1&4* isoform presence on heavy polysomes increased relative to free, 40S-60S-80S

monosome, and light polysome fractions relative to control (Ctrl) shRNA (Fig. 5a). *Elavl4-v2* did not shift significantly (Fig. 5b) and again *Elavl4-v3* unreliably amplified. This aligns with increased translation of *Ren* via the *Elavl4-v1&4* 5' UTR in *Celf1* knockdown (Fig. 5c). In addition, we found increased translation of *Ren* via the *Elavl4-v3* 5' UTRs in *Celf1* knockdown (Fig. 5c). Silencing *Celf1* in N2a cells elevated steady-state levels of *Elavl4*
**Fig. 4** Celf1S translationally repressed *Elavl4-v3* and *-v4* via their 5' UTRs. **a** Relative mRNA levels determined by qRT-PCR of *pan-Elavl4* after either Control (Ctrl), *Celf1L*, or *Celf1S* overexpression (OE) in N2a cells (n = 6 transfections). Data represent the mean and SEM. Data normalized to *Gapdh* and then to Ctrl. Statistics: one-way ANOVA. n.s. = not significant. **b** (left) Western blot for Ctrl, *Celf1L*, or *Celf1S* OE in N2a cells (n = 6 transfections). Immunoblot of pan-Elavl4 (upper panel). Gapdh is the loading control (lower panel). (right) Densitometry quantification of Western blots. Data represent the mean, SD. Data normalized to Gapdh and then to Ctrl. Statistics: one-way ANOVA. *p* < 0.05. **c** Schematic of 3' UTR luciferase construct (left of bar graph). *Elavl4* isoforms' common 3' UTR downstream of *Luciferase (Luc)*, and cotransfected with either Ctrl OE/Celf1S into N2a cells. Translation assayed by measuring LUC RLU over *Luc* mRNA levels (n = 9 transfections). Data represent the mean and SEM. Statistics: Student's t-test. n.s. = not significant. **d** Schematic of 5' UTR Renilla construct (left of bar graph). N2a cells were cotransfected with Ctrl OE or *Celf1S* and *Elavl4-v2* 5' UTR upstream of *Renilla* (*Ren*). Same assay as in **c**; REN measured over *Ren* mRNA (n = 3 transfections). Data represent the mean and SEM. Statistics: Student's t-test. n.s. = not significant. **e, f** Deletion strategy of Celf1-regulatory regions in 5' UTRs of *Elavl4-v3* (**e**) and *Elavl4-v1&4* (**f**). *Elavl4-v3* or *Elavl4-v1&4* 5' UTRs Renilla constructs; same strategy as in **d**. Dashed lines represent deleted regions. Numbers indicate positions in UTR, 5' → 3'. n = 4 transfections for each experiment. bp = base pairs. Putative Celf1 regulatory sequence in *Elavl4* 5' UTRs indicated under their respective bar graphs. Data represent the mean and SEM. Statistics: Student's t-test and one-way ANOVA with Tukey's (**e, f**). *p* < 0.05; **p* < 0.01; ***p* < 0.001; ****p* < 0.0001; n.s. = not significant. **g, h** Celf1 RIP of *Elavl4-v3* (**g**) and *Elavl4-v1&4* (**h**) 5' UTRs. UTRs have WT, mutated (Mut), or deleted (Del) *Celf1* binding sequences. 18s was used as negative control. Data represent the mean and SEM; normalized to *Gapdh*. n = 3 RIPs from 3 transfections; Statistics: Two-way ANOVA with Tukey's. *p* < 0.05, **p* < 0.01, ***p* < 0.001.

a

F = Free M = 40S-60S-80S LP = Light polysomes HP = Heavy polysomes

Elavl4-v1&4

<figure><img src="image_31.png" /></figure>

<figure><img src="image_32.png" /></figure>

<figure><img src="image_33.png" /></figure>

<figure><img src="image_34.png" /></figure>

b

Elavl4-v2

<figure><img src="image_35.png" /></figure>

<figure><img src="image_36.png" /></figure>

<figure><img src="image_37.png" /></figure>

<figure><img src="image_38.png" /></figure>

d

WT

■

Celf1 HET

e

<figure><img src="image_39.png" /></figure>

<figure><img src="image_40.png" /></figure>

IB:

Celf1

IB:

Gapdh

IB:

E13

IB:

E13

IB:

E13

IB:

Gapdh

IB:

Gapdh

IB:

Gapdh

WT

■

Celf1 HET

f

■

Celf1 KO

c

O

Ctrl shRNA

△

Celf1 shRNA

<figure><img src="image_41.png" /></figure>

**Fig. 5** Celf1 deletion derepressed *Elavl4-v3* and *Elavl4-v4* translation. **a, b** Polysome fractionation of N2a cells transfected with either Ctrl shRNA or *Celf1* shRNA. Difference between free (F), 40S-60S-80S monosome (M), light polysome (LP), and heavy polysome (HP) fractions are represented for *Elavl4-v1&4* (**a**) and *Elavl4-v2* (**b**). This graph plots mean difference for each condition; a higher mean difference means a greater abundance in the minuend (n = 3 fractionations, each fractionation was from one 10 cm plate transfection). Statistics: Welch's t-test. ***p* < 0.001. *p* < 0.05, n.s. = not significant. **c**, N2a cells cotransfected with Ctrl shRNA or *Celf1* shRNA plasmids and *Elavl4-v3* or *Elavl4-v1&4-5' UTRs* cloned upstream of Renilla. Same assay as in Fig. 4. Multiple T-tests. *p* < 0.05, ***p* < 0.001. **d** qRT-PCR determined levels of different *Elavl4* isoforms in E13 WT and *Celf1* HET neocortices. Data represent the mean and SEM. N = 3 separate spins with 6 brains per spin. Data normalized to *Gapdh*. **e** (left) Western blot for WT and *Celf1* HET neocortices at E13 (n = 9 neocortices). (right) Densitometry quantification of Western blots. Data represent the mean and SD. Data normalized to *Gapdh* and then to WT. Statistics: Student's t-test. *p* < 0.05. **f** Neocortical VZ at E13 of WT (n = 6 animals) and *Celf1* KO (n = 2 animals). IHC for pan-Elavl4 (green). DAPI shown in blue. Scale bar: 40 µm.

mRNAs (Supplementary Fig. 6d, f, h), but there was no change in splicing nor in decay (Supplementary Figs. 5b, 6c, e, g). Collectively, our results indicate that a decrease in *Celf1* derepresses *Elavl4* mRNA translation.

To determine the effect of *Celf1* on *Elavl4* in vivo, we compared E13 WT mice to *Celf1* heterozygous (*Celf1* HET) littermates. *Celf1* HETs were analyzed because global deletion of

*Celf1* in mice reduces viability in a background-strain dependent manner³⁶,³⁷. While steady state levels of *Elavl4-v2* did not change, *Elavl4-v1&4* and *-v3* mRNA decreased in *Celf1* HET (Fig. 5d), contrasting the increase in *Elavl4* protein expression at E13 (Fig. 5e). Low amounts of mRNAs paralleled by a high level of corresponding protein expression has been reported in multiple systems³,¹⁰,¹¹,³⁸,³⁹. From the *Celf1* KO neocortices that we
<figure><img src="image_42.png" /><figcaption>**Fig. 6** *Elavl4-v3* OE lacking its 5' UTR in neural stem cells promotes LL identity. **a** Schematic: Neocortical IUE at E13 (left) with CAG-GFP and either Ctrl or *Elavl4* constructs were analyzed at E17 for colocalization between GFP and neocortical layer markers as indicated in the chart on the right. Nctx = neocortex, UL = upper layers, LL = lower layers, IC = intracortically projecting. **b** Representative images of E17 IUE Nctx, corpus callosum, and striatum. GFP + axons (green) are indicated by arrowheads. Scale bar: 200 µm. **c** Representative images of immunostained E17 IUE Nctx electroporated with GFP and either Ctrl (n = 3 animals), *Nestin:Elavl4-v3* (n = 4 animals), or *Nestin:Elavl4-v4* (n = 4 animals) OE. IHC for GFP (green) and cortical neuron identity markers (red) listed in **a**. Scale bar: 20 µm. **d** Quantification of colocalization between layer markers and GFP, normalized to Ctrl. Data represent the mean of sections and SEM. Statistics: one-way ANOVA with Tukey's. \*p < 0.05, \*\*p < 0.01, \*\*\*p < 0.001.</figcaption></figure>

managed to obtain, Elavl4 protein was enriched in the VZ when compared to WT at E13 (Fig. 5f), suggesting that Celf1 suppresses *Elavl4* mRNA translation in vivo during early neocortical neurogenesis.

*Elavl4-v3* promotes lower layer (LL) identity when not repressed in RG. To determine the consequence of premature *Elavl4* derepression in early neocortices, we cloned *Elavl4-v3*

and -v4 without their native 5' UTRs downstream of the RG-specific *Nestin* promoter. Ctrl, *Elavl4-v3*, and -v4 OE constructs were in *utero* electroporated (IUE) with CAG-GFP into E13 RG (Fig. 6a). GFP+ axons projected appropriately for the developmental stage in each condition (Fig. 6b). IHC of transfected neocortices at E17 revealed that when *Elavl4-v3* was not translationally repressed in RG, there was an increase in Tle4+ and an upward trend in Ctip2+ LL identity compared to Ctrl (Fig. 6c, d).
This finding is in line with our previous report that postnatal Tle4<sup>+</sup> neuronal identity is reduced by *Elavl4* deletion⁹. There was no effect on Cdp<sup>+</sup> upper layer (UL) neuronal identity by E17 when *Elavl4-v3* is OE in early RG (Fig. 6d). However, premature *Elavl4-v4* expression in RG resulted in a slight increase of Satb2<sup>+</sup> intracortically (IC) projecting neurons and a trending decrease of Ctip2<sup>+</sup> LL neurons when compared to Ctrl or *Nestin:Elavl4-v3*, respectively (Fig. 6c, d). These data indicate that premature expression of *Elavl4-v3* in RG results in a misbalance of neocortical glutamatergic neurons.

**Elavl4-v3 or -v4 with native 5' UTRs promote UL identity.** We next tested whether neuronal identities are disrupted when *Elavl4* isoforms are aberrantly expressed in neurons. *Elavl4-v3* and *-v4* were cloned under the *Cdk5r* promoter which is selectively expressed in early post-mitotic neurons⁴⁰. Ctrl, *Cdk5r:Elavl4-v3*, or *Cdk5r:Elavl4-v4* OE were IUE with CAG-GFP into RG at E13. Analysis at P0 revealed strikingly different results than the OE experiments in RG. When either *Elavl4-v3* or *-v4* are OE in early neurons, there was an increased number of UL Cdp<sup>+</sup> and IC Satb2<sup>+</sup> neuronal identities compared to Ctrl (Fig. 7a-c, g). These data suggest that the aberrant expression of *Elavl4* in post-mitotic neurons promotes UL identities.

To determine if these effects were due to the lack of the native 5' UTRs, we added isoform-specific 5' UTRs to the cell-type specific promoter-driven constructs used thus far. Constructs were IUE with CAG-GFP into E13 RG. Appropriate expression of each 5' UTR construct was confirmed using a specific FISH probe (Supplementary Fig. 7a). Surprisingly, while the UL Cdp+ identity was not affected by *Elavl4-v3* OE in RG (Figs. 6d and 7g), the addition of the native 5' UTR in RG increased this subpopulation (Fig. 7b, g; *Nestin:5' UTR-Elavl4-v3*). UL Cdp+ and IC Satb2+ identities were increased after OE of either *v-3* or *-v4* isoform in neurons. This effect was potentiated by the presence of 5' UTR in RG for *Elavl4-v3* (*Nestin:5' UTR-Elavl4-v3*). Meanwhile, the presence of the native *Elavl4-v4* 5' UTR potentiated the increase of UL Cdp+ and IC Satb2+ identities only when expressed in neurons (*Cdk5r:5' UTR-Elavl4-v4*), but not in RG (*Nestin:5' UTR-Elavl4-v4*). These results indicate that *Elavl4* isoforms promote UL and IC identity when the 5' UTR is present.

**LLs are sensitive to the presence of distinct *Elavl4* 5' UTRs.** We next turned our focus to the effect of 5' UTRs on LL identities. LL Ctip2+ identity decreased regardless of 5' UTR presence when *Elavl4-v4* was expressed in neurons (Fig. 7d, g; *Cdk5r:Elavl4-v4* and *Cdk5r:5' UTR-Elavl4-v4*). *Elavl4-v3* OE in RG decreased LL Ctip2+ identity only when its native 5' UTR was present (Fig. 7d, g; *Nestin:5' UTR-Elavl4-v3*). *Elavl4-v4* OE in neurons also decreased LL Ctip2+ identity (Fig. 7d, g; *Cdk5r:Elavl4-v4*). However, this isoform again proved to be especially sensitive to the presence of its 5' UTR. When the *Elavl4-v4* is OE in neurons with its 5' UTR, the LL Ctip2+ identity was normalized (Fig. 7d, g; *Cdk5r:5' UTR-Elavl4-v4*). Moreover, when the *Elavl4-v3* is introduced in the RG with its 5' UTR, the LL Ctip2+ identity increased (Fig. 7d, g; *Nestin:5' UTR-Elavl4-v3*).

The LL Tle4<sup>+</sup> subpopulation was likewise sensitive to isoform and 5' UTR presence. Even though both *Elavl4-v3* and *-v4* promoted Tle4<sup>+</sup> identity even when expressed in RG (Figs. 6d and 7g; *Nestin:Elavl4-v3* and *Nestin:Elavl4-v4*), *Elavl4-v3* decreased this identity in the presence of its 5' UTR in both RG and neurons (Fig. 7e, g; *Nestin:5' UTR-Elavl4-v3* and *Cdk5r:5' UTR-Elavl4-v3*). In contrast, by the presence of *Elavl4-v4* 5' UTR in either RG or neurons, Tle4<sup>+</sup> identities were unaffected (Fig. 7e, g; *Nestin:5' UTR-Elavl4-v4* or *Cdk5r:5' UTR-Elavl4-v4*). This indicates that

translational regulation of *Elavl4-v4* leads to the appropriate balance of the Tle4<sup>+</sup> identity (Fig. 7g). Surprisingly, the different constructs also appeared to have differential effects on migration (Fig. 7a, top row, Supplementary Fig. 9a, b).

GFP+ axons in all experimental conditions crossed the corpus callosum (Fig. 7f) and followed the white-matter tract above the VZ/SVZ⁴¹. However, differences were observed in the contralateral hemisphere of the IUE (Fig. 7f). We observed axons turning and growing into the opposite hemisphere when the 5' UTR was present in the OE constructs. This defect was particularly robust when the *Nestin* promoter was driving *Elavl4-v3* (*Nestin:5' UTR-Elavl4-v3*) and when *Cdk5r* was driving *Elavl4-v4* (*Cdk5r:5' UTR-Elavl4-v4*) (Fig. 7f). This finding suggests the excessive outgrowth when overexpression was restricted to the appropriate cell types (*Elavl4-v3* in RG and *Elavl4-v4* in neurons). When compared to Ctrl OE, fewer GFP+ axons reached the striatum in all conditions except when *Elavl4-v4* OE was regulated by its 5' UTR (Fig. 7f; *Nestin:5' UTR-Elavl4-v4* and *Cdk5r:5' UTR-Elavl4-v4*). It should also be noted that GFP+ cells from FACS-sorted transgenic lines co-express *Cdp* and *Tle4* mRNA (Supplementary Fig. 8c), supporting the importance of translational control and corroborating recent findings that RG are transcriptionally primed¹⁵. These data corroborate the idea that translational regulation of distinct *Elavl4* isoforms in early post-mitotic neurons strongly impacts the development of glutamatergic neurons to the risk of neuroanatomical pathologies.

**Celf1 regulates the balance of glutamatergic neurons.** Since Celf1 translationally repressed *Elavl4* isoforms through their 5' UTRs, we tested if mutating Celf1 binding sites in *Elavl4* 5' UTRs abolished the defects observed in Fig. 7. We IUE CAG-GFP and either Ctrl OE, 5' *UTR-Elavl4-v3*, 5' *UTR-Elavl4-v4*, or Celf1 binding site (BS)-mutated constructs (*Celf1BSmut:5' UTR-Elavl4-v3* and *Celf1BSmut:5'UTR-Elavl4-v4*) at E13 and analyzed transfected neocortices at P0. The mutation of Celf1 BS in 5' UTRs of *Elavl4* isoforms abolished the effects on UL Cdp+ and LL Ctip2+ identities induced by native 5' UTRs (Supplementary Fig. 9c, d).

To assess Celf1's possible role on identities, we first analyzed Celf1 HET and KO neocortices at E13 when *Elavl4* protein is increased. Celf1 deletion increased Cdp+ and decreased Ctip2+ identity (Fig. 8a, Supplementary Fig. 7c-f), mirroring *Elavl4* OE phenotypes. To interrogate the role of Celf1 cell-autonomously, we performed IUE electroporation of Ctrl shRNA+CAG-GFP or Celf1 shRNA+CAG-RFP into separate embryos of the same litter and performed in vitro analysis using a previously reported approach¹³ (Supplementary Fig. 8a). After three days in vitro (DIV) we found that Cdp+ identity increased (Fig. 8b) and Ctip2+ identity decreased in Celf1 knockdown (Supplementary Fig. 8b). These in vitro findings substantiated the observation in the Celf1 HET and KO neocortices (Fig. 8a, Supplementary Fig. 7c-f). Notably, the identity changes observed in Celf1-depleted neocortices and neurons is comparable to *Elavl4-v3* and *-v4* OE without their 5' UTRs (Fig. 7).

The possible Celf1 HET effects on identities are independent of the cell cycle, Tbr2+ progenitors, or neurogenesis (Supplementary Fig. 7b). However, we found a higher retention of RG Pax6+ progenitors in Celf1 HET neocortex (Supplementary Fig. 7b). Since *Elavl4* protein is increased in Celf1 HETs, the Pax6+ retention corroborates previous findings where *Elavl4* decrease resulted in fewer progenitors²⁸,²⁹. Overall, these data suggest that Celf1 acts on the balance of glutamatergic neurons through *Elavl4*.

**Celf1 regulates development of the neocortex via *Elavl4*.** The data thus far suggest that Celf1 reduction leads to a selective neuronal identity misbalance, possibly induced by *Elavl4* increase
<figure><img src="image_43.png" /><figcaption>**Fig. 7** Native 5' UTRs of *Elavl4* isoforms dictate neocortical glutamatergic identities. **a** (top) PO representative confocal images of E13 IUE neocortices electroporated with CAG-GFP and either Ctrl (n = 6 animals), Cdk5r:Elavl4-v3 (n = 5 animals), Cdk5r:5' UTR-Elavl4-v3 (n = 4 animals), Nestin:5' UTR-Elavl4-v3 (n = 3 animals), Cdk5r:Elavl4-v4 (n = 6 animals), Cdk5r:5' UTR-Elavl4-v4 (n = 3 animals), Nestin:5' UTR-Elavl4-v4 (n = 3 animals) OE plasmids. IHC for GFP (green), and cortical layer identity markers Cdp, Satb2, Ctip2, or Tle4 (all red). Scale bar: 20 µm. DAPI shown in blue. **b-e** Quantification of colocalization between layer markers (arrowheads in A) and GFP normalized to controls. Data represent the mean of sections and SEM. Statistics: one-way ANOVA with a post-hoc Tukey's. \*p < 0.05, \*\*p < 0.01, \*\*\*p < 0.001, \*\*\*\*p < 0.0001. All significance stars that are located directly above bars are in reference to the Ctrl, otherwise indicated with line above. **f** PO representative confocal images of E13 IUE neocortices electroporated with GFP and either Ctrl (n = 6 animals), Cdk5r:Elavl4-v3 (n = 5 animals), Cdk5r:5' UTR-Elavl4-v3 (n = 4 animals), Nestin:5' UTR-Elavl4-v3 (n = 3 animals), Cdk5r:Elavl4-v4 (n = 6 animals), Cdk5r:5' UTR-Elavl4-v4 (n = 3 animals), Nestin:5' UTR-Elavl4-v4 (n = 3 animals) OE plasmids IHC for GFP (green). Arrows indicate GFP+ axons passing through corpus callosum and striatum ipsilateral to IUE site. Arrowheads indicate GFP+ axons growing into neocortex contralateral to IUE site. DAPI shown in blue. Scale bar: 200 µm. **g** Summary table of results in Fig. 6 (E17; gray) and 7a-7e (PO; colored). Increased colocalization is indicated with a green up arrow and decreased colocalization is indicated with a red down arrow.</figcaption></figure>
<figure><img src="image_44.png" /><figcaption>b</figcaption></figure>

<figure><img src="image_45.png" /><figcaption>c</figcaption></figure>

<figure><img src="image_46.png" /><figcaption>d</figcaption></figure>

<figure><img src="image_47.png" /><figcaption>e</figcaption></figure>

<figure><img src="image_48.png" /><figcaption>f</figcaption></figure>

<figure><img src="image_49.png" /></figure>

<figure><img src="image_50.png" /><figcaption>g</figcaption></figure>

<figure><img src="image_51.png" /></figure>

<figure><img src="image_52.png" /><figcaption>h</figcaption></figure>

<figure><img src="image_53.png" /></figure>

<figure><img src="image_54.png" /></figure>

**Fig. 8 Celf1 regulates development of glutamatergic neurons in developing neocortices through Elavl4.** a Neocortical E13 VZ of WT (n = 6 animals) and Celf1 KO (n = 2 animals). IHC for Cdp (red). DAPI shown in blue. Scale bar: 100 µm. **b** (left) Representative IHC of three days in vitro primary neocortical neuronal cultures. Transfected neurons are expressing GFP (Ctrl shRNA) or RFP (Celf1 shRNA) and neuronal identity marker Cdp (white). Scale bar: 40 µm. (right) Quantification of Cdp colocalization. Data represent the mean and SEM. N = 1 litter (14 embryos). Statistics: Students t-test. ***p < 0.001.
**c** Schematic: Celf1 translationally represses Elavl4-v3, which decreases Cdp+ identities. Once Celf1 is silenced, Elavl4-v3 translation is derepressed, increasing Cdp+ identities. **d** PO representative confocal images of E13 IUE neocortices electroporated with CAG-GFP and either Ctrl shRNA (n = 5 animals) or Celf1 shRNA (n = 4 animals) at E13. Ventricular zone/subventricular zone (VZ/SVZ) shown underneath. GFP is green. DAPI shown in blue. Scale bar: 100 µm. UL = upper layers, LL = lower layers. **e** (left) PO representative confocal images of transfected neocortices in **d**, immunostained for cortical neuronal identity markers Cdp and Ctip2 in CP. Arrowheads indicate colocalization. Scale bar: 40 µm. (right) Quantification of Cdp and Ctip2 colocalization. Data represent the mean and SEM. Statistics: Student's t-test. ***p < 0.001. **f** (left) Representative confocal images of transfected neocortices in **e**, immunostained for cortical neuronal identity marker Cdp in VZ/SVZ. Arrowheads indicate colocalization. Scale bar: 40 µm. (right) Quantification of Cdp colocalization. Data represent the mean of sections and SEM. Statistics: Student's t-test. ***p < 0.001. **g** (left) PO representative confocal images of E13 IUE neocortices electroporated with CAG-GFP and either Ctrl shRNA (n = 5 animals), Celf1 shRNA and Celf1 OE (n = 3 animals), or Celf1 shRNA and Elavl4-v3 shRNA (n = 3 animals). IHC for GFP (green) and layer marker Cdp (red). Arrowheads indicate colocalization. Scale bar: 40 µm. UL = upper layer. (right) Quantification of Cdp colocalization. Data represent the mean of sections and SEM. Statistics: one-way ANOVA with Tukey's. *p* < 0.05, ***p* < 0.001, n.s. = not significant. **h** (left) PO representative confocal images of E13 IUE neocortices electroporated with GFP and either Ctrl (n = 4 animals), Celf1 OE (n = 3 animals), or Celf1 OE and 5' UTR-Elavl4-v3 (n = 3 animals) OE plasmids IHC for GFP (green), and cortical layer identity markers Cdp (red). Arrowheads indicate colocalization. DAPI shown in blue. Scale bar: 40 µm. (right) Quantification of Cdp colocalization. Data represent the mean and SEM. Statistics: one-way ANOVA with Tukey's. ***p* < 0.001. **i, j** Representative confocal images of the corpus callosum, contralateral neocortex, and striatum of PO neocortices that underwent IUE with GFP and either Ctrl shRNA (n = 5 animals) or Celf1 shRNA (n = 4 animals) (**i**) or either Ctrl OE (n = 3 animals), Celf1 OE (n = 3 animals), or Celf1 OE and 5' UTR-Elavl4-v3 (**j**). Scale bar is 200 µm.
(Fig. 8c). We next IUE E13 RG with CAG-GFP and either Ctrl or Celf1 shRNA and analyzed neocortices at P0. Celf1 silencing increased the number of transfected cells within the VZ/SVZ region (Fig. 8d, bottom), possibly due to a migration defect. Upon examination of neurons that reached the CP, there was no difference in UL Cdp+ identity between Ctrl and Celf1 shRNA, while LL Ctip2+ identity in the CP decreased (Fig. 8e). However, when we analyzed GFP+ cells in the VZ, UL Cdp+ identity increased when Celf1 was silenced (Fig. 8f). Co-expression of Celf1 shRNA and Celf1 OE decreased the UL Cdp+ identity (Fig. 8g), suggesting the specificity of Celf1 shRNAs.

To test if increased UL Cdp+ identity after Celf1 silencing is due to *Elavl4-v3* (Fig. 8c), we obtained specific *Elavl4-v3* shRNAs (Supplementary Fig. 4f). *Elavl4-v3* shRNA decreased UL Cdp+ identities and this was rescued by co-IUE with human *ELAVL4-V3* OE that is resistant to the mouse-specific shRNA (Supplementary Fig. 8f). Remarkably, silencing *Elavl4-v3* in *Celf1*-depleted neurons decreased UL Cdp+ identities (Fig. 8g) initially induced by *Celf1* silencing (Fig. 8a, b, f). Taken together, these data suggest that *Elavl4-v3* is downstream to *Celf1* in the regulation of neocortical neuronal identities and their development.

To further investigate the pathway, we IUE E13 neocortices with CAG-GFP and either Ctrl or Celf1s OE and analyzed the transfected neocortices at P0. Celf1s OE decreased Cdp+ identity (Fig. 8h). The Celf1s OE effect on the UL Cdp+ identity was reversed when 5' UTR-Elavl4-v3 (Supplementary Fig. 4e) was co-IUE (Fig. 8h). Interestingly, co-IUE of 5' UTR-Elavl4-v3 and Celf1s OE rescued UL Cdp+ identity up to control levels. This rescue further indicates a suppressive upstream role of Celf1s on *Elavl4* phenotypes.

When Celf1 was silenced in vivo, Tle4+ identity decreased (Supplementary Fig. 8e). Celf1s OE in Celf1 shRNA cells partially rescued Tle4+ identity and silencing of *Elavl4-v3* in Celf1 shRNA cells reversed the Tle4+ identities (Supplementary Fig. 8e), further corroborating that *Elavl4-v3* is downstream to Celf1. While Celf1 silencing can promote Tle4+ identity in vitro, these data indicate that Tle4+ identity is indirectly influenced by Celf1 in vivo, possibly non cell-autonomously.

When Celf1 was silenced, GFP+ axons were found to cross the CC and accumulate above the VZ/SVZ. Meanwhile, a considerably reduced number of GFP+ axons reached the striatum (Fig. 8i). This axonal distribution matched the increase in Cdp+ and decrease in Ctip2+ neurons projecting through callosum or subcortically, respectively (Fig. 8b, e, f, Supplementary Fig. 8b). Remarkably, Celf1s OE induced an aberrant GFP+ tract to terminate ipsilaterally to the IUE instead of crossing the corpus callosum (Fig. 8j). The abnormal ipsilateral tract was rescued when Celf1s was co-OE with 5' UTR-Elavl4-v3 (Fig. 8j). All in all, these data indicate that Celf1 regulates development of glutamatergic neurons through Elavl4. Given that Celf1 OE at E13 affected connectivity, we IUE E16.5 neocortices with CAG-RFP+ Ctrl OE in one hemisphere and CAG-GFP+ Celf1s OE in the contralateral hemisphere and analyzed neocortices at P7. Remarkably, the hemisphere with Celf1s OE had ectopic GFP+ subcortical projections descending through the striatum and reaching the internal capsule (Supplementary Fig. 9e). In the contralateral RFP+ Ctrl hemisphere no axons were observed in the thalamus or the internal capsule. These data further suggest that Celf1s OE can disrupt callosal development and induce ectopic subcortical projections from IC-projecting neurons.

**CELF1 and ELAVL4 are expressed in human ventral RG.** Given the striking impact of Celf1-mediated regulation on identity and connectivity of glutamatergic neurons, we examined if

the 5' UTRs of mouse *Elavl4* isoforms were conserved in human. At least fifteen human *ELAVL4* isoforms exist (https://genome.ucsc.edu/), suggesting an increase in the complexity of translational control throughout evolution. The 5' UTR of mouse *-v1&4* was highly conserved across species, including the region of Celf1 regulatory site (Fig. 9a, b). Meanwhile, the 5' UTRs of mouse *-v3* and *Elavl3* (HuC) were less conserved (Fig. 9b), suggesting evolutionary divergence of the ELAVL4 protein.

We next wanted to determine if the CELF1 and ELAVL4 protein expression patterns in human developing neocortices resemble those we found in mice. We found conserved changes in ELAVL4 and CELF1 expression in the VZ and SVZ as human neocortical neurogenesis proceeds (Fig. 9c). ELAVL4 protein was enriched in the ventral RG of the VZ at later stages of neurogenesis, 17 and 20 post-conceptual weeks (PCW), which was decreased at 10 and 13 PCW. ELAVL4 mRNA has been reported to be expressed ventral RG and as a novel human intermediate progenitor marker in human developing neocortices between 14 and 16 PCW<sup>42,43</sup>. In agreement with this, we found ELAVL4-protein expression in the inner subventricular zone (ISVZ) and outer subventricular zone (OSVZ) of all analyzed stages, but also scattered across the periventricular fiber rich zone (PVFRZ) positioned above VZ and ISVZ<sup>44,45</sup> (Fig. 9c). The PVFRZ has been suggested to contain contralateral callosal fibers in human<sup>44</sup>, which is in agreement with murine contralateral callosal fibers positioned above VZ/SVZ that we observed in Figs. 7f and 8i, j. In addition, we found ELAVL4-protein to be expressed in ventral RG after 17 PCW, at later stages of human neurogenesis when CELF1 expression decreases. We found CELF1 to be enriched in ventral RG in the VZ during 10 and 13 PCW and decreased by 17 and 20 PCW (Fig. 9c). Thus, CELF1 had a contrasting temporal expression pattern to ELAVL4 in the VZ (Fig. 9c), similar to what we observed in mouse neocortex (Fig. 3d). In addition to confirming the relevance of our findings for human neocortical development, our findings identify markers of early and late human ventral RG.

*ELAVL4* mRNA expression was present in the VZ at all stages (Fig. 9d), including early on when ELAVL4 protein is low (Fig. 9c). The discrepancy between *ELAVL4* mRNA and ELAVL4 protein expression is similar to what we observe in mouse early RG (Figs. 1–3). These findings demonstrate conserved expression patterns and suggest at least partially conserved ELAVL4/CELF1 regulatory mechanism between early and late stages of neurogenesis in both mouse and human.

**Multiple mRNAs associated with neurodevelopmental disorders are translationally regulated in developing neocortices.** Given the conservation of the expression patterns of ELAVL4 and CELF1 in ventral RG during neocortical development in humans, we sought to explore the disease implications of the translational regulatory mechanisms. Initially, evolutionary constraint scores were computed by the Exome Aggregation Consortium (ExAC) as a measure of intolerance to damaging mutations<sup>46</sup>. We found that RBPs are highly constrained, including the members of the ELAVL family (Supplementary Fig. 10). Next, we performed enrichment analyses for published mRNA targets of the neuronal ELAVLs (ELAVL1, CELF1, and CELF4). Various degrees of enrichment were found in risk genes for neurodevelopmental disorders (NDDs) as defined by the DDG2P database, intellectual disability (ID; in OMIM), epileptic encephalopathy (EE; in OMIM) and ASD<sup>47</sup> (Fig. 9e, Supplementary Fig. 10a). We also found several translationally regulated mRNAs to be targets of an RBP associated with ID and ASD, Fragile X mental retardation protein (FMRP)<sup>48,49</sup> (Supplementary Data 3). The human orthologues of the translationally-regulated mRNAs during
<figure><img src="image_55.png" /><figcaption>b</figcaption></figure>

<figure><img src="image_56.png" /><figcaption>c</figcaption></figure>

<figure><img src="image_57.png" /><figcaption>d</figcaption></figure>

<figure><img src="image_58.png" /><figcaption>e</figcaption></figure>

<figure><img src="image_59.png" /><figcaption>f</figcaption></figure>

<figure><img src="image_60.png" /></figure>

neocortical development (Fig. 1, Supplementary Data 1) were also found to be constrained (Supplementary Fig. 10b).

Given that these genes are under strong selective pressure, we conjectured that they might be essential for physiological neurodevelopment and vulnerable targets in NDDs. To explore this hypothesis, we cross referenced the list of differentially translated mRNAs at E16 with risk genes for NDDs in the

DDG2P database; 61 of the 850 unique human orthologs derepressed at E16 are NDD genes (57 OMIM disease genes); 62 of the 1,070 genes that are repressed at E16 were NDD genes (59 in OMIM) (Supplementary Data 2). The mouse orthologs of genes associated with EE (e.g., *GRIN1* and *GABRA1*), ID (e.g., *MEF2C*), and ASD (e.g., *CHD8*, *MYT1L*), often in co-morbidity, were amongst mRNAs showing robust change in their polysome
**Fig. 9 ELAVL4 and CELF1 have opposing expression in developing human neocortical VZ.** a Alignment of mouse *Elavl4* 5' UTRs (top). Conservation of UTRs shown below. b The plot shows the probability (y-axis) that each nucleotide in the 5' UTR (x-axis) belongs to a conserved element for *Elavl4-v1/4* (pink), *Elavl4-v3* (light blue) and *Elavl3* (green). Each dot represents a base pair. Celf1 regulatory binding sites (BS)-are indicated by the black bar/shaded area. nt, nucleotide. c Representative confocal images of human ventricular zone (VZ), inner subventricular zone (ISVZ), periventricular fiber rich zone (PVFRZ), and outer subventricular zone (OSVZ) at 10.5, 13, 15, 17, and 20 PCW in the frontal lobe. IHC for PAX6 (green) and either CELF1 (red) or ELAVL4 (red), indicated above image. DAPI shown in blue. d Representative images of in situ hybridization (ISH) for *Elavl4* on human VZ, ISVZ, PVFRZ, OSVZ at 11, 13, 15, 17, and 20 PCW in the frontal lobe. Scale bar: 100 µm. e The heatmap shows the degree of enrichment for mRNA targets of neuronal ELAVL5s (nELAVL5s), ELAVL1, CELF1 and CELF4 across neurodevelopmental disorders genes defined by the Developmental Disorders Genotype-Phenotype Database (NDD_DDG2P), epileptic encephalopathy genes in OMIM (EE_OMIM), intellectual disability genes in OMIM (ID_OMIM) and ASD risk genes (PMID: 31981491). Color coding indicates -log10(p-value) after 10,000 permutations. Number of overlapping genes is indicated in each cell. f Examples of top translationally regulated genes implicated with ASD, ID, and EE. Note that some of these genes can manifest with co-morbidities.

association (Fig. 9f; Supplementary Data 2). Collectively, these findings indicate that genes subjected to translational control during neocorticogenesis, including RBPs, are under strong selective pressure and might contribute to the vulnerability leading to NDDs.

Discussion

The data contained herein demonstrate that one mechanism to achieve the incredible diversity of glutamatergic neurons is translational control by specific RBPs acting on alternative 5' UTRs of corresponding mRNAs. We demonstrate that identity, subpopulation balance, and connectivity of neocortical glutamatergic neurons are dependent on translational control of *Elavl4* isoforms at the 5' UTRs by Celf1 (Fig. 10). Previously, selective isoform-specific protein synthesis (mRNA translation) and the role of timed releases of translational breaks in any system was poorly understood.

To our knowledge, our study is the first genome-wide report of isoforms that are translationally repressed and derepressed during neocortical neurogenesis; many of these translationally-regulated genes turn out to be implicated in the etiology of NDDs. Our findings reveal that a robust set of mRNAs and their isoforms are poised to be translated rapidly at an appropriate time point. This sets up not only a highly plastic environment, but also a vulnerable one during cellular development.

The pattern of timed translational repression/derepression we revealed applies to *Elavl4* isoforms. The diverse effects of translational control on *Elavl4* isoforms within early and late RG and neurons is due to their alternative 5' UTRs. The presence of specific 5' UTRs in *Elavl4* isoforms reversed LL subpopulation-specific glutamatergic phenotypes that were seen with the OE of unregulated *Elavl4*. The differential identity acquisition seen in our cell subtype-specific experiments may further rely on available cell-specific transcripts that RBPs selectively act upon; as 5' UTR-driven expression favors expression in one subtype identity over another. *Elavl4* isoforms also teach us that the significance of an mRNA cannot be guessed by its abundance. While *Elavl4-v3* and *Elavl4-v1&v4* are heavily regulated by Celf1, *Elavl4-v2* is less influenced by Celf1. Yet, *Elavl4-v2* is the isoform with the highest abundance. Besides differences in their 5' UTRs and expression patterns, *Elavl4-v3* and -*v4* exhibit important differences in their “linker” regions positioned between the RNA binding domains, which was found to determine its association with polysomes⁵⁰,⁵¹. Therefore, both the UTR and CDS may account for the differential phenotypes that are reflected in the isoform-specific experiments of our study.

Our screen led us to investigate another RBP, Celf1, which acts as an upstream regulator of *Elavl4* at its alternative 5' UTRs. We found that Celf1 suppressed translation of distinct *Elavl4* isoforms by binding their alternative 5' UTRs. Mutation of Celf1 binding sites also abolished the subpopulation phenotype observed upon OE of *Elavl4* isoforms with their native 5' UTRs.

Given the number of common binding motifs associated with other RBPs in our screen, Celf1 is likely only one example of an RBP that is an isoform-specific translation regulator and that has broad impacts across cell types during neocortical neurogenesis. Nevertheless, our findings underscore the potential diversity in expression that alternative splicing and precise dynamics of selective isoform translation provides.

We also report for the first time that Celf1 silencing itself disrupted the identity, subpopulation balance, and connectivity of glutamatergic neurons. This phenotype can be reversed by silencing of the translationally derepressed *Elavl4* isoform, *Elavl4-v3*. The OE of Celf1 also leads to disrupted development of glutamatergic neurons that can be corrected by overexpressing *Elavl4-v3* bearing its native 5' UTR. These findings reaffirm that *Elavl4* is downstream to Celf1 and demonstrate that both an upstream regulator and its mRNA target with native regulatory sites work together to give us a balanced neocortex. It is reasonable to expect that other regulators like miRNAs or other RBPs will also compete with Celf1 to not only bind, but also to regulate *Elavl4* expression²⁰,⁵².

There is immense potential for increased understanding of NDDs by studying the mechanisms that control translation during development. Our study has revealed the magnitude of translationally repressed and derepressed mRNAs and their associated isoforms. Disrupted functions and mutations of RBPs are linked with several NDDs associated with abnormal neocortical development. This is particularly interesting given the various neuroanatomical aberrations caused by Celf1 and *Elavl4* disruptions. Indeed, prenatal deletion of *Elavl4* has been associated with disrupted development of neocortical circuits, misbalance of glutamatergic identities, seizure susceptibility, and autism-like repetitive behaviors⁹,¹⁹,⁵¹,⁵³,⁵⁴. Another paradigmatic example of an RBP associated with NDD is Fragile X syndrome, which is caused by loss of the FMRP and the translational control it executes⁴⁹,⁵⁵⁻⁵⁷. FMRP⁴⁹,⁵⁸, ELAVL1¹³, neuronal ELAV, Celf1 and Celf4⁵⁹ (Fig. 9e), Celf6⁶⁰ and CPEB4⁶¹ target mRNAs are encoded by ASD risk genes⁶²⁻⁶⁴. We found that both CELF1 and ELAVL4 are expressed in developing prenatal human neocortices suggesting their roles and vulnerability to disruption during development. RBPs can thus have diverse effects in different NDDs by causing global misregulation of their mRNAs. This suggests a previously overlooked mechanism contributing to the manifestation of comorbid conditions. Future studies are needed to elucidate the complete picture of collaborative and competitive differences in the control of isoform-specific 5' UTRs, which we can now appreciate as critical in fully understanding how neocortical RG and neurons develop.

Methods

**Animals.** Experiments involving animals were carried out in accordance with Rutgers University Medical School's Institutional Animal Care and Use Committee. For timed pregnancies, adult pregnant female CD-1 mice were purchased from
<figure><img src="image_61.png" /><figcaption>**Fig. 10 Summary graphic.** Development of glutamatergic neurons is regulated by Celf1 translational repression of *Elavl4* through its 5' UTRs.</figcaption></figure>

Charles River Laboratories. The day of vaginal plug discovery was considered embryonic day 0.5 (E0.5). The *Elavl4-GFP* transgenic line was obtained from GENSAT (www.gensat.org) and previously reported⁵. *Elavl4* KO, Celf1 KO, *Nestin-GFP*, and *Tbr2-GFP* transgenic mice were described previously⁹,²⁸,³⁰⁻³³,³⁶. For all embryonic experiments, mice of either gender were used and gender was not determined. Gender was also not determined in human fetal tissue.

**Postmortem human brain tissue.** All experiments on human tissue were carried out in accordance to protocols approved by the Institutional Review Board (IRB) of the Ethical Committee of the School of Medicine, University of Zagreb and were carried out in accordance with the Declaration of Helsinki 2000. Fetal brain specimens are part of the Zagreb Brain Collection, received from clinical hospitals affiliated with the School of Medicine, University of Zagreb⁴⁵. Tissue was fixed with a post-mortem delay of less than 6 h. Post-conceptual weeks (PCW) were determined according to the clinical pregnancy records and crown-rump length.

**Polysome fractionation.** Polysome fractionation of neocortices was performed as previously described¹⁰. Briefly, embryonic neocortices were dissected in ice-cold HBSS media and immediately frozen on dry ice (n = 6 brains per spin). Tissue was homogenized in polysome extraction buffer (PEB) and the input was layered on top of a 10–50% sucrose density gradient in thin wall polypropylene centrifuge tubes (Cat. no. 347357), with surplus input being stored for further processing. Gradiations were spun on Sorvall Discovery M120SE in Sorvall S55S rotor with swinging buckets (#18507) at 39,000 RPM for 50 min. Fractions were collected using Fluorinert FC-40 (Cat. no. F9755) and a Brandel fraction collector (no. 621140007), with UV absorbance recorded at 254 nm (Brandel UA-6). Samples were stored at –80 °C until further processing.

N2a cell cultures were grown in 100 mm plates to 75-85% confluency and transfected with 15 µg of plasmid using Lipofectamine 2000 in a 2:1 ratio

(Lipofectamine:DNA). A 15 min cyclohexamide treatment (0.1 mg/ml) was conducted at the end of the 48-hour transfection. Polysome fractionation of frozen neocortices and cell pellets started with lysis in PEB, freshly supplemented with EDTA-free protease inhibitor (Santa Cruz Biotechnology; no. sc-29131), RNaseOut (Invitrogen; no. 100000840), 20 mM DTT (Invitrogen; no. NP0009), and 0.1 mg/ mL cyclohexamide (Santa Cruz Biotechnology; no. sc-3508A). Lysis was conducted on ice for 10 minutes, pipetting up and down every two minutes. Lysate was spun for 10 min at 4 °C and 5 K RPM; the resultant supernatant was spun for 5 min at 4 °C and 14k RPM. Samples were measured using a Nanodrop under the A₂₉₉ wavelength to normalize loading to total protein content. Sample was loaded onto a 10-50% sucrose gradient. Gradiations were prepared the evening before in 2 ml polyallomer tubes (Beckman Coulter; no. 347357) and left to stabilize to a continuous gradient overnight at 4 °C. Gradiations were spun for 120 min using a Thermo Fisher Sorvall MX 120+ micro-ultracentrifuge and the Sorvall S-55-S swinging bucket rotor at 4 °C. Fractions were collected until the presence of Fluorinert FC-40 (Sigma; no. F9755) was detected. Samples were stored at –80 °C until further processing.

RNA was isolated using TRIzol reagent and mostly according to the manufacturer's protocol briefly: the nucleic acid from each fraction was isolated using a 3:1 ratio of TRIzol:sample. The mixture was washed with chloroform twice to isolate and purify the aqueous phase. Total RNA was precipitated with one volume of isopropanol/glycogen (Thermo Fisher; no. R0551). The pellet was washed in 75% ethanol and the resulting pellet was dried and resolubilized with Molecular Biology Grade water (Corning; no. 46-000-CV) supplemented with RNaseOut for ~30 min, vortexing frequently. DNase treatment was conducted according to the manufacturer's protocol (Invitrogen; no. AM1907); DNase was inactivated with phenol:chloroform (Fisher Scientific; no. BP1754I). RNA was further purified by chloroform extraction and precipitated with 3 M sodium acetate (G-Biosciences, no. R010). Samples were stored at –80 °C in 100% ethanol until further processing.
qRTs were conducted as described in the relevant methods subsection. Polysome inputs were diluted 1:25 for *Elav4-v1/4* and *Elav4-v2* assays and normalized to *Gapdh*. Fractions were diluted 1:5 for *Elav4-v1/4* and *Elav4-v2*. *Elav4-v3* assays were run with undiluted sample. Non-template controls were used to determine background levels of amplification and verify fraction amplification. Polysome fractionation was recorded with analog and digital traces recorded by UV-vis at 260 nm and rRNA 18 S was tested for in each fraction. CELF1 short OE and KD were confirmed by qRT using a *Celf1-pan* TaqMan probe. Ct values were used in comparison analysis. We first applied the closure operator so that the components added to 1.0⁶⁵. Differences in polysome fractions (e.g., HP - F) were tested between conditions using Welch's *t*-test. Nominal significance levels are indicated in the figures.

**RNA-seq, GO, and KEGG analysis.** RNA-seq was performed at RUCDR Infinite Biologics™ using Illumina technology and analyzed as described previously¹⁰. Briefly, RNA was isolated from input and polysome fractions by using TRIzol LS (Life Technologies; no. 10296028) following the manufacturer's protocol. Next, equal volumes of RNA extracted from fractions 4–7 were pooled together for analysis of “40S–60S–80S-associated cytoplasmic RNA” and equal volumes of RNA extracted from fractions 9–12 were pooled together for “polysome-associated cytoplasmic RNA”. Three biological replicates were analyzed for input, 40S–60S–80S, and polysome RNA at E13 and E16 in WT (18 total samples). Sequencing libraries were prepared by using the Illumina TruSeq RNA Sample Preparation Kit v2 according to the manufacturer's protocol. Libraries were quantified by using the Library Quantification Kit Illumina/Universal (KAPA Biosystems) and then diluted and symmetrically pooled. We performed 2 × 75-bp paired-end sequencing using the Illumina HiSeq2500. Sequencing data have been deposited in the Gene Expression Omnibus (GEO) database, www.ncbi.nlm.nih.gov/geo (accession no. GSE77674).

Results were aligned with the mm10 mouse genome found in the University of California Santa Cruz transcript map (Illumina iGenomes) using TopHat (v 2.1.1), and comparisons between groups were made in Cufflinks (v 2.2.1) and cummerbund (v 2.22.0)⁶⁶. Significant differences were judged using a 5% false discovery rate (FDR). The detection of up- or down-regulated mRNAs is more sensitive to lower expression levels in the polysome fraction than in total cellular RNA, so differences in polysome mRNAs were assessed whether these mRNAs were detected in total or not. Relatively few translationally derepressed mRNAs were below the detection threshold in total RNA (97 out of 971 mRNAs). Lists of regulated genes were assessed for enrichment of functional groups or pathways using g-Profiler²⁶. Transcript isoforms that were significantly increased or decreased (from E13 to E16) in polysomes, but unchanged in total RNA were used to prepare unique Ensembl Gene ID lists (using the biomaRt package in R)⁶⁷ for upload to the g-Profiler site. Significantly enriched gene ontology groups, as represented by the best term per parent ontology chosen by p-value, are shown (Supplementary Fig. 1). Similarly, the best terms from a KEGG pathway search are included.

**Identifying common RNA-binding motifs.** Transcript isoforms that were significantly up- or down-regulated in E13 vs. E16 polysomes and unregulated in total RNA were selected and used to download FASTA-formatted cDNA sequences using biomaRt. For comparison, all known cDNA FASTA sequences were downloaded. Motifs enriched in either the up- or down-regulated polysome lists compared with all cDNAs were identified using Homer v 4.9⁶⁸, with the “-rna” flag.

**RNA-binding proteins.** RNA-binding proteins (RBPs) were selected from the census curated by Gerstberger et al.⁶⁹ filtered for RBPs with mRNA as consensus RNA target and “established” as supporting evidence as (205 unique genes).

**Primary antibodies.** The following primary antibodies and dilutions were used: mouse anti-Elav4 (anti-HuD; Santa Cruz Biotechnology (SCB), 1:1,000, sc-28299), rabbit anti-Pax6 (Biolegend previously Covance, 1:250, 901301 previously PRB-278), rabbit anti-Tbr2 (Abcam, 1:250, ab23345), rabbit anti-Cdp (SCB, 1:250, sc13024), rat anti-Bcl11b (SCB, 1:250, sc-98514), mouse anti-Celf1 (SCB, 1:100, sc20003), rat anti-Ctip2 (Abcam, 1:250, ab18465), mouse anti-Tle4 (SCB, 1:250, sc-365406), chicken anti-GFP (Aves, 1:1000, GFP-1020), mouse anti-Gapdh (Millipore Sigma, 1:2000, MAB374), mouse anti-Satb2 (Abcam, 1:250, ab51502), rabbit anti-Nestin (Sigma Aldrich, 1:100, N5413), chicken anti-N-Cadherin (Takara, 1:250, M110 Clone NCD-2), rabbit anti-NeuN (Millipore Sigma, 1:500, ABN78), rabbit anti-pH3 (Millipore Sigma, 1:1000, 06-570), rat anti-Brdu/Cldu (Novus, NB500-169), mouse anti-Brdu/IdU (BD Biosciences, 1:100, B44), and rabbit anti-Renilla (Thermo Fisher, 1:1000, PA1-180). Appropriate species-specific Donkey secondary antibodies were used at a 1:250 dilution and were obtained from Jackson ImmunoResearch.

**Plasmids.** Elav4 isoform specific overexpression plasmids were commercially obtained from Genecopoeia: *Elav4-v1* (Prod ID: EX-Mm30154-M46), *Elav4-v2* (Prod ID: EX-Mm26736-M46), *Elav4-v3* (Prod ID: EX-Mm29523-M46), and *Elav4-v4* (Prod ID: EX-Mm29524-M46). *Celf1* isoform specific overexpression constructs were obtained from Origene: *Celf1 Long* (Variant 1; Cat# MC217278)

and *Celf1 Short* (Variant 2; Cat# MC216621). *Celf1* shRNAs were commercially obtained from Origene; *Celf1* mouse shRNA Plasmid (Locus ID 13046) (Cat# TF514753B), and Genecopoeia; *Celf1 Variant 1 Long* (Prod ID; MSH031592-mU6) and *Celf1 Variant 2 Short* (Prod ID: MSH043146-mU6). *Elav4-v3* shRNA was custom made (Origene, Cat# HC121119A, TTTAAGAGAAGAGTCGAAGCGCT GCGAGA).

The pCdk5r-Fezf2-ires-GFP plasmid was generously shared as a gift from Paula Arlotta. Briefly, to clone the *Elav4-v3* (Prod ID: EX-Mm29523-M46) and *Elav4-v4* (Prod ID: EX-Mm29524-M46) overexpression under the Cdk5r promoter, the HpaI and XhoI were used to digest the pCdk5r vector and the In-Fusion HD Cloning Plus kit (Clontech, Cat. no. 638909) was used to insert *Elav4-v3* or *-v4* clone. The new construct lacks the original cloning sites. The pNestin-EGFP plasmid was obtained from Addgene (# 38777). The *Elav4-v3* and *Elav4-v4* plasmids above were digested by SacI and SpeI to remove the CMV promoter, and the In-Fusion HD Cloning Plus kit was used to insert the Nestin promoter of pNestin-EGFP plasmid into *Elav4-v3* and *Elav4-v4* vectors. Clones were confirmed through sequencing.

**Elav4 UTR reporter constructs.** The full-length 5' UTRs of *Elav4* variants were amplified from wild type mouse genomic DNA. The *Elav4-v3* 5' UTR sequence (1-97 bp) was amplified with the forward primer 5'-AAAAActtaagGCGCGGGAC CACGTGAGAA-3' and the reverse primer 5'- AAAAAccatggCTTCGCGGAGGC GGGGT-3'. The PCR products were digested with AfIII and NcoI and ligated upstream of the Renilla gene in the LightSwitch 5'UTR reporter vector (SwitchGear Genomics, #S690005). Four *Elav4-v4* 5' UTR subclones (1–341, 1–143, 144–341, and 144–248 bp) were similarly amplified and cloned. All primer sequences can be found in Supplementary Data 4. All other *Elav4-v3* and *Elav4-v4* 5' UTR subclones were synthesized by IDT with sticky ends of the above-mentioned restriction sites for direct ligation. The *Elav4* 3' UTR was amplified with Q5 High-Fidelity DNA Polymerase (NEB) and cloned into the XhoI and NotI sites of psiCheck-2 vector (Promega), downstream of the Renilla luciferase gene. All UTR clones were confirmed with sequencing.

Translationally-relevant binding sites in *Elav4* 5'UTRs were determined by translation assays using subclones described above. These 12 bp sites in the *Elav4-v3* and *Elav4-v4* 5'UTRs were mutated by changing the WT bases to their complements. For reporter gene constructs, full length *Elav4* 5'UTRs with binding site mutated (BSmut) were synthesized by IDT bearing AfIII and NcoI sites at the two ends, and then subcloned upstream of the Renilla CDS in the pLightswitch vector above. The same strategy was used for overexpression constructs, except that the synthesized 5' UTR fragment contains flanking sequence from CMV-*Elav4*-OE CDS vectors (Origene) for subsequent subcloning to SacI and Bsu361 sites. All clones were confirmed by sequencing.

**Elav4 isoform diagram.** The mouse *Elav4* isoform diagrams were retrieved from AceView.

**Celf1 isoforms protein alignment.** Celf1L and Celf1S protein sequences were obtained from USCs genome browser. Protein sequences were aligned using Clustal Omega.

**Immunohistochemistry on mouse brain tissue.** Standard methods were used for IHC as described previously¹⁰. Briefly, embryonic brains and brains from postnatal day 0 (P0) were dissected and fixed by immersion for 6 to 12 h in 4% formaldehyde (PFA, Sigma-Aldrich, Cat. no. 158127) in phosphate buffered saline (PBS, Corning cellgro, Cat. no. 21-040-CV), pH 7.4, as described in Ref. ¹⁰ (Embryonic brains for 6–8 h, P0 brains for 10–12 h). After fixation, the brains were washed three times in 1x PBS and stored in 30% sucrose for later use. Brains were coronally-sectioned at 70–80 µm using a Leica VT1000S vibratome. Free-floating sections were then incubated in blocking solution (normal donkey serum (Jackson ImmunoResearch, Cat. no. 017-000-121), albumin (Biomatik, Cat. no. A2134), 0.2% Glycine (BDH, Cat. no. BDH4156), 0.2% l-lysine (Sigma, Cat. no. L5501), and 0.4% Triton (Sigma, X-100)) gently shaking for 2–4 h at room temperature. Following incubation in blocking solution, the sections were immediately transferred to a solution of primary antibody resuspended in blocking solution +0.4% Triton X-100 and incubated, gently rocking, overnight at 4 °C. The following day, sections were washed 3 times in 1x PBS before being transferred to a secondary antibody solution diluted in blocking solution without Triton X-100 for 1–2 h gently shaking at room temperature. Next, sections were washed again for three times in 1x PBS before being incubated in 1 µg/ml of DAPI (Fisher Scientific, Cat. no. D1306) for 10 min. Finally, sections were washed two more times in 1x PBS before being mounted with Vectashield mounting media (Vector Laboratories, H1000).

**Immunohistochemistry on human brain tissue.** Tissue was fixed in 4% PFA (Biognost, cat. no. FNB4) for up to 48 h, dissected coronally in three blocks, embedded in paraffin (Merck, cat. no. 107300) and sectioned on a microtome (Leica, SM2000R) at 20 µm thick sections. Prior to immunohistochemistry, a standard process of deparaffinization was performed in a series of xylol and alcohol. After four washes in 1x PBS, blocking solution containing 1% BSA and
0.5% Triton X-100 in PBS was applied on sections for 1 h. Blocking solution was replaced with primary antibodies (mouse anti-Elav4, 1:1000; rabbit anti-Pax6, mouse anti-Celf1, 1:200) diluted in blocking solutions and kept overnight at 4 °C. After incubation, sections were 3× washed in 1× PBS and appropriate secondary antibodies (donkey Alexa Fluor 546 (cat. no. A10036) and Alexa Fluor 488 (cat. no. A21206), diluted 1:1000, Thermo Fisher Scientific) were applied for 2 h. Following 3× washes in 1× PBS, DAPI (Sigma Aldrich, St. Louis, USA) was applied for 1 min according to manufacturer instructions. Sections were covered in Aqueous Mounting Medium (DAKO, Carpinteria, CA, USA).

**FISH on mouse brain tissue.** IUE of CAG-GFP and either 5'UTR-Elav4-v3 or 5' UTR-Elav4-v4 under Cdk5r or Nestin promoters was done at E13. Brains were fixed after 24 (Nestin driven constructs) or 36 (Cdk5r driven constructs) hours with RNAse-free 4% PFA for 8 h. From the transfected brains, 70 µm thick sections were permeabilized with 500 µl of PBS-T (1X PBS with 0.5% Triton X-100) in a 24-well plate for 10 min at RT. After a rinse with 1X PBS, the slices were rehydrated in a solution of 10% formamide + 2X SSC for 10 min at RT with 2 ng/µl of a DNA probe (5'UTR-specific 90nt-cy3 or cy5, IDT). Sections were then hybridized in 250 µl of hybridization solution (50% formamide, 5X SSC, 5X Denhardt's solution, 500 ng/µl Salmon Sperm DNA, 250 ng/µl Yeast tRNA) overnight at 37 °C using the HB-100 Hybridizer (UVP Laboratory Products). After two washes for 30 min at 30 °C using 10% formamide + 2X SSC, sections were stained 10 min at room temperature with DAPI and mounted⁷⁰.

**In situ hybridization on human brain tissue.** Fixed and paraffin-embedded coronal sections were mounted on slides, rinsed in 2 × SSC, prehybridized in hybridization buffer (HB) at 45 °C in hybridization oven for 1 h. Hybridization was performed in 100 ng/ml of digoxigenin-labeled anti-sense hElav4-pan-DIG (IDT) or sense probes in HB, at 45 °C overnight. Following hybridization, the sections were washed in 2 × SSC, 2 × SSC/50% formamide, 0.1 × SSC/50% formamide and 0.1 × SSC, at 55 °C. DIG-labeled signal detection was performed with anti-DIG-AP antibody, as recommended by the manufacturer (Roche). Colorimetric detection was done using nitroblue tetrazolium (NBT) and 5-bromo-4-chloro-3-indolyl-phosphate (BCIP) color developing substrate (Promega). Color development reaction was stopped by transferring sections into 10 mM Tris-HCl, pH 8.0, and finally in 1 mM EDTA. Sections were analyzed and images taken with high-resolution slide scanner NanoZoomer 2.0RS (Hamamatsu).

**Neurogenesis analysis.** Animals were bred for timed pregnancies to E13. In addition to standard handling of embryonic tissue, intraperitoneal injections were administered to pregnant dams to allow for thymidine analog incorporation. Cldu (Sigma; no. C6891), dissolved in 1 × PBS, was injected (50 mg/kg) at E12, 24 h prior to extraction of pups. Idu (Sigma; no. I7125), dissolved in 7 mM NaOH in 1 × PBS and made fresh daily, was injected (50 mg/kg) 1 h before brain fixation. Once brain sections were cut at 70 µm thickness, to allow for Cldu/Idu antibody penetration, antigen retrieval was performed. Sections were treated with 1 M HCl for 15 min shaking at room temperature, followed by a 15 min stationary treatment with 2 M HCl. Acid was washed off with 1 × PBS, 4 washes 5 min each. Our standard IHC protocol was followed thereafter. Cortex was imaged with a 20x objective. Files were opened in Fiji (NIH; ImageJ ver. 2.0.0-rc-43/1.51r) and all cells in cortex were marked using the Cell Counter plugin (license GPLv3). Cells were counted if they were between the pial surface and the VZ, and if they were positive for DAPI staining. Positive staining in individual channels was determined first. Then colocalizations were determined by the presence of two different markers on the same cell. Two-tailed unpaired t-tests were used to test comparisons for significant differences.

**Confocal imaging.** All images except from human were taken with an Olympus BX61WI confocal microscope using 10 ×, 20×, or 60× objectives and processed using Fluoview FV-1000. Human IHC was analyzed using Leica TCS SP8 X FLIM confocal microscope. All representative images and images used in analysis were taken with the same confocal settings per experiment to allow for accurate comparisons of fluorescent intensity. Sections were imaged from the top of the cortical plate to the VZ and merged in the photo editing software Gimp2 and Fiji⁷¹. For quantifications of binned images, a rectangle approximately a third of the width of the image was drawn from CC to the tops of the apical tufts at the pia. The rectangle with split into ten equal bins from CC to pia. Cells were counted per condition as appropriate for the experiment, either GFP+ or GFP+ colocalized with a layer marker. DAPI was used to confirm cell nuclei.

**RNA immunoprecipitation.** Six of E13 or E16 neocortices were dissected, pooled together, and considered to be one biological sample. Three biological samples were used for RNA immunoprecipitation (RIP) experiments. Celf1 RNA immunoprecipitation validated antibody (Millipore, Cat. no. 03-104) was used coupled to EZ-Magna RIP RNA-Binding Protein Immunoprecipitation Kit (Millipore, Cat. no. 17-701) and qRT-PCR analysis. For each target tested, three qRT-PCR technical replicates were done from each immunoprecipitation of each biological sample. To examine Celf1 binding to Elav4 5' UTR, the Elav4-v3 (or v4) 5' UTR-Renilla + Celf1S OE constructs (0.2 and 1 µg/12-well plate) were co-transfected into N2a

cells, which were then incubated for 3 days. 10⁷ cells were used for one RIP experiment of IgG or Celf1 and followed by the Millipore RIP kit. The effect of Celf1 binding to Elav4-v3 (or v4) 5' UTR was measured with a Renilla TaqMan probe. 3–4 replicates of qRT-PCR were performed for each RIP.

**Western blot.** Neocortical tissue and N2a cells were lysed in Tissue Protein Extraction Reagent (Thermo Scientific, Cat. no. 78510) and protein concentrations were determined using the Pierce 660 nm reagent (Thermo Fisher, Cat. no. 22660) on a Nanodrop ND-1000 spectrophotometer. Protein samples were analyzed using the NuPAGE system (Life Technologies) per manufacturer's instructions and with 4-12% Bis-Tris gels (Cat. no. NP0335/6) and transferred to nitrocellulose membranes (GVS Life Sciences, Cat. no. 1215471). The resulting membranes were blocked for 1 h at room temperature with 5% nonfat dried milk (VWR, Cat. no. M203). and 10% fetal bovine serum (FBS) (Gemini, 900-108) in PBS with 0.4% Triton-X-100 (PBST). Primary antibody was added and the blots were incubated overnight at 4 °C gently shaking, followed by three washes in PBST. Secondary antibody diluted in 10% FBS in PBST was then added for 1 h. Development of blots was performed with ChemiGlow West Chemiluminescence Substrate Kit (Proteinsimple, Cat. no. 60-12596-00-2) and visualized with Gbox (Syngene). Quantification of band intensity was performed with GeneTools (Syngene) and Fiji⁷¹. The experimental protein of interest was normalized to Gapdh levels on the same blot. At least three biological replicates were performed for all analyses and an ANOVA or Student's t-test was used to compare averages, as appropriate. *p* < 0.05 was considered significant.

**Laser capture micro dissection.** Whole brains were dissected from developing embryos at E13 or E16 in an RNAse-free environment and flash frozen. Brains were then sectioned into 60 µm slices on a cryostat at -20 °C and placed onto Molecular Machines & Industries (MMI) RNAse Free Membrane Slides (MMI AG, Cat. no. 50102). The sections were rinsed with nuclease-free water and immersed for 1 min in 95% ethanol. Dried slides were subjected to microdissection using an MMI SmartCut Plus microscope. Neocortical compartments were determined by the cytoarchitecture of cresyl violet-stained tissue. Brains underwent laser capture microdissection of the ventricular zone and cortical plate. RNA was harvested from compartmentalized tissue using the RNAqueous-micro kit LCM protocol (Ambion). RNA from at least three brains was analyzed per experiment (n = 3) using quantitative real-time PCR (qRT-PCR, see below). Significant changes between conditions were analyzed using the ΔΔCt method. A student's t-test or ANOVA was performed to evaluate significant changes and *p* < 0.05 was considered significant.

**Preparation and FACS sorting of GFP+ cells.** WT female CD1 mice were crossed with Nestin-GFP or Tbr2-GFP males for timed pregnancies. The morning of plug observation was considered to be E0.5. Embryos were taken at E13 or E16 and observed under a dissecting microscope for GFP expression. GFP+ neocortices were dissected, pooled together, and single cell suspension per pregnant animal was prepared as described¹³. GFP+ cells were then FACS-sorted at the local Rutgers University Facility, spun down, and immediately homogenized in TRIzol (Invitrogen) for further RNA and protein isolation. Average values for cells obtained (from three separate spins): E13 Nestin-GFP, 348 K cells; E16 Nestin-GFP, 1.73 M; E13 Tbr2-GFP, 128 K; E16 Tbr2-GFP, 525 K. Sample purity was determined through ROI determination as well as post-sort visual examination on a fluorescent microscope at the facility.

**Quantitative real-time PCR (qRT-PCR).** RNA was isolated using TRIzol (Invitrogen) per manufacturer instructions and as previously described¹⁰. After ethanol was removed from the precipitate, the mRNA pellet was resuspended in 20–30 µL of water. Residual DNA was removed by incubating the mRNA with Turbo DNAse in DNAse buffer, and then the DNAse was inactivated with the DNAse inactivation reagent (Invitrogen, AM1907). The mRNA concentration was measured using a Nanodrop ND-1000 spectrophotometer. For LCM samples, 10–25 ng of total RNA was used per reaction. For all other reactions, 50 ng per reaction was used. The Applied Biosystems StepOne Real-Time System and reagents (TaqMan RNA-to-Ct 1-Step Kit, Thermo Fisher, #4392938) were used to perform qRT-PCR. The results were analyzed using the ΔΔCt method with Gapdh as a normalization control unless otherwise specified. For qRT-PCR of Elav4 variants, commercially available TaqMan probes (Thermo Fisher) were obtained.

**In utero electroporation.** In utero electroporation (IUE) was performed at E12, E13, or E16 and analyzed at E13, E17, P0, or P7 as previously described⁹. Co-transfections accomplished with approximately 1 µL of mix containing 3.5–4 µg/µL of the control vector or the vector of interest along with 4 µg/µL CAG-GFP reporter in a 4:1 ratio (construct:reporter). For rescue experiments, 8 µg/µL of the shRNA and overexpression vectors were combined and injected to achieve the final concentration of the OE experiments (~3.5–4 µg of construct, ~800–900 ng of reporter, per µL). For each IUE, at least three transfected neocortices were used in experiments and this is indicated in figure legends with the “n” value. Two to five sections from each transfected per staining were used in quantifications.
**Cell culture/transfections.** Neuroblastoma N2a cells (ATCC, Cat. no. CCL-131) were grown in Dulbecco's modified Eagle medium containing 10% FBS (Gemini, Cat. no. 900-108), 1% GlutaMAX (Gibco, 35050-061), 1% sodium pyruvate (Gibco, Cat. no. 11360-070), and 1% penicillin streptomycin (Corning, Cat. no. 30-001-CI). TrypLE Express (Gibco, Cat. no. 12604-021) dissociation reagent was used for regular maintenance. Lipofectamine 2000 (Invitrogen, Cat. no. 11668-019) was used to perform transfections per the manufacturer's protocol.

**Renilla luciferase reporter assay.** N2a cells were seeded in 12-well plates and co-transfected the next day with 0.2 µg of either of the *Elavl4* 5' UTR-Renilla reporter vectors and 1 µg of either the Ctrl OE, *Celf1S* OE, or *Celf1L* OE vectors. Forty-eight hours after transfection, cells were trypsinized with TrypLE Express (Gibco, Cat. no. 12604-021) and divided equally into two fractions. One fraction of the cells was spun down and lysed at room temperature for 15 min with 250 µL of 1x Passive Lysis Buffer (Promega, in Ref# E1910 Kit). The other fraction of cells was used for RNA isolation using TRIzol (Invitrogen). All protein lysates were quantified with the Pierce 660 nm Protein Assay reagent (Thermo Fisher) on a NanoDrop spectrophotometer and normalized to the same concentration with lysis buffer for equal loading. To measure the Renilla luminescence (REN), 10 µL of normalized protein lysate was mixed with 50 µL of working solution made from the LightSwitch Luciferase Assay Reagent (SwitchGear Genomics, Cat. no. 32035). REN readings were immediately performed on a GloMax 20/20 luminometer (Promega). The levels of Renilla mRNA (*Ren*) with 50 ng of Turbo DNase (Invitrogen, Cat. no. AM2238) treated RNA fraction were quantified by qRT-PCR using the custom Taqman probe Rensp (Thermo Fisher). The REN luminescence readings were then normalized to *Ren* mRNA levels for each sample. Each assay condition was performed with at least three biological replicates. ANOVA or Student t-test statistical comparisons of REN/*Ren* between experimental conditions were performed. *P* < 0.05 was considered significant. The same approach was used for experiments with *Celf1* knockdown except that cells were harvested 72 h after co-transfection of reporter plasmids with the control shRNA or *Celf1* shRNA plasmid. Similar experiments were also conducted for *Elavl4* 3' UTR reporter gene assay, except that the Dual-Luciferase Reporter Assay kit (Promega) was used for determining the Firefly Luciferase luminescence.

**mRNA decay assay.** N2a cells were seeded in 100 mm dishes and transfected the next day with 30 µg of either the Ctrl OE, *Celf1S* OE, or *Celf1L* OE plasmids. Twenty-four hours later the bulk transfected were reseeded evenly into 35 mm dishes and rested overnight. The cells were then treated with Actinomycin D (Sigma, Cat. no. A9415) at a final concentration of 5 µg/mL and incubated for 15 min to stop transcription before starting the time point. Cells were then harvested at 0, 2, 4 and 6 h for RNA isolation and qRT-PCR analysis of mRNA decay. The same approach was used for experiments with *Celf1* knockdown (where the control shRNA or *Celf1* shRNA plasmids were bulk transfected) except that cells were left in 35 mm dishes for 48 h after replating and before Actinomycin D treatment.

**Primary neuronal cultures.** Alternating pups were co-electroporated with Ctrl shRNA and CAG-GFP or *Celf1* shRNA and CAG-RFP and primary neuronal cultures were performed as described¹³.

**Constraint scores.** Constraint metrics were obtained from the Exome Aggregation Consortium (ExAC), Cambridge, MA (URL: http://exac.broadinstitute.org) [May 2016]. We used the missense Z-score, computed as deviation of observed number of missense variants from the expected number, as a measure of intolerance to missense variation⁴⁶. We used the probability of loss-of-function intolerance (pLI) as a measure of intolerance to loss-of-function mutations⁴⁶. For the partitioning of genes in pLI tranches, deviation from the distribution observed for all genes was evaluated using Kolgomorov-Smirnov test (ks.test, R).

**Enrichment analyses.** We extracted the 1216 unique PAR-CLIP targets of ELAVL1 (Lebedeva et al., 2011⁷², Supplementary Data 2, genes described as "Conservative targets", PMID:21723171), 924 unique top HITS-CLIP targets of nELAVLs (Scheckel et al., 2016⁷³, Supplementary File 1, PMID: 26894958), 431 human orthologs of the top 490 iCLIP targets of CELF4 (Wagnon et al., 2012⁷⁴, Supplementary Data 1, PMID:2320943), 1439 unique HITS-CLIP targets of CELF1 (Le Tonquèze et al.⁷⁵, Supplementary Data 3, genes described as "inLETONQUEZE=1", PMID:27222809). We also used 915 unique genes in the Developmental Disorders Genotype-Phenotype Database (DDG2P, https://decipher.sanger.ac.uk/info/ddg2p) with "brain" included in organ specificity (referred to as NDD_DDG2P), 61 unique gene associated with EE in the Online Mendelian Inheritance in Man (OMIM) database (EE_OMIM), 142 unique genes associated with ID in OMIM (ID_OMIM), and 102 unique genes associated with ASD (ASD_Satterstrom) (PMID: 31981491). For the enrichment analyses, we used ELAVL1, nELAVLs, CELF4 and CELF1 genes are input sets, the NDD_DDG2P, EE_OMIM, ID_OMIM, and ASD_Satterstrom as target sets, and 19,584 unique human genes as background set. For each comparison between the input and target sets, we first constructed the empirical distribution by sampling the same number of genes as in the input set from the background set 10,000 times. The P value was then computed by calculating the number of sampled gene lists that had at

least as many overlapping genes with the target sets as the input set, divided by 10,000 iterations.

The 915 NDD_DDG2P genes were also crossed with the 850 human orthologs of the 917 unique mouse genes de-repressed at E16, and with the 1,070 human orthologs of the 1,093 unique mouse genes repressed at E16. Overlapping genes were further mapped to OMIM.

**Quantification.** Cell counting was done in double blind fashion where neither the person imaging nor quantifying knew the experimental condition. The percentage of colocalization between GFP and the glutamatergic identity markers of interest was quantified in Fiji and is presented as the fraction of total GFP+ cells. For the cell migration analysis, 10 equal sized bins from either (1) the top of the cortical plate to bottom of the cortical plate or (2) from the top of the cortical plate to the bottom of the VZ were drawn for each experiment. Binning conditions were kept constant across an experiment and the number of GFP+ cells in each bin was counted to be presented as the fraction of total GFP+ cells.

**Statistical analysis.** Graph bars represent mean and ±SEM or mean and ±SD as noted in figure legends. The number (n) of replicates for each experiment is also noted in figure legends. Appropriate statistical tests were used (either unpaired Student's t-test or one way ANOVA Tukey's post-hoc test for multiple comparisons), also noted in figure legends. Statistical significance is considered achieved when *p* < 0.05 and is reported as: \* *p* < 0.05, \*\* *p* < 0.01, \*\*\* *p* < 0.001, \*\*\*\* *p* < 0.0001.

**Reporting summary.** Further information on research design is available in the Nature Research Reporting Summary linked to this article.

## Data availability

The unprocessed sequencing data have been deposited in the NIH GEO with accession number: GSE77647. Processed data used for the analyses are included in Supplementary Data files 1–3. The source data underlying Figs. 1g–j, 2b, 3b, 3e, 4b, 4d, 4g, 4h, 5a–e, 6d, 7b–e, 8e–h and Supplementary Figs. 1a, 2b, c, 4c–e, 5a, 5c–e, 6a–h, 7b, c, 8c, 9b, c are provided as a Source Data file.

Received: 4 October 2018; Accepted: 3 March 2020;
Published online: 03 April 2020

## References

1. Kwan, K. Y., Sestan, N. & Anton, E. S. Transcriptional co-regulation of neuronal migration and laminar identity in the neocortex. *Development* **139**, 1535–1546 (2012).

2. Pilaz, L. & Silver, D. Post-transcriptional regulation in corticogenesis: how RNA-binding proteins help build the brain. *Wiley Interdiscip. Rev. RNA* **6**, 501–515 (2015).

3. Popovitchenko, T. & Rasin, M. R. Transcriptional and post-transcriptional mechanisms of the development of neocortical lamination. *Front. Neuroanat.* **11**, 102 (2017).

4. Breunig, J., Haydar, T. & Rakic, P. Neural stem cells: historical perspective and future prospects. *Neuron* **70**, 614–625 (2011).

5. DeBoer, E., Kraushar, M., Hart, R. & Rasin, M. Post-transcriptional regulatory elements and spatiotemporal specification of neocortical stem cells and projection neurons. *Neuroscience* **248**, 499–528 (2013).

6. Lodato, S. & Arlotta, P. Generating neuronal diversity in the mammalian cerebral cortex. *Annu. Rev. Cell Dev. Biol.* **31**, 699–720 (2015).

7. Lennox, A., Mao, H. & Silver, D. RNA on the brain: emerging layers of post-transcriptional regulation in cerebral cortex development. *Wiley Interdiscip. Rev. Dev. Biol.* **7**, e290 (2018).

8. Zhang, X. et al. Cell-type-specific alternative splicing governs cell fate in the developing cerebral cortex. *Cell* **166**, 1147–1162 (2016).

9. DeBoer, E. et al. Prenatal deletion of the RNA-binding protein HuD Disrupts postnatal cortical circuit maturation and behavior. *J. Neurosci.* **34**, 3674–3686 (2014).

10. Kraushar, M. L. et al. Temporally defined neocortical translation and polysome assembly are determined by the RNA-binding protein Hu antigen R. *Proc. Natl Acad. Sci. USA* **111**, E3815–E3824 (2014).

11. Kraushar, M. L. et al. Thalamic WNT3 secretion spatiotemporally regulates the neocortical ribosome signature and mRNA translation to specify neocortical cell subtypes. *J. Neurosci.* **35**, 10911–10926 (2015).

12. Yang, G., Smibert, C., Kaplan, D. & Miller, F. An eIF4E1/4E-T complex determines the genesis of neurons from precursors by translationally repressing a proneurogenic transcription program. *Neuron* **84**, 723–739 (2014).

13. Popovitchenko, T. et al. The RNA binding protein HuR determines the differential translation of autism-associated FoxP subfamily members in the developing neocortex. *Sci. Rep.* **6**, 28998 (2016).
14. Chau, K. F. et al. Downregulation of ribosome biogenesis during early forebrain development. *Elife* 7, e36998 (2018).

15. Zahr, S. et al. Translational repression complex in developing mammalian neural stem cells that regulates neuronal specification. *Neuron* **97**, 520–537 (2018).

16. Jackson, R., Hellen, C. & Pestova, T. The mechanism of eukaryotic translation initiation and principles of its regulation. *Nat. Rev. Mol. Cell Biol.* **11**, 113–127 (2010).

17. Hinnebusch, A. G., Ivanov, I. P. & Sonenberg, N. Translational control by 5'-untranslated regions of eukaryotic mRNAs. *Science* **352**, 1413–1416 (2016).

18. Tushev, G. et al. Alternative 3'UTRs modify the localization, regulatory potential, stability, and plasticity of mRNAs in neuronal compartments. *Neuron* **98**, 495–511 (2018).

19. Simone, L. & Keene, J. Mechanisms coordinating ELAV/Hu mRNA regulons. *Curr. Opin. Genet. Dev.* **23**, 35–43 (2013).

20. Gardiner, A. S., Twiss, J. L. & Perrone-Bizzozero, N. I. Competing interactions of RNA-binding proteins, microRNAs, and their targets control neuronal development and function. *Biomolecules* **5**, 2903–2918 (2015).

21. Theil, K., Herzog, M. & Rajewsky, N. Post-transcriptional regulation by 3' UTRs can be masked by regulatory elements in 5' UTRs. *Cell Rep.* **22**, 3217–3226 (2018).

22. Ayoub, A. et al. Transcriptional programs in transient embryonic zones of the cerebral cortex defined by high-resolution mRNA sequencing. *Proc. Natl Acad. Sci. USA* **108**, 14950–14955 (2011).

23. Liu, L. et al. Competition between RNA-binding proteins CELF1 and HuR modulates MYC translation and intestinal epithelium renewal. *Mol. Biol. Cell* **26**, 1797–1810 (2015).

24. Timchenko, N., Lu, A., Welm, X. & Timchenko, L. CUG repeat binding protein (CUGBP1) interacts with the 5' region of C/EBP mRNA and regulates translation of C/EBP isoforms. *Nucleic Acids Res.* **27**, 4517–4525 (1999).

25. Lee, J. & Cooper, T. Pathogenic mechanisms of myotonic dystrophy. *Biochem. Soc. Trans.* **37**, 1281–1286 (2009).

26. Reimand, J. et al. g:Profiler—a web server for functional interpretation of gene lists (2016 update). *Nucleic Acids Res.* **44**, W83–W89 (2016).

27. Kraushar, M. L., Popovitchenko, T., Volk, N. L. & Rasin, M. R. The frontier of RNA metamorphosis and ribosome signature in neocortical development. *Int. J. Dev. Neurosci.* **55**, 131–139 (2016).

28. Akamatsu, W. et al. The RNA-binding protein HuD regulates neuronal cell identity and maturation. *Proc. Natl. Acad. Sci. USA* **102**, 4625–4630 (2005).

29. Wang, F. et al. Positive feedback between RNA-binding protein HuD and transcription factor SATB1 promotes neurogenesis. *Proc. Natl Acad. Sci. USA* **112**, E4995–E5004 (2015).

30. Mignone, J. L., Kukekov, V., Chiang, A., Steindler, D. & Enikolopov, G. Neural stem and progenitor cells in nestin-GFP transgenic mice. *J. Comp. Neurol.* **469**, 311–324 (2004).

31. Mignone, J., Peunova, N. & Enikolopov, G. Nestin-based reporter transgenic mouse lines. *Methods Mol. Biol.* **1453**, 7–14 (2016).

32. Arnold, S., Sugnaseelan, J., Groszer, M., Srinivas, S. & Robertson, E. Generation and analysis of a mouse line harboring GFP in the Eomes/Tbr2 locus. *Genesis* **47**, 775–781 (2009).

33. Arnold, S. J. et al. The T-box transcription factor Eomes/Tbr2 regulates neurogenesis in the cortical subventricular zone. *Genes Dev.* **22**, 2479–2484 (2008).

34. Tsuda, K. et al. Structural basis for the sequence-specific RNA-recognition mechanism of human CUG-BP1 RRM3. *Nucleic Acids Res.* **37**, 5151–5166 (2009).

35. Teplova, M., Song, J., Gaw, H. Y., Teplov, A. & Patel, D. J. Structural insights into RNA recognition by the alternate-splicing regulator CUG-binding protein 1. *Structure* **18**, 1364–1377 (2010).

36. Kress, C., Gautier-Courteille, C., Osborne, H., Babinet, C. & Paillard, L. Inactivation of CUG-BP1/CELF1 causes growth, viability, and spermatogenesis defects in mice. *Mol. Cell. Biol.* **27**, 1146–1157 (2006).

37. Cibois, M., Boulanger, G., Audic, Y., Paillard, L. & Gautier-Courteille, C. Inactivation of the Celf1 gene that encodes an RNA-binding protein delays the first wave of spermatogenesis in mice. *PLoS ONE* **7**, e46337 (2012).

38. Carlyle, B. C. et al. A multiregional proteomic survey of the postnatal human brain. *Nat. Neurosci.* **20**, 1787–1795 (2017).

39. Xu, R. et al. OLIG2 drives abnormal neurodevelopmental phenotypes in human iPSC-based organoid and chimeric mouse models of down syndrome. *Cell Stem Cell* **24**, 908–926 (2019).

40. Rouaux, C. & Arlotta, P. Direct lineage reprogramming of post-mitotic callosal neurons into corticofugal neurons in vivo. *Nat. Cell Biol.* **15**, 214–221 (2013).

41. Fenlon, L. & Richards, L. Contralateral targeting of the corpus callosum in normal and pathological brain function. *Trends Neurosci.* **38**, 264–272 (2015).

42. Pollen, A. A. et al. Molecular identity of human outer radial glia during cortical development. *Cell* **163**, 55–67 (2015).

43. Nowakowski, T. J. et al. Spatiotemporal gene expression trajectories reveal developmental hierarchies of the human cortex. *Science* **358**, 1318–1323 (2017).

44. Zunic Isasegi, I. et al. Interactive histogenesis of axonal strata and proliferative zones in the human fetal cerebral wall. *Brain Struct. Funct.* **223**, 3919–3943 (2018).

45. Judaš, M. et al. The Zagreb collection of human brains: a unique, versatile, but underexploited resource for the neuroscience community. *Ann. N. Y. Acad. Sci.* **1225**, E105–E130 (2011).

46. Lek, M. et al. Exome Aggregation Consortium. Analysis of protein-coding genetic variation in 60,706 humans. *Nature* **536**, 285–291 (2016).

47. Satterstrom, F. K. et al. Large-scale exome sequencing study implicates both developmental and functional changes in the neurobiology of autism. *Cell* **180**, 568–584 (2020).

48. Kwan, K. Y. et al. Species-dependent posttranscriptional regulation of NOS1 by FMRP in the developing cerebral cortex. *Cell* **149**, 899–911 (2012).

49. Darnell, J. et al. FMRP stalls ribosomal translocation on mRNAs linked to synaptic function and autism. *Cell* **146**, 247–261 (2011).

50. Fukao, A. et al. The ELAV protein HuD stimulates cap-dependent translation in a poly(A)- and eIF4A-dependent manner. *Mol. Cell* **36**, 1007–1017 (2009).

51. Hayashi, S., Yano, M., Igarashi, M., Okano, H. J. & Okano, H. Alternative role of HuD splicing variants in neuronal differentiation. *J. Neurosci. Res.* **93**, 399–409 (2014).

52. Yokoi, S. et al. 3'UTR length-dependent control of SynGAP isoform α2 mRNA by FUS and ELAV-like proteins promotes dendritic spine maturation and cognitive function. *Cell Rep.* **20**, 3071–3084 (2017).

53. Zybura-Broda, K. et al. HuR (Elavl1) and HuB (Elavl2) stabilize matrix metalloproteinase-9 mRNA during seizure-induced Mmp-9 expression in neurons. *Front. Neurosci.* **12**, 1–15 (2018).

54. Darnell, J. & Klann, E. The translation of translational control by FMRP: therapeutic targets for FXS. *Nat. Neurosci.* **16**, 1530–1536 (2013).

55. De Rubeis, S. et al. CYFIP1 coordinates mRNA translation and cytoskeleton remodeling to ensure proper dendritic spine formation. *Neuron* **79**, 1169–1182 (2013).

56. De Rubeis, S. et al. Synaptic, transcriptional and chromatin genes disrupted in autism. *Nature* **515**, 209–215 (2014).

57. De Rubeis, S. et al. CYFIP1 coordinates mRNA translation and cytoskeleton remodeling to ensure proper dendritic spine formation. *Neuron* **79**, 1169–1182 (2013).

58. De Rubeis, S. et al. Synaptic, transcriptional and chromatin genes disrupted in autism. *Nature* **515**, 209–215 (2014).

59. Dasgupta, T. & Ladd, A. N. The importance of CELF control: molecular and biological roles of the CUG-BP, Elav-like family of RNA-binding proteins. *Wiley Interdiscip. Rev. RNA* **3**, 104–121 (2012).

60. Dougherty, J. D. et al. The disruption of Celf6, a gene identified by translational profiling of serotonergic neurons, results in autism-related behaviors. *J. Neurosci.* **33**, 2732–2753 (2013).

61. Parras, A. et al. Autism-like phenotype and risk gene mRNA deadenylation by CPEB4 mis-splicing. *Nature* **560**, 441–446 (2018).

62. Geschwind, D. H. & Levitt, P. Autism spectrum disorders: developmental disconnection syndromes. *Curr. Opin. Neurobiol.* **17**, 103–111 (2007).

63. State, M. W. & Levitt, P. The conundrums of understanding genetic risks for autism spectrum disorders. *Nat. Neurosci.* **14**, 1499–1506 (2011).

64. Sestan, N. & State, M. W. Lost in translation: traversing the complex path from genomics to therapeutics in autism spectrum disorder. *Neuron* **100**, 406–423 (2018).

65. Aitchison, J. The Statistical Analysis of Compositional Data Monographs on Statistics and Applied Probability. (Chapman & Hall Ltd., London, UK, 1986).

66. Trapnell, C. et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. *Nat. Protoc.* **7**, 562–578 (2012).

67. Durinck, S., Spellman, P., Birney, E. & Huber, W. Mapping identifiers for the integration of genomic datasets with the R/Bioconductor package biomaRt. *Nat. Protoc.* **4**, 1184–1191 (2009).

68. Heinz, S. et al. Simple Combinations of Lineage-Determining Transcription Factors Prime cis-Regulatory Elements Required for Macrophage and B Cell Identities. *Mol. Cell* **38**, 576–589 (2010).

69. Gerstberger, S., Hafner, M. & Tuschl, T. A census of human RNA-binding proteins. *Nat. Rev. Genet.* **15**, 829–845 (2014).

70. Park, S. W., Kuroda, M. I., & Park, Y. Regulation of histone H4 Lys16 acetylation by predicted alternative secondary structures in roX noncoding RNAs. *Mol Cell Biol.* **28**, 4952–4962 (2008).

71. Schindelin, J. et al. Fiji: an open-source platform for biological-image analysis. *Nat. Methods* **9**, 676–682 (2012).

72. Lebedeva, S., et al. Transcriptome-wide Analysis of Regulatory Interactions of the RNA-Binding Protein HuR. *Mol. Cell* **43**, 340–352 (2011).
73. Scheckel, C. et al., Regulatory Consequences of Neuronal ELAV-like Protein Binding to Coding and Non-Coding RNAs in Human Brain. *Elife*. **5**, e10421 (2016).

74. Wagnon, JL. et al., CELF4 Regulates Translation and Local Abundance of a Vast Set of mRNAs, Including Genes Associated With Regulation of Synaptic Function. *PLoS Genet*. **8**, e1003067 (2012).

75. Le Tonquèze O, Gschloessl B, Legagneux V, Paillard L, & Audic Y. Identification of CELF1 RNA targets by CLIP-seq in human HeLa cells. *Genomics Data* **8**, 97–103 (2016)

## Acknowledgements

This work was supported by National Institutes of Health (NIH) grants (NS064303 and NS075367) and Robert Wood Johnson Medical School start-up funds (M.R.R.). This work was also supported in part by the Assistant Secretary of Defense for Health Affairs endorsed by the Department of Defense, through the Epilepsy Research Program under Award No. W81XQH-18-1-0338 to H.Z. and M.R.R. Opinions, interpretations, conclusions and recommendations are those of the author and are not necessarily endorsed by the Department of Defense. The U.S. Army Medical Research Acquisition Activity, 820 Chandler Street, Fort Detrick MD 21702-5014 is the awarding and administering acquisition office. In addition, this work was supported in part by the "Research Cooperability" Program of the Croatian Science Foundation funded by the European Union from the European Social Fund under the Operational Programme Efficient Human Resources 2014-2020 PSZ-2019-02-4710 to Z.K. and M.R.R. This publication was also supported by Adris Foundation (ZK), Croatian Science Foundation projects CSF-IP-09-2014-4517 (IK) and co-financed by the Scientific Centre of Excellence for Basic, Clinical and Translational Neuroscience (project "Experimental and clinical research of hypoxic-ischemic damage in perinatal and adult brain"); GA KK01.1.1.01.0007 funded by the European Union through the European Regional Development Fund), and NIH grant NS089578 to HZ. N.F.P. received support from the Aresty Undergraduate Research Fellowship and the Thermo Fisher Scientific Antibody Scholarship Program. J.D.S. was funded by an Undergraduate Research Award from the Rutgers Honors College-New Brunswick. M.L.K. was supported by a NIH NRSA F30 MD/PhD Fellowship (1F30MH106220). S.D.R. is supported by the Beatrice and Samuel A. Seaver Foundation and the Fascitelli Research Scholar Award. This study was supported by the Germany's Excellence Strategy (CIBSS-EXC-2189- Project ID 390939984) to S.J.A. We thank Ishan Khosla for help with text editing, and Erik DeBoer, Aditi Dubey, Althea Stillman, Elizaveta Muravieva, Jack Pike, Faizan Rafi, and Midori Lofton for technical support, Dr. Paola Arlotta (Harvard University) for sharing the Cdk5r construct, Dr. Grigori Enikolopov (Stony Brook University) for *Nestin-GFP* mice, Drs. Lori Covey (Rutgers University) and Mengqing Xiang (Rutgers University) for luminometer usage, and Dr. Paul Copeland for microultracentrifuge usage.

## Author contributions

Conceptualization, M.R.R., T.P., and N.F.P.; Investigation, T.P., Y.P., N.F.P., X.L., Z.K., Y.L., J.D.S., I.S. M.L.K., N.L.V., S.M.P., H.R.S.W., D.L., S.S.K., A.W., and M.R.R.; Formal analysis R.P.H., S.D.R., N.F.P., T.P., X.L., Y.P., S.B., and M.R.R.; Writing—Original Draft, M.R.R., T.P., and N.F.P.; Writing—Review and Editing, T.P., N.F.P., X.L., Y.P., I.S., H.O., S.J.A. R.P.H., S.D.R., H.Z., Z.K., H.O., and M.R.R.; Funding acquisition, M.R.R.; Resources, M.S., H.Z., S.J.A., W.A., H.O.; and L.P. S.D.R., R.P.H. and I.K. provided expertise and feedback.

## Competing interests

The authors declare no competing interests.

## Additional information

Supplementary information is available for this paper at https://doi.org/10.1038/s41467-020-15412-8.

Correspondence and requests for materials should be addressed to M.-R.R.

Peer review information *Nature Communications* thanks the anonymous reviewer(s) for their contribution to the peer review of this work.

Reprints and permission information is available at http://www.nature.com/reprints

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_62.png" /></figure>

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2020