WDFY3 mutation alters laminar position and morphology of cortical neurons
Molecular Autism volume 13, Article number: 27 (2022)
Proper cerebral cortical development depends on the tightly orchestrated migration of newly born neurons from the inner ventricular and subventricular zones to the outer cortical plate. Any disturbance in this process during prenatal stages may lead to neuronal migration disorders (NMDs), which can vary in extent from focal to global. Furthermore, NMDs show a substantial comorbidity with other neurodevelopmental disorders, notably autism spectrum disorders (ASDs). Our previous work demonstrated focal neuronal migration defects in mice carrying loss-of-function alleles of the recognized autism risk gene WDFY3. However, the cellular origins of these defects in Wdfy3 mutant mice remain elusive and uncovering it will provide critical insight into WDFY3-dependent disease pathology.
Here, in an effort to untangle the origins of NMDs in Wdfy3lacZ mice, we employed mosaic analysis with double markers (MADM). MADM technology enabled us to genetically distinctly track and phenotypically analyze mutant and wild-type cells concomitantly in vivo using immunofluorescent techniques.
We revealed a cell autonomous requirement of WDFY3 for accurate laminar positioning of cortical projection neurons and elimination of mispositioned cells during early postnatal life. In addition, we identified significant deviations in dendritic arborization, as well as synaptic density and morphology between wild type, heterozygous, and homozygous Wdfy3 mutant neurons in Wdfy3-MADM reporter mice at postnatal stages.
While Wdfy3 mutant mice have provided valuable insight into prenatal aspects of ASD pathology that remain inaccessible to investigation in humans, like most animal models, they do not a perfectly replicate all aspects of human ASD biology. The lack of human data makes it indeterminate whether morphological deviations described here apply to ASD patients or some of the other neurodevelopmental conditions associated with WDFY3 mutation.
Our genetic approach revealed several cell autonomous requirements of WDFY3 in neuronal development that could underlie the pathogenic mechanisms of WDFY3-related neurodevelopmental conditions. The results are also consistent with findings in other ASD animal models and patients and suggest an important role for WDFY3 in regulating neuronal function and interconnectivity in postnatal life.
During prenatal neurogenesis, newly born neurons are deployed from proliferative compartments surrounding the ventricles towards the surface of the brain where they will settle into their proper laminae and nuclei to form functional circuits [1,2,3,4]. Any disturbance of this tightly orchestrated process can result in neuronal migration disorders (NMDs), birth defects with often devastating consequences for affected individuals [5,6,7,8,9]. In the cerebral cortex, the primary site of many NMDs, projection neurons are generated in the proliferative layers of the ventricular and subventricular zones (VZ, SVZ) from radial glial cells (RGCs) and migrate radially to the pial surface along radial glial fibers [2, 10, 11]. The earliest born neurons building the cortical plate become the deepest layer of the cortex, and ensuing waves of newly born neurons migrate past their predecessors forming layer upon layer below the superficial marginal zone. NMDs can range from relatively mild cases such as focal cortical dysplasia with only limited migration anomalies [12,13,14] to crippling diseases like lissencephaly, which display a reduction in gyri and a thickening of gray matter over a large area or the entire cerebral cortex often leading to death in early childhood [15, 16].
Interestingly, there is a pronounced comorbidity of NMDs with autism spectrum disorders (ASDs), a fact that may contribute to the high prevalence of epilepsy in ASDs with co-diagnoses reaching up to 40% of autism cases [17, 18]. Indeed, evidence from postmortem analysis of human cases [19,20,21,22,23] and more recently from ASD mouse models [24, 25] strongly suggests that focal cortical dysplasia is a common feature of the autistic brain, further strengthening the notion that a subset of autism cases is rooted in dysregulations of prenatal neurogenesis.
We previously reported on the generation of Wdfy3 mutant mice and on the neurodevelopmental anomalies associated with Wdfy3 loss-of-function . WDFY3 is an autism risk gene and causative in a range of other often comorbid neurodevelopmental disorders, including neurodevelopmental delay, intellectual disability, and hypotonia [26,27,28,29,30,31,32]. Depending on the specific allele, heterozygous WDFY3 loss in humans is typically associated with macrocephaly, a feature faithfully modeled in heterozygous and homozygous mutant mice. Wdfy genes (1–4) encode a small family of four neuronally expressed, intracellular molecules associated with vesicular transport. WDFY3 is a member of the BEACH (beige and CHS proteins) protein family and contains, in addition to the BEACH domain, a PH domain, five WD40 domains, and a C-terminal FYVE (Fab1/YOTB/Vac1/EEA1) domain enabling its integration into vesicular membranes [32, 33]. Human WDFY3 has been shown to act as an autophagosomal scaffolding protein required for the selective recruitment and degradation of macromolecular components such as aggregation-prone proteins. Untangling the cellular causes underlying the morphological abnormalities in Wdfy3 mutant mice, we identified an essential role for WDFY3 in regulating a subset of RGC divisions by promoting asymmetric differentiative over symmetric proliferative divisions. Consequently, WDFY3 loss promotes self-renewing divisions in RGCs, transient expansion of the progenitor pool, increased neurogenesis, and larger brains. In addition, homozygous Wdfy3 mutation in mice results in focal neuronal migration defects in which deep-layer neurons are typically mispositioned within upper layers.
Radial migration of newly born cortical neurons is regulated by factors both intrinsic to the migrating cells and/or extrinsic, determined by the environment they navigate. Cell autonomous effectors may involve cytoskeletal proteins that control cellular motility, such as neuronal tubulins  and their regulators [34, 35], while cell non-autonomous effects may involve chemoattractants expressed by surrounding cells, such as reelin . Earlier findings of WDFY3’s cytosolic and possibly vesicular localization within dividing progenitors suggested cell autonomous control over neuronal migration, but in the absence of functional evidence remained speculative.
Here, we employed mosaic analysis with double markers (MADM) [37, 38] to trace the laminar positioning of cortical projection neurons depending on Wdfy3 genotype and quantitatively assess mispositioned cells to endogenous Wdfy3 loss. In addition, we addressed the question of whether Wdfy3 mutation has effects on neuronal morphology and circuit integration by evaluating dendritic arborization, spine density, and morphology. Our results demonstrate a predominantly cell autonomous role of Wdfy3 in regulating either aspect of neuronal laminar and circuit integration.
Animal husbandry and genotyping
Animals were housed in Plexiglas cages (55 × 33 × 19 cm) and maintained in facilities approved by the Association for Assessment and Accreditation of Laboratory Animal Care (AAALAC) International under standard laboratory conditions (21 ± 2 °C; 55 ± 5% humidity) on a 12-h light/dark cycle, with ad libitum access to both water and standard rodent chow. Animal handling protocols were approved by the University of California at Davis Institutional Animal Care and Use Committee overseen by the AAALAC International accreditation program (latest accreditation in February 14, 2020) and in compliance with the ARRIVE  and NIH guidelines .
Mice carrying the Wdfy3lacZ (Wdfy3tm1a(KOMP)Mbp) allele were generated and genotyped as previously described  and maintained on C57BL/6NJ background. To generate Wdfy3lacZ-MADM-5 mice, we crossed Wdfy3+/lacZ with homozygous MADM-5GT/GT mice [38, 41]. Subsequently, compound heterozygous offspring was crossed with homozygous MADM-5TG/TG; Emx1Cre mice. The resulting offspring was genotyped for the presence of Wdfy3lacZ, Emx1Cre, GT, and TG cassettes and once pups positive for all markers identified they were processed for histological analysis.
Perfusion, brain collection, tissue preservation, and sectioning
Mice collected for histological analysis were transcardially perfused with phosphate-buffered saline (PBS) followed by 4% paraformaldehyde (PFA) in PBS using a medical pump at a rate of 1 mL/min and immersed in 4% PFA/PBS overnight at 4 °C. The next day skulls were washed with PBS, brains excavated, cryopreserved through gradual sucrose immersion (15%, 30% in PBS), embedded in O.C.T. compound (Fisher Healthcare), flash frozen on dry ice, and then transferred to − 80 °C until sectioning. O.C.T.-embedded brains were coronally sectioned on a Leica CM-1950 cryostat at a thickness of 30 μm. Sections were transferred onto glass slides (Thermo Fisher Scientific, Waltham, MA) and placed on a slide warmer (Premiere XH-2004) at 40 °C for a minimum of 36 h. Slides from the warmer were then placed in a slide box and frozen at − 80 °C until used.
Alternatively, dissected brains used for analysis of cell morphology were transferred into PBS/0.1% NaN3 until embedded in 4% agarose/PBS for immediate vibratome sectioning. Agarose-embedded brains were submerged in ice-cold PBS and coronally sectioned at 200 μm on a Leica VT-1000S vibratome. Sections were then transferred into 24-well plates with 1 mL of PBS on ice and stored at 4 °C until use.
Slide-mounted sections were fixed with 4% PFA/PBS, followed by three washes in PBS for 5 min each. For antigen retrieval, slides were submerged in 1 × working concentration of Diva Decloaker solution (Biocare Medical, Pacheco, CA) while heated in a Decloaking Pressure Chamber at 90 °C for 30 min. Following this treatment, slides were washed in PBS three times for 5 min and tissue blocked with 10% donkey serum in PBS with 1% Triton X-100 (PBST) for 1 h at room temperature. After being washed with PBS, sections were incubated with chosen primary antibodies for 18 h at 4 °C, subsequently washed five times for 10 min each with PBS, and secondary antibodies applied for 2 h at room temperature. A list of all antibodies used in this study is provided in Table 1. After secondary incubation, slides were washed five times for 10 min each in PBS and for 5 min submerged in a 0.1% 4′,6-diamidino-2-phenylindole (DAPI) (Thermo Fisher Scientific) to create a nuclear counterstain. After being washed twice 10 min each with PBS, sections were covered with Fluoromount-G (Thermo Fisher Scientific), mounted with cover glass (Thermo Fisher Scientific), and then allowed to dry in a dark ventilated area for a minimum of 20 h before being imaged.
Agarose-embedded sections (200 μm) were immunostained free floating on an orbital shaker (Labnet) at 30 rpm using the process described above while omitting antigen retrieval. Transfer of sections into washes, blocking, and incubation liquids was done using a paintbrush. Antibody incubation times were extended to account for increased tissue thickness compared with cryosections with primary incubation applied for 64 h at 4 °C and secondary incubation for 18 h at 4 °C. Sections were finally washed 6 times for 45 min each.
To increase clarity and microscope laser penetration, post immunostaining 200-μm-thick sections used for cell morphological analysis were immersed overnight at 4 °C in ultrafast optical clearing method (FOCM) solution; 30% (wt/vol) urea, 20% (wt/vol) D-sorbitol, and 5% (wt/vol) glycerol were dissolved in dimethyl sulfoxide (DMSO) . The next day the sections were removed from the FOCM solution, transferred to microscope slides (Thermo Fisher Scientific), and surrounded with Fluoromount-G (Thermo Fisher Scientific) before coverslipping. Mounted slides were set to dry in the dark for a minimum of three days before imaging.
Microscopy and image processing
To assess genotype distribution and laminar positioning, fluorescent confocal microscopy and image acquisition was performed on an inverted Nikon C1 microscope with associated software, typically at 10 × magnification. Images used for cell morphological analysis, including dendritic arbor and spine analysis, were obtained on a Nikon A1 microscope with associated software. In either case, we imaged sections under varying acquisition parameters for best possible comparisons of the two fluorophores (GFP and tdT). Subsequently, images were uploaded and cells were counted in NIS-Elements. Cellular morphology was analyzed in FIJI using the Simple Neurite Tracer plugin where further adjustments ensured even brightness and detail of GFP and tdT in grayscale.
All data were processed, analyzed, and graphical figures created in GraphPad Prism 9. To remove suspected outliers, data were initially processed with the ROUT outlier test (Q value = 1%). Applicable data were then tested for normality, variance, and standard deviation before analysis. Lamination analysis was done using Fisher’s exact test, population analysis was done using an unpaired Student’s t test, and Sholl profile analysis was statistically evaluated by two-way analysis of variance (ANOVA), followed by a Tukey’s multiple comparison test. Analysis of bouton density was done using a one-way ANOVA with multiple comparisons, while bouton subtype distribution analysis was done using a two-way ANOVA with multiple comparisons. When applicable and reasonable to add, data are reported as a mean with standard deviation. Bar graphs with replicates are reported as mean with standard error of the mean. Results were considered to be statistically significant if p ≤ 0.05. Individual data points largely correspond to replicates of one brain. The extent of significance between groups is indicated with one, two, three, or four asterisks if p values were equal to or less than 0.05, 0.01, 0.001, or 0.0001, respectively.
Wdfy3 mutant cells get disproportionally lost during early postnatal development in MADM brains
To test the efficacy of the Wdfy3-MADM system, we obtained brains of mice carrying all four transgenic modifications (Wdfy3lacZ, MADM-5GT, MADM-5TG, and Emx1Cre) and immunolabeled for endogenous tdT and eGFP expression. We opted to analyze three male and three female brains each at two developmental stages, postnatal day (P)8, a time point at which developmental neurogenesis is completed and all cortical layers are formed and distinguished, and P30, a time point in which cell death associated with activity-dependent neuron elimination has largely faded . After imaging, sections of both stages revealed the expected sparse labeling of neurons and astrocytes within all cortical layers with cells distributed across labels associated with the anticipated recombination of fluorophore encoding cassettes. We detected red tdT+ cells expected to be wild type (WT or + / +), green homozygous mutant GFP+ cells (lacZ/lacZ), and yellow heterozygous tdT+/GFP+ cells (+ /lacZ) (Fig. 1A, B). We proceeded to count all cells in ten sections each of 3 male and 3 female brains per stage and converted the added numbers of each genotype to ratios of a whole. In total, we counted at P8 685 tdT+ cells, 1,330 tdT+/GFP+ cells, and 658 GFP+ cells and at P30 614 tdT+ cells, 1,236 tdT+/GFP+ cells, and 262 GFP+ cells. Genotype populations were compared across each other and between the two time points (P8 and P30). In P8 brains, the distribution of GFP labeled homozygous mutant neurons was ~ 24%, likely due to cell death significantly decreasing to ~ 12% in P30 brains (0.2368 ± 0.029 at P8 and 0.1216 ± 0.028 at P30; p < 0.0001) (Fig. 1C). This relative decrease in mutant cells at P30 was accompanied by a significant relative increase in both WT (0.2634 ± 0.026 at P8 and 0.2949 ± 0.010 at P30; p = 0.02) and + /lacZ cells (0.4998 ± 0.024 at P8 and 0.5835 ± 0.026 at P30; p = 0.0002) (Fig. 1C). No significant differences between the sexes were observed.
To unveil direct evidence of increased apoptotic rates of lacZ/lacZ neurons, we performed cleaved CAS3 immunofluorescence at P8 with the aim to establish ratios of tdT+/CAS3+, GFP+/CAS3+, or tdT+/GFP+/CAS3+ cells and derive from these ratios statistically significant differences between genotypes. Labeling 40 sections of six brains with this approach and analyzing a total of ~ 1,800 neurons of all genotypes, we did not discover a single cell in which the GFP or tdT signal would overlap with the CAS3 signal. The sparsity of MADM-labeled cells in combination with the rarity of cell death in the postnatal brain presented a barrier to discovery, an inherent limitation of the approach used here (Additional file 1: Fig. S1).
Cortical lamination errors occur more frequently in Wdfy3 homozygous mutant cells
Realizing that early postnatal development is associated with the disproportionate loss of lacZ/lacZ neurons, we proceeded to examine whether this loss is associated with the laminar positioning of these neurons. We performed immunofluorescent labeling targeting endogenous tdT and eGFP expression produced as a result of the MADM system, and colabeled for cortical layer markers TBR1, CTIP2, and BRN2, each at a time. TBR1 is predominantly expressed by layer VI neurons, CTIP2 by layer V neurons, and BRN2 by layer II/III neurons, respectively. The approach allowed us to identify neurons of either fluorescent marker and genotype and, if positive for one of the cortical layer markers, assess whether they were also correctly positioned within their respective layer.
Analyzing a total number of 730 cells of six brains (3 male, 3 female) at P8, mutant neurons in Wdfy3-MADM mice were mispositioned more often than their WT counterpart. TBR1+ (VI) mutant cells were correctly positioned in 46% of instances compared to ~ 91% of WT cells (WT 90.66%; lacZ/lacZ, 46%; p < 0.0001) (Fig. 2A, G). Mutant CTIP2+ cells (V) were also significantly less often correctly positioned compared with WT cells (WT, 78.79%; lacZ/lacZ, 29.76%; p < 0.0001) (Fig. 2C, H). The same observation was made with BRN2+ cells (II/III) with mutant neurons being significantly less often correctly located compared with WT neurons (WT, 89.78%; lacZ/lacZ, 53.91%; p < 0.0001) (Fig. 2E, I). Examining layer-specific positioning of 735 neurons in six Wdfy3-MADM brains at P30 (3 male, 3 female), we found similar, but less exaggerated discrepancies between WT and lacZ/lacZ neurons in P8 brains. TBR1+ mutant cells were aligned in the correct layer in ~ 63% of instances, compared to 87% for their WT counterpart (WT, 86.81%; lacZ/lacZ, 62.82%; p = 0.0003) (Fig. 2B, G). CTIP2+ mutant cells showed correct layer position ~ 23% of the time compared with of ~ 65% WT cells (WT, 65.06%; lacZ/lacZ, 23.38%; p < 0.0001) (Fig. 2D, H). Lastly, BRN2+ mutant neurons were placed in the correct layer ~ 68% of the time, compared with ~ 89% of the time for WT cells (WT, 88.52%; lacZ/lacZ, 67.98%; p < 0.0001) (Fig. 2E, I).
Comparing the two stages (P8 and P30) with each other, the proportion of correctly aligned TBR1+ mutant neurons increased from ~ 46% at P8 to ~ 63% at P30 (p = 0.0338). Comparably, BRN2+ mutant neurons with correct layer position increased from ~ 54% in P8 brains, to ~ 68% in P30 brains (p = 0.0095), while CTIP2+ mutant neurons did not show significant differences between the two developmental time points (p = 0.3787). Also, no significant differences were present for WT neurons of any marker between the two stages and no significant deviations between the sexes were recorded. Combining the data of all three markers together, ~ 45% of mutant neurons showed correct laminar placement in P8 brains, increasing to ~ 59% in P30 brains (P8, 44.87%; P30, 58.58%; p = 0.0003). This relative increase likely reflects the selective elimination of mispositioned lacZ/lacZ neurons during the examined postnatal period.
Wdfy3 mutant pyramidal neurons show decreased dendritic arborization
Neurodevelopmental disorders, including autism, have been linked with changes in neuronal morphology [44, 45]. In the MADM system, sparse labeling and the fact that cells are entirely filled by the either fluorescent protein (eGFP or tdT) provided us with the option to examine cellular morphological parameters and compare them between genotypes. We focused in particular on the complexity of dendritic arbors and the density of dendritic spines or boutons. To capture labeled neurons to the greatest possible extent, 200-μm-thick sections of three male and three female agarose-embedded Wdfy3-MADM brains were prepared and immunostained for eGFP and tdT applying a free-floating protocol. Subsequently, sections were cleared by means of the ultrafast optical clearing method (FOCM)  and mounted on microscope slides for deep imaging of labeled cortical pyramidal neurons. Acquired z-stacks were traced to create Sholl profiles of 81 upper layer (II/III) and 81 deep-layer (V/VI) cortical projection neurons of either genotype.
For either upper layer or deep-layer neurons, both WT and heterozygous cells were more complex than homozygous mutant cells, while no noteworthy differences were observed between sexes. For each genotype, the mean number of intersections peeked at ~ 100 μm from the soma and significant differences between genotypes were mostly confined to distances of 15 μm to 200 μm from the soma. The number of maximum intersections for deep-layer WT neurons at ~ 100 μm was in average 18 while lacZ/lacZ neurons showed a peak average of 13.7 intersections (p < 0.05—0.0001). Heterozygous neurons were slightly but not significantly more complex than WT with an average of 19.5 intersections at ~ 100 μm distance from the soma (+ /lacZ vs. lacZ/lacZ, p < 0.05—0.001) (Fig. 3D). In upper layers, WT and + /lacZ neurons showed near identical Sholl profiles with peak intersections at ~ 100 μm being on average 22.3, while lacZ/lacZ neurons showed a highest average of 15.7 intersections at ~ 100 μm from the soma (WT or + /lacZ vs. lacZ/lacZ, p < 0.05—0.0001) (Fig. 3E).
Wdfy3 mutation increases synaptic spine density
To assess dendritic spine density and subtype distribution, 200-μm sections of four male and three female brains were immunolabeled and cleared via FOCM as described above. From acquired z-stacks of 42 labeled cells of each genotype, 10-μm (in z-depth) regions of dendrites were then extracted and dendritic spines manually counted while also tracking morphological features associated with widely described subtypes, including filopodia, long thin, thin, stubby, mushroom, and branched .
A decrease in Wdfy3 gene dosage led to an overall increase in spine density. In WT neurons, spine density averaged 0.29 spines per 1 μm. In heterozygous neurons, density significantly rose to an average of 0.42 spines per μm. In lacZ/lacZ cells, we recorded an even higher density of 0.54 spines per μm, while no significant sex- or layer-specific differences (II/III vs. V) were observed in any genotype leading us to combine all counts. (WT, 0.2907 ± 0.1111; + /lacZ, 0.4188 ± 0.1004; lacZ/lacZ 0.5409 ± 0.1609; WT vs. + /lacZ p = 0.0286; WT vs. lacZ/lacZ p = 0.0011; + /lacZ vs. lacZ/lacZ; p = 0.0287.)
Dendritic spines can be distinguished by morphological criteria that associate with their state of maturation and function. A well-replicated rank order of maturation distinguishes spines that are filopodia, long thin, thin, stubby, mushroom shaped, and branched . While assessing overall spine density, we traced these different subtypes and converted counts to a percentage of a whole to analyze distribution across genotypes. No significant differences between genotypes were detected from a two-way analysis of variance (ANOVA), followed by a Tukey’s multiple comparison test (Fig. 4H). Irrespective of genotype, the most prevalent spine types were long thin, thin, and mushroom, comprising ~ 90% of all counted spines. In summary, we found Wdfy3 mutation to affect spine density in a gene dosage-responsive manner, but to have no effect on spine subtype distribution.
Human WDFY3 mutation will result in neurodevelopmental disorders, such as intellectual disability, neurodevelopmental delay, and most frequently ASD associated with macrocephaly [26,27,28,29,30, 32]. While the megalencephalic phenotype has been faithfully replicated in both heterozygous and homozygous Wdfy3 mutant mice and its origins traced to altered RGC divisions [24, 32], the developmental origins of the focal neuronal migration defects of homozygous mutant mice remain uncertain. Migration defects of cortical projection neurons can have either intrinsic causes or be the result of disrupted cell–environment interactions [4, 47, 48]. For instance, predominantly intrinsic mechanisms underlie lissencephaly caused by mutations in LIS1 (PAFAH1B1) or DCX [49, 50]. LIS1 interacts with cytoplasmic dynein, a molecular motor required for microtubule-based transport [51, 52]. Lis1 knockdown in mice will block interkinetic nuclear migration, a requirement for both RGC divisions and the motility of neuronal precursors that will derive from them . Similarly, mutations in neuronal tubulin TUBA3 (Tuba1 in mice) will cause classic lissencephaly . Migration defects rooted in the loss of extrinsic signals have been extensively studied in reeler mice that carry a mutation in reelin (Reln) [54, 55], a secreted glycoprotein produced by Cajal–Retzius cells occupying the marginal zone during development . Reln signaling will activate an intracellular cascade of pathways that will promote migration and correct positioning of cortical neurons critical components of which are Notch signaling  and cadherin-mediated cell adhesion .
During mouse cortical neurogenesis, WDFY3 is predominantly expressed by dividing RGCs, the pia, and diffusely in the intermediate zone . By using Emx1Cre as a driver, this study sought to isolate cell autonomous vs. non-cell autonomous effects of WDFY3 loss on radial migration by conditionally limiting mutation to cortical projection neurons and their progenitors. Our results strongly suggest the predominance of intrinsic WDFY3 activity guiding migration and possibly survival of cortical neurons in the developing cortex. Indeed, while at P8 WT cells overwhelmingly localize within the cortical layers specific to their marker identity, lacZ/lacZ cells do so at significantly reduced frequency, typically at only half the rate of WT cells. At P30, we found the ratio of lacZ/lacZ neurons that are correctly positioned within layers VI (TBR1+) and II/III (BRN2+) to be elevated compared to P8. This finding very likely points to corrective mechanisms that are designed to eliminate misplaced and/or inadequately circuit-integrated neurons from the developing brain via programmed cell death. Indeed, TUNEL analysis of mouse parietal neocortical fields 1, 3, and 40 (largely corresponding to the area examined in the present study) revealed cell death to peak at P4, but also to continue well into the third postnatal week . These findings are supported by earlier studies that recorded in rodents an average of 30% neuronal loss during the first postnatal month . More recent work, however, describes a lesser, approximately 12% loss of excitatory neocortical neurons, confined to the period between P2 and P5 . In light of our own data, it is conceivable that this period may be extended in case pathological dysregulations require expanded corrections. That early postnatal neuronal apoptosis is controlled by electrical activity appears also strongly supported by elegant work that associated electrical activity patterns with neuronal death rates in two distinct cortical areas (M1 and S1) . Thus, the suggestion that misplaced neurons may be especially targeted for removal as their incorporation into physiologically functional networks is less likely, plausibly explains our observations of progressively fewer lacZ/lacZ neurons located outside their correct layers.
Homozygous mutant neurons display reduced dendritic arborization compared with the other genotypes. While interpretation of this finding remains speculative, possible explanations may be the loss of WDFY3-mediated functions required for cellular physiology and dendritic function. As previously reported, WDFY3 supports optimal cellular bioenergetics via mitophagy and glycophagy [62, 63]. An alternative explanation may center around a possible WDFY3 requirement for cytoskeletal organization, as widespread axonal defects in homozygous mutants suggest (unpublished observations and ). Interestingly, decreased dendritic complexity in Wdfy3 mutant neurons agrees with findings in other ASD risk genes [65,66,67,68,69], has been observed postmortem in the hippocampi of autism cases , and may be a point of convergence between ASD and Rett syndrome .
Progressive Wdfy3 gene loss results in concomitant increase in synaptic spine density. This finding agrees with previous reports of added dendritic spine density in layer V pyramidal neurons in the temporal lobe of ASD cases . Apparently, augmented spine density in ASD is the result of reduced developmental spine pruning precipitated by hyperactivated mTOR signaling and resulting autophagy impairment . Weakened autophagic function appears to be the molecular cause of reduced spine elimination in ASD cases, as well as in Tsc1/2- and Atg7-deficient cortical projection neurons, as analysis of mutant mice confirmed, drawing intriguing parallels to WDFY3 and its well-established role in regulating selective macroautophagy. While previous results did not show diminished autophagic flux in brains of Wdfy3 mutant mice, these studies were either conducted in embryonic  or perinatal stages of homozygous mutants  and, due to perinatal lethality of homozygous Wdfy3 mutation, only in heterozygous adults . Thus, subtle changes in autophagy function caused by the synergy of homozygous Wdfy3 loss and progressed age may have eluded detection. Of course, it is possible that greater spine density in Wdfy3 mutant cells is caused by autophagy-independent mechanisms that currently remain elusive.
Wdfy3-MADM produced the predicted sparse labeling with WT red and lacZ/lacZ green cells being generated in the same ratio, while + /lacZ yellow cells were twice as frequent as homozygous genotypes. The rates of labeled cells are consistent with past findings and the expectation that yellow cells being born out of G2-Z events, G1 recombination, and postmitotic recombination (G0) will amount to approximately double as many cells as either red or green cells [37, 38, 41]. The fact that yellow heterozygous cells are born at different timepoints poses difficulties in comparing phenotypic features of this genotype to the other genotypes. However, we noted that at least with respect to spine density heterozygous cells assumed an intermediate phenotype between WT and homozygous mutant cells. This finding indicated dosage sensitivity which may be relevant to human disease pathogenesis that is typically associated with haploinsufficiency. In these situations, additional, conceivably environmental, factors may be responsible for full expression of the autistic phenotype [72,73,74,75].
Whether migration defects are present in typically heterozygous WDFY3 patients remains uncertain as of now. While MRI data for some patients are available with no evidence of focal cortical dysplasia , limitations in resolution cannot provide certainty that small-scale lamination defects do not occur and may influence the clinical presentation in these patients. Our work in mice found no mispositioned neurons in heterozygous mutants, but is also facing the constraint that only two hypomorphic alleles have been examined in this respect, much fewer than the number of known human alleles (49 alleles listed in denovo-db and 13 alleles reported in ). Potential lamination defects have not been analyzed in other available Wdfy3 mutant mouse lines .
Wdfy3 mutant mice have been a valuable tool in understanding prenatal aspects of ASD pathology that remain inaccessible to investigation in humans. However, like most animal models they are not a perfect model in replicating all aspects of human ASD biology, as ASD-relevant behavioral deviations in these mice remain inconclusive as of yet. It is also uncertain whether alterations in neuronal morphology described here apply to ASD patients, considering the greater complexity of the human brain and the lack of patient data.
In humans, WDFY3 is a recognized autism risk gene associated with macrocephaly. Wdfy3 mutations in mice produce a complex phenotype of cortical maldevelopment including focal neuronal migration defects of excitatory neurons. Using MADM technology in this study, we arrive at the conclusion that these migration errors are driven by cell autonomous mechanisms that act directly on mutant neurons and their progenitors rather than impaired cell–environment interactions. In addition, we discovered alterations in mutant neuronal morphology consistent with other reports documenting decreased dendritic arbor complexity but increased spine density in autism mouse models and human cases. Our findings further underline the validity of Wdfy3 mutant mice as ASD models and point to the significance of changes at the chemical synapse in ASD neuropathology.
Availability of data and materials
Datasets analyzed and experimental tools such as transgenic mice used in the current study are available from the corresponding author on reasonable request.
Analysis of variance
Autism spectrum disorders
Green fluorescence protein
Ultrafast optical clearing method
Mosaic analysis with double markers
Tandem dimer Tomato
WD repeat and FYVE domain containing 3
Nadarajah B, Parnavelas JG. Modes of neuronal migration in the developing cerebral cortex. Nat Rev Neurosci. 2002;3(6):423–32.
Ayala R, Shu T, Tsai LH. Trekking across the brain: the journey of neuronal migration. Cell. 2007;128(1):29–43.
Kwan KY, Sestan N, Anton ES. Transcriptional co-regulation of neuronal migration and laminar identity in the neocortex. Development. 2012;139(9):1535–46.
Hippenmeyer S. Molecular pathways controlling the sequential steps of cortical projection neuron migration. Adv Exp Med Biol. 2014;800:1–24.
Gleeson JG, Walsh CA. Neuronal migration disorders: from genetic diseases to developmental mechanisms. Trends Neurosci. 2000;23(8):352–9.
Gupta A, Tsai LH, Wynshaw-Boris A. Life is a journey: a genetic look at neocortical development. Nat Rev Genet. 2002;3(5):342–55.
Bielas S, et al. Cortical neuronal migration mutants suggest separate but intersecting pathways. Annu Rev Cell Dev Biol. 2004;20:593–618.
Guerrini R, Filippi T. Neuronal migration disorders, genetics, and epileptogenesis. J Child Neurol. 2005;20(4):287–99.
Guerrini R, Parrini E. Neuronal migration disorders. Neurobiol Dis. 2010;38(2):154–66.
Kriegstein AR, Noctor SC. Patterns of neuronal migration in the embryonic cortex. Trends Neurosci. 2004;27(7):392–9.
Wu Q, et al. The dynamics of neuronal migration. Adv Exp Med Biol. 2014;800:25–36.
Tassi L, et al. Focal cortical dysplasia: neuropathological subtypes, EEG, neuroimaging and surgical outcome. Brain. 2002;125(Pt 8):1719–32.
Krsek P, et al. Different features of histopathological subtypes of pediatric focal cortical dysplasia. Ann Neurol. 2008;63(6):758–69.
Kabat J, Krol P. Focal cortical dysplasia - review. Pol J Radiol. 2012;77(2):35–43.
Kato M, Dobyns WB. Lissencephaly and the molecular basis of neuronal migration. Hum Mol Genet, 2003. 12 Spec No 1: p. R89–96.
Fry AE, Cushion TD, Pilz DT. The genetics of lissencephaly. Am J Med Genet C Semin Med Genet. 2014;166C(2):198–210.
Danielsson S, et al. Epilepsy in young adults with autism: a prospective population-based follow-up study of 120 individuals diagnosed in childhood. Epilepsia. 2005;46(6):918–23.
Mouridsen SE, Rich B, Isager T. A longitudinal study of epilepsy and other central nervous system diseases in individuals with and without a history of infantile autism. Brain Dev. 2011;33(5):361–6.
Bailey A, et al. A clinicopathological study of autism. Brain. 1998;121(Pt 5):889–905.
Hutsler JJ, Love T, Zhang H. Histological and magnetic resonance imaging assessment of cortical layering and thickness in autism spectrum disorders. Biol Psychiatry. 2007;61(4):449–57.
Wegiel J, et al. The neuropathology of autism: defects of neurogenesis and neuronal migration, and dysplastic changes. Acta Neuropathol. 2010;119(6):755–70.
Casanova MF, et al. Focal cortical dysplasias in autism spectrum disorders. Acta Neuropathol Commun. 2013;1:67.
Stoner R, et al. Patches of disorganization in the neocortex of children with autism. N Engl J Med. 2014;370(13):1209–19.
Orosco LA, et al. Loss of Wdfy3 in mice alters cerebral cortical neurogenesis reflecting aspects of the autism pathology. Nat Commun. 2014;5:4692.
Choi GB, et al. The maternal interleukin-17a pathway in mice promotes autism-like phenotypes in offspring. Science. 2016;351(6276):933–9.
Iossifov I, et al. De novo gene disruptions in children on the autistic spectrum. Neuron. 2012;74(2):285–99.
Iossifov I, et al. The contribution of de novo coding mutations to autism spectrum disorder. Nature. 2014;515(7526):216–21.
Wang T, et al. De novo genic mutations among a Chinese autism spectrum disorder cohort. Nat Commun. 2016;7:13316.
Stessman HA, et al. Targeted sequencing identifies 91 neurodevelopmental-disorder risk genes with autism and developmental-disability biases. Nat Genet. 2017;49(4):515–26.
Yuen RKC, et al. Whole genome sequencing resource identifies 18 new candidate genes for autism spectrum disorder. Nat Neurosci. 2017;20(4):602–11.
Guo H, et al. Inherited and multiple de novo mutations in autism/developmental delay risk genes suggest a multifactorial model. Mol Autism. 2018;9:64.
Le Duc, D., et al., Pathogenic WDFY3 variants cause neurodevelopmental disorders and opposing effects on brain size. Brain, 2019.
Keays DA, et al. Mutations in alpha-tubulin cause abnormal neuronal migration in mice and lissencephaly in humans. Cell. 2007;128(1):45–57.
Gleeson JG, et al. Doublecortin, a brain-specific gene mutated in human X-linked lissencephaly and double cortex syndrome, encodes a putative signaling protein. Cell. 1998;92(1):63–72.
Magen D, et al. Autosomal recessive lissencephaly with cerebellar hypoplasia is associated with a loss-of-function mutation in CDK5. Hum Genet. 2015;134(3):305–14.
Hong SE, et al. Autosomal recessive lissencephaly with cerebellar hypoplasia is associated with human RELN mutations. Nat Genet. 2000;26(1):93–6.
Zong H, et al. Mosaic analysis with double markers in mice. Cell. 2005;121(3):479–92.
Contreras X, et al. A genome-wide library of MADM mice for single-cell genetic mosaic analysis. Cell Rep. 2021;35(12):109274.
Percie-du-Sert N, et al. The ARRIVE guidelines 2.0: Updated guidelines for reporting animal research. J Cereb Blood Flow Metab. 2020;40(9):1769–77.
Council, N.R., Committee for the Update of the Guide for the Care and Use of Laboratory animals. Guide for the Care and Use of Laboratory Animals. Guide for the Care and Use of Laboratory Animals. 2011, Washington (DC): National Academies Press (US) National Academy of Sciences.
Amberg N, Hippenmeyer S. Genetic mosaic dissection of candidate genes in mice using mosaic analysis with double markers. STAR Protoc. 2021;2(4):100939.
Zhu X, et al. Ultrafast optical clearing method for three-dimensional imaging with cellular resolution. Proc Natl Acad Sci USA. 2019;116(23):11480–9.
Verney C, et al. Independent controls for neocortical neuron production and histogenetic cell death. Dev Neurosci. 2000;22(1–2):125–38.
Zoghbi HY. Postnatal neurodevelopmental disorders: meeting at the synapse? Science. 2003;302(5646):826–30.
Martinez-Cerdeno V. Dendrite and spine modifications in autism and related neurodevelopmental disorders in patients and animal models. Dev Neurobiol. 2017;77(4):393–404.
Risher WC, et al. Rapid Golgi analysis method for efficient and unbiased classification of dendritic spines. PLoS ONE. 2014;9(9):e107591.
Evsyukova I, Plestant C, Anton ES. Integrative mechanisms of oriented neuronal migration in the developing brain. Annu Rev Cell Dev Biol. 2013;29:299–353.
Hansen AH, Hippenmeyer S. Non-cell-autonomous mechanisms in radial projection neuron migration in the developing cerebral cortex. Front Cell Dev Biol. 2020;8:574382.
Dobyns WB, et al. Lissencephaly. A human brain malformation associated with deletion of the LIS1 gene located at chromosome 17p13. JAMA. 1993;270(23):2838–42.
Pilz DT, et al. LIS1 and XLIS (DCX) mutations cause most classical lissencephaly, but different patterns of malformation. Hum Mol Genet. 1998;7(13):2029–37.
Faulkner NE, et al. A role for the lissencephaly gene LIS1 in mitosis and cytoplasmic dynein function. Nat Cell Biol. 2000;2(11):784–91.
Smith DS, et al. Regulation of cytoplasmic dynein behaviour and microtubule organization by mammalian Lis1. Nat Cell Biol. 2000;2(11):767–75.
Tsai JW, et al. LIS1 RNA interference blocks neural stem cell division, morphogenesis, and motility at multiple stages. J Cell Biol. 2005;170(6):935–45.
D’Arcangelo G, et al. A protein related to extracellular matrix proteins deleted in the mouse mutant reeler. Nature. 1995;374(6524):719–23.
Hirotsune S, et al. The reeler gene encodes a protein with an EGF-like motif expressed by pioneer neurons. Nat Genet. 1995;10(1):77–83.
Ogawa M, et al. The reeler gene-associated antigen on Cajal-Retzius neurons is a crucial molecule for laminar organization of cortical neurons. Neuron. 1995;14(5):899–912.
Hashimoto-Torii K, et al. Interaction between Reelin and Notch signaling regulates neuronal migration in the cerebral cortex. Neuron. 2008;60(2):273–84.
Franco SJ, et al. Reelin regulates cadherin function via Dab1/Rap1 to control neuronal migration and lamination in the neocortex. Neuron. 2011;69(3):482–97.
Nikolic M, Gardner HA, Tucker KL. Postnatal neuronal apoptosis in the cerebral cortex: physiological and pathophysiological mechanisms. Neuroscience. 2013;254:369–78.
Wong FK, et al. Pyramidal cell regulation of interneuron survival sculpts cortical networks. Nature. 2018;557(7707):668–73.
Blanquie O, et al. Electrical activity controls area-specific expression of neuronal apoptosis in the mouse developing cerebral cortex. Elife. 2017;6.
Napoli E, et al. Beyond autophagy: a novel role for autism-linked Wdfy3 in brain mitophagy. Sci Rep. 2018;8(1):11348.
Napoli E, et al. Wdfy3 regulates glycophagy, mitophagy, and synaptic plasticity. J Cereb Blood Flow Metab. 2021;41(12):3213–31.
Dragich JM, et al. Autophagy linked FYVE (Alfy/WDFY3) is required for establishing neuronal connectivity in the mammalian brain. Elife. 2016. 5.
Berman RF, et al. Abnormal dendrite and spine morphology in primary visual cortex in the CGG knock-in mouse model of the fragile X premutation. Epilepsia. 2012;53(Suppl 1):150–60.
de Anda FC, et al. Autism spectrum disorder susceptibility gene TAOK2 affects basal dendrite formation in the neocortex. Nat Neurosci. 2012;15(7):1022–31.
Anderson GR, et al. Candidate autism gene screen identifies critical role for cell-adhesion molecule CASPR2 in dendritic arborization and spine development. Proc Natl Acad Sci USA. 2012;109(44):18120–5.
Khatri N, et al. The autism protein Ube3A/E6AP remodels neuronal dendritic arborization via caspase-dependent microtubule destabilization. J Neurosci. 2018;38(2):363–78.
Bahry JA, et al. An autism-associated de novo mutation in GluN2B destabilizes growing dendrites by promoting retraction and pruning. Front Cell Neurosci. 2021;15:692232.
Raymond GV, Bauman ML, Kemper TL. Hippocampus in autism: a Golgi analysis. Acta Neuropathol. 1996;91(1):117–9.
Tang G, et al. Loss of mTOR-dependent macroautophagy causes autistic-like synaptic pruning deficits. Neuron. 2014;83(5):1131–43.
Chaste P, Leboyer M. Autism risk factors: genes, environment, and gene-environment interactions. Dialogues Clin Neurosci. 2012;14(3):281–92.
Stamou M, et al. Neuronal connectivity as a convergent target of gene x environment interactions that confer risk for Autism Spectrum Disorders. Neurotoxicol Teratol. 2013;36:3–16.
Tordjman S, et al. Gene x Environment interactions in autism spectrum disorders: role of epigenetic mechanisms. Front Psychiatry. 2014;5:53.
Courchesne E, et al. The ASD Living Biology: from cell proliferation to clinical phenotype. Mol Psychiatry. 2019;24(1):88–107.
We would like to thank Angelo Harlan de Crescenzo for early contributions to this project.
This study was funded by NIMH R21MH115347 to KSZ. KSZ is further supported by Shriners Hospitals for Children.
Ethics approval and consent to participate
Animals were handled in accordance with protocols approved by the University of California at Davis Institutional Animal Care and Use Committee overseen by the AAALAC International accreditation program (latest accreditation in February 14, 2020) and in compliance with the ARRIVE guidelines (Animal Research: Reporting in Vivo Experiments).
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The original version of this article was revised: the omitted 4th author (Alexios A. Panoutsopoulos) was added to the author group.
About this article
Cite this article
Schaaf, Z.A., Tat, L., Cannizzaro, N. et al. WDFY3 mutation alters laminar position and morphology of cortical neurons. Molecular Autism 13, 27 (2022). https://doi.org/10.1186/s13229-022-00508-3
- Cerebral cortex
- Neuronal migration
- Excitatory neurons
- Dendritic spines