NOVA2-mediated RNA regulation is required for axonal pathfinding during development
Abstract
The neuron specific RNA-binding proteins NOVA1 and NOVA2 are highly homologous alternative splicing regulators. NOVA proteins regulate at least 700 alternative splicing events in vivo, yet relatively little is known about the biologic consequences of NOVA action and in particular about functional differences between NOVA1 and NOVA2. Transcriptome-wide searches for isoform-specific functions, using NOVA1 and NOVA2 specific HITS-CLIP and RNA-seq data from mouse cortex lacking either NOVA isoform, reveals that NOVA2 uniquely regulates alternative splicing events of a series of axon guidance related genes during cortical development. Corresponding axonal pathfinding defects were specific to NOVA2 deficiency: Nova2-/- but not Nova1-/- mice had agenesis of the corpus callosum, and axonal outgrowth defects specific to ventral motoneuron axons and efferent innervation of the cochlea. Thus we have discovered that NOVA2 uniquely regulates alternative splicing of a coordinate set of transcripts encoding key components in cortical, brainstem and spinal axon guidance/outgrowth pathways during neural differentiation, with severe functional consequences in vivo.
https://doi.org/10.7554/eLife.14371.001eLife digest
The first step of producing a protein involves the DNA of a gene being copied to form a molecule of RNA. This RNA molecule can often be processed to create several different “messenger” RNAs (mRNAs), each of which are used to produce a different protein by a process known as alternative splicing. A class of proteins that bind to RNA molecules controls alternative splicing. These “splicing factors” ensure that the right protein variant is produced at the right time and in the right place to carry out the appropriate activity.
Many genes that play important roles in the nervous system have been reported to undergo alternative splicing to generate different protein variants. However, it is unclear whether alternative splicing is important for controlling how the nervous system develops, during which time the neurons connect to the cells that they will communicate with. Forming these connections involves part of the neuron, called the axon, growing along a precise path through the nervous system to reach its destination.
Two RNA-binding proteins called NOVA1 and NOVA2 are produced exclusively in the central nervous system, where they regulate a number of actions including alternative splicing. So far, no differences in the roles of NOVA1 and NOVA2 have been identified, and relatively little is known about their actions in the brain.
Saito et al. have addressed these missing puzzle pieces by combining RNA analysis methods with an analysis of the structure of the nervous system of mice that lack either NOVA1 or NOVA2. This approach identified where NOVA1 and NOVA2 bind on mRNAs, and showed that the mRNAs are processed in different ways in the developing mouse brain depending on which form of the NOVA protein is bound to it.
Further analysis of the data revealed that NOVA2, and not NOVA1, regulates splicing in a series of RNA molecules that help to guide axons to the correct locations in the developing mouse brain. A related study by Leggere et al. also reported on the role that NOVA proteins play in the alternative splicing of one of these genes, called Dcc.
Saito et al. also found defects in the nervous systems of the mice that lacked NOVA2 that only occurred in these mice and resulted from certain axons being unable to follow the correct path to their target cells. These led to major defects, such as agenesis of the corpus callosum (a complete lack of connection between the right and left sides of the brain). Further defects affected how specific subsets of motor neurons connect to muscles and how cochlear neurons in the brainstem connect to the inner ear. The next steps are to explore how the processing of RNA molecules by NOVA2 causes these defects, and to assess whether these actions relate to developmental brain disorders in humans.
https://doi.org/10.7554/eLife.14371.002Introduction
During central nervous system (CNS) development, a neuron extends its axon through a complex yet precise path to reach its final destination by sensing extracellular signals called guidance cues. These cues are sensed by the growth cone, a motile structure at the extending axon edge, and they control growth cone motility through directed cytoskeletal remodeling. Netrins, slits, semaphorins, and ephrins are the major classic guidance cues and elicit attractive or repulsive responses in growth cones via specific receptors (Brose et al., 1999; Cheng et al., 1995; Drescher et al., 1995; Fan and Raper, 1995; Kapfhammer and Raper, 1987; Kennedy et al., 1994; Kidd et al., 1999; Serafini et al., 1994). An important aspect of axon guidance is the spatial and temporal control of response to the guidance cues. For example, the spinal cord commissural axon reaching the midline senses netrin-1, secreted from the floorplate as a chemoattractive cue; however, once it has crossed the floorplate, this cue becomes repulsive (Kennedy et al., 1994; Kidd et al., 1998; Tessier-Lavigne et al., 1988; Zou et al., 2000). Furthermore, the spatiotemporally restricted expression of Robo3 alternative splicing isoforms in spinal cord commissure axons are essential for the switching of the growth cone response to the axon guidance cues (Chen et al., 2008), indicating that spatiotemporally regulated protein isoform expression and diversity is crucial to establish proper neuronal networks.
Alternative splicing and alternative polyadenylation can produce multiple messenger RNAs (mRNAs) possessing distinct coding and regulatory sequences from a single gene. The regulated processes that generate such mRNA diversity are orchestrated by RNA-binding proteins (RBPs). In the nervous system, alternative splicing has many important roles, including controlling the spatial and temporal expression of protein isoforms that are necessary for neurodevelopment and the modification of synaptic plasticity (Li et al., 2007; Licatalosi et al., 2008; Ule and Darnell, 2006). Significantly, human genetic studies have indicated that RNA misregulation resulting from defects in RBP expression and function are linked to numerous diseases, including Fragile X syndrome, spinal muscular atrophy, spinocerebellar ataxias, motoneuron disease and others (Cooper et al., 2009; Darnell, 2010; Lukong et al., 2008).
NOVA1 and NOVA2, RBPs initially identified as targets in autoimmune motor neuron disease (Buckanovich et al., 1993; Darnell and Posner, 2003), are RNA-binding splicing regulators specifically expressed in neurons in the CNS (Buckanovich et al., 1996; Yang et al., 1998). NOVA1 and NOVA2 harbor three KH-type RNA binding domains, and in vitro RNA selection (Buckanovich and Darnell, 1997; Jensen et al., 2000b; Yang et al., 1998) and X-ray crystallography (Lewis et al., 2000; Teplova et al., 2011) demonstrate that NOVA1 and NOVA2 bind directly to RNA sequences harboring YCAY motifs. The development of RNA:RBP crosslinking and immunoprecipitation (CLIP) (Ule et al., 2003; 2005a) followed by high-throughput sequencing (HITS-CLIP) (Licatalosi et al., 2008) methods has enabled identification of genome-wide RBP:RNA interaction maps in vivo. The bioinformatic integration of HITS-CLIP binding data and functional outcome revealed by exon junction microarray and RNA-seq data sets led to the conclusion that NOVA regulates alternative splicing following the rule that the position of NOVA binding to pre-mRNA predicts its action to enhance or inhibit alternate exon inclusion (Licatalosi et al., 2008; Ule et al., 2005b; 2006; Wu et al., 2013, Zhang et al., 2010). These analyses revealed ~700 NOVA1/NOVA2 target alternate exons and allowed to predict the biological process regulated by NOVA1/NOVA2. However, these studies have found no biochemical actions that were unique to NOVA1 versus NOVA2 paralogues, reflecting their nearly identical KH-type RNA binding domains.
In the present work, we find that NOVA2 uniquely regulates a series of alternative splicing events of axon guidance related genes, including deleted in colorectal carcinoma (Dcc), Roundabout, Axon Guidance Receptor, Homolog 2 (Robo2), Slit homolog 2 (Slit2), and EPH Receptor A5 (Epha5). These findings derived from combining transcriptome-wide NOVA1 and NOVA2 specific HITS-CLIP with RNA-seq analysis in the cortex of mice lacking either Nova1 or Nova2. Nova2 deficiency results in some common phenotypes, such as severe growth retardation, as previously reported in Nova1-/- mice (Jensen et al., 2000a), as well as unique phenotypes, including agenesis of the corpus callosum (ACC) in Nova2-/- mice. Unexpectedly, given that motoneurons express both NOVA1 and NOVA2, a subset of motoneuron axons directed to the ventral diaphragm show outgrowth defects in Nova2-/- but not Nova1-/- mice. In addition, the efferent innervation to the cochlea which is derived during development from a ventral subpopulation of facial motoneurons, is normal in Nova1-/- mice, reduced in Nova2-/- mice, and completely stalled in the absence of both Nova1 and Nova2, suggesting that alternative splicing regulator NOVA2 plays critical roles in axon pathfinding context in mammals and that NOVA1 may have a cooperative role in this process. Based on these observations from genome-wide and histological analyses of Nova2-/- mice, we conclude that NOVA2-mediated RNA regulation is essential for CNS development by regulating neural networks wiring.
Results
Generation of Nova2 null mice and growth retardation in Nova2 null mice
We utilized targeted mutagenesis to disrupt Nova2 genetic function in mice. A targeting cassette was designed to replace the first exon of Nova2 gene (Figure 1A). The genotype of the resulting progenies was confirmed both by Southern blot (Figure 1B) and PCR (Figure 1C). The progeny displayed the expected Mendelian ratio of mice for the heterozygous and homozygous Nova2 mutation. To confirm that homozygous Nova2 mutant mice did not express NOVA2 protein, Western blot analysis was carried out on protein extracts from P0 cortex (Figure 1D) of wild-type, Nova2-/-, and Nova1-/- mice, using anti-NOVA1, anti-NOVA2, and antisera from a patient with paraneoplastic opsoclonus-myoclonus ataxia (POMA), which recognizes all NOVA protein species (Yang et al., 1998). Expression of NOVA2 protein isoforms was absent in Nova2-/-mice (Figure 1D).
Interestingly, more than one protein isoform was absent in protein extracts from Nova2 null mice: a single 50–52 kDa NOVA2 protein species corresponding to the predicted size of the previously described Nova2 ORF, and a series of previously described (Ule et al., 2005a; Eom et al., 2013) but uncharacterized protein isoforms of ~70 kDa recognized by anti-NOVA2 antibody and POMA antisera, which we named SuperNOVA2. We found three putative initiator ATG sequences upstream of the known Nova2 start codon on the 5’ UTR of the transcript. The mobility of SuperNOVA2 on electrophoresis was comparable with the protein product expressed from SuperNOVA2 plasmid containing full-length sequence of the 5’ UTR (Figure 1E), confirming that SuperNOVA2 translation started from an upstream codon in the Nova2 5’ UTR. Using anti-NOVA1 or SuperNOVA2/NOVA2 specific antibodies, we found that NOVA1 protein levels were significantly increased, by ~40% in Nova2-/- mice cortex (Figure 1D) but NOVA1 anatomic expression pattern was similar to wild-type mice sections (data not shown). Thus NOVA1 has partially overlapping distributions with SuperNOVA2/NOVA2 (hereinafter we defined SuperNOVA2/NOVA2 as NOVA2), is upregulated by unknown means in the absence of Nova2, and may have partially redundant functions with NOVA2 in the cortex.
Nova2 null animals were born indistinguishable from littermates but failed to thrive, demonstrating progressive motor dysfunction and overt motor weakness, and they died an average of 14–18 days after birth (Figure 1—figure supplement 1). There was no apparent phenotype in Nova2 heterozygotes, although we previously found that when assessed by electroencephalography, 6 month old Nova2+/- (and Nova1+/-) mice had frequent synchronous cortical interictal discharges consistent with epilepsy (Eom et al., 2013). In conclusion, we demonstrate that the Nova2 gene encodes for multiple protein isoforms, including an identified ~70 kDa SuperNOVA2 isoform, and that Nova2 is necessary for post-natal survival.
Genome-wide NOVA1- and NOVA2-RNA interaction maps
The NOVA1/2 splicing-regulatory network and its target alternative splicing events have been previously defined by utilizing NOVA2 and pan-NOVA HITS-CLIP data, splicing sensitive microarray data, NOVA1/2-binding motifs, and integrative modeling (Licatalosi et al., 2008; Ule et al., 2005b; 2006; Zhang et al., 2010), yet these data sets were integrated data sets from multiple developmental stages, CNS regions, and Nova1/Nova2 genotypes. Relatively little is known about the biologic consequences of NOVA1 and NOVA2 action and about the functional differences between NOVA1 and NOVA2 in a specified CNS region and at a particular developmental stage.
To identify NOVA1 and NOVA2 direct RNA targets in E18.5 cortex, a transcriptome-wide library of NOVA1-RNA and NOVA2-RNA interactions were generated by HITS-CLIP from E18.5 wild-type mice cortex and sequenced. NOVA1 and NOVA2 CLIP reads were filtered and aligned to the mouse genome (mm9). A total of 2,318,553 unique NOVA2 CLIP tags were obtained from three biological replicates (483,446, 836,678, and 998,429 unique tags, respectively), with a total of 139,007 clusters from at least two or three biologic replicates (defined as biological complexity of two (BC2)). A total of 910,342 unique NOVA1 CLIP tags were obtained from three biological replicates (319,904, 292,709, and 297,729 unique tags, respectively), with a total of 54,546 clusters (BC2). To identify sequence features most commonly associated with either NOVA1-RNA or NOVA2-RNA interactions, 4 nucleotide (nt) long motifs were counted in BC2 NOVA1 and NOVA2 clusters and in a set of shuffled control sequences as a control. Each of the top five ranking tetramers in NOVA1 and NOVA2 clusters were the same (UCAU, CAUC, CAUU, UCAC, and UUCA) (Figure 2A) and coincided with the previously identified YCAY NOVA binding motif. Enrichment of the YCAY motif around NOVA1 and NOVA2 clusters was comparable between NOVA1 and NOVA2 and greatly enriched within NOVA1 and NOVA2 CLIP clusters (Figure 2B), suggesting that NOVA1 and NOVA2 recognized the same YCAY RNA sequence. NOVA1 and NOVA2 CLIP tag positions on known NOVA regulated alternatively spliced exons (Zhang et al., 2010) were comparable (Figure 2C and D). For example, alternative splicing of Agrn minigene containing exon 31–34 were regulated by NOVA1 and NOVA2 in the same splicing direction, in this case mediating the inclusion of the exon (Figure 2—figure supplement 1), indicating that NOVA1 and NOVA2 have similar RNA interaction profiles on NOVA1/2 target alternative splicing sites and are capable of regulating the same alternative splicing events.
Twenty-six percent or eighteen percent, respectively, of the NOVA2 or NOVA1 BC2 clusters consisted of a minimum of 10 tags (defined as peak height (PH); NOVA2; 36,591 of 139,007, NOVA1; 9988 of 54,546). 80.1% of NOVA1 clusters were shared with NOVA2 cluster positions, while 22.1% of NOVA2 clusters overlapped with NOVA1 (Figure 2E). Interestingly, the genomic distribution of NOVA2 clusters and NOVA1 showed different distributions within their target transcripts (Figure 2F and Figure 2—figure supplement 2). 77.7% of NOVA2 clusters were located in introns and 13.9% in exons (5’ UTR, CDS, and 3’ UTR) compared to 45.4% of NOVA1 clusters that were in introns and 45.1% in exons (Figure 2F). 84.9% of NOVA2 specific BC2 clusters (PH>=10) were enriched in introns, while 74.1% of NOVA1 specific BC2 clusters (PH>=10) were in exons (Figure 2—figure supplement 2), indicating that NOVA2 binds preferentially to introns than exons, suggesting that NOVA2 may play a greater nuclear role than NOVA1, and demonstrating that RNA-interaction profiles on a genome-wide scale are different between NOVA homologues.
The distribution of NOVA2 throughout the brain mirrored previous immunohistochemical and in situ hybridization data (Yang et al., 1998) showed that NOVA2 was expressed at high levels in cortex and hippocampus, and at lower levels in midbrain and spinal cord, where NOVA1 was expressed at high levels in a generally reciprocal fashion, with low levels in the cortex and relatively high levels in the midbrain and spinal cord (Figure 2—figure supplement 3). The NOVA2 expressed cell in the cortical plate of neocortex was ubiquitously distributed at comparable expression level, yet NOVA1 was expressed in the specified cell types. Taken together, the HITS-CLIP and immunohistochemical data suggest that NOVA1 and NOVA2 perform unique biological functions in different brain areas and cell types, and that in those few cortical neurons expressing NOVA1, NOVA2 might be expected to have some redundant activity, while the reciprocal may not so often be the case.
NOVA2-dependent alternative splicing changes of axon guidance related genes
To identify NOVA1 or NOVA2 target transcripts whose abundance or alternative splicing was changed in either Nova1-/- or Nova2-/- cortex, paired-end 125 nucleotide RNA-seq libraries were prepared from littermate E18.5 wild-type, Nova1-/-, and Nova2-/- mice cortex. A total of 225,054,059 and 213,596,515 unique reads were obtained from three biological replicates of wild-type and Nova1-/- mice, respectively, and a total of 190,282,241 and 180,200,564 unique reads were obtained from three biological replicates of wild-type and Nova2-/- mice, respectively. These experiments showed that the abundance of 60 transcripts were significantly changed in Nova2-/- mice compared to only 2 transcripts which showed steady-state changes in Nova1-/- mice compared to wild-type littermates (FDR<0.05 and |log2FC|>=1) (Figure 3A), consistent with lower levels of NOVA1 expression in cortex. More changes were seen in analyzing splicing-dependent changes, and this pattern of differential NOVA2 and NOVA1 effects was also seen. We integrated NOVA1 and NOVA2 specific HITS-CLIP data, RNA-seq data, and NOVA binding YCAY motifs to define 1991 NOVA2 and 58 NOVA1-mediated alternative splicing events in the mouse cortex (Figure 3B and Figure 3—figure supplement 1). Interestingly, there was no correlation between delta I (ΔI) of wild-type versus Nova1-/- and wild-type versus Nova2-/- (R = 0.09) (Figure 3C). Taken together, this data shows that Nova1 and Nova2 have different RNA regulatory networks in developing cortical neurons; while the quantitative differences likely relate to the greater and more unique cellular expression of NOVA2 relative to NOVA1 in cortex, they also support the suggestion that the two NOVA homologues may regulate unique biological processes.
-
Figure 3—source data 1
- https://doi.org/10.7554/eLife.14371.011
-
Figure 3—source data 2
- https://doi.org/10.7554/eLife.14371.012
We explored the biologic pathways regulated by NOVA1 and NOVA2 with gene ontology (GO) analysis on the transcripts harboring NOVA-regulated alternative exons. Among the KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway terms associated with NOVA2-regulated genes (FDR<0.05) were axon guidance signaling (FDR=0.0039) and adherens junctions signaling (FDR = 0.017) (Figure 3—figure supplement 2 and 3). In the GO terms associated with NOVA2-regulated genes (FDR<0.05), axonal projection related terms (cell morphogenesis, neuron projection morphogenesis, axonogenesis, cell morphogenesis involved in neuron differentiation, and cell projection morphogenesis) were also enriched (Figure 3—source data 1). NOVA2-dependent alternative splicing events of axon guidance related genes were validated by semi-quantitative RT-PCR with RNA prepared from E18.5 wild-type Nova2-/-, and Nova1-/- cortex (Figure 4 and Figure 4—figure supplement 1). The alternative splicing events (Dcc exon 17, Slit2 exon 28b, Robo2 exons 6b and 21, Epha5 exon 7, Arhgef12 exon 4, Ppp3cb exon 10b, Neo1 exon 26, and Rock1 exon 27b) were significantly changed in Nova2-/- but not in Nova1-/- mice cortex, coinciding with RNA-seq results. No association was detected in the NOVA1-regulated alternative exons GO or KEGG pathways terms, likely due to the smaller number of NOVA1-regulated transcripts in cortex.
RNA-seq analysis in the E18.5 midbrain and hindbrain of wild-type and Nova1-/- mice, where NOVA1 is more abundantly expressed, identified that 119 of 1991 NOVA2-dependent alternative splicing events were changed in Nova1-/- mice (|ΔI|>=0.1, FDR<0.1) (Figure 3—figure supplement 4A). There was no correlation between ΔI of wild-type versus Nova1-/- (mid- and hindbrain) and wild-type versus Nova2-/- (cortex) (R=0.20) (Figure 3—figure supplement 4B) as well as identified between ΔI of wild-type versus Nova1-/- (cortex) and wild-type versus Nova2-/- (cortex). In the midbrain and hindbrain of Nova1-/- mice, only one alternative splicing event (Robo2 exon 6b) was significantly changed in the NOVA2-regulated genes list of KEGG axon guidance pathway but it was smaller change when compared with ΔI of wild-type versus Nova2-/- (cortex) (Figure 3—figure supplement 4C). Taken together, these data suggest that NOVA2 regulates a distinct set of transcripts related, in particular, to axon guidance, in the E18.5 mouse cortex.
NOVA2-mediated alternative splicing switch of axon guidance genes during mouse cortical development
To pursue the potential connection between NOVA2 and axonal guidance, we examined NOVA2 expression during cortical mouse brain development. We found that NOVA2 expression level detected by immunofluorescence was high in the cortical plate (CP) and subplate (SP) that consists of post-mitotic neurons at E18.5 (Figure 2—figure supplement 3). Similarly a high-expression of NOVA2 in CP and SP was also observed at E12.5, E14.5, and E16.5 (Figure 5A), indicating that cortical NOVA2 expression level was progressively up-regulated during neural differentiation from neural progenitor cells (NPC).
To assess whether NOVA2 regulated alternative splicing is also regulated in a developmentally regulated manner, we compared expression of alternatively spliced NOVA2 target exons in E12.5 wild-type, E18.5 wild-type, and E18.5 Nova2-/- cortex (Figure 5B). This revealed that NOVA2-target splicing events were developmentally regulated between E12.5 and E18.5, such that splicing in E18.5 Nova2-/- cortex reverted to patterns seen earlier in E12.5 wild-type cortex. Based on these observations, we conclude that NOVA2 switches the alternative splicing patterns of a series of axon guidance genes during mouse cortical development.
Axon guidance defects in the brain of Nova2 but not Nova1 null mice
The observations that NOVA2 regulates alternative splicing events in transcripts encoding axon guidance signaling factors led us to test whether axon guidance itself is affected in Nova2 null mice. When coronal and horizontal sections of E18.5 mice brains were double-immunostained for L1, an axonal marker, and TAG-1, a corticofugal axon marker, we discovered that Nova2-/- mice, but not wild-type littermates, had agenesis of the corpus callosum (ACC) (14/14 mice) (Figure 6 and Figure 6—figure supplement 1A). ACC was in all Nova2-/- mice serial sections examined (Figure 6—figure supplement 1B), was seen along with Probst bundles (abnormal collections of cells characteristically seen in patients with ACC; arrows in Figure 6A–o and Figure 6—figure supplement 1A–l), and was also clearly evident in sections stained with the pioneer axon marker Neuropilin-1 (Figure 6—figure supplement 1C). Interestingly, Nova1-/- mice showed normal CC formation (Figure 6—figure supplement 1D). The anterior commissure axons were observed in both Nova2-/- and Nova1-/- mice (Figure 6—figure supplement 2).
To independently corroborate these observations, corpus callosal commissural axons were visualized with DiI anterograde tracer placed into the cingulate cortex. No CC commissural axons were detected that crossed the midline in Nova2-/- brain, however, were clearly detected in anterior and posterior sections of wild-type brain at P0 (Figure 6B). Corpus callosal commissural axons from the lateral neocortex also failed to cross the midline in Nova2-/- mice (Figure 6—figure supplement 1E). Interestingly, corticothalamic axons of both wild-type and Nova2-/- mice terminated in a similar thalamic region at E18.5 (arrows in Figure 6B), indicating that only specific aspects of axon guidance are disrupted in Nova2 null mice. These observations suggest that Nova2 but not Nova1 regulate a series of alternative splicing events that are necessary for forming the CC, and demonstrate that NOVA2 and NOVA1 have different functions in CC axon guidance pathways.
To understand the axon guidance phenotype in more detail, we compared the pathway of L1 positive axons traversing the upper cortical layer in the cortex of Nova2-/- and wild-type mice at E16.5 and E18.5 (E16.5; Figure 7 and E18.5; Figure 6A). In wild-type mice, L1 positive axons passing between the subplate and subventricular zone were detected as one wide bundle, whereas in Nova2-/- mice L1 positive axon routes were separated into two bundles (Figure 7A). The L1 positive axons detected in the deeper cortical layer in Nova2-/- appeared normal, in that they passed through similar cortical layers as did axons in wild-type mice. In contrast, separated axons passing through the upper L1 positive axonal pathway, which were Netrin-G1 (NTNG1) positive axons, were only seen in Nova2-/- mice and were passed along a subplate path enriched in NURR1 (a subplate marker), as revealed by L1/NTNG1 or L1/NURR1 double-immunostaining (Figure 7B and Figure 7C). These results indicate that only a portion of the L1 positive axonal path, which is NTNG1 positive, is specifically affected by the absence of Nova2. We conclude that NOVA2 plays a role as a modifier for a unique set of neuron-specific axons in the developing cortex.
Motor innervation defects in the ventral diaphragm in NOVA2 dependent manner
We previously reported that NOVA1/NOVA2 regulate the alternative splicing of the agrin Z exons (Ruggiu et al., 2009). Mice that lack both Nova family members fail to cluster AChR at neuromuscular junctions, including those of the diaphragm. Given the role for NOVA2 in axon guidance in the developing cortex (Figure 6, Figure 6—figure supplement 1, and Figure 7), we revisited the role of NOVA1 and NOVA2 on motoneuron axon guidance. The axons of peripheral motor neurons innervating the diaphragm muscle were visualized by immunostaining of neurofilaments at E18.5 (Figure 8A) and at E14.5 and E16.5 (Figure 8B) and the innervation percentage of muscle in each quadrant of diaphragm at E18.5 was quantified (Figure 8C). Although the motor innervation patterns and percentages of the dorsal diaphragm were normal in Nova2-/- and Nova1/Nova2-double knockout (dKO), those of the ventral diaphragm were severely abnormal with approximately 60–85% of the muscle uninnervated in either Nova2-/- and Nova1/Nova2 dKO mice (p<0.001). Similar abnormalities were observed in Nova2-/- and Nova1/Nova2 dKO mice at E14.5 and E16.5 (Figure 8B). These innervation defects were specific to Nova2-deficient mice, as defects in ventral diaphragmatic innervation were never observed in wild-type or Nova1-/- mice. We cannot rule out a role of NOVA1 in the axon guidance of a subset of motoneurons given the slight decrease in ventral innervation in the Nova1/2 double-KO compared to the Nova2-/- mice. Nevertheless, taken together, these data demonstrate that in addition to what was observed in CC axons, Nova2 is also essential for phrenic nerve innervation of the ventral but not dorsal diaphragm.
Auditory efferent innervation defect in NOVA2 dependent manner
The inner ear is innervated by two types of fibers, afferent innervation consisting of both vestibular and cochlear (spiral ganglion neurons; Mao et al., 2014) and from efferent fibers derived from rhombomere 4 of the hindbrain (Simmons et al., 2011). Efferents are guided by and extend along afferents to reach the sensory epithelia of the ear arriving at approximately the same time (Fritzsch et al., 1998; Ma et al., 2000). Efferents are a unique population of ventral brainstem neurons derived from facial motoneurons that project to the inner ear in contrast to facial muscle fibers or glands (Karis et al. 2001). Some transcription factors have been hypothesized to be relevant for guiding inner ear efferents (Duncan and Fritzsch, 2013), however, essentially nothing is known about how these neurons diverge from facial motoneurons to selectively reach the ear to innervate hair cells and afferent processes (Simmons et al., 2011). Dye tracing analysis comparing with control mice at E14.5 (Figure 9A,A’) showed that Nova2 (Figure 9C,C’) but not Nova1 (Figure 9B,B’) was essential for proper progression of efferent growth along afferents. Moreover, although efferent neurons expressed both NOVA1 (Figure 9—figure supplement 1I–J) and NOVA2 (Figure 9—figure supplement 1F–H), Nova1 by itself was not capable of maintaining the function. In the absence of both Nova1/Nova2, the efferents stalled when they reached vestibular ganglion neurons (Figure 9D,D’) while afferent fibers that also expressed both NOVA1 and NOVA2 (Figure 9—figure supplement 1A–E, K) reached the cochlea normally. Our data suggested that Nova2 played a crucial role for efferent guidance relative to Nova1, but both cooperate for normal efferent fiber extension. Interestingly, vestibular efferent innervation which segregated from cochlear efferents at E14.5 (Bruce et al., 1997) was unaffected by the absence of Nova1 or Nova2 alone (see green efferents reaching the posterior canal crista (PC) in Figure 9B,C) but was also completely stalled in Nova1/Nova2 double knockouts (green label, Figure 9D). Postnatal Nova2-/- mice vestibular innervation was also normal (not shown).
The cochlea has a tonotopic map where the sensory hair cells at each position along the organ are most sensitive to a particular frequency, forming a continuous gradient from high frequency in the base to low frequency in the apex (Howard et al., 1988; Hudspeth, 1989; Liberman, 1982). Dye tracing analysis showed that the loss of efferent innervation was frequency specific; it was almost completely lost in the base of the Nova2-/- cochlea and only partially lost towards the apex (compare Figure 9A’ and C’ ). A more detailed analysis was possible at E18.5 (Figure 9E–L), when the medial olivocochlear (MOC) efferent axons expand through the tunnel of Corti, prior to their reaching the outer hair cells (OHC) at P2 and establishing functional synapses at P7 (Simmons et al., 2011). These axons stained strongly with the Thy-1, YFPJ line that also has high expression in motor neurons, demonstrating the motor origin of auditory efferents. Thy-1 positive axons reached the inner hair cells (IHC) at E18.5 in control (Figure 9E,I) and Nova1-/- (Figure 9F,J) mice but this innervation was reduced in the Nova2-/- (Figure 9G,K) mice, and this was especially evident in the base (Figure 9K). In the absence of both Nova1/Nova2 there were no labeled axons reaching the cochlea (Figure 9H,L), showing that efferents remain stalled. Imaging of efferent innervation at the onset of hearing (P14) showed a clear reduction in axons reaching IHC and specially OHC in the base when comparing Nova2-/- mice (Figure 10C,D,G,H) to controls (Figure 10A,B,E,F). At this age we also found a reduction in efferent synapse formation as demonstrated by reduced labeling with the pre-synaptic marker synaptophysin (blue label at the base of IHC and OHC). Whole mount apical cochlea (Segment 2, see map in Figure 10O) at P10 showed a reduction in axons reaching IHC and an almost complete absence of axons reaching OHC (Figure 10K,L) compared to controls (Figure 10I,J).
These defects were also evident in apical cholinergic innervation (Dalian et al., 2001; Maison et al., 2003) after the onset of hearing (P15-P25; Figure 10—figure supplement 1A–C). At these stages, we were able to demonstrate that the loss of axonal innervation was functionally significant, as assessed by auditory brainstem responses (ABRs) and distortion products otoacustic emissions (DPOAEs) recordings (Figure 10M–N) In addition, an apparent increase in afferent innervation was found in Nova2-/- mice by neurofilament immunodetection (Figure 10A–L and Figure 10—figure supplement 1D–E). The difference was significant in the apex as determined by quantitation of IHC ribbon synapses (Meyer et al., 2009) (Figure 10—figure supplement 1F–J). Taken together, these data demonstrate a role for NOVA in cochlear innervation, efferent pathfinding, and normal hearing function in vivo.
Discussion
In the present work, we provide evidence for a novel and unique role for the tissue specific splicing regulator NOVA2 as an axon pathfinding modifier in cortical CC axons, motoneuron, and auditory efferents. In prior studies, integration of transcriptome-wide HITS-CLIP analysis, exon junction microarray data sets, and bioinformatics led to the identification of a large set of NOVA1/2 RNA functional interaction sites and major biological pathways which NOVA1/2 orchestrates in vivo (Licatalosi et al., 2008; Ule et al., 2005b; Ule and Darnell, 2006; Zhang et al., 2010). Given identical binding characteristics in shared transcripts (e.g. Agrin, others in Figure 2), we infer the following: 1) specificity is determined primarily by biochemical interactions: kD of NOVA KH domains binding to accessible YCAY motifs; this might include differential levels of NOVA1/NOVA2 expression within a cell; 2) differential binding is a function of differences in cell expression and/or differences in cell biology (e.g. subcellular distribution of NOVA1 and NOVA2). Evidence for differences in cell expression is clear in the cortex. Differences in actions within cells such as motoneurons that express both NOVA1 and NOVA2 isoforms could arise from different levels of NOVA1 and NOVA2 within the neurons, or due to differential localization of NOVA1 and NOVA2 within the motoneuron, as suggested by the differences in CLIP binding (intron versus 3’ UTR) of the two proteins (Figure 2F). The development of single cell-type CLIP would be able to address this issue. Prior genomic and biochemical studies have found no actions that were unique to either NOVA1 or NOVA2 paralogues. The present study identifies NOVA2 as playing a unique role in axonal pathfinding in embryonic development by combining histological analysis with genome wide NOVA1 and NOVA2 specific HITS-CLIP and RNA-seq analysis in Nova1-/- and Nova2-/- mice. These data build upon previously identified roles for NOVA1 and NOVA2 in early development as modifiers of neuronal migration in late-generated cortical neurons and Purkinje neurons through control of splicing of Dab1 exon 7b/7c (Yano et al., 2010), for Nova1/Nova2 as essential factors for AChR clustering at the NMJ through the regulation of Agrin Z exons (Ruggiu et al., 2009), and for Nova2 in the induction of long term potentiation of slow inhibitory, but not excitatory, post synaptic potential currents after coincidence detection in CA1 hippocampal neurons (Huang et al., 2005). The current study provides the first evidence that the formation of corpus callosum in brain, proper innervation of ventral diaphragm, and efferent innervation to the cochlea at developing embryonic stage all require specific actions mediated by NOVA2 (Figure 11).
Immunohistochemical and DiI tracer analysis showed that Nova2 deficiency results in severe ACC (Figure 6, and Figure 6—figure supplement 1). There are over 50 mouse genes known to be required for the formation of the corpus callosum (reviewed in Paul et al., 2007; Edwards et al., 2014). NOVA2 HITS-CLIP and RNA-seq analysis revealed that 11 NOVA2 target alternative splicing events on 9 genes (Dcc, Robo2, Epha5, Slit2, Ank2, Dclk1, Enah, Cask, Ptprs) intersected with the genes associated with mouse ACC phenotypes and human ACC syndromes (Figure 6—figure supplement 3), suggesting that the ACC phenotype in Nova2-/- mice may be caused by splicing dysregulation of multiple coordinately NOVA2-regulated transcripts, consistent with previous observations that NOVA proteins regulate subsets of transcripts that mediate coordinate biologic functions (Ule et al., 2005b). Given that Netrin-1 null mice and Dcc null mice shows similar ACC phenotype to Nova2 null mice (Fazeli et al., 1997; Serafini et al., 1996), we tried to rescue the ACC phenotype with the DCC-long isoform, that is decreased in the cortex of Nova2-/- mice, by in utero electroporation (data not shown). The expression of DCC-long (long isoform of exon 17) isoform in cingulate cortex at E14.5 was not sufficient to form the commissure corpus callosum, indicating that the DCC-long isoform was not sufficient to rescue the ACC phenotype, and suggesting the possibility that the ACC phenotype in Nova2-/- mice either does not involve the DCC-long isoform or is caused by the dysregulation of multiple transcriptomes.
We found that in the cortex, NOVA2 regulates alternative splicing events of two netrin receptors: Dcc exon 17 and Neo1 exon 27 in developmentally regulated manner (Figure 4, Figure 4—figure supplement 1, and Figure 5). In the spinal cord, Dcc exon 17 is regulated by NOVA1/2 and required for the spinal commissural neuron development (see accompanying paper by Leggere et al.). Recent structure analysis of netrin-1/NEO1 and netrin-1/DCC-short (short isoform of exon 17) revealed that netrin-1/NEO1 form a 2:2 heterotetrameric netrin-1/NEO1 complex and netrin-1/DCC-short have a continuous netrin-1/DCC assembly (Xu et al., 2014). Netrin-1 recognizes the fourth and fifth fibronectin type III domains (termed FN4 and FN5, respectively) of both NEO1 and DCC. Interestingly, the NOVA2-regulated Dcc exon 17 is located in the linker region between the FN4 and FN5 domains, and this splicing switch from DCC-short to DCC-long isoform is enough to support the architectural switch of netrin-1/DCC complex from the continuous assembly to the 2:2 assembly (Xu et al., 2014). This suggests that NOVA2-dependent splicing regulation results in architecturally distinct netrin-1/DCC complexes that elicit distinct downstream signaling. Changes in NOVA2-dependent regulation of Dcc exon 17 splicing during development (Figure 5) further indicate that NOVA2-regulated splicing controls the ability of cortical and spinal neurons to utilize distinct netrin-1/DCC assemblies to respond to a variety of cellular demands.
NEO1, another member of the Dcc family, shares ~50% amino acid identity with DCC (Vielmetter et al., 1994) and has four known alternative splicing sites. The NEO1 extracellular region is a receptor for repulsive guidance molecule (RGM) family molecules in addition to netrins. Upon ligand binding to the extracellular domain of NEO1 and DCC, their intracellular domains can activate multiple downstream signal transduction pathways, that modulate chemotropic axon guidance (Bradford et al., 2009). The inclusion of Neo1 exon 27 inserts a 53 amino acid sequence into the intracellular region. Inclusion is increased in E18.5 Nova2-/- cortex, suggesting that NOVA2-dependent alternative splicing of Neo1 exon 27 regulates downstream signal transduction. It is well-established that DCC receptors mediate chemo-attraction upon netrin-1 binding, while the UNC5 receptor alone or heterodimerized with DCC elicits chemo-repulsive activity in neuronal tissues (Hong et al., 1999). Interestingly, UNC5B interacts with NEO1 as a co-receptor for RGMa and associates with leukemia-associated guanine nucleotide exchange factor (LARG, also termed as ARHGEF12) to mediate the signal transduction leading to growth cone collapse (Hata et al., 2009). Although the abundance of Arhgef12 RNA in cortex was comparable between E18.5 wild-type and Nova2-/- mice (log2FC = 0.10, FDR = 0.78), Arhgef12 exon 4 inclusion was significantly increased in Nova2-/- mice (Figure 4—figure supplement 1), suggesting the possibility that NOVA2 modulates RGM signal transduction mediated by NEO1, UNC5b, and ARHGEF12 through coordinate regulation of alternative splicing of those RNAs.
Interestingly, exon 7–8 deletion of Srrm4/nSR100 gene, an alternative splicing regulator, mice shows partial ACC phenotype and neurite outgrowth defect in diaphragm (Quesnel-Vallieres et al., 2015) in addition to hearing defects (see below), indicating that alternative splicing regulation mediated by RBPs (e.g. NOVA2, SRRM4/nSR100) is required for the formation of corpus callosum and neurite outgrowth into diaphragm.
Our analyses were focused on axon guidance related genes but do not rule out the possibilities that axonal pathfinding defect phenotypes in Nova2-/- were caused by deficits of multiple biological processes. Several relevant processes were enriched in GO and KEGG pathway analysis, including protein modification, phosphorylation, cellular protein metabolic process, synaptic transmission, chromatin modification, adherens junction, and cell morphogenesis involved in neuron differentiation. Our previous work has showed that the proliferation of cortical neural progenitor cells was not disrupted in Nova2-/- (Yano et al., 2010). RNA-seq data showed that the transcriptome abundance of the neuronal differentiation markers were comparative between wild-type and Nova2-/- mice (Figure 3—source data 2), indicating that the differentiation from neuronal progenitor to neuron is comparable between wild-type and Nova2-/- mice.
Despite the important role of alternative splicing in the formation of the cochlear tonotopic gradient (Navaratnam et al., 1997; Rosenblatt et al., 1997; Miranda-Rottmann et al., 2010) and in hair cell development and function (Kollmar et al., 1997; Liu et al., 2007; Webb et al., 2011) to date only two alternative splicing factors have been implicated in hearing: Srrm4/nSR100 (Nakano et al., 2012) and Sfswap (Moayedi et al., 2014). Both are expressed in hair cells necessary for the development of the sensory epithelium, in contrast with Nova1 and Nova2, which are the first neuron-specific splicing factors described to regulate cochlear innervation.
Interestingly auditory efferents are a ventral population of rhombomere 4 facial brachial motoneurons that segregate during development (Karis et al., 2001) and while facial innervation is normal, the ventral auditory efferents are selectively affected in the Nova2-/- mice (Figure 9A–D and data not shown). This specificity is similar to the observation that dorsal motoneuron innervation of the diaphragm is normal while ventral innervation (Figure 8) is affected in the Nova2-/- mouse, suggesting a more general ventral specialization of axonal guidance regulation by NOVA2. These observations are reminiscent of those seen with defects in GATA-3 innervation of facial motor neurons (Karis et al., 2001) and dorsal/ventral pathfinding of limb motoneurons regulated by EphA4 (Helmbacher et al., 2000), suggesting the possibility of an integrated role for NOVA2 in dorsal/ventral axonal pathfinding.
To assess OHC function we measured DPOAEs (Lonsbury-Martin and Martin 1990) in Nova2-/- mice. There was a significant reduction in DPOAEs in the apex of Nova2-/- mice cochlea (Up to 10 dB SPL, Figure 10N) where a significant reduction in cholinergic innervation was also found (Figure 10—figure supplement 1A–C). Nevertheless the reduction in efferent innervation is comparatively more severe in the base of the cochlea (Figure 9 and 10A–H) suggesting that OHCs sensitive to low frequency sounds are more dependent on Nova2 controlled efferent innervation. Direct testing of efferent function by electrical activation of these nerves has been shown to reduce the baseline DPOAE intensity (Liberman et al., 1996; Vetter et al., 1999). Preliminary testing in Nova2-/- mice showed that this inhibitory function was not affected (Stéphane S Maison, personal communication), suggesting that the residual efferent innervation was able to sustain the function of further reducing the DPOAEs.
Only an indirect mouse model of efferent innervation defect has been available to date: a mutation in the α9α10 hair cell nicotinic receptor (Taranda et al., 2009) in which the DPOAEs threshold levels are reduced by approximately 7–14 dB SPL (a reduction in 10 dB SPL is perceived approximately as half the intensity). In addition, the reduction in DPOAEs found in the Nova2-/- mice could explain the ABR threshold elevation at the same frequency range (Figure 10M) (Sun and Kim 1999; Taranda et al., 2009). The opposite possibility, that the defect evidenced by ABRs is responsible for the reduction in DPOAEs is unlikely since it has been reported that even with 80% loss of IHC (with intact OHC) the DPOAEs are only reduced by 10–20 dB SPL (Trautwein et al., 1996). It has also been recently shown that afferent innervation to the OHC has no effect on DPOAEs (Froud et al., 2015). Despite the reduction in ABRs which measure afferent pathway function, we found an increase in afferent innervation (Figure 10A–L) quantified by synaptic ribbon count in Nova2-/- mice (Figure 10—figure supplement 1E–J) although further study will be necessary to determine if this synapses are functional. Increase in afferent innervation could be a compensatory mechanism to decreased lateral olivocochlear (LOC) efferents, which form synapses on IHC afferents, or a direct effect of the loss of Nova2 expression in afferent neurons through an ephrin-dependent pathway (Defourny et al., 2013). An interesting possibility is that the transient synapses formed between MOC efferents and IHC during development (Simmons et al., 2011) have an important role in IHC development that is impaired in Nova2-/- mice.
Taken together, our observations suggest that protein/RNA diversity provided by NOVA2-mediated RNA regulation is required for proper axon pathfinding and formation of complex synapses/neural networks, particularly in dorsal/ventral choices, and that alternative splicing switches mediated by NOVA2 may regulate key developmental steps in mammalian biology and pathogenesis of neurological diseases.
Material and methods
Generation of Nova2 null mice
Request a detailed protocolTwo fragments from the Nova2 genomic locus of 2.2 and 6 kb, flanking a 1.5 kb DNA fragment harboring the exon containing the initiating methionine, were cloned into a targeting vector. An IRES-Cre-FRT-Neo-FRT targeting cassette was inserted into the Nova2 locus upstream of the initiator ATG in the first known coding exon of the Nova2 gene. Linearized targeting vector plasmid was electroporated into ES, and G418-resistant clones (150 ug/mL G418 for 24 hr) harboring homologous recombinants were screened by Southern blot with a 5’ Nova2 genomic probe and injected into mouse blastocysts to produce germline transmitted chimeras. The neomycin cassette was excised by breeding germline-transformed mice with transgenic animals expressing Flp recombinase under the control of the CMV promoter. The original Nova2+/- (and Nova1+/-) mice obtained in a B57B/6 background were backcrossed for 10 generations to CD-1 and FVB strains and to the thy1-YFP line J (YFPJ) strain (Jackson labs Stock 003709) which has high expression of YFP in motor neurons (Feng et al., 2000) or crossed to the R26R strain (Jackson lab stock 003309) for Nova2 tissue-wide expression analysis. Most experiments were done in CD-1 background, but mice in the inbred FVB background showed no phenotypic differences. FVB mice were used for the hearing tests reported here because they present normal hearing thresholds (Zheng et al., 1999) and less inter individual differences. All procedures in mice were performed in compliance with protocols approved by the Institutional Animal Care and Use Committee (IACUC) of the Rockefeller University or the Comité de déontologie de l’expérimentation sur les animaux (CDEA) of the Univeristy of Montreal.
Antibodies
Primary antibodies used for immunohistochemistory and western blot were as follows; goat anti-NOVA2 (C-16) (sc-10546, Santa Cruz), rabbit anti-NOVA1 [EPR13847] (ab183024, abcam), human anti-pan NOVA (anti-Nova paraneoplastic human serum), rabbit anti-PTBP2 (Polydorides et al., 2000), rat anti-L1 (MAB5272, Millipore), goat anti-TAG1/Contactin-2 (AF4439, R&D systems), goat anti-NetrinG1a (AF1166, R&D systems), rabbit anti-NURR1 (M-196) (sc-5568, Santa Cruz), rabbit anti-neurofilament (AB1981, Chemicon), mouse anti-NF200 (N5389, SIGMA), mouse anti-ctbp2 (612044, BD), rabbit anti-myoVI (sc-50461, Santa Cruz), and goat anti-neuropilin-1 (AF566, R&D systems). Anti-NOVA1 and anit-NOVA2 antibody specificity for immunohistochemistory and immunoprecipitation was confirmed (Figure 1—figure supplement 2).
Immunohistochemistory
Request a detailed protocolE16.5, E18.5 mice brains were fixed with 4% praraformaldehyde (PFA)/PBS at 4 degrees overnight, and sequentially replaced to 15% sucrose/PBS and 30% sucrose/PBS at 4 degrees for cryo-protection, then embedded in OCT compound. Frozen brains were sliced into 80 μm thick sections on a cryostat (CM3050S, LEICA). Brain sections and whole embryos fixed with 4% PFA/PBS were subjected to floating or whole-mount immunohistochemistry. These samples were washed three times with PBS at room temperature (RT), incubated with 0.2% Triton X-100/PBS at R.T., blocked with 1.5% normal donkey serum (NDS)/PBS at RT, and then incubated overnight at 4 degrees with primary antibodies or Alexa-dye conjugated phalloidin in 1.5% NDS/PBS followed by incubation with Alexa 488, 555, or 647 conjugated donkey secondary antibodies (1:1000) in 1.5% NDS/PBS. Images of immunostained specimens were collected by BZ-X700 (KEYENCE), fluoscence microscopes Axioplan 2 (Zeiss), or an inverted TCS SP8 laser scanning confocal microscope (LEICA) at The Rockefeller University Bio-Imaging Resource Center.
Cochlea tissue processing and imaging
Request a detailed protocolEmbryos were obtained by C-section of time-mated dams and fixed by overnight immersion in 4%PFA. Postnatal mice were intracardialy perfused with 4% PFA, cochleae immediately removed and perfused through the round window with a 1 ml syringe. After post-fixation 2:30 hr at 4°C cochleae were washed in PBS, directly dehydrated in 30% sucrose at 4°C or previously decalcified for 24–48 hr in 120 mM EDTA in PBS. For beta-galactosidase expression analysis the decalcified cochleae were perfused through the round window with M-1 embedding matrix (Thermo Scientific) and frozen in a block of the same matrix, 40 µm cryosections were incubated in staining solution (4 mM potassium ferricyanide, 4 mM potassium ferrocyanide, 1 mg/ml x-gal) dissolved in wash buffer (2 mM magnesium chloride, 0,01% sodium deoxycholate, 0,1% NP40 and 0.2 mM buffer phosphate pH 8.0) overnight at 37°C, washed in wash buffer 3 hr at 4°C, post-fixed again in 4% PFA and washed in PBS before mounting. For AChE staining the non-decalcified cochleae were dissected by carefully chipping the bone away and incubated in a solution containing 0.5 mg/ml acetylcholine iodide as explained elsewhere (Willott, 2001). For whole mount immunostainings the decalcified cochleae were dissected cutting into half turns (See Figure 10O), permeabilized and blocked in 0.1% TX100, 5% donkey serum in PBS. Primary and secondary antibodies were incubated sequentially 24 hr at 4°C. Washes were done in 0.05% TX100. For cryosections the decalcified cochleae perfused through the round window with M-1 embedding matrix (Thermo Scientific) or the fixed embryos were cut into 14–20 µm sections, permeabilized and blocked in 0.1% Triton-X100, 5% donkey serum in PBS. Primary antibodies were incubated 24 hr at 4°C and secondary antibodies for 2 hr at room temperature. Washes were done in PBS. Alexa labeled Phalloidin (1:100) was added to the secondary antibody mix. Images were collected using a Zeiss LSM 510 confocal microscope at The Rockefeller University Bio-Imaging Resource Center and a Zeiss LSM 700 confocal microscope at the IRIC core of the University of Montreal.
Axon tracing
Request a detailed protocolFor DiI trace analysis, DiIC18(3) crystal (D3911, molecular probes) were placed in cingulate cortex or lateral neocortex of P0 brain. After 8–21 days at 37 degrees in 4% PFA/PBS to allow dye diffusion, samples were cryo-protected with 30% sucrose/PBS at 4 degrees, embedded in OCT compound, and cut into 100 μm sections on cryostat. Counterstaining was with 1 μg/mL DAPI (4’,6-Diamidino-2-phenylindole dihydrochloride). The images of immunostained specimens were collected by a fluoscence microscope Axioplan 2 (Zeiss). Inner ear afferents and efferents were labeled with lipophylic (NeuroVue) dyes (Molecular Targeting Technologies; MTTI). These dyes were placed into specific areas for selective labeling of nerve fibers (Fritzsch et al., 2005; Duncan et al., 2015) Dyes were placed into the cerebellum and the alar plate of rhombomere 2 to label the afferent fibers to the ear. Dye was also placed into the basal plate of rhombomere 4 for efferent nerve labeling. The preparations were then incubated in 4% PFA at 36 degrees C for three days. After dye diffusion the ear was dissected out and placed on a glass slide with glycerol and a cover slip on top to be imaged with a Leica SP5 confocal microscope.
Western blot and semi-quantitative RT-PCR
Request a detailed protocolWestern blotting and RT-PCR analyses were performed as described elsewhere (Yano et al., 2010). For quantification, three biological replicates were used and quantified from obtained images by using ImageJ free software (http://rsb.info.nih.gov/ij/index.html).
Hearing tests
Request a detailed protocolABRs were recorded from deeply anesthetized mice. The positive needle electrode was inserted subdermally at the vertex, the negative electrode was placed beneath the pinna of the left ear, and the ground electrode was located on the hind leg. ABRs were evoked by tone bursts ranging from 4–64 kHz and were produced by a closed-field electrostatic speaker connected to a driver (EC‑1 and ED‑1; Tucker-Davis Technologies). The 5 ms tone bursts were presented 33.3 times per second; their 0.5 ms onsets and offsets were tapered with a squared cosine function. The speaker’s audio output was transmitted into the ear through a custom acoustic assembly. Sound-pressure levels were measured with a calibrated microphone and preamplifier connected to a conditioning amplifier (4939‑A‑011 and 2690‑A‑0S1, Brüel and Kjær). The response was amplified x10,000 and bandpass filtered at 0.3‑3 kHz (P55, Natus Neurology Inc.). The amplified response was then digitally sampled at 10 μs intervals with a data acquisition device (PCIe-6353, National Instruments) controlled by custom software (LabVIEW 2010, National Instruments). The responses to 1000 bursts were averaged at each intensity level to determine the threshold, defined as the lowest level at which a response peak was distinctly and reproducibly present. For each frequency, sound-pressure level was decreased from 80 dB SPL in 5 dB steps until threshold was reached. DPOAEs were then elicited with an acoustic coupler that allowed for recording ear-canal sound pressure levels via a probe tube concentrically situated within the common sound-delivery tube. The acoustic coupler consisted of two electrostatic speakers (EC‑1, Tucker-Davis Technologies) to generate primary tones and a miniature microphone (EK‑23103, Knowles) to measure ear-canal sound pressure (coupler from Mike Ravicz, EPL, Eye&Ear inf, Boston). The speakers and the microphone were both calibrated using the calibrated microphone described above. The 2f1 − f2 distortion product was measured with f2 = 6–54 kHz, f2/f1 = 1.2, and L1 = L2 = 55‑75 dB SPL. The acoustic signal was amplified by a preamplifier (ER‑10B+, Etymotic Research) and the sound pressure measured in the ear canal was digitally sampled at 10 μs intervals with the data-acquisition system described above. Each frequency pair was presented for 1 s. After fast Fourier transforms had been computed and averaged over ten consecutive traces, the amplitudes of the 2f1 − f2 distortion product and the surrounding noise floor (+/- 100 Hz of 2f1 − f2) for each frequency pair were determined, a procedure requiring 17 s of data acquisition and processing time.
HITS-CLIP
Request a detailed protocolNOVA1 and NOVA2 CLIP was performed on E18.5 wild-type cortex using three biological replicates, as described elsewhere (Licatalosi et al., 2008). High-throughput sequencing was performed at the Rockefeller University Genome Resource Center. Sequence tags were aligned to the mouse genome (mm9) by novoalign. Unique tags were collected by eliminating PCR duplicates.
RNA-seq and analysis
Request a detailed protocolE18.5 mouse cortex RNA from wild-type and Nova2-/- littermate and E18.5 mouse cortex and midbrain/hindbrain RNA from wild-type and Nova1-/- littermate (three biological replicates for each genotype) was prepared using Trizol (ambion) and ribosomal RNA was removed from 1 μg RNA using Ribo-Zero rRNA removal Kit (epicentre). Standard RNA-seq libraries were prepared using TruSeq RNA Sample Preparation Kit v2 (illumina) following manufacturer’s instructions. High-throughput sequencing was performed on Hi-seq 2500 (illumina) to obtain 125 nucleotide paired-end reads at New York Genome Center. Reads were aligned to the mouse genome (mm10) using OLego (Wu et al., 2013). The downstream splicing analysis was performed as described previously (Yan et al., 2015) and then liftovered to mm9. Differential expression analysis was performed with the R package edgeR (Robinson et al., 2009) available on Bioconductor 3.0 (Gentleman et al., 2004).
Accession numbers
Request a detailed protocolRNA-seq and HITS-CLIP data have been deposited in GEO under accession number GSE69711 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE69711).
Data availability
-
Nova HITS-CLIP and RNA-Seq in mouse cortexPublicly available at NCBI Gene Expression Omnibus (accession no: GSE69711).
References
-
Netrin-1: diversity in developmentThe International Journal of Biochemistry & Cell Biology 41:487–493.https://doi.org/10.1016/j.biocel.2008.03.014
-
The development of vestibulocochlear efferents and cochlear afferents in miceInternational Journal of Developmental Neuroscience 15:671–692.https://doi.org/10.1016/S0736-5748(96)00120-7
-
The onconeural antigen Nova-1 is a neuron-specific RNA-binding protein, the activity of which is inhibited by paraneoplastic antibodiesJournal of Neuroscience 16:1114–1122.
-
The neuronal RNA binding protein Nova-1 recognizes specific RNA targets in vitro and in vivoMolecular and Cellular Biology 17:3194–3201.https://doi.org/10.1128/MCB.17.6.3194
-
Usherin is required for maintenance of retinal photoreceptors and normal development of cochlear hair cellsProceedings of the National Academy of Sciences of the United States of America 104:4413–4418.https://doi.org/10.1073/pnas.0610950104
-
Handbook of Mouse Auditory Research (CRC Press)189–204, Cochlear hair cell densities and inner-ear staining techniques, Handbook of Mouse Auditory Research (CRC Press), 10.1201/9781420038736.ch13.
-
Paraneoplastic syndromes involving the nervous systemThe New England Journal of Medicine 349:1543–1554.https://doi.org/10.1056/NEJMra023009
-
RNA regulation in neurologic disease and cancerCancer Research and Treatment 42:125–129.https://doi.org/10.4143/crt.2010.42.3.125
-
Nova autoregulation reveals dual functions in neuronal splicingThe EMBO Journal 24:1608–1620.https://doi.org/10.1038/sj.emboj.7600630
-
In Situ Hybridization Methods339–352, Combining whole-mount in situ hybridization with neuronal tracing and immunohistochemistry, In Situ Hybridization Methods, Springer New York, 99, 10.1007/978-1-4939-2303-8_17.
-
The combined effects of trkB and trkC mutations on the innervation of the inner earInternational Journal of Developmental Neuroscience 16:493–505.https://doi.org/10.1016/S0736-5748(98)00043-4
-
Unc5B associates with LARG to mediate the action of repulsive guidance moleculeThe Journal of Cell Biology 184:737–750.https://doi.org/10.1083/jcb.200807029
-
Targeting of the EphA4 tyrosine kinase receptor affects dorsal/ventral pathfinding of limb motor axonsDevelopment 127:3313–3324.
-
Mechanoelectrical transduction by hair cellsAnnual Review of Biophysics and Biophysical Chemistry 17:99–124.https://doi.org/10.1146/annurev.bb.17.060188.000531
-
The tetranucleotide UCAY directs the specific recognition of RNA by the Nova K-homology 3 domainProceedings of the National Academy of Sciences of the United States of America 97:5740–5745.https://doi.org/10.1073/pnas.090553997
-
Collapse of growth cone structure on contact with specific neurites in cultureJournal of Neuroscience 7:201–212.
-
Transcription factor GATA-3 alters pathway selection of olivocochlear neurons and affects morphogenesis of the earThe Journal of Comparative Neurology 429:615–630.https://doi.org/10.1002/1096-9861(20010122)429:4<615::AID-CNE8>3.0.CO;2-F
-
Hair cell-specific splicing of mRNA for the alpha1D subunit of voltage-gated Ca2+ channels in the chicken's cochleaProceedings of the National Academy of Sciences of the United States of America 94:14889–14893.https://doi.org/10.1073/pnas.94.26.14889
-
Neuronal regulation of alternative pre-mRNA splicingNature Reviews. Neuroscience 8:819–831.https://doi.org/10.1038/nrn2237
-
The cochlear frequency map for the cat: labeling auditory-nerve fibers of known characteristic frequencyThe Journal of the Acoustical Society of America 72:1441–1449.https://doi.org/10.1121/1.388677
-
The ipsilaterally evoked olivocochlear reflex causes rapid adaptation of the 2f1-f2 distortion product otoacoustic emissionThe Journal of the Acoustical Society of America 99:3572–3584.https://doi.org/10.1121/1.414956
-
RNA-binding proteins in human genetic diseaseTrends in Genetics : TIG 24:416–425.https://doi.org/10.1016/j.tig.2008.05.004
-
Neurogenin 1 null mutant ears develop fewer, morphologically normal hair cells in smaller sensory epithelia devoid of innervationJournal of the Association for Research in Otolaryngology 1:129–143.https://doi.org/10.1007/s101620010017
-
Olivocochlear innervation in the mouse: immunocytochemical maps, crossed versus uncrossed contributions, and transmitter colocalizationThe Journal of Comparative Neurology 455:406–416.https://doi.org/10.1002/cne.10490
-
Tuning of synapse number, structure and function in the cochleaNature Neuroscience 12:444–453.https://doi.org/10.1038/nn.2293
-
Agenesis of the corpus callosum: genetic, developmental and functional aspects of connectivityNature Reviews Neuroscience 8:287–299.https://doi.org/10.1038/nrn2107
-
A brain-enriched polypyrimidine tract-binding protein antagonizes the ability of Nova to regulate neuron-specific alternative splicingProceedings of the National Academy of Sciences 97:6350–6355.https://doi.org/10.1073/pnas.110128397
-
Rescuing Z+ agrin splicing in Nova null mice restores synapse formation and unmasks a physiologic defect in motor neuron firingProceedings of the National Academy of Sciences 106:3513–3518.https://doi.org/10.1073/pnas.0813112106
-
BookDevelopment of the inner ear efferent systemIn: Ryugo DK, Fay RR, Popper AN, editors. Auditory and Vestibular Efferents. New York, NY: Springer Verlag. pp. 187–216.https://doi.org/10.1007/978-1-4419-7070-1_7
-
Adaptation of 2f1-2f2 distortion product otoacoustic emission in young-adult and old CBA and C57 miceThe Journal of the Acoustical Society of America 105:3399–3409.https://doi.org/10.1121/1.424668
-
Nova regulates brain-specific splicing to shape the synapseNature Genetics 37:844–852.https://doi.org/10.1038/ng1610
-
RNA binding proteins and the regulation of neuronal synaptic plasticityCurrent Opinion in Neurobiology 16:102–110.https://doi.org/10.1016/j.conb.2006.01.003
-
Boca RatonHandbook of mouse auditory research: from behavior to molecular biology, Boca Raton, CRC Press, 10.1201/9781420038736.
-
OLego: fast and sensitive mapping of spliced mRNA-Seq reads using small seedsNucleic Acids Research 41:5149–5163.https://doi.org/10.1093/nar/gkt216
-
Systematic discovery of regulated and conserved alternative exons in the mammalian brain reveals NMD modulating chromatin regulatorsProceedings of the National Academy of Sciences of the United States of America 112:3445–3450.https://doi.org/10.1073/pnas.1502849112
-
The neuronal RNA-binding protein Nova-2 is implicated as the autoantigen targeted in POMA patients with dementiaProceedings of the National Academy of Sciences of the United States of America 95:13254–13259.https://doi.org/10.1073/pnas.95.22.13254
Article and author information
Author details
Funding
Japan Society for the Promotion of Science (JSPS postdoctoral fellowship for research abroad)
- Yuhki Saito
Naito Foundation (Subsidy for Inter-Institute Researches)
- Yuhki Saito
Hearing Health Foundation
- Soledad Miranda-Rottmann
National Institutes of Health (R01EY024261)
- Harald J Junge
Boettcher Foundation
- Zhe Chen
Canadian Institutes of Health Research (RNI00109)
- Roberto Araya
Howard Hughes Medical Institute
- A J Hudspeth
- Robert B Darnell
National Institutes of Health (R01 NS069473)
- Robert B Darnell
National Institutes of Health (R01 NS081706)
- Robert B Darnell
National Institutes of Health (R01 NS034389)
- Robert B Darnell
Office of Vice President for Research (P30 DC010362)
- Jeremy S Duncan
- Bernd Fritzsch
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Claudia Scheckel, Hun-Way Hwang, Michael J Moore, and Jennifer C Darnell for critical discussion and reading of manuscript and Shengdong Ke for help with bioinformatics processing of HITS-CLIP and RNA-seq data sets, and are grateful to all other members of Robert B Darnell’s laboratory for kind support. We thank M Charles Liberman and the members of his laboratory for technical help and valuable discussion. This work was supported by grants from the National Institutes of Health (NS034389 and NS081706) to RBD. SM’s research was supported in part by a grant of the Hearing Health Foundation and by CIHR grant RNI00109 to RA. ZC was supported by the Boettcher Foundation. HJJ was supported by the NIH (R01EY024261). YS was supported by a JSPS postdoctoral fellowship for research abroad and The Naito Foundation Subsidy for Inter-Institute Researches. BF and JAD were supported by P30 DC010362 and the Office of Vice President for Research at the University of Iowa. RBD and AJH are Investigators of the Howard Hughes Medical Institute.
Ethics
Animal experimentation: This studies were performed in compliance with protocols (#14678 and #07069) approved by the Institutional Animal Care and Use Committee (IACUC) of the Rockefeller University or with protocols (13-185 and 15-002) approved by the Comité de déontologie de l'expérimentation sur les animaux (CDEA) of the Univeristy of Montreal.
Copyright
© 2016, Saito et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 4,455
- views
-
- 954
- downloads
-
- 94
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Neuroscience
Neural implants have the potential to restore lost sensory function by electrically evoking the complex naturalistic activity patterns of neural populations. However, it can be difficult to predict and control evoked neural responses to simultaneous multi-electrode stimulation due to nonlinearity of the responses. We present a solution to this problem and demonstrate its utility in the context of a bidirectional retinal implant for restoring vision. A dynamically optimized stimulation approach encodes incoming visual stimuli into a rapid, greedily chosen, temporally dithered and spatially multiplexed sequence of simple stimulation patterns. Stimuli are selected to optimize the reconstruction of the visual stimulus from the evoked responses. Temporal dithering exploits the slow time scales of downstream neural processing, and spatial multiplexing exploits the independence of responses generated by distant electrodes. The approach was evaluated using an experimental laboratory prototype of a retinal implant: large-scale, high-resolution multi-electrode stimulation and recording of macaque and rat retinal ganglion cells ex vivo. The dynamically optimized stimulation approach substantially enhanced performance compared to existing approaches based on static mapping between visual stimulus intensity and current amplitude. The modular framework enabled parallel extensions to naturalistic viewing conditions, incorporation of perceptual similarity measures, and efficient implementation for an implantable device. A direct closed-loop test of the approach supported its potential use in vision restoration.
-
- Neuroscience
Brain water homeostasis not only provides a physical protection, but also determines the diffusion of chemical molecules key for information processing and metabolic stability. As a major type of glia in brain parenchyma, astrocytes are the dominant cell type expressing aquaporin water channel. How astrocyte aquaporin contributes to brain water homeostasis in basal physiology remains to be understood. We report that astrocyte aquaporin 4 (AQP4) mediates a tonic water efflux in basal conditions. Acute inhibition of astrocyte AQP4 leads to intracellular water accumulation as optically resolved by fluorescence-translated imaging in acute brain slices, and in vivo by fiber photometry in mobile mice. We then show that aquaporin-mediated constant water efflux maintains astrocyte volume and osmotic equilibrium, astrocyte and neuron Ca2+ signaling, and extracellular space remodeling during optogenetically induced cortical spreading depression. Using diffusion-weighted magnetic resonance imaging (DW-MRI), we observed that in vivo inhibition of AQP4 water efflux heterogeneously disturbs brain water homeostasis in a region-dependent manner. Our data suggest that astrocyte aquaporin, though bidirectional in nature, mediates a tonic water outflow to sustain cellular and environmental equilibrium in brain parenchyma.