Abstract
RNA helicases regulate RNA metabolism, but their substrate specificity and in vivo function remain largely unknown. We isolated spontaneous mutant zebrafish that exhibit an abnormal dorsal bend at the beginning of tactile-evoked escape swimming. Similar behavioral defects were observed in zebrafish embryos treated with strychnine, which blocks glycine receptors (GlyRs), suggesting that the abnormal motor response in mutants may be attributable to a deficit in glycinergic synaptic transmission. We identified a missense mutation in the gene encoding RNA helicase Dhx37. In Dhx37 mutants, ribosomal RNA levels were unchanged, whereas GlyR α1, α3, and α4a subunit mRNA levels were decreased due to a splicing defect. We found that Dhx37 can interact with GlyR α1, α3, and α4a transcripts but not with the GlyR α2 subunit mRNA. Overexpression of GlyR α1, α3, or α4a subunits in Dhx37-deficient embryos restored normal behavior. Conversely, antisense-mediated knockdown of multiple GlyR α subunits in wild-type embryos was required to recapitulate the Dhx37 mutant phenotype. These results indicate that Dhx37 is specifically required for the biogenesis of a subset of GlyR α subunit mRNAs, thereby regulating glycinergic synaptic transmission and associated motor behaviors. To our knowledge, this is the first identification of pathologically relevant substrates for an RNA helicase.
Introduction
RNA helicases belong to a large DNA/RNA helicase superfamily and are classified into Upf1-like, Ski2-like, RIG-I-like, NS3/NPH-II, DEAD-box, and DEAH-box subfamilies (Bleichert and Baserga, 2007). These RNA helicases take part in various aspects of RNA metabolism. In yeast, the DEAH-box protein Dhr1p regulates the synthesis of ribosomal RNA (rRNA) and is involved in the preribosomal complex (Colley et al., 2000; Grandi et al., 2002). However, the biological function of the vertebrate homolog Dhx37 has not been investigated. In vivo physiological functions of RNA helicases remain largely unknown, with the exception of a few involved in disease (Moreira et al., 2004; Boon et al., 2007; Pena et al., 2007). Preference for certain substrate RNAs was reported only for RNA helicase A (Hartman et al., 2006). Whether other RNA helicases have specific targets remains unknown.
Zebrafish represent a useful model for motor study. At 2 d postfertilization (dpf), tactile stimulation evokes an escape behavior, which consists of an initial turn and subsequent swimming (Saint-Amant and Drapeau, 1998). Forward genetics has identified zebrafish mutants as being defective in this escape behavior (Granato et al., 1996). For example, bandoneon mutants carry mutations in glrbb, encoding the glycine receptor (GlyR) βb subunit. Since GlyR α and β subunits form heteromeric pentamers, with the β subunit being essential for the synaptic clustering of GlyRs (Lynch, 2004), the GlyR βb mutants exhibited a loss of inhibitory glycinergic synaptic transmission. Due to the consequent loss of the reciprocal inhibition between the left and right sides of the spinal cord, motor neurons activate simultaneously on both sides (Grillner, 2003; Fetcho et al., 2008), which results in bilateral muscle activation and thus dorsal flexure of the body in GlyR βb mutants (Hirata et al., 2005). The same motor deficit is observed in zebrafish embryos treated with a high concentration of strychnine, a GlyR antagonist (Granato et al., 1996; Hirata et al., 2005; McDearmid et al., 2006). On the other hand, the application of strychnine at a lower concentration causes a transient dorsal bend followed by swimming.
In this study, we characterized a zebrafish mutant, dhx37nig1, that exhibits a tactile-evoked dorsal bend, followed by swimming. Positional cloning revealed a missense mutation in the DEAH-box RNA helicase Dhx37. In dhx37nig1 mutants, mRNA levels of selected GlyR α subunit genes were decreased, whereas rRNA levels were unchanged. Indeed, our RNA analyses showed that Dhx37 binds to a subset of GlyR transcripts and regulates their splicing. Overexpression of selected recombinant GlyR α subunits restored the normal touch response in Dhx37-depleted embryos. Thus, Dhx37 plays an essential role in the biogenesis of GlyR α subunit mRNAs and is indispensable for normal escape behavior in vertebrates.
Materials and Methods
Animals.
Zebrafish were bred and raised according to the guidelines set forth by the National Institute of Genetics. The dhx37nig1 mutation was found in zebrafish breeding stock obtained from a pet shop. Animals from either sex were used in this study.
Video recording.
Embryonic behaviors were recorded using a high-speed camera at 200 frames per second (HAS-220, Ditect), as previously described (Hirata et al., 2005).
Electrophysiology.
The dissection protocols for in vivo patch recordings have been described previously (Hirata et al., 2005). To record spontaneous glycinergic currents, 1 μm tetrodotoxin (TTX; Sigma), 10 μm 6-cyano-7-nitroquinoxaline-2,3-dione (CNQX; Sigma), 50 μm aminophosphonovalerate (APV; Sigma), 10 μm bicuculline (Sigma), and 20 μm d-tubocurarine (Sigma) were bath applied.
Mapping, cloning, mRNA rescue, and antisense knockdown.
A mutant carrier fish was crossed with a WIK fish for meiotic mapping. Cloning, mRNA rescue, and antisense knockdown were performed as described previously (Hirata et al., 2005). The following primers and morpholino oligonucleotides (MOs) were used. We performed injection of MOs and/or RNAs in three independent experimental trials and obtained consistent results: zDhx37: ATGGGCAGATTGAGAAAGAAACA, TTAAAAGTCTAGGTTTAGAAGGTCAGA; Dhx37 MO1: TGTTTCTTTCTCAATCTGCCCATGG; Dhx37 MO2: ATCAAGTGTTTTACCTTGTTGCGGA; GlyR α1 MO: ACAAATAAATCCCGAGTGCGAACAT; GlyR α2 MO: GGAGTTTGACCGAAGGGCGAGTCAT; GlyR α3 MO: CAGAGCGCAGTCTTTTCCGAGTCAT; GlyR α4a MO: AAATCCTTATGACCTGAGGGAGCAT; GlyR α4b MO: C GATCCTCCAGATCACAGAAAACAT; GlyR βa MO: G TGCTGACTTACGAGGGACAGATGA; GlyR βb MO: AAATGCACTTACGAAGGACAGTAA; and control MO (random sequence with no target): CCTCTTACCTCAGTTACAATTTATA.
To pharmacologically increase glycinergic transmission for behavioral rescue, NFPS (N[3-(4′-fluorophenyl)-3-(4′-phenylphenoxy)propyl]sarcosine, Santa Cruz Biotechnology), a blocker of Gly transporter 1 (GlyT1), was injected into mutants as described previously (Cui et al., 2005).
RNA analysis.
Northern blotting was performed using a full coding sequence as a probe. Quantitative PCR (qPCR) was performed using a KAPA SYBR Fast qPCR Kit (Kapa Biosystems). The following PCR primers were used: GlyR α1: CTCTCTTCCCAAGGTCTCG, GCCTCGTCCTCCTTCAG; GlyR α2: CTGTACAGCATCAGGCTGAC, TGGTGTATCCGAAGCTCTCC; GlyR α3: GCAGCTGGAGAGTTTTGGTTAC, GCATGTGAACTTGCCTGTGTTG; GlyR α4a: GAATGTGCTTTACAGCATCAGGC, CGTTCATGGTGTAGCCAAAGC; GlyR α4b: GTATAGCATCAGACTCACGCTG, CAGATCATTCATAGTGTAGCCAAAGC; GlyR βa: CGGCCGAATTTCAAAGGAATC, GCAGGAAGATATTCACACGATAATCCATT; GlyR βb: GTTCTCATCAGCATGAGGTTGTC, GTCATCTGTGGTGTAACCAAAGC; Gephyrin a: GAGTGTTTATGAAGCCAGGTTTGC, CTGACACTGGATTACCTGGAAG; Gephyrin b: GGATGCTCTCCAGACCTG, GATGAACTGGAAGCATTCCTGTG; GlyT1: CTACAGAAATGGAGGAGGTGC, CATTCCATAGCCCACACCTTTAAAC; GlyT2: GGTCGGAATCTCTTACATTTATGG, GTACAGACTAAGTGCCAAAATGAATGTG; vesicular inhibitory amino acid transporter (VIAAT): AACGCGATCCAGGGCATG, GAGGATTTTGCCTGTGTAGCAG; and Actb1: CCTTCCTTCCTGGGTATGG, GGGGGAGCAATGATCTTGATC.
Immunostaining.
Immunostaining of spinal cord sections was performed as described previously (Hirata et al., 2005). The following antibodies were used: anti-GlyRα (1:1000; mAb4a, IgG1, Synaptic Systems); anti-gephyrin (1:1000; clone 45, IgG1, Synaptic Systems); anti-acetyl tubulin (1:1000; 6-11B-1, IgG2b, Sigma); Alexa Fluor 488-conjugated anti-mouse IgG1; Alexa Fluor 488-conjugated anti-mouse IgG2b (1:1000; Invitrogen). Quantitative analysis of immunolabeling signals was done using ImageJ (NIH).
Immunoprecipitation and Western blotting.
Immunoprecipitation and Western blots were performed as described previously (Hirata et al., 2005). Anti-FLAG affinity gel and mouse IgG agarose (Sigma) were used for immuneprecipitation and its control. Anti-GlyRα (1:2000), anti-gephyrin (1:2000), anti-acetyl tubulin (1:2000), anti-FLAG (1:1000; M2, Sigma), and peroxidase-conjugated anti-mouse IgG (1:1000; Invitrogen) were used in Can Get Signal immunoreaction enhancer solution (Toyobo). The intensity of bands was quantified using Multi Gauge (Fujifilm).
Results
dhx37nig1 mutants display abnormal escape behavior due to a defect in glycinergic synaptic transmission
A recessive mutation, nig1 (later referred to as dhx37nig1), was identified in our breeding stock of zebrafish. At 2 dpf, tactile stimulation of wild-type zebrafish induced an escape response that consists of an initial turn to the lateral side and subsequent swimming (Movie 1, Clip 1). The initial turn is executed by activation of the trunk muscles on one side, and the subsequent swimming is mediated by alternating muscle contractions on either side of the trunk. By contrast, most nig1 mutants (92%, 44/48 mutants) responded to touch with a dorsal bend of the body followed by swimming (Movie 1, Clip 2). The rest of the mutants (8%, 4/48 mutants) exhibited a strong dorsal bend without subsequent swimming (Movie 1, Clip 3). Since a dorsal flexure of the trunk is typically seen in embryos treated with strychnine, a GlyR antagonist, we assayed the touch response following exposure to a low or high dose of strychnine in the bath and classified the response into the following three groups: level 1 (normal), exhibiting a lateral turn and subsequent swimming (Movie 1, Clip 4); level 2 (mild), exhibiting a dorsal bend followed by swimming of >2 cm (Movie 1, Clip 5); and level 3 (severe), exhibiting a dorsal bend without escape swimming (<2 cm; Movie 1, Clip 6). In the absence of strychnine, all of the wild-type embryos exhibited a normal response (Fig. 1A,G). Application of 30 μm strychnine induced level 2 (21%, 17/81) and level 3 (38%, 31/81) responses in wild-type embryos (Fig. 1B). This strychnine-induced motor deficit was dose dependent, because treatment with 70 μm strychnine worsened the response (level 3: 88%, 57/65; Fig. 1C). Without strychnine, most mutant embryos exhibited the level 2 phenotype (Fig. 1D). Treatment with 30 μm strychnine induced severe motor deficits (level 3: 78%, 40/51; Fig. 1E), as did treatment with 70 μm strychnine (level 3: 87%, 48/55; Fig. 1F). These results suggest that glycinergic transmission is compromised in nig1 mutants. Although nig1 mutants did not show apparent developmental defects, they became thinner and died at 7–10 dpf, possibly from an ineffective motor response and feeding difficulties.
To see whether glycinergic transmission is affected in nig1 mutants, we measured miniature glycinergic currents in motor neurons at 2 dpf. In wild-type zebrafish, spontaneous glycinergic synaptic currents in the presence of TTX, CNQX, APV, bicuculline, and d-tubocurarine, which block Na+ channel, AMPA, NMDA, GABAA, and nACh receptors, respectively, were observed at 0.39 ± 0.17 Hz (n = 7; Fig. 2A), with a frequency comparable to a previous report (0.33 ± 0.07 Hz; Hirata et al., 2005). However, the spontaneous events were less frequently seen in mutants (0.08 ± 0.05 Hz, n = 7, p < 0.01; Fig. 2B), suggesting that glycinergic synaptic transmission is indeed decreased in mutants. GlyT1 is known to remove glycine from the synaptic cleft, and the blockade of GlyT1 enhances glycinergic transmission (Cui et al., 2005; Eulenburg et al., 2005; Mongeon et al., 2008). We therefore injected NFPS, an inhibitor of GlyT1, into the ventricles of mutant embryos. Although 28% of NFPS-injected mutants (5/18 mutants) became immobile, 17% of the mutants (3/18 mutants) responded to touch without dorsal bend at the onset of the escape Movie 1, Clip 7). These results indicate that glycinergic synaptic transmission is compromised in nig1 mutants.
To determine whether nig1 mutants have morphological defects at glycinergic synapses, the distribution of gephyrin, GlyRs, and acetyl tubulin was assessed by immunolabeling of spinal cord sections. Gephyrin and acetyl tubulin were found within the lateral region of the mutant spinal cord as in the wild type (Fig. 2C,D). Synaptic localization of GlyR clusters was also found at the lateral spinal cord in wild type but with lower intensity in mutants (Fig. 2E,F). Quantitative analysis confirmed that GlyR signals in mutant spinal cord was lower (0.56 ± 0.11, n = 6, p < 0.05; Fig. 2G) compared with wild-type zebrafish (n = 6). Western blotting of whole embryo protein extracts with anti-GlyRα also showed that GlyR α subunit levels were decreased in mutants (0.37 ± 0.15, n = 6, p < 0.01; Fig. 2H,I) compared with wild-type zebrafish (n = 6), whereas the amount of gephyrin protein (1.06 ± 0.21, n = 6) and acetyl tubulin protein (1.10 ± 0.14, n = 6) was unchanged. These results, along with those of our behavioral assays and electrophysiology testing, indicate that glycinergic transmission is compromised due to a reduction of GlyR subunit levels in nig1 mutants.
nig1 encodes the DEAH-box RNA helicase Dhx37
We next performed genetic mapping of nig1 locus and identified a missense mutation (L489P) in an RNA helicase, DEAH-box protein 37 (Dhx37; GenBank accession number AB739007). To address whether the L489P mutation is responsible for dhx37nig1 mutant phenotype, we injected control (green fluorescent protein) RNA, wild-type Dhx37 RNA, or Dhx37 L489P RNA into recently fertilized embryos from dhx37nig1 heterozygous carriers and assayed touch responses. A quarter of control RNA-injected progeny (29%, 22/75) and Dhx37 (L489P) RNA-injected progeny (27%, 20/75) displayed the mutant response, whereas most of the progeny injected with wild-type Dhx37 RNA (94%, 61/65) exhibited normal escape behavior, and only 6% (4/65) displayed an abnormal response. The molecular identification of dhx37 as the causative gene was also supported by knockdown of Dhx37 using antisense MOs. Antisense Dhx37 MO1, which blocks translation of Dhx37, was injected into wild-type embryos. Most MO1-injected embryos (82%, 71/87) displayed an abnormal touch response, whereas all control MO-injected embryos exhibited a normal response. Similarly, injection of Dhx37 MO2, which interferes with splicing of dhx37 mRNAs also induced a mutant phenotype (91%, 84/92). These RNA rescue and antisense knockdown experiments confirm that the dhx37nig1 mutation is responsible for the abnormal dorsal bend in the mutant escape behavior.
Dhx37 regulates biogenesis of GlyR α subunit mRNAs
To investigate why a defect in the RNA helicase Dhx37 affects glycinergic synapses, we measured the expression of key genes involved in glycinergic transmission. These included five GlyR α subunits (α1, α2, α3, α4a, and α4b), two β subunits (βa and βb), two isoforms of gephyrin (a and b), two glycine transporters (GlyT1 and GlyT2), and VIAAT. Our qPCR revealed that GlyR α1, α3, α4a, and βa subunits and gephyrin b transcripts were decreased in mutants (α1: 0.74 ± 0.09, n = 6, p < 0.05; α3: 0.74 ± 0.08, n = 6, p < 0.05; α4a: 0.76 ± 0.07, n = 6, p < 0.05; βa: 0.62 ± 0.05, n = 6, p < 0.01; gephyrin b: 0.62 ± 0.08, n = 5, p < 0.01; Fig. 3A).
Some RNA helicases are involved in pre-mRNA splicing. Since qPCR can detect transcripts that are not fully spliced (i.e., could have a retained intron at another junction distant to the qPCR site), this method does not accurately measure the levels of fully spliced mRNAs. To quantify the amount of fully spliced mRNA, we performed Northern blotting and confirmed that GlyR α1, α3, and α4a subunits and gephyrin b mRNAs were decreased in dhx37nig1 mutants (α1: 0.34 ± 0.17, n = 3, p < 0.05; α3: 0.26 ± 0.12, n = 3, p < 0.05; α4a: 0.24 ± 0.11, n = 3, p < 0.01; gephyrin b: 0.37 ± 0.15, n = 3, p < 0.05; Fig. 3B). The level of actin β1 mRNA was unchanged in mutants (1.03 ± 0.14, n = 3). We also tried to detect GlyR βa mRNA, but specific bands were not detectable even in wild-type fish, presumably because the expression level is very low at this stage (Hirata et al., 2005). The amounts of 28S and 18S rRNA were comparable between wild-type and mutant fish (28S rRNA: 0.98 ± 0.07, n = 3; 18S rRNA: 1.04 ± 0.13, n = 3). Our qPCR and Northern blot analyses indicate that mRNA levels for GlyR and gephyrin are decreased in dhx37nig1 mutants, whereas rRNA biogenesis remains unchanged.
We next assessed the splicing of GlyR and gephyrin genes by RT-PCR, revealing that unspliced transcripts for the GlyR α4a subunit are increased in dhx37nig1 mutants (Fig. 3C). Similarly, unspliced gephyrin b products were upregulated. By contrast, splicing defects in actin β1 transcripts were not observed in dhx37nig1 mutants. We also attempted to detect unspliced transcripts for GlyR α1 and α3 subunit genes, but they were not accessible, because all of the introns are too long to amplify alongside normal short products. These results suggest that Dhx37 is involved in pre-mRNA splicing and that Dhx37 has substrate specificity for certain transcripts important for glycinergic synapse function.
The substrate specificity of Dhx37 was further investigated using an RNA immunoprecipitation assay. HEK293 cells were cotransfected with GlyR cDNAs that contained 5′- and 3′-UTR sequences together with a FLAG-tagged Dhx37 expression vector. Cell extracts were immuneprecipitated with anti-FLAG or with control IgG and subjected to RT-PCR. FLAG-Dhx37 immunoprecipitates contained GlyR α1, α3, and α4a subunit transcripts, but not α2 subunit transcripts (Fig. 3D). Taken together, this assay clearly indicates that Dhx37 has substrate specificity for selected GlyR α subunit transcripts.
Normal touch response is restored by overexpression of GlyR α subunits in dhx37nig1 mutants
Since RNA analysis revealed that the expression of GlyR α1, α3, α4a, and βa subunits and gephyrin b was decreased in dhx37nig1 mutants, we investigated whether these GlyR subunits/gephyrin isoforms were also responsible for the behavioral deficiency. Capped RNAs for various GlyR subunits and gephyrin were synthesized in vitro and injected into wild-type embryos along with Dhx37 MO2 in three independent experimental trials. Injection of MO2 alone induced aberrant escape behavior, but coinjection of GlyR α1, α3, or α4a subunit RNAs with MO2 increased the rate of normal behavior (α1: 51%, 39/77; α3: 56%, 34/61; α4a: 51%, 35/69; Fig. 4A). Mixed application of GlyR α1, α3, and α4a subunit RNAs with MO2 also restored the normal escape (56%, 18/32). On the other hand, neither GlyR α2 nor α4b subunit RNAs restored the normal touch response (α2: 13%, 10/75; α4b: 8%, 5/66 showed normal responses). Likewise, embryos injected with MO2 and either GlyR βa or βb subunit RNA, or gephyrin a or b RNA failed to exhibit recovery of the touch response (data not shown). Thus, the shortage of GlyR α subunit, but not of GlyR β or gephyrin, is responsible for the abnormal escape response in dhx37nig1 mutants.
To assess which α subunit (α1, α3, or α4a) is primarily responsible for the GlyR α subunit depletion in mutants, we knocked down each GlyR subunit in wild-type embryos. Upon tactile stimulation, all control MO-injected embryos exhibited a lateral turn and swimming (Fig. 4B). However, injection of the GlyR α1 MO caused the mild level 2 phenotype in some fish (9%, 5/58). Injection of the GlyR α3 MO increased the ratio of level 2 phenotype (15%, 13/85), while the GlyR α4a MO had a stronger effect in a greater percentage of fish (71%, 69/97). Knockdown of the multiple GlyR α subunits had an additive effect—mixed application of two MOs decreased the touch response compared with single-MO injections (α1 + α3: 30.2%, 16/53; α1 + α4a: 82%, 81/99; α3 + α4a: 88%, 83/94 showed level 2 and 3 responses). Moreover, the escape behavior of embryos injected with all three MOs (α1, α3, and α4a) was severely impaired (level 3 response, 76%, 48/63). Although βa subunit morphants did not exhibit any behavioral defects, knockdown of GlyR βb led to the severe level 3 response (85%, 41/48), similar to the triple knockdown of α1, α3, and α4a subunits.
Discussion
In summary, we have identified a zebrafish dhx37nig1 mutant that displays an abnormal dorsal bend at the beginning of the escape response. GlyR α1, α3, and α4a subunit mRNAs were downregulated and mis-spliced in dhx37nig1 mutants. The depletion of GlyR α1, α3, and α4a subunits is responsible for the motor deficits observed, since these could be improved by restoring expression of these GlyR α subunits. Mutations in human DHX37 may cause startle disease/hyperekplexia associated with defective glycinergic transmission.
Substrate specificity of Dhx37
The physiological significance of RNA helicases has been characterized in only a few cases. For example, mutations in the RNA helicase senataxin are responsible for neurodegenerative diseases (Moreira et al., 2004). Retinitis pigmentosa is caused by a dysfunction of pre-mRNA splicing factor 8 (Boon et al., 2007; Pena et al., 2007). However, the pathologically relevant substrate of these RNA helicases remains unknown. We demonstrated that Dhx37 interacts with selected GlyR α subunit mRNAs and specifically regulates biogenesis of these transcripts. Does Dhx37 have the other targets? Most embryos injected with MO2 and Dhx37 RNA displayed a normal touch response (90%, 28/31), while embryos injected with MO2 and mixture of GlyR α1, α3, and α4a subunit RNAs showed motor recovery at a lower efficiency (56%, 18/32). This suggests that Dhx37 has targets other than GlyR subunit mRNAs.
Subunit composition of GlyRs
The initial characterization of purified GlyRs suggested that GlyRs were composed of 3α and 2β subunits (Langosch et al., 1988), but recent functional and structural studies revised the stoichiometry to 2α3β (Grudzinska et al., 2005; Dutertre et al., 2012; Yang et al., 2012). Although our knockdown experiments indicate that α4a appears to be the major GlyR α subunit involved in the escape behavior, the α1 and α3 subunits also contribute. We also showed that the βb but not βa subunit is essential for touch responses. This is consistent with our previous finding that the βb subunit is expressed by 1 dpf, while the expression of βa is very low until 3 dpf (Hirata et al., 2005). Taken together, it is likely that GlyR α1, α3, and α4a subunits coassemble with the GlyR βb subunit, either in different receptor complexes or in heteromeric combinations, and that knockdown of all three α subunits is required to mimic the severe GlyR βb-deficient phenotype.
Footnotes
This work was supported by Grants-in-Aid for Scientific Research (23115719 and 25920008) from the Ministry of Education, Culture, Sports, Science, and Technology of Japan, and by Takeda Science Foundation to H.H.; and by Action Medical Research (1966) to R.J.H. We thank Megumi Takahashi, and Drs. Toshinobu Fujiwara, Norimasa Iwanami, Masashi Tanimoto, Shin Takagi, Yoichi Oda, Koichi Kawakami, Sean E. Low, and John Y. Kuwada for preliminary findings, materials, and discussion.
The authors declare no competing financial interests.
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.
- Correspondence should be addressed to Hiromi Hirata. Center for Frontier Research, National Institute of Genetics, 1111 Yata, Mishima 411-8540, Japan. hihirata{at}nig.ac.jp
This article is freely available online through the J Neurosci Author Open Choice option.