Adolescent Tuning of Association Cortex in Human Structural Brain Networks

The heritability and structural correlates of resting-state fMRI complexity

The complexity of fMRI signals quantifies temporal dynamics of spontaneous neural activity, which has been increasingly recognized as providing important insights into cognitive functions and psychiatric disorders. However, its heritability and structural underpinnings are not well understood. Here, we utilize multi-scale sample entropy to extract resting-state fMRI complexity in a large healthy adult sample from the Human Connectome Project. We show that fMRI complexity at multiple time scales is heritable in broad brain regions. Heritability estimates are modest and regionally variable. We relate fMRI complexity to brain structure including surface area, cortical myelination, cortical thickness, subcortical volumes, and total brain volume. We find that surface area is negatively correlated with fine-scale complexity and positively correlated with coarse-scale complexity in most cortical regions, especially the association cortex. Most of these correlations are related to common genetic and environmental effects. We also find positive correlations between cortical myelination and fMRI complexity at fine scales and negative correlations at coarse scales in the prefrontal cortex, lateral temporal lobe, precuneus, lateral parietal cortex, and cingulate cortex, with these correlations mainly attributed to common environmental effects. We detect few significant associations between fMRI complexity and cortical thickness. Despite the non-significant association with total brain volume, fMRI complexity exhibits significant correlations with subcortical volumes in the hippocampus, cerebellum, putamen, and pallidum at certain scales. Collectively, our work establishes the genetic basis and structural correlates of resting-state fMRI complexity across multiple scales, supporting its potential application as an endophenotype for psychiatric disorders.

Neural mechanism of non-adaptive cognitive emotion regulation in patients with non-suicidal self-injury

BACKGROUND

The incidence of non-suicidal self-injury (NSSI) has been on the rise in recent years. Studies have shown that people with NSSI have difficulties in emotion regulation and cognitive control. In addition, some studies have investigated the cognitive emotion regulation of people with NSSI which found that they have difficulties in cognitive emotion regulation, but there was a lack of research on cognitive emotion regulation strategies and related neural mechanisms.

METHODS

This study included 117 people with NSSI (age = 19.47 ± 5.13, male = 17) and 84 non-NSSI participants (age = 19.86 ± 4.14, male = 16). People with NSSI met the DSM-5 diagnostic criteria, and non-NSSI participants had no mental or physical disorders. The study collected all participants' data of Cognitive Emotion Regulation Questionnaire (CERQ) and functional magnetic resonance imaging (fMRI) to explore the differences in psychological performance and brain between two groups. Afterwards, Machine learning was used to select the found differential brain regions to obtain the highest correlation regions with NSSI. Then, Allen's Human Brain Atlas database was used to compare with the information on the abnormal brain regions of people with NSSI to find the genetic information related to NSSI. In addition, gene enrichment analysis was carried out to find the related pathways and specific cells that may have differences.

RESULTS

The differences between NSSI participants and non-NSSI participants were as follows: positive refocusing (t = -4.74, p < 0.01); refocusing on plans (t = -4.11, p < 0.01); positive reappraisal (t = -9.22, p < 0.01); self-blame (t = 6.30, p < 0.01); rumination (t = 3.64, p < 0.01); catastrophizing (t = 9.10, p < 0.01), and blaming others (t = 2.52, p < 0.01), the precentral gyrus (t = 6.04, pFDR < 0.05) and the rolandic operculum (t = -4.57, pFDR < 0.05). Rolandic operculum activity was negatively correlated with blaming others (r = -0.20, p < 0.05). Epigenetic results showed that excitatory neurons (p < 0.01) and inhibitory neurons (p < 0.01) were significant differences in two pathways, "trans-synaptic signaling" (p < -log108) and "modulation of chemical synaptic transmission" (p < -log108) in both cells.

CONCLUSIONS

People with NSSI are more inclined to adopt non-adaptive cognitive emotion regulation strategies. Rolandic operculum is also abnormally active. Abnormal changes in the rolandic operculum of them are associated with non-adaptive cognitive emotion regulation strategies. Changes in the excitatory and inhibitory neurons provide hints to explore the abnormalities of the neurological mechanisms at the cellular level of them. Trial registration number NCT04094623.

Quantitative susceptibility mapping of multiple system atrophy and Parkinson's disease correlates with neurotransmitter reference maps

BACKGROUND

Multiple system atrophy (MSA) and Parkinson's disease (PD) are neurodegenerative disorders characterized by α-synuclein pathology, disrupted iron homeostasis and impaired neurochemical transmission. Considering the critical role of iron in neurotransmitter synthesis and transport, our study aims to identify distinct patterns of whole-brain iron accumulation in MSA and PD, and to elucidate the corresponding neurochemical substrates.

METHODS

A total of 122 PD patients, 58 MSA patients and 78 age-, sex-matched health controls underwent multi-echo gradient echo sequences and neurological evaluations. We conducted voxel-wise and regional analyses using quantitative susceptibility mapping to explore MSA or PD-specific alterations in cortical and subcortical iron concentrations. Spatial correlation approaches were employed to examine the topographical alignment of cortical iron accumulation patterns with normative atlases of neurotransmitter receptor and transporter densities. Furthermore, we assessed the associations between the colocalization strength of neurochemical systems and disease severity.

RESULTS

MSA patients exhibited increased susceptibility in the striatal, midbrain, cerebellar nuclei, as well as the frontal, temporal, occipital lobes, and anterior cingulate gyrus. In contrast, PD patients displayed elevated iron levels in the left inferior occipital gyrus, precentral gyrus, and substantia nigra. The excessive iron accumulation in MSA or PD correlated with the spatial distribution of cholinergic, noradrenaline, glutamate, serotonin, cannabinoids, and opioid neurotransmitters, and the degree of this alignment was related to motor deficits.

CONCLUSIONS

Our findings provide evidence of the interaction between iron accumulation and non-dopamine neurotransmitters in the pathogenesis of MSA and PD, which inspires research on potential targets for pharmacotherapy.

Electroencephalography connectome changes in chronic insomnia disorder are correlated with neurochemical signatures

STUDY OBJECTIVES

This study aimed to investigate the alterations in resting-state electroencephalography (EEG) global brain connectivity (GBC) in patients with chronic insomnia disorder (CID) and to explore the correlation between macroscale connectomic variances and microscale neurotransmitter distributions.

METHODS

We acquired 64-channel EEG from 35 female CID patients and 34 healthy females. EEG signals were source-localized using individual brain anatomy and orthogonalized to mitigate volume conduction. Correlation coefficients between band-limited source-space power envelopes of the DK 68 atlas were computed and averaged across regions to determine specific GBC values. A support vector machine (SVM) classifier utilizing GBC features was employed to differentiate CID patients from controls. We further used Neurosynth and a 3D atlas of neurotransmitter receptors/transporters to assess the cognitive functions and neurotransmitter landscape associated with CID cortical abnormality maps, respectively.

RESULTS

CID patients exhibited elevated GBC within the medial prefrontal cortex and limbic cortex, particularly at the gamma carrier frequency, compared to controls (pFDR < .05). GBC patterns were found to effectively distinguish CID patients from controls with a precision of 90.8% in the SVM model. The cortical abnormality maps were significantly correlated with meta-analytic terms like "cognitive control" and "emotion regulation." Notably, GBC patterns were associated with neurotransmitter profiles (pspin < .05), with neurotransmitter systems such as norepinephrine, dopamine, and serotonin making significant contributions.

CONCLUSIONS

This work characterizes the EEG connectomic profile of CID, facilitating the cost-effective clinical translation of EEG-derived markers. Additionally, the linkage between GBC patterns and neurotransmitter distribution offers promising avenues for developing targeted treatment strategies for CID.

Obesity and the cerebral cortex: Underlying neurobiology in mice and humans

Obesity is a major modifiable risk factor for Alzheimer's disease (AD), characterized by progressive atrophy of the cerebral cortex. The neurobiology of obesity contributions to AD is poorly understood. Here we show with in vivo MRI that diet-induced obesity decreases cortical volume in mice, and that higher body adiposity associates with lower cortical volume in humans. Single-nuclei transcriptomics of the mouse cortex reveals that dietary obesity promotes an array of neuron-adverse transcriptional dysregulations, which are mediated by an interplay of excitatory neurons and glial cells, and which involve microglial activation and lowered neuronal capacity for neuritogenesis and maintenance of membrane potential. The transcriptional dysregulations of microglia, more than of other cell types, are like those in AD, as assessed with single-nuclei cortical transcriptomics in a mouse model of AD and two sets of human donors with the disease. Serial two-photon tomography of microglia demonstrates microgliosis throughout the mouse cortex. The spatial pattern of adiposity-cortical volume associations in human cohorts interrogated together with in silico bulk and single-nucleus transcriptomic data from the human cortex implicated microglia (along with other glial cells and subtypes of excitatory neurons), and it correlated positively with the spatial profile of cortical atrophy in patients with mild cognitive impairment and AD. Thus, multi-cell neuron-adverse dysregulations likely contribute to the loss of cortical tissue in obesity. The dysregulations of microglia may be pivotal to the obesity-related risk of AD.

