Abstract
Introduction The cerebellum, most known for its role in motor control, exerts a key role in cognition. Multiple lines of evidence across human functional, lesion and animal data point to a role of the cerebellum, in particular of Crus I, Crus II and Lobule VIIB, in cognitive function. However, whether cerebellar substrates pertaining to distinct facets of cognitive function exist is not known.
Methods We analyzed structural neuroimaging data from the Healthy Brain Network (HBN). Cerebellar parcellation was performed via a standard validated automated segmentation pipeline (CERES) with stringent visual quality check (n = 662 subjects retained from initial n = 1452). We used data-driven canonical correlation analyses (CCA) to examine regional gray matter volumetric (GMV) differences in association to cognitive function assessed with the NIH Toolbox Cognition Domain (NIH-TB). Our multivariate analyses accounted for psychopathology severity, age, sex, scan location and intracranial volume.
Results Multivariate CCA uncovered a significant correlation between two components entailing a latent cognitive canonical variate composed of NIH-TB subscales and the brain canonical variate (cerebellar regions’ GMV and intracranial volume, ICV). A bootstrapping and a permutation procedure ensured the results are statistically significant and the CCA model, stable. The identified components correspond to only partly shared cerebellar -cognitive function relationship with a first map encompassing cognitive flexibility (r=0.89) and speed of processing (r=0.65) associated with regional gray matter volume in Crus II (r=0.57) and Lobule X (r=0.59) and a second map including the Crus I (r=0.49) and Lobule VI (r=0.49) associated with cognitive control (r=-0.51). Working memory associations were similarly present in both these maps (Crus II, Lobule X, Crus I and Lobule VI) for the first (r=0.52) and second (r=0.51) component.
Discussion Our results show evidence in favor of structural sub-specialization in the cerebellum, independently of psychopathology contributions to cognitive function and brain structure. Overall, these findings highlight a prominent role for the human cerebellum in cognitive function for flexible and stable adaptive behavior.
Introduction
The cerebellum is a fascinating infratentorial brain structure with a pivotal role in human cognition (Schmahmann et al. 2019). Research on the cerebellum has been traditionally limited to its role in motor control, even though the majority of the cerebellar cortex is not involved in motor action planning or execution (King et al. 2019). Cerebellar lesions across different diagnostic entities are associated with a diverse palette of cognitive deficits including disturbances of executive function such as planning, set-shifting, working memory and verbal fluency (Schmahmann and Sherman 1998).
A recent study linked cerebellar anatomy to cognitive functioning and found that anatomical features predicted both general cognitive function and psychopathology (Moberget et al. 2019). However, regional cerebellar morphometry differences relating to general cognitive function in psychosis (Moussa-Tooks et al. 2022) or in autism (Laidi et al. 2022) were no different from controls. One reason for the discrepancy could be the difference in these studies’ approaches (dimensional vs. case-control comparisons). These studies looked into general cognitive function which may fail to capture anatomy-cognition links in the cerebellum because the relationships possibly concern specific cognitive domains. Studies are starting to progressively elucidate the functional organization of the cerebellum (King et al. 2019; Buckner et al. 2011; Guell et al. 2018). Yet, a finer-grained nuanced investigation of the distinct facets of cognition is currently lacking leaving open the question of whether a structural cerebellar subspecialization exists with respect to cognitive abilities. Previous evidence showed that brain lesions in Crus I and Crus II, VIIB (Stoodley et al. 2016; Stoodley and Schmahmann 2010; Kansal et al. 2017) and to a lower extent VIIIA and VI are associated with executive function performance. A seminal study in patients with cerebellar degeneration showed that distinct components of cognitive function (e.g. executive function, working memory, perceptual processing and so on) relate differently to cerebellar topography (Kansal et al. 2017). Nevertheless, even though lesion studies are informative they have limitations, and a large-scale examination of the cerebellar mapping of distinct cognitive components is lacking. Furthermore, it is worth noting that psychopathology severity (e.g., levels of anxiety, depression and so on) has not been systematically accounted for in the reviewed studies examining associations with cognitive function. Psychopathology severity impacts cognitive function (Chavez-Baldini et al. 2021) but also brain structure properties particularly in the developing brain (Patel et al. 2022; Mattoni, Wilson, and Olino 2021; Romer, Ren, and Pizzagalli 2023). Recently, cerebellar structure has been shown to be linked to both general cognitive function and psychopathology (Moberget et al. 2019). However, how cognitive mapping in cerebellar anatomy is represented independently of psychopathology contributions is not fully elucidated. Such investigation is of key interest with respect to the ensuing potential for clinical (e.g. neuromodulation) translation (Yao et al. 2022). In sum, a major gap remains in the current understanding of cerebellar contributions to cognitive function and psychopathology: it is not clear whether sub-specializations in cerebellar anatomy pertaining to components of distinct cognitive function exist and whether such differences can be observed independently of psychopathology severity.
In the current investigation, we examined for the first time how cerebellar regional anatomy may support cognitive function, capitalizing on a large dataset of transdiagnostic population employing a dimensional approach in agreement with the Research Domain Criteria framework (Cuthbert 2014). Our aim was to first outline gray matter volume changes in the cerebellum across distinct facets of cognitive function (e.g., executive function, working memory, cognitive flexibility, processing speed). Considering we were interested in cohorts with detailed cognitive phenotyping, we used the Healthy Brain Network (HBN) (Alexander et al. 2017), a landmark transdiagnostic mental health neuroimaging and behavioral dataset in a few thousand children and adolescents. The HBN includes predominantly unmedicated children and teenagers allowing us to examine cognitive cerebellar correlates unconfounded by chronic psychotropic consumption. Because spurious results can arise from quality control issues regarding neuroimaging scans, a rigorous quality assessment with visual inspection of all images is key to ensure the robustness of the results (Laidi et al. 2022). We used a data-driven multivariate canonical correlation analysis model (CCA) to evaluate the association between cerebellar anatomy and cognitive phenotype. Importantly, we used both permutation testing and bootstrapping to assess the significance and the robustness of our results.
Methods
Subjects
In the current investigation, we used data coming from an openly shared dataset, the Healthy Brain Network (HBN) project (Alexander et al. 2017). The HBN is a transdiagnostic dataset of neuroimaging and psychopathological assessments from a cohort of psychiatric or at-risk population of children and adolescents (5-21 years) (Alexander et al. 2017). Participants with severe neurological disorder or acute psychotic episodes are excluded in this cohort. In our study, considering our focus was on neurocognitive functioning, we excluded subjects with an intellectual deficiency (age-corrected IQ below 70, as measured with the Wechsler Adult Intelligence Scale (WASI-II) or the Wechsler Intelligence Scale for Children (WISC-V) (Na and Burns 2016). The full clinical assessment of the HBN cohort is described elsewhere in depth (Alexander et al. 2017).
Assessments
We were motivated to examine distinct contributions of the specific cognitive function aspects to the cerebellar anatomy. In the HBN database, cognition is quantified via the NIH Toolbox Cognition domain (Weintraub et al. 2013). The four subscales of the NIH Toolbox Cognition domain used in the HBN are detailed as follows. 1) The NIH Flanker assesses inhibitory executive control and attention and requires participants to focus on a target stimulus and ignore flanking stimuli. 2) NIH Card Sort assesses cognitive flexibility and requires participants to apply one rule to two target pictures (e.g., matching by color) and then another (e.g., matching by shape). 3) NIH List assesses working memory function and requires participants to sequence visually and orally presented stimuli e.g. by size. 4) NIH Pattern Comparison Processing Speed Test assesses processing speed by requiring participants to compare two side-by-side pictures (same vs. different). In all subscales, higher scores mean better ability. We used the standardized normative scores for each subscale. Psychopathology severity is quantified with the Child Behavior Checklist (CBCL) (Achenbach and Rescorla 2001), a widely employed scale that measures emotional, behavioral and social problems in children and teenagers of 1.5 - 18 years old. These scales have a mean t-score of 50 with a standard deviation of 10. A t-score ≤ 64 indicates non-clinical symptoms, a t-score between 65 and 69 indicates problems rated high enough to be of concern but not overtly deviant, a t-score ≥70 indicates clinical symptoms (Achenbach and Rescorla 2001; Carta et al. 2020).
MRI acquisition
Acquisition of MRI scans was done in three sites in New York City: Staten Island, Rutgers University and Cornell Brain Imaging Center. The specific details of each acquisition protocol are as follows: Staten Island images were acquired on a 1.5T Siemens Avanto (TR = 2730 ms, TE = 1.64 ms, flip angle = 7°, slice number= 176, voxel dimensions = 1.0 × 1.0 × 1.0 mm3). Rutgers University images were acquired on a 3T Siemens Tim Trio (TR = 2500 ms, TE = 3.15 ms, flip angle = 8°, slice number= 224, voxel dimensions = 0.8 × 0.8 × 0.8 mm3). Cornell Brain Imaging Center images were acquired on a Siemens Prisma 3T MRI (TR = 2500 ms, TE = 3.15 ms, flip angle = 8°, slice number= 224, voxel dimensions = 0.8 × 0.8 × 0.8 mm3)
MRI processing
All subjects were processed using the CERES pipeline (Romero et al. 2017). This fully automated method relies on an multi-atlas patch-based strategy that has been compared with manual tracing and performs well compared to other segmentation methods (Carass et al. 2018). All structural T1 MRIs were processed by PC on a high computing performance cluster in Bordeaux, France. The CERES pipeline follows the parcellation protocol described in (Park et al. 2014) (), which provides a parcellation of the cerebellum and gray matter volumes for all cerebellar Lobules, with the exception of the cerebellar vermis, which is included in every Lobule. Moreover, the CERES pipeline provides a mask of intracranial volume (ICV) (Manjón et al. 2014).
Quality control
The quality control procedure was done in two steps: 1) visual inspection of the raw T1 images 2) visual inspection of the images issued from the parcellation procedure in every slice for each spatial plan of the cerebellum by an expert rater (YE) - blind to the clinical features of each participant. We identified subjects with non-cerebellar voxels labeled as voxels belonging to the cerebellum, and vice versa, and subjects with parcellation errors within the cerebellar lobules. The same procedure has been applied previously (Laidi et al. 2022). No images with parcellation defects were included in further analyses. After the preprocessing and quality control step we excluded 602 individuals after visual inspection of the raw T1 images (280 individuals), parcellation errors (279 individuals), low IQ (43 individuals) and incomplete psychometric scores (60 individuals). A final sample of 662 individuals was included in the subsequent neuroimaging analyses. A summary of the repartition of the excluded subjects can be found in supplementary figure S1.
Statistical analyses
Canonical correlation analysis
We performed our multivariate analyses with scikit-learn library (Pedregosa et al. 2011). We employed a regularized kernel Canonical Correlation Analysis (CCA), using an open-source python pyrcca package (Bilenko and Gallant 2016), as a multivariate approach to evaluate the association between cerebellar anatomy (component with anatomical features “X”) and cognitive phenotype (component with clinical scores of interest “Y”). In brief, CCA solves the canonical spaces in which the maximal correlation of projected datasets occurs, not preassuming the directionality of the relationship between datasets (Hardoon et al. 2007). One documented disadvantage of CCA is overfitting to noise correlation of the datasets. To overcome this limitation, we implemented the algorithm used by the pyrcca toolbox. This algorithm constrains the number of components and find the optimal regularization parameters in a data-driven manner using a 10-fold cross validation approach (Bilenko and Gallant 2016). Here, we investigated a range of components ([2, 3, 4]) and a range of parameters ([0.0001, 0.01, 1, 100]). Such step resulted in the best number of components of 2 and best regularization parameter of 0.0001 that were used for the subsequent analyses.
In our analysis, the anatomical component included our regions of interest, namely the anterior lobe (Lobules I – V), Lobule VI, Crus I, Crus II, Lobule VIIB, VIIIA, IX and X and the ICV, after regressing out the effect of scan location, age and sex. We included in the cognitive component the standardized scores of the subscales of NIH toolbox, namely List subscale indexing working memory, Card subscale indexing cognitive flexibility, Flanker’s subscale indexing cognitive control and Processing subscale indexing processing speed. We then computed the correlation between the two canonical components for clinical and anatomical features. We repeated these analyses by including psychopathology severity quantified by the total t-score of the Child Behavior Checklist (CBCL) (Mazefsky et al. 2011).
Assessment of statistical significance and model stability
To assess statistical significance, we used non-parametric permutation testing (Ericson and Zinoviev 2001). Permutation testing involves random rearrangement of samples without replacement in order to estimate the population distribution and in turn, test the null hypothesis. Thus, p-value would be defined as the proportion of permuted samples that test statistically higher than our observed sample. In this study, unless specified otherwise, the threshold of significance was set to p=0.05 corresponding to the r value higher than the r value of the 95 percentile in a 10’ 000 random permutation test.
Crucially, we assessed the model stability through bootstrapping analysis (Carpenter and Bithell 2000). Bootstrapping is used to create a sampling distribution by repeatedly taking random samples with replacement from the original sample. We then performed CCA on each bootstrapped sample and collected these results to perform summary statistics (mean and confidence intervals). The average estimated from these multiple random samples may be used to infer results regarding the robustness of the CCA results for the original sample. Here, we used a bootstrap of 10’000 random samples with replacement which showed a normal distribution of canonical correlations. The results were considered stable if the 95% confidence interval of the bootstrap distribution would not include zero.
Results
Study sample: demographics and psychopathology
The demographic characteristics of our sample can be found in Table 1. In sum, the cohort included in our study (n=662), had a mean age of 10.5 years [range 5.82 – 17.74 years old] and were predominantly males (58%).
Next, we investigated the relationship between cognitive function assessed with the NIH-TB subscales and psychopathology severity quantified by the CBCL total t-score. We found positive correlations between the individual cognitive subscales scores with small effect sizes (Figure 2), indicating that these subscales measure little overlapping constructs. However, there were no associations of the cognitive subscales with psychopathology severity (CBCL) indicating no significant impact of psychopathology on distinct components of cognitive function.
Cerebellar correlates of cognitive function
We then turned to the brain to examine whether we could uncover latent neural and behavioral dimensions to our data with a multivariate CCA approach (including age, sex, scan location, total intracranial volume, cognition (NIH-TB subscales). Such an approach allowed us to identify two significant correlation between the first cognitive canonical variate (the subscales of the NIH toolbox Cognition domain) and the first brain canonical variate (regional cerebellar gray matter volume and intracranial volume) at r=.22, as well as between the second clinical canonical variate (the subscales of the NIH toolbox Cognition domain) and the second brain canonical variate (regional cerebellar gray matter volume and intracranial volume) at r=.16. To assess significance, we conducted permutation testing with 10’000 tests at 95% that revealed that the correlation of both components was significant (statistical threshold for component 1: r = .20, and for component 2: r = .15) (Figure 2).
Next we performed a bootstrapping analysis to examine the stability of our multivariate CCA model (Sauerbrei and Royston 2007). Our results showed that our CCA results are stable and non-zero within the [5:95] confidence interval of the results generated by bootstrapping analysis.
Regarding the first component, cognitive flexibility (indexed by the NIH card subscale), processing speed (indexed by NIH processing scale) and working memory (indexed by the NIH list subscale) loaded the most on the clinical canonical variate at 0.89 (large effect), 0.65 and 0.52 (moderate effects), respectively. Gray matter volume in the Crus II and Lobule X loaded with a moderate effect size on the brain canonical variate at 0.57 and 0.59, respectively (Figure 3). Regarding the second component, working memory (indexed by the NIH list subscale) and cognitive control (indexed by NIH Flanker subscale) loaded with a moderate effect on the cognitive canonical variate at 0.51 and -0.56 respectively. Gray matter volume in the Crus I and Lobule VI loaded also moderately on the brain canonical variate at 0.49 each. Overall, these results
We then asked whether there the differences in cerebellar volumes observed with cognitive function could be explained by psychopathology severity. To address this, we stratified the total CBCL t-score in the established normative cut-offs of the CBCL with a t-score of ≤ 64 for non-clinical symptoms, a t-score between 65 and 69 for borderline individuals with risk for problem behaviors, and a t-score ≥ 70 for clinical symptoms (Achenbach and Rescorla 2001; Carta et al. 2020). We then computed the effect size (Cohen’s d) of differences in cerebellar volume using standardized mean differences accounting for age, sex, scan location and intracranial volume. We found no significant differences in cerebellar volumes across these categories (Figure 3D) showing a lack of psychopathology effects in cerebellar anatomy. Furthermore, we included the CBCL total t-scores and the NIH subscales in the same CCA model with permutation and bootstrapping. This analysis did not impact our main cognitive-cerebellum CCA results remaining both statistically significant and stable (Supplementary material, Figure S3). We asked then whether specific dimensions of psychopathology could drive cerebellar changes and performed CCA analyses with the subscales of the CBCL. These results were not statistically significant and did not survive bootstrapping (Supplementary material, Figure S4). In sum, our results indicated that cognitive function but not psychopathology severity drove the observed cerebellar changes.
Discussion
Cognitive neuroscience is only beginning to unravel the role of the cerebellum in higher “supratentorial” cognitive functions. Despite historically being framed as a “motor control” brain region, extensive human neuroimaging and lesion evidence has suggested a cerebellar role in cognitive function. Furthermore, previous evidence uncovered a correspondence of cerebellar anatomy with general cognition and psychopathology. Therefore, we set out to examine for the first time the cerebellar topography in connection to specific components of cognition. Our multivariate analyses (CCA) outlined how different components of cognitive function map onto cerebellar morphometry independently of psychopathology severity in support of the cerebellar cognitive and affective syndrome (Schmahmann and Sherman 1998; Jacobi et al. 2021). In particular, we showed only partly shared cerebellar maps of cognitive function (Figure 3, A-C): a first map encompassing cognitive flexibility (large effect size) and speed of processing (moderate effect size) associated with regional gray matter volume in Crus II and Lobule X and a second map including the Crus I and Lobule VI associated with cognitive control (moderate effect size). Working memory associations were similarly present in both these maps (Crus II, Lobule X, Crus I and Lobule VI) with similar moderate effect sizes. These results account for psychopathology severity and other confounds and suggest that such correspondence between cerebellar anatomy may go across transdiagnostic boundaries. Crucially, permutation testing and bootstrapping analyses showed that these relationships are significant and our CCA model is robust and stable (Sauerbrei and Royston 2007).
Our findings highlight an association of the cerebellar Crus II and Lobule X structure with cognitive flexibility with a large effect size. Importantly, we show that cerebellar gray matter volume is linked to cognitive flexibility abilities in a dimensional across-diagnostic categories manner. While theoretical accounts have posited a role for the cerebellum in the flexible coordination of cognitive processes (Cognitive Dysmetria theory), strikingly, no large-scaled evidence for such contribution existed prior to our work. Animal research evidence indicated that hemi-cerebellectomized animals are unable to flexibly switch to a new set of rules, despite having intact motor responses (De Bartolo et al. 2009; Dickson et al. 2017). A few small-sample human studies (Ben-Soussan et al. 2015; Kansal et al. 2017; Badaly et al. 2022; Kühn et al. 2012; Paradiso et al. 1997; Parker et al. 2008; Koppelmans et al. 2017; Bauer et al. 2009; Bernard and Seidler 2013) have underlined a cerebellar role in mediating cognitive flexibility. The cerebellar correspondence of cognitive flexibility fits well with evidence from cerebellar lesions (Stoodley and Schmahmann 2010; Argyropoulos et al. 2020; Stoodley et al. 2016) and theoretical accounts of cognitive dysmetria theory (Andreasen, Paradiso, and O’Leary 1998). According to this theory, the cerebellum plays a key role in coordinating different cognitive and affective processes, quite similarly to its role in motor coordination. Impairments in cognitive flexibility seem to be prevalent in a variety of psychiatric disorders across the lifespan (O’Donnell et al. 2017; Dotson et al. 2020; Goodall et al. 2018; Geurts, Corbett, and Solomon 2009; Verdejo-Garcia et al. 2015) and they represent a potentially important pharmacological (Doss et al. 2021), psychotherapeutic (Johnco, Wuthrich, and Rapee 2014), or neuromodulation (de Boer et al. 2021) target.
Our results point to two cognitive-anatomical maps both related to working memory in agreement with previous findings (Hayter, Langdon, and Ramnani 2007; Guell, Gabrieli, and Schmahmann 2018; Marvel and Desmond 2016; Ashida et al. 2019; Begue et al. 2022). Furthermore, we show a first component characterizing the positive association of cognitive flexibility and speed of processing with regional gray matter volume in Crus II and Lobule X and a second, separate, component that captures the relationship of Crus I and Lobule VI with cognitive control (negative association). Speed of processing (the ability to quickly process information), working memory (the ability to hold and manipulate information during short periods of time) and cognitive flexibility (the ability to switch rapidly between mental states and tasks) are interconnected cognitive capacities that are important for flexible behavior. On the other hand, cognitive stability - the ability to maintain stable cognitive representations is also important for consistent and reliable behavior, and requires less task-switching and more working memory capacity. Previous literature has focused on supratentorial prefronto-striatal networks of working memory, cognitive flexibility and stability (Cools and D’Esposito 2010; Cools 2016; Westbrook and Braver 2016). In this circuitry, cognitive flexibility and cognitive stability have been shown to rely on dopaminergic signaling in the striatum, as demonstrated in human PET neuroimaging (Klostermann et al. 2012; Rieckmann et al. 2011; Landau et al. 2009) and task-related studies (Berry et al. 2016; Samanez-Larkin et al. 2013), and in the dorsolateral prefrontal cortices (Borwick et al. 2020). In line with these findings, Braver and Westbrook (Westbrook and Braver 2016) attributed a key role to dopaminergic neurons in the flexible coordination of cognitive processes for goal-directed behavior. In other words, during goal-directed behavior, phasic dopamine release promotes updating of these representations and therefore cognitive flexibility, whereas tonic dopamine release maintains the stability of cognitive representations. In light of our results, we speculatively propose that, similarly, dopaminergic signaling may underlie cerebellar participation in cognitive flexibility or stability either through direct local dopaminergic signaling in the cerebellum or via indirect (e.g., distant) cerebellar prefrontal closed loops as part of “the rich club” (Watson et al. 2014). Indeed, high levels of dopamine have been found in postmortem cerebellum of humans, rats, and monkeys reviewed in (Flace et al. 2021). Further research is needed to examine whether and how dopaminergic signaling underlies cerebellar correlates of cognitive function.
Our study has several strengths. First, we overcome previous shortcomings of case-control studies by endorsing a dimensional approach that captures phenotypic gradients in a large cohort. To our knowledge, this is the largest study to date to ever examine cerebellar contributions to distinct cognitive components including cognitive flexibility. Second, we employ a validated pipeline (CERES) with careful and stringent quality control to ensure optimal preprocessing and avoid spurious results. In the current work, we have included only high-quality imaging data surviving a stringent visual quality check (e.g., only 662 of the images have passed the quality check of the initial n = 1452 subjects, Figure S1) using the same quality - control protocol employed previously (Laidi et al. 2022). Third, this study provides statistically significant (permutation testing) and stable (bootstrapping) data-driven results in the largest transdiagnostic sample to date to ascertain the significance of the results and the stability of the CCA model. Permutation testing allowed us to establish that the obtained results are statistically significant. However, obtaining significant results does not exclude the possibility of having a random sampling error (e.g., the sample does not represent the general population). Such a possibility can be ruled out by bootstrapping: the results of the original cohort can be compared with bootstrapped samples, allowing us to examine the robustness of the model. Here, our results survived bootstrapping and we can confidently state that our results are not due to a random sampling error. Regarding limitations, our cross-sectional design and correlational analyses prevent inferences on the causal nature of the observed changes in the cerebellum. Furthermore, the findings may not generalize to older populations considering that in this cohort only young developing individuals (mean age of 10.5 years) were included. Future follow-up studies should examine how these relationships are expressed in adult cohorts.
Our work links cerebellar morphometry to distinct components of cognitive function including cognitive flexibility. These functions are observed to be altered in psychiatric disorders such as schizophrenia, depression, autism, and obsessive-compulsive disorders (Dajani and Uddin 2015), all of which also are shown to have cerebellar aberrations (Phillips et al. 2015). Given the recent advance in cerebellar non-invasive brain stimulation and its association with neuroimaging (Brady et al. 2019; Yao et al. 2022), our work opens the perspective of cerebellar targeting across different psychiatric diagnoses for cognitive improvement. Overall, our results elucidate for the first time the cerebellar anatomical circuitry supporting inter-individual differences in cognitive function and highlight a prominent role for the human cerebellum in distinct aspects of cognition for flexible adaptive behavior.
Data Availability
data are available on request from the authors.
Supplementary material
Supplementary Material 1. Linear univariate analyses
We also examined post-hoc linear univariate relationships of cerebellar regional gray matter to cognitive abilities measured by NIH toolbox subscales controlling for age, ICV, and psychopathology severity (indexed by CBCL total scores) as continuous variables and sex, scan location as categorical variables (Figure S2). We found a main positive effect of NIH List subscale indexing working memory in gray matter volumes of the anterior lobe (p=0. 015), Lobules VI (p = 0.0008), Crus I (p = 0.0008), Crus II (p =0.001), VIIB (p =0.0008), VIIIA (p=0.003) and X (p = 0.0008). Further, we found a main positive effect of NIH Card indexing cognitive flexibility in Crus II (p =.001), VIIB (p=.015), VIIIA (p=0.03) and Lobule X (p=.004). Cognitive control (indexed by NIH Flanker’s subscore) correlated with gray matter volume in the anterior lobe gray matter volume (p=0.04) and Crus II (p=0.03) whereas processing speed (indexed by NIH processing subscore) positively correlated to gray matter volume in Crus II (p=0.04) and Lobule X (p= 0.04) (all p-values are FDR corrected).