Spatial and temporal mapping of de novo mutations in schizophrenia to a fetal prefrontal cortical network.

Genes disrupted in schizophrenia may be revealed by de novo mutations in affected persons from otherwise healthy families. Furthermore, during normal brain development, genes are expressed in patterns specific to developmental stage and neuroanatomical structure. We identified de novo mutations in persons with schizophrenia and then mapped the responsible genes onto transcriptome profiles of normal human brain tissues from age 13 weeks gestation to adulthood. In the dorsolateral and ventrolateral prefrontal cortex during fetal development, genes harboring damaging de novo mutations in schizophrenia formed a network significantly enriched for transcriptional coexpression and protein interaction. The 50 genes in the network function in neuronal migration, synaptic transmission, signaling, transcriptional regulation, and transport. These results suggest that disruptions of fetal prefrontal cortical neurogenesis are critical to the pathophysiology of schizophrenia. These results also support the feasibility of integrating genomic and transcriptome analyses to map critical neurodevelopmental processes in time and space in the brain.

In vivo whole-cortex marker of excitation-inhibition ratio indexes cortical maturation and cognitive ability in youth.

A balanced excitation-inhibition ratio (E/I ratio) is critical for healthy brain function. Normative development of cortex-wide E/I ratio remains unknown. Here, we noninvasively estimate a putative marker of whole-cortex E/I ratio by fitting a large-scale biophysically plausible circuit model to resting-state functional MRI (fMRI) data. We first confirm that our model generates realistic brain dynamics in the Human Connectome Project. Next, we show that the estimated E/I ratio marker is sensitive to the gamma-aminobutyric acid (GABA) agonist benzodiazepine alprazolam during fMRI. Alprazolam-induced E/I changes are spatially consistent with positron emission tomography measurement of benzodiazepine receptor density. We then investigate the relationship between the E/I ratio marker and neurodevelopment. We find that the E/I ratio marker declines heterogeneously across the cerebral cortex during youth, with the greatest reduction occurring in sensorimotor systems relative to association systems. Importantly, among children with the same chronological age, a lower E/I ratio marker (especially in the association cortex) is linked to better cognitive performance. This result is replicated across North American (8.2 to 23.0 y old) and Asian (7.2 to 7.9 y old) cohorts, suggesting that a more mature E/I ratio indexes improved cognition during normative development. Overall, our findings open the door to studying how disrupted E/I trajectories may lead to cognitive dysfunction in psychopathology that emerges during youth.

Genomic loci influence patterns of structural covariance in the human brain.

Normal and pathologic neurobiological processes influence brain morphology in coordinated ways that give rise to patterns of structural covariance (PSC) across brain regions and individuals during brain aging and diseases. The genetic underpinnings of these patterns remain largely unknown. We apply a stochastic multivariate factorization method to a diverse population of 50,699 individuals (12 studies and 130 sites) and derive data-driven, multi-scale PSCs of regional brain size. PSCs were significantly correlated with 915 genomic loci in the discovery set, 617 of which are newly identified, and 72% were independently replicated. Key pathways influencing PSCs involve reelin signaling, apoptosis, neurogenesis, and appendage development, while pathways of breast cancer indicate potential interplays between brain metastasis and PSCs associated with neurodegeneration and dementia. Using support vector machines, multi-scale PSCs effectively derive imaging signatures of several brain diseases. Our results elucidate genetic and biological underpinnings that influence structural covariance patterns in the human brain.

ASLPrep: a platform for processing of arterial spin labeled MRI and quantification of regional brain perfusion.

Arterial spin labeled (ASL) magnetic resonance imaging (MRI) is the primary method for noninvasively measuring regional brain perfusion in humans. We introduce ASLPrep, a suite of software pipelines that ensure the reproducible and generalizable processing of ASL MRI data.

Psychosis superspectrum I: Nosology, etiology, and lifespan development.