Reconfigurations of cortical manifold structure during reward-based motor learning

Adaptive motor behavior depends on the coordinated activity of multiple neural systems distributed across the brain. While the role of sensorimotor cortex in motor learning has been well established, how higher-order brain systems interact with sensorimotor cortex to guide learning is less well understood. Using functional MRI, we examined human brain activity during a reward-based motor task where subjects learned to shape their hand trajectories through reinforcement feedback. We projected patterns of cortical and striatal functional connectivity onto a low-dimensional manifold space and examined how regions expanded and contracted along the manifold during learning. During early learning, we found that several sensorimotor areas in the dorsal attention network exhibited increased covariance with areas of the salience/ventral attention network and reduced covariance with areas of the default mode network (DMN). During late learning, these effects reversed, with sensorimotor areas now exhibiting increased covariance with DMN areas. However, areas in posteromedial cortex showed the opposite pattern across learning phases, with its connectivity suggesting a role in coordinating activity across different networks over time. Our results establish the neural changes that support reward-based motor learning and identify distinct transitions in the functional coupling of sensorimotor to transmodal cortex when adapting behavior.

Heterogeneity of morphometric similarity networks in health and schizophrenia

Introduction

Morphometric similarity is a recently developed neuroimaging phenotype of inter-regional connectivity by quantifying the similarity of a region to other regions based on multiple MRI parameters. Altered average morphometric similarity has been reported in psychotic disorders at the group level, with considerable heterogeneity across individuals. We used normative modeling to address cross-sectional and longitudinal inter-individual heterogeneity of morphometric similarity in health and schizophrenia.

Methods

Morphometric similarity for 62 cortical regions was obtained from baseline and follow-up T1-weighted scans of healthy individuals and patients with chronic schizophrenia. Cortical regions were classified into seven predefined brain functional networks. Using Bayesian Linear Regression and taking into account age, sex, image quality and scanner, we trained and validated normative models in healthy controls from eleven datasets (n = 4310). Individual deviations from the norm (z-scores) in morphometric similarity were computed for each participant for each network and region at both timepoints. A z-score ≧ than 1.96 was considered supra-normal and a z-score ≦ -1.96 infra-normal. As a longitudinal metric, we calculated the change over time of the total number of infra- or supra-normal regions per participant.

Results

At baseline, patients with schizophrenia had decreased morphometric similarity of the default mode network and increased morphometric similarity of the somatomotor network when compared with healthy controls. The percentage of patients with infra- or supra-normal values for any region at baseline and follow-up was low (<6%) and did not differ from healthy controls. Mean intra-group changes over time in the total number of infra- or supra-normal regions were small in schizophrenia and healthy control groups (<1) and there were no significant between-group differences.

Conclusions

In a case-control setting, a decrease of morphometric similarity within the default mode network may be a robust finding implicated in schizophrenia. However, normative modeling suggests that significant reductions and changes over time of regional morphometric similarity are evident only in a minority of patients.

Topological abnormalities of the morphometric similarity network of the cerebral cortex in schizophrenia

A morphometric similarity (MS) network can be constructed using multiple magnetic resonance imaging parameters of each cortical region. An MS network can be used to assess the similarity between cortical regions. Although MS networks can detect microstructural alterations and capture connections between histologically similar cortical areas, the influence of schizophrenia on the topological characteristics of MS networks remains unclear. We obtained T1- and diffusion-weighted images of 239 healthy controls and 190 individuals with schizophrenia to construct the MS network. Group comparisons of the mean MS of the cortical regions and subnetworks were performed. The strengths of the connections between the cortical regions and the global and nodal network indices were compared between the groups. Clinical associations with the network indices were tested using Spearman's rho. Compared with healthy controls, individuals with schizophrenia had significant group differences in the mean MS of several cortical regions and subnetworks. Individuals with schizophrenia had both superior and inferior strengths of connections between cortical regions compared with those of healthy controls. We observed regional abnormalities of the MS network in individuals with schizophrenia regarding lower centrality values of the pars opercularis, superior frontal, and superior temporal areas. Specific nodal network measures of the right pars opercularis and left superior temporal areas were associated with illness duration in individuals with schizophrenia. We identified regional abnormalities of the MS network in schizophrenia with the left superior temporal area possibly being a key region in topological organization and cortical connections.

Neuroimaging epicenters as potential sites of onset of the neuroanatomical pathology in schizophrenia

Schizophrenia lacks a clear definition at the neuroanatomical level, capturing the sites of origin and progress of this disorder. Using a network-theory approach called epicenter mapping on cross-sectional magnetic resonance imaging from 1124 individuals with schizophrenia, we identified the most likely "source of origin" of the structural pathology. Our results suggest that the Broca's area and adjacent frontoinsular cortex may be the epicenters of neuroanatomical pathophysiology in schizophrenia. These epicenters can predict an individual's response to treatment for psychosis. In addition, cross-diagnostic similarities based on epicenter mapping over of 4000 individuals diagnosed with neurological, neurodevelopmental, or psychiatric disorders appear to be limited. When present, these similarities are restricted to bipolar disorder, major depressive disorder, and obsessive-compulsive disorder. We provide a comprehensive framework linking schizophrenia-specific epicenters to multiple levels of neurobiology, including cognitive processes, neurotransmitter receptors and transporters, and human brain gene expression. Epicenter mapping may be a reliable tool for identifying the potential onset sites of neural pathophysiology in schizophrenia.

Association between hearing ability and cortical morphology in the elderly: multiparametric mapping, cognitive relevance, and neurobiological underpinnings

BACKGROUND

Hearing impairment is a common condition in the elderly. However, a comprehensive understanding of its neural correlates is still lacking.

METHODS

We recruited 284 elderly adults who underwent structural MRI, magnetic resonance spectroscopy, audiometry, and cognitive assessments. Individual hearing abilities indexed by pure tone average (PTA) were correlated with multiple structural MRI-derived cortical morphological indices. For regions showing significant correlations, mediation analyses were performed to examine their role in the relationship between hearing ability and cognitive function. Finally, the correlation maps between hearing ability and cortical morphology were linked with publicly available connectomic gradient, transcriptomic, and neurotransmitter maps.

FINDINGS

Poorer hearing was related to cortical thickness (CT) reductions in widespread regions and gyrification index (GI) reductions in the right Area 52 and Insular Granular Complex. The GI in the right Area 52 mediated the relationship between hearing ability and executive function. This mediating effect was further modulated by glutamate and N-acetylaspartate levels in the right auditory region. The PTA-CT correlation map followed microstructural connectomic hierarchy, were related to genes involved in certain biological processes (e.g., glutamate metabolic process), cell types (e.g., excitatory neurons and astrocytes), and developmental stages (i.e., childhood to young adulthood), and covaried with dopamine receptor 1, dopamine transporter, and fluorodopa. The PTA-GI correlation map was related to 5-hydroxytryptamine receptor 2a.

INTERPRETATION

Poorer hearing is associated with cortical thinning and folding reductions, which may be engaged in the relationship between hearing impairment and cognitive decline in the elderly and have different neurobiological substrates.

FUNDING

See the Acknowledgements section.

Unraveling the molecular relevance of brain phenotypes: A comparative analysis of null models and test statistics

Correlating transcriptional profiles with imaging-derived phenotypes has the potential to reveal possible molecular architectures associated with cognitive functions, brain development and disorders. Competitive null models built by resampling genes and self-contained null models built by spinning brain regions, along with varying test statistics, have been used to determine the significance of transcriptional associations. However, there has been no systematic evaluation of their performance in imaging transcriptomics analyses. Here, we evaluated the performance of eight different test statistics (mean, mean absolute value, mean squared value, max mean, median, Kolmogorov-Smirnov (KS), Weighted KS and the number of significant correlations) in both competitive null models and self-contained null models. Simulated brain maps (n = 1,000) and gene sets (n = 500) were used to calculate the probability of significance (Psig) for each statistical test. Our results suggested that competitive null models may result in false positive results driven by co-expression within gene sets. Furthermore, we demonstrated that the self-contained null models may fail to account for distribution characteristics (e.g., bimodality) of correlations between all available genes and brain phenotypes, leading to false positives. These two confounding factors interacted differently with test statistics, resulting in varying outcomes. Specifically, the sign-sensitive test statistics (i.e., mean, median, KS, Weighted KS) were influenced by co-expression bias in the competitive null models, while median and sign-insensitive test statistics were sensitive to the bimodality bias in the self-contained null models. Additionally, KS-based statistics produced conservative results in the self-contained null models, which increased the risk of false negatives. Comprehensive supplementary analyses with various configurations, including realistic scenarios, supported the results. These findings suggest utilizing sign-insensitive test statistics such as mean absolute value, max mean in the competitive null models and the mean as the test statistic for the self-contained null models. Additionally, adopting the confounder-matched (e.g., coexpression-matched) null models as an alternative to standard null models can be a viable strategy. Overall, the present study offers insights into the selection of statistical tests for imaging transcriptomics studies, highlighting areas for further investigation and refinement in the evaluation of novel and commonly used tests.

