Widespread Autonomic Physiological Coupling Across the Brain-Body Axis.

The brain is closely attuned to visceral signals from the body's internal environment, as evidenced by the numerous associations between neural, hemodynamic, and peripheral physiological signals. We show that these brain-body co-fluctuations can be captured by a single spatiotemporal pattern. Across several independent samples, as well as single-echo and multi-echo fMRI data acquisition sequences, we identify widespread co-fluctuations in the low-frequency range (0.01 - 0.1 Hz) between resting-state global fMRI signals, neural activity, and a host of autonomic signals spanning cardiovascular, pulmonary, exocrine and smooth muscle systems. The same brain-body co-fluctuations observed at rest are elicited by arousal induced by cued deep breathing and intermittent sensory stimuli, as well as spontaneous phasic EEG events during sleep. Further, we show that the spatial structure of global fMRI signals is maintained under experimental suppression of end-tidal carbon dioxide (PETCO2) variations, suggesting that respiratory-driven fluctuations in arterial CO2 accompanying arousal cannot explain the origin of these signals in the brain. These findings establish the global fMRI signal as a significant component of the arousal response governed by the autonomic nervous system.

Multimodal neural correlates of childhood psychopathology.

Complex structural and functional changes occurring in typical and atypical development necessitate multidimensional approaches to better understand the risk of developing psychopathology. Here, we simultaneously examined structural and functional brain network patterns in relation to dimensions of psychopathology in the Adolescent Brain Cognitive Development dataset. Several components were identified, recapitulating the psychopathology hierarchy, with the general psychopathology ( p ) factor explaining most covariance with multimodal imaging features, while the internalizing, externalizing, and neurodevelopmental dimensions were each associated with distinct morphological and functional connectivity signatures. Connectivity signatures associated with the p factor and neurodevelopmental dimensions followed the sensory-to-transmodal axis of cortical organization, which is related to the emergence of complex cognition and risk for psychopathology. Results were consistent in two separate data subsamples, supporting generalizability, and robust to variations in analytical parameters. Our findings help in better understanding biological mechanisms underpinning dimensions of psychopathology, and could provide brain-based vulnerability markers.

Longitudinal changes in brain asymmetry track lifestyle and disease.

Human beings may have evolved the largest asymmetries of brain organization in the animal kingdom. Hemispheric left-vs-right specialization is especially pronounced in our species-unique capacities. Yet, brain asymmetry features appear to be strongly shaped by non-genetic influences. We hence charted the largest longitudinal brain-imaging adult resource, yielding evidence that brain asymmetry changes continuously in a manner suggestive of neural plasticity. In the UK Biobank population cohort, we demonstrate that asymmetry changes show robust associations across 959 distinct phenotypic variables spanning 11 categories. We also find that changes in brain asymmetry over years co-occur with changes among specific lifestyle markers. Finally, we reveal relevance of brain asymmetry changes to major disease categories across thousands of medical diagnoses. Our results challenge the tacit assumption that asymmetrical neural systems are highly conserved throughout adulthood.

Ventral attention network connectivity is linked to cortical maturation and cognitive ability in childhood.

The human brain experiences functional changes through childhood and adolescence, shifting from an organizational framework anchored within sensorimotor and visual regions into one that is balanced through interactions with later-maturing aspects of association cortex. Here, we link this profile of functional reorganization to the development of ventral attention network connectivity across independent datasets. We demonstrate that maturational changes in cortical organization link preferentially to within-network connectivity and heightened degree centrality in the ventral attention network, whereas connectivity within network-linked vertices predicts cognitive ability. This connectivity is associated closely with maturational refinement of cortical organization. Children with low ventral attention network connectivity exhibit adolescent-like topographical profiles, suggesting that attentional systems may be relevant in understanding how brain functions are refined across development. These data suggest a role for attention networks in supporting age-dependent shifts in cortical organization and cognition across childhood and adolescence.

Intrinsic structural covariation links cerebellum subregions to the cerebral cortex.

The human cerebellum is increasingly recognized to be involved in non-motor and higher-order cognitive functions. Yet, its ties with the entire cerebral cortex have not been holistically studied in a whole-brain exploration with a unified analytical framework. Here, we characterized dissociable cortical-cerebellar structural covariation patterns based on regional gray matter volume (GMV) across the brain in n=38,527 UK Biobank participants. Our results invigorate previous observations in that important shares of cortical-cerebellar structural covariation are described as i) a dissociation between the higher-level cognitive system and lower-level sensorimotor system, as well as ii) an anticorrelation between the visual-attention system and advanced associative networks within the cerebellum. We also discovered a novel pattern of ipsilateral, rather than contralateral, cerebral-cerebellar associations. Furthermore, phenome-wide association assays revealed key phenotypes, including cognitive phenotypes, lifestyle, physical properties, and blood assays, associated with each decomposed covariation pattern, helping to understand their real-world implications. This systems neuroscience view paves the way for future studies to explore the implications of these structural covariations, potentially illuminating new pathways in our understanding of neurological and cognitive disorders.

