The coding region of FOXG1’s N-terminus is a mutation hotspot
The human FOXG1 gene possesses three prominently high GC-rich regions in the coding sequences for the N-terminus of the FOXG1 protein, seven C (7C, c.250–256), seven G (7G, c.454–460), and six G (6G, c.501–506) (Fig. 1a). Correspondingly, FS-causing frameshift variants are over-represented in the N-terminal region of the FOXG1 protein prior to the DNA-binding Forkhead domain. In the FOXG1 Syndrome Research Registry, which enrolled 256 individuals with FS, frameshift variants are the most frequent, accounting for 40.6% of participants (104 out of 256 cases) (Fig. 1b). 81% of frameshift variants (84 out of 104 cases) occur before the DNA-binding domain (Fig. 1c). 64.4% of frameshift cases (67 out of 104 cases) were found in the 7C, 7G, and 6G regions, indicating that these C/G-rich regions are prone to the mutations that cause FS.
We reviewed the clinical data of eight individuals with FOXG1 variants at c.250–256 (7C) or c.454–460 (7G) positions (Fig. 1d,e). Over 50% of individuals with variants in these regions showed microcephaly, strabismus, movement disorder, hypotonia, epilepsy, sleep disturbance, and abnormalities of the corpus callosum agenesis (Fig. 1e). Five of eight patients also presented with delayed myelination detected via magnetic resonance imaging of the brain (Fig. 1d).
Generation of mice carrying Q84Pfs allele
To understand the pathophysiology of FS, we generated FS mouse model bearing a specific FS-causing FOXG1 variant. The knock-in mice carry c.250dupC (p.Q84Pfs*31, herein referred to as Q84Pfs), equivalent to human c.256dupC (p.Q86Pfs*35, referred to as Q86Pfs), representing one of the most prevalent variants in FS (Fig. 1a, 2a,b).
Foxg1Q84Pfs/Q84Pfs homozygous (Q84Pfs-Homo) mice exhibited a marked loss of forebrain tissue and died soon after birth (Fig. 2c, Supplementary Fig. 1). They also showed craniofacial defects, such as a short frontonasal process, and underdeveloped eye morphology (Supplementary Fig. 1). In contrast, Foxg1Q84Pfs/+ heterozygous (Q84Pfs-Het) mice survived into adulthood.
Structural brain deficits in Q84Pfs-Het mice
Compared to wild-type (WT) littermate controls, adult Q84Pfs-Het mice showed smaller brains, shortened corpus callosum, and profound septum defects (Fig. 2d,e). The hippocampus size was reduced, and the malformation of the dentate gyrus was notable (Fig. 2d). Thus, Q84Pfs-Het mice recapitulate the brain structural deficits in human FS with N-terminal frameshift variants (Fig. 1d,e).
RNA-seq analysis identified the dysregulated gene sets in Q84Pfs-Het cortex
To identify the dysregulated genes contributing to the deficits of Q84Pfs-Het brains, we performed RNA-seq in P1 cortices of Q84Pfs-Het and littermate WT control mice and compared their gene expression patterns. 222 genes were significantly changed in the Q84Pfs-Het cortex relative to control WT cortex, representing differentially expressed genes (DEGs) (Fig. 3a, Supplementary Dataset). Interestingly, genes expressed in the cortical upper-layer projection neurons were downregulated, whereas deep-layer projection neuronal genes were upregulated (Fig. 3b), suggesting that the microcephaly of Q84Pfs-Het cortex is at least partly attributed to the reduction of the upper-layer, not deep-layer. Many genes involved in GABAergic interneuron development were also downregulated in Q84Pfs-Het cortex (Fig. 3b), indicating defects in the development of cortical interneurons.
Approximately 30% of DEGs in Q84Pfs-Het cortex were associated with Foxg1-bound ChIP-seq peaks in the developing cortex 11 regardless of up- or down-regulation (Fig. 3c), suggesting that Foxg1 regulates its direct target genes both negatively and positively. Consistently, the HOMER motif analysis27 revealed that the top two significantly enriched motifs in Foxg1-bound peaks annotated to either up- or down-regulated genes are the FOX motif and E-box. E-box serves as the binding site for basic helix-loop-helix (bHLH) and Zbtb18 (Rp58) transcription factors, both of which were shown to collaborate with Foxg1 11,28 (Fig. 3d). The NF1 site was also significantly associated with both up- and down-regulated Foxg1-target genes (Fig. 3d). These data suggest that a subset of dysregulated genes in Q84Pfs-Het mice are direct target genes of Foxg1 and that Foxg1 cooperates with ZBTB18, bHLH transcription factor, and NF1 in the developing cortex.
The dysregulated pathways in Q84Pfs-Het cortex
To gain insights into molecular and cellular pathways leading to defects of Q84Pfs-Het mice, we performed the gene set enrichment analysis (GSEA) of DEGs using the Database for Annotation, Visualization and Integrated Discovery (DAVID) 29 and Enrichr 30.
In tissue analysis, DEGs were enriched in neuronal epitheliums, prefrontal cortex, cerebral cortex, cingulate gyrus, and superior frontal gyrus (Fig. 3e), consistent with the use of cortex for RNA-seq. In cell type analysis, the upregulated genes were significantly associated with cortical layer 6 glutamatergic neurons and, interestingly, oligodendrocytes, in which Foxg1 function remains unclear (Fig. 3f). The downregulated genes were most highly enriched for the “interneuron: embryonic prefrontal cortex” (Fig. 3f). Thus, our bioinformatic analyses indicate that the three major cell types in the neonatal cortex, excitatory projection neurons, inhibitory GABAergic interneurons, and the oligodendrocyte lineage, have a deficiency in Q84Pfs-Het mice.
Consistent with the known function of Foxg1 10–14, DEGs of Q84Pfs-Het cortex were enriched for the genes controlling neuronal projection development, neuronal cell body, axon, and transcription regulations (Fig. 3g,h). Interestingly, the gene involved in the negative regulation of cell proliferation was downregulated in Q84Pfs-Het cortex (Fig. 3h), suggesting that progenitor proliferation may increase in Q84Pfs-Het cortex unlike the heterozygous mice carrying Foxg1-null allele (global Foxg1-Het)2,5,8,9.
Intriguingly, many synaptic genes were significantly dysregulated in Q84Pfs-Het cortex. The upregulated genes were enriched for biological process (BP) terms of synaptic vesicle transport, synaptic vesicle exocytosis and endocytosis, and presynaptic membrane assembly, and cellular component (CC) terms of synapse and presynapse (Fig. 3g,h). Further, the downregulated genes were also strongly associated with the CC terms of synapse and postsynaptic density (Fig. 3g,h). Our results underline the new role of Foxg1 in synapse development. It is also noteworthy that the upregulated and downregulated genes were more strongly associated with presynapse and postsynaptic density, respectively.
Remarkably, both up- and down-regulated genes were enriched for translation and ribosome (Fig. 3g,h), suggesting defects in translation in Q84Pfs-Het brains. Other enriched terms, such as locomotory behavior and memory, are related to FS symptoms, such as intellectual disability and hyperkinetic-dyskinetic movements 19–23.
Together, our results highlight the new role of Foxg1 in the oligodendrocyte development and synapse formation and function.
Developmental defects of neural progenitors and excitatory and inhibitory neurons in Q84Pfs-Het cortex
We next assessed if the dysregulated genes in Q84Pfs-Het brains led to cortical defects suggested by our GO analyses. Interestingly, although Q84Pfs-Het cortex was thinner than the WT cortex by E16 (Fig. 4a,b), Pax6+ neural progenitor cells (NPCs) and Tbr2+ intermediate progenitor cells (IPCs) increased in Q84Pfs-Het cortex at E16 and P0 (Fig. 4a-f). Consistently, phospho-histone H3+ proliferating cells in the ventricular zone were also significantly increased in Q84Pfs-Het cortex (Fig. 4a-f). Combined with the downregulation of genes for negative regulation of cell proliferation (Fig. 3h), these results indicate that progenitor proliferation was augmented in Q84Pfs-Het mice in contrast to the reduced progenitor proliferation in global Foxg1-Het mice 2,5,8,9.
The production of Cux1+ upper-layer neurons was delayed in E16 Q84Pfs-Het cortex, as shown by a lack of Cux1+ neurons in the superficial area of the cortex (Fig. 5a). Cux1+ upper-layer neurons remained significantly reduced in P30 Q84Pfs-Het cortex (Fig. 5b,e), indicating that the delayed upper-layer neuron generation is not compensated at the later stages. The deep-layer neurons showed a trend of increase in number, but the thickness of the deep-layer did not significantly increase (Fig. 5a,b,d). Combined with reduced cortex thickness (Fig. 5. a-c), these results show that mainly the reduction of the upper-layer accounts for the thinner cortex.
Dlx1+ GABAergic interneurons were substantially decreased in Q84Pfs-Het cortex (Fig. 5f,g), consistent with downregulation of genes involved in cerebral cortex tangential migration (i.e., GABAergic interneuron migration deficits) (Fig. 3f,h).
Q84Pfs-Het mice displayed two prominent axon projection deficits in the cortex, matching the dysregulation of neuron projection development genes (Fig. 3g,h). First, a substantial fraction of callosal axons was stalled at the midline, forming the Probst bundle (Fig. 5h). Second, Q84Pfs-Het showed misprojected L1+ axon bundles crossing the cortical plate toward the superficial area. These misrouted L1+ axons expressed Ntng1 (netrin G1) (Fig. 5i), indicating the defects of thalamocortical axon guidance in Q84Pfs-Het cortex.
Together, Q84Pfs-Het cortex exhibited deficits in progenitor proliferation, axon projection, layer formation, and interneuron migration.
Myelination defects in Q84Pfs-Het cortex
Although delayed myelination is one of the most common features of FS brains (Fig. 1d), the role of Foxg1 in the oligodendrocyte lineage development remains unclear. To test if the FS mouse model recapitulates this prominent white matter phenotype in humans, we monitored oligodendrocyte differentiation and myelination in Q84Pfs-Het mice at P30, the active myelination phase. Interestingly, Olig2+ oligodendrocyte lineage cells, many of which include oligodendrocyte progenitor cells (OPCs), significantly increased in Q84Pfs-Het cortex (Fig. 6a,b), consistent with the association of upregulated genes with oligodendrocyte (Fig. 3f). Despite of the increased Olig2+ cells, Mbp+ myelinated area was substantially reduced in Q84Pfs-Het brains (Fig. 6c,d). Moreover, Q84Pfs-Het brains exhibited a marked reduction of the arborization and complexity of myelination patterns (Fig. 6c,e). Our results indicate that the oligodendrocyte development process was disrupted, and the myelination steps, including the microstructural organization of myelination, was impaired in Q84Pfs-Het brains, providing the pathological mechanisms underlying white matter deficiency in FS.
Dysregulated genes in Q84Pfs-Het cortex were associated with motor dysfunction and autistic-like behaviors
To ask if the dysregulated genes in the Q84Pfs-Het cortex are connected to specific human conditions, we performed GSEA with the DisGeNET database, which contains collections of genes and variants associated with human diseases 31 (Fig. 7a). The upregulated genes were significantly associated with impaired social interactions, Alzheimer’s disease, ASD, Parkinson’s disease, and autonomic nervous system disorders (Fig. 7a). The downregulated genes were enriched for dystonia, myoclonic encephalopathy, myoclonus, and ataxia, as well as major depressive disorder. Together, Q84Pfs-Het transcriptome changes strongly linked to impairments in movement, social interactions, and autonomic nervous system function in Q84Pfs-Het mice.
Given that both FS and Huntington’s disease (HD) exhibit many neurological symptoms, including abnormal involuntary movements 22,32, we compared the dysregulated genes between Q84Pfs mice and HD using the database of HD molecular signatures (HDSigDB, hdinhd.org) (Fig. 7b). This analysis revealed the striking resemblance between Q84Pfs-Het mice and HD mouse models 33,34. The high-ranking categories that resemble the up- and down-regulated genes in Q84Pfs-Het mice were the up- and down-regulated genes in the hippocampus of the HD mouse model Q175 33, respectively (Fig. 7b). Notably, the two categories among the top five gene sets for the upregulated genes were oligodendrocyte progenitor cell (OPC) genes (Fig. 7b), consistent with our finding that Q84Pfs-Het cortex showed the increased Olig2+ oligodendrocyte lineage cells (Fig. 6a,b).
Collectively, the transcriptome changes in Q84Pfs-Het cortex are closely connected to the clinical features of FS, such as movement disorders, autism-like behaviors, and social impairment.
Q84Pfs-Het mice showed movement deficits, repetitive behaviors, and prolonged behavior arrest
To investigate whether Q84Pfs-Het mice exhibit behavior phenotypes corresponding to the above molecular and cellular changes, we performed behavioral assessment tests at adolescence stage P30 and adult stages P60 and P90.
There was no significant difference in body weights of Q84Pfs-Het and WT mice at these stages (Fig. 7c). In the wire hang test that evaluates the motor function and muscle strength, Q84Pfs-Het mice showed significantly reduced hang time, which worsened as mice aged (Fig. 7d). In the open field test that monitors locomotor activity and exploratory behaviors, Q84Pfs-Het mice moved significantly less than WT mice at P60 and P90 (Fig. 7e). At P30, Q84Pfs-Het mice showed a tendency for travel distance reduction over time compared to WT mice, but the total travel distance did not differ significantly (Fig. 7e). These results indicate deficits in muscle strength and movement in Q84Pfs-Het mice.
Q84Pfs-Het mice exhibited a strikingly extended behavior arrest, defined by paused locomotion for longer than 3 min at one episode during free exploration of the arena (Fig. 7f, Supplementary Fig. 2 and videos). Throughout these episodes of remarkably pronounced stops, Q84Pfs-Het mice showed a high degree of postural control, typically positioning their bodies facing the center of the arena, and no apparent overt movements, such as visual survey, rearing, or digging. The proportion of Q84Pfs-Het mice that exhibited the extended behavior arrest was 11%, 22%, and 47% at P30, P60, and P90, respectively, whereas WT mice did not show behavior arrest without purposeful movements under the same condition (Fig. 7f). Notably, Q84Pfs-Het mice without prolonged behavior arrests still showed a decreased travel distance, suggesting that the extended arrest is not the only contributing factor to movement reduction.
Next, we assessed ASD-related behaviors. The center time in the open field test was markedly decreased as Q84Pfs-Het mice aged from P30 to P90, relative to WT mice (Fig. 7g), indicating increased anxiety levels. Q84Pfs-Het mice exhibited a significant increase in repetitive grooming behavior at all three stages, with a tendency for bigger differences from WT mice as they age (Fig. 7h). These results suggest heightened repetitive behavior in Q84Pfs-Het mice. Intriguingly, in the marble-burying test, Q84Pfs-Het mice buried more marbles at P30 but fewer marbles at the adult stages than WT mice at P60 and P90 (Fig. 7i), suggesting age-dependent changes in the marble burying behavior.
Together, our data uncovered that Q84Pfs-Het mice have motor coordination defects, increased anxiety levels, remarkably prolonged behavior arrest, and repetitive behaviors, all observed in FS 19–23 and strongly linked to molecular and cellular changes of Q84Pfs-Het brains.