Brain structural covariances in the ageing brain in the UK Biobank

The morphologic properties of brain regions co-vary or correlate with each other. Here we investigated the structural covariances of cortical thickness and subcortical volumes in the ageing brain, along with their associations with age and cognition, using cross-sectional data from the UK Biobank (N = 42,075, aged 45-83 years, 53% female). As the structural covariance should be estimated in a group of participants, all participants were divided into 84 non-overlapping, equal-sized age groups ranging from the youngest to the oldest. We examined 84 cortical thickness covariances and subcortical covariances. Our findings include: (1) there were significant differences in the variability of structural covariance in the ageing process, including an increased variance, and a decreased entropy. (2) significant enrichment in pairwise correlations between brain regions within the occipital lobe was observed in all age groups; (3) structural covariance in older age, especially after the age of around 64, was significantly different from that in the youngest group (median age 48 years); (4) sixty-two of the total 528 pairs of cortical thickness correlations and 10 of the total 21 pairs of subcortical volume correlations showed significant associations with age. These trends varied, with some correlations strengthening, some weakening, and some reversing in direction with advancing age. Additionally, as ageing was associated with cognitive decline, most of the correlations with cognition displayed an opposite trend compared to age associated patterns of correlations.

Network enrichment significance testing in brain-phenotype association studies

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

Neural correlates of obesity across the lifespan

Associations between brain and obesity are bidirectional: changes in brain structure and function underpin over-eating, while chronic adiposity leads to brain atrophy. Investigating brain-obesity interactions across the lifespan can help better understand these relationships. This study explores the interaction between obesity and cortical morphometry in children, young adults, adults, and older adults. We also investigate the genetic, neurochemical, and cognitive correlates of the brain-obesity associations. Our findings reveal a pattern of lower cortical thickness in fronto-temporal brain regions associated with obesity across all age cohorts and varying age-dependent patterns in the remaining brain regions. In adults and older adults, obesity correlates with neurochemical changes and expression of inflammatory and mitochondrial genes. In children and older adults, adiposity is associated with modifications in brain regions involved in emotional and attentional processes. Thus, obesity might originate from cognitive changes during early adolescence, leading to neurodegeneration in later life through mitochondrial and inflammatory mechanisms.

Functional connectivity development along the sensorimotor-association axis enhances the cortical hierarchy

Human cortical maturation has been posited to be organized along the sensorimotor-association axis, a hierarchical axis of brain organization that spans from unimodal sensorimotor cortices to transmodal association cortices. Here, we investigate the hypothesis that the development of functional connectivity during childhood through adolescence conforms to the cortical hierarchy defined by the sensorimotor-association axis. We tested this pre-registered hypothesis in four large-scale, independent datasets (total n = 3355; ages 5-23 years): the Philadelphia Neurodevelopmental Cohort (n = 1207), Nathan Kline Institute-Rockland Sample (n = 397), Human Connectome Project: Development (n = 625), and Healthy Brain Network (n = 1126). Across datasets, the development of functional connectivity systematically varied along the sensorimotor-association axis. Connectivity in sensorimotor regions increased, whereas connectivity in association cortices declined, refining and reinforcing the cortical hierarchy. These consistent and generalizable results establish that the sensorimotor-association axis of cortical organization encodes the dominant pattern of functional connectivity development.

Opposite changes in morphometric similarity of medial reward and lateral non-reward orbitofrontal cortex circuits in obesity

Obesity has a profound impact on metabolic health thereby adversely affecting brain structure and function. However, the majority of previous studies used a single structural index to investigate the link between brain structure and body mass index (BMI), which hinders our understanding of structural covariance between regions in obesity. This study aimed to examine the relationship between macroscale cortical organization and BMI using novel morphometric similarity networks (MSNs). The individual MSNs were first constructed from individual eight multimodal cortical morphometric features between brain regions. Then the relationship between BMI and MSNs within the discovery sample of 434 participants was assessed. The key findings were further validated in an independent sample of 192 participants. We observed that the lateral non-reward orbitofrontal cortex (lOFC) exhibited decoupling (i.e., reduction in integration) in obesity, which was mainly manifested by its decoupling with the cognitive systems (i.e., DMN and FPN) while the medial reward orbitofrontal cortex (mOFC) showed de-differentiation (i.e., decrease in distinctiveness) in obesity, which was mainly represented by its de-differentiation with the cognitive and attention systems (i.e., DMN and VAN). Additionally, the lOFC showed de-differentiation with the visual system in obesity, while the mOFC showed decoupling with the visual system and hyper-coupling with the sensory-motor system in obesity. As an important first step in revealing the role of underlying structural covariance in body mass variability, the present study presents a novel mechanism that underlies the reward-control interaction imbalance in obesity, thus can inform future weight-management approaches.

Brain activation and connection across resting and motor-task states in patients with generalized tonic-clonic seizures

AIMS

Motor abnormalities have been identified as one common symptom in patients with generalized tonic-clonic seizures (GTCS) inspiring us to explore the disease in a motor execution condition, which might provide novel insight into the pathomechanism.

METHODS

Resting-state and motor-task fMRI data were collected from 50 patients with GTCS, including 18 patients newly diagnosed without antiepileptic drugs (ND_GTCS) and 32 patients receiving antiepileptic drugs (AEDs_GTCS). Motor activation and its association with head motion and cerebral gradients were assessed. Whole-brain network connectivity across resting and motor states was further calculated and compared between groups.

RESULTS

All patients showed over-activation in the postcentral gyrus and the ND_GTCS showed decreased activation in putamen. Specifically, activation maps of ND_GTCS showed an abnormal correlation with head motion and cerebral gradient. Moreover, we detected altered functional network connectivity in patients within states and across resting and motor states by using repeated-measures analysis of variance. Patients did not show abnormal connectivity in the resting state, while distributed abnormal connectivity in the motor-task state. Decreased across-state network connectivity was also found in all patients.

CONCLUSION

Convergent findings suggested the over-response of activation and connection of the brain to motor execution in GTCS, providing new clues to uncover motor susceptibility underlying the disease.

Early life stress and functional network topology in children

Brain networks are continuously modified throughout development, yet this plasticity can also make functional networks vulnerable to early life stress. Little is currently known about the effect of early life stress on the functional organization of the brain. The current study investigated the association between environmental stressors and network topology using data from the Adolescent Brain Cognitive DevelopmentSM (ABCD®) Study. Hierarchical modeling identified a general factor of environmental stress, representing the common variance across multiple stressors, as well as four subfactors including familial dynamics, interpersonal support, neighborhood SES deprivation, and urbanicity. Functional network topology metrics were obtained using graph theory at rest and during tasks of reward processing, inhibition, and affective working memory. The general factor of environmental stress was associated with less specialization of networks, represented by lower modularity at rest. Local metrics indicated that general environmental stress was also associated with less efficiency in the subcortical-cerebellar and visual networks while showing greater efficiency in the default mode network at rest. Subfactors of environmental stress were associated with differences in specialization and efficiency in select networks. The current study illustrates that a wide range of stressors in a child's environment are associated with differences in brain network topology.

Altered intersubject functional variability of brain white-matter in major depressive disorder and its association with gene expression profiles