Differential spatial working memory-related functional network reconfiguration in young and older adults.

Functional brain networks have preserved architectures in rest and task; nevertheless, previous work consistently demonstrated task-related brain functional reorganization. Efficient rest-to-task functional network reconfiguration is associated with better cognition in young adults. However, aging and cognitive load effects, as well as contributions of intra- and internetwork reconfiguration, remain unclear. We assessed age-related and load-dependent effects on global and network-specific functional reconfiguration between rest and a spatial working memory (SWM) task in young and older adults, then investigated associations between functional reconfiguration and SWM across loads and age groups. Overall, global and network-level functional reconfiguration between rest and task increased with age and load. Importantly, more efficient functional reconfiguration associated with better performance across age groups. However, older adults relied more on internetwork reconfiguration of higher cognitive and task-relevant networks. These reflect the consistent importance of efficient network updating despite recruitment of additional functional networks to offset reduction in neural resources and a change in brain functional topology in older adults. Our findings generalize the association between efficient functional reconfiguration and cognition to aging and demonstrate distinct brain functional reconfiguration patterns associated with SWM in aging, highlighting the importance of combining rest and task measures to study aging cognition.

Brain networks identified by functional connectivity (FC) have preserved architectures from rest to task and across task demands. Higher similarity, implying more efficient network reconfiguration, was associated with better cognition and task performance in young adults. To examine how it may be influenced by aging, we compared whole-brain and network-level FC similarities between resting-state and spatial working memory fMRI in young and older adults. At whole-brain level and higher order cognitive networks, older adults evidenced less efficient network reconfiguration from rest to task than young adults. Importantly, more efficient reconfiguration was associated with better accuracy. This relationship relied more on internetwork connections in older adults. Despite reduced neural resources compared to young, maintaining efficient network updating still contributes to better cognition at older age.

A network correspondence toolbox for quantitative evaluation of novel neuroimaging results.

Decades of neuroscience research has shown that macroscale brain dynamics can be reliably decomposed into a subset of large-scale functional networks, but the specific spatial topographies of these networks and the names used to describe them can vary across studies. Such discordance has hampered interpretation and convergence of research findings across the field. To address this problem, we have developed the Network Correspondence Toolbox (NCT) to permit researchers to examine and report spatial correspondence between their novel neuroimaging results and sixteen widely used functional brain atlases, consistent with recommended reporting standards developed by the Organization for Human Brain Mapping. The atlases included in the toolbox show some topographical convergence for specific networks, such as those labeled as default or visual. Network naming varies across atlases, particularly for networks spanning frontoparietal association cortices. For this reason, quantitative comparison with multiple atlases is recommended to benchmark novel neuroimaging findings. We provide several exemplar demonstrations using the Human Connectome Project task fMRI results and UK Biobank independent component analysis maps to illustrate how researchers can use the NCT to report their own findings through quantitative evaluation against multiple published atlases. The NCT provides a convenient means for computing Dice coefficients with spin test permutations to determine the magnitude and statistical significance of correspondence among user-defined maps and existing atlas labels. The NCT also includes functionality to incorporate additional atlases in the future. The adoption of the NCT will make it easier for network neuroscience researchers to report their findings in a standardized manner, thus aiding reproducibility and facilitating comparisons between studies to produce interdisciplinary insights.

High-effect gene-coding variants impact cognition, mental well-being, and neighborhood safety substrates in brain morphology.

Our genetic makeup, together with environmental and social influences, shape our brain's development. Yet, the imaging genetics field has struggled to integrate all these modalities to investigate the interplay between genetic blueprint, environment, human health, daily living skills and outcomes. Hence, we interrogated the Adolescent Brain Cognitive Development (ABCD) cohort to outline the effects of rare high-effect genetic variants on brain architecture and corresponding implications on cognitive, behavioral, psychosocial, and socioeconomic traits. Specifically, we designed a holistic pattern-learning algorithm that quantitatively dissects the impacts of copy number variations (CNVs) on brain structure and 962 behavioral variables spanning 20 categories in 7,657 adolescents. Our results reveal associations between genetic alterations, higher-order brain networks, and specific parameters of the family well-being (increased parental and child stress, anxiety and depression) or neighborhood dynamics (decreased safety); effects extending beyond the impairment of cognitive ability or language capacity, dominantly reported in the CNV literature. Our investigation thus spotlights a far-reaching interplay between genetic variation and subjective life quality in adolescents and their families.

