Permutation inference for canonical correlation analysis

Dimensional and Categorical Solutions to Parsing Depression Heterogeneity in a Large Single-Site Sample

BACKGROUND

Recent studies have reported significant advances in modeling the biological basis of heterogeneity in major depressive disorder, but investigators have also identified important technical challenges, including scanner-related artifacts, a propensity for multivariate models to overfit, and a need for larger samples with more extensive clinical phenotyping. The goals of the current study were to evaluate dimensional and categorical solutions to parsing heterogeneity in depression that are stable and generalizable in a large, single-site sample.

METHODS

We used regularized canonical correlation analysis to identify data-driven brain-behavior dimensions that explain individual differences in depression symptom domains in a large, single-site dataset comprising clinical assessments and resting-state functional magnetic resonance imaging data for 328 patients with major depressive disorder and 461 healthy control participants. We examined the stability of clinical loadings and model performance in held-out data. Finally, hierarchical clustering on these dimensions was used to identify categorical depression subtypes.

RESULTS

The optimal regularized canonical correlation analysis model yielded 3 robust and generalizable brain-behavior dimensions that explained individual differences in depressed mood and anxiety, anhedonia, and insomnia. Hierarchical clustering identified 4 depression subtypes, each with distinct clinical symptom profiles, abnormal resting-state functional connectivity patterns, and antidepressant responsiveness to repetitive transcranial magnetic stimulation.

CONCLUSIONS

Our results define dimensional and categorical solutions to parsing neurobiological heterogeneity in major depressive disorder that are stable, generalizable, and capable of predicting treatment outcomes, each with distinct advantages in different contexts. They also provide additional evidence that regularized canonical correlation analysis and hierarchical clustering are effective tools for investigating associations between functional connectivity and clinical symptoms.

Brain Functional Connectome Defines a Transdiagnostic Dimension Shared by Cognitive Function and Psychopathology in Preadolescents

BACKGROUND

Cognitive function and general psychopathology are two important classes of human behavior dimensions that are individually related to mental disorders across diagnostic categories. However, whether these two transdiagnostic dimensions are linked to common or distinct brain networks that convey resilience or risk for the development of psychiatric disorders remains unclear.

METHODS

The current study is a longitudinal investigation with 11,875 youths from the Adolescent Brain Cognitive Development (ABCD) Study at ages 9 to 10 years at the onset of the study. A machine learning approach based on canonical correlation analysis was used to identify latent dimensional associations of the resting-state functional connectome with multidomain behavioral assessments including cognitive functions and psychopathological measures. For the latent resting-state functional connectivity factor showing a robust behavioral association, its ability to predict psychiatric disorders was assessed using 2-year follow-up data, and its genetic association was evaluated using twin data from the same cohort.

RESULTS

A latent functional connectome pattern was identified that showed a strong and generalizable association with the multidomain behavioral assessments (5-fold cross-validation: ρ = 0.68-0.73 for the training set [n = 5096]; ρ = 0.56-0.58 for the test set [n = 1476]). This functional connectome pattern was highly heritable (h2 = 74.42%, 95% CI: 56.76%-85.42%), exhibited a dose-response relationship with the cumulative number of psychiatric disorders assessed concurrently and at 2 years post-magnetic resonance imaging scan, and predicted the transition of diagnosis across disorders over the 2-year follow-up period.

CONCLUSIONS

These findings provide preliminary evidence for a transdiagnostic connectome-based measure that underlies individual differences in the development of psychiatric disorders during early adolescence.

Data-driven, generalizable prediction of adolescent sleep disturbances in the multisite Adolescent Brain Cognitive Development Study

STUDY OBJECTIVES

Sleep disturbances are common in adolescence and associated with a host of negative outcomes. Here, we assess associations between multifaceted sleep disturbances and a broad set of psychological, cognitive, and demographic variables using a data-driven approach, canonical correlation analysis (CCA).

METHODS

Baseline data from 9093 participants from the Adolescent Brain Cognitive Development (ABCD) Study were examined using CCA, a multivariate statistical approach that identifies many-to-many associations between two sets of variables by finding combinations for each set of variables that maximize their correlation. We combined CCA with leave-one-site-out cross-validation across ABCD sites to examine the robustness of results and generalizability to new participants. The statistical significance of canonical correlations was determined by non-parametric permutation tests that accounted for twin, family, and site structure. To assess the stability of the associations identified at baseline, CCA was repeated using 2-year follow-up data from 4247 ABCD Study participants.

RESULTS

Two significant sets of associations were identified: (1) difficulty initiating and maintaining sleep and excessive daytime somnolence were strongly linked to nearly all domains of psychopathology (r2 = 0.36, p < .0001); (2) sleep breathing disorders were linked to BMI and African American/black race (r2 = 0.08, p < .0001). These associations generalized to unseen participants at all 22 ABCD sites and were replicated using 2-year follow-up data.

CONCLUSIONS

These findings underscore interwoven links between sleep disturbances in early adolescence and psychological, social, and demographic factors.

Neuroanatomical, transcriptomic, and molecular correlates of math ability and their prognostic value for predicting learning outcomes

Foundational mathematical abilities, acquired in early childhood, are essential for success in our technology-driven society. Yet, the neurobiological mechanisms underlying individual differences in children's mathematical abilities and learning outcomes remain largely unexplored. Leveraging one of the largest multicohort datasets from children at a pivotal stage of knowledge acquisition, we first establish a replicable mathematical ability-related imaging phenotype (MAIP). We then show that brain gene expression profiles enriched for candidate math ability-related genes, neuronal signaling, synaptic transmission, and voltage-gated potassium channel activity contributed to the MAIP. Furthermore, the similarity between MAIP gene expression signatures and brain structure, acquired before intervention, predicted learning outcomes in two independent math tutoring cohorts. These findings advance our knowledge of the interplay between neuroanatomical, transcriptomic, and molecular mechanisms underlying mathematical ability and reveal predictive biomarkers of learning. Our findings have implications for the development of personalized education and interventions.

Subthalamic nucleus-language network connectivity predicts dopaminergic modulation of speech function in Parkinson's disease