Major depressive disorder (MDD) is a clinically heterogeneous disorder. Its mechanism is still unknown. Although the altered intersubject variability in functional connectivity (IVFC) within gray-matter has been reported in MDD, the alterations to IVFC within white-matter (WM-IVFC) remain unknown. Based on the resting-state functional MRI data of discovery (145 MDD patients and 119 healthy controls [HCs]) and validation cohorts (54 MDD patients, and 78 HCs), we compared the WM-IVFC between the two groups. We further assessed the meta-analytic cognitive functions related to the alterations. The discriminant WM-IVFC values were used to classify MDD patients and predict clinical symptoms in patients. In combination with the Allen Human Brain Atlas, transcriptome-neuroimaging association analyses were further conducted to investigate gene expression profiles associated with WM-IVFC alterations in MDD, followed by a set of gene functional characteristic analyses. We found extensive WM-IVFC alterations in MDD compared to HCs, which were associated with multiple behavioral domains, including sensorimotor processes and higher-order functions. The discriminant WM-IVFC could not only effectively distinguish MDD patients from HCs with an area under curve ranging from 0.889 to 0.901 across three classifiers, but significantly predict depression severity (r = 0.575, p = 0.002) and suicide risk (r = 0.384, p = 0.040) in patients. Furthermore, the variability-related genes were enriched for synapse, neuronal system, and ion channel, and predominantly expressed in excitatory and inhibitory neurons. Our results obtained good reproducibility in the validation cohort. These findings revealed intersubject functional variability changes of brain WM in MDD and its linkage with gene expression profiles, providing potential implications for understanding the high clinical heterogeneity of MDD.

Multimodal Associations of FKBP5 Methylation with Emotion-Regulatory Brain Circuits

BACKGROUND

Understanding the biological processes underlying individual differences in emotion regulation and stress responsivity is a key challenge for translational neuroscience. The gene FKBP5 is a core regulator in molecular stress signaling that is implicated in the development of psychiatric disorders. Yet it remains unclear how FKBP5 DNA methylation (DNAm) in peripheral blood relates to individual differences in measures of neural structure and function, and their relevance to daily-life stress responsivity.

METHODS

Here, we characterize multimodal correlates of FKBP5 DNAm by combining epigenetic data with neuroimaging and Ambulatory Assessment in a sample of 395 healthy individuals.

RESULTS

First, we show that FKBP5 demethylation as a psychiatric risk factor relates to an anxiety-associated reduction of gray matter volume in the ventromedial prefrontal cortex (vmPFC), a brain area that is involved in emotion regulation and mental health risk and resilience. This effect of epigenetic upregulation of FKBP5 on neuronal structure is more pronounced where FKBP5 is epigenetically downregulated at baseline. Leveraging 208 functional MRI scans during a well-established emotion processing task we find that FKBP5 DNAm in peripheral blood is associated with functional difference of prefrontal-limbic circuits modulating affective responsivity to daily stressors, which we measured using ecological momentary assessment in daily life.

CONCLUSIONS

Overall, we demonstrate how FKBP5 contributes to interindividual differences in neural and real-life affect regulation via structural and functional changes in prefrontal-limbic brain circuits.

A Method to Estimate Longitudinal Change Patterns in Functional Network Connectivity of the Developing Brain Relevant to Psychiatric Problems, Cognition, and Age

Aim: To develop an approach to evaluate multiple overlapping brain functional change patterns (FCPs) in functional network connectivity (FNC) and apply to study developmental changes in brain function. Introduction: FNC, the network analog of functional connectivity (FC), is commonly used to capture the intrinsic functional relationships among brain networks. Ongoing research on longitudinal changes of intrinsic FC across whole-brain functional networks has proven useful for characterizing age-related changes, but to date, there has been little focus on capturing multivariate patterns of FNC change with brain development. Methods: In this article, we introduce a novel approach to evaluate multiple overlapping FCPs by utilizing FNC matrices. We computed FNC matrices from the large-scale Adolescent Brain Cognitive Development data using fully automated spatially constrained independent component analysis (ICA). We next evaluated changes in these patterns for a 2-year period using a second-level ICA on the FNC change maps. Results: Our proposed approach reveals several highly structured (modular) FCPs and significant results including strong brain FC between visual and sensorimotor domains that increase with age. We also find several FCPs that are associated with longitudinal changes of psychiatric problems, cognition, and age in the developing brain. Interestingly, FCP cross-covariation, reflecting coupling between maximally independent FCPs, also shows significant differences between upper and lower quartile loadings for longitudinal changes in age, psychiatric problems, and cognition scores, as well as baseline age in the developing brain. FCP patterns and results were also found to be highly reliable based on analysis of data collected in a separate scan session. Conclusion: In sum, our results show evidence of consistent multivariate patterns of functional change in emerging adolescents and the proposed approach provides a useful and general tool to evaluate covarying patterns of whole-brain functional changes in longitudinal data. Impact statement In this article, we introduce a novel approach utilizing functional network connectivity (FNC) matrices to estimate multiple overlapping brain functional change patterns (FCPs). The findings demonstrate several well-structured FCPs that exhibit significant changes for a 2-year period, particularly in the functional connectivity between the visual and sensorimotor domains. In addition, we discover several FCPs that are associated with psychopathology, cognition, and age. Finally, our proposed approach for studying age-related FCPs represents a pioneering method that provides a valuable tool for assessing interconnected patterns of whole-brain functional changes in longitudinal data and may be useful to study change over time with applicability to many other areas, including the study of longitudinal changes within diagnostic groups, treatment effects, aging effects, and more.

General and specific patterns of cortical gene expression as spatial correlates of complex cognitive functioning

Gene expression varies across the brain. This spatial patterning denotes specialised support for particular brain functions. However, the way that a given gene's expression fluctuates across the brain may be governed by general rules. Quantifying patterns of spatial covariation across genes would offer insights into the molecular characteristics of brain areas supporting, for example, complex cognitive functions. Here, we use principal component analysis to separate general and unique gene regulatory associations with cortical substrates of cognition. We find that the region-to-region variation in cortical expression profiles of 8235 genes covaries across two major principal components: gene ontology analysis suggests these dimensions are characterised by downregulation and upregulation of cell-signalling/modification and transcription factors. We validate these patterns out-of-sample and across different data processing choices. Brain regions more strongly implicated in general cognitive functioning (g; 3 cohorts, total meta-analytic N = 39,519) tend to be more balanced between downregulation and upregulation of both major components (indicated by regional component scores). We then identify a further 29 genes as candidate cortical spatial correlates of g, beyond the patterning of the two major components (|β| range = 0.18 to 0.53). Many of these genes have been previously associated with clinical neurodegenerative and psychiatric disorders, or with other health-related phenotypes. The results provide insights into the cortical organisation of gene expression and its association with individual differences in cognitive functioning.

Spatial transcriptomic patterns underlying amyloid-β and tau pathology are associated with cognitive dysfunction in Alzheimer's disease

Amyloid-β (Aβ) and tau proteins accumulate within distinct neuronal systems in Alzheimer's disease (AD). Although it is not clear why certain brain regions are more vulnerable to Aβ and tau pathologies than others, gene expression may play a role. We study the association between brain-wide gene expression profiles and regional vulnerability to Aβ (gene-to-Aβ associations) and tau (gene-to-tau associations) pathologies by leveraging two large independent AD cohorts. We identify AD susceptibility genes and gene modules in a gene co-expression network with expression profiles specifically related to regional vulnerability to Aβ and tau pathologies in AD. In addition, we identify distinct biochemical pathways associated with the gene-to-Aβ and the gene-to-tau associations. These findings may explain the discordance between regional Aβ and tau pathologies. Finally, we propose an analytic framework, linking the identified gene-to-pathology associations to cognitive dysfunction in AD at the individual level, suggesting potential clinical implication of the gene-to-pathology associations.

Molecular signatures of attention networks

Attention network theory proposes three distinct types of attention-alerting, orienting, and control-that are supported by separate brain networks and modulated by different neurotransmitters, that is, norepinephrine, acetylcholine, and dopamine. Here, we explore the extent of cortical, genetic, and molecular dissociation of these three attention systems using multimodal neuroimaging. We evaluated the spatial overlap between fMRI activation maps from the attention network test (ANT) and cortex-wide gene expression data from the Allen Human Brain Atlas. The goal was to identify genes associated with each of the attention networks in order to determine whether specific groups of genes were co-expressed with the corresponding attention networks. Furthermore, we analyzed publicly available PET-maps of neurotransmitter receptors and transporters to investigate their spatial overlap with the attention networks. Our analyses revealed a substantial number of genes (3871 for alerting, 6905 for orienting, 2556 for control) whose cortex-wide expression co-varied with the activation maps, prioritizing several molecular functions such as the regulation of protein biosynthesis, phosphorylation, and receptor binding. Contrary to the hypothesized associations, the ANT activation maps neither aligned with the distribution of norepinephrine, acetylcholine, and dopamine receptor and transporter molecules, nor with transcriptomic profiles that would suggest clearly separable networks. Independence of the attention networks appeared additionally constrained by a high level of spatial dependency between the network maps. Future work may need to reconceptualize the attention networks in terms of their segregation and reevaluate the presumed independence at the neural and neurochemical level.

