Abstract
Loss-of-function mutations in Reelin and DAB1 signaling pathways disrupt proper neuronal positioning in the cerebral neocortex and hippocampus, but the underlying molecular mechanisms remain elusive. Here, we report that heterozygous yotari mice harboring a single autosomal recessive yotari mutation of Dab1 exhibited a thinner neocortical layer 1 than wild-type mice on postnatal day (P)7. However, a birth-dating study suggested that this reduction was not caused by failure of neuronal migration. In utero electroporation-mediated sparse labeling revealed that the superficial layer neurons of heterozygous yotari mice tended to elongate their apical dendrites within layer 2 than within layer 1. In addition, the CA1 pyramidal cell layer in the caudo-dorsal hippocampus was abnormally split in heterozygous yotari mice, and a birth-dating study revealed that this splitting was caused mainly by migration failure of late-born pyramidal neurons. Adeno-associated virus (AAV)-mediated sparse labeling further showed that many pyramidal cells within the split cell had misoriented apical dendrites. These results suggest that regulation of neuronal migration and positioning by Reelin-DAB1 signaling pathways has unique dependencies on Dab1 gene dosage in different brain regions.
Significance Statement
DAB1 is a cytoplasmic adaptor protein essential for transmission of the extracellular Reelin signal to cytoplasmic proteins that regulate cortical development. In this study, we found that Dab1 is haplosufficient for the regulation of neuronal migration but haploinsufficient for control of layer 1 thickness in the cerebral neocortex. Alternatively, the migration of a subpopulation of hippocampal pyramidal neurons is sensitive to Dab1 gene haploinsufficiency. This study suggests that neural development in the cerebral neocortex and hippocampus are differentially sensitive to Dab1 gene dose.
Introduction
In the mammalian cerebral neocortex, excitatory neurons are born mainly in the ventricular and subventricular zones and migrate radially toward the pial surface by passing through preexisting neuronal layers. Therefore, excitatory cortical neurons are aligned in an “early-born deep” and “late-born superficial” sequence (Cooper, 2008; Honda et al., 2011; Ohtaka-Maruyama and Okado, 2015). Radially migrating neurons express receptors for Reelin, a large secreted glycoprotein mainly produced by Cajal-Retzius cells in the marginal zone that is critical for the positioning of excitatory neurons (D’Arcangelo et al., 1995; Hirotsune et al., 1995; Ogawa et al., 1995). Activation of the Reelin receptors apolipoprotein E receptor 2 (APOER2) and/or very low-density lipoprotein receptor (VLDLR) stimulates Src-mediated tyrosine phosphorylation of the nucleocytoplasmic-shuttling adaptor protein DAB1 (Howell et al., 1997, 2000; Trommsdorff et al., 1999; Honda and Nakajima, 2006, 2016), which in turn interacts with a host of molecular pathways controlling neuronal migration (Assadi et al., 2003; Bock et al., 2003; Pramatarova et al., 2003; Ballif et al., 2004; Chen et al., 2004; Huang et al., 2004; Feng et al., 2007; Franco et al., 2011; Jossin and Cooper, 2011; Sekine et al., 2012; Matsunaga et al., 2017). Thus, the null mutation of Reelin (Falconer, 1951; Montgomery et al., 1994; Bar et al., 1995; D’Arcangelo et al., 1995), the double null mutation of Apoer2 and Vldlr (D’Arcangelo et al., 1999; Hiesberger et al., 1999; Trommsdorff et al., 1999), and Dab1 null mutation (Howell et al., 1997; Sheldon et al., 1997; Yoneshima et al., 1997; Kojima et al., 2000) result in developmental defects in the cerebral neocortex, hippocampus, and cerebellum. In cerebral neocortex, for instance, Reelin, or Dab1 null mutation results in the absence of layer 1 and severe neuronal migration failure (Hamburgh, 1963; Caviness, 1982; D’Arcangelo, 2005). Similarly, excessive cytosolic DAB1 protein causes severe migration failure in wild-type mice (Honda and Nakajima, 2016), suggesting that the regulation of cytosolic DAB1 protein amount is important for the migration of excitatory neurons. Naturally occurring deletions, duplications, and loss-of-function mutations of genes sometimes cause gene dosage effects on function (Veitia and Potier, 2015). For diploid organisms, the protein product from one allele is often sufficient to allow normal growth, development, and physiological function; however, a small number of genes exhibit dosage sensitivity, leading to deficits in the heterozygous state (termed haploinsufficiency). In general, genes encoding enzymes tend to be dose-insensitive, whereas genes that encode proteins with structural, transcriptional, or regulatory functions are more likely to be dose-sensitive (Kondrashov and Koonin, 2004; Morrill and Amon, 2019). DAB1 is a versatile protein that mediates multiple functions, but the underlying signaling pathways for many of these functions have not been fully elucidated. Therefore, it is important to dissect DAB1 pathways and determine the gene-dose sensitivity of each. Further, among the previous studies related to the Reelin-DAB1 signaling pathway, there are several reports showing that attenuation of Reelin signaling pathway caused defects in the cerebral cortical development (Herrick and Cooper, 2002; Hack et al., 2007; Kohno et al., 2015; de Frutos et al., 2016; Hirota et al., 2018; Hirota and Nakajima, 2020) and might be associated with the pathogenesis of a number of neuropsychiatric disorders including schizophrenia (Impagnatiello et al., 1998; Guidotti et al., 2000), bipolar disorder (Guidotti et al., 2000), autism spectrum disorder (Fatemi et al., 2001, 2005), and Alzheimer’s disease (Chin et al., 2007). However, it is not clear whether there is Dab1 gene-dosage dependency for each reported phenotype.
To examine whether Dab1 exhibits haploinsufficiency or haplosufficiency in various neurodevelopmental functions, we compared neocortical and hippocampal structures between wild-type and heterozygous yotari mice harboring a single Dab1 null mutation (Sheldon et al., 1997; Yoneshima et al., 1997; Kojima et al., 2000). We found that heterozygous yotari mice at postnatal day (P)7 exhibited reduced neocortical layer 1 thickness. In addition, some neurons located in the superficial layer tended to elongate their dendrites ectopically within layer 2. Also, we found that the late-born CA1 pyramidal cell layer in the caudal hippocampus is abnormally split in heterozygous yotari mice. These findings suggest that the Dab1 gene dosage is important for maintenance of the layer 1 thickness and dendrite guidance of superficial-layer neurons in the cerebral neocortex and for the migration of late-born pyramidal neurons in the caudal hippocampus.
Materials and Methods
Mice
Heterozygous yotari or wild-type females were mated with heterozygous yotari males (RRID:IMSR_RBRC05456). Noon on the day the vaginal plug was found was designated as embryonic day (E)0.5. All animal procedures were performed in accordance with the Keio University Animal Care and Use Committee’s regulations, Japanese Government Law Concerning the Protection and Control of Animals, and Japanese Government Notification of Feeding and Safekeeping of Animals.
Plasmids
A fluorescent protein expression vector, pCAGGS-RG, capable of Cre-dependent switching from DsRed to enhanced green fluorescent protein (EGFP) expression, and a Cre expression vector, pDCX-Cre, expressing Cre under control of the Doublecortin (DCX) promoter were used for sparse labeling of neurons as has been shown before (Kohno et al., 2015). To construct an adeno-associated virus (AAV)-DCX-Cre vector, the DCX promoter was obtained from Dcx4kb-EGFP (X. Wang et al., 2007; kindly provided by Q. Lu, Beckman Research Institute of the City of Hope), Cre cDNA was obtained from pCXN-Cre (Koresawa et al., 2000; kindly provided by S. Miyagawa, Osaka University), and the AAV plasmid backbone were obtained from AAV-CAG-GFP (a kind gift from K. Svoboda, Addgene plasmid #28014; http://n2t.net/addgene:28014; RRID:Addgene_28014). To construct AAV-CAG-Double-floxed Inverse Orientation (DIO)-EGFP and AAV-CAG-DIO-mCherry vectors, the DIO cassette sequence information was obtained from the double-floxed inverse ChR2-EYFP vector (Sohal et al., 2009), and the sequence was synthesized by Genscript. This DIO cassette was then inserted into the AAV-CAG backbone obtained from AAV-CAG-GFP, yielding AAV-CAG-DIO1. The cDNAs encoding EGFP and mCherry were amplified by PCR from pEGFP-N1 (Clontech) and pmCherry (Clontech), respectively, and inserted into AAV-CAG-DIO1 in reverse orientation with respect to the CAG promoter to allow Cre-dependent expression of the reporter genes. The adenovirus helper plasmid pAdΔF6 and vector pAAV2/9 were kindly provided by the University of Pennsylvania Penn Vector Core.
Bromodeoxyuridine (BrdU) labeling
Bromodeoxyuridine (Merck) was first dissolved in sterile phosphate-buffered saline (PBS) at 10 mg/ml, and the solution was administered to pregnant dams by intraperitoneal injection at 50 μg BrdU/g body weight on gestation days 12.5, 14.5, and 16.5.
Immunohistochemistry
Mouse pups treated as indicated were anesthetized using isoflurane and brains fixed by transcardial perfusion with 4% paraformaldehyde (PFA). Brains were then isolated, postfixed in PFA at 4°C for 2 h, washed with PBS, incubated sequentially in 10%, 20%, and 30% sucrose/PBS for cryopreservation, embedded in a mixture of optimal cutting temperature (OCT) compound (Sakura Finetek) and 30% sucrose/PBS solution, and cryosectioned in the coronal plane at 20 μm thick. For layer marker staining, sections were washed three times (5 min/wash) in PBS containing 0.1% Tween 20 (PBStw) at room temperature (RT), autoclaved at 105°C for 5 min in 0.01 m citrate buffer (pH 6.0) for antigen retrieval, then incubated sequentially with 5% bovine serum albumin (BSA)/PBStw at RT for 5 min and primary antibodies against BRN2 (Santa Cruz Biotechnology, sc-6029; RRID:AB_2167385), RORB (Perseus Proteomics, PP-N7927-00; RRID:AB_1964364), or TBR1 (abcam, ab31940; RRID:AB_2200219) in 5% BSA/PBStw at RT for 1 h. For visualization of immunolabeling, sections were washed three times in PBStw and incubated first with Alexa 488-labeled secondary antibody against mouse IgG (Thermo Fisher Scientific, A-21202; RRID:AB_141607) or rabbit mouse IgG (Thermo Fisher Scientific, A-11008; RRID:AB_143165) with 4′,6-diamidino-2-phenylindole (DAPI; Life Technologies) at RT for 1 h. After three washes in PBS, sections were mounted with PermaFluor Aqueous Mounting Medium (Thermo Fisher Scientific). For double staining with BrdU and either BRN2 or CTIP2 (official symbol is BCL11B) antibodies, sections were autoclaved, blocked with 5% BSA/PBStw, incubated with goat anti-BRN2 (Santa Cruz, sc-6029; RRID:AB_2167385) or rat anti-CTIP2 (abcam, ab18465; RRID:AB_2064130) in 5% BSA/PBStw, washed, and incubated with Alexa 555-labeled anti-goat IgG (Thermo Fisher Scientific, A-21432; RRID:AB_2535853) or Alexa 594-labeled anti-rat IgG (Thermo Fisher Scientific, A-21209; RRID:AB_2535795). Sections were rinsed with PBStw, fixed in 4% PFA for 10 min, washed three times in PBStw (5 min/wash) to remove trace PFA, incubated with 2N HCl at 37°C for 30 min, rewashed three times in PBStw (5 min/wash) to remove HCl, blocked in 5% BSA/PBStw for 5 min, incubated with mouse anti-BrdU (BD Biosciences, 347580; RRID:AB_10015219) for 1 h, washed again three times in PBStw, and incubated with an Alexa 488-labeled secondary antibody against mouse IgG (Thermo Fisher Scientific, A-21202; RRID:AB_141607) containing DAPI for 1 h. After three washes in PBStw, sections were embedded in PermaFluor. For DAPI staining alone, cryosections were washed in PBStw for 5 min, incubated with PBStw containing DAPI for 5 min, washed with PBStw, and embedded in PermaFluor. For immunohistochemical staining of mCherry in the hippocampus, OCT-embedded brains were sectioned in the coronal plane at 100 μm thick and floated on PBS in a 12-well plates. Selected sections were transferred to 24-well plates containing PBS, and treated with PBS containing 5% BSA and 0.5% Triton X-100 (blocking/permeabilization buffer) for 30 min to 1 h at RT on a shaker. After removal of the blocking/permeabilization buffer, sections were treated with primary antibody against RFP (Rockland, 600-401-379; RRID:AB_2209751) at 4°C overnight on a shaker, washed three times with blocking/permeabilization buffer, and incubated with Alexa Fluor 555-labeled secondary antibody against rabbit IgG (Thermo Fisher Scientific, A-31 572; RRID:AB_162543) containing DAPI for 2–3 h on a shaker. After three additional washes in PBS, the sections were embedded in PermaFluor. Confocal images were acquired using an Olympus FV1000 confocal microscope (Olympus), Leica SP8 confocal microscope (Leica microsystems), or Zeiss LSM700 confocal microscope (Zeiss).
Nissl staining
For Nissl staining, OCT-embedded brains obtained from P0 or P7 pups were sectioned in the coronal plane at 20 μm thick, and the sections floated on a small droplet of distilled water (DW) placed on a glass slide. The excess DW under the section was removed with a pipette and allowed to air dry. Then, sections were washed three times with PBS (5 min/wash), dehydrated in an ascending ethanol gradient (70%, 80%, 90%, 95%, 99%, and 100%, each for 1 min), and then rehydrated in a descending ethanol gradient (100%, 99%, 95%, 90%, 80%, and 70%, each for 1 min). Sections were treated with 0.1% cresyl violet solution for 30 min, washed with DW for 5 min, destained with 0.005% acetic acid containing absolute ethanol for 30 s, and washed with DW for 5 min. Stained sections were dehydrated in the same ascending ethanol gradient, cleared with 100% xylene, and mounted with Entellan mounting medium (Merck).
In utero electroporation
In utero electroporation was conducted as previously described (Tabata and Nakajima, 2001, 2008). We first prepared an anesthetic mixture of three drugs in PBS, 75 ng/ml medetomidine hydrochloride (Domitor; Nippon Zenyaku Kogyo), 400 ng/ml midazolam (Sandoz), and 0.5 mg/ml butorphanol (Vetorphale; Meiji Seika Pharma). Pregnant mice at E16.5 were deeply anesthetized by the three-drug mixture at 10 μl/body weight (g), and the uterus was exteriorized. Plasmid DNA solution containing 2.5 μg/μl pCAGGS-RG and 200 ng/μl pDCX-Cre supplemented with 0.1% fast green at 1:10 v/v was injected into the lateral ventricle of embryos using a glass capillary. Then, while holding the head of individual embryos with a tweezer-type electrode (CUY21; NEPA Gene) from outside the uterine wall, 35-V electric pulses were applied four times, each for 50 ms at 950-ms intervals, using an electroporator (NEPA21; NEPA Gene). The uterus was then pushed back into the abdominal cavity, and 75 ng/ml atipamezole hydrochloride (Antisedan; Nippon Zenyaku kogyo) solution dissolved in PBS was injected at the same volume as the three-drug anesthetic mixture. Electroporated mice allowed to continue gestation and fixed by transcardial perfusion with 4% PFA at P7. Cryosections were prepared as described above.
Recombinant adeno-associated virus (rAAV) preparation
Recombinant adeno-associated viruses were prepared as described previously (Burger and Nash, 2016). A mixture of plasmid DNAs including (1) AAV-DCX-Cre, pAdΔF6, and pAAV2/9, (2) AAV-CAG-DIO-EGFP, pAdΔF6, and pAAV2/9, or (3) AAV-CAG-DIO-mCherry, pAdΔF6, and pAAV2/9, was used to transfect AAVpro 293T cells (Takara Bio) seeded on twenty 15-cm dishes using polyethyleneimine “Max” (PEI-Max; Polysciences). After 2 d, cells were collected and lysed by three freeze–thaw treatments and lysates treated with Benzonase nuclease (Millipore). The rAAVs were purified from the lysates by ultracentrifugation using a discontinuous gradient of iodixanol (OptiPrep; Abbott Diagnostics Technologies AS). Banded rAAVs were collected from the 40% iodixanol layer, and the iodixanol buffer exchange to PBS and concentration of the rAAV was performed using a Vivaspin 20 ultrafiltration unit (100 kDa molecular weight cut off; Sartorius). For determination of the genomic titer, rAAVs were treated sequentially with DNaseI and Proteinase K, and the quantity of viral DNAs was determined by quantitative PCR using a StepOnePlus real-time PCR system (Applied Biosystems).
rAAV injection
New-born pups (P0) were anesthetized on ice, and injected with AAV-DCX-Cre at ∼6 × 108 vector genomes (vg)/μl, AAV-CAG-DIO-EGFP at ∼2 × 108 vg/μl, and AAV-CAG-DIO-mCherry at ∼2 × 108 vg/μl supplemented with 1:10 (v/v) 0.1% fast green into the lateral ventricles at 3 μl per ventricle using glass capillary connected to a microinjector. The rAAV-injected pups were placed on a warming plate until they regained normal color and full typical new-born activity, and then returned to their mothers. On P7, rAAV-injected mice were fixed by transcardial perfusion of 4% PFA and 100-μm-thick brain cryosections were prepared in the coronal plane as described above.
Measurement of layer 1, cortical wall, and relative layer 1 thickness
Layer thicknesses were measured on digital images of brain coronal sections taken from the dorsolateral region of the rostral (somatosensory and auditory), middle (visual, parietal, and auditory), and caudal (visual and auditory) cortex using Fiji software (RRID:SCR_002285; Schindelin et al., 2012). To measure neocortical layer 1 thickness, a straight line was drawn along the shortest radial distance from the pial surface to the ventricular surface (Fig. 2A,H; Extended Data Fig. 2-1). Intersecting points with the pial surface, layer 1–layer 2 border, and ventricular surface were designated A, B, and C, respectively. First, layer 1 thickness between points A and B, and cortical wall thickness between points A and C were measured. Then, a point A + 1 located 100 μm to the right of point A was drawn on the pial surface. Similarly, points A + 2 to A + 5 were drawn on the pial surface to the right of A + 1, while points A − 1 to A − 5 were drawn on the pial surface to the left of point A. Straight lines were drawn toward the ventricular surface from points A + 1 to A + 5 and A − 1 to A − 5. These lines were measured as layer 1 thickness and cortical wall thickness, respectively, using Fiji software. Measurements were performed using six P7 brains from six different litters. Relative layer 1 thickness was calculated by dividing layer 1 thickness by cortical wall thickness and multiplying by 100. Measurement results were compared between groups by independent samples Student’s t test using the R language environment (R Foundation for Statistical Computing; https://www.r-project.org/; RRID:SCR_001905) as described in Table 1.
Extended Data Figure 1-1
Nissl staining of the cerebral cortex at P0. Nissl-stained coronal sections of cerebral cortex from wild-type (+/+) and heterozygous yotari (+/yot) mice obtained at P0. Scale bar: 1 mm. Download Figure 1-1, EPS file.
Extended Data Figure 2-1
Schematic illustrations showing the three brain regions along the rostro-caudal axis and the 11 positions analyzed in each brain region. Schematic illustrations of the three brain regions analyzed for layer 1 and cortical wall thickness, the rostral region containing the somatosensory and auditory cortices (Rostral), middle brain region containing the visual, parietal, and auditory cortices (Middle), and the caudal brain region containing the visual and auditory cortices (Caudal). In the lower panel, the 11 positions analyzed in each brain region (A, A + 1–5, A − 1–5) are shown by dashed lines. S: somatosensory cortex; A: auditory cortex; V: visual cortex, P: parietal cortex. Download Figure 2-1, EPS file.
Measurements of relative cell position and distance from the ventricular surface
Distances from the ventricular surface to marker-positive cells (D1; see Fig. 3B), distances from the ventricular surface to the pial surface through each of these cells (D2; see Fig. 3B), and distances from the ventricular surface to BrdU/BRN2 double-positive cells (Fig. 4) were measured using Fiji software. Relative position was then calculated as D1/D2 × 10. These measurements were performed on three brains obtained from three different litters, and results were compared between groups by independent samples Student’s t test, unpaired Welch’s t test, or Mann–Whitney U test using the R language environment as described in Table 1.
Extended Data Figure 3-1
The neocortical region from layer 2 to the ventricular surface of heterozygous yotari mice is expanded and the number of cells is increased but there is no change in the cell density. A, Schematic showing the dorsolateral cerebral neocortex. The measurement area was from the boundary of layer 1 (L1) and layer 2 (L2) to the ventricular surface (VS; 200-μm width, dark gray area). PS: pial surface; L6: layer 6. B–D, Sections of the rostral brain region containing the somatosensory and auditory cortex (Rostral), middle brain region containing the visual, parietal, and auditory cortex (Middle), and caudal brain region containing the visual and auditory cortex (Caudal) from six wild-type (W) and heterozygous yotari mice (He) each, all from different litters. Confocal images were then taken of the somatosensory, parietal, and visual cortex, respectively. Numbers of DAPI-positive cells within the measurement region are shown in C and cell density in D on combination dot and scatter plots. The black dots with a black bar indicate the mean and SD, respectively, and each gray dot is a raw datapoint. The total numbers of analyzed mice are shown below each graph. Statistical significance was evaluated by independent samples Student’s t test (B–D: Middle) or Mann–Whitney U test (D: Rostral and Caudal). The p-values are labeled on the graph. *p < 0.05, ** p < 0.01, ns (not significant) p > 0.05. Download Figure 3-1, EPS file.
Extended Data Figure 3-2
Heterozygous yotari mice exhibit an expanded RORB-positive cell layer compared to wild-type mice. A, Schematic illustration showing how layer thickness was defined as the distance from the 20th to 80th percentile of layer marker-positive cells. B, Thicknesses of TBR1-positive, RORB-positive, or BRN2-positive cell layers measured on coronal sections of the somatosensory (Rostral), parietal (Middle), and visual (Caudal) cortex from three wild-type (W) and three heterozygous yotari mice (He), each of different litters. The black dot and bar in the leftmost dot plot indicate the mean and SD, respectively, each gray dot indicates the mean from one brain, and each black dot in the scatter plot is a datapoint obtained from the rostral, middle, or caudal brain regions. The total number of analyzed coronal brain sections is shown at the bottom of the graph. Groups were compared by independent samples Student’s t test. The p-values are labeled on the graph. **p < 0.01, ns (not significant) p > 0.05. Download Figure 3-2, EPS file.
Measurement of neocortical layer thickness by the distribution of layer marker-positive cells
For quantitative comparisons, we defined layer thickness as the distance between the 20th and 80th percentile of layer marker-positive cells (Extended Data Fig. 3-2A). The markers used were TBR1, RORB, or BRN2. Measurement was performed using coronal sections of somatosensory, parietal, and visual cortices obtained from three brains from different litters. Then the average thickness of each brain along the rostro-caudal axis was calculated. Results were compared between groups by independent samples Student’s t test using the R language environment.
Measurement of pyramidal neuron apical dendrite angle in cerebral neocortex and hippocampus
To measure the apical dendrite angle of pyramidal neurons in the cerebral neocortex, a perpendicular line was drawn from the pial surface to the center of the target cell body, from which another straight line was drawn to the first branch point of the apical dendrite (Fig. 7D,E). The angle between these two straight lines was measured using Fiji software and defined as the apical dendrite angle. To measure the apical dendrite angle of hippocampal pyramidal neurons, a line perpendicular to the pyramidal cell layer was drawn to the center of the target cell body and another from the target cell body center to the first branch point of the apical dendrite (Fig. 6H). The angle between the lines was then measured using Fiji software. Measurements were conducted on more than three brains, all obtained from different litters. Angles were compared between groups by the Steel–Dwass test using the R language environment as described in Table 1.
Measurement of apical dendrite length and number of branching points within neocortical layer 1 or layers 2
Sparsely labeled EGFP-positive cells using electroporation as described above were selected, and their apical dendrite length and number of apical dendrite branching points within layer 1 and layer 2 were measured using Fiji software, respectively. Brains were obtained from more than three mice, each from different litters, and results compared between groups by independent samples Student’s t test or Mann–Whitney U test using the R language environment as described in Table 1.
Measurement of the total apical dendrite length and total number of branching points in the hippocampus
Total apical dendrite length and branch point number were determined using Fiji software. Measurements were performed using more than three brains derived from different litters, and results compared between groups by Tukey–Kramer test using the R language environment as described in Table 1.
Experimental design and statistical analysis
Experimental design details, sample sizes, and results of statistical analyses are presented in the corresponding Materials and Methods and Results sections, figure legends, and Table 1. Mice of either sex were used throughout the study. Exact p-values are listed on Table 1 and also labeled on figures. A p < 0.05 was considered significant for all tests.
Results
Heterozygous yotari mutant mice exhibit a thinner neocortical layer 1 than wild-type mice and a split CA1 pyramidal cell layer in the caudal hippocampus
To examine whether neuronal positioning is disrupted by reduced Dab1 gene expression, we compared the cytoarchitecture of the cerebral neocortex and hippocampus between wild-type and heterozygous yotari mutant mice using Nissl staining at P0 and P7. There were no obvious differences in cortical layer structure between genotypes at P0 (Extended Data Fig. 1-1), but layer 1 was significantly thinner in both the rostral brain region containing somatosensory and auditory cortices (Fig. 1, Rostral) and caudal brain region containing visual and auditory cortices (Fig. 1, Caudal) of P7 heterozygous yotari mutant mice compared with P7 wild-type mice.
To quantitatively evaluate potential differences in neocortical layer thickness between wild-types and heterozygous yotari mutant mice, we stained nuclei in coronal brain sections with DAPI (Fig. 2B–G,B’–G’) and measured layer 1 thickness (Fig. 2I) and cortical wall thickness (Fig. 2K; see Materials and Methods). Consistent with Nissl staining, layer 1 thickness (Fig. 2I) of the rostral brain region containing somatosensory and auditory cortices (Fig. 2; Extended Data Fig. 2-1, Rostral), middle brain region containing visual, parietal, and auditory cortices (Fig. 2; Extended Data Fig. 2-1, Middle), and caudal brain region containing visual and auditory cortices (Fig. 2; Extended Data Fig. 2-1, Caudal) were significantly reduced in yotari heterozygous mice compared with wild-type mice. Moreover, relative layer 1 thickness, calculated as the ratio of layer 1 thickness to cortical wall thickness, was also reduced in rostral, middle, and caudal regions of heterozygous yotari mouse neocortex (Fig. 2J). In the hippocampus as well, layer structure was abnormal in heterozygous yotari mice, but was manifested by an abnormal split in the caudal CA1 pyramidal cell layer (Fig. 1C,C’’).
The RORB-positive layer 4 but not TBR1-positive and BRN2-positive layers 5/6 and 2 are expanded in the neocortex of heterozygous yotari mice
The reduced layer 1 thickness observed among heterozygous yotari mice by Nissl and DAPI staining may result from disruption of the cortical layer structure and/or positional changes in cortical neurons (as seen in reeler and yotari homozygous mutants; Howell et al., 1997; Sheldon et al., 1997; Yoneshima et al., 1997; Trommsdorff et al., 1999). To distinguish between these possibilities, we examined whether relative position or distance from the ventricular surface of cortical neurons is differed between genotypes by staining coronal brain sections with antibodies specifically targeting layer 5/6 neurons (TBR1), layer 4 neurons (RORB), or layer 2/3 neurons (BRN2), and measured individual cell positions (Fig. 3). While there were no significant differences in relative position of cells (Fig. 3C), scatter plots revealed that BRN2-positive layer 2/3 neurons tended to be located farther from the ventricular surface in heterozygous yotari mice than wild-type mice (Fig. 3D), suggesting that layer 2/3 is shifted by cell dispersion into layer 1.
We then examined whether this reduction in layer 1 thickness was accompanied by changes in other layers. Briefly, we defined a 200-μm wide measurement area from the layer 1–2 boundary to the ventricular surface (Extended Data Fig. 3-1A) and examined layer 2/3, 4, and 5/6 area (Extended Data Fig. 3-1B), cell number (Extended Data Fig. 3-1C), and cell density (Extended Data Fig. 3-1D) in somatosensory (Rostral), parietal (Middle), and visual (Caudal) cortices. We found that the measured area was significantly larger in the middle and caudal brain regions of heterozygous yotari mice (Extended Data Fig. 3-1B). There was also a significant increase in cell number within the middle brain region of heterozygous yotari mice, and a numerical increase in both rostral and caudal brain regions that did not reach statistical significance (Extended Data Fig. 3-1C). In contrast, there were no differences in cell density between wild-type and heterozygous yotari mutant mice (Extended Data Fig. 3-1D).
We also examined the thickness of each layer by measuring the distributions of neurons immunopositive for layer-specific markers (Extended Data Fig. 3-2). For comparison, we defined the layer thickness as the length along which the 20th–80th percentile population of layer marker-positive cells is distributed (Extended Data Fig. 3-2A) as it was difficult to clearly determine the boundaries of layers based on the distribution of the all layer marker positive cells. This analysis revealed no significant difference in TBR1-positive layer 5/6 thickness or BRN2-positive layer 2/3 thickness, but RORB-positive layer 4 thickness was greater in heterozygous yotari mice than wild-type mice (Extended Data Fig. 3-2B). These results suggest that a half reduction in the Dab1 gene dosage causes a decrease in layer 1 thickness as well as an increase in the area of the neocortical region from layer 2 to the ventricular surface and cell number with an expansion of RORB-positive layer 4 thickness.
To further explore whether the reduction of DAB1 protein causes superficial positional shift of the superficial layer neurons, late-born neurons were labeled by BrdU at E16.5 and visualized at P9 (Fig. 4A) together with BRN2 immunostaining in somatosensory (Rostral), parietal (Middle), and visual (Caudal) cortices (Fig. 4B,C). There was no significant genotype difference in either mean relative position of BrdU/BRN2 double-positive neurons (Fig. 4B) or mean distance from the ventricular surface (Fig. 4C). Similarly, there were no significant differences in the number of E16.5-BrdU-labeled cells, BRN2-positive cells, and BrdU/BRN2 double-positive cells in a 200-μm-wide region from the layers 1–2 boundary to the ventricular surface between wild-type and heterozygous yotari mice (Extended Data Fig. 4-1).
Extended Data Figure 4-1
There were no significant changes in number and cell density of E16.5-BrdU-labeled cells, BRN2-positive cells, and BrdU/BRN2 double-positive cells between wild-type and heterozygous yotari mice. The numbers of BrdU/BRN2 double-positive cells (A), E16.5-BrdU-labeled cells (C), and BRN2-positive cells (E) were counted in a 200-μm-wide region between the layer 1–2 boundary and ventricular surface (shown in Extended Data Fig. 3-1A). Corresponding cell densities were then calculated (B, D, F). Measurements were performed in the somatosensory (rostral), parietal (middle), and the visual (caudal) cortex of three wild-type (W) and three yotari heterozygous mice (He), each from different litters. The large black dots and black bars in the graph represent the mean and SD, respectively. The gray dots indicate the mean from one brain and each small black dot in the scatter plot is a raw datapoint obtained from the rostral, middle, or caudal brain region. The total number of analyzed brain sections is shown at the bottom of the graph. Groups were compared by Mann–Whitney U test (A) or independent samples Student’s t test (B–F). The p-values are labeled on the graph. ns (not significant) p > 0.05. Download Figure 4-1, EPS file.
Collectively, these results confirmed that a reduction in Dab1 gene expression does not severely disrupt the basic layer structure or layer-specific spatial distribution of neocortical neurons. In addition, these findings suggest that expression from a single Dab1 gene is sufficient for the radial migration of late-born neurons, and that the reduction in layer 1 thickness is not caused by disruption of radial migration of the neocortical neurons.
Abnormal splitting of the CA1 pyramidal layer in heterozygous yotari mice is caused mainly by migration failure of late-born neurons
Nissl staining revealed an abnormally split CA1 cell layer in the caudal hippocampus of P7 heterozygous yotari mice (Fig. 1C,C’’). To investigate the identity of the abnormally located cells and to determine their birthdates, mice were injected with BrdU at E12.5, E14.5, and E16.5, and brains were fixed for co-staining with the CA1 neurons marker CTIP2 at P9 (Fig. 5A). Subsequent analysis revealed that almost all abnormally located cells were CTIP2-positive and most were produced in the late stages of development as shown by greater BrdU staining at E16.5 compared with earlier times (Fig. 5A,B).
We then examined whether splitting of the CTIP2-positive pyramidal cell layer is caused by failure of neuronal migration. According to the GFP-labeling study, E16.5-born neurons labeled with GFP by in utero electroporation migrate within the stratum oriens (SO) around P0–P3 and reach the top of the pyramidal cell layer at P4 (Kitazawa et al., 2014). Thus, we examined the positioning of E16.5 BrdU-labeled and CTIP2-positive cells at P3. We found that more E16.5 BrdU-labeled CTIP2-positive cells were distributed in the SO of heterozygous yotari mice than wild-type mice at P3 (Fig. 5D). These results suggest that abnormal splitting of the pyramidal cell layer occurs during migration within the SO and is already started at P3 (region shown with brackets in Fig. 5C +/yot).
Mispositioned hippocampal pyramidal neurons in heterozygous yotari mice have shorter and misoriented apical dendrites
Previous studies have shown that null mutation of the Dab1 gene impairs pyramidal neuron development and reduces dendritic complexity in the hippocampus (Niu et al., 2004; Matsuki et al., 2008). Here, we examined whether reduced Dab1 also results in abnormal apical dendrite morphology by sparse in utero labeling of hippocampal pyramidal neurons using a DIO reporter system (Sohal et al., 2009; Fig. 6A). Briefly, a reporter gene (GFP or mCherry) is inserted within the DIO cassette in reverse orientation with respect to the promoter, thereby preventing translation under normal conditions but allowing translation in the presence of Cre recombinase as the reporter gene orientation is switched. Further, to induce DIO-mediated reporter gene expression specifically in neurons, Cre expression was controlled by a promoter for Doublecortin (Dcx), a gene thought to be expressed predominantly in newly generated neurons, including hippocampal pyramidal neurons (X. Wang et al., 2007; Yoo et al., 2011). Thus, we co-injected rAAV vectors carrying (1) a Cre recombinase driven by a Dcx promoter (AAV-DCX-Cre), (2) DIO-EGFP driven by a CAG promoter (AAV-CAG-DIO-EGFP), or (3) DIO-mCherry driven by a CAG promoter (AAV-CAG-DIO-mCherry) to wild-type and heterozygous yotari mice at P0, and prepared brain sections at P7. For sparse labeling, concentrations of AAV-CAG-DIO-EGFP and AAV-CAG-DIO-mCherry were appropriately reduced based on prior examinations. We then compared apical dendrite morphology among (1) mCherry-positive cells distributed within the pyramidal cell layer of wild-type mice (Fig. 6B–C’), (2) mCherry-positive cells normally distributed within the pyramidal cell layer of heterozygous yotari mice (Fig. 6E,E’, yellow arrowheads), and (3) mCherry-positive cells ectopically distributed within the SO of heterozygous yotari mice (Fig. 6D–E’, cyan arrowheads). Total apical dendrite length was significantly shorter among ectopically positioned neurons in heterozygous yotari mice compared with wild-types (Fig. 6F). Total apical dendritic length was also numerically reduced among normally positioned neurons in heterozygous yotari mice compared with wild types, although the difference did not reach statistical significance. In addition, the total number of branching points was numerically reduced among normally positioned and ectopically positioned neurons of heterozygous yotari mice compared with wild-types (Fig. 6G). Furthermore, many of the primary apical dendrites in ectopically positioned neurons of heterozygous yotari mice were orientated in abnormal directions, although the average did not differ significantly between genotypes (Fig. 6H).
Neocortical layer 2 pyramidal neurons in heterozygous yotari mice tend to extend their apical dendrites within layer 2 rather than in layer 1
In addition to a migration defect, null mutation, or knock-down of Dab1 causes hypomorphic growth and misorientation of apical dendrites in the cerebral neocortex (Niu et al., 2004; Olson et al., 2006; Franco et al., 2011). Layer 1 is composed mainly of the apical dendrites projecting from pyramidal neurons in layers 2/3 and 5, suggesting that the thinner layer 1 documented here may stem from fewer apical dendrites reaching into the superficial cortex. To simultaneously observe positioning of neurons and the fine dendritic morphology, we introduced a plasmid DNA solution containing pCAGGS-RG, which expresses DsRed under the control of a CAG promoter, and pDCX-Cre into late-born superficial neurons by in utero electroporation at E16.5 (Fig. 7A). Since the DsRed cassette is flanked by two loxP sites, in the presence of Cre, the DsRed cassette is removed and EGFP begins to be transcribed. By reducing the concentration of pDCX-Cre, we sparsely labeled neurons with EGFP. For quantitative analysis of apical dendrite morphology, we first measured the angle between the first branch point (the end of the main shaft of the apical dendrite) and the perpendicular line passing through the center of the cell body (Fig. 7D,E) but found no significant difference between genotypes (Fig. 7E). Similarly, there were no genotype differences in total length and total number of branching points (Fig. 7F,H). However, we found significant differences both in the total length of apical dendrites and the total number of branching points when measured separately within layers 1 and 2 (Fig. 7G,I). Those apical dendrites within layer 1 were significantly shorter and less morphologically complex in heterozygous yotari mutant mice than wild-type mice, while those within layer 2 were longer and more branched in heterozygous yotari mutant mice. These results suggest that a reduction of DAB1 protein disrupts the normal guidance of apical dendrites.
Discussion
We demonstrate that Dab1 is haploinsufficient for the formation of the cerebral neocortex and hippocampus, although specific developmental processes are differentially sensitive to Dab1 gene dose. Neocortical layer 1 is thinner in heterozygous yotari mutant mice than in wild-type mice (Figs. 1 and 2), but layer marker staining and BrdU birth-dating indicated that half reduction of Dab1 does not disrupt radial neuronal migration and layer formation (Figs. 3 and 4). Rather, sparse labeling revealed that layer 2 neurons tend to elongate apical dendrites within layer 2 rather than layer 1 (Fig. 7). Heterozygous yotari mutant mice also exhibit splitting of the CA1 pyramidal cell layer in the caudal hippocampus (Figs. 1 and 5), and BrdU birth-dating studies suggests that this developmental abnormality is caused, at least in part, by migration defects of late-born pyramidal neurons (Fig. 5). Further, sparse labeling revealed that mispositioned neurons in the SO tend to have short and abnormally orientated dendrites (Fig. 6). Thus, reduction of Dab1 gene expression has distinct effects on development of the cerebral neocortex and hippocampus.
Homozygous null mutation of the Dab1 gene causes severe migration failure both in the cerebral neocortex and hippocampus (Howell et al., 1997; Yoneshima et al., 1997; Franco et al., 2011; Blume et al., 2017). However, half reduction of Dab1 gene expression did not cause migration failure in the cerebral neocortex, and impaired only late-born pyramidal neuron migration in the caudal hippocampus. This regional difference may be attributable to distinct migration modes. Early-born pyramidal neurons in the cerebral neocortex migrate radially by somal translocation, while late-born neurons follow a series of steps including multipolar migration, locomotion, and terminal translocation (Noctor et al., 2001; Tabata and Nakajima, 2003; Tabata et al., 2009; Sekine et al., 2011; Ohtaka-Maruyama and Okado, 2015; Kon et al., 2017). In contrast, hippocampal pyramidal neurons migrate radially by a climbing mode after a multipolar stage (Kitazawa et al., 2014; Hayashi et al., 2015). These various migration modes may have different DAB1 protein sensitivities (requirements), although it is still unclear why only the caudal hippocampus is affected as migration modes are similar across the hippocampal rostral–caudal axis. It may be that a subpopulation of pyramidal neurons in caudal hippocampus, possibly comprising a specific sublayer, requires both Dab1 gene copies for proper migration in the caudal hippocampus (Soltesz and Losonczy, 2018).
Homozygous null mutation of the Dab1 gene also causes complete loss of cerebral neocortex layer 1 (Sweet et al., 1996; Howell et al., 1997; Ware et al., 1997; Yoneshima et al., 1997). In contrast to this severe phenotype, half reduction of Dab1 gene resulted only in a thinner layer 1 throughout the rostral–caudal axis (Figs. 1 and 2). This abnormality appears to arise after radial neuronal migration is complete, as no obvious anatomic change in neocortical structure was observed at P0 (Extended Data Fig. 1-1), migration failure was not observed (Fig. 4), and layer-specific marker staining showed no positional changes in TBR1-positive, RORB-positive, or BRN2-positive cells (Fig. 3). These effects of heterozygous Dab1 yotari mutation further suggest two distinct signaling pathways for excitatory neurons with differential dependencies on Dab1 gene dosage: a less sensitive pathway for radial neuronal migration and a sensitive pathway for the control of layer 1 thickness.
Several other mutant mouse lines exhibit abnormal and ectopic positioning of cells in cortical layer 1, including Dab1p45/– (Herrick and Cooper, 2002), Reelin ΔC-KI (Kohno et al., 2015), Vldlr KO (Hack et al., 2007; Hirota and Nakajima, 2020), Apoer2 KO (Hirota et al., 2018), DeltaNp73Cre;ROSA26stop-dt-a (de Frutos et al., 2016), Emx1Cre;Wnt3adt-a (de Frutos et al., 2016), Emx1-Cre;Pcdh-γfcon3/fcon3 (Garrett et al., 2012), Plexin A2/A4 double KO (Hatanaka et al., 2019), and Sema6a KO (Hatanaka et al., 2019) mice. Many of these lines harbor mutations in the Reelin signaling pathway. The Dab1p45 gene encodes a hypomorphic alternative splice variant of Dab1, and Dab1p45/– mice also show thinning of layer 1 and splitting of the CA1 pyramidal layer (Herrick and Cooper, 2002). Since the Dab1p45 allele was created by the insertion of Dab1p45 cDNA, only the DAB1 p45 protein is expressed in this mouse (Herrick and Cooper, 2002). Therefore, the Dab1p45/– condition is distinct from the heterozygous null condition, and analysis of heterozygous Dab1 mutant mice such as yotari could be particularly valuable for elucidating the pathomechanisms of human diseases associated with heterozygous Dab1 loss-of-function mutations. While it is yet unknown whether haploinsufficiency of the Dab1 gene has similar effects on human cortical and hippocampal development, haploinsufficiency of the related Reelin gene is associated with various neuropsychiatric disorders (Guidotti et al., 2016; Ishii et al., 2016; Hattori and Kohno, 2021). Given the critical role of Dab1 in the canonical Reelin signaling pathway (Honda et al., 2011; Jossin, 2020), the findings of this study could provide potential explanations for human diseases.
Similar to heterozygous yotari mice, Reelin ΔC-KI mouse lacking the C-terminal region of Reelin, which is thought to be required for efficient DAB1 phosphorylation (Nakano et al., 2007), gradually show thinning of layer 1 and apical dendrites with abnormal orientations that tend to avoid growing into the layer 1 (Kohno et al., 2015). The DeltaNp73Cre;ROSA26stop-dt-a and the Emx1Cre;Wnt3adt-a mouse lines also underexpress Reelin protein and show reduced layer 1 thickness and cortical neurons with decreased apical dendrite tuft complexity (de Frutos et al., 2016). Since the Reelin receptor VLDLR is predominantly expressed in the superficial cerebral neocortex (Hirota et al., 2015), Vldlr KO would result in deficient Reelin signaling within superficial cortical layers, and indeed Vldlr KO also results in ectopic positioning of CUX1-positive or CUX2-positive neurons in layer 1 (Hack et al., 2007; Hirota and Nakajima, 2020). Therefore, it appears that partial attenuation of the Reelin-Dab1 signaling pathway does not inhibit radial neuronal migration in the cerebral neocortex but causes a decrease in the layer 1 thickness.
We suggest that these deficits account at least in part for the reduction in layer 1 thickness. Layer 1 is composed of the apical dendrites of excitatory neurons within layers 2/3 and 5, local interneurons, and Cajal-Retzius cells, thalamic axons, and axons of excitatory neurons projecting from other layers and regions of the cerebral neocortex (Solari and Stoner, 2011; Harris and Shepherd, 2015; Genescu and Garel, 2021; Gesuita and Karayannis, 2021; Schuman et al., 2021). However, there is a particularly strong association between apical dendrites and layer 1 thickness as evidenced by (1) the tendency of apical dendrites projecting from sparsely labeled superficial layer neurons in heterozygous yotari mice to elongate within layer 2 but not layer 1 (Fig. 7), (2) the reduced layer 1 thickness and associated apical dendrite maldevelopment among Reelin C-terminal deletion mutant mice (Kohno et al., 2015), Reelin-producing Cajal-Retzius cell-ablated mice (de Frutos et al., 2016), and Dab1 knock-down or knock-out experiments (Olson et al., 2006; Franco et al., 2011; Sekine et al., 2011, 2012), and (3) defects in apical dendrite morphology observed in other mutants with abnormal layer 1 reduction, such as Emx1-Cre;Pcdh-γfcon3/fcon3 (Garrett et al., 2012), PlexinA2/A4 double KO, and Sema6a KO mice (Hatanaka et al., 2019). Thus, abnormalities in the apical dendrites of neocortical pyramidal neurons may explain reduced layer 1 thickness. It is also possible that the increase in cell number within other neocortical layers (Extended Data Fig. 3-1) leads to a reduction in layer 1 thickness. The causal relationship between the increase in cell number and the expansion of the neocortical area, and why the cell number increased, are important subjects for future study.
In addition to abnormal development of neocortical layer 2 apical dendrites, we also observed a decrease in the total length of apical dendrites projecting from ectopic pyramidal neurons of the caudal hippocampus as well as misoriented apical dendrites among subset of neurons. While it is possible that the half reduction of DAB1 protein directly caused these phenotypes (Niu et al., 2004; Matsuki et al., 2008), it is also possible that these abnormalities are secondary to migration failure or that neurons ectopically distributed in the SO region are still migrating and these misoriented structures are leading processes rather than dendrites. Regardless of mechanism, the resulting maldevelopment in the hippocampus may contribute to a myriad of neuropsychiatric diseases exhibiting deficits in learning and memory formation (Hainmueller and Bartos, 2020).
Similar to the Dab1 gene, the Reelin gene is also known to exhibit haploinsufficiency (Liu et al., 2001; Carboni et al., 2004; Qiu et al., 2006; Nullmeier et al., 2011; Bouamrane et al., 2016; L. Wang et al., 2019), and reduced Reelin levels have been observed in the brains of patients with schizophrenia (Impagnatiello et al., 1998; Guidotti et al., 2000), bipolar disorder (Guidotti et al., 2000), autism (Fatemi et al., 2001, 2005), and Alzheimer’s disease (Chin et al., 2007). Moreover, heterozygous reeler mice and Reelin ΔC-KI mice exhibit similar neurochemical, behavioral, and cognitive abnormalities to those observed in autism or schizophrenia (Costa et al., 2001; Tremolizzo et al., 2002; Carboni et al., 2004; Nullmeier et al., 2011; Sakai et al., 2016; Sánchez-Hidalgo et al., 2022). The anatomic abnormalities we observed in heterozygous yotari mice may also be related to behavioral defects because layer 1 and hippocampus are essential for higher-order cognitive functions (Hainmueller and Bartos, 2020; Schuman et al., 2021). Consistent with this possibility, behavioral abnormalities observed in schizophrenia, such as hyperactivity, decreased anxiety-like behavior, and impaired working memory, have been observed in dorsal forebrain-specific Dab1 conditional knock-out mice (Imai et al., 2017). These findings suggest that therapeutic enhancement of Reelin signaling could prevent or ameliorate the cognitive deficits of schizophrenia and other neurodevelopmental disorders (Guidotti et al., 2016; Ishii et al., 2016; Hattori and Kohno, 2021). Indeed, administration of Reelin to mouse models of cognitive dysfunction such as Angelman syndrome and schizophrenia could prevent or ameliorate the cognitive deficits. (Hethorn et al., 2015; Ishii et al., 2015; Ibi et al., 2020; Sawahata et al., 2021). If a similar neuropsychiatric abnormality is observed because of Dab1 haploinsufficiency, it may be possible to ameliorate the defect by activating the Dab1 signaling through Reelin supplementation.
In conclusion, while the radial migration of neocortical neurons can proceed normally when Dab1 is reduced to half, Dab1 is haploinsufficient for maintaining layer 1 thickness and for the migration of late-born pyramidal neurons in the caudal hippocampus. Further studies are required to explore why these developmental events are differentially dependent on Dab1 expression level. We suggest that elucidation of the Reelin signaling pathways may provide clues to neuropsychiatric disease pathogenesis and potential therapeutic strategies.
Acknowledgments
Acknowledgments: We thank Dr. Jun-ichi Miyazaki (Osaka University, Osaka, Japan) for providing the pCAGGS, Dr. Q. Lu (Beckman Research Institute of the City of Hope) for the Dcx4kb-EGFP, Dr. S. Miyagawa (Osaka University) for the pCXN-Cre, Karel Svoboda (Howard Hughes Medical Institute) for AAV-CAG-GFP, University of Pennsylvania Penn Vector Core for pAdΔF6 and pAAV2/9, technicians of Nakajima laboratory, especially Yuki Nomura for her help in analyzing data, and all the members of the Nakajima laboratory for their helpful discussion and valuable comments.
Footnotes
The authors declare no competing financial interests.
This work was supported by JSPS KAKENHI Grant Numbers JP21K06413, JP18K06508, JP24700357, JP16H06482, JP20H05688, JP22K19365, and JP20K06670; Asahi Glass Foundation; Takeda Science Foundation; Keio Gijuku Fukuzawa Memorial Fund for the Advancement of Education and Research; Keio Gijuku Academic Development Funds; Gifu Pharmaceutical University Research Fund; and Toukai Foundation for Technology.
This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license, which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.