Speech impediments are a prominent yet understudied symptom of Parkinson's disease (PD). While the subthalamic nucleus (STN) is an established clinical target for treating motor symptoms, these interventions can lead to further worsening of speech. The interplay between dopaminergic medication, STN circuitry, and their downstream effects on speech in PD is not yet fully understood. Here, we investigate the effect of dopaminergic medication on STN circuitry and probe its association with speech and cognitive functions in PD patients. We found that changes in intrinsic functional connectivity of the STN were associated with alterations in speech functions in PD. Interestingly, this relationship was characterized by altered functional connectivity of the dorsolateral and ventromedial subdivisions of the STN with the language network. Crucially, medication-induced changes in functional connectivity between the STN's dorsolateral subdivision and key regions in the language network, including the left inferior frontal cortex and the left superior temporal gyrus, correlated with alterations on a standardized neuropsychological test requiring oral responses. This relation was not observed in the written version of the same test. Furthermore, changes in functional connectivity between STN and language regions predicted the medication's downstream effects on speech-related cognitive performance. These findings reveal a previously unidentified brain mechanism through which dopaminergic medication influences speech function in PD. Our study sheds light into the subcortical-cortical circuit mechanisms underlying impaired speech control in PD. The insights gained here could inform treatment strategies aimed at mitigating speech deficits in PD and enhancing the quality of life for affected individuals.

Multimodal neuroimaging correlates of physical-cognitive covariation in Chilean adolescents. The Cogni-Action Project

Health-related behaviours have been related to brain structural features. In developing settings, such as Latin America, high social inequality has been inversely associated with several health-related behaviours affecting brain development. Understanding the relationship between health behaviours and brain structure in such settings is particularly important during adolescence when critical habits are acquired and ingrained. In this cross-sectional study, we carry out a multimodal analysis identifying a brain region associated with health-related behaviours (i.e., adiposity, fitness, sleep problems and others) and cognitive/academic performance, independent of socioeconomic status in a large sample of Chilean adolescents. Our findings suggest that the relationship between health behaviours and cognitive/academic performance involves a particular brain phenotype that could play a mediator role. These findings fill a significant gap in the literature, which has largely focused on developed countries and raise the possibility of promoting healthy behaviours in adolescence as a means to influence brain structure and thereby cognitive/academic achievement, independently of socioeconomic factors. By highlighting the potential impact on brain structure and cognitive/academic achievement, policymakers could design interventions that are more effective in reducing health disparities in developing countries.

Leveraging multivariate approaches to advance the science of early-life adversity

Since the landmark Adverse Childhood Experiences (ACEs) study, adversity research has expanded to more precisely account for the multifaceted nature of adverse experiences. The complex data structures and interrelated nature of adversity data require robust multivariate statistical methods, and recent methodological and statistical innovations have facilitated advancements in research on childhood adversity. Here, we provide an overview of a subset of multivariate methods that we believe hold particular promise for advancing the field's understanding of early-life adversity, and discuss how these approaches can be practically applied to explore different research questions. This review covers data-driven or unsupervised approaches (including dimensionality reduction and person-centered clustering/subtype identification) as well as supervised/prediction-based approaches (including linear and tree-based models and neural networks). For each, we highlight studies that have effectively applied the method to provide novel insight into early-life adversity. Taken together, we hope this review serves as a resource to adversity researchers looking to expand upon the cumulative approach described in the original ACEs study, thereby advancing the field's understanding of the complexity of adversity and related developmental consequences.

Neurobiologically realistic neural network enables cross-scale modeling of neural dynamics

Fundamental principles underlying computation in multi-scale brain networks illustrate how multiple brain areas and their coordinated activity give rise to complex cognitive functions. Whereas brain activity has been studied at the micro- to meso-scale to reveal the connections between the dynamical patterns and the behaviors, investigations of neural population dynamics are mainly limited to single-scale analysis. Our goal is to develop a cross-scale dynamical model for the collective activity of neuronal populations. Here we introduce a bio-inspired deep learning approach, termed NeuroBondGraph Network (NBGNet), to capture cross-scale dynamics that can infer and map the neural data from multiple scales. Our model not only exhibits more than an 11-fold improvement in reconstruction accuracy, but also predicts synchronous neural activity and preserves correlated low-dimensional latent dynamics. We also show that the NBGNet robustly predicts held-out data across a long time scale (2 weeks) without retraining. We further validate the effective connectivity defined from our model by demonstrating that neural connectivity during motor behaviour agrees with the established neuroanatomical hierarchy of motor control in the literature. The NBGNet approach opens the door to revealing a comprehensive understanding of brain computation, where network mechanisms of multi-scale activity are critical.

Multivariate Relationships Between Health Outcomes and Health System Performance Indicators: An Integrated Factor Analysis With Canonical Correlations

OBJECTIVES

This study aimed to investigate the relationships between sets of variables related to health system performance indicators and health outcomes.

METHODS

The relationships between a set of health outcomes and a set of health system performance indicators of a developing country were examined using multivariate statistical analysis techniques. A combinative strategy of explanatory factor analysis and the canonical correlation coefficient was used to define linear structural relationships between study variables. Province-based data were gathered from2 official statistical records of the Turkish Statistical Institute for the year 2019. Life expectancy at birth, infant mortality rate, and crude death rate were accepted as health outcome indicators.

RESULTS

The explanatory factor analysis indicated 2 independent variable groups, namely (1) health-related human resources and capacity and (2) health service utilization characteristics. The results of the canonical correlation analysis illustrated good performance to define sparse linear combinations of the 2 groups of variables. There existed strong positive correlations between health outcomes and health-related human resources and capacity indicators (rc = 0.83; P < .001) and health service utilization indicators (rc = 0.59; P < .001).

CONCLUSIONS

The results of this study support the view that there is a linear and strong positive relationship between health outcomes and health-related human resources and capacity indicators. Further studies will combine big data analytics with multivariate statistical analysis techniques by studying large health system performance data sets.

Modeling Climate Change Effects on Genetic Diversity of an Endangered Horse Breed Using Canonical Correlations

The historical increase in the occurrence of extreme weather events in Spain during the last thirty years makes it a perfect location for the evaluation of climate change. Modeling the effects of climate change on domestic animals' genetic diversity may help to anticipate challenging situations. However, animal populations' short life cycle and patent lack of historical information during extended periods of time drastically compromise the evaluation of climate change effects. Locally adapted breeds' gene pool is the base for their improved resilience and plasticity in response to climate change's extreme climatic conditions. The preservation of these domestic resources offers selection alternatives to breeders who seek such improved adaptability. The Spanish endangered autochthonous Hispano-Arabian horse breed is perfectly adapted to the conditions of the territory where it was created, developed, and widespread worldwide. The possibility to trace genetic diversity in the Hispano-Arabian breed back around seven decades and its global ubiquity make this breed an idoneous reference subject to act as a model for other international populations. Climate change's shaping effects on the genetic diversity of the Hispano-Arabian horse breed's historical population were monitored from 1950 to 2019 and evaluated. Wind speed, gust speed, or barometric pressure have greater repercussions than extreme temperatures on genetic diversity. Extreme climate conditions, rather than average modifications of climate, may push breeders/owners to implement effective strategies in the short to medium term, but the effect will be plausible in the long term due to breed sustainability and enhanced capacity of response to extreme climate events. When extreme climatic conditions occur, breeders opt for mating highly diverse unrelated individuals, avoiding the production of a large number of offspring. People in charge of domestic population conservation act as catalyzers of the regulatory changes occurring during breeds' climate change adaptive process and may identify genes conferring their animals with greater adaptability but still maintaining enhanced performance. This model assists in determining how owners of endangered domestic populations should plan their breeding strategies, seeking the obtention of animals more resilient and adapted to climate-extreme conditions. This efficient alternative is focused on the obtention of increased profitability from this population and in turn ensuring their sustainability.