Combined Connectomics, MAPT Gene Expression, and Amyloid Deposition to Explain Regional Tau Deposition in Alzheimer Disease

OBJECTIVE

We aimed to test whether region-specific factors, including spatial expression patterns of the tau-encoding gene MAPT and regional levels of amyloid positron emission tomography (PET), enhance connectivity-based modeling of the spatial variability in tau-PET deposition in the Alzheimer disease (AD) spectrum.

METHODS

We included 685 participants (395 amyloid-positive participants within AD spectrum and 290 amyloid-negative controls) with tau-PET and amyloid-PET from 3 studies (Alzheimer's Disease Neuroimaging Initiative, 18 F-AV-1451-A05, and BioFINDER-1). Resting-state functional magnetic resonance imaging was obtained in healthy controls (n = 1,000) from the Human Connectome Project, and MAPT gene expression from the Allen Human Brain Atlas. Based on a brain-parcellation atlas superimposed onto all modalities, we obtained region of interest (ROI)-to-ROI functional connectivity, ROI-level PET values, and MAPT gene expression. In stepwise regression analyses, we tested connectivity, MAPT gene expression, and amyloid-PET as predictors of group-averaged and individual tau-PET ROI values in amyloid-positive participants.

RESULTS

Connectivity alone explained 21.8 to 39.2% (range across 3 studies) of the variance in tau-PET ROI values averaged across amyloid-positive participants. Stepwise addition of MAPT gene expression and amyloid-PET increased the proportion of explained variance to 30.2 to 46.0% and 45.0 to 49.9%, respectively. Similarly, for the prediction of patient-level tau-PET ROI values, combining all 3 predictors significantly improved the variability explained (mean adjusted R2 range across studies = 0.118-0.148, 0.156-0.196, and 0.251-0.333 for connectivity alone, connectivity plus MAPT expression, and all 3 modalities combined, respectively).

INTERPRETATION

Across 3 study samples, combining the functional connectome and molecular properties substantially enhanced the explanatory power compared to single modalities, providing a valuable tool to explain regional susceptibility to tau deposition in AD. ANN NEUROL 2024;95:274-287.

Structural connectome architecture shapes the maturation of cortical morphology from childhood to adolescence

Cortical thinning is an important hallmark of the maturation of brain morphology during childhood and adolescence. However, the connectome-based wiring mechanism that underlies cortical maturation remains unclear. Here, we show cortical thinning patterns primarily located in the lateral frontal and parietal heteromodal nodes during childhood and adolescence, which are structurally constrained by white matter network architecture and are particularly represented using a network-based diffusion model. Furthermore, connectome-based constraints are regionally heterogeneous, with the largest constraints residing in frontoparietal nodes, and are associated with gene expression signatures of microstructural neurodevelopmental events. These results are highly reproducible in another independent dataset. These findings advance our understanding of network-level mechanisms and the associated genetic basis that underlies the maturational process of cortical morphology during childhood and adolescence.

Modular subgraphs in large-scale connectomes underpin spontaneous co-fluctuation events in mouse and human brains

Previous studies have adopted an edge-centric framework to study fine-scale network dynamics in human fMRI. To date, however, no studies have applied this framework to data collected from model organisms. Here, we analyze structural and functional imaging data from lightly anesthetized mice through an edge-centric lens. We find evidence of "bursty" dynamics and events - brief periods of high-amplitude network connectivity. Further, we show that on a per-frame basis events best explain static FC and can be divided into a series of hierarchically-related clusters. The co-fluctuation patterns associated with each cluster centroid link distinct anatomical areas and largely adhere to the boundaries of algorithmically detected functional brain systems. We then investigate the anatomical connectivity undergirding high-amplitude co-fluctuation patterns. We find that events induce modular bipartitions of the anatomical network of inter-areal axonal projections. Finally, we replicate these same findings in a human imaging dataset. In summary, this report recapitulates in a model organism many of the same phenomena observed in previously edge-centric analyses of human imaging data. However, unlike human subjects, the murine nervous system is amenable to invasive experimental perturbations. Thus, this study sets the stage for future investigation into the causal origins of fine-scale brain dynamics and high-amplitude co-fluctuations. Moreover, the cross-species consistency of the reported findings enhances the likelihood of future translation.

The Variegation of Human Brain Vulnerability to Rare Genetic Disorders and Convergence With Behaviorally Defined Disorders

BACKGROUND

Diverse gene dosage disorders (GDDs) increase risk for psychiatric impairment, but characterization of GDD effects on the human brain has so far been piecemeal, with few simultaneous analyses of multiple brain features across different GDDs.

METHODS

Here, through multimodal neuroimaging of 3 aneuploidy syndromes (XXY [total n = 191, 92 control participants], XYY [total n = 81, 47 control participants], and trisomy 21 [total n = 69, 41 control participants]), we systematically mapped the effects of supernumerary X, Y, and chromosome 21 dosage across a breadth of 15 different macrostructural, microstructural, and functional imaging-derived phenotypes (IDPs).

RESULTS

The results revealed considerable diversity in cortical changes across GDDs and IDPs. This variegation of IDP change underlines the limitations of studying GDD effects unimodally. Integration across all IDP change maps revealed highly distinct architectures of cortical change in each GDD along with partial coalescence onto a common spatial axis of cortical vulnerability that is evident in all 3 GDDs. This common axis shows strong alignment with shared cortical changes in behaviorally defined psychiatric disorders and is enriched for specific molecular and cellular signatures.

CONCLUSIONS

Use of multimodal neuroimaging data in 3 aneuploidies indicates that different GDDs impose unique fingerprints of change in the human brain that differ widely depending on the imaging modality that is being considered. Embedded in this variegation is a spatial axis of shared multimodal change that aligns with shared brain changes across psychiatric disorders and therefore represents a major high-priority target for future translational research in neuroscience.

Age-dependent alterations in the coordinated development of subcortical regions in adolescents with social anxiety disorder

Subcortical brain regions play essential roles in the pathology of social anxiety disorder (SAD). While adolescence is the peak period of SAD, the relationships between altered development of the subcortical regions during this period and SAD are still unclear. This study investigated the age-dependent alterations in structural co-variance among subcortical regions and between subcortical and cortical regions, aiming to reflect aberrant coordination during development in the adolescent with SAD. High-resolution T1-weighted images were obtained from 76 adolescents with SAD and 67 healthy controls (HC), ranging from 11 to 17.9 years. Symptom severity was evaluated with the Social Anxiety Scale for Children (SASC) and the Depression Self Rating Scale for Children (DSRS-C). Structural co-variance and sliding age-window analyses were used to detect age-dependent group differences in inter-regional coordination patterns among subcortical regions and between subcortical and cortical regions. The volume of the striatum significantly correlated with SAD symptom severity. The SAD group exhibited significantly enhanced structural co-variance among key regions of the striatum (putamen and caudate). While the co-variance decreased with age in healthy adolescents, the co-variance in SAD adolescents stayed high, leading to more apparent group differences in middle adolescence. Moreover, the striatum's mean structural co-variance with cortical regions decreased with age in HC but increased with age in SAD. Adolescents with SAD suffer aberrant developmental coordination among the key regions of the striatum and between the striatum and cortical regions. The degree of incoordination is age-dependent, which may represent a neurodevelopmental trait of SAD.

Transcriptional patterns of the cortical Morphometric Inverse Divergence in first-episode, treatment-naïve early-onset schizophrenia

Early-onset Schizophrenia (EOS) is a profoundly progressive psychiatric disorder characterized by both positive and negative symptoms, whose pathogenesis is influenced by genes, environment and brain structure development. In this study, the MIND (Morphometric Inverse Divergence) network was employed to explore the relationship between morphological similarity and specific transcriptional expression patterns in EOS patients. This study involved a cohort of 187 participants aged between 7 and 17 years, consisting of 97 EOS patients and 90 healthy controls (HC). Multiple morphological features were used to construct the MIND network for all participants. Furthermore, we explored the associations between MIND network and brain-wide gene expression in EOS patients through partial least squares (PLS) regression, shared genetic predispositions with other psychiatric disorders, functional enrichment of PLS weighted genes, as well as transcriptional signature assessment of cell types, cortical layers, and developmental stages. The MIND showed similarity differences in the orbitofrontal cortex, pericalcarine cortex, lingual gyrus, and multiple networks in EOS patients compared to HC. Moreover, our exploration revealed a significant overlap of PLS2 weighted genes linking to EOS-related MIND differences and the dysregulated genes reported in other psychiatric diseases. Interestingly, genes correlated with MIND changes (PLS2-) exhibited a significant enrichment not only in metabolism-related pathways, but also in specific astrocytes, cortical layers (specifically layer I and III), and posterior developmental stages (late infancy to young adulthood stages). However, PLS2+ genes were primarily enriched in synapses signaling-related pathways and early developmental stages (from early-mid fetal to neonatal early infancy) but not in special cell types or layers. These findings provide a novel perspective on the intricate relationship between macroscopic morphometric structural abnormalities and microscopic transcriptional patterns during the onset and progression of EOS.

Single-subject cortical morphological brain networks: Phenotypic associations and neurobiological substrates

Although single-subject morphological brain networks provide an important way for human connectome studies, their roles and origins are poorly understood. Combining cross-sectional and repeated structural magnetic resonance imaging scans from adults, children and twins with behavioral and cognitive measures and brain-wide transcriptomic, cytoarchitectonic and chemoarchitectonic data, this study examined phenotypic associations and neurobiological substrates of single-subject morphological brain networks. We found that single-subject morphological brain networks explained inter-individual variance and predicted individual outcomes in Motor and Cognition domains, and distinguished individuals from each other. The performance can be further improved by integrating different morphological indices for network construction. Low-moderate heritability was observed for single-subject morphological brain networks with the highest heritability for sulcal depth-derived networks and higher heritability for inter-module connections. Furthermore, differential roles of genetic, cytoarchitectonic and chemoarchitectonic factors were observed for single-subject morphological brain networks. Cortical thickness-derived networks were related to the three factors with contributions from genes enriched in membrane and transport related functions, genes preferentially located in supragranular and granular layers, overall thickness in the molecular layer and thickness of wall in the infragranular layers, and metabotropic glutamate receptor 5 and dopamine transporter; fractal dimension-, gyrification index- and sulcal depth-derived networks were only associated with the chemoarchitectonic factor with contributions from different sets of neurotransmitter receptors. Most results were reproducible across different parcellation schemes and datasets. Altogether, this study demonstrates phenotypic associations and neurobiological substrates of single-subject morphological brain networks, which provide intermediate endophenotypes to link molecular and cellular architecture and behavior and cognition.

Subcortical influences on the topology of cortical networks align with functional processing hierarchies

fMRI of the human brain reveals spatiotemporal patterns of functional connectivity (FC), forming distinct cortical networks. Lately, subcortical contributions to these configurations are receiving renewed interest, but investigations rarely focus explicitly on their effects on cortico-cortical FC. Here, we employ a straightforward multivariable approach and graph-theoretic tools to assess subcortical impact on topological features of cortical networks. Given recent evidence showing that structures like the thalamus and basal ganglia integrate input from multiple networks, we expect increased segregation between cortical networks after removal of subcortical effects on their FC patterns. We analyze resting state data of young and healthy participants (male and female; N = 100) from the human connectome project. We find that overall, the cortical network architecture becomes less segregated, and more integrated, when subcortical influences are accounted for. Underlying these global effects are the following trends: 'Transmodal' systems become more integrated with the rest of the network, while 'unimodal' networks show the opposite effect. For single nodes this hierarchical organization is reflected by a close correspondence with the spatial layout of the principal gradient of FC (Margulies et al., 2016). Lastly, we show that the limbic system is significantly less coherent with subcortical influences removed. The findings are validated in a (split-sample) replication dataset. Our results provide new insight regarding the interplay between subcortex and cortical networks, by putting the integrative impact of subcortex in the context of macroscale patterns of cortical organization.

Resilient functioning is associated with altered structural brain network topology in adolescents exposed to childhood adversity

Childhood adversity is one of the strongest predictors of adolescent mental illness. Therefore, it is critical that the mechanisms that aid resilient functioning in individuals exposed to childhood adversity are better understood. Here, we examined whether resilient functioning was related to structural brain network topology. We quantified resilient functioning at the individual level as psychosocial functioning adjusted for the severity of childhood adversity in a large sample of adolescents (N = 2406, aged 14-24). Next, we examined nodal degree (the number of connections that brain regions have in a network) using brain-wide cortical thickness measures in a representative subset (N = 275) using a sliding window approach. We found that higher resilient functioning was associated with lower nodal degree of multiple regions including the dorsolateral prefrontal cortex, the medial prefrontal cortex, and the posterior superior temporal sulcus (z > 1.645). During adolescence, decreases in nodal degree are thought to reflect a normative developmental process that is part of the extensive remodeling of structural brain network topology. Prior findings in this sample showed that decreased nodal degree was associated with age, as such our findings of negative associations between nodal degree and resilient functioning may therefore potentially resemble a more mature structural network configuration in individuals with higher resilient functioning.

Differences in whole-brain metabolism are associated with the expression of genes related to neurovascular unit integrity and synaptic plasticity in temporal lobe epilepsy

PURPOSE

Temporal lobe epilepsy (TLE) is a common, polygenic epilepsy syndrome that involves glucose hypometabolism in the epileptogenic zone. However, the transcriptional and cellular signatures underlying the metabolism in TLE remain unclear.

METHODS

In this retrospective study, 2-[18F]-fluoro-2-deoxy-D-glucose ([18F]FDG) positron emission tomography (PET) scans of TLE patients (n = 104) who underwent anterior temporal lobectomy were consecutively collected between 2016 and 2021. The transcriptional profiles of TLE risk genes across the brain were identified by the gene expression analyses from six TLE patients and twelve postmortem donors (six from the Allen Human Brain Atlas). Integrating the neuroimaging and transcriptomic data, we examined the relationship between the expression of TLE-associated genes and metabolic alterations in TLE. Furthermore, we performed functional enrichment analyses of the genes with higher weight in partial least squares regression using Metascape.

RESULTS

A total of 104 patients with TLE (mean age 29 ± 9 years, 50% male) and 30 healthy controls (HCs) (mean age 31 ± 6 years, 53% male) were enrolled. Compared to that of HCs, patients with TLE showed hypometabolism in the temporal lobes and adjacent structures but hypermetabolism in the thalamus and basal ganglia. The cortical map of inter-group differences in cerebral metabolism was spatially correlated with the expression of a weighted combination of genes enriched in ontology terms and pathways related to neurovascular unit (NVU) integrity and synaptic plasticity.

DISCUSSION

Our findings, combined with the analysis of neuroimaging and transcriptional data, suggest that genes related to NVU integrity and synaptic plasticity may drive alterations to brain metabolism that mediate the genetic risk of TLE.

Neural variability in three major psychiatric disorders

Across the major psychiatric disorders (MPDs), a shared disruption in brain physiology is suspected. Here we investigate the neural variability at rest, a well-established behavior-relevant marker of brain function, and probe its basis in gene expression and neurotransmitter receptor profiles across the MPDs. We recruited 219 healthy controls and 279 patients with schizophrenia, major depressive disorder, or bipolar disorders (manic or depressive state). The standard deviation of blood oxygenation level-dependent signal (SDBOLD) obtained from resting-state fMRI was used to characterize neural variability. Transdiagnostic disruptions in SDBOLD patterns and their relationships with clinical symptoms and cognitive functions were tested by partial least-squares correlation. Moving beyond the clinical sample, spatial correlations between the observed patterns of SDBOLD disruption and postmortem gene expressions, Neurosynth meta-analytic cognitive functions, and neurotransmitter receptor profiles were estimated. Two transdiagnostic patterns of disrupted SDBOLD were discovered. Pattern 1 is exhibited in all diagnostic groups and is most pronounced in schizophrenia, characterized by higher SDBOLD in the language/auditory networks but lower SDBOLD in the default mode/sensorimotor networks. In comparison, pattern 2 is only exhibited in unipolar and bipolar depression, characterized by higher SDBOLD in the default mode/salience networks but lower SDBOLD in the sensorimotor network. The expression of pattern 1 related to the severity of clinical symptoms and cognitive deficits across MPDs. The two disrupted patterns had distinct spatial correlations with gene expressions (e.g., neuronal projections/cellular processes), meta-analytic cognitive functions (e.g., language/memory), and neurotransmitter receptor expression profiles (e.g., D2/serotonin/opioid receptors). In conclusion, neural variability is a potential transdiagnostic biomarker of MPDs with a substantial amount of its spatial distribution explained by gene expressions and neurotransmitter receptor profiles. The pathophysiology of MPDs can be traced through the measures of neural variability at rest, with varying clinical-cognitive profiles arising from differential spatial patterns of aberrant variability.

The genetic relationships between brain structure and schizophrenia

Genetic risks for schizophrenia are theoretically mediated by genetic effects on brain structure but it has been unclear which genes are associated with both schizophrenia and cortical phenotypes. We accessed genome-wide association studies (GWAS) of schizophrenia (N = 69,369 cases; 236,642 controls), and of three magnetic resonance imaging (MRI) metrics (surface area, cortical thickness, neurite density index) measured at 180 cortical areas (N = 36,843, UK Biobank). Using Hi-C-coupled MAGMA, 61 genes were significantly associated with both schizophrenia and one or more MRI metrics. Whole genome analysis with partial least squares demonstrated significant genetic covariation between schizophrenia and area or thickness of most cortical regions. Genetic similarity between cortical areas was strongly coupled to their phenotypic covariance, and genetic covariation between schizophrenia and brain phenotypes was strongest in the hubs of structural covariance networks. Pleiotropically associated genes were enriched for neurodevelopmental processes and positionally concentrated in chromosomes 3p21, 17q21 and 11p11. Mendelian randomization analysis indicated that genetically determined variation in a posterior cingulate cortical area could be causal for schizophrenia. Parallel analyses of GWAS on bipolar disorder, Alzheimer's disease and height showed that pleiotropic association with MRI metrics was stronger for schizophrenia compared to other disorders.

Diverse functional interaction driven by control-default network hubs supports creative thinking

Complex cognitive processes, like creative thinking, rely on interactions among multiple neurocognitive processes to generate effective and innovative behaviors on demand, for which the brain's connector hubs play a crucial role. However, the unique contribution of specific hub sets to creative thinking is unknown. Employing three functional magnetic resonance imaging datasets (total N = 1,911), we demonstrate that connector hub sets are organized in a hierarchical manner based on diversity, with "control-default hubs"-which combine regions from the frontoparietal control and default mode networks-positioned at the apex. Specifically, control-default hubs exhibit the most diverse resting-state connectivity profiles and play the most substantial role in facilitating interactions between regions with dissimilar neurocognitive functions, a phenomenon we refer to as "diverse functional interaction". Critically, we found that the involvement of control-default hubs in facilitating diverse functional interaction robustly relates to creativity, explaining both task-induced functional connectivity changes and individual creative performance. Our findings suggest that control-default hubs drive diverse functional interaction in the brain, enabling complex cognition, including creative thinking. We thus uncover a biologically plausible explanation that further elucidates the widely reported contributions of certain frontoparietal control and default mode network regions in creativity studies.

Network Enrichment Significance Testing in Brain-Phenotype Association Studies

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

Eigenmode-based approach reveals a decline in brain structure-function liberality across the human lifespan

While brain function is supported and constrained by the underlying structure, the connectome-based link estimated by current approaches is either relatively moderate or accompanied by high model complexity, with the essential principles underlying structure-function coupling remaining elusive. Here, by proposing a mapping method based on network eigendecomposition, we present a concise and strong correspondence between structure and function. We show that the explanation of functional connectivity can be significantly improved by incorporating interactions between different structural eigenmodes. We also demonstrate the pronounced advantage of the present mapping in capturing individual-specific information with simple implementation. Applying our methodology to the human lifespan, we find that functional diversity decreases with age, with functional interactions increasingly dominated by the leading functional mode. We also find that structure-function liberality weakens with age, which is driven by the decreases in functional components that are less constrained by anatomy, while the magnitude of structure-aligned components is preserved. Overall, our work enhances the understanding of structure-function coupling from a collective, connectome-oriented perspective and promotes a more refined identification of functional portions relevant to human aging, holding great potential for mechanistic insights into individual differences associated with cognition, development, and neurological disorders.

Association between gene expression and functional-metabolic architecture in Parkinson's disease

Gene expression plays a critical role in the pathogenesis of Parkinson's disease (PD). How gene expression profiles are correlated with functional-metabolic architecture remains obscure. We enrolled 34 PD patients and 25 age-and-sex-matched healthy controls for simultaneous 18 F-FDG-PET/functional MRI scanning during resting state. We investigated the functional gradients and the ratio of standard uptake value. Principal component analysis was used to further combine the functional gradients and glucose metabolism into functional-metabolic architecture. Using partial least squares (PLS) regression, we introduced the transcriptomic data from the Allen Institute of Brain Sciences to identify gene expression patterns underlying the affected functional-metabolic architecture in PD. Between-group comparisons revealed significantly higher gradient variation in the visual, somatomotor, dorsal attention, frontoparietal, default mode, and subcortical network (pFDR  < .048) in PD. Increased FDG-uptake was found in the somatomotor and ventral attention network while decreased FDG-uptake was found in the visual network (pFDR  < .008). Spatial correlation analysis showed consistently affected patterns of functional gradients and metabolism (p = 2.47 × 10-8 ). PLS analysis and gene ontological analyses further revealed that genes were mainly enriched for metabolic, catabolic, cellular response to ions, and regulation of DNA transcription and RNA biosynthesis. In conclusion, our study provided genetic pathological mechanism to explain imaging-defined brain functional-metabolic architecture of PD.

Single-subject cortical morphological brain networks across the adult lifespan

Age-related changes in focal cortical morphology have been well documented in previous literature; however, how interregional coordination patterns of the focal cortical morphology reorganize with advancing age is not well established. In this study, we performed a comprehensive analysis of the topological changes in single-subject morphological brain networks across the adult lifespan. Specifically, we constructed four types of single-subject morphological brain networks for 650 participants (aged from 18 to 88 years old), and characterized their topological organization using graph-based network measures. Age-related changes in the network measures were examined via linear, quadratic, and cubic models. We found profound age-related changes in global small-world attributes and efficiency, local nodal centralities, and interregional similarities of the single-subject morphological brain networks. The age-related changes were mainly embodied in cortical thickness networks, involved in frontal regions and highly connected hubs, concentrated on short-range connections, characterized by linear changes, and susceptible to connections between limbic, frontoparietal, and ventral attention networks. Intriguingly, nonlinear (i.e., quadratic or cubic) age-related changes were frequently found in the insula and limbic regions, and age-related cubic changes preferred long-range morphological connections. Finally, we demonstrated that the morphological similarity in cortical thickness between two frontal regions mediated the relationship between age and cognition measured by Cattell scores. Taken together, these findings deepen our understanding of adaptive changes of the human brain with advancing age, which may account for interindividual variations in behaviors and cognition.

Functional connectome fingerprinting across the lifespan

Systematic changes have been observed in the functional architecture of the human brain with advancing age. However, functional connectivity (FC) is also a powerful feature to detect unique "connectome fingerprints," allowing identification of individuals among their peers. Although fingerprinting has been robustly observed in samples of young adults, the reliability of this approach has not been demonstrated across the lifespan. We applied the fingerprinting framework to the Cambridge Centre for Ageing and Neuroscience cohort (n = 483 aged 18 to 89 years). We found that individuals are "fingerprintable" (i.e., identifiable) across independent functional MRI scans throughout the lifespan. We observed a U-shape distribution in the strength of "self-identifiability" (within-individual correlation across modalities), and "others-identifiability" (between-individual correlation across modalities), with a decrease from early adulthood into middle age, before improving in older age. FC edges contributing to self-identifiability were not restricted to specific brain networks and were different between individuals across the lifespan sample. Self-identifiability was additionally associated with regional brain volume. These findings indicate that individual participant-level identification is preserved across the lifespan despite the fact that its components are changing nonlinearly.

High-amplitude network co-fluctuations linked to variation in hormone concentrations over the menstrual cycle

Many studies have shown that the human endocrine system modulates brain function, reporting associations between fluctuations in hormone concentrations and brain connectivity. However, how hormonal fluctuations impact fast changes in brain network organization over short timescales remains unknown. Here, we leverage a recently proposed framework for modeling co-fluctuations between the activity of pairs of brain regions at a framewise timescale. In previous studies we showed that time points corresponding to high-amplitude co-fluctuations disproportionately contributed to the time-averaged functional connectivity pattern and that these co-fluctuation patterns could be clustered into a low-dimensional set of recurring "states." Here, we assessed the relationship between these network states and quotidian variation in hormone concentrations. Specifically, we were interested in whether the frequency with which network states occurred was related to hormone concentration. We addressed this question using a dense-sampling dataset (N = 1 brain). In this dataset, a single individual was sampled over the course of two endocrine states: a natural menstrual cycle and while the subject underwent selective progesterone suppression via oral hormonal contraceptives. During each cycle, the subject underwent 30 daily resting-state fMRI scans and blood draws. Our analysis of the imaging data revealed two repeating network states. We found that the frequency with which state 1 occurred in scan sessions was significantly correlated with follicle-stimulating and luteinizing hormone concentrations. We also constructed representative networks for each scan session using only "event frames"-those time points when an event was determined to have occurred. We found that the weights of specific subsets of functional connections were robustly correlated with fluctuations in the concentration of not only luteinizing and follicle-stimulating hormones, but also progesterone and estradiol.

Myelination and excitation-inhibition balance synergistically shape structure-function coupling across the human cortex

Recent work has demonstrated that the relationship between structural and functional connectivity varies regionally across the human brain, with reduced coupling emerging along the sensory-association cortical hierarchy. The biological underpinnings driving this expression, however, remain largely unknown. Here, we postulate that intracortical myelination and excitation-inhibition (EI) balance mediate the heterogeneous expression of structure-function coupling (SFC) and its temporal variance across the cortical hierarchy. We employ atlas- and voxel-based connectivity approaches to analyze neuroimaging data acquired from two groups of healthy participants. Our findings are consistent across six complementary processing pipelines: 1) SFC and its temporal variance respectively decrease and increase across the unimodal-transmodal and granular-agranular gradients; 2) increased myelination and lower EI-ratio are associated with more rigid SFC and restricted moment-to-moment SFC fluctuations; 3) a gradual shift from EI-ratio to myelination as the principal predictor of SFC occurs when traversing from granular to agranular cortical regions. Collectively, our work delivers a framework to conceptualize structure-function relationships in the human brain, paving the way for an improved understanding of how demyelination and/or EI-imbalances induce reorganization in brain disorders.