This review describes the Hierarchical Taxonomy of Psychopathology (HiTOP) model of psychosis-related psychopathology, the psychosis superspectrum. The HiTOP psychosis superspectrum was developed to address shortcomings of traditional diagnoses for psychotic disorders and related conditions including low reliability, arbitrary boundaries between psychopathology and normality, high symptom co-occurrence, and heterogeneity within diagnostic categories. The psychosis superspectrum is a transdiagnostic dimensional model comprising two spectra-psychoticism and detachment-which are in turn broken down into fourteen narrow components, and two auxiliary domains-cognition and functional impairment. The structure of the spectra and their components are shown to parallel the genetic structure of psychosis and related traits. Psychoticism and detachment have distinct patterns of association with urbanicity, migrant and ethnic minority status, childhood adversity, and cannabis use. The superspectrum also provides a useful model for describing the emergence and course of psychosis, as components of the superspectrum are relatively stable over time. Changes in psychoticism predict the onset of psychosis-related psychopathology, whereas changes in detachment and cognition define later course. Implications of the superspectrum for genetic, socio-environmental, and longitudinal research are discussed. A companion review focuses on neurobiology, treatment response, and clinical utility of the superspectrum, and future research directions.

Polygenic risk of social isolation behavior and its influence on psychopathology and personality.

Social isolation has been linked to a range of psychiatric issues, but the behavioral component that drives it is not well understood. Here, a genome-wide associations study (GWAS) was carried out to identify genetic variants that contribute specifically to social isolation behavior (SIB) in up to 449,609 participants from the UK Biobank. 17 loci were identified at genome-wide significance, contributing to a 4% SNP-based heritability estimate. Using the SIB GWAS, polygenic risk scores (PRS) were derived in ALSPAC, an independent, developmental cohort, and used to test for association with self-reported friendship scores, comprising items related to friendship quality and quantity, at age 12 and 18 to determine whether genetic predisposition manifests during childhood development. At age 18, friendship scores were associated with the SIB PRS, demonstrating that the genetic factors can predict related social traits in late adolescence. Linkage disequilibrium (LD) score correlation using the SIB GWAS demonstrated genetic correlations with autism spectrum disorder (ASD), schizophrenia, major depressive disorder (MDD), educational attainment, extraversion, and loneliness. However, no evidence of causality was found using a conservative Mendelian randomization approach between SIB and any of the traits in either direction. Genomic Structural Equation Modeling (SEM) revealed a common factor contributing to SIB, neuroticism, loneliness, MDD, and ASD, weakly correlated with a second common factor that contributes to psychiatric and psychotic traits. Our results show that SIB contributes a small heritable component, which is associated genetically with other social traits such as friendship as well as psychiatric disorders.

Sex differences in the functional topography of association networks in youth.

Prior work has shown that there is substantial interindividual variation in the spatial distribution of functional networks across the cerebral cortex, or functional topography. However, it remains unknown whether there are sex differences in the topography of individualized networks in youth. Here, we leveraged an advanced machine learning method (sparsity-regularized non-negative matrix factorization) to define individualized functional networks in 693 youth (ages 8 to 23 y) who underwent functional MRI as part of the Philadelphia Neurodevelopmental Cohort. Multivariate pattern analysis using support vector machines classified participant sex based on functional topography with 82.9% accuracy (P < 0.0001). Brain regions most effective in classifying participant sex belonged to association networks, including the ventral attention, default mode, and frontoparietal networks. Mass univariate analyses using generalized additive models with penalized splines provided convergent results. Furthermore, transcriptomic data from the Allen Human Brain Atlas revealed that sex differences in multivariate patterns of functional topography were spatially correlated with the expression of genes on the X chromosome. These results highlight the role of sex as a biological variable in shaping functional topography.

Immune status of the murine 22q11.2 deletion syndrome model.

Mice modeling the hemizygous deletion of chromosome 22q11.2 (22qMc) have been utilized to address various clinical phenotypes associated with the disease, including cardiac malformations, altered neural circuitry, and behavioral deficits. Yet, the status of the T cell compartment, an important clinical concern among 22q11.2 deletion syndrome (22qDS) patients, has not been addressed. While infancy and early childhood in 22qDS are associated with deficient T cell numbers and thymic hypoplasia, which can be severe in a small subset of patients, studies suggest normalization of the T cell counts by adulthood. We found that adult 22qMc do not exhibit thymic hypoplasia or altered thymic T cell development. Our findings that immune cell counts and inflammatory T cell activation are unaffected in 22qMc lend support to the hypothesis that human 22qDS immunodeficiencies are secondary to thymic hypoplasia, rather than intrinsic effects due to the deletion. Furthermore, the 22q11.2 deletion does not impact the differentiation capacity of T cells, nor their activity and response during inflammatory activation. Thus, 22qMc reflects the T cell compartment in adult 22qDS patients, and our findings suggest that 22qMc may serve as a novel model to address experimental and translational aspects of immunity in 22qDS.