The Crucial Role of Breeder and Dog Owners Associations in Safeguarding Genetic Heritage of Endangered Balearic Dog Breeds: Gender Preference and Registry Adscription

This study delves into the complex relationships between indigenous dog breeds in the Balearic Isles and their human counterparts, specifically breeders and owners. Using Canonical Correlation Analysis, the research examines variables such as breed registries and the number of breeders/owners, uncovering significant correlations within registries. For example, an increase in female auxiliary registrations corresponds to a decline in foundational registrations, indicating shifts in breed documentation dynamics. Similarly, a rise in definitive female registrations coincides with a decrease in foundational female registrations, suggesting increased pedigree awareness across generations. Beyond registries, the study explores the correlation between breeders/owners and various initial records, highlighting that a notable increase in breeders positively influences initial registrations, definitive totals, and overall counts, underscoring their crucial role in early breed stages. Gender preferences in registrations are noted, with a historical bias towards female entries during foundational stages gradually shifting in favor of males in definitive registrations. In conclusion, the research underscores the interconnected roles of breeders, owners, and comprehensive registries in preserving genetic diversity among Balearic dog breeds, emphasizing the need for ongoing efforts to address gaps in genealogical data for a more accurate understanding of breed dynamics.

A multiscale characterization of cortical shape asymmetries in early psychosis

