Chd8 mutation in oligodendrocytes alters microstructure and functional connectivity in the mouse brain
Molecular Brain volume 13, Article number: 160 (2020)
CHD8 encodes a chromatin-remodeling factor and is one of the most recurrently mutated genes in individuals with autism spectrum disorder (ASD). Although we have recently shown that mice heterozygous for Chd8 mutation manifest myelination defects and ASD-like behaviors, the detailed mechanisms underlying ASD pathogenesis have remained unclear. Here we performed diffusion tensor imaging (DTI) and resting-state functional magnetic resonance imaging (rsfMRI) in oligodendrocyte lineage-specific Chd8 heterozygous mutant mice. DTI revealed that ablation of Chd8 specifically in oligodendrocytes of mice was associated with microstructural changes of specific brain regions including the cortex and striatum. The extent of these changes in white matter including the corpus callosum and fornix was correlated with total contact time in the reciprocal social interaction test. Analysis with rsfMRI revealed changes in functional brain connectivity in the mutant mice, and the extent of such changes in the cortex, hippocampus, and amygdala was also correlated with the change in social interaction. Our results thus suggest that changes in brain microstructure and functional connectivity induced by oligodendrocyte dysfunction might underlie altered social interaction in mice with oligodendrocyte-specific CHD8 haploinsufficiency.
Autism spectrum disorder (ASD) encompasses a range of neurodevelopmental disorders characterized by deficits in social interaction and communication as well as by restricted and repetitive behaviors. Structural and functional alterations in the brain of individuals with ASD have been identified by magnetic resonance imaging (MRI). For example, volumetric differences in several brain regions including the cortex, striatum, amygdala, and cerebellum have been detected in such individuals [1, 2]. Diffusion tensor imaging (DTI), a specific type of MRI, has also revealed microstructural changes in gray and white matter that are associated with ASD . Moreover, resting-state functional MRI (rsfMRI) has uncovered a complex pattern of both hypo- and hyperconnectivity across the brain of ASD patients [4, 5]. Whether the alterations detected by such imaging studies contribute to the symptoms of ASD has remained to be demonstrated, however.
Mutations in the gene for chromodomain helicase DNA-binding protein 8 (CHD8) have been identified as a highly penetrant risk factor for ASD by exome-sequencing studies [6,7,8,9]. CHD8 is a member of the CHD family of enzymes that belong to the SNF2 superfamily of ATP-dependent chromatin remodelers, and it has been shown to regulate the expression of developmental and ASD risk genes including those related to synaptic function and neurogenesis [10,11,12,13]. In humans, CHD8 mutations give rise to core ASD phenotypes as well as macrocephaly and gastrointestinal complaints . MRI of individuals with CHD8 mutations has revealed white matter defects, ventriculomegaly, and unspecified abnormalities [14,15,16]. Furthermore, we and others have found that heterozygous deletion of Chd8 in mice results in altered social behavior, increased anxiety-like behavior, and cognitive deficits reminiscent of those seen in humans with CHD8 mutations [17,18,19,20,21]. Chd8 mutant mice are thus a powerful tool for investigation of the pathogenesis of ASD.
MRI is a key method for translational research to bridge the gap between patients and animal models. In addition to MRI evaluation of ASD patients, recent studies have identified structural and functional alterations across the brain in several mouse models of ASD. One such mouse model, the BTBR T + tf/J mouse, manifests volumetric and structural changes in the striatum, nucleus accumbens, thalamus, and amygdala that were found to be correlated with changes in sociability and repetitive behaviors [22, 23]. Several genetic mouse models of ASD also manifest altered functional connectivity among various brain regions [24,25,26,27,28,29]. MRI studies in Chd8 heterozygous mutant mice have thus suggested the presence of microstructural changes in several brain regions as well as functional overconnectivity in the cortico-hippocampal network [19, 20]. Although we recently found that CHD8 haploinsufficiency in oligodendrocytes results in myelination defects and delayed action potential propagation as well as recapitulates certain behavioral phenotypes of Chd8 heterozygous mutant mice , whether oligodendrocyte dysfunction induced by Chd8 mutation gives rise to such structural and functional changes in the brain has remained unclear.
We have now searched for structural and functional changes in the brain of oligodendrocyte lineage-specific Chd8 heterozygous mutant mice with the use of MRI. DTI revealed an altered microstructure in some regions of the cortex and striatum in the mutant mice. Furthermore, rsfMRI showed that oligodendrocyte lineage-specific ablation of Chd8 affected functional connectivity between several brain regions related to social behavior, and that these changes were correlated with the deficit apparent in the reciprocal social interaction test. Our results thus uncover an apparent relation among brain structural and functional alterations as well as a behavioral change induced by oligodendrocyte dysfunction due to Chd8 mutation.
Materials and methods
Generation of Chd8LF/F mice was described previously . A pair of loxP sites flanking exons 11 to 13 of Chd8 in Chd8L+/F mice was deleted by Cre recombinase. Olig1-Cre transgenic mice were described previously . Chd8L+/F mice were crossed with Olig1-Cre heterozygous mice to produce Olig1-Cre/Chd8L+/F mice. Offspring were backcrossed onto the C57BL/6J line for at least nine generations. All experiments were performed with male mice. Mice were genotyped by polymerase chain reaction-based analysis of genomic DNA with primers for Chd8 (5′-CCCAAAAGACCAAATCAAACAAAC-3′, 5′-CCATAGGCTGAAGAACCGTAATTG-3′, and 5′-AGGCTTAGAAACCCGTCGAG-3′) and Cre (5′-AGGTTCGTTCACTCATGGA-3′ and 5′-TCGACCAGTTTAGTTACCC-3′). All animals were maintained under specific pathogen-free conditions.
Reciprocal social interaction test
Two mice at 8 to 10 weeks of age and of identical genotypes that had previously been housed in different cages were placed together in a box (40 by 40 by 30 cm) and allowed to explore freely for 10 min (n = 12 mice per genotype). Images were captured at a rate of three frames per second. Analysis was performed automatically with the use of ImageSI software, which was developed by Tsuyoshi Miyakawa at Fujita Health University and is based on ImageJ [30, 32].
rsfMRI analysis in awake mice
rsfMRI was performed on awake mice (n = 12 mice per genotype) as previously described [33,34,35,36]. Male mice at 6 weeks of age were anesthetized with 2% isoflurane. A custom-made acrylic head bar (3 by 3 by 27 mm) was mounted along the sagittal suture of the exposed skull with the use of dental cement (Super-Bond C&B; Sun Medical, Shiga, Japan). The exposed skull was covered with the dental acrylic, and the animal was returned to its cage and allowed to recover. Mice were acclimated to a mock rsfMRI environment for 2 h per day for at least 10 days before being subjected to rsfMRI in the awake state. Under such conditions, body weight, fecal weight, the electromyogram, and the electrocardiogram of the mice resembled those of unrestrained control mice by the 6th day of habituation . The reciprocal social interaction test was performed before rsfMRI recording with a 7.0-T MRI apparatus equipped with actively shielded gradients at a maximum strength of 700 mT/m (Biospec 70/16; Bruker, Ettlingen, Germany), a cryogenic quadrature radio frequency surface probe (CryoProbe Z120046, Bruker), and ParaVision 6.01 software interface (Bruker). Animal respiration was monitored during the performance of rsfMRI. There was a gap of ~ 1 mm between the top of the head bar holder and the CryoProbe . B0 homogeneity was ensured with the use of automatic iterative FASTMAP methods (ParaVision 6.01). For construction of a reference for brain anatomy, high-resolution T2-weighted images of the whole brain were acquired from mice under anesthesia with 2% isoflurane with the use of a Rapid Acquisition with Relaxation Enhancement (RARE) method and with the following parameters: effective time to echo (TE), 27 ms; time to repetition (TR), 2500 ms; number of averages, 10; RARE factor, 8; matrix size, 182 by 182; field of view (FOV), 12.8 by 12.8; spatial resolution, 0.07 by 0.07 mm; slice thickness, 0.45 mm; number of slices, 30; and slice gap, 0.05 mm. The blood oxygen level-dependent (BOLD) signal was acquired with the use of a gradient echo–echo planar imaging method with the following parameters: TE, 15 ms; TR, 2000 ms; number of averages, 1; segment, 1; band width, 250 kHz; matrix size, 64 by 64; FOV, 12.8 by 12.8; spatial resolution, 0.2 by 0.2 mm; slice thickness, 0.7 mm; number of slices, 20; slice gap, 0.05 mm; flip angle, 70°; and dummy scan, 0. The scan was repeated 300 times in 10 min. Preprocessing and statistical analysis were performed with SPM12 (Welcome Trust Centre for Neuroimaging, London, UK) and in-house software written in MATLAB. Image preprocessing was first performed individually for each animal. The reference images obtained by RARE were segmented with tissue probability maps of gray matter, white matter, and cerebrospinal fluid . Time-series images of BOLD-rsfMRI were realigned to correct for residual head motion and slice timing, and were coregistered to the reference structural images. All images were then spatially normalized and coregistered to the standard structural brain . Finally, the images were resliced to a resolution of 0.2 by 0.2 by 0.2 mm and smoothed with a full width at half maximum (FWHM) Gaussian kernel of 0.6 mm. Functional connectivity was calculated with time-series images of BOLD signals with functional connectivity toolbox (CONN, available at https://www.nitrc.org/projects/conn) and in-house software written in MATLAB. The ART-based scrubbing method in CONN (movement threshold, 2; rotation threshold, 1.15) was used to remove motion artifacts. Signals from white matter and cerebrospinal fluid were also regressed out. Data were detrended and temporally filtered at 0.01 to 0.1 Hz. The global mean signal was regressed out from the preprocessed time-series data in order to reduce global noise and nonneural signal correlations. Analysis of rsfMRI data without global signal regression (GSR) was also performed. Seed-based functional connectivities were computed with the use of regions of interest (ROIs) of 18 bilateral anatomic regions of the cortex divided into left and right hemispheres (MO, SS, AUD, VIS, ACA, PL, ILA, ORB, AI) or 11 bilateral brain loci related to social behavior (PL, ILA, ORB, DG, CA1, SUB, BLA, CP, ACB, VTA, PAG) that were defined by our flexible annotation atlas of the mouse brain based on the Allen Mouse Brain Atlas . Pearson’s correlation coefficients for the time courses of two ROIs were calculated after denoising and despiking and were transformed into Fisher’s z-score. Regions of white matter and cerebrospinal fluid were excluded. An 18 by 18 or 11 by 11 connectivity matrix was then generated for the two groups.
To provide some indication of the overall signal level in our imaging data, we calculated the temporal signal-to-noise ratio (tSNR) of BOLD signal fluctuations as the mean voxel value across 194 time points divided by the temporal standard deviation of the voxel signal . The averaged tSNR value was 70.8 ± 6.6 (n = 24), which is similar to that for a previous study , suggesting that the loss of tSNR is unlikely in our study.
DTI acquisition and analysis
The method for DTI acquisition was performed as previously described but with some modifications for the in vivo DTI scans [40, 41]. A DTI-spin echo protocol was used for the acquisition of DTI images of the whole brain of mice (n = 12 mice per genotype) under anesthesia with 2% isoflurane with the following parameters: TR, 2800 ms; TE, 17.5 ms; number averaged, 1; segment, 8; band width, 250 kHz; matrix size, 75 by 75; FOV, 13.5 by 13.5 mm; spatial resolution, 0.18 by 0.18 mm; slice thickness, 0.65 mm; number of slices, 21; slice gap, 0.05 mm; b-value, 1000 s/mm2; image number with b = 0 s/mm2, 6; number of encoding directions of the motion probing gradient, 30; diffusion time δ, 3 ms; Δ, 9 ms; acquisition time, 806 s. The method for DTI analysis was as described previously [40, 41]. The analysis was performed with DSI studio (http://dsi-studio.labsolver.org). A diffusion tensor was modeled at each voxel of the DTI image. Scalar anisotropy and diffusivity maps were obtained from the resulting diffusion tensor eigenvalues (λ1, λ2, λ3), which captured the size of the longest, middle, and shortest axes of the ellipsoid. Fractional anisotropy (FA), radial diffusivity (RD), mean diffusivity (MD), and axial diffusivity (AD) were calculated from standard formulae [42, 43]. An ROI was set with the use of MarsBar (MRC Cognition and Brain Sciences Unit, Cambridge, UK) for significantly altered clusters (voxel-level P < 0.01, with a cluster-level false discovery rate [FDR] correction of P < 0.05) including MO, SS, acc and pcc, ac, aCP and pCP, CP, ACA, CLA, fx, PALd, fi, and MRN. The same ROI was used to measure FA, RD, MD, and AD values with MarsBar.
General protocol for behavioral tests and MRI analysis
Chd8 mutant (Olig1-Cre/Chd8L+/F) or control (Olig1-Cre/Chd8+/+) mice were subjected to the reciprocal social interaction test and MRI analysis on the same day. The procedures were performed in the following sequence: social interaction test, awake rsfMRI analysis, DTI acquisition, and T2-weighted imaging.
All bar graphs show the mean + SEM. Statistical analysis was performed with MATLAB (MathWorks), JMP (SAS Institute, Cary, NC, USA), and Excel (Microsoft, Redmond, WA, USA) software. The Benjamini–Hochberg FDR correction was applied to correct for multiple comparisons within the results of each analysis. Correlations among brain microstructure, functional connectivity, and behavior were evaluated with Pearson’s correlation coefficient and with correction for multiple comparisons. Uncorrected P values are shown together with asterisks (*P < 0.05, **P < 0.01, ***P < 0.001) to indicate the significance level after correction by the FDR.
Oligodendrocyte lineage-specific Chd8 heterozygous mutant mice manifest microstructural alterations in the brain
We performed DTI analysis to examine microstructural changes in axons and the myelin sheath throughout the brain in oligodendrocyte lineage–specific Chd8 heterozygous mutant (Olig1-Cre/Chd8L+/F) and control (Olig1-Cre/Chd8+/+) mice. The significantly altered regions of FA, MD, RD, and AD in Olig1-Cre/Chd8L+/F mice relative to control mice were represented in coronal maps (voxel-level P < 0.01, with a cluster-level FDR correction of P < 0.05) (Fig. 1a and Additional file 1: Table S1). Compared with control mice, Olig1-Cre/Chd8L+/F mice showed a decrease in FA and an increase in MD and RD in the motor cortex (MO), somatosensory cortex (SS), anterior cingulate cortex (ACA), and anterior commissure (ac). The FA values in MO, SS, ACA, and ac of Olig1-Cre/Chd8L+/F mice were significantly lower than those of control mice (Fig. 1b). The MD value in ac, but not those in MO, SS, and ACA, was significantly higher in Olig1-Cre/Chd8L+/F mice than in control animals (Fig. 1c). The RD values in MO, SS, ACA, and ac were significantly higher in Olig1-Cre/Chd8L+/F mice than in control mice (Fig. 1d). Furthermore, the AD value in ac was significantly higher in the mutant mice compared with control mice, whereas those in MO, SS, and ACA did not differ significantly between the two genotypes (Fig. 1e). These results suggested that the microstructure of some regions of the cortex and ac were altered in Olig1-Cre/Chd8L+/F mice.
We also evaluated social behavior in Olig1-Cre/Chd8L+/F and control mice before MRI analysis (Additional file 3: Fig. S1), thereby allowing comparison of performance in the reciprocal social interaction test with MRI measurements in the same animals. We previously found that Olig1-Cre/Chd8L+/F mice manifest altered social interaction including an increased contact time during the reciprocal social interaction test . The brain regions showing a significant correlation between total contact time in the behavioral test and DTI indices (FA, MD, RD, and AD) were represented in coronal maps (voxel-level P < 0.01, with a cluster-level FDR correction of P < 0.05) (Fig. 2a and Additional file 1: Table S1). The FA values in the anterior and posterior corpus callosum (acc and pcc), anterior and posterior caudoputamen (aCP and pCP), fornix (fx), dorsal pallidum (PALd), fimbria (fi), and midbrain reticular nucleus (MRN) were negatively correlated with social interaction (Fig. 2b, c, and Table 1). The RD values in acc, pcc, aCP, pCP, fx, claustrum (CLA), PALd, fi, and MRN were positively correlated with social interaction (Fig. 2b, c, and Table 1). These findings suggested that the structural changes in white matter and striatal pathways were correlated with behavioral change in Olig1-Cre/Chd8L+/F mice.
Oligodendrocyte lineage-specific Chd8 heterozygous mutant mice manifest altered functional brain connectivity
Given that altered inter- and intrahemispheric functional connectivity has been detected in individuals with ASD , we next assessed long-range functional brain connectivity by measurement of the correlation of BOLD signals across brain regions in awake Olig1-Cre/Chd8L+/F and control mice with rsfMRI. Functional connectivity between pairs of cortical ROIs was compared between Olig1-Cre/Chd8L+/F and control mice. ROIs of 18 bilateral anatomic regions in the left and right hemispheres of the cortex were defined from our flexible annotation atlas based on the Allen Brain Atlas , and we constructed correlation matrices of functional connectivity between these ROIs in Olig1-Cre/Chd8L+/F and control mice (Fig. 3a). Whereas the extent of functional connectivity between right agranular insular cortex (AI) and left or right infralimbic cortex (ILA) was significantly increased in Olig1-Cre/Chd8L+/F mice compared with control mice (Fig. 3b), the extent of connectivity of corresponding ROIs between the right and left hemispheres did not differ significantly between the two genotypes (Fig. 3c). To examine the possible effect of GSR on these results, we also analyzed rsfMRI data without such regression. Without GSR, we did not detect any significant differences between the two genotypes with regard to the extent of functional connectivity between the right and left hemispheres (Additional file 3: Fig. S2).
To investigate further the effects of oligodendrocyte lineage-specific Chd8 heterozygous mutation on functional connectivity, we evaluated the connectivity between prelimbic cortex (PL), ILA, orbitofrontal cortex (ORB), dentate gyrus (DG), cornu ammonis 1 (CA1), subiculum (SUB), basolateral amygdala (BLA), CP, nucleus accumbens (ACB), ventral tegmental area (VTA), and periaqueductal gray (PAG), all of which are related to social behavior [45, 46]. We constructed correlation matrices of functional connectivity between these ROIs in Olig1-Cre/Chd8L+/F and control mice (Fig. 4a). The correlation between PL and DG, ILA and PAG, or ORB and PAG was increased in Olig1-Cre/Chd8L+/F mice, whereas the negative correlation between PL and BLA was diminished in Olig1-Cre/Chd8L+/F mice compared with control mice (Fig. 4b). We also found that the extent of functional connectivity between DG and CA1 was reduced in Olig1-Cre/Chd8L+/F mice compared with control mice (Fig. 4b). Although these differences did not survive multiple corrections, uncorrected statistics supported the alterations of functional connectivity in Olig1-Cre/Chd8L+/F mice. Similar results were obtained without GSR for the correlation between ORB and PAG and between PL and BLA (Additional file 3: Fig. S3). To evaluate the possible effect of head motion on functional connectivity, we calculated the framewise displacement from six parameters of head motion, which were estimated during the realignment process , and then examined the relation between framewise displacement and functional connectivity. Functional connectivity between PL and BLA did not show a significant correlation with head motion (Additional file 3: Fig. S4), suggesting that head motion during the acquisition of rsfMRI data did not affect the functional connectivity values. These results thus suggested that functional brain connectivity related to social behavior was altered in Olig1-Cre/Chd8L+/F mice.
Altered functional brain connectivity is correlated with structural and behavioral changes
We examined whether the structural changes detected in the brain of Olig1-Cre/Chd8L+/F mice were associated with the extent of functional connectivity. The medial prefrontal cortex (mPFC) including PL and ILA projects to the dorsomedial portions of CP in rodents , and structural changes in CP were apparent in Olig1-Cre/Chd8L+/F mice (Fig. 1). The RD value in CP was positively correlated with the extent of functional connectivity between PL and CP in Olig1-Cre/Chd8L+/F and control mice, although these differences did not survive multiple corrections (Fig. 5a). We also detected a correlation of functional connectivity between several other regions with structural changes in CP (Table 2). The fx is the major association pathway connecting the hippocampus with other regions including mPFC , and we detected a correlation between structural changes of fx and social interaction (Fig. 2c). We also found that FA and RD values in fx were correlated with the extent of functional connectivity between PL and CA1 among several other brain regions (Fig. 5b and Table 2).
We finally investigated whether altered functional connectivity in Olig1-Cre/Chd8L+/F mice was correlated with altered social interaction. The extent of functional connectivity between PL and the hippocampus including CA1, DG, and SUB was negatively correlated with social interaction, whereas the extent of functional connectivity between BLA and either PL or ORB was positively correlated with social interaction (Fig. 5c, d). We also detected a significant correlation between the extent of functional connectivity between several other regions (with or without GSR) and total interaction time during the reciprocal social interaction test (Fig. 5d, Table 2, and Additional file 3: Fig. S5). Among these results, the correlation between the extent of functional connectivity in PL-CA1, PL-SUB, and ORB-BLA and social interaction survived FDR corrections. Examination of the relation between functional connectivity and social interaction within the same genotype revealed that almost all significant correlations observed in the analysis of both genotypes were also apparent in Olig1-Cre/Chd8L+/F mice (Additional file 2: Table S2). These results suggested that altered functional connectivity in Olig1-Cre/Chd8L+/F mice is associated with both brain structural and behavioral changes.
Alterations in brain microstructure and functional connectivity are common characteristics of individuals with ASD as well as of mouse models of this disorder. To investigate which brain regions are altered in association with specific genetic mutations and whether these alterations contribute to autism-related behavior is key to understanding the etiology of ASD. We have now shown that oligodendrocyte-specific ablation of Chd8 results in changes to morphology and functional connectivity across the brain. We also found that altered functional connectivity between several brain regions was correlated with structural and behavioral changes.
Altered social interaction is a prominent characteristic of ASD. Three independent groups have reproducibly found that Chd8 heterozygous knockout mice manifest an increased contact time during the reciprocal social interaction test as well as altered social novelty preference during the three-chamber test [17, 18, 20]. Although increased contact time is not a typical behavioral characteristic of other ASD model mice, it may be a characteristic of social behavior in Chd8 mutant mice. Given that Chd8 mutation restricted to cells of the oligodendrocyte lineage recapitulated this behavioral phenotype , this mouse model appears to be well suited to analysis of functional brain connectivity underlying social behavior. Olig1-Cre/Chd8L+/F mice in the present study did not show an increased contact time in the reciprocal social interaction test (Additional file 3: Fig. S1), probably because of differences related to sample size, head bar attachment, surgery, habituation procedure, and environmental factors. However, our results have revealed that changes in brain structure or functional connectivity are correlated with social interaction in Olig1-Cre/Chd8L+/F mice.
Our DTI data revealed an increase in RD values and a reduction in FA values in MO, SS, ACA, and ac in Olig1-Cre/Chd8L+/F mice compared with control mice. Whereas FA is a summary measure of microstructural changes, RD reflects the structural status of myelin as well as axons . Given that both Chd8 heterozygous mutant mice and Olig1-Cre/Chd8L+/F mice manifest a reduced thickness of the myelin sheath , the altered FA and RD values of Olig1-Cre/Chd8L+/F mice might be attributable to myelination defects in these animals. It is of note that regions showing microstructural changes in Olig1-Cre/Chd8L+/F mice include the corticostriatal pathway, which plays a role in cognition, motivation, and reward processing  and is implicated in ASD pathogenesis . Striatal dysfunction has thus previously been demonstrated in Chd8 heterozygous mutant mice . Defects in synaptic transmission by medium spiny neurons have also been detected in Chd8 heterozygous mutant mice, and Chd8 ablation in ACB, but not in the dorsal striatum, was found to facilitate acquired motor learning . Microstructural changes due to oligodendrocyte defects might thus contribute to changes in neural activity and behavior in these mice. Although we have previously shown with the use of electron microscopy that the myelin sheath is thinner in cc of Olig1-Cre/Chd8L+/F mice , we did not detect any significant differences in FA or RD values in cc, likely as a result of the relatively small effect of Chd8 mutation on myelination or of technical limitations such as those related to the resolution or signal-to-noise ratio of DTI. We previously showed that myelination defects are present throughout the brain of Chd8 mutant mice including cortical regions , suggesting that myelination is also impaired in the regions that show alterations by DTI in the present study. Structural alterations in cc that are correlated with behavioral changes may be indicative of mild axonal or myelination defects.
The performance of behavioral analysis before MRI allowed us to assess the correlation between behavioral phenotype and structural or functional changes in the brain. Abnormal social behavior is the most prominent phenotype of individuals with ASD and is also observed in both Chd8 heterozygous mutant mice and Olig1-Cre/Chd8L+/F mice [17, 18, 20, 30]. Consistent with the previous finding that FA in cc and fx was negatively correlated with autism symptom severity as measured by the Autism Diagnostic Interview-Revised (ADI-R) in individuals with ASD , we found that FA and RD values in several brain regions including cc and fx were correlated with total contact time in the reciprocal social interaction test in Olig1-Cre/Chd8L+/F and control mice. Whereas microstructural changes in the cortex and striatum were thus detected in Olig1-Cre/Chd8L+/F mice, it remains unclear whether these changes are attributable to alterations in axons and dendrites in addition to myelination defects. Further analysis will be required to identify the mechanisms underlying such structural changes in several brain regions.
Action potential propagation is delayed as a result of myelination defects in Olig1-Cre/Chd8L+/F mice . However, it had remained unclear whether this delay affects neural networks in the brain of these mice. Given that Chd8 mutation results in a reduced conduction velocity in the cc connecting the left and right hemispheres , we investigated interhemispheric functional connectivity in Olig1-Cre/Chd8L+/F mice. However, we detected little effect of such Chd8 mutation on inter- or intrahemispheric functional connectivity, suggesting that the reduction in conduction velocity in the cc of Olig1-Cre/Chd8L+/F mice is not sufficient to influence functional connectivity in the examined regions. In contrast, we found that synchronization of rsfMRI signals between several brain regions related to ASD was altered in the mutant mice. Changes in functional connectivity in the mPFC, hippocampus, and amygdala were recently detected in individuals with ASD from four independent cohorts , with the amygdala having previously been shown to play a key role in social behavior and emotional information processing [53, 54]. The amygdala has reciprocal connections with multiple brain areas including the prelimbic prefrontal cortex, and our rsfMRI data showed that alterations in BLA connectivity in oligodendrocyte lineage-specific Chd8 heterozygous mutant mice were significantly correlated with total contact time in the reciprocal social interaction test. We also found that functional connectivity between PL and CA1 was associated with alterations in FA in fx and the behavioral deficit in the mutant mice. Excitatory neurons in the ventral hippocampus project to mPFC through fx, and these connections are related to social behavior [55, 56]. We speculate that delayed action potential propagation due to a myelination defect in fx might affect postsynaptic responses in mPFC and thereby give rise to altered social interaction in Olig1-Cre/Chd8L+/F mice. It is also possible that compensatory changes as a result of oligodendrocyte dysfunction during development might influence behavioral characteristics of Olig1-Cre/Chd8L+/F mice. For example, Nogo receptor ligand, a myelin glycoprotein expressed in oligodendrocytes, inhibits neurite outgrowth . It remains to be determined whether behavioral deficits in Olig1-Cre/Chd8L+/F mice depend on myelination defects or are secondary effects of oligodendrocyte dysfunction.
We applied GSR to eliminate undesired global noise and bias and to detect negative correlation. We performed awake rsfMRI, and movement-related noise is thought to be eliminated with the use of GSR. We were able to detect many changes in functional connectivity. We believe that the results obtained with GSR are thus more reliable in the present study. However, given the possibility that differences might be diminished with the use of GSR , we also provide the results of functional connectivity analysis without GSR (Additional file 3: Figs. S2, S3, and S5).
Whereas Chd8 heterozygous mutant mice were previously found to manifest increased connectivity in entorhinal, retrosplenial, and auditory cortical as well as posterior hippocampal areas , we detected no such alterations in these brain regions of Olig1-Cre/Chd8L+/F mice. This difference might be explained by a difference in affected cell types of the mutant mice (Chd8+/− versus Olig1-Cre/Chd8L+/F mice) or in experimental conditions (anesthetized or awake mice during rsfMRI). Furthermore, recent studies with various mouse models of ASD have detected changes in functional connectivity in many brain regions including PFC, MO, SS, ACA, AI, retrosplenial cortex, ventral hippocampus, CP, amygdala, hindbrain, thalamic nuclei, and medial-dorsal thalamus [24,25,26,27,28,29]. These ASD model mice as well as our mice overlap to only a small extent in their signatures of affected brain networks. Further studies will be necessary to investigate the causal relation between altered functional connectivity and ASD-like behaviors, as well as to elucidate how each genetic mutation gives rise to such connectivity changes in individuals with ASD or in model animals. Collectively, our results indicate that oligodendrocyte dysfunction induced by Chd8 mutation influences neuroanatomic and functional connectivity throughout the brain. They may thus be of translational value for further ASD research as well as provide insight into neural mechanisms underlying ASD.
Availability of data and materials
The data sets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Anterior cingulate cortex
Anterior corpus callosum
Agranular insular cortex
Autism spectrum disorder
Blood oxygen level-dependent
Chromodomain helicase DNA-binding protein 8
Diffusion tensor imaging
False discovery rate
Field of view
Global signal regression
Medial prefrontal cortex
Magnetic resonance imaging
Midbrain reticular nucleus
Posterior corpus callosum
Region of interest
Resting-state functional MRI
Effective time to echo
Time to repetition
Temporal signal-to-noise ratio
Ventral tegmental area
Amaral DG, Schumann CM, Nordahl CW. Neuroanatomy of autism. Trends Neurosci. 2008;31(3):137–45.
Stigler KA, McDonald BC, Anand A, Saykin AJ, McDougle CJ. Structural and functional magnetic resonance imaging of autism spectrum disorders. Brain Res. 2011;1380:146–61.
Travers BG, Adluru N, Ennis C, Trompdo PM, Destiche D, Doran S, Bigler ED, Lange N, Lainhart JE, Alexander AL. Diffusion tensor imaging in autism spectrum disorder: a review. Autism Res. 2012;5(5):289–313.
Di Martino A, Zuo XN, Kelly C, Grzadzinski R, Mennes M, Schvarcz A, Rodman J, Lord C, Castellanos FX, Milham MP. Shared and distinct intrinsic functional network centrality in autism and attention-deficit/hyperactivity disorder. Biol Psychiatry. 2013;74(8):623–32.
Uddin LQ, Supekar K, Menon V. Reconceptualizing functional brain connectivity in autism from a developmental perspective. Front Hum Neurosci. 2013;7:458.
O’Roak BJ, Vives L, Girirajan S, Karakoc E, Krumm N, Coe BP, Levy R, Ko A, Lee C, Smith JD, et al. Sporadic autism exomes reveal a highly interconnected protein network of de novo mutations. Nature. 2012;485(7397):246–50.
O’Roak BJ, Vives L, Fu W, Egertson JD, Stanaway IB, Phelps IG, Carvill G, Kumar A, Lee C, Ankenman K, et al. Multiplex targeted sequencing identifies recurrently mutated genes in autism spectrum disorders. Science. 2012;338(6114):1619–22.
Neale BM, Kou Y, Liu L, Ma’ayan A, Samocha KE, Sabo A, Lin CF, Stevens C, Wang LS, Makarov V, et al. Patterns and rates of exonic de novo mutations in autism spectrum disorders. Nature. 2012;485(7397):242–5.
Talkowski ME, Rosenfeld JA, Blumenthal I, Pillalamarri V, Chiang C, Heilbut A, Ernst C, Hanscom C, Rossin E, Lindgren AM, et al. Sequencing chromosomal abnormalities reveals neurodevelopmental loci that confer risk across diagnostic boundaries. Cell. 2012;149(3):525–37.
Nishiyama M, Oshikawa K, Tsukada Y, Nakagawa T, Iemura S, Natsume T, Fan Y, Kikuchi A, Skoultchi AI, Nakayama KI. CHD8 suppresses p53-mediated apoptosis through histone H1 recruitment during early embryogenesis. Nat Cell Biol. 2009;11(2):172–82.
Nishiyama M, Skoultchi AI, Nakayama KI. Histone H1 recruitment by CHD8 is essential for suppression of the Wnt-beta-catenin signaling pathway. Mol Cell Biol. 2012;32(2):501–12.
Sugathan A, Biagioli M, Golzio C, Erdin S, Blumenthal I, Manavalan P, Ragavendran A, Brand H, Lucente D, Miles J, et al. CHD8 regulates neurodevelopmental pathways associated with autism spectrum disorder in neural progenitors. Proc Natl Acad Sci USA. 2014;111(42):E4468-4477.
Cotney J, Muhle RA, Sanders SJ, Liu L, Willsey AJ, Niu W, Liu W, Klei L, Lei J, Yin J, et al. The autism-associated chromatin modifier CHD8 regulates other autism risk genes during human neurodevelopment. Nat Commun. 2015;6:6404.
Bernier R, Golzio C, Xiong B, Stessman HA, Coe BP, Penn O, Witherspoon K, Gerdts J, Baker C, Vulto-van Silfhout AT, et al. Disruptive CHD8 mutations define a subtype of autism early in development. Cell. 2014;158(2):263–76.
Wang J, Liu J, Gao Y, Wang K, Jiang K. Autism spectrum disorder early in development associated with CHD8 mutations among two Chinese children. BMC Pediatr. 2018;18(1):338.
Douzgou S, Liang HW, Metcalfe K, Somarathi S, Tischkowitz M, Mohamed W, Kini U, McKee S, Yates L, Bertoli M, et al. The clinical presentation caused by truncating CHD8 variants. Clin Genet. 2019;96(1):72–84.
Katayama Y, Nishiyama M, Shoji H, Ohkawa Y, Kawamura A, Sato T, Suyama M, Takumi T, Miyakawa T, Nakayama KI. CHD8 haploinsufficiency results in autistic-like phenotypes in mice. Nature. 2016;537(7622):675–9.
Platt RJ, Zhou Y, Slaymaker IM, Shetty AS, Weisbach NR, Kim JA, Sharma J, Desai M, Sood S, Kempton HR, et al. Chd8 mutation leads to autistic-like behaviors and impaired striatal circuits. Cell Rep. 2017;19(2):335–50.
Gompers AL, Su-Feher L, Ellegood J, Copping NA, Riyadh MA, Stradleigh TW, Pride MC, Schaffler MD, Wade AA, Catta-Preta R, et al. Germline Chd8 haploinsufficiency alters brain development in mouse. Nat Neurosci. 2017;20(8):1062–73.
Suetterlin P, Hurley S, Mohan C, Riegman KLH, Pagani M, Caruso A, Ellegood J, Galbusera A, Crespo-Enriquez I, Michetti C, et al. Altered neocortical gene expression, brain overgrowth and functional over-connectivity in Chd8 haploinsufficient mice. Cereb Cortex. 2018;28(6):2192–206.
Jung H, Park H, Choi Y, Kang H, Lee E, Kweon H, Roh JD, Ellegood J, Choi W, Kang J, et al. Sexually dimorphic behavior, neuronal activity, and gene expression in Chd8-mutant mice. Nat Neurosci. 2018;21(9):1218–28.
Ellegood J, Babineau BA, Henkelman RM, Lerch JP, Crawley JN. Neuroanatomical analysis of the BTBR mouse model of autism using magnetic resonance imaging and diffusion tensor imaging. Neuroimage. 2013;70:288–300.
Dodero L, Damiano M, Galbusera A, Bifone A, Tsaftsaris SA, Scattoni ML, Gozzi A. Neuroimaging evidence of major morpho-anatomical and functional abnormalities in the BTBR T+TF/J mouse model of autism. PLoS ONE. 2013;8(10):e76655.
Tsurugizawa T, Tamada K, Ono N, Karakawa S, Kodama Y, Debacker C, Hata J, Okano H, Kitamura A, Zalesky A, et al. Awake functional MRI detects neural circuit dysfunction in a mouse model of autism. Sci Adv. 2020;6(6):eaav4520.
Michetti C, Caruso A, Pagani M, Sabbioni M, Medrihan L, David G, Galbusera A, Morini M, Gozzi A, Benfenati F, et al. The knockout of synapsin II in mice impairs social behavior and functional connectivity generating an ASD-like phenotype. Cereb Cortex. 2017;27(10):5014–23.
Chelini G, Zerbi V, Cimino L, Grigoli A, Markicevic M, Libera F, Robbiati S, Gadler M, Bronzoni S, Miorelli S, et al. Aberrant somatosensory processing and connectivity in mice lacking engrailed-2. J Neurosci. 2019;39(8):1525–38.
Pagani M, Bertero A, Liska A, Galbusera A, Sabbioni M, Barsotti N, Colenbier N, Marinazzo D, Scattoni ML, Pasqualetti M, et al. Deletion of autism risk gene Shank3 disrupts prefrontal connectivity. J Neurosci. 2019;39(27):5299–310.
Haberl MG, Zerbi V, Veltien A, Ginger M, Heerschap A, Frick A. Structural-functional connectivity deficits of neocortical circuits in the Fmr1 (-/y) mouse model of autism. Sci Adv. 2015;1(10):e1500775.
Liska A, Bertero A, Gomolka R, Sabbioni M, Galbusera A, Barsotti N, Panzeri S, Scattoni ML, Pasqualetti M, Gozzi A. Homozygous loss of autism-risk gene CNTNAP2 results in reduced local and long-range prefrontal functional connectivity. Cereb Cortex. 2018;28(4):1141–53.
Kawamura A, Katayama Y, Nishiyama M, Shoji H, Tokuoka K, Ueta Y, Miyata M, Isa T, Miyakawa T, Hayashi-Takagi A, et al. Oligodendrocyte dysfunction due to Chd8 mutation gives rise to behavioral deficits in mice. Hum Mol Genet. 2020;29(8):1274–91.
Lu QR, Sun T, Zhu Z, Ma N, Garcia M, Stiles CD, Rowitch DH. Common developmental requirement for Olig function indicates a motor neuron/oligodendrocyte connection. Cell. 2002;109(1):75–86.
Hattori S, Okumura Y, Takao K, Yamaguchi Y, Miyakawa T. Open source code for behavior analysis in rodents. Neuropsychopharmacol Rep. 2019;39(1):67–9.
Yoshida K, Mimura Y, Ishihara R, Nishida H, Komaki Y, Minakuchi T, Tsurugizawa T, Mimura M, Okano H, Tanaka KF, et al. Physiological effects of a habituation procedure for functional MRI in awake mice using a cryogenic radiofrequency probe. J Neurosci Methods. 2016;274:38–48.
Komaki Y, Hikishima K, Shibata S, Konomi T, Seki F, Yamada M, Miyasaka N, Fujiyoshi K, Okano HJ, Nakamura M, et al. Functional brain mapping using specific sensory-circuit stimulation and a theoretical graph network analysis in mice with neuropathic allodynia. Sci Rep. 2016;6:37802.
Matsubayashi K, Nagoshi N, Komaki Y, Kojima K, Shinozaki M, Tsuji O, Iwanami A, Ishihara R, Takata N, Matsumoto M, et al. Assessing cortical plasticity after spinal cord injury by using resting-state functional magnetic resonance imaging in awake adult mice. Sci Rep. 2018;8(1):14406.
Takata N, Sugiura Y, Yoshida K, Koizumi M, Hiroshi N, Honda K, Yano R, Komaki Y, Matsui K, Suematsu M, et al. Optogenetic astrocyte activation evokes BOLD fMRI response with oxygen consumption without neuronal activity modulation. Glia. 2018;66(9):2013–23.
Hikishima K, Komaki Y, Seki F, Ohnishi Y, Okano HJ, Okano H. In vivo microscopic voxel-based morphometry with a brain template to characterize strain-specific structures in the mouse brain. Sci Rep. 2017;7(1):85.
Takata N, Sato N, Komaki Y, Okano H, Tanaka KF. Flexible annotation atlas of the mouse brain: combining and dividing brain structures of the Allen Brain Atlas while maintaining anatomical hierarchy. BioRxiv. 2020.
Murphy K, Bodurka J, Bandettini PA. How long to scan? The relationship between fMRI temporal signal to noise ratio and necessary scan duration. Neuroimage. 2007;34(2):565–74.
Abe Y, Komaki Y, Seki F, Shibata S, Okano H, Tanaka KF. Correlative study using structural MRI and super-resolution microscopy to detect structural alterations induced by long-term optogenetic stimulation of striatal medium spiny neurons. Neurochem Int. 2019;125:163–74.
Yano R, Hata J, Abe Y, Seki F, Yoshida K, Komaki Y, Okano H, Tanaka KF. Quantitative temporal changes in DTI values coupled with histological properties in cuprizone-induced demyelination and remyelination. Neurochem Int. 2018;119:151–8.
Le Bihan D, Mangin JF, Poupon C, Clark CA, Pappata S, Molko N, Chabriat H. Diffusion tensor imaging: concepts and applications. J Magn Reson Imaging. 2001;13(4):534–46.
Alexander AL, Lee JE, Lazar M, Field AS. Diffusion tensor imaging of the brain. Neurotherapeutics. 2007;4(3):316–29.
Anderson JS, Druzgal TJ, Froehlich A, DuBray MB, Lange N, Alexander AL, Abildskov T, Nielsen JA, Cariello AN, Cooperrider JR, et al. Decreased interhemispheric functional connectivity in autism. Cereb Cortex. 2011;21(5):1134–46.
Barak B, Feng G. Neurobiology of social behavior abnormalities in autism and Williams syndrome. Nat Neurosci. 2016;19(6):647–55.
Allsop SA, Vander Weele CM, Wichmann R, Tye KM. Optogenetic insights on the relationship between anxiety-related behaviors and social deficits. Front Behav Neurosci. 2014;8:241.
Power JD, Barnes KA, Snyder AZ, Schlaggar BL, Petersen SE. Spurious but systematic correlations in functional connectivity MRI networks arise from subject motion. Neuroimage. 2012;59(3):2142–54.
Gabbott PL, Warner TA, Jays PR, Salway P, Busby SJ. Prefrontal cortex in the rat: projections to subcortical autonomic, motor, and limbic centers. J Comp Neurol. 2005;492(2):145–77.
Mulder AB, Arts MP, Lopes da Silva FH. Short- and long-term plasticity of the hippocampus to nucleus accumbens and prefrontal cortex pathways in the rat, in vivo. Eur J Neurosci. 1997;9(8):1603–11.
Haber SN. Corticostriatal circuitry. Dialogues Clin Neurosci. 2016;18(1):7–21.
Balsters JH, Mantini D, Wenderoth N. Connectivity-based parcellation reveals distinct cortico-striatal connectivity fingerprints in Autism Spectrum Disorder. Neuroimage. 2018;170:412–23.
Holiga S, Hipp JF, Chatham CH, Garces P, Spooren W, D'Ardhuy XL, Bertolino A, Bouquet C, Buitelaar JK, Bours C, et al. Patients with autism spectrum disorders display reproducible functional connectivity alterations. Sci Transl Med. 2019;11(481).
Felix-Ortiz AC, Tye KM. Amygdala inputs to the ventral hippocampus bidirectionally modulate social behavior. J Neurosci. 2014;34(2):586–95.
Felix-Ortiz AC, Burgos-Robles A, Bhagat ND, Leppla CA, Tye KM. Bidirectional modulation of anxiety-related and social behaviors by amygdala projections to the medial prefrontal cortex. Neuroscience. 2016;321:197–209.
Phillips ML, Robinson HA, Pozzo-Miller L. Ventral hippocampal projections to the medial prefrontal cortex regulate social memory. Elife. 2019;8:e44182.
Sun Q, Li X, Li A, Zhang J, Ding Z, Gong H, Luo Q. Ventral hippocampal-prefrontal interaction affects social behavior via parvalbumin positive neurons in the medial prefrontal cortex. iScience. 2020;23(3):100894.
Wang KC, Koprivica V, Kim JA, Sivasankaran R, Guo Y, Neve RL, He Z. Oligodendrocyte-myelin glycoprotein is a Nogo receptor ligand that inhibits neurite outgrowth. Nature. 2002;417(6892):941–4.
Gotts SJ, Saad ZS, Jo HJ, Wallace GL, Cox RW, Martin A. The perils of global signal regression for group comparisons: a case study of Autism Spectrum Disorders. Front Hum Neurosci. 2013;7:356.
We thank C. Yamada, M. Goto, and M. Tanaka for technical assistance; laboratory members for discussion; T. Miyakawa (Fujita Health University) for providing the ImageSI application; and A. Ohta for help with preparation of the manuscript. A.K. was supported by a fellowship from the Japan Society for the Promotion of Science (JSPS).
This study was supported by Japan Society for the Promotion of Science, KAKENHI Grants (18H05215, 19H05220, and 16K07032) from the Ministry of Education, Culture, Sports, Science, and Technology of Japan.
Ethics approval and consent to participate
All methods were carried out in accordance with relevant guidelines and regulations, and the study was approved by the animal ethics committees of Kyushu University and the Central Institute for Experimental Animals.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Kawamura, A., Abe, Y., Seki, F. et al. Chd8 mutation in oligodendrocytes alters microstructure and functional connectivity in the mouse brain. Mol Brain 13, 160 (2020). https://doi.org/10.1186/s13041-020-00699-x
- Autism spectrum disorder
- Diffusion tensor imaging
- Functional connectivity