Introduction
22q11.2 copy number variants (CNVs) are rare, recurrent structural genetic variants that have profound impacts on neurodevelopment and confer increased risk for several developmental neuropsychiatric disorders, including autism spectrum disorder (ASD) and attention deficit/hyperactivity disorder (ADHD) [1,2,3,4]. The 22q11.2 deletion (22qDel) is estimated to occur in one in ~3700 live births [1], where 90% of cases are estimated to be due to de novo events [5]. The 22qDel is also associated with a range of medical and neurobehavioral outcomes, including congenital heart defects, immune deficiencies, and intellectual disability [1, 6, 7], and confers a markedly increased risk of psychosis [1,2,3,4]. The 22q11.2 duplication (22qDup) is estimated to occur in one in ~1600 live births [1] and, conversely, confers a putative reduced risk of psychosis [4, 8, 9]. Moreover, the 22qDup is more likely to be inherited rather than de novo, and may have a more variable clinical phenotype presentation [10, 11] compared to the 22qDel. Despite their divergent risk for psychotic illness, both 22qDel and 22qDup are associated with elevated rates of ASD, ADHD, and intellectual disability compared to controls [1, 3, 12], and disrupt dosage-sensitive genes that play fundamental roles during neurodevelopment. Thus, in contrast to research on individuals with behaviorally defined psychiatric disorders, studying a molecularly confirmed CNV like the 22q11.2 locus provides a genetically traceable and hypothesis-driven probe of brain alterations in individuals with a high risk of developing brain disorders. Moreover, reciprocal CNVs (i.e., gain and loss of genomic material at the same locus) allows for investigation of converging and diverging neurodevelopmental pathways of brain alterations by examining gene dosage effects on brain phenotypes.
22q11.2 CNVs have been associated with medium to large differences in several measures of brain structure, including grey and white matter [13,14,15,16], suggesting altered cellular architecture. However, the complex neuronal composition and dispersion of the cortical mantle of the brain makes it challenging to study the cellular architecture in MRI-derived measures of grey matter, such as cortical thickness. In contrast, the white matter consists of predominantly axons, which can be probed using diffusion MRI (dMRI) in vivo. As water diffusion is restricted or hindered by biological tissue, dMRI can be used to generate measures of white matter microstructure that is influenced by the cellular architecture along white matter tracts [17,18,19]. For instance, water molecules in an environment of dense and coherent white matter axons will diffuse along the primary direction of the axons, which has been commonly measured by diffusion tensor imaging (DTI)-derived metrics such as fractional anisotropy (FA, i.e., standardized value of anisotropic diffusion) and mean diffusivity (i.e., measure of diffusion in all directions).
A previous consortium-based DTI-study of 22qDel carriers reported an overall pattern of higher FA and lower diffusivity across most white matter tracts in 22qDel carriers compared to controls [13], which could reflect an increase in the cumulative cellular membrane circumference, increase in smaller tortuous axons [13], and/or higher density of axons with disproportionately small diameters [13, 20] compared to controls. Indeed, other studies have suggested other lines of evidence for increased axonal density in 22qDel carriers, such as excessive prenatal overgrowth of thalamocortical axons in a cellular model of the 22qDel [21] and a higher frequency of deep layer cortical neurons in a post-mortem examination of a 3-month-old infant with 22qDel [22]. In contrast to 22qDel, only one prior small dMRI study has included 22qDup carriers [23]. Despite the modest sample size, the results indicated widespread lower FA and higher diffusivity in 22qDup carriers compared to controls [23]. Although the biological underpinnings of these differences are unknown, the opposing pattern in 22qDel and 22qDup carriers suggests a gene dosage effect, which may reflect divergent effects of the 22q11.2 CNV on the cellular architecture of white matter tracts. Indeed, the 22q11.2 locus is a hotspot for genes critical to white matter development, including axonal morphology and function in white matter such as regulation of axonal diameter, growth and branching (e.g., ZDHHC8 [24, 25], TBX1 [26], RTN4R/Nogo-66 [27], TXNRD2 [28]), as well as dendritic morphology and neuronal excitability (e.g., DGCR8 [29, 30]). In addition, both 22qDel and 22qDup affect expression of genes outside the 22q11.2 locus [21, 29, 31, 32], including genes that influence axonal morphology [21]. However, to our knowledge, no previous studies have directly compared axonal density in white matter tracts between 22qDel and 22qDup carriers to determine if a gene dosage effect exists. Moreover, such effects may be region-specific due to the unique developmental trajectories of white matter fibers. For example, limbic, commissural and projection fibers develop and mature earlier than association fibers [33,34,35,36]. Indeed, it has been postulated that the nature of the white matter alterations in 22qDel carriers may differ between commissural and association fibers [13], potentially influenced by the maturational patterns of these fiber tracts. Specifically, the presence of commissural fibers are observed in the first trimester to second trimester [37,38,39,40,41,42], whereas the emergence of association fibers occurs later during the second trimester [38, 39, 42, 43]. The number of axons peak around birth for commissural fibers [44, 45], followed by axonal elimination[45, 46], axon growth and myelination from birth through adolescence [42, 47, 48]. The association fibers are the least mature fibers during infancy [49] and are the last to mature, as indicated by developmental changes in DTI measures well into adulthood [36]. Thus, it is plausible that the 22q11.2 CNV may have maturation-specific effects on commissural and association tracts.
Conventional DTI measures are inherently nonspecific to the underlying white matter cellular architecture [18, 50, 51] and can be influenced by alterations to axonal density and/or dispersion that can be driven by early (e.g., axonal proliferation or axonal pruning) or later (e.g., axonal diameter expansion) neurodevelopmental events. In contrast, neurite orientation dispersion and density imaging (NODDI) is a dMRI method that is better suited to test the hypothesis of altered axonal density in 22q11.2 CNV carriers, as it provides measures of the biophysical properties of white matter tracts, including axonal density (i.e., intracellular volume fraction, ICVF), axonal dispersion (i.e., orientation dispersion index, ODI), and free water diffusion/cerebrospinal fluid (i.e., isotropic water diffusion, ISO). It is important to note that NODDI generates metrics that provide indirect proxies of underlying microstructural properties and cannot directly disentangle the specific biological events (e.g., along-axon diameter variations and/or glial changes) that drive these signal differences [52]. Nevertheless, as a novel technique that disentangles the composite signal of FA derived from DTI, NODDI metrics allow for more precise and biologically meaningful inferences about the underlying white matter architecture. Moreover, NODDI-derived measures have been found to correlate with their histological counterparts [53,54,55] and to be altered in several neuropsychiatric disorders, including psychosis and ASD [56]. Thus, examining NODDI measures in 22q11.2 CNV carriers may be informative for detecting converging neurobiological patterns with idiopathic developmental neuropsychiatric disorders. In addition, ICVF has been found to be more sensitive to age-related differences in white matter microstructure compared to DTI-derived metrics during development [57]. As such, this metric can yield novel insight into age-specific neurobiological events that underlie changes in white matter microstructure [58,59,60], such as prenatal proliferation of axons[44, 45] and prolonged periods of axonal diameter expansion across childhood and adolescence [58]. Thus, by mapping the neurodevelopmental trajectories of the axonal architecture across white matter regions, we can obtain important insight into the maturational mechanisms that underlie altered white matter observed in 22qDel and 22qDup carriers.
Here, we aim to characterize the effects of 22q11.2 CNVs on axonal density and axonal dispersion in white matter regions. We hypothesize gene-dosage effects on axonal density in white matter regions such that 22qDel carriers will exhibit higher and 22qDup carriers will exhibit lower axonal density, respectively. We also aim to investigate spatiotemporal effects of the 22q11.2 CNV on axonal architecture by examining age-related trajectories using mixed cross-sectional and longitudinal data, and by leveraging a large dataset of typically developing controls from the Human Connectome Project (HCP) [61, 62]. We reasoned that alterations to prenatal neurodevelopmental processes will manifest as robust group differences in NODDI measures, whereas alterations to later neurodevelopmental processes will be reflected by altered age-related trajectories in NODDI measures across childhood and adolescence.
Methods
Participants
A total of 964 participants with 1004 scans were included in the analysis, including 50 22qDel carriers, 24 22qDup carriers and 890 age and sex matched typically developing controls (Table 1). The 22qDel carriers were recruited from either the University of California, Los Angeles (UCLA) or Stanford University, whereas the 22qDup carriers were recruited from UCLA only (supplementary note 1). The control group included unrelated healthy control participants recruited from UCLA and Stanford University, and healthy control participants derived from the HCP, including HCP-Development [61] and HCP-Aging [62] to establish robust developmental trajectories from childhood to adulthood (see sTable 1–2 for baseline characteristics across scanner sites, sTable 3 for binned age groups, sFigures 1–2).
Diffusion MRI acquisition, preprocessing and quality control
Multi-shell dMRI was acquired using a 3 T Siemens Prisma scanner with the following parameters: TE = 89.2 ms, TR = 3230 ms, flip angle = 78°, slice thickness = 1.5 mm, voxel size = 1.5 mm isotropic obtained with 7 b = 0 s/mm2 images, and 3 b = 200 s/mm2, 6 b = 500 s/mm2, 46 b = 1500 s/mm2 and 46 b = 3000 s/mm2 diffusion weighted volumes. The scans were acquired with both anterior to posterior and posterior to anterior phase encoding directions, resulting in a total of 216 volumes, following the HCP acquisition protocol [63]. Manual inspections and statistical analyses were conducted to exclude participants with excessive head movement (see supplementary note 1 for details). We also included dMRI data derived from the HCP-D (https://www.humanconnectome.org/study/hcp-lifespan-development) and HCP-A (https://www.humanconnectome.org/study/hcp-lifespan-aging) as described elsewhere [61, 62]. All the raw dMRIs were denoised with LPCA and corrected for Gibbs ringing artifacts using DiPy [64, 65]; susceptibility induced distortions, head movement, and eddy currents using TOPUP and EDDY in FSL [66, 67]. Afterwards, we ran the diffusion preprocessing pipeline from HCP (https://github.com/Washington-University/HCPpipelines). NODDI was fitted using the dMIPY package (https://github.com/AthenaEPI/dmipy). The NODDI derived scalar maps were registered to the ENIGMA-DTI template [68] and average skeletonized ROI-based measures were calculated based on the Johns Hopkins University White Matter atlas [69] following the ENIGMA-DTI protocol (https://enigma.ini.usc.edu/protocols/dti-protocols/) [68], including anterior corona radiata (ACR), anterior limb of internal capsule (ALIC), body of corpus callosum (BCC), corpus callosum (CC), cingulum cingulate gyrus (CGC), cingulum hippocampus (CGH), corona radiata (CR), corticospinal tract (CST), extreme/external capsule (EC), fornix (FX), fornix (crus)/stria terminalis (FXST), genu of corpus callosum (GCC), internal capsule (IC), posterior corona radiata (PCR), posterior limb of the internal capsule (PLIC), posterior thalamic radiation (PTR), retrolenticular part of the internal capsule (RLIC), splenium of corpus callosum (SCC), superior corona radiata (SCR), superior fronto-occipital fasciculus (SFO), superior longitudinal fasciculus (SLF), sagittal stratum (SS), tapetum of the corpus callosum (TAP), uncinate fasciculus (UNC) (see sFigure 3 for a brain map). Due to repeated measures across scanner sites, we harmonized all of the NODDI measures using longCombat [70]. Due to the focus on white matter regions, which mostly consist of myelinated axons [71], we define ICVF as a measure of axonal density, ODI as a measure of axonal dispersion and ISO as a measure of free water diffusivity [51].
Statistical analyses
Group comparisons
We used generalized additive mixed models (GAMM) in R [72] to test for group differences in ICVF, ISO and ODI across 27 variables (24 regions of interest (ROIs) and 3 average measures). Here, we compared 22qDel carriers to controls, 22qDup carriers to controls, and 22qDel carriers to 22qDup carriers. In addition, we ran gene dosage analyses on NODDI measures by treating copy number as a continuous variable (22qDel = 1, control = 2, and 22qDup = 3). All group analyses were adjusted for the smoothed age effect using cubic regression splines, sex, and repeated measures through the inclusion of participant ID as random intercept. All GAMM analyses were fitted using default arguments, where k = -1 is set for the smoothed age terms and the use of restricted maximum likelihood when fitting all models. The standardized beta is reported and used as a measure of effect size. To account for multiple testing, we applied FDR correction across 75 (i.e., dependent variables) * 4 (i.e., three group analyses and one gene dosage analysis) = 300 comparisons. Finally, we conducted several sensitivity analyses to test the robustness of the results by including 1) participants from the same scanner site only, and by adjusting for 2) cerebral white matter volume, 3) intracranial volume or 4) cerebrospinal fluid volume (see supplementary note 2 for details).
Developmental trajectories
To examine the developmental trajectories of the NODDI measures, we examine the smoothed effect of age, with cubic regression splines using GAMM, in controls, 22qDel carriers and 22qDup carriers. All age-related analyses were conducted on site-harmonized data, and adjusted for the fixed effect of sex and the random effect of subject ID. First, we present the smoothed effect of age for the control group, along with the corresponding FDR-corrected p-values specific to the control group analysis, to establish typical developmental trajectories across white matter regions. The effective degrees of freedom (edf) represent the complexity of the smoothed age term, where values close to 1 indicate a linear relationship and higher values indicate non-linear relationships between chronological age and the response variable. We then compared the developmental trajectories of 22qDel and 22qDup carriers to the control group. A significant interaction term between the smoothed-age term and 22qDel or 22qDup indicates a difference in developmental trajectory compared to controls. To account for multiple testing, we applied FDR-correction across all the interaction terms (i.e., age*22qDel and age*22qDup).
Results
Gene dosage effects of the 22q11.2 CNV on intracellular volume fraction
Mean ICVF was significantly greater in 22qDel carriers compared to controls, and significantly greater across all regions except for FX and CGH. In contrast, 22qDup carriers showed significantly lower mean ICVF, and lower ICVF across all white matter regions compared to controls. Reflecting a similar pattern, there was a negative gene dosage effect of the 22q11.2 CNV (i.e., decreasing ICVF with increasing copy number) on mean ICVF and regional ICVF across all white matter regions (Fig. 1, left panel, sTable 3). Sensitivity analyses revealed that scanner site, cerebral white matter volume, intracranial volume, and cerebrospinal fluid volume did not influence these results (sFigure 4, sTable 5–8).
The control group and the 22qDel group are used as reference groups for the group comparisons. Gene dosage effects are based on a continuous variable of copy number (i.e., 22qDel = 1, control = 2, and 22qDup = 3). E.g., Significant negative gene dosage effects indicate that an increased number of copies of genes in the 22q11.2 genomic locus is associated with lower values of the diffusion measures. ACR anterior corona radiata, ALIC anterior limb of internal capsule, BCC body of corpus callosum, CC corpus callosum, CGC cingulum cingulate gyrus, CGH cingulum hippocampus, CR corona radiata, CST corticospinal tract, EC extreme/external capsule, FX fornix. FXST fornix (crus)/stria terminalis, GCC genu of corpus callosum, IC internal capsule, PCR posterior corona radiata, PLIC posterior limb of the internal capsule, PTR posterior thalamic radiation, RLIC retrolenticular part of the internal capsule, SCC splenium of corpus callosum, SCR superior corona radiata, SFO superior fronto-occipital fasciculus, SLF superior longitudinal fasciculus, SS sagittal stratum, TAP tapetum of the corpus callosum, UNC uncinate fasciculus.
Higher mean, but regionally variable effects, of 22q11.2 CNV on orientation dispersion index
Both 22qDel and 22qDup carriers showed greater mean ODI compared to controls; however, there were region-specific differences. 22qDel carriers had higher ODI in limbic (i.e., FXST, CGH), commissural (i.e., GCC) and association (i.e., SS, SLF, SFO, EC) fiber tracts, but also directionally variable effects on projection tracts (i.e., higher ODI in ACR and PTR; lower ODI in IC, PLIC, ALIC and CST) compared to controls. 22qDup carriers showed higher ODI in projection (i.e., IC, PLIC, ALIC, CST), commissural (i.e., BCC) and association (i.e., CGC, UNC, EC) fiber tracts, but also lower ODI in PCR, where projections fibers predominate, compared to controls (Fig. 1, middle panel, sTable 3). The sensitivity analyses revealed the same pattern, albeit with fewer significant results; opposing effects of the 22qDel and 22qDup on ODI in projection (i.e., PTR, PLIC, IC) and association regions (i.e., SS and SLF) yielded the most robust effects (sFigure 5, sTable 5–8).
Reduced free water diffusion in 22qDup carriers
There was no significant difference in mean ISO between 22qDel and controls. Regionally, however, 22qDel carriers exhibited higher ISO in a few projection tracts (i.e., PCR, ACR, PTR), and directionally variable effects on limbic (i.e., higher in FX; lower in FXST) commissural fibers (i.e., higher in TAP; lower in SCC) compared to controls (Fig. 1, right panel). Sensitivity analyses revealed only higher ISO in the FX and PTR in 22qDel compared to controls after adjusting for intracranial volume, white matter volume (significant for FX and PTR), and cerebrospinal fluid volume (significant for FX; sFigure 3, sTable 3). 22qDup carriers showed lower mean ISO compared to controls, primarily driven by lower ISO in white matter regions predominated by projection fibers (RLIC, CR, PCR, SCR, ACR, PTR) and association fibers (UNC, SS, SLF, SFO, Fig. 1, right panel). However, similar to the 22qDel carriers, 22qDup carriers also exhibited higher ISO in FX and lower ISO in FXST. Overall, sensitivity analyses yielded the same pattern of significant group-level differences between 22qDup and controls, except for the FX, UNC, RLIC and ACR that were no longer significant (sFigure 6, sTable 5–8).
22q11.2 CNV carriers show typical developmental trajectories in axonal density and dispersion
In controls, there were significant age effects on ICVF, ODI and ISO across all white matter regions from childhood to adulthood, except for ODI in the TAP (sTable 9). Most of the measures showed a non-linear age effect across white matter regions, except for ODI in ALIC, CST, EC, SFO, SLF and SS, which showed a linear age effect indicated by their effective degrees of freedom (i.e., edf = 1). The age-related trajectories of 22qDel and 22qDup carriers did not statistically differ from the age-related trajectories of the control group for the global and regional measures of ICVF (Fig. 2), ODI (Fig. 3) or ISO (Fig. 4) after adjusting for multiple comparisons (sTable 10). Follow-up analyses also revealed that several group differences, most prominently in ICVF, are already evident in childhood and persist through adolescence and into adulthood (sFigure 7, sTable 11).
ACR anterior corona radiata, ALIC anterior limb of internal capsule, BCC body of corpus callosum, CC corpus callosum, CGC cingulum cingulate gyrus, CGH cingulum hippocampus, CR corona radiata, CST corticospinal tract, EC extreme/external capsule, FX fornix. FXST fornix (crus)/stria terminalis, GCC genu of corpus callosum, IC internal capsule, PCR posterior corona radiata, PLIC posterior limb of the internal capsule, PTR posterior thalamic radiation, RLIC retrolenticular part of the internal capsule, SCC splenium of corpus callosum, SCR superior corona radiata, SFO superior fronto-occipital fasciculus, SLF superior longitudinal fasciculus, SS sagittal stratum, TAP tapetum of the corpus callosum, UNC uncinate fasciculus.
ACR anterior corona radiata, ALIC anterior limb of internal capsule, BCC body of corpus callosum, CC corpus callosum, CGC cingulum cingulate gyrus, CGH cingulum hippocampus, CR corona radiata, CST corticospinal tract, EC extreme/external capsule, FX fornix. FXST fornix (crus)/stria terminalis, GCC genu of corpus callosum, IC internal capsule, PCR posterior corona radiata, PLIC posterior limb of the internal capsule, PTR posterior thalamic radiation, RLIC retrolenticular part of the internal capsule, SCC splenium of corpus callosum, SCR superior corona radiata, SFO superior fronto-occipital fasciculus, SLF superior longitudinal fasciculus, SS sagittal stratum, TAP tapetum of the corpus callosum, UNC uncinate fasciculus.
ACR anterior corona radiata, ALIC anterior limb of internal capsule, BCC body of corpus callosum, CC corpus callosum, CGC cingulum cingulate gyrus, CGH cingulum hippocampus, CR corona radiata, CST corticospinal tract, EC extreme/external capsule, FX fornix. FXST fornix (crus)/stria terminalis, GCC genu of corpus callosum, IC internal capsule, PCR posterior corona radiata, PLIC posterior limb of the internal capsule, PTR posterior thalamic radiation, RLIC retrolenticular part of the internal capsule, SCC splenium of corpus callosum, SCR superior corona radiata, SFO superior fronto-occipital fasciculus, SLF superior longitudinal fasciculus, SS sagittal stratum, TAP tapetum of the corpus callosum, UNC uncinate fasciculus.
Discussion
Hitherto, this is the largest study to examine biophysical dMRI-derived measures in brain white matter among 22q11.2 CNV carriers. The results showed widespread higher ICVF in 22qDel and lower ICVF in 22qDup compared to controls, where gene dosage at the 22q11.2 locus was associated with lower ICVF across all white matter regions, suggesting a strong effect of the 22q11.2 CNV on axonal density. Further, we found higher mean axonal dispersion (ODI for both 22qDel and 22qDup compared to controls, but also directionally variable effects of the 22q11.2 CNV on ODI across white matter regions; overall, having more copies of 22q11.2 genes was associated with lower ODI in regions where association, projection and limbic fibers predominate. However, we also observed that more copies at the 22q11.2 genomic locus were associated with higher ODI in IC and CST regions. In addition, we observed lower ISO in 22qDup carriers compared to controls in several white matter tracts where projection and association fibers predominate. Finally, we find that 22qDel and 22qDup carriers showed similar developmental trajectories to controls in all NODDI measures across childhood and adolescence.
The current study extends on the previous findings of higher FA in 22qDel [13] and lower FA in 22qDup [23] relative to typically developing controls, by focusing on measures of two of the main contributors to FA: axonal density and axonal dispersion [51]. Across regions, we found robust evidence for higher axonal density in 22qDel carriers and lower axonal density in 22qDup carriers compared to controls, suggesting that axonal density is the main contributor to differences in FA alterations in white matter microstructure among 22q11.2 CNV carriers. It should be noted, however, that some white matter regions have previously been found to show lower (i.e., FXST, SLF, EC) or no difference in FA (i.e., RLIC, UNC, SS, PTR, SFO) in 22qDel carriers compared to controls [13]. Here, we find that all these regions, except for the RLIC and UNC, were characterized by both higher axonal density and higher axonal dispersion in 22qDel carriers compared to controls. Thus, it seems that higher axonal dispersion counteracts the influence of higher ICVF, yielding lower or similar FA relative to controls in 22qDel carriers.
The greatest dosage effects of the 22q11.2 genomic region on ICVF were found in commissural tracts and the corona radiata, which consist of fibers in the corpus callosum and projection fibers from motor cortex, respectively. Such strong opposing effects on white matter architecture may be driven by mechanisms that are directly related to the deletion and duplication of genes within the 22q11.2 locus. However, the exact mechanisms underlying the dosage effect on these tracts are unknown, but it should be noted that the corpus callosum has been found to be larger in children with the 22qDel [73], which contrasts with the overall pattern of lower brain volume in 22qDel carriers [14, 15]. As corpus callosum is made up by subregions where the size of the axons typically differ in humans, e.g., the relative density of thick and thin axons (i.e., the axon diameter distribution) differs across the corpus callosum [74,75,76], it seems plausible that axonal overgrowth or axonal diameter expansion may underlie alterations to axonal density of commissural fibers and corpus callosum volume. Indeed, deletion of genes within the 22q11.2 locus (e.g., ZDHHC8 and RTN4R/Nogo-66) has been linked to axonal growth and branching in mice models of the 22qDel [24, 27]. Others have shown that heterozygous deletion of the TBX1 among mice can result in relatively more axons with small diameters and no large diameters in the fimbria compared to wild-type mice [26], which is in line with the predictions from previous diffusion MRI studies on 22q11.2 deletion carriers [13, 20].
Our results also revealed higher ODI in 22qDel carriers for a few regions that are predominated by association fibers. Lower FA has previously been found in these white matter regions in 22qDel compared to controls [13]. While lower FA can indicate both lower axonal density and/or higher axonal dispersion/crossing fibers [51], our result of increased ODI indicates that this is primarily driven by higher axonal dispersion in association tracts in 22qDel carriers relative to controls. Moreover, long association fibers typically originate from upper cortical layers and show protracted development compared to limbic and projection fibers [33,34,35,36], which may indicate that the white matter fibers originating from upper cortical layers show different morphology compared to fibers originating from deeper cortical layers. This interpretation aligns with several lines of previous research, including a mouse model of 22qDel showing fewer neural progenitors [77], shorter axons, and less dendritic branching in neurons from the upper cortical layers relative to wild type mice [28]. The latter finding was observed after deletion of the TXNRD2 gene, which also resulted in fewer long distance cortical connections and more local connections, which may yield higher axonal dispersion for long distance/association fibers in 22qDel carriers (e.g., lower alignment and/or rerouting through compensation of alternative or less direct paths). In addition, a post-mortem study reported lower frequency of neurons in the upper cortical layers in a 22qDel carrier [22], and a neuroimaging study reported altered axonal morphology (i.e., utilizing ultra-strong dMRI gradients sensitive to axonal morphology) and lower white matter volume in association tracts in 22qDel carriers compared to controls [20], supporting this interpretation.
Overall, the results point towards a divergent pattern of white matter microstructure between 22qDel and 22qDup carriers; however, there were a few points of convergence. Such diverging and converging effects are of particular interest in research on genomic structural variants, where gene dosage effects are often observed for neuroimaging derived features [16, 23, 78] but where the clinical phenotypic profile may converge (e.g., increased risk of ASD for both 22qDel and 22qDup carriers [12]). In the current study, we observed that a few regions exhibit similar directions of effect for both 22qDel and 22qDup, including higher ODI in EC, higher ISO in FX and lower ISO in FXST. These results indicate higher axonal dispersion in EC in both 22qDel and 22qDup, whereas the differences in FX and FXST may be driven by CSF partial volume in the voxels at the fornix and adjacent areas. White matter alterations in the external capsule have also been associated with idiopathic ASD (i.e., lower FA in twins with autism compared to control twins) [79] and symptom severity, and it is one of the white matter regions most susceptible to environmental factors in individuals with ASD [79]. Moreover, a recent pilot study found within-subject white matter changes in the external capsule among toddlers with ASD after a behavioral intervention that improved verbal and communication skills [80], possibly implicating a role of the external capsule in individuals with or at risk for speech and communication problems. Our results suggest altered axonal circuitry in the external capsule may be implicated in the behavioral phenotypes of ASD commonly observed in 22q11.2 CNV carriers, although future research on the relationship between external capsule white matter structure and autistic phenotypes in 22q11.2 CNV carriers is required to understand this mechanism.
The 22qDel is also one of the greatest genetic risk factors for psychosis [4, 81]. A previous neuroimaging study has found evidence for lower axonal density in commissural and association fibers and higher axonal density and free water in the anterior thalamic radiation, predominated by projection fibers, among individuals with first episode psychosis[56]. The differences in thalamic radiation converge with our results of higher ICVF and ISO in the PTR among 22qDel carriers compared to controls. A recent study reported thalamocortical axonal overgrowth in 22qDel thalamocortical organoids, likely due to elevated FOXP2 expression as a trans-effect of the 22qDel [21]. In addition, hyperconnectivity in thalamic circuitry has been reported among 22qDel carriers compared to controls [82,83,84], which has also been hypothesized as a potential biomarker for psychosis [85]. It is also noteworthy that the 22qDup, carriers of a putative protective genetic factor for psychosis [4, 8, 9], showed the opposite pattern in our study, reflected by lower ICVF and ISO in the PTR. This may suggest that a higher axonal density of thalamic projection fibers, possibly due to thalamocortical axonal overgrowth, is implicated in the etiology of psychosis.
Finally, despite widespread group differences between 22q11.2 CNV carriers and controls, we did not find strong support for 22q11.2 CNV-dependent differences in the developmental trajectories of NODDI measures. Thus, the results do not support alterations in late maturational changes in axonal architecture, such as cumulative axonal diameter expansion. These findings contrast with the notably altered developmental pattern of resting-state functional connectivity observed in 22qDel carriers compared to controls [82, 86]. Concomitant with these findings, a parallel investigation in a mouse model of 22qDel found evidence for developmental changes in dendritic spine density [86]. Taken together, these results suggest that changes in synapses and spine density may reshape brain circuits during adolescence, while altered axonal packing, as reflected by ICVF, may be established earlier. Still, it is important to note that several factors may contribute to these contrasting findings, including limited statistical power to detect interactions, high within-group variability, and/or imperfect biophysical specificity of the NODDI measures. Moreover, the observed alterations in axonal density may instead be reflective of atypical prenatal axonogenesis or lack of early axonal pruning in 22qDel and 22qDup carriers. This is in line with previous research on 22qDel-derived cortical spheroids exhibiting excessive prenatal overgrowth of thalamocortical axons [21] and neural progenitors that stay longer in the cell cycling state [87], which may be smaller in size, as indicated by reduced neurospheres derived from human induced pluripotent stem cells from 22qDel carriers [88]. To speculate, the gene dosage effect of 22q11.2 genomic region on axonal density observed in our study may be reflective of a differential effect on gene expression levels of genes that are involved in axonogenesis, possibly through trans-regulatory mechanisms such as the mediating effect of FOXP2 gene expression levels on 22qDel thalamic neurons [21]. However, to our knowledge, no studies have utilized cellular models of the reciprocal 22qDup; such future studies are warranted to determine differential gene dosage effects on axonal development.
Strengths and limitations
The current study includes a unique sample across a wide age-range consisting of both 22qDel and 22qDup carriers and a large typically developing control sample with multi-shell dMRI data, allowing for gene dosage and developmental analyses on metrics derived from advanced neuroimaging techniques. The use of advanced dMRI measures provides novel insight into the axonal architecture of the white matter microstructure in 22q11.2 CNV carriers, insight that cannot be achieved through conventional DTI measures [18, 51]. We also performed several sensitivity analyses demonstrating the robustness of the results of the group level analyses by analyzing individuals derived from the same scanner site only and adjusting for brain volumetric differences. The results revealed largely the same pattern, albeit with some group differences becoming attenuated due to reduction in sample size. It is also important to note that small sample sizes—due to the population prevalence of these CNVs—limit our ability to statistically detect small effects. Despite this, previous studies have found robust differences in white matter microstructure with smaller sample sizes than our study [20], which remains the largest study on biophysical white matter measures in 22q11.2 CNV carriers to date. Here, we also utilized GAMMs to capture nonlinear age effects, account for repeated measures and avoid biased group comparisons, especially important in instances where sample sizes are modest and age distributions differ between groups. Although we did not find evidence for altered developmental trajectories in 22q11.2 CNV carriers, we urge caution when interpreting these results, as the interaction effects may be more subtle in magnitude than the corresponding main effects and thus require greater power to detect. Thus, we cannot exclude the possibility of subtly altered developmental trajectories in NODDI measures among 22q11.2 CNV carriers compared to controls. However, because robust group differences in white matter microstructure between 22qDel carriers and controls are already evident in early adolescence [20], and neither our analyses nor a previous multi-site DTI study of 22qDel vs controls [13] show notable age-by-group interactions, there is currently no clear evidence that these group differences are driven by large maturational shifts during adolescence. Future studies with larger sample sizes, dense longitudinal data and methods that examine differences in smoothed confidence intervals (e.g., [16, 89]) may shed new light on possible age-dependent differences across the lifespan in 22q11.2 CNV carriers. Finally, despite the advantages of NODDI, it is still subject to the limitations inherent to dMRI, which is used to infer tissue microstructure by measuring water diffusion within a voxel of tissue. In addition, although NODDI measures have been shown to reflect their histological counterparts [53,54,55], it is important to note that biophysical models are emerging techniques that remain under active development [51, 90,91,92,93,94]. Our results also raise new questions. For instance, it is unclear if the altered axonal density in 22q11.2 CNV carriers is a consequence of an absolute difference in total neurons in the brain, altered axonal diameter distribution, or other mechanisms that alter space requirements in white matter tissue. Future studies utilizing both in-vitro cellular models and post-mortem examinations of 22q11.2 CNV carriers are warranted.
Conclusions
The results of the current study provide new insights into the underlying neurobiology of the altered white matter microstructure in 22q11.2 CNV carriers, linking copy numbers at the 22q11.2 genomic locus to altered axonal density. However, we do not find evidence for altered age-related changes in axonal density or dispersion across childhood and adolescence, possibly reflecting atypical axonogenesis and/or axonal pruning during early neurodevelopment in 22q11.2 CNV carriers. Future studies are warranted to characterize the morphological features of white matter axons across fetal development, and to connect these developmental alterations to neuropsychiatric phenotypes.
Data availability
Phenotypic data are publicly available from the National Institute of Mental Health Data Archives: https://nda.nih.gov/edit_collection.html?id=2414. Neuroimaging data included in this study are available from the corresponding author upon reasonable request and subject to a data-sharing agreement. Data from the Human Connectome Project database are publicly available to qualified researchers upon registration and agreement to data use terms (https://www.humanconnectome.org).
References
-
Olsen L, Sparsø T, Weinsheimer SM, Dos Santos MBQ, Mazin W, Rosengren A, et al. Prevalence of rearrangements in the 22q11.2 region and population-based risk of neuropsychiatric and developmental disorders in a Danish population: a case-cohort study. Lancet Psychiatry. 2018;5:573–80.
-
Schneider M, Debbané M, Bassett AS, Chow EWC, Fung WLA, van den Bree MBM, et al. Psychiatric disorders from childhood to adulthood in 22q11.2 deletion syndrome: results from the international consortium on brain and behavior in 22q11.2 deletion syndrome. Am J Psychiatry. 2014 June;171:627–39.
-
Hoeffding LK, Trabjerg BB, Olsen L, Mazin W, Sparsø T, Vangkilde A, et al. Risk of psychiatric disorders among individuals with the 22q11.2 deletion or duplication: a danish nationwide, register-based study. JAMA Psychiatry. 2017;74:282–90.
-
Marshall CR, Howrigan DP, Merico D, Thiruvahindrapuram B, Wu W, Greer DS, et al. Contribution of copy number variants to schizophrenia from a genome-wide study of 41,321 subjects. Nat Genet. 2017;49:27–35.
-
Emanuel BS, McDonald-McGinn D, Saitta SC, Zackai EH. CHAPTER 2 – The 22q11.2 deletion syndrome. Adv Pediatr. 2001 ;48:39–73.
-
Boot E, Óskarsdóttir S, Loo JCY, Crowley TB, Orchanian-Cheff A, Andrade DM, et al. Updated clinical practice recommendations for managing adults with 22q11.2 deletion syndrome. Genet Med [Internet]. 2023;25. https://www.gimjournal.org/article/S1098-3600(22)01028-0/fulltext.
-
McDonald-McGinn DM, Sullivan KE, Marino B, Philip N, Swillen A, Vorstman JAS, et al. 22q11.2 deletion syndrome. Nat Rev Dis Primer. 2015;1:15071.
-
Rees E, Kirov G, Sanders A, Walters JTR, Chambert KD, Shi J, et al. Evidence that duplications of 22q11.2 protect against schizophrenia. Mol Psychiatry. 2014;19:37–40.
-
Owen MJ, Legge SE, Rees E, Walters JTR, O’Donovan MC. Genomic findings in schizophrenia and their implications. Mol Psychiatry. 2023;28:3638–47.
-
Ou Z, Berg JS, Yonath H, Enciso VB, Miller DT, Picker J, et al. Microduplications of 22q11.2 are frequently inherited and are associated with variable phenotypes. Genet Med. 2008;10:267–77.
-
Portnoï MF. Microduplication 22q11.2: A new chromosomal syndrome. Eur J Med Genet. 2009;52:88–93.
-
Lin A, Vajdi A, Kushan-Wells L, Helleman G, Hansen LP, Jonas RK, et al. Reciprocal copy number variations at 22q11.2 produce distinct and convergent neurobehavioral impairments relevant for schizophrenia and autism spectrum disorder. Biol Psychiatry. 2020;88:260–72.
-
Villalón-Reina JE, Martínez K, Qu X, Ching CRK, Nir TM, Kothapalli D, et al. Altered white matter microstructure in 22q11.2 deletion syndrome: a multisite diffusion tensor imaging study. Mol Psychiatry. 2020;25:2818–31.
-
Sun D, Ching CRK, Lin A, Forsyth JK, Kushan L, Vajdi A, et al. Large-scale mapping of cortical alterations in 22q11.2 deletion syndrome: Convergence with idiopathic psychosis and effects of deletion size. Mol Psychiatry. 2020;25:1822–34.
-
Ching CRK, Gutman BA, Sun D, Villalon Reina J, Ragothaman A, Isaev D, et al. Mapping subcortical brain alterations in 22q11.2 deletion syndrome: effects of deletion size and convergence with idiopathic neuropsychiatric illness. Am J Psychiatry. 2020;177:589–600.
-
Schleifer CH, O’Hora KP, Fung H, Xu J, Robinson TA, Wu AS, et al. Effects of gene dosage and development on subcortical nuclei volumes in individuals with 22q11.2 copy number variations. Neuropsychopharmacology. 2024;49:1024–32.
-
Mori S, Zhang J. Principles of diffusion tensor imaging and its applications to basic neuroscience research. Neuron. 2006;51:527–39.
-
Pierpaoli C, Jezzard P, Basser PJ, Barnett A, Di Chiro G. Diffusion tensor MR imaging of the human brain. Radiology. 1996;201:637–48.
-
Beaulieu C. The basis of anisotropic water diffusion in the nervous system – a technical review. NMR Biomed. 2002;15:435–55.
-
Raven EP, Veraart J, Kievit RA, Genc S, Ward IL, Hall J, et al. In vivo evidence of microstructural hypo-connectivity of brain white matter in 22q11.2 deletion syndrome. Mol Psychiatry. 2023;28:4342–52.
-
Shin D, Kim CN, Ross J, Hennick KM, Wu SR, Paranjape N, et al. Thalamocortical organoids enable in vitro modeling of 22q11.2 microdeletion associated with neuropsychiatric disorders. Cell Stem Cell. 2024;31:421–.e8.
-
Wu P, Teot L, Murdoch G, Monaghan-Nichols AP, McFadden K. Neuropathology of 22q11 deletion syndrome in an infant. Pediatr Dev Pathol. 2014;17:386–92.
-
Seitz-Holland J, Lyons M, Kushan L, Lin A, Villalon-Reina JE, Cho KIK, et al. Opposing white matter microstructure abnormalities in 22q11.2 deletion and duplication carriers. Transl Psychiatry. 2021;11:1–11.
-
Mukai J, Tamura M, Fénelon K, Rosen AM, Spellman TJ, Kang R, et al. Molecular substrates of altered axonal growth and brain connectivity in a mouse model of schizophrenia. Neuron. 2015;86:680–95.
-
Mukai J, Dhilla A, Drew LJ, Stark KL, Cao L, MacDermott AB, et al. Palmitoylation-dependent neurodevelopmental deficits in a mouse model of 22q11 microdeletion. Nat Neurosci. 2008;11:1302–10.
-
Hiramoto T, Sumiyoshi A, Yamauchi T, Tanigaki K, Shi Q, Kang G, et al. Tbx1, a gene encoded in 22q11.2 copy number variant, is a link between alterations in fimbria myelination and cognitive speed in mice. Mol Psychiatry. 2022;27:929–38.
-
Budel S, Padukkavidana T, Liu BP, Feng Z, Hu F, Johnson S, et al. Genetic variants of Nogo-66 receptor with possible association to schizophrenia block myelin inhibition of axon growth. J Neurosci. 2008;28:13161–72.
-
Fernandez A, Meechan DW, Karpinski BA, Paronett EM, Bryan CA, Rutz HL, et al. Mitochondrial dysfunction leads to cortical under-connectivity and cognitive impairment. Neuron. 2019;102:1127–.e3.
-
Khan TA, Revah O, Gordon A, Yoon SJ, Krawisz AK, Goold C, et al. Neuronal defects in a human cellular model of 22q11.2 deletion syndrome. Nat Med. 2020;26:1888–98.
-
Stark KL, Xu B, Bagchi A, Lai WS, Liu H, Hsu R, et al. Altered brain microRNA biogenesis contributes to phenotypic deficits in a 22q11-deletion mouse model. Nat Genet. 2008;40:751–60.
-
Lin A, Forsyth JK, Hoftman GD, Kushan-Wells L, Jalbrzikowski M, Dokuru D, et al. Transcriptomic profiling of whole blood in 22q11.2 reciprocal copy number variants reveals that cell proportion highly impacts gene expression. Brain Behav Immun – Health. 2021;18:100386.
-
Nehme R, Pietiläinen O, Artomov M, Tegtmeyer M, Valakh V, Lehtonen L, et al. The 22q11.2 region regulates presynaptic gene-products linked to schizophrenia. Nat Commun. 2022;13:3690.
-
Huang H, Zhang J, Wakana S, Zhang W, Ren T, Richards LJ, et al. White and gray matter development in human fetal, newborn and pediatric brains. Neuroimage. 2006;33:27–38.
-
Kochunov P, Williamson DE, Lancaster J, Fox P, Cornell J, Blangero J, et al. Fractional anisotropy of water diffusion in cerebral white matter across the lifespan. Neurobiol Aging. 2012;33:9–20.
-
Wilson S, Pietsch M, Cordero-Grande L, Price AN, Hutter J, Xiao J, et al. Development of human white matter pathways in utero over the second and third trimester. Proc Natl Acad Sci USA. 2021;118:e2023598118.
-
Lebel C, Beaulieu C. Longitudinal development of human brain wiring continues from childhood into adulthood. J Neurosci. 2011;31:10937–47.
-
Huang H, Xue R, Zhang J, Ren T, Richards LJ, Yarowsky P, et al. Anatomical characterization of human fetal brain development with diffusion tensor magnetic resonance imaging. J Neurosci. 2009;29:4263–73.
-
Vasung L, Raguz M, Kostovic I, Takahashi E. Spatiotemporal relationship of brain pathways during human fetal development using high-angular resolution diffusion MR imaging and histology. Front Neurosci [Internet]. 2017;11. Available from: https://www.frontiersin.org/journals/neuroscience/articles/10.3389/fnins.2017.00348/full.
-
Horgos B, Mecea M, Boer A, Szabo B, Buruiana A, Stamatian F, et al. White matter dissection of the fetal brain. Front Neuroanat. 2020;14:584266.
-
Rakic P, Yakovlev PI. Development of the corpus callosum and cavum septi in man. J Comp Neurol. 1968;132:45–72.
-
Ren T, Anderson A, Shen WB, Huang H, Plachez C, Zhang J, et al. Imaging, anatomical, and molecular analysis of callosal formation in the developing human fetal brain. Anat Rec A Discov Mol Cell Evol Biol. 2006;288A:191–204.
-
Ouyang M, Dubois J, Yu Q, Mukherjee P, Huang H. Delineation of early brain development from fetuses to infants with diffusion MRI and beyond. Neuroimage. 2019;185:836–50.
-
Takahashi E, Folkerth RD, Galaburda AM, Grant PE. Emerging cerebral connectivity in the human fetal brain: an MR tractography study. Cereb Cortex. 2012;22:455–64.
-
Kostović I, Jovanov-Milošević N. The development of cerebral connections during the first 20–45 weeks’ gestation. Semin Fetal Neonatal Med. 2006;11:415–22.
-
LaMantia AS, Rakic P. Axon overproduction and elimination in the corpus callosum of the developing rhesus monkey. J Neurosci. 1990;10:2156–75.
-
LaMantia AS, Rakic P. Axon overproduction and elimination in the anterior commissure of the developing rhesus monkey. J Comp Neurol. 1994;340:328–36.
-
Paus T. Growth of white matter in the adolescent brain: Myelin or axon? Brain Cogn. 2010;72:26–35.
-
Benes FM. Myelination of cortical-hippocampal relays during late adolescence. Schizophr Bull. 1989;15:585–93.
-
Kulikova S, Hertz-Pannier L, Dehaene-Lambertz G, Buzmakov A, Poupon C, Dubois J. Multi-parametric evaluation of the white matter maturation. Brain Struct Funct. 2015;220:3657–72.
-
Jones DK, Knösche TR, Turner R. White matter integrity, fiber count, and other fallacies: The do’s and don’ts of diffusion MRI. Neuroimage. 2013;73:239–54.
-
Zhang H, Schneider T, Wheeler-Kingshott CA, Alexander DC. NODDI: Practical in vivo neurite orientation dispersion and density imaging of the human brain. Neuroimage. 2012;61:1000–16.
-
Lee HH, Yaros K, Veraart J, Pathan JL, Liang FX, Kim SG, et al. Along-axon diameter variation and axonal orientation dispersion revealed with 3D electron microscopy: implications for quantifying brain white matter microstructure with histology and diffusion MRI. Brain Struct Funct. 2019;224:1469–88.
-
Grussu F, Schneider T, Tur C, Yates RL, Tachrount M, Ianuş A, et al. Neurite dispersion: a new marker of multiple sclerosis spinal cord pathology? Ann Clin Transl Neurol. 2017;4:663–79.
-
Jespersen SN, Bjarkam CR, Nyengaard JR, Chakravarty MM, Hansen B, Vosegaard T, et al. Neurite density from magnetic resonance diffusion measurements at ultrahigh field: Comparison with light microscopy and electron microscopy. Neuroimage. 2010;49:205–16.
-
Sato K, Kerever A, Kamagata K, Tsuruta K, Irie R, Tagawa K, et al. Understanding microstructure of the brain by comparison of neurite orientation dispersion and density imaging (NODDI) with transparent mouse brain. Acta Radiol Open. 2017;6:2058460117703816.
-
Kraguljac NV, Guerreri M, Strickland MJ, Zhang H. Neurite orientation dispersion and density imaging in psychiatric disorders: a systematic literature review and a technical note. Biol Psychiatry Glob Open Sci. 2022;3:10–21.
-
Genc S, Malpas CB, Holland SK, Beare R, Silk TJ. Neurite density index is sensitive to age related differences in the developing brain. Neuroimage. 2017;148:373–80.
-
Genc S, Raven EP, Drakesmith M, Blakemore SJ, Jones DK. Novel insights into axon diameter and myelin content in late childhood and adolescence. Cereb Cortex N Y NY. 2023;33:6435–48.
-
Fan Q, Tian Q, Ohringer NA, Nummenmaa A, Witzel T, Tobyne SM, et al. Age-related alterations in axonal microstructure in the corpus callosum measured by high-gradient diffusion MRI. Neuroimage. 2019;191:325–36.
-
Benes FM, Turtle M, Khan Y, Farol P. Myelination of a key relay zone in the hippocampal formation occurs in the human brain during childhood, adolescence, and adulthood. Arch Gen Psychiatry. 1994;51:477–84.
-
Somerville LH, Bookheimer SY, Buckner RL, Burgess GC, Curtiss SW, Dapretto M, et al. The lifespan human connectome project in development: a large-scale study of brain connectivity development in 5–21 year olds. Neuroimage. 2018;183:456.
-
Bookheimer SY, Salat DH, Terpstra M, Ances BM, Barch DM, Buckner RL, et al. The lifespan human connectome project in aging: an overview. Neuroimage. 2019;185:335–48.
-
Harms MP, Somerville LH, Ances BM, Andersson J, Barch DM, Bastiani M, et al. Extending the human connectome project across ages: imaging protocols for the lifespan development and aging projects. Neuroimage. 2018;183:972–84.
-
Garyfallidis E, Brett M, Amirbekian B, Rokem A, Van Der Walt S, Descoteaux M, et al. Dipy, a library for the analysis of diffusion MRI data. Front Neuroinformatics [Internet]. 2014;8. https://www.frontiersin.org/articles/10.3389/fninf.2014.00008.
-
Manjón JV, Coupé P, Concha L, Buades A, Collins DL, Robles M. Diffusion weighted image denoising using overcomplete local PCA. PLoS ONE. 2013;8:e73021.
-
Andersson JLR, Skare S, Ashburner J. How to correct susceptibility distortions in spin-echo echo-planar images: application to diffusion tensor imaging. Neuroimage. 2003;20:870–88.
-
Andersson JLR, Sotiropoulos SN. An integrated approach to correction for off-resonance effects and subject movement in diffusion MR imaging. Neuroimage. 2016;125:1063–78.
-
Jahanshad N, Kochunov P, Sprooten E, Mandl RC, Nichols TE, Almassy L, et al. Multi-site genetic analysis of diffusion images and voxelwise heritability analysis: A pilot project of the ENIGMA–DTI working group. Neuroimage. 2013;81:455–69.
-
Mori S, Oishi K, Jiang H, Jiang L, Li X, Akhter K, et al. Stereotaxic white matter atlas based on diffusion tensor imaging in an ICBM template. Neuroimage. 2008;40:570–82.
-
Beer JC, Tustison NJ, Cook PA, Davatzikos C, Sheline YI, Shinohara RT, et al. Longitudinal ComBat: A method for harmonizing longitudinal multi-scanner imaging data. Neuroimage. 2020;220:117129.
-
Pakkenberg B, Pelvig D, Marner L, Bundgaard MJ, Gundersen HJG, Nyengaard JR, et al. Aging and the human neocortex. Exp Gerontol. 2003;38:95–9.
-
Wood S. Generalized additive models: an introduction with R. Chapman and Hall/CRC; 2017.
-
Shashi V, Francis A, Hooper SR, Kranz PG, Zapadka M, Schoch K, et al. Increased corpus callosum volume in children with chromosome 22q11.2 deletion syndrome is associated with neurocognitive deficits and genetic polymorphisms. Eur J Hum Genet. 2012;20:1051–7.
-
Garic D, Yeh FC, Graziano P, Dick AS. In vivo restricted diffusion imaging (RDI) is sensitive to differences in axonal density in typical children and adults. Brain Struct Funct. 2021;226:2689–705.
-
Aboitiz F, Scheibel AB, Fisher RS, Zaidel E. Fiber composition of the human corpus callosum. Brain Res. 1992;598:143–53.
-
Liewald D, Miller R, Logothetis N, Wagner HJ, Schüz A. Distribution of axon diameters in cortical white matter: an electron-microscopic study on three human brains and a macaque. Biol Cybern. 2014;108:541–57.
-
Meechan DW, Tucker ES, Maynard TM, LaMantia AS. Diminished dosage of 22q11 genes disrupts neurogenesis and cortical development in a mouse model of 22q11 deletion/DiGeorge syndrome. Proc Natl Acad Sci. 2009;106:16434–45.
-
Lin A, Ching CRK, Vajdi A, Sun D, Jonas RK, Jalbrzikowski M, et al. Mapping 22q11.2 gene dosage effects on brain morphometry. J Neurosci. 2017;37:6183–99.
-
Hegarty JP, Monterrey JC, Tian Q, Cleveland SC, Gong X, Phillips JM, et al. A twin study of altered white matter heritability in youth with autism spectrum disorder. J Am Acad Child Adolesc Psychiatry. 2024;63:65–79.
-
Saaybi S, AlArab N, Hannoun S, Saade M, Tutunji R, Zeeni C, et al. Pre- and post-therapy assessment of clinical outcomes and white matter integrity in autism spectrum disorder: pilot study. Front Neurol. 2019;10:877.
-
Singh T, Poterba T, Curtis D, Akil H, Al Eissa M, Barchas JD, et al. Rare coding variants in ten genes confer substantial risk for schizophrenia. Nature. 2022;604:509–16.
-
Schleifer C, Lin A, Kushan L, Ji JL, Yang G, Bearden CE, et al. Dissociable disruptions in thalamic and hippocampal resting-state functional connectivity in youth with 22q11.2 deletions. J Neurosci. 2019;39:1301–19.
-
Mancini V, Zöller D, Schneider M, Schaer M, Eliez S. Abnormal development and dysconnectivity of distinct thalamic nuclei in patients with 22q11.2 deletion syndrome experiencing auditory hallucinations. Biol Psychiatry Cogn Neurosci Neuroimaging. 2020;5:875–90.
-
Moreau CA, Kumar K, Harvey A, Huguet G, Urchs SGW, Schultz LM, et al. Brain functional connectivity mirrors genetic pleiotropy in psychiatric conditions. Brain. 2022;146:1686–96.
-
Steullet P. Thalamus-related anomalies as candidate mechanism-based biomarkers for psychosis. Schizophr Res. 2020;226:147–57.
-
Alvino FG, Gini S, Minetti A, Pagani M, Sastre-Yagüe D, Barsotti N, et al. Synaptic-dependent developmental dysconnectivity in 22q11.2 deletion syndrome. Sci Adv. 2025;11:eadq2807.
-
Rao SB, Brundu F, Chen Y, Sun Y, Zhu H, Shprintzen RJ, et al. Aberrant pace of cortical neuron development in brain organoids from patients with 22q11.2 deletion syndrome and schizophrenia [Internet]. bioRxiv; 2023 https://www.biorxiv.org/content/10.1101/2023.10.04.557612v1.
-
Toyoshima M, Akamatsu W, Okada Y, Ohnishi T, Balan S, Hisano Y, et al. Analysis of induced pluripotent stem cells carrying 22q11.2 deletion. Transl Psychiatry. 2016;6:e934–e934.
-
Jalbrzikowski M, Lin A, Vajdi A, Grigoryan V, Kushan L, Ching CRK, et al. Longitudinal trajectories of cortical development in 22q11.2 copy number variants and typically developing controls. Mol Psychiatry. 2022;27:4181–90.
-
Jelescu IO, Budde MD. Design and validation of diffusion MRI models of white Matter. Front Phys [Internet]. 2017: https://www.frontiersin.org/journals/physics/articles/10.3389/fphy.2017.00061/full.
-
Guerrero JM, Adluru N, Bendlin BB, Goldsmith HH, Schaefer SM, Davidson RJ, et al. Optimizing the intrinsic parallel diffusivity in NODDI: An extensive empirical evaluation. PLoS ONE. 2019;14:e0217118.
-
Jelescu IO, Veraart J, Fieremans E, Novikov DS. Degeneracy in model parameter estimation for multi-compartmental diffusion in neuronal tissue. NMR Biomed. 2016;29:33–47.
-
Lampinen B, Szczepankiewicz F, Mårtensson J, van Westen D, Sundgren PC, Nilsson M. Neurite density imaging versus imaging of microscopic anisotropy in diffusion MRI: A model comparison using spherical tensor encoding. Neuroimage. 2017;147:517–31.
-
Kamiya K, Hori M, Aoki S. NODDI in clinical research. J Neurosci Methods. 2020;346:108908.
Acknowledgements
We thank Sergiu Pasca, M.D., for contributing to the study. This work used computational and storage services associated with the Hoffman2 Cluster which is operated by the UCLA Office of Advanced Research Computing’s Research Technology Group and servers operated by the USC Mark and Mary Stevens Neuroimaging and Informatics Institute (INI). Data from the UCLA and Stanford cohorts were supported by grants from the MCHRI Uytengsu-Hamilton 22q11 Neuropsychiatry Research Award (UH22QEXTFY22-04) and NIH 1R21MH116473-01A1, R37MH085953, 1U01MH119736-01 (to CEB), R01MH100900 (to JH), and R01MH129858 (to CEB, PMT, IES, OAA). JEVR acknowledges funding from NIH (T32AG058507). IES and RB acknowledge funding from South-Eastern Norway Regional Health Authority (#2020060). Research reported in this publication was supported by the National Institute Of Mental Health of the National Institutes of Health under Award Number U01MH109589/ U01AG052564 and by funds provided by the McDonnell Center for Systems Neuroscience at Washington University in St. Louis. The HCP-Development 2.0 Release data used in this report came from DOI: 10.15154/1520708. The HCP-Aging 2.0 Release data used in this report came from DOI: 10.15154/1520707.
Ethics declarations
Competing interests
OAA is a consultant for Cortechs.ai and Precision Health and has received speaker’s honoraria from Lundbeck, Janssen, Otsuka, Lilly and Sunovion. The remaining authors have no conflicts of interest to declare.
Ethics approval and consent to participate
All participants gave verbal and written informed consent to participate in the study, after study procedures were fully explained. Participants under 18 years of age provided written assent and their parent/guardian provided written consent. All study procedures and documents conform to the standards of the Declaration of Helsinki and were approved by the University of California Los Angeles (UCLA) (IRB#10-00107) and Stanford University Institutional Review Boards (Protocol ID:27795, IRB Number: 350).
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Boen, R., Villalón-Reina, J.E., Kushan, L. et al. Gene dosage effects of 22q11.2 copy number variants on in-vivo measures of white matter axonal density and dispersion. Mol Psychiatry (2026). https://doi.org/10.1038/s41380-026-03489-4
-
Received:
-
Revised:
-
Accepted:
-
Published:
-
Version of record:
-
DOI: https://doi.org/10.1038/s41380-026-03489-4