Network enrichment significance testing in brain-phenotype association studies.

Functional networks often guide our interpretation of spatial maps of brain-phenotype associations. However, methods for assessing enrichment of associations within networks of interest have varied in terms of both scientific rigor and underlying assumptions. While some approaches have relied on subjective interpretations, others have made unrealistic assumptions about spatial properties of imaging data, leading to inflated false positive rates. We seek to address this gap in existing methodology by borrowing insight from a method widely used in genetics research for testing enrichment of associations between a set of genes and a phenotype of interest. We propose network enrichment significance testing (NEST), a flexible framework for testing the specificity of brain-phenotype associations to functional networks or other sub-regions of the brain. We apply NEST to study enrichment of associations with structural and functional brain imaging data from a large-scale neurodevelopmental cohort study.

Source-based morphometry reveals structural brain pattern abnormalities in 22q11.2 deletion syndrome.

22q11.2 deletion syndrome (22q11DS) is the most frequently occurring microdeletion in humans. It is associated with a significant impact on brain structure, including prominent reductions in gray matter volume (GMV), and neuropsychiatric manifestations, including cognitive impairment and psychosis. It is unclear whether GMV alterations in 22q11DS occur according to distinct structural patterns. Then, 783 participants (470 with 22q11DS: 51% females, mean age [SD] 18.2 [9.2]; and 313 typically developing [TD] controls: 46% females, mean age 18.0 [8.6]) from 13 datasets were included in the present study. We segmented structural T1-weighted brain MRI scans and extracted GMV images, which were then utilized in a novel source-based morphometry (SBM) pipeline (SS-Detect) to generate structural brain patterns (SBPs) that capture co-varying GMV. We investigated the impact of the 22q11.2 deletion, deletion size, intelligence quotient, and psychosis on the SBPs. Seventeen GMV-SBPs were derived, which provided spatial patterns of GMV covariance associated with a quantitative metric (i.e., loading score) for analysis. Patterns of topographically widespread differences in GMV covariance, including the cerebellum, discriminated individuals with 22q11DS from healthy controls. The spatial extents of the SBPs that revealed disparities between individuals with 22q11DS and controls were consistent with the findings of the univariate voxel-based morphometry analysis. Larger deletion size was associated with significantly lower GMV in frontal and occipital SBPs; however, history of psychosis did not show a strong relationship with these covariance patterns. 22q11DS is associated with distinct structural abnormalities captured by topographical GMV covariance patterns that include the cerebellum. Findings indicate that structural anomalies in 22q11DS manifest in a nonrandom manner and in distinct covarying anatomical patterns, rather than a diffuse global process. These SBP abnormalities converge with previously reported cortical surface area abnormalities, suggesting disturbances of early neurodevelopment as the most likely underlying mechanism.

Similarity-based multimodal regression.

To better understand complex human phenotypes, large-scale studies have increasingly collected multiple data modalities across domains such as imaging, mobile health, and physical activity. The properties of each data type often differ substantially and require either separate analyses or extensive processing to obtain comparable features for a combined analysis. Multimodal data fusion enables certain analyses on matrix-valued and vector-valued data, but it generally cannot integrate modalities of different dimensions and data structures. For a single data modality, multivariate distance matrix regression provides a distance-based framework for regression accommodating a wide range of data types. However, no distance-based method exists to handle multiple complementary types of data. We propose a novel distance-based regression model, which we refer to as Similarity-based Multimodal Regression (SiMMR), that enables simultaneous regression of multiple modalities through their distance profiles. We demonstrate through simulation, imaging studies, and longitudinal mobile health analyses that our proposed method can detect associations between clinical variables and multimodal data of differing properties and dimensionalities, even with modest sample sizes. We perform experiments to evaluate several different test statistics and provide recommendations for applying our method across a broad range of scenarios.

QSIPrep: an integrative platform for preprocessing and reconstructing diffusion MRI data.

Diffusion-weighted magnetic resonance imaging (dMRI) is the primary method for noninvasively studying the organization of white matter in the human brain. Here we introduce QSIPrep, an integrative software platform for the processing of diffusion images that is compatible with nearly all dMRI sampling schemes. Drawing on a diverse set of software suites to capitalize on their complementary strengths, QSIPrep facilitates the implementation of best practices for processing of diffusion images.