Lower motor performance is linked with poor sleep quality, depressive symptoms, and grey matter volume alterations.

Motor performance (MP) is essential for functional independence and well-being, particularly in later life. However, the relationship between behavioural aspects such as sleep quality and depressive symptoms, which contribute to MP, and the underlying structural brain substrates of their interplay remains unclear. This study used three population-based cohorts of younger and older adults (n=1,950) from the Human Connectome Project-Young Adult (HCP-YA), HCP-Aging (HCP-A), and enhanced Nathan Kline Institute-Rockland sample (eNKI-RS). Several canonical correlation analyses were computed within a machine learning framework to assess the associations between each of the three domains (sleep quality, depressive symptoms, grey matter volume (GMV)) and MP. The HCP-YA analyses showed progressively stronger associations between MP and each domain: depressive symptoms (unexpectedly positive, r=0.13, SD=0.06), sleep quality (r=0.17, SD=0.05), and GMV (r=0.19, SD=0.06). Combining sleep and depressive symptoms significantly improved the canonical correlations (r=0.25, SD=0.05), while the addition of GMV exhibited no further increase (r=0.23, SD=0.06). In young adults, better sleep quality, mild depressive symptoms, and GMV of several brain regions were associated with better MP. This was conceptually replicated in young adults from the eNKI-RS cohort. In HCP-Aging, better sleep quality, fewer depressive symptoms, and increased GMV were associated with MP. Robust multivariate associations were observed between sleep quality, depressive symptoms and GMV with MP, as well as age-related variations in these factors. Future studies should further explore these associations and consider interventions targeting sleep and mental health to test the potential effects on MP across the lifespan.

Sex differences in the functional network underpinnings of psychotic-like experiences in children.

Psychotic-like experiences (PLEs) include a range of sub-threshold symptoms that resemble aspects of psychosis but do not necessarily indicate the presence of psychiatric illness. These experiences are highly prevalent in youth and are associated with developmental disruptions across social, academic, and emotional domains. While not all youth who report PLEs develop psychosis, many develop other psychiatric illnesses during adolescence and adulthood. As such, PLEs are theorized to represent early markers of poor mental health. Here, we characterized the similarities and differences in the neurobiological underpinnings of childhood PLEs across the sexes using a large sample from the ABCD Study (n=5,260), revealing sex-specific associations between functional networks connectivity and PLEs. We find that although the networks associated with PLEs overlap to some extent across the sexes, there are also crucial differences. In females, PLEs are associated with dispersed cortical and non-cortical connections, whereas in males, they are primarily associated with functional connections within limbic, temporal parietal, somato/motor, and visual networks. These results suggest that early transdiagnostic markers of psychopathology may be distinct across the sexes, further emphasizing the need to consider sex in psychiatric research as well as clinical practice.

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.

Protective role of parenthood on age-related brain function in mid- to late-life.

The experience of parenthood can profoundly alter one's body, mind, and environment, yet we know little about the long-term associations between parenthood and brain function and aging in adulthood. Here, we investigate the link between number of children parented (parity) and age on brain function in 19,964 females and 17,607 males from the UK Biobank. In both females and males, increased parity was positively associated with functional connectivity, particularly within the somato/motor network. Critically, the spatial topography of parity-linked effects was inversely correlated with the impact of age on functional connectivity across the brain for both females and males, suggesting that a higher number of children is associated with patterns of brain function in the opposite direction to age-related alterations. These results indicate that the changes accompanying parenthood may confer benefits to brain health across the lifespan, highlighting the importance of future work to understand the associated mechanisms.

A multidimensional investigation of sleep and biopsychosocial profiles with associated neural signatures.