A phylogenetically-conserved axis of thalamocortical connectivity in the human brain

The thalamus enables key sensory, motor, emotive, and cognitive processes via connections to the cortex. These projection patterns are traditionally considered to originate from discrete thalamic nuclei, however recent work showing gradients of molecular and connectivity features in the thalamus suggests the organisation of thalamocortical connections occurs along a continuous dimension. By performing a joint decomposition of densely sampled gene expression and non-invasive diffusion tractography in the adult human thalamus, we define a principal axis of genetic and connectomic variation along a medial-lateral thalamic gradient. Projections along this axis correspond to an anterior-posterior cortical pattern and are aligned with electrophysiological properties of the cortex. The medial-lateral axis demonstrates phylogenetic conservation, reflects transitions in neuronal subtypes, and shows associations with neurodevelopment and common brain disorders. This study provides evidence for a supra-nuclear axis of thalamocortical organisation characterised by a graded transition in molecular properties and anatomical connectivity.

Robust estimation of cortical similarity networks from brain MRI

Structural similarity is a growing focus for magnetic resonance imaging (MRI) of connectomes. Here we propose Morphometric INverse Divergence (MIND), a new method to estimate within-subject similarity between cortical areas based on the divergence between their multivariate distributions of multiple MRI features. Compared to the prior approach of morphometric similarity networks (MSNs) on n > 11,000 scans spanning three human datasets and one macaque dataset, MIND networks were more reliable, more consistent with cortical cytoarchitectonics and symmetry and more correlated with tract-tracing measures of axonal connectivity. MIND networks derived from human T1-weighted MRI were more sensitive to age-related changes than MSNs or networks derived by tractography of diffusion-weighted MRI. Gene co-expression between cortical areas was more strongly coupled to MIND networks than to MSNs or tractography. MIND network phenotypes were also more heritable, especially edges between structurally differentiated areas. MIND network analysis provides a biologically validated lens for cortical connectomics using readily available MRI data.