Domain-specific associations between psychopathology and neurocognitive functioning.

BACKGROUND: Neurocognitive dysfunction is a transdiagnostic finding in psychopathology, but relationships among cognitive domains and general and specific psychopathology dimensions remain unclear. This study aimed to examine associations between cognition and psychopathology dimensions in a large youth cohort. METHOD: The sample (N = 9350; age 8-21 years) was drawn from the Philadelphia Neurodevelopmental Cohort. Data from structured clinical interviews were modeled using bifactor confirmatory factor analysis (CFA), resulting in an overall psychopathology ('p') factor score and six orthogonal psychopathology dimensions: dysphoria/distress, obsessive-compulsive, behavioral/externalizing, attention-deficit/hyperactivity, phobias, and psychosis. Neurocognitive data were aggregated using correlated-traits CFA into five factors: executive functioning, memory, complex cognition, social cognition, and sensorimotor speed. We examined relationships among specific and general psychopathology dimensions and neurocognitive factors. RESULTS: The final model showed both overall and specific associations between cognitive functioning and psychopathology, with acceptable fit (CFI = 0.91; TLI = 0.90; RMSEA = 0.024; SRMR = 0.054). Overall psychopathology and most psychopathology dimensions were negatively associated with neurocognitive functioning (phobias [p < 0.0005], behavioral/externalizing [p < 0.0005], attention-deficit/hyperactivity [p < 0.0005], psychosis [p < 0.0005 to p < 0.05]), except for dysphoria/distress and obsessive-compulsive symptoms, which were positively associated with complex cognition (p < 0.05 and p < 0.01, respectively). CONCLUSION: By modeling a broad range of cognitive and psychopathology domains in a large, diverse sample of youth, we found aspects of neurocognitive functioning shared across clinical phenotypes, as well as domain-specific patterns. Findings support transdiagnostic examination of cognitive performance to parse variability in the link between neurocognitive functioning and clinical phenotypes.

Immune and oxidative stress biomarkers in pediatric psychosis and psychosis-risk: Meta-analyses and systematic review.