Psychosis has often been linked to abnormal cortical asymmetry, but prior results have been inconsistent. Here, we applied a novel spectral shape analysis to characterize cortical shape asymmetries in patients with early psychosis across different spatial scales. We used the Human Connectome Project for Early Psychosis dataset (aged 16-35), comprising 56 healthy controls (37 males, 19 females) and 112 patients with early psychosis (68 males, 44 females). We quantified shape variations of each hemisphere over different spatial frequencies and applied a general linear model to compare differences between healthy controls and patients with early psychosis. We further used canonical correlation analysis to examine associations between shape asymmetries and clinical symptoms. Cortical shape asymmetries, spanning wavelengths from about 22 to 75 mm, were significantly different between healthy controls and patients with early psychosis (Cohen's d = 0.28-0.51), with patients showing greater asymmetry in cortical shape than controls. A single canonical mode linked the asymmetry measures to symptoms (canonical correlation analysis r = 0.45), such that higher cortical asymmetry was correlated with more severe excitement symptoms and less severe emotional distress. Significant group differences in the asymmetries of traditional morphological measures of cortical thickness, surface area, and gyrification, at either global or regional levels, were not identified. Cortical shape asymmetries are more sensitive than other morphological asymmetries in capturing abnormalities in patients with early psychosis. These abnormalities are expressed at coarse spatial scales and are correlated with specific symptom domains.

Static and dynamic fMRI-derived functional connectomes represent largely similar information

Functional connectivity (FC) of blood oxygen level-dependent (BOLD) fMRI time series can be estimated using methods that differ in sensitivity to the temporal order of time points (static vs. dynamic) and the number of regions considered in estimating a single edge (bivariate vs. multivariate). Previous research suggests that dynamic FC explains variability in FC fluctuations and behavior beyond static FC. Our aim was to systematically compare methods on both dimensions. We compared five FC methods: Pearson's/full correlation (static, bivariate), lagged correlation (dynamic, bivariate), partial correlation (static, multivariate), and multivariate AR model with and without self-connections (dynamic, multivariate). We compared these methods by (i) assessing similarities between FC matrices, (ii) by comparing node centrality measures, and (iii) by comparing the patterns of brain-behavior associations. Although FC estimates did not differ as a function of sensitivity to temporal order, we observed differences between the multivariate and bivariate FC methods. The dynamic FC estimates were highly correlated with the static FC estimates, especially when comparing group-level FC matrices. Similarly, there were high correlations between the patterns of brain-behavior associations obtained using the dynamic and static FC methods. We conclude that the dynamic FC estimates represent information largely similar to that of the static FC.

When Age Tips the Balance: a Dual Mechanism Affecting Hemispheric Specialization for Language

Aging engenders neuroadaptations, generally reducing specificity and selectivity in functional brain responses. Our investigation delves into the functional specialization of brain hemispheres within language-related networks across adulthood. In a cohort of 728 healthy adults spanning ages 18 to 88, we modeled the trajectories of inter-hemispheric asymmetry concerning the principal functional gradient across 37 homotopic regions of interest (hROIs) of an extensive language network, known as the Language-and-Memory Network. Our findings reveal that over two-thirds of Language-and-Memory Network hROIs undergo asymmetry changes with age, falling into two main clusters. The first cluster evolves from left-sided specialization to right-sided tendencies, while the second cluster transitions from right-sided asymmetry to left-hemisphere dominance. These reversed asymmetry shifts manifest around midlife, occurring after age 50, and are associated with poorer language production performance. Our results provide valuable insights into the influence of functional brain asymmetries on language proficiency and present a dynamic perspective on brain plasticity during the typical aging process.

Disruptions of Hierarchical Cortical Organization in Early Psychosis and Schizophrenia

BACKGROUND

The cerebral cortex is organized hierarchically along an axis that spans unimodal sensorimotor to transmodal association areas. This hierarchy is often characterized using low-dimensional embeddings, termed gradients, of interregional functional coupling estimates measured with resting-state functional magnetic resonance imaging. Such analyses may offer insights into the pathophysiology of schizophrenia, which has been frequently linked to dysfunctional interactions between association and sensorimotor areas.

METHODS

To examine disruptions of hierarchical cortical function across distinct stages of psychosis, we applied diffusion map embedding to 2 independent functional magnetic resonance imaging datasets: one comprising 114 patients with early psychosis and 48 control participants, and the other comprising 50 patients with established schizophrenia and 121 control participants. Then, we analyzed the primary sensorimotor-to-association and secondary visual-to-sensorimotor gradients of each participant in both datasets.

RESULTS

There were no significant differences in regional gradient scores between patients with early psychosis and control participants. Patients with established schizophrenia showed significant differences in the secondary, but not primary, gradient compared with control participants. Gradient differences in schizophrenia were characterized by lower within-network dispersion in the dorsal attention (false discovery rate [FDR]-corrected p [pFDR] < .001), visual (pFDR = .003), frontoparietal (pFDR = .018), and limbic (pFDR = .020) networks and lower between-network dispersion between the visual network and other networks (pFDR < .001).

CONCLUSIONS

These findings indicate that differences in cortical hierarchical function occur along the secondary visual-to-sensorimotor axis rather than the primary sensorimotor-to-association axis as previously thought. The absence of differences in early psychosis suggests that visual-sensorimotor abnormalities may emerge as the illness progresses.

Investigating the spontaneous brain activities of patients with subjective cognitive decline and mild cognitive impairment: an amplitude of low-frequency fluctuation functional magnetic resonance imaging study

Background

Subjective cognitive decline (SCD) and mild cognitive impairment (MCI) are neurodegenerative processing stages of Alzheimer's disease (AD). Cognitive decline is thought to manifest in intrinsic brain activity changes, but research results yielded conflicting and few studies have explored the roles of brain regions in cognitive decline, and sensitivity of the cognitive field to changes in the altered intrinsic brain activity.

Methods

In this cross-sectional study, 158 elderly participants were recruited from the memory clinic of the First Affiliated Hospital of Nanjing Medical University from July 2019 to May 2021, and grouped into SCD (n=73), MCI (n=46), and normal controls (NC) (n=39). The amplitude of low-frequency fluctuation (ALFF) was calculated and evaluated among the groups. Then canonical correlation analysis (CCA) was conducted to investigate the associations between imaging outcomes and cognitive behaviors.

Results

Neuropsychological tests in different cognitive dimensions and ALFF values of the prefrontal, parietal, and temporal gyrus, were significantly different (P<0.05) among the three groups, with no appreciable decline in daily activity. The changes in intrinsic activities were closely related to the decline in cognitive function (R=0.73, P=0.002). ALFF values in the left middle occipital gyrus, right middle frontal gyrus, left superior frontal gyrus, left angular gyrus, and superior temporal gyrus played significant roles in the analysis, while the Montreal Cognitive Assessment (MoCA) and Auditory-Verbal Learning Test scores were found to be more sensitive to changes in ALFF values.

Conclusions

Spontaneous brain activity is a stable imaging biomarker of cognitive impairment. ALFF changes of the prefrontal, occipital, left angular, and temporal gyrus were sensitive to identifying cognitive decline, and the scores of the Auditory-Verbal Learning Test and MoCA could predict the abnormal intrinsic activities.

Multivariate links between the developmental timing of adversity exposure and white matter tract integrity in adulthood

Early-life adversity is pervasive worldwide and represents a potent risk factor for increased mental health burden across the lifespan. However, there is substantial individual heterogeneity in associations between adversity exposure, neurobiological changes, and mental health problems. Accounting for key features of adversity such as the developmental timing of exposure may clarify associations between adversity, neurodevelopment, and mental health. The present study leverages sparse canonical correlation analysis to characterize modes of covariation between age of adversity exposure and the integrity of white matter tracts throughout the brain in a sample of 107 adults. We find that adversity exposure during middle childhood (ages 5-6 and 8-9 in particular) is consistently linked with alterations in white matter tract integrity, such that tracts supporting sensorimotor functions display higher integrity in relation to adversity exposure while tracts supporting cortico-cortical communication display lower integrity. Further, latent patterns of tract integrity linked with adversity experienced across preschool age and middle childhood (ages 4-9) were associated with trauma-related symptoms in adulthood. Our findings underscore that adversity exposure may differentially affect white matter in a function- and developmental-timing specific manner and suggest that adversity experienced between ages 4-9 may shape the development of global white matter tracts in ways that are relevant for adult mental health.

Autism Is Associated With Interindividual Variations of Gray and White Matter Morphology

BACKGROUND

Although many studies have explored atypicalities in gray matter (GM) and white matter (WM) morphology of autism, most of them relied on unimodal analyses that did not benefit from the likelihood that different imaging modalities may reflect common neurobiology. We aimed to establish brain patterns of modalities that differentiate between individuals with and without autism and explore associations between these brain patterns and clinical measures in the autism group.

METHODS

We studied 183 individuals with autism and 157 nonautistic individuals (age range, 6-30 years) in a large, deeply phenotyped autism dataset (EU-AIMS LEAP [European Autism Interventions-A Multicentre Study for Developing New Medications Longitudinal European Autism Project]). Linked independent component analysis was used to link all participants' GM volume and WM diffusion tensor images, and group comparisons of modality shared variances were examined. Subsequently, we performed univariate and multivariate brain-behavior correlation analyses to separately explore the relationships between brain patterns and clinical profiles.

RESULTS

One multimodal pattern was significantly related to autism. This pattern was primarily associated with GM volume in bilateral insula and frontal, precentral and postcentral, cingulate, and caudate areas and co-occurred with altered WM features in the superior longitudinal fasciculus. The brain-behavior correlation analyses showed a significant multivariate association primarily between brain patterns that involved variation of WM and symptoms of restricted and repetitive behavior in the autism group.

CONCLUSIONS

Our findings demonstrate the assets of integrated analyses of GM and WM alterations to study the brain mechanisms that underpin autism and show that the complex clinical autism phenotype can be interpreted by brain covariation patterns that are spread across the brain involving both cortical and subcortical areas.

Biomarkers of professional cybersportsmen: Event related potentials and cognitive tests study

"Faster, higher, stronger" is the motto of any professional athlete. Does that apply to brain dynamics as well? In our paper, we performed a series of EEG experiments on Visually Evoked Potentials and a series of cognitive tests-reaction time and visual search, with professional eSport players in Counter-Strike: Global Offensive (CS:GO) and novices (control group) in order to find important differences between them. EEG data were studied in a temporal domain by Event-Related Potentials (ERPs) and in a frequency domain by Variational Mode Decomposition. The EEG analysis showed that the brain reaction of eSport players is faster (P300 latency is earlier on average by 20-70 ms, p < 0.005) and stronger (P300 peak amplitude is higher on average by 7-9 mkV, p < 0.01). Professional eSport players also exhibit stronger stimulus-locked alpha-band power. Besides, the Spearman correlation analysis showed a significant correlation between hours spend in CS:GO and mean amplitude of P200 and N200 for the professional players. The comparison of cognitive test results showed the superiority of the professional players to the novices in reaction time (faster) and choice reaction time-faster reaction, but similar correctness, while a significant difference in visual search skills was not detected. Thus, significant differences in EEG signals (in spectrograms and ERPs) and cognitive test results (reaction time) were detected between the professional players and the control group. Cognitive tests could be used to separate skilled players from novices, while EEG testing can help to understand the skilled player's level. The results can contribute to understanding the impact of eSport on a player's cognitive state and associating eSport with a real sport. Moreover, the presented results can be useful for evaluating eSport team members and making training plans.

The Cellular Underpinnings of the Human Cortical Connectome

The functional properties of the human brain arise, in part, from the vast assortment of cell types that pattern the cortex. The cortical sheet can be broadly divided into distinct networks, which are further embedded into processing streams, or gradients, that extend from unimodal systems through higher-order association territories. Here, using transcriptional data from the Allen Human Brain Atlas, we demonstrate that imputed cell type distributions are spatially coupled to the functional organization of cortex, as estimated through fMRI. Cortical cellular profiles follow the macro-scale organization of the functional gradients as well as the associated large-scale networks. Distinct cellular fingerprints were evident across networks, and a classifier trained on post-mortem cell-type distributions was able to predict the functional network allegiance of cortical tissue samples. These data indicate that the in vivo organization of the cortical sheet is reflected in the spatial variability of its cellular composition.

Static and dynamic functional connectomes represent largely similar information

Functional connectivity (FC) of blood-oxygen-level-dependent (BOLD) fMRI time series can be estimated using methods that differ in sensitivity to the temporal order of time points (static vs. dynamic) and the number of regions considered in estimating a single edge (bivariate vs. multivariate). Previous research suggests that dynamic FC explains variability in FC fluctuations and behavior beyond static FC. Our aim was to systematically compare methods on both dimensions. We compared five FC methods: Pearson's/full correlation (static, bivariate), lagged correlation (dynamic, bivariate), partial correlation (static, multivariate) and multivariate AR model with and without self-connections (dynamic, multivariate). We compared these methods by (i) assessing similarities between FC matrices, (ii) by comparing node centrality measures, and (iii) by comparing the patterns of brain-behavior associations. Although FC estimates did not differ as a function of sensitivity to temporal order, we observed differences between the multivariate and bivariate FC methods. The dynamic FC estimates were highly correlated with the static FC estimates, especially when comparing group-level FC matrices. Similarly, there were high correlations between the patterns of brain-behavior associations obtained using the dynamic and static FC methods. We conclude that the dynamic FC estimates represent information largely similar to that of the static FC.

Dysconnection and cognition in schizophrenia: A spectral dynamic causal modeling study

Schizophrenia (SZ) is a severe mental disorder characterized by failure of functional integration (aka dysconnection) across the brain. Recent functional connectivity (FC) studies have adopted functional parcellations to define subnetworks of large-scale networks, and to characterize the (dys)connection between them, in normal and clinical populations. While FC examines statistical dependencies between observations, model-based effective connectivity (EC) can disclose the causal influences that underwrite the observed dependencies. In this study, we investigated resting state EC within seven large-scale networks, in 66 SZ and 74 healthy subjects from a public dataset. The results showed that a remarkable 33% of the effective connections (among subnetworks) of the cognitive control network had been pathologically modulated in SZ. Further dysconnection was identified within the visual, default mode and sensorimotor networks of SZ subjects, with 24%, 20%, and 11% aberrant couplings. Overall, the proportion of discriminative connections was remarkably larger in EC (24%) than FC (1%) analysis. Subsequently, to study the neural correlates of impaired cognition in SZ, we conducted a canonical correlation analysis between the EC parameters and the cognitive scores of the patients. As such, the self-inhibitions of supplementary motor area and paracentral lobule (in the sensorimotor network) and the excitatory connection from parahippocampal gyrus to inferior temporal gyrus (in the cognitive control network) were significantly correlated with the social cognition, reasoning/problem solving and working memory capabilities of the patients. Future research can investigate the potential of whole-brain EC as a biomarker for diagnosis of brain disorders and for neuroimaging-based cognitive assessment.

Overlapping brain correlates of superior cognition among children at genetic risk for Alzheimer's disease and/or major depressive disorder

Early life adversity (ELA) tends to accelerate neurobiological ageing, which, in turn, is thought to heighten vulnerability to both major depressive disorder (MDD) and Alzheimer's disease (AD). The two conditions are putatively related, with MDD representing either a risk factor or early symptom of AD. Given the substantial environmental susceptibility of both disorders, timely identification of their neurocognitive markers could facilitate interventions to prevent clinical onset. To this end, we analysed multimodal data from the Adolescent Brain and Cognitive Development study (ages 9-10 years). To disentangle genetic from correlated genetic-environmental influences, while also probing gene-adversity interactions, we compared adoptees, a group generally exposed to substantial ELA, with children raised by their biological families via genetic risk scores (GRS) from genome-wide association studies. AD and MDD GRSs predicted overlapping and widespread neurodevelopmental alterations associated with superior fluid cognition. Specifically, among adoptees only, greater AD GRS were related to accelerated structural maturation (i.e., cortical thinning) and higher MDD GRS were linked to delayed functional neurodevelopment, as reflected in compensatory brain activation on an inhibitory control task. Our study identifies compensatory mechanisms linked to MDD risk and highlights the potential cognitive benefits of accelerated maturation linked to AD vulnerability in late childhood.

Serotonin regulation of behavior via large-scale neuromodulation of serotonin receptor networks

Although we understand how serotonin receptors function at the single-cell level, what role different serotonin receptors play in regulating brain-wide activity and, in turn, human behavior, remains unknown. Here, we developed transcriptomic-neuroimaging mapping to characterize brain-wide functional signatures associated with specific serotonin receptors: serotonin receptor networks (SRNs). Probing SRNs with optogenetics-functional magnetic resonance imaging (MRI) and pharmacology in mice, we show that activation of dorsal raphe serotonin neurons differentially modulates the amplitude and functional connectivity of different SRNs, showing that receptors' spatial distributions can confer specificity not only at the local, but also at the brain-wide, network level. In humans, using resting-state functional MRI, SRNs replicate established divisions of serotonin effects on impulsivity and negative biases. These results provide compelling evidence that heterogeneous brain-wide distributions of different serotonin receptor types may underpin behaviorally distinct modes of serotonin regulation. This suggests that serotonin neurons may regulate multiple aspects of human behavior via modulation of large-scale receptor networks.

Multimodal and multidomain lesion network mapping enhances prediction of sensorimotor behavior in stroke patients

Beyond the characteristics of a brain lesion, such as its etiology, size or location, lesion network mapping (LNM) has shown that similar symptoms after a lesion reflects similar dis-connectivity patterns, thereby linking symptoms to brain networks. Here, we extend LNM by using a multimodal strategy, combining functional and structural networks from 1000 healthy participants in the Human Connectome Project. We apply multimodal LNM to a cohort of 54 stroke patients with the aim of predicting sensorimotor behavior, as assessed through a combination of motor and sensory tests. Results are two-fold. First, multimodal LNM reveals that the functional modality contributes more than the structural one in the prediction of sensorimotor behavior. Second, when looking at each modality individually, the performance of the structural networks strongly depended on whether sensorimotor performance was corrected for lesion size, thereby eliminating the effect that larger lesions generally produce more severe sensorimotor impairment. In contrast, functional networks provided similar performance regardless of whether or not the effect of lesion size was removed. Overall, these results support the extension of LNM to its multimodal form, highlighting the synergistic and additive nature of different types of network modalities, and their corresponding influence on behavioral performance after brain injury.

Canonical Correlation Analysis and Partial Least Squares for Identifying Brain-Behavior Associations: A Tutorial and a Comparative Study

Canonical correlation analysis (CCA) and partial least squares (PLS) are powerful multivariate methods for capturing associations across 2 modalities of data (e.g., brain and behavior). However, when the sample size is similar to or smaller than the number of variables in the data, standard CCA and PLS models may overfit, i.e., find spurious associations that generalize poorly to new data. Dimensionality reduction and regularized extensions of CCA and PLS have been proposed to address this problem, yet most studies using these approaches have some limitations. This work gives a theoretical and practical introduction into the most common CCA/PLS models and their regularized variants. We examine the limitations of standard CCA and PLS when the sample size is similar to or smaller than the number of variables. We discuss how dimensionality reduction and regularization techniques address this problem and explain their main advantages and disadvantages. We highlight crucial aspects of the CCA/PLS analysis framework, including optimizing the hyperparameters of the model and testing the identified associations for statistical significance. We apply the described CCA/PLS models to simulated data and real data from the Human Connectome Project and Alzheimer's Disease Neuroimaging Initiative (both of n > 500). We use both low- and high-dimensionality versions of these data (i.e., ratios between sample size and variables in the range of ∼1-10 and ∼0.1-0.01, respectively) to demonstrate the impact of data dimensionality on the models. Finally, we summarize the key lessons of the tutorial.

Shared brain and genetic architectures between mental health and physical activity

Physical activity is correlated with, and effectively treats various forms of psychopathology. However, whether biological correlates of physical activity and psychopathology are shared remains unclear. Here, we examined the extent to which the neural and genetic architecture of physical activity and mental health are shared. Using data from the UK Biobank (N = 6389), we applied canonical correlation analysis to estimate associations between the amplitude and connectivity strength of subnetworks of three major neurocognitive networks (default mode, DMN; salience, SN; central executive networks, CEN) with accelerometer-derived measures of physical activity and self-reported mental health measures (primarily of depression, anxiety disorders, neuroticism, subjective well-being, and risk-taking behaviors). We estimated the genetic correlation between mental health and physical activity measures, as well as putative causal relationships by applying linkage disequilibrium score regression, genomic structural equational modeling, and latent causal variable analysis to genome-wide association summary statistics (GWAS N = 91,105-500,199). Physical activity and mental health were associated with connectivity strength and amplitude of the DMN, SN, and CEN (r's ≥ 0.12, p's < 0.048). These neural correlates exhibited highly similar loading patterns across mental health and physical activity models even when accounting for their shared variance. This suggests a largely shared brain network architecture between mental health and physical activity. Mental health and physical activity (including sleep) were also genetically correlated (|rg| = 0.085-0.121), but we found no evidence for causal relationships between them. Collectively, our findings provide empirical evidence that mental health and physical activity have shared brain and genetic architectures and suggest potential candidate subnetworks for future studies on brain mechanisms underlying beneficial effects of physical activity on mental health.

Canonical correlation analysis as a statistical method to relate underwater acoustic propagation and ocean fluctuations

Numerical models are currently used to understand how environmental fluctuations impact acoustic propagation. Such a process can be tedious in complex fluctuating environments. This letter proposes a complementary approach based upon canonical correlation analysis (CCA) to determine statistical relationships between two sets of observed acoustic and oceanographic variables. It is shown, as an example, how CCA puts forward the impact of external and internal tide on shallow water propagation. Results are consistent with the physical understanding of tide impact on acoustic propagation. They encourage the use of CCA for complex studies, in particular, for environments fluctuating under several environmental phenomena.

Network hub centrality and working memory performance in schizophrenia

Cognitive impairment, and working memory deficits in particular, are debilitating, treatment-resistant aspects of schizophrenia. Dysfunction of brain network hubs, putatively related to altered neurodevelopment, is thought to underlie the cognitive symptoms associated with this illness. Here, we used weighted degree, a robust graph theory metric representing the number of weighted connections to a node, to quantify centrality in cortical hubs in 29 patients with schizophrenia and 29 age- and gender-matched healthy controls and identify the critical nodes that underlie working memory performance. In both patients and controls, elevated weighted degree in the default mode network (DMN) was generally associated with poorer performance (accuracy and reaction time). Higher degree in the ventral attention network (VAN) nodes in the right superior temporal cortex was associated with better performance (accuracy) in patients. Degree in several prefrontal and parietal areas was associated with cognitive performance only in patients. In regions that are critical for sustained attention, these correlations were primarily driven by between-network connectivity in patients. Moreover, a cross-validated prediction analysis showed that a linear model using a summary degree score can be used to predict an individual’s working memory accuracy (r = 0.35). Our results suggest that schizophrenia is associated with dysfunctional hubs in the cortical systems supporting internal and external cognition and highlight the importance of topological network analysis in the search of biomarkers for cognitive deficits in schizophrenia.

Psychological Resilience and Neurodegenerative Risk: A Connectomics-Transcriptomics Investigation in Healthy Adolescent and Middle-Aged Females

Adverse life events can inflict substantial long-term damage, which, paradoxically, has been posited to stem from initially adaptative responses to the challenges encountered in one's environment. Thus, identification of the mechanisms linking resilience against recent stressors to longer-term psychological vulnerability is key to understanding optimal functioning across multiple timescales. To address this issue, our study tested the relevance of neuro-reproductive maturation and senescence, respectively, to both resilience and longer-term risk for pathologies characterised by accelerated brain aging, specifically, Alzheimer's Disease (AD). Graph theoretical and partial least squares analyses were conducted on multimodal imaging, reported biological aging and recent adverse experience data from the Lifespan Human Connectome Project (HCP). Availability of reproductive maturation/senescence measures restricted our investigation to adolescent (N =178) and middle-aged (N=146) females. Psychological resilience was linked to age-specific brain senescence patterns suggestive of precocious functional development of somatomotor and control-relevant networks (adolescence) and earlier aging of default mode and salience/ventral attention systems (middle adulthood). Biological aging showed complementary associations with the neural patterns relevant to resilience in adolescence (positive relationship) versus middle-age (negative relationship). Transcriptomic and expression quantitative trait locus data analyses linked the neural aging patterns correlated with psychological resilience in middle adulthood to gene expression patterns suggestive of increased AD risk. Our results imply a partially antagonistic relationship between resilience against proximal stressors and longer-term psychological adjustment in later life. They thus underscore the importance of fine-tuning extant views on successful coping by considering the multiple timescales across which age-specific processes may unfold.

The individuality of shape asymmetries of the human cerebral cortex

Asymmetries of the cerebral cortex are found across diverse phyla and are particularly pronounced in humans, with important implications for brain function and disease. However, many prior studies have confounded asymmetries due to size with those due to shape. Here, we introduce a novel approach to characterize asymmetries of the whole cortical shape, independent of size, across different spatial frequencies using magnetic resonance imaging data in three independent datasets. We find that cortical shape asymmetry is highly individualized and robust, akin to a cortical fingerprint, and identifies individuals more accurately than size-based descriptors, such as cortical thickness and surface area, or measures of inter-regional functional coupling of brain activity. Individual identifiability is optimal at coarse spatial scales (~37 mm wavelength), and shape asymmetries show scale-specific associations with sex and cognition, but not handedness. While unihemispheric cortical shape shows significant heritability at coarse scales (~65 mm wavelength), shape asymmetries are determined primarily by subject-specific environmental effects. Thus, coarse-scale shape asymmetries are highly personalized, sexually dimorphic, linked to individual differences in cognition, and are primarily driven by stochastic environmental influences.

A hierarchical Bayesian model to find brain-behaviour associations in incomplete data sets

Canonical Correlation Analysis (CCA) and its regularised versions have been widely used in the neuroimaging community to uncover multivariate associations between two data modalities (e.g., brain imaging and behaviour). However, these methods have inherent limitations: (1) statistical inferences about the associations are often not robust; (2) the associations within each data modality are not modelled; (3) missing values need to be imputed or removed. Group Factor Analysis (GFA) is a hierarchical model that addresses the first two limitations by providing Bayesian inference and modelling modality-specific associations. Here, we propose an extension of GFA that handles missing data, and highlight that GFA can be used as a predictive model. We applied GFA to synthetic and real data consisting of brain connectivity and non-imaging measures from the Human Connectome Project (HCP). In synthetic data, GFA uncovered the underlying shared and specific factors and predicted correctly the non-observed data modalities in complete and incomplete data sets. In the HCP data, we identified four relevant shared factors, capturing associations between mood, alcohol and drug use, cognition, demographics and psychopathological measures and the default mode, frontoparietal control, dorsal and ventral networks and insula, as well as two factors describing associations within brain connectivity. In addition, GFA predicted a set of non-imaging measures from brain connectivity. These findings were consistent in complete and incomplete data sets, and replicated previous findings in the literature. GFA is a promising tool that can be used to uncover associations between and within multiple data modalities in benchmark datasets (such as, HCP), and easily extended to more complex models to solve more challenging tasks.

Functional Connectivity in Antipsychotic-Treated and Antipsychotic-Naive Patients With First-Episode Psychosis and Low Risk of Self-harm or Aggression: A Secondary Analysis of a Randomized Clinical Trial

IMPORTANCE

Altered functional connectivity (FC) is a common finding in resting-state functional magnetic resonance imaging (rs-fMRI) studies of people with psychosis, yet how FC disturbances evolve in the early stages of illness, and how antipsychotic treatment influences these disturbances, remains unknown.

OBJECTIVE

To investigate longitudinal FC changes in antipsychotic-naive and antipsychotic-treated patients with first-episode psychosis (FEP).

DESIGN, SETTING, AND PARTICIPANTS

This secondary analysis of a triple-blind, randomized clinical trial was conducted over a 5-year recruitment period between April 2008 and December 2016 with 59 antipsychotic-naive patients with FEP receiving either a second-generation antipsychotic or a placebo pill over a treatment period of 6 months. Participants were required to have low suicidality and aggression, to have a duration of untreated psychosis of less than 6 months, and to be living in stable accommodations with social support. Both FEP groups received intensive psychosocial therapy. A healthy control group was also recruited. Participants completed rs-fMRI scans at baseline, 3 months, and 12 months. Data were analyzed from May 2019 to August 2020.

INTERVENTIONS

Resting-state functional MRI was used to probe brain FC. Patients received either a second-generation antipsychotic or a matched placebo tablet. Both patient groups received a manualized psychosocial intervention.

MAIN OUTCOMES AND MEASURES

The primary outcomes of this analysis were to investigate (1) FC differences between patients and controls at baseline; (2) FC changes in medicated and unmedicated patients between baseline and 3 months; and (3) associations between longitudinal FC changes and clinical outcomes. An additional aim was to investigate long-term FC changes at 12 months after baseline. These outcomes were not preregistered.

RESULTS

Data were analyzed for 59 patients (antipsychotic medication plus psychosocial treatment: 28 [47.5%]; mean [SD] age, 19.5 [3.0] years; 15 men [53.6%]; placebo plus psychosocial treatment: 31 [52.5%]; mean [SD] age, 18.8 [2.7]; 16 men [51.6%]) and 27 control individuals (mean [SD] age, 21.9 [1.9] years). At baseline, patients showed widespread functional dysconnectivity compared with controls, with reductions predominantly affecting interactions between the default mode network, limbic systems, and the rest of the brain. From baseline to 3 months, patients receiving placebo showed increased FC principally within the same systems; some of these changes correlated with improved clinical outcomes (canonical correlation analysis R = 0.901; familywise error-corrected P = .005). Antipsychotic exposure was associated with increased FC primarily between the thalamus and the rest of the brain.

CONCLUSIONS AND RELEVANCE

In this secondary analysis of a clinical trial, antipsychotic-naive patients with FEP showed widespread functional dysconnectivity at baseline, followed by an early normalization of default mode network and cortical limbic dysfunction in patients receiving placebo and psychosocial intervention. Antipsychotic exposure was associated with FC changes concentrated on thalamocortical networks.

TRIAL REGISTRATION

ACTRN12607000608460.

Interindividual Signatures of fMRI Temporal Fluctuations

The complexity and variability of human brain activity, such as quantified from Functional Magnetic Resonance Imaging (fMRI) time series, have been widely studied as potential markers of healthy and pathological states. However, the extent to which fMRI temporal features exhibit stable markers of inter-individual differences in brain function across healthy young adults is currently an open question. In this study, we draw upon two widely used time-series measures-a nonlinear complexity measure (sample entropy; SampEn) and a spectral measure of low-frequency content (fALFF)-to capture dynamic properties of resting-state fMRI in a large sample of young adults from the Human Connectome Project. We observe that these two measures are closely related, and that both generate reproducible patterns across brain regions over four different fMRI runs, with intra-class correlations of up to 0.8. Moreover, we find that both metrics can uniquely differentiate subjects with high identification rates (ca. 89%). Canonical correlation analysis revealed a significant relationship between multivariate brain temporal features and behavioral measures. Overall, these findings suggest that regional profiles of fMRI temporal characteristics may provide stable markers of individual differences, and motivate future studies to further probe relationships between fMRI time series metrics and behavior.

Baseline Functional Connectivity in Resting State Networks Associated with Depression and Remission Status after 16 Weeks of Pharmacotherapy: A CAN-BIND Report

Understanding the neural underpinnings of major depressive disorder (MDD) and its treatment could improve treatment outcomes. So far, findings are variable and large sample replications scarce. We aimed to replicate and extend altered functional connectivity associated with MDD and pharmacotherapy outcomes in a large, multisite sample. Resting-state fMRI data were collected from 129 patients and 99 controls through the Canadian Biomarker Integration Network in Depression. Symptoms were assessed with the Montgomery-Åsberg Depression Rating Scale (MADRS). Connectivity was measured as correlations between four seeds (anterior and posterior cingulate cortex, insula and dorsolateral prefrontal cortex) and all other brain voxels. Partial least squares was used to compare connectivity prior to treatment between patients and controls, and between patients reaching remission (MADRS ≤ 10) early (within 8 weeks), late (within 16 weeks), or not at all. We replicated previous findings of altered connectivity in patients. In addition, baseline connectivity of the anterior/posterior cingulate and insula seeds differentiated patients with different treatment outcomes. The stability of these differences was established in the largest single-site subsample. Our replication and extension of altered connectivity highlighted previously reported and new differences between patients and controls, and revealed features that might predict remission prior to pharmacotherapy. Trial registration:  ClinicalTrials.gov: NCT01655706.

Shared and Anxiety-Specific Pediatric Psychopathology Dimensions Manifest Distributed Neural Correlates

BACKGROUND

Imaging research has not yet delivered reliable psychiatric biomarkers. One challenge, particularly among youth, is high comorbidity. This challenge might be met through canonical correlation analysis designed to model mutual dependencies between symptom dimensions and neural measures. We mapped the multivariate associations that intrinsic functional connectivity manifests with pediatric symptoms of anxiety, irritability, and attention-deficit/hyperactivity disorder (ADHD) as common, impactful, co-occurring problems. We evaluate the replicability of such latent dimensions in an independent sample.

METHODS

We obtained ratings of anxiety, irritability, and ADHD, and 10 minutes of resting-state functional magnetic resonance imaging data, from two independent cohorts. Both cohorts (discovery: n = 182; replication: n = 326) included treatment-seeking youth with anxiety disorders, with disruptive mood dysregulation disorder, with ADHD, or without psychopathology. Functional connectivity was modeled as partial correlations among 216 brain areas. Using canonical correlation analysis and independent component analysis jointly we sought maximally correlated, maximally interpretable latent dimensions of brain connectivity and clinical symptoms.

RESULTS

We identified seven canonical variates in the discovery and five in the replication cohort. Of these canonical variates, three exhibited similarities across datasets: two variates consistently captured shared aspects of irritability, ADHD, and anxiety, while the third was specific to anxiety. Across cohorts, canonical variates did not relate to specific resting-state networks but comprised edges interconnecting established networks within and across both hemispheres.

CONCLUSIONS

Findings revealed two replicable types of clinical variates, one related to multiple symptom dimensions and a second relatively specific to anxiety. Both types involved a multitude of broadly distributed, weak brain connections as opposed to strong connections encompassing known resting-state networks.

Multimodal Imaging Brain Markers in Early Adolescence Are Linked with a Physically Active Lifestyle

The World Health Organization promotes physical exercise and a healthy lifestyle as means to improve youth development. However, relationships between physical lifestyle and human brain development are not fully understood. Here, we asked whether a human brain-physical latent mode of covariation underpins the relationship between physical activity, fitness, and physical health measures with multimodal neuroimaging markers. In 50 12-year old school pupils (26 females), we acquired multimodal whole-brain MRI, characterizing brain structure, microstructure, function, myelin content, and blood perfusion. We also acquired physical variables measuring objective fitness levels, 7 d physical activity, body mass index, heart rate, and blood pressure. Using canonical correlation analysis, we unravel a latent mode of brain-physical covariation, independent of demographics, school, or socioeconomic status. We show that MRI metrics with greater involvement in this mode also showed spatially extended patterns across the brain. Specifically, global patterns of greater gray matter perfusion, volume, cortical surface area, greater white matter extra-neurite density, and resting state networks activity covaried positively with measures reflecting a physically active phenotype (high fit, low sedentary individuals). Showing that a physically active lifestyle is linked with systems-level brain MRI metrics, these results suggest widespread associations relating to several biological processes. These results support the notion of close brain-body relationships and underline the importance of investigating modifiable lifestyle factors not only for physical health but also for brain health early in adolescence.SIGNIFICANCE STATEMENT An active lifestyle is key for healthy development. In this work, we answer the following question: How do brain neuroimaging markers relate with young adolescents' level of physical activity, fitness, and physical health? Combining advanced whole-brain multimodal MRI metrics with computational approaches, we show a robust relationship between physically active lifestyles and spatially extended, multimodal brain imaging-derived phenotypes. Suggesting a wider effect on brain neuroimaging metrics than previously thought, this work underlies the importance of studying physical lifestyle, as well as other brain-body relationships in an effort to foster brain health at this crucial stage in development.