Atypical functional connectivity hierarchy in Rolandic epilepsy

Functional connectivity hierarchy is an important principle in the process of brain functional organization and an important feature reflecting brain development. However, atypical brain network hierarchy organization in Rolandic epilepsy have not been systematically investigated. We examined connectivity alteration with age and its relation to epileptic incidence, cognition, or underlying genetic factors in 162 cases of Rolandic epilepsy and 117 typically developing children, by measuring fMRI multi-axis functional connectivity gradients. Rolandic epilepsy is characterized by contracting and slowing expansion of the functional connectivity gradients, highlighting the atypical age-related change of the connectivity hierarchy in segregation properties. The gradient alterations are relevant to seizure incidence, cognition, and connectivity deficit, and development-associated genetic basis. Collectively, our approach provides converging evidence for atypical connectivity hierarchy as a system-level substrate of Rolandic epilepsy, suggesting this is a disorder of information processing across multiple functional domains, and established a framework for large-scale brain hierarchical research.

Progressive trajectories of schizophrenia across symptoms, genes, and the brain

BACKGROUND

Schizophrenia is characterized by complex psychiatric symptoms and unclear pathological mechanisms. Most previous studies have focused on the morphological changes that occur over the development of the disease; however, the corresponding functional trajectories remain unclear. In the present study, we aimed to explore the progressive trajectories of patterns of dysfunction after diagnosis.