OBJECTIVE: While genetic and cohort studies suggest immune and reduction/oxidation (redox) alterations occur in psychosis, less is known about potential alterations in children and adolescents. METHODS: We conducted a systematic review to identify immune and redox biomarker studies in children and adolescents (mean age ≤ 18 years old) across the psychosis spectrum: from psychotic like experiences, which are common in children, to threshold psychotic disorders like schizophrenia. We conducted meta-analyses when at least three studies measured the same biomarker. RESULTS: The systematic review includes 38 pediatric psychosis studies. The meta-analyses found that youth with threshold psychotic disorders had higher neutrophil/lymphocyte ratio (Hedge's g = 0.40, 95 % CI 0.17 - 0.64), tumor necrosis factor (Hedge's g = 0.38, 95 % CI 0.06 - 0.69), C-reactive protein (Hedge's g = 0.38, 95 % CI 0.05 - 0.70), interleukin-6 (Hedge's g = 0.35; 95 % CI 0.11 - 0.64), and total white blood cell count (Hedge's g = 0.29, 95 % CI 0.12 - 0.46) compared to youth without psychosis. Other immune and oxidative stress meta-analytic findings were very heterogeneous. CONCLUSION: Results from several studies are consistent with the hypothesis that signals often classified as "proinflammatory" are elevated in threshold pediatric psychotic disorders. Data are less clear for immune markers in subthreshold psychosis and redox markers across the subthreshold and threshold psychosis spectrum. Immune and redox biomarker intervention studies are lacking, and research investigating interventions targeting the immune system in threshold pediatric psychosis is especially warranted.

Traumatic stress load and stressor reactivity score associated with accelerated gray matter maturation in youths indexed by normative models.

Understanding how traumatic stress affects typical brain development during adolescence is critical to elucidate underlying mechanisms related to both maladaptive functioning and resilience after traumatic exposures. The current study aimed to map deviations from normative ranges of brain gray matter for youths with traumatic exposures. For each cortical and subcortical gray matter region, normative percentiles of variations were established using structural MRI from typically developing youths without any traumatic exposure (n = 245; age range = 8-23) from the Philadelphia Neurodevelopmental Cohort (PNC). The remaining PNC participants with neuroimaging data (n = 1129) were classified as either within the normative range (5-95%), delayed (>95%) or accelerated (<5%) maturational ranges for each region using the normative model. An averaged quantile regression index was calculated across all regions. Mediation models revealed that high traumatic stress load was positively associated with poorer cognitive functioning and greater psychopathology, and these associations were mediated by accelerated gray matter maturation. Furthermore, higher stressor reactivity scores, which represent a less resilient response under traumatic stress, were positively correlated with greater acceleration of gray matter maturation (r = 0.224, 95% CI = [0.17, 0.28], p < 0.001), suggesting that more accelerated maturation was linked to greater stressor response regardless of traumatic stress load. We conclude that traumatic stress is a source of deviation from normative brain development associated with poorer cognitive functioning and more psychopathology in the long run.

Alprazolam modulates persistence energy during emotion processing in first-degree relatives of individuals with schizophrenia: a network control study.

Schizophrenia is marked by deficits in facial affect processing associated with abnormalities in GABAergic circuitry, deficits also found in first-degree relatives. Facial affect processing involves a distributed network of brain regions including limbic regions like amygdala and visual processing areas like fusiform cortex. Pharmacological modulation of GABAergic circuitry using benzodiazepines like alprazolam can be useful for studying this facial affect processing network and associated GABAergic abnormalities in schizophrenia. Here, we use pharmacological modulation and computational modeling to study the contribution of GABAergic abnormalities toward emotion processing deficits in schizophrenia. Specifically, we apply principles from network control theory to model persistence energy - the control energy required to maintain brain activation states - during emotion identification and recall tasks, with and without administration of alprazolam, in a sample of first-degree relatives and healthy controls. Here, persistence energy quantifies the magnitude of theoretical external inputs during the task. We find that alprazolam increases persistence energy in relatives but not in controls during threatening face processing, suggesting a compensatory mechanism given the relative absence of behavioral abnormalities in this sample of unaffected relatives. Further, we demonstrate that regions in the fusiform and occipital cortices are important for facilitating state transitions during facial affect processing. Finally, we uncover spatial relationships (i) between regional variation in differential control energy (alprazolam versus placebo) and (ii) both serotonin and dopamine neurotransmitter systems, indicating that alprazolam may exert its effects by altering neuromodulatory systems. Together, these findings provide a new perspective on the distributed emotion processing network and the effect of GABAergic modulation on this network, in addition to identifying an association between schizophrenia risk and abnormal GABAergic effects on persistence energy during threat processing.

Rare coding variants as risk modifiers of the 22q11.2 deletion implicate postnatal cortical development in syndromic schizophrenia.

22q11.2 deletion is one of the strongest known genetic risk factors for schizophrenia. Recent whole-genome sequencing of schizophrenia cases and controls with this deletion provided an unprecedented opportunity to identify risk modifying genetic variants and investigate their contribution to the pathogenesis of schizophrenia in 22q11.2 deletion syndrome. Here, we apply a novel analytic framework that integrates gene network and phenotype data to investigate the aggregate effects of rare coding variants and identified modifier genes in this etiologically homogenous cohort (223 schizophrenia cases and 233 controls of European descent). Our analyses revealed significant additive genetic components of rare nonsynonymous variants in 110 modifier genes (adjusted P = 9.4E-04) that overall accounted for 4.6% of the variance in schizophrenia status in this cohort, of which 4.0% was independent of the common polygenic risk for schizophrenia. The modifier genes affected by rare coding variants were enriched with genes involved in synaptic function and developmental disorders. Spatiotemporal transcriptomic analyses identified an enrichment of coexpression between modifier and 22q11.2 genes in cortical brain regions from late infancy to young adulthood. Corresponding gene coexpression modules are enriched with brain-specific protein-protein interactions of SLC25A1, COMT, and PI4KA in the 22q11.2 deletion region. Overall, our study highlights the contribution of rare coding variants to the SCZ risk. They not only complement common variants in disease genetics but also pinpoint brain regions and developmental stages critical to the etiology of syndromic schizophrenia.

Psychosis brain subtypes validated in first-episode cohorts and related to illness remission: results from the PHENOM consortium.

Using machine learning, we recently decomposed the neuroanatomical heterogeneity of established schizophrenia to discover two volumetric subgroups-a 'lower brain volume' subgroup (SG1) and an 'higher striatal volume' subgroup (SG2) with otherwise normal brain structure. In this study, we investigated whether the MRI signatures of these subgroups were also already present at the time of the first-episode of psychosis (FEP) and whether they were related to clinical presentation and clinical remission over 1-, 3-, and 5-years. We included 572 FEP and 424 healthy controls (HC) from 4 sites (Sao Paulo, Santander, London, Melbourne) of the PHENOM consortium. Our prior MRI subgrouping models (671 participants; USA, Germany, and China) were applied to both FEP and HC. Participants were assigned into 1 of 4 categories: subgroup 1 (SG1), subgroup 2 (SG2), no subgroup membership ('None'), and mixed SG1 + SG2 subgroups ('Mixed'). Voxel-wise analyses characterized SG1 and SG2 subgroups. Supervised machine learning analyses characterized baseline and remission signatures related to SG1 and SG2 membership. The two dominant patterns of 'lower brain volume' in SG1 and 'higher striatal volume' (with otherwise normal neuromorphology) in SG2 were identified already at the first episode of psychosis. SG1 had a significantly higher proportion of FEP (32%) vs. HC (19%) than SG2 (FEP, 21%; HC, 23%). Clinical multivariate signatures separated the SG1 and SG2 subgroups (balanced accuracy = 64%; p < 0.0001), with SG2 showing higher education but also greater positive psychosis symptoms at first presentation, and an association with symptom remission at 1-year, 5-year, and when timepoints were combined. Neuromorphological subtypes of schizophrenia are already evident at illness onset, separated by distinct clinical presentations, and differentially associated with subsequent remission. These results suggest that the subgroups may be underlying risk phenotypes that could be targeted in future treatment trials and are critical to consider when interpreting neuroimaging literature.

Validation of the cognitive section of the Penn computerized adaptive test for neurocognitive and clinical psychopathology assessment (CAT-CCNB).

BACKGROUND: The Penn Computerized Neurocognitive Battery is an efficient tool for assessing brain-behavior domains, and its efficiency was augmented via computerized adaptive testing (CAT). This battery requires validation in a separate sample to establish psychometric properties. METHODS: In a mixed community/clinical sample of N = 307 18-to-35-year-olds, we tested the relationships of the CAT tests with the full-form tests. We compared discriminability among recruitment groups (psychosis, mood, control) and examined how their scores relate to demographics. CAT-Full relationships were evaluated based on a minimum inter-test correlation of 0.70 or an inter-test correlation within at least 0.10 of the full-form correlation with a previous administration of the full battery. Differences in criterion relationships were tested via mixed models. RESULTS: Most tests (15/17) met the minimum criteria for replacing the full-form with the updated CAT version (mean r = 0.67; range = 0.53-0.80) when compared to relationships of the full-forms with previous administrations of the full-forms (mean r = 0.68; range = 0.50-0.85). Most (16/17) CAT-based relationships with diagnostics and other validity criteria were indistinguishable (interaction p > 0.05) from their full-form counterparts. CONCLUSIONS: The updated CNB shows psychometric properties acceptable for research. The full-forms of some tests should be retained due to insufficient time savings to justify the loss in precision.

Associations between neighborhood socioeconomic status, parental education, and executive system activation in youth.

Socioeconomic status (SES) can impact cognitive performance, including working memory (WM). As executive systems that support WM undergo functional neurodevelopment during adolescence, environmental stressors at both individual and community levels may influence cognitive outcomes. Here, we sought to examine how SES at the neighborhood and family level impacts task-related activation of the executive system during adolescence and determine whether this effect mediates the relationship between SES and WM performance. To address these questions, we studied 1,150 youths (age 8-23) that completed a fractal n-back WM task during functional magnetic resonance imaging at 3T as part of the Philadelphia Neurodevelopmental Cohort. We found that both higher neighborhood SES and parental education were associated with greater activation of the executive system to WM load, including the bilateral dorsolateral prefrontal cortex, posterior parietal cortex, and precuneus. The association of neighborhood SES remained significant when controlling for task performance, or related factors like exposure to traumatic events. Furthermore, high-dimensional multivariate mediation analysis identified distinct patterns of brain activity within the executive system that significantly mediated the relationship between measures of SES and task performance. These findings underscore the importance of multilevel environmental factors in shaping executive system function and WM in youth.