Developmental trajectory of cortical and subcortical structural asymmetry in youths
This study used ABCD datasets, comprising over 11,000 young participants aged from 10 to 14 years (SM Table 1). Asymmetry indices (AI), commonly defined as 2*(L-R)/(L + R), were utilized to evaluate the asymmetry in the cerebral structures (Kong et al., 2018), where “L” and “R” signifies the structural attributes of the left and right hemisphere, respectively. To ensure comparability of asymmetry indices across different age, we normalized the asymmetry measures at all three time points using the sum of left and right hemispheric attributes (L + R) at baseline (age 10).
Figure 2A displays the group-averaged asymmetry patterns in 3 ages. The asymmetry patterns of cortical surface area (CSA), cortical thickness (CT) and subcortical volume (Vol) are relatively stable, converge to those observed in the large adult sample (Kong et al., 2018). CSA mainly showed left-lateralization in sensorimotor, auditory area, anterior and posterior cingulate, and right lateralization in frontotemporal and inferior parietal cortex. CT mainly showed left-lateralization in anterior and dorsal prefrontal cortex and para-hippocampal gyrus, and right-lateralization in occipito-temporal areas. For subcortical structures, leftward lateralization in volume is observed in the lateral ventricle (LatVent), pallidum, ventral diencephalon (DC), and thalamus, while rightward lateralization occurs in the hippocampus, amygdala, caudate, pallidum, and inferior lateral ventricle (InfLatVent).
To explore the developmental trend of structural asymmetry, we utilized linear mixed models to delineate the effects of age on structural asymmetry (using R package lme4), where gender, pubertal stage, race, handedness and intracranial volume (ITV) were treated as fixed effects while family ID nested in image acquisition sites and subject IDs were separately accounted for as two random effects. Multiple comparisons were conducted using the False Discovery Rate (FDR) method (34 * 2 cortical asymmetry measures and 12 subcortical asymmetry measures, totaling 80 comparisons). We found that as age increases, motor area, medial prefrontal and temporo-occipital regions showed a left-lateralized trend in CSA, while the parietal lobule and frontal triangular gyrus showed a right-lateralized trend. For CT, the frontal and medial temporal lobes become more left-lateralized, while the parieto-occipital and lateral orbitofrontal lobes become more right-lateralized with age. For subcortex, thalamus, putamen and caudate exhibit increased left-lateralization with age, whereas pallidum tends to become more right-lateralized with age (Fig. 2B). Overall, the most pronounced developmental effect in asymmetry are observed in subcortical region volume (the smallest p = 1e-22, t = 9.77, thalamus), followed by CT (middle orbital frontal cortex, t = 5.93, p = 3e-9; inferior parietal lobe, t=-5.88, p = 4e-9) whereas CSA demonstrates relatively small age effects (inferior parietal lobe, t = 4.11, p = 3.9e-5) (Fig. 2B and SM Table 2).
We further investigated the correlation between regional CSA/CT/subcortical volume with age to understand how structural changes with age may affect asymmetry development (Fig. 2C). We noted a general decreasing trend in regional CSA/CT from 10 to 14 years, which suggest the presence of neural pruning (Juraska & Drzewiecki, 2020; Spear, 2013). This result indicates that the development of cortical asymmetry in pre-adolescence may be primarily driven by more pruning on one hemisphere compared to the other.
Relationship between asymmetry, cognition/psychiatric phenotypes and early factors
Linear mixed models were utilized to assess the relationships between baseline structural asymmetry and cognitive functions cognition [fluid, crystal and total intelligence from NIH Toolbox Cognitive Battery, Average reaction time of Little Man Task (LMT, measuring visual-spatial processing), and WISC-V Matrix Reasoning Total Raw Score (WISCV TRS)] and psychiatric symptoms [internalization, externalization, and total problem scores from Child Behavior Checklist (CBCL), and common mental disorder diagnoses from parent-reported KSADS] (B. Casey et al., 2018), with FDR correction. Baseline data were selected to maximize sample size (N = 11,362). Alongside common covariates described above, the models include additional fixed-effect covariates of age, while grouping variables (nested image acquisition sites and family) were treated as random effects in ABCD association analyses.
We observed a series of associations between asymmetry and behavioral phenotypes (Fig. 3A-B, SM Table 3). Specifically, for cortical structures, cognitive abilities (particularly crystallized intelligence) was negatively associated with asymmetry of lateral orbital frontal cortex (LOrbF) (CSA, NIH fluid intelligence: t=-4.68, p = 3e-6; NIH crystallized intelligence, t=-5.27, p = 1e-7; NIH total intelligence: t=-6.03, p = 2e-9) and superior frontal cortex (SF) (CT, NIH fluid: t=-3.67, p = 2e-4; NIH crystal: t=-5.1, p = 3e-7; NIH total: t=-5.13, p = 3e-7), and was positively related with para-hippocampus (PaH) (CT, NIH fluid: t = 3.17, p = 2e-3; NIH crystal: t = 5.31, p = 1e-7; NIH total: t = 4.51, p = 7e-6; WISCV TRS, t = 3.61, p = 3e-4, see Fig. 4B). The CT asymmetry of insula, fusiform gyrus and isthmic cingulate also positively correlated with intelligence. For subcortical structures, asymmetry of thalamus correlated positively with cognitive abilities including NIH total intelligence (volume, t = 2.57, p = 0.01) and fluid intelligence (t = 3.41, p = 0.0007). Additionally, CBCL total problems were positively associated with asymmetry of LOrbF (CSA, t = 3.57, p = 0.0004) and middle cingulate (CT, t = 3.02, p = 0.0026), and was negatively associated with asymmetry in the inferior lateral ventricles (volume, t=-3.25, p = 0.0012) and thalamus (volume, t=-2.84, p = 0.0047). For psychiatric diagnosis, only decreased ventral diencephalon volume was observed in post-traumatic stress disorder (PTSD, t=-3.47, p = 5e-4).
To understand the function of regions whose asymmetry was associated with behaviroal phenotypes, we calculated the correlation between t-maps of crystallized intelligence (with CSA/CT asymmetry, Fig. 3B) and meta-analysis maps from Neurosynth (https://www.neurosynth.org/) corresponding to 127 mental terms (Hansen et al., 2022). The t-map of CSA asymmetry-cognition association was positively correlated with meta-maps of memory and language (speech perception) while negatively correlated with decision-making (uncertainty, reward expectation). The t-map of CT asymmetry-cognition association was positively correlated with visual processing, imagination, and consciousness, and negatively with memory (autobiographical, retrieval), language, and sociability (social cognition, semantic memory, context), see Fig. 3C.
We furthermore investigated whether regions in left or right hemisphere are responsible for the observed asymmetry pattern that was related to cognition. We focused on 7 regions whose asymmetry was most strongly correlated with crystallized intelligence, using a linear mixed model (FDR q < 0.001) to assess group (top 30% vs. bottom 30% in crystallized intelligence), hemisphere (left vs. right), and interaction effects (Supplemental Method). In high-intelligence group, changed asymmetry occurred alongside increased CSA, CT, and volume, compared to the low-intelligence group. For instance, bilateral thalamic volume in the high-intelligence group was larger than that in the low-intelligence group (t = 22.2, p = 1.2e-105), with the right thalamus being larger than the left (t = 105, p < 0.0001). The interaction effect was also significant (t = 5.51, p = 3.72e-8), indicating that the left > right asymmetry in the high-intelligence group was increased compared to the low-intelligence group (Fig. 4 and SM Table 9).
To explore the impact of environmental factors on brain asymmetry, three latent factors defined by Group Factor Analysis (GFA) based on 22 baseline variables from ABCD 2.0.1 (Gonzalez et al., 2019) were used for association analysis: general socioeconomic status (SES) factor (parent income/education, housing security and neighborhood safety), social factor (adolescents' perceived social support from family and school) and perinatal health factor (positively associated with birth weight/gestational age, negatively with maternal age/prenatal medical conditions). We found that perinatal factors were negatively associated with laterality of SF (CT, t=-6.33, p = 3e-10), rostral middle frontal gyrus and anterior/isthmus cingulate, and with asymmetry of LOrbF (CSA, t=-6.59, p = 5e-11) and transverse/inferior temporal cortex. Conversely, perinatal factors were positively associated with asymmetry of inferior parietal cortex (CT) and thalamus, amygdala and ventral diencephalon (volume). SES factors show a positive correlation with asymmetry of thalamus (volume, t = 3.83, p = 0.0001), while social factors did not exhibit significant correlation (Fig. 5A).
Perinatal factors shape cognitive abilities via corpus callosum-mediated structural asymmetry
To further investigate the neurobiological (macrostructural) and environmental factors shaping asymmetry and whether such factors impact cognitive phenotypes, we conducted a series of pathway analyses using longitudinal data. All models were build using the lavaan package of R, and all covariates in association analysis using linear mixed models were regressed out from CC and AIs before modeling.
Considering the influence of the corpus callosum (CC) on brain asymmetry (Hinkley et al., 2016; Karolis, Corbetta, & Thiebaut de Schotten, 2019; Tzourio-Mazoyer, 2016), we first examined the relationship between CC integrity and brain asymmetry. We used the first principal component (CC PC1) from 4 diffusion imaging indices which explains 50.17% variance, correlates positively with microstructural integrity, and increases with age (t = 38.2, p = 1.66e-300) (Supplemental Method, SM Table 2&4). Cross-lagged panel models (CLPMs) revealed that CC integrity at age 10 significantly influenced structural asymmetry at age 12. Specifically, CC PC1 correlated positively with asymmetry of orbitofrontal cortex (CSA, z = 7.61, p = 2.7e-14), frontal cortex (CT, especially superior frontal gyrus, z = 8.45, p < 0.0001), anterior cingulate (CT), superior temporal cortex (CT), pallidum (volume, z = 10.9, p < 0.0001), and amygdala (volume, z = 11.4, p < 0.0001), and correlated negatively with asymmetry of insula (CT, z=-9.62, p < 0.0001), parietal and temporal cortex (CT), and thalamus (volume), see Fig. 5A. Interestingly, this pattern was inversely related to that between asymmetry and crystallized intelligence (Pearson r=-0.47, p = 1e-5). Asymmetry at age 10 did not significantly affect CC at age 12, indicating a unidirectional effect of CC on asymmetry during pre-adolescence (SM Tables 5–6).
We then examined the cascading relationship between four factors by longitudinal structural equation models (SEMs): (1) perinatal factors, which correlated most with asymmetry among early factors, (2) CC PC1 at age 10, (3) CT asymmetry at age 12, and (4) crystallized intelligence at age 12. Our findings show that better perinatal health is associated with lower CC PC1 levels at age 10 (β = 0.02, p = 0.04), which correlates positively with asymmetry in superior frontal gyrus (CT, β = 3e-3, p < 0.0001) and pallidum (volume, β = 0.01, p < 0.0001) and negatively correlates with asymmetry of para-hippocampus (CT, β = 3e-4, p = 1e-6) at age 12, ultimately leading to enhanced crystallized intelligence at age 12 (Fig. 5B and SM Table 7–8).
Biological processes that modulate asymmetry development
To elucidate the genetic underpinnings of structural brain asymmetry both at baseline and throughout its longitudinal development in youth, we conducted Genome-Wide Association Study (GWAS) using whole-genome sequencing data from the ABCD study. We focus on baseline (bl) structural asymmetry, denoted as AIbl (10 years old), and the rate of asymmetry change with age, denoted as ΔAI [(12-year follow-up − 10-year baseline)/age gap]. Our analysis revealed that ΔAI exhibited a more extensive correlation with single nucleotide polymorphisms (SNPs) (p < 1e-8) across the brain compared to AIbl (Fig. 6A-B). Subsequently, all significant SNPs were mapped to gene loci based on the Genome Reference Consortium Human Build 37 (GRCh37) database, enabling the identification of gene sets significantly associated with AIbl and ΔAI. Enrichment analysis based on Functional Mapping and Annotation (FUMA) platform (https://fuma.ctglab.nl/) (Watanabe, Taskesen, Van Bochoven, & Posthuma, 2017) show that both baseline and longitudinal structural asymmetry are related to synapse-related processes, neuron projection, neurogenesis, and membrane transport (Fig. 6C).