Sleep is essential for optimal functioning and health. Interconnected to multiple biological, psychological and socio-environmental factors (i.e., biopsychosocial factors), the multidimensional nature of sleep is rarely capitalized on in research. Here, we deployed a data-driven approach to identify sleep-biopsychosocial profiles that linked self-reported sleep patterns to inter-individual variability in health, cognition, and lifestyle factors in 770 healthy young adults. We uncovered five profiles, including two profiles reflecting general psychopathology associated with either reports of general poor sleep or an absence of sleep complaints (i.e., sleep resilience) respectively. The three other profiles were driven by the use of sleep aids and social satisfaction, sleep duration and cognitive performance, and sleep disturbance linked to cognition and mental health. Furthermore, identified sleep-biopsychosocial profiles displayed unique patterns of brain network organization. In particular, somatomotor network connectivity alterations were involved in the relationships between sleep and biopsychosocial factors. These profiles can potentially untangle the interplay between individuals' variability in sleep, health, cognition and lifestyle - equipping research and clinical settings to better support individual's well-being.

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 non-invasively 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 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 association cortex) is linked to better cognitive performance. This result is replicated across North American (8.2 to 23.0 years old) and Asian (7.2 to 7.9 years 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.

Gradients of Brain Organization: Smooth Sailing from Methods Development to User Community.

Multimodal neuroimaging grants a powerful in vivo window into the structure and function of the human brain. Recent methodological and conceptual advances have enabled investigations of the interplay between large-scale spatial trends - or gradients - in brain structure and function, offering a framework to unify principles of brain organization across multiple scales. Strong community enthusiasm for these techniques has been instrumental in their widespread adoption and implementation to answer key questions in neuroscience. Following a brief review of current literature on this framework, this perspective paper will highlight how pragmatic steps aiming to make gradient methods more accessible to the community propelled these techniques to the forefront of neuroscientific inquiry. More specifically, we will emphasize how interest for gradient methods was catalyzed by data sharing, open-source software development, as well as the organization of dedicated workshops led by a diverse team of early career researchers. To this end, we argue that the growing excitement for brain gradients is the result of coordinated and consistent efforts to build an inclusive community and can serve as a case in point for future innovations and conceptual advances in neuroinformatics. We close this perspective paper by discussing challenges for the continuous refinement of neuroscientific theory, methodological innovation, and real-world translation to maintain our collective progress towards integrated models of brain organization.

A multidimensional investigation of sleep and biopsychosocial profiles with associated neural signatures.

Sleep is essential for optimal functioning and health. Interconnected to multiple biological, psychological and socio-environmental factors (i.e., biopsychosocial factors), the multidimensional nature of sleep is rarely capitalized on in research. Here, we deployed a data-driven approach to identify sleep-biopsychosocial profiles that linked self-reported sleep patterns to inter-individual variability in health, cognition, and lifestyle factors in 770 healthy young adults. We uncovered five profiles, including two profiles reflecting general psychopathology associated with either reports of general poor sleep or an absence of sleep complaints (i.e., sleep resilience) respectively. The three other profiles were driven by sedative-hypnotics-use and social satisfaction, sleep duration and cognitive performance, and sleep disturbance linked to cognition and mental health. Furthermore, identified sleep-biopsychosocial profiles displayed unique patterns of brain network organization. In particular, somatomotor network connectivity alterations were involved in the relationships between sleep and biopsychosocial factors. These profiles can potentially untangle the interplay between individuals' variability in sleep, health, cognition and lifestyle - equipping research and clinical settings to better support individual's well-being.

Organization of the human cerebral cortex estimated within individuals: networks, global topography, and function.

The cerebral cortex is populated by specialized regions that are organized into networks. Here we estimated networks from functional MRI (fMRI) data in intensively sampled participants. The procedure was developed in two participants (scanned 31 times) and then prospectively applied to 15 participants (scanned 8-11 times). Analysis of the networks revealed a global organization. Locally organized first-order sensory and motor networks were surrounded by spatially adjacent second-order networks that linked to distant regions. Third-order networks possessed regions distributed widely throughout association cortex. Regions of distinct third-order networks displayed side-by-side juxtapositions with a pattern that repeated across multiple cortical zones. We refer to these as supra-areal association megaclusters (SAAMs). Within each SAAM, two candidate control regions were adjacent to three separate domain-specialized regions. Response properties were explored with task data. The somatomotor and visual networks responded to body movements and visual stimulation, respectively. Second-order networks responded to transients in an oddball detection task, consistent with a role in orienting to salient events. The third-order networks, including distinct regions within each SAAM, showed two levels of functional specialization. Regions linked to candidate control networks responded to working memory load across multiple stimulus domains. The remaining regions dissociated across language, social, and spatial/episodic processing domains. These results suggest that progressively higher-order networks nest outward from primary sensory and motor cortices. Within the apex zones of association cortex, there is specialization that repeatedly divides domain-flexible from domain-specialized regions. We discuss implications of these findings, including how repeating organizational motifs may emerge during development.NEW & NOTEWORTHY The organization of cerebral networks was estimated within individuals with intensive, repeat sampling of fMRI data. A hierarchical organization emerged in each individual that delineated first-, second-, and third-order cortical networks. Regions of distinct third-order association networks consistently exhibited side-by-side juxtapositions that repeated across multiple cortical zones, with clear and robust functional specialization among the embedded regions.

Macroscale intrinsic dynamics are associated with microcircuit function in focal and generalized epilepsies.

Epilepsies are a group of neurological disorders characterized by abnormal spontaneous brain activity, involving multiscale changes in brain functional organizations. However, it is not clear to what extent the epilepsy-related perturbations of spontaneous brain activity affect macroscale intrinsic dynamics and microcircuit organizations, that supports their pathological relevance. We collect a sample of patients with temporal lobe epilepsy (TLE) and genetic generalized epilepsy with tonic-clonic seizure (GTCS), as well as healthy controls. We extract massive temporal features of fMRI BOLD time-series to characterize macroscale intrinsic dynamics, and simulate microcircuit neuronal dynamics used a large-scale biological model. Here we show whether macroscale intrinsic dynamics and microcircuit dysfunction are differed in epilepsies, and how these changes are linked. Differences in macroscale gradient of time-series features are prominent in the primary network and default mode network in TLE and GTCS. Biophysical simulations indicate reduced recurrent connection within somatomotor microcircuits in both subtypes, and even more reduced in GTCS. We further demonstrate strong spatial correlations between differences in the gradient of macroscale intrinsic dynamics and microcircuit dysfunction in epilepsies. These results emphasize the impact of abnormal neuronal activity on primary network and high-order networks, suggesting a systematic abnormality of brain hierarchical organization.

MRI economics: Balancing sample size and scan duration in brain wide association studies.

A pervasive dilemma in neuroimaging is whether to prioritize sample size or scan duration given fixed resources. Here, we systematically investigate this trade-off in the context of brain-wide association studies (BWAS) using resting-state functional magnetic resonance imaging (fMRI). We find that total scan duration (sample size × scan duration per participant) robustly explains individual-level phenotypic prediction accuracy via a logarithmic model, suggesting that sample size and scan duration are broadly interchangeable. The returns of scan duration eventually diminish relative to sample size, which we explain with principled theoretical derivations. When accounting for fixed costs associated with each participant (e.g., recruitment, non-imaging measures), we find that prediction accuracy in small-scale BWAS might benefit from much longer scan durations (>50 min) than typically assumed. Most existing large-scale studies might also have benefited from smaller sample sizes with longer scan durations. Both logarithmic and theoretical models of the relationships among sample size, scan duration and prediction accuracy explain well-predicted phenotypes better than poorly-predicted phenotypes. The logarithmic and theoretical models are also undermined by individual differences in brain states. These results replicate across phenotypic domains (e.g., cognition and mental health) from two large-scale datasets with different algorithms and metrics. Overall, our study emphasizes the importance of scan time, which is ignored in standard power calculations. Standard power calculations inevitably maximize sample size at the expense of scan duration. The resulting prediction accuracies are likely lower than would be produced with alternate designs, thus impeding scientific discovery. Our empirically informed reference is available for future study design: WEB_APPLICATION_LINK.

DeepResBat: deep residual batch harmonization accounting for covariate distribution differences.

Pooling MRI data from multiple datasets requires harmonization to reduce undesired inter-site variabilities, while preserving effects of biological variables (or covariates). The popular harmonization approach ComBat uses a mixed effect regression framework that explicitly accounts for covariate distribution differences across datasets. There is also significant interest in developing harmonization approaches based on deep neural networks (DNNs), such as conditional variational autoencoder (cVAE). However, current DNN approaches do not explicitly account for covariate distribution differences across datasets. Here, we provide mathematical results, suggesting that not accounting for covariates can lead to suboptimal harmonization. We propose two DNN-based covariate-aware harmonization approaches: covariate VAE (coVAE) and DeepResBat. The coVAE approach is a natural extension of cVAE by concatenating covariates and site information with site- and covariate-invariant latent representations. DeepResBat adopts a residual framework inspired by ComBat. DeepResBat first removes the effects of covariates with nonlinear regression trees, followed by eliminating site differences with cVAE. Finally, covariate effects are added back to the harmonized residuals. Using three datasets from three continents with a total of 2787 participants and 10085 anatomical T1 scans, we find that DeepResBat and coVAE outperformed ComBat, CovBat and cVAE in terms of removing dataset differences, while enhancing biological effects of interest. However, coVAE hallucinates spurious associations between anatomical MRI and covariates even when no association exists. Future studies proposing DNN-based harmonization approaches should be aware of this false positive pitfall. Overall, our results suggest that DeepResBat is an effective deep learning alternative to ComBat. Code for DeepResBat can be found here: https://github.com/ThomasYeoLab/CBIG/tree/master/stable_projects/harmonization/An2024_DeepResBat.