METHODS

Eighty-six patients with schizophrenia and 120 healthy controls were recruited as the discovery dataset. Based on multiple functional indicators of resting-state brain functional magnetic resonance imaging, we conducted a duration-sliding dynamic analysis framework to investigate trajectories in association with disease progression. Neuroimaging findings were associated with clinical symptoms and gene expression data from the Allen Human Brain Atlas database. A replication cohort of patients with schizophrenia from the University of California, Los Angeles, was used as the replication dataset for the validation analysis.

RESULTS

Five stage-specific phenotypes were identified. A symptom trajectory was characterized by positive-dominated, negative ascendant, negative-dominated, positive ascendant, and negative surpassed stages. Dysfunctional trajectories from primary and subcortical regions to higher-order cortices were recognized; these are associated with abnormal external sensory gating and a disrupted internal excitation-inhibition equilibrium. From stage 1 to stage 5, the importance of neuroimaging features associated with behaviors gradually shifted from primary to higher-order cortices and subcortical regions. Genetic enrichment analysis identified that neurodevelopmental and neurodegenerative factors may be relevant as schizophrenia progresses and highlighted multiple synaptic systems.

CONCLUSIONS

Our convergent results indicate that progressive symptoms and functional neuroimaging phenotypes are associated with genetic factors in schizophrenia. Furthermore, the identification of functional trajectories complements previous findings of structural abnormalities and provides potential targets for drug and non-drug interventions in different stages of schizophrenia.

Functional and Structural Longitudinal Change Patterns in Adolescent Brain

Functional magnetic resonance imaging (fMRI) and structural magnetic resonance imaging (sMRI) are two widely used techniques to analyze longitudinal brain functional and structural change in adolescents. Although longitudinal changes in intrinsic functional and structural changes have been studied separately, most studies focus on univariate change rather than estimating multivariate patterns of functional network connectivity (FNC) and gray matter (GM) changes with increased age. To analyze whole-brain structural and functional changes with increased age, we suggest two complementary techniques (1: linking of functional change pattern (FCP) to voxel-wise ∆GM and 2: the connection between FCP and structural change pattern (SCP)). In this study, we apply our approaches to the functional and GM data from the large-scale Adolescent Brain and Cognitive Development (ABCD) data. We find a significant correlation between FCP and voxel-wise ∆GM for two components. We also investigate the links between FCP and SCP and hypothesize that functional connectivity and GM continue to exhibit linked changes during adolescence.Clinical Relevance- This work captures the whole-brain functional and structural change patterns link by introducing two complementary techniques.

Spatial Correspondence of LSD-Induced Variations on Brain Functioning at Rest With Serotonin Receptor Expression

BACKGROUND

Lysergic acid diethylamide (LSD) is an atypical psychedelic compound that exerts its effects through pleiotropic actions, mainly involving 1A/2A serotoninergic (5-HT) receptor subtypes. However, the mechanisms by which LSD promotes a reorganization of the brain's functional activity and connectivity are still partially unknown.

METHODS

Our study analyzed resting-state functional magnetic resonance imaging data acquired from 15 healthy volunteers undergoing LSD single-dose intake. A voxelwise analysis investigated the alterations of the brain's intrinsic functional connectivity and local signal amplitude induced by LSD or by a placebo. Quantitative comparisons assessed the spatial overlap between these 2 indices of functional reorganization and the topography of receptor expression obtained from a publicly available collection of in vivo, whole-brain atlases. Finally, linear regression models explored the relationships between changes in resting-state functional magnetic resonance imaging and behavioral aspects of the psychedelic experience.

RESULTS

LSD elicited modifications of the cortical functional architecture that spatially overlapped with the distribution of serotoninergic receptors. Local signal amplitude and functional connectivity increased in regions belonging to the default mode and attention networks associated with high expression of 5-HT2A receptors. These functional changes correlate with the occurrence of simple and complex visual hallucinations. At the same time, a decrease in local signal amplitude and intrinsic connectivity was observed in limbic areas, which are dense with 5-HT1A receptors.

CONCLUSIONS

This study provides new insights into the neural processes underlying the brain network reconfiguration induced by LSD. It also identifies a topographical relationship between opposite effects on brain functioning and the spatial distribution of different 5-HT receptors.