An expanding manifold in transmodal regions characterizes adolescent reconfiguration of structural connectome organization

The relationship between demographic factors and brain hierarchical changes following antidepressant treatment in patients remitted from depression

To investigate the associations between demographic factors and brain hierarchical changes following successful selective serotonin reuptake inhibitor (SSRI) treatment, 57 major depressive disorder (MDD) patients who achieved remission after a 12-week SSRI treatment and 39 healthy controls (HCs) were recruited. MDD patients underwent diffusion tensor imaging (DTI) scans before treatment and after the 12-week SSRI treatment. Depression severity was evaluated with the Hamilton Rating Scale for Depression (HAMD) using the total score and the subscales: retardation, cognitive impairment, anxiety, and sleep disturbance. All HCs also underwent DTI scans after enrollment. Building on gradient mapping techniques, we developed a set of measures to quantify the dispersion within functional communities and also studied demographic-relevant differences in the three-dimensional gradient space of remitted MDD patients. We defined the Z-scores of the gradients in the pre-treatment group relative to the HC group as the disease pattern, post-treatment group relative to the HC group as the recovery pattern. The results showed that the disease pattern of depression is associated with age, as older age groups exhibit more severe impairments in depression. A significant difference was detected in the dispersion of the frontoparietal network (FPN) between pre-treatment and post-treatment patients. With the moderating effect of the age of onset, the dispersion of the FPN was related to the improvement in cognitive impairment, the dorsal attention network (DAN) was related to the improvement in retardation symptoms. Our findings help clinicians be alert to the role of demographic effects on clinical efficacy when treating depressed patients.

Changes in cortical manifold structure following stroke and its relation to behavioral recovery in the male macaque

Stroke, a major cause of disability, disrupts brain function and motor skills. Previous research has mainly focused on reorganization of the motor system post-stroke, but the effects on other brain areas and their influence on recovery is poorly understood. Here, we use functional neuroimaging in a nonhuman primate model (23 male Cynomolgus Macaques), we explore how ischemic stroke affects whole-brain cortical architecture and its relation to spontaneous behavioral recovery. By projecting patterns of cortical functional connectivity onto a low-dimensional manifold space, we find that several regions in both sensorimotor cortex and higher-order transmodal cortex exhibit significant shifts in their manifold embedding from pre- to post-stroke. Furthermore, we observe that changes in default mode and limbic network regions, and not preserved sensorimotor cortical regions, are associated with animal behavioral recovery post-stroke. These results establish the whole-brain functional changes associated with stroke, and suggest an important role for higher-order transmodal cortex in post-stroke outcomes.

Prognostic model for predicting Alzheimer's disease conversion using functional connectome manifolds

BACKGROUND

Early detection of Alzheimer's disease (AD) is essential for timely management and consideration of therapeutic options; therefore, detecting the risk of conversion from mild cognitive impairment (MCI) to AD is crucial during neurodegenerative progression. Existing neuroimaging studies have mostly focused on group differences between individuals with MCI (or AD) and cognitively normal (CN), discarding the temporal information of conversion time. Here, we aimed to develop a prognostic model for AD conversion using functional connectivity (FC) and Cox regression suitable for conversion event modeling.

METHODS

We developed a prognostic model using a large-scale Alzheimer's Disease Neuroimaging Initiative dataset, and it was validated using external data obtained from the Open Access Series of Imaging Studies. We considered individuals who were initially CN or had MCI but progressed to AD and those with MCI with no progression to AD during the five-year follow-up period. As the exact conversion time to AD is unknown, we inferred this information using imputation approaches. We generated cortex-wide principal FC gradients using manifold learning techniques and computed subcortical-weighted manifold degrees from baseline functional magnetic resonance imaging data. A penalized Cox regression model with an elastic net penalty was adopted to define a risk score predicting the risk of conversion to AD, using FC gradients and clinical factors as regressors.

RESULTS

Our prognostic model predicted the conversion risk and confirmed the role of imaging-derived manifolds in the conversion risk. The brain regions that largely contributed to predicting AD conversion were the heteromodal association and visual cortices, as well as the caudate and hippocampus. Our risk score based on Cox regression was consistent with the expected disease trajectories and correlated with positron emission tomography tracer uptake and symptom severity, reinforcing its clinical usefulness. Our findings were validated using an independent dataset. The cross-sectional application of our model showed a higher risk for AD than that for MCI, which correlated with symptom severity scores in the validation dataset.

CONCLUSION

We proposed a prognostic model predicting the risk of conversion to AD. The associated risk score may provide insights for early intervention in individuals at risk of AD conversion.

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.

Exploring nonlinear dynamics in brain functionality through phase portraits and fuzzy recurrence plots

Much of the complexity and diversity found in nature is driven by nonlinear phenomena, and this holds true for the brain. Nonlinear dynamics theory has been successfully utilized in explaining brain functions from a biophysics standpoint, and the field of statistical physics continues to make substantial progress in understanding brain connectivity and function. This study delves into complex brain functional connectivity using biophysical nonlinear dynamics approaches. We aim to uncover hidden information in high-dimensional and nonlinear neural signals, with the hope of providing a useful tool for analyzing information transitions in functionally complex networks. By utilizing phase portraits and fuzzy recurrence plots, we investigated the latent information in the functional connectivity of complex brain networks. Our numerical experiments, which include synthetic linear dynamics neural time series and a biophysically realistic neural mass model, showed that phase portraits and fuzzy recurrence plots are highly sensitive to changes in neural dynamics and can also be used to predict functional connectivity based on structural connectivity. Furthermore, the results showed that phase trajectories of neuronal activity encode low-dimensional dynamics, and the geometric properties of the limit-cycle attractor formed by the phase portraits can be used to explain the neurodynamics. Additionally, our results showed that the phase portrait and fuzzy recurrence plots can be used as functional connectivity descriptors, and both metrics were able to capture and explain nonlinear dynamics behavior during specific cognitive tasks. In conclusion, our findings suggest that phase portraits and fuzzy recurrence plots could be highly effective as functional connectivity descriptors, providing valuable insights into nonlinear dynamics in the brain.

Hippocampal connectivity patterns echo macroscale cortical evolution in the primate brain

While the hippocampus is key for human cognitive abilities, it is also a phylogenetically old cortex and paradoxically considered evolutionarily preserved. Here, we introduce a comparative framework to quantify preservation and reconfiguration of hippocampal organisation in primate evolution, by analysing the hippocampus as an unfolded cortical surface that is geometrically matched across species. Our findings revealed an overall conservation of hippocampal macro- and micro-structure, which shows anterior-posterior and, perpendicularly, subfield-related organisational axes in both humans and macaques. However, while functional organisation in both species followed an anterior-posterior axis, we observed a marked reconfiguration in the latter across species, which mirrors a rudimentary integration of the default-mode-network in non-human primates. Here we show that microstructurally preserved regions like the hippocampus may still undergo functional reconfiguration in primate evolution, due to their embedding within heteromodal association networks.

Connectome reorganization associated with temporal lobe pathology and its surgical resection

Network neuroscience offers a unique framework to understand the organizational principles of the human brain. Despite recent progress, our understanding of how the brain is modulated by focal lesions remains incomplete. Resection of the temporal lobe is the most effective treatment to control seizures in pharmaco-resistant temporal lobe epilepsy (TLE), making this syndrome a powerful model to study lesional effects on network organization in young and middle-aged adults. Here, we assessed the downstream consequences of a focal lesion and its surgical resection on the brain's structural connectome, and explored how this reorganization relates to clinical variables at the individual patient level. We included adults with pharmaco-resistant TLE (n = 37) who underwent anterior temporal lobectomy between two imaging time points, as well as age- and sex-matched healthy controls who underwent comparable imaging (n = 31). Core to our analysis was the projection of high-dimensional structural connectome data-derived from diffusion MRI tractography from each subject-into lower-dimensional gradients. We then compared connectome gradients in patients relative to controls before surgery, tracked surgically-induced connectome reconfiguration from pre- to postoperative time points, and examined associations to patient-specific clinical and imaging phenotypes. Before surgery, individuals with TLE presented with marked connectome changes in bilateral temporo-parietal regions, reflecting an increased segregation of the ipsilateral anterior temporal lobe from the rest of the brain. Surgery-induced connectome reorganization was localized to this temporo-parietal subnetwork, but primarily involved postoperative integration of contralateral regions with the rest of the brain. Using a partial least-squares analysis, we uncovered a latent clinical imaging signature underlying this pre- to postoperative connectome reorganization, showing that patients who displayed postoperative integration in bilateral fronto-occipital cortices also had greater preoperative ipsilateral hippocampal atrophy, lower seizure frequency and secondarily generalized seizures. Our results bridge the effects of focal brain lesions and their surgical resections with large-scale network reorganization and interindividual clinical variability, thus offering new avenues to examine the fundamental malleability of the human brain.

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.

Tau follows principal axes of functional and structural brain organization in Alzheimer's disease

Alzheimer's disease (AD) is a brain network disorder where pathological proteins accumulate through networks and drive cognitive decline. Yet, the role of network connectivity in facilitating this accumulation remains unclear. Using in-vivo multimodal imaging, we show that the distribution of tau and reactive microglia in humans follows spatial patterns of connectivity variation, the so-called gradients of brain organization. Notably, less distinct connectivity patterns ("gradient contraction") are associated with cognitive decline in regions with greater tau, suggesting an interaction between reduced network differentiation and tau on cognition. Furthermore, by modeling tau in subject-specific gradient space, we demonstrate that tau accumulation in the frontoparietal and temporo-occipital cortices is associated with greater baseline tau within their functionally and structurally connected hubs, respectively. Our work unveils a role for both functional and structural brain organization in pathology accumulation in AD, and supports subject-specific gradient space as a promising tool to map disease progression.

Shifts in structural connectome organization in the limbic and sensory systems of patients with episodic migraine

Migraine is a complex neurological condition characterized by recurrent headaches, which is often accompanied by various neurological symptoms. Magnetic resonance imaging (MRI) is a powerful tool for investigating whole-brain connectivity patterns; however, systematic assessment of structural connectome organization has rarely been performed. In the present study, we aimed to examine the changes in structural connectivity in patients with episodic migraines using diffusion MRI. First, we computed structural connectivity using diffusion MRI tractography, after which we applied dimensionality reduction techniques to the structural connectivity and generated three low-dimensional eigenvectors. We subsequently calculated the manifold eccentricity, defined as the Euclidean distance between each data point and the center of the data in the manifold space. We then compared the manifold eccentricity between patients with migraines and healthy controls, revealing significant between-group differences in the orbitofrontal cortex, temporal pole, and sensory/motor regions. Between-group differences in subcortico-cortical connectivity further revealed significant changes in the amygdala, accumbens, and caudate nuclei. Finally, supervised machine learning effectively classified patients with migraines and healthy controls using cortical and subcortical structural connectivity features, highlighting the importance of the orbitofrontal and sensory cortices, in addition to the caudate, in distinguishing between the groups. Our findings confirmed that episodic migraine is related to the structural connectome changes in the limbic and sensory systems, suggesting its potential utility as a diagnostic marker for migraine.

Age-related changes of dopamine D1 and D2 receptors expression in parvalbumin-positive cells of the orbitofrontal and prelimbic cortices of mice

Dopamine (DA) plays a pivotal role in reward processing, cognitive functions, and emotional regulation. The prefrontal cortex (PFC) is a critical brain region for these processes. Parvalbumin-positive (PV+) neurons are one of the major classes of inhibitory GABAergic neurons in the cortex, they modulate the activity of neighboring neurons, influencing various brain functions. While DA receptor expression exhibits age-related changes, the age-related changes of these receptors in PV+ neurons, especially in the PFC, remain unclear. To address this, we investigated the expression of DA D1 (D1R) and D2 (D2R) receptors in PV+ neurons within the orbitofrontal (OFC) and prelimbic (PrL) cortices at different postnatal ages (P28, P42, P56, and P365). We found that the expression of D1R and D2R in PV+ neurons showed both age- and region-related changes. PV+ neurons in the OFC expressed a higher abundance of D1 than those in the PrL, and those neurons in the OFC also showed higher co-expression of D1R and D2R than those in the PrL. In the OFC and PrL, D1R in PV+ neurons increased from P28 and reached a plateau at P42, then receded to express at P365. Meanwhile, D2R did not show significant age-related changes between the two regions except at P56. These results showed dopamine receptors in the prefrontal cortex exhibit age- and region-specific changes, which may contribute to the difference of these brain regions in reward-related brain functions.

Connectome architecture shapes large-scale cortical alterations in schizophrenia: a worldwide ENIGMA study

Schizophrenia is a prototypical network disorder with widespread brain-morphological alterations, yet it remains unclear whether these distributed alterations robustly reflect the underlying network layout. We tested whether large-scale structural alterations in schizophrenia relate to normative structural and functional connectome architecture, and systematically evaluated robustness and generalizability of these network-level alterations. Leveraging anatomical MRI scans from 2439 adults with schizophrenia and 2867 healthy controls from 26 ENIGMA sites and normative data from the Human Connectome Project (n = 207), we evaluated structural alterations of schizophrenia against two network susceptibility models: (i) hub vulnerability, which examines associations between regional network centrality and magnitude of disease-related alterations; (ii) epicenter mapping, which identifies regions whose typical connectivity profile most closely resembles the disease-related morphological alterations. To assess generalizability and specificity, we contextualized the influence of site, disease stages, and individual clinical factors and compared network associations of schizophrenia with that found in affective disorders. Our findings show schizophrenia-related cortical thinning is spatially associated with functional and structural hubs, suggesting that highly interconnected regions are more vulnerable to morphological alterations. Predominantly temporo-paralimbic and frontal regions emerged as epicenters with connectivity profiles linked to schizophrenia's alteration patterns. Findings were robust across sites, disease stages, and related to individual symptoms. Moreover, transdiagnostic comparisons revealed overlapping epicenters in schizophrenia and bipolar, but not major depressive disorder, suggestive of a pathophysiological continuity within the schizophrenia-bipolar-spectrum. In sum, cortical alterations over the course of schizophrenia robustly follow brain network architecture, emphasizing marked hub susceptibility and temporo-frontal epicenters at both the level of the group and the individual. Subtle variations of epicenters across disease stages suggest interacting pathological processes, while associations with patient-specific symptoms support additional inter-individual variability of hub vulnerability and epicenters in schizophrenia. Our work outlines potential pathways to better understand macroscale structural alterations, and inter- individual variability in schizophrenia.

Diffusion time-related structure-function coupling reveals differential association with inter-individual variations in body mass index

Body mass index (BMI) is an indicator of obesity, and recent neuroimaging studies have demonstrated that inter-individual variations in BMI are associated with altered brain structure and function. However, the mechanism underlying the alteration of structure-function correspondence according to BMI is under-investigated. In this study, we studied structural and functional connectivity derived from diffusion MRI tractography and inter-regional correlations of functional MRI time series, respectively. We combined the structural and functional connectivity information using the Riemannian optimization approach. First, the low-dimensional principal eigenvectors (i.e., gradients) of the structural connectivity were generated by applying diffusion map embedding with varying diffusion times. A transformation was identified so that the structural and functional embeddings share the same coordinate system, and subsequently, the functional connectivity matrix was simulated. Then, we generated gradients from the simulated functional connectivity matrix. We found the most apparent cortical hierarchical organization differentiating between low-level sensory and higher-order transmodal regions in the middle of the diffusion time, indicating that the hierarchical organization of the brain may reflect the intermediate mechanisms of mono- and polysynaptic communications. Associations between the functional gradients and BMI were strongest when the hierarchical structure was the most evident. Moreover, the gradient-BMI association map was related to the microstructural features, and the findings indicated that the BMI-related structure-function coupling was significantly associated with brain microstructure, particularly in higher-order transmodal areas. Finally, transcriptomic association analysis revealed the potential biological underpinnings specifying gene enrichment in the striatum, hypothalamus, and cortical cells. Our findings provide evidence that structure-function correspondence is strongly coupled with BMI when hierarchical organization is the most apparent and that the associations are related to the multiscale properties of the brain, leading to an advanced understanding of the neural mechanisms related to BMI.

Structural connectome alterations between individuals with autism and neurotypical controls using feature representation learning

Autism spectrum disorder is one of the most common neurodevelopmental conditions associated with sensory and social communication impairments. Previous neuroimaging studies reported that atypical nodal- or network-level functional brain organization in individuals with autism was associated with autistic behaviors. Although dimensionality reduction techniques have the potential to uncover new biomarkers, the analysis of whole-brain structural connectome abnormalities in a low-dimensional latent space is underinvestigated. In this study, we utilized autoencoder-based feature representation learning for diffusion magnetic resonance imaging-based structural connectivity in 80 individuals with autism and 61 neurotypical controls that passed strict quality controls. We generated low-dimensional latent features using the autoencoder model for each group and adopted an integrated gradient approach to assess the contribution of the input data for predicting latent features during the encoding process. Subsequently, we compared the integrated gradient values between individuals with autism and neurotypical controls and observed differences within the transmodal regions and between the sensory and limbic systems. Finally, we identified significant associations between integrated gradient values and communication abilities in individuals with autism. Our findings provide insights into the whole-brain structural connectome in autism and may help identify potential biomarkers for autistic connectopathy.

Connectome-wide structure-function coupling models implicate polysynaptic alterations in autism

Autism spectrum disorder (ASD) is one of the most common neurodevelopmental diagnoses. Although incompletely understood, structural and functional network alterations are increasingly recognized to be at the core of the condition. We utilized multimodal imaging and connectivity modeling to study structure-function coupling in ASD and probed mono- and polysynaptic mechanisms on structurally-governed network function. We examined multimodal magnetic resonance imaging data in 80 ASD and 61 neurotypical controls from the Autism Brain Imaging Data Exchange (ABIDE) II initiative. We predicted intrinsic functional connectivity from structural connectivity data in each participant using a Riemannian optimization procedure that varies the times that simulated signals can unfold along tractography-derived personalized connectomes. In both ASD and neurotypical controls, we observed improved structure-function prediction at longer diffusion time scales, indicating better modeling of brain function when polysynaptic mechanisms are accounted for. Prediction accuracy differences (∆prediction accuracy) were marked in transmodal association systems, such as the default mode network, in both neurotypical controls and ASD. Differences were, however, lower in ASD in a polysynaptic regime at higher simulated diffusion times. We compared regional differences in ∆prediction accuracy between both groups to assess the impact of polysynaptic communication on structure-function coupling. This analysis revealed that between-group differences in ∆prediction accuracy followed a sensory-to-transmodal cortical hierarchy, with an increased gap between controls and ASD in transmodal compared to sensory/motor systems. Multivariate associative techniques revealed that structure-function differences reflected inter-individual differences in autistic symptoms and verbal as well as non-verbal intelligence. Our network modeling approach sheds light on atypical structure-function coupling in autism, and suggests that polysynaptic network mechanisms are implicated in the condition and that these can help explain its wide range of associated symptoms.

Identifying subgroups of eating behavior traits unrelated to obesity using functional connectivity and feature representation learning

Eating behavior is highly heterogeneous across individuals and cannot be fully explained using only the degree of obesity. We utilized unsupervised machine learning and functional connectivity measures to explore the heterogeneity of eating behaviors measured by a self-assessment instrument using 424 healthy adults (mean ± standard deviation [SD] age = 47.07 ± 18.89 years; 67% female). We generated low-dimensional representations of functional connectivity using resting-state functional magnetic resonance imaging and estimated latent features using the feature representation capabilities of an autoencoder by nonlinearly compressing the functional connectivity information. The clustering approaches applied to latent features identified three distinct subgroups. The subgroups exhibited different levels of hunger traits, while their body mass indices were comparable. The results were replicated in an independent dataset consisting of 212 participants (mean ± SD age = 38.97 ± 19.80 years; 35% female). The model interpretation technique of integrated gradients revealed that the between-group differences in the integrated gradient maps were associated with functional reorganization in heteromodal association and limbic cortices and reward-related subcortical structures such as the accumbens, amygdala, and caudate. The cognitive decoding analysis revealed that these systems are associated with reward- and emotion-related systems. Our findings provide insights into the macroscopic brain organization of eating behavior-related subgroups independent of obesity.

Longitudinal development of the human white matter structural connectome and its association with brain transcriptomic and cellular architecture

From childhood to adolescence, the spatiotemporal development pattern of the human brain white matter connectome and its underlying transcriptomic and cellular mechanisms remain largely unknown. With a longitudinal diffusion MRI cohort of 604 participants, we map the developmental trajectory of the white matter connectome from global to regional levels and identify that most brain network properties followed a linear developmental trajectory. Importantly, connectome-transcriptomic analysis reveals that the spatial development pattern of white matter connectome is potentially regulated by the transcriptomic architecture, with positively correlated genes involve in ion transport- and development-related pathways expressed in excitatory and inhibitory neurons, and negatively correlated genes enriches in synapse- and development-related pathways expressed in astrocytes, inhibitory neurons and microglia. Additionally, the macroscale developmental pattern is also associated with myelin content and thicknesses of specific laminas. These findings offer insights into the underlying genetics and neural mechanisms of macroscale white matter connectome development from childhood to adolescence.

Compressed primary-to-transmodal gradient is accompanied with subcortical alterations and linked to neurotransmitters and cellular signatures in major depressive disorder

Major depressive disorder (MDD) has been shown to involve widespread changes in low-level sensorimotor and higher-level cognitive functions. Recent research found that a primary-to-transmodal gradient could capture a cortical hierarchical organization ranging from perception and action to cognition in healthy subjects, but a prominent gradient dysfunction in MDD patients. However, whether and how this cortical gradient is linked to subcortical impairments and whether it is reflected in the microscale neurotransmitter systems and cell type-specific transcriptional signatures remain largely unknown. Data were acquired from 323 MDD patients and 328 sex- and age-matched healthy controls derived from the REST-meta-MDD project, and the human brain neurotransmitter systems density maps and gene expression data were drawn from two publicly available datasets. We investigated alterations of the primary-to-transmodal gradient in MDD patients and their correlations with clinical symptoms of depression and anxiety, as well as  their paralleled subcortical impairments. The correlations between MDD-related gradient alterations and densities of the neurotransmitter systems and gene expression information were assessed, respectively. The results demonstrated that MDD patients had a compressed primary-to-transmodal gradient accompanied by paralleled alterations in subcortical regions including the caudate, amygdala, and thalamus. The case-control gradient differences were spatially correlated with the densities of the neurotransmitter systems including the serotonin and dopamine receptors, and meanwhile with gene expression enriched in astrocytes, excitatory and inhibitory neuronal cells. These findings mapped the paralleled subcortical impairments in cortical hierarchical organization and also helped us understand the possible molecular and cellular substrates of the co-occurrence of high-level cognitive impairments with low-level sensorimotor abnormalities in MDD.

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.

Macroscale Thalamic Functional Organization Disturbances and Underlying Core Cytoarchitecture in Early-Onset Schizophrenia

BACKGROUND AND HYPOTHESIS

Schizophrenia is a polygenetic mental disorder with heterogeneous positive and negative symptom constellations, and is associated with abnormal cortical connectivity. The thalamus has a coordinative role in cortical function and is key to the development of the cerebral cortex. Conversely, altered functional organization of the thalamus might relate to overarching cortical disruptions in schizophrenia, anchored in development.

STUDY DESIGN

Here, we contrasted resting-state fMRI in 86 antipsychotic-naive first-episode early-onset schizophrenia (EOS) patients and 91 typically developing controls to study whether macroscale thalamic organization is altered in EOS. Employing dimensional reduction techniques on thalamocortical functional connectome (FC), we derived lateral-medial and anterior-posterior thalamic functional axes.

STUDY RESULTS

We observed increased segregation of macroscale thalamic functional organization in EOS patients, which was related to altered thalamocortical interactions both in unimodal and transmodal networks. Using an ex vivo approximation of core-matrix cell distribution, we found that core cells particularly underlie the macroscale abnormalities in EOS patients. Moreover, the disruptions were associated with schizophrenia-related gene expression maps. Behavioral and disorder decoding analyses indicated that the macroscale hierarchy disturbances might perturb both perceptual and abstract cognitive functions and contribute to negative syndromes in patients.

CONCLUSIONS

These findings provide mechanistic evidence for disrupted thalamocortical system in schizophrenia, suggesting a unitary pathophysiological framework.

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.

Functional and microstructural plasticity following social and interoceptive mental training

The human brain supports social cognitive functions, including Theory of Mind, empathy, and compassion, through its intrinsic hierarchical organization. However, it remains unclear how the learning and refinement of social skills shapes brain function and structure. We studied if different types of social mental training induce changes in cortical function and microstructure, investigating 332 healthy adults (197 women, 20-55 years) with repeated multimodal neuroimaging and behavioral testing. Our neuroimaging approach examined longitudinal changes in cortical functional gradients and myelin-sensitive T1 relaxometry, two complementary measures of cortical hierarchical organization. We observed marked changes in intrinsic cortical function and microstructure, which varied as a function of social training content. In particular, cortical function and microstructure changed as a result of attention-mindfulness and socio-cognitive training in regions functionally associated with attention and interoception, including insular and parietal cortices. Conversely, socio-affective and socio-cognitive training resulted in differential microstructural changes in regions classically implicated in interoceptive and emotional processing, including insular and orbitofrontal areas, but did not result in functional reorganization. Notably, longitudinal changes in cortical function and microstructure predicted behavioral change in attention, compassion and perspective-taking. Our work demonstrates functional and microstructural plasticity after the training of social-interoceptive functions, and illustrates the bidirectional relationship between brain organisation and human social skills.

Connectopic mapping techniques do not reflect functional gradients in the brain

Functional gradients, in which response properties change gradually across a brain region, have been proposed as a key organising principle of the brain. Recent studies using both resting-state and natural viewing paradigms have indicated that these gradients may be reconstructed from functional connectivity patterns via "connectopic mapping" analyses. However, local connectivity patterns may be confounded by spatial autocorrelations artificially introduced during data analysis, for instance by spatial smoothing or interpolation between coordinate spaces. Here, we investigate whether such confounds can produce illusory connectopic gradients. We generated datasets comprising random white noise in subjects' functional volume spaces, then optionally applied spatial smoothing and/or interpolated the data to a different volume or surface space. Both smoothing and interpolation induced spatial autocorrelations sufficient for connectopic mapping to produce both volume- and surface-based local gradients in numerous brain regions. Furthermore, these gradients appeared highly similar to those obtained from real natural viewing data, although gradients generated from real and random data were statistically different in certain scenarios. We also reconstructed global gradients across the whole-brain - while these appeared less susceptible to artificial spatial autocorrelations, the ability to reproduce previously reported gradients was closely linked to specific features of the analysis pipeline. These results indicate that previously reported gradients identified by connectopic mapping techniques may be confounded by artificial spatial autocorrelations introduced during the analysis, and in some cases may reproduce poorly across different analysis pipelines. These findings imply that connectopic gradients need to be interpreted with caution.

Whole-brain functional gradients reveal cortical and subcortical alterations in patients with episodic migraine

Migraine is a type of headache with multiple neurological symptoms. Prior neuroimaging studies in patients with migraine based on functional magnetic resonance imaging have found regional as well as network-level alterations in brain function. Here, we expand on prior studies by establishing whole-brain functional connectivity patterns in patients with migraine using dimensionality reduction techniques. We studied functional brain connectivity in 50 patients with episodic migraine and sex- and age-matched healthy controls. Using dimensionality reduction techniques that project high-dimensional functional connectivity onto low-dimensional representations (i.e., eigenvectors), we found significant between-group differences in the eigenvectors between patients with migraine and healthy controls, particularly in the sensory/motor and limbic cortices. Furthermore, we assessed between-group differences in subcortical connectivity with subcortical weighted manifolds defined by subcortico-cortical connectivity multiplied by cortical eigenvectors and revealed significant alterations in the amygdala. Finally, leveraging supervised machine learning, we moderately predicted headache frequency using cortical and subcortical functional connectivity features, again indicating that sensory and limbic regions play a particularly important role in predicting migraine frequency. Our study confirmed that migraine is a hierarchical disease of the brain that shows alterations along the sensory-limbic axis, and therefore, the functional connectivity in these areas could be a useful marker to investigate migraine symptomatology.

Correspondence of functional connectivity gradients across human isocortex, cerebellum, and hippocampus

Gradient mapping is an important technique to summarize high dimensional biological features as low dimensional manifold representations in exploring brain structure-function relationships at various levels of the cerebral cortex. While recent studies have characterized the major gradients of functional connectivity in several brain structures using this technique, very few have systematically examined the correspondence of such gradients across structures under a common systems-level framework. Using resting-state functional magnetic resonance imaging, here we show that the organizing principles of the isocortex, and those of the cerebellum and hippocampus in relation to the isocortex, can be described using two common functional gradients. We suggest that the similarity in functional connectivity gradients across these structures can be meaningfully interpreted within a common computational framework based on the principles of predictive processing. The present results, and the specific hypotheses that they suggest, represent an important step toward an integrative account of brain function.

BrainStat: A toolbox for brain-wide statistics and multimodal feature associations

Analysis and interpretation of neuroimaging datasets has become a multidisciplinary endeavor, relying not only on statistical methods, but increasingly on associations with respect to other brain-derived features such as gene expression, histological data, and functional as well as cognitive architectures. Here, we introduce BrainStat - a toolbox for (i) univariate and multivariate linear models in volumetric and surface-based brain imaging datasets, and (ii) multidomain feature association of results with respect to spatial maps of post-mortem gene expression and histology, task-based fMRI meta-analysis, as well as resting-state fMRI motifs across several common surface templates. The combination of statistics and feature associations into a turnkey toolbox streamlines analytical processes and accelerates cross-modal research. The toolbox is implemented in both Python and MATLAB, two widely used programming languages in the neuroimaging and neuroinformatics communities. BrainStat is openly available and complemented by an expandable documentation.

Distinct patterns of cortical manifold expansion and contraction underlie human sensorimotor adaptation

Sensorimotor learning is a dynamic, systems-level process that involves the combined action of multiple neural systems distributed across the brain. Although much is known about the specialized cortical systems that support specific components of action (such as reaching), we know less about how cortical systems function in a coordinated manner to facilitate adaptive behavior. To address this gap, our study measured human brain activity using functional MRI (fMRI) while participants performed a classic sensorimotor adaptation task and used a manifold learning approach to describe how behavioral changes during adaptation relate to changes in the landscape of cortical activity. During early adaptation, areas in the parietal and premotor cortices exhibited significant contraction along the cortical manifold, which was associated with their increased covariance with regions in the higher-order association cortex, including both the default mode and fronto-parietal networks. By contrast, during Late adaptation, when visuomotor errors had been largely reduced, a significant expansion of the visual cortex along the cortical manifold was associated with its reduced covariance with the association cortex and its increased intraconnectivity. Lastly, individuals who learned more rapidly exhibited greater covariance between regions in the sensorimotor and association cortices during early adaptation. These findings are consistent with a view that sensorimotor adaptation depends on changes in the integration and segregation of neural activity across more specialized regions of the unimodal cortex with regions in the association cortex implicated in higher-order processes. More generally, they lend support to an emerging line of evidence implicating regions of the default mode network (DMN) in task-based performance.

Age differences in the functional architecture of the human brain

The intrinsic functional organization of the brain changes into older adulthood. Age differences are observed at multiple spatial scales, from global reductions in modularity and segregation of distributed brain systems, to network-specific patterns of dedifferentiation. Whether dedifferentiation reflects an inevitable, global shift in brain function with age, circumscribed, experience-dependent changes, or both, is uncertain. We employed a multimethod strategy to interrogate dedifferentiation at multiple spatial scales. Multi-echo (ME) resting-state fMRI was collected in younger (n = 181) and older (n = 120) healthy adults. Cortical parcellation sensitive to individual variation was implemented for precision functional mapping of each participant while preserving group-level parcel and network labels. ME-fMRI processing and gradient mapping identified global and macroscale network differences. Multivariate functional connectivity methods tested for microscale, edge-level differences. Older adults had lower BOLD signal dimensionality, consistent with global network dedifferentiation. Gradients were largely age-invariant. Edge-level analyses revealed discrete, network-specific dedifferentiation patterns in older adults. Visual and somatosensory regions were more integrated within the functional connectome; default and frontoparietal control network regions showed greater connectivity; and the dorsal attention network was more integrated with heteromodal regions. These findings highlight the importance of multiscale, multimethod approaches to characterize the architecture of functional brain aging.

Multiscale neural gradients reflect transdiagnostic effects of major psychiatric conditions on cortical morphology

It is increasingly recognized that multiple psychiatric conditions are underpinned by shared neural pathways, affecting similar brain systems. Here, we carried out a multiscale neural contextualization of shared alterations of cortical morphology across six major psychiatric conditions (autism spectrum disorder, attention deficit/hyperactivity disorder, major depression disorder, obsessive-compulsive disorder, bipolar disorder, and schizophrenia). Our framework cross-referenced shared morphological anomalies with respect to cortical myeloarchitecture and cytoarchitecture, as well as connectome and neurotransmitter organization. Pooling disease-related effects on MRI-based cortical thickness measures across six ENIGMA working groups, including a total of 28,546 participants (12,876 patients and 15,670 controls), we identified a cortex-wide dimension of morphological changes that described a sensory-fugal pattern, with paralimbic regions showing the most consistent alterations across conditions. The shared disease dimension was closely related to cortical gradients of microstructure as well as neurotransmitter axes, specifically cortex-wide variations in serotonin and dopamine. Multiple sensitivity analyses confirmed robustness with respect to slight variations in analytical choices. Our findings embed shared effects of common psychiatric conditions on brain structure in multiple scales of brain organization, and may provide insights into neural mechanisms of transdiagnostic vulnerability.

Adolescent development of multiscale structural wiring and functional interactions in the human connectome

Adolescence is a time of profound changes in the physical wiring and function of the brain. Here, we analyzed structural and functional brain network development in an accelerated longitudinal cohort spanning 14 to 25 y (n = 199). Core to our work was an advanced in vivo model of cortical wiring incorporating MRI features of corticocortical proximity, microstructural similarity, and white matter tractography. Longitudinal analyses assessing age-related changes in cortical wiring identified a continued differentiation of multiple corticocortical structural networks in youth. We then assessed structure-function coupling using resting-state functional MRI measures in the same participants both via cross-sectional analysis at baseline and by studying longitudinal change between baseline and follow-up scans. At baseline, regions with more similar structural wiring were more likely to be functionally coupled. Moreover, correlating longitudinal structural wiring changes with longitudinal functional connectivity reconfigurations, we found that increased structural differentiation, particularly between sensory/unimodal and default mode networks, was reflected by reduced functional interactions. These findings provide insights into adolescent development of human brain structure and function, illustrating how structural wiring interacts with the maturation of macroscale functional hierarchies.

Dissociable multi-scale patterns of development in personalized brain networks

The brain is organized into networks at multiple resolutions, or scales, yet studies of functional network development typically focus on a single scale. Here, we derive personalized functional networks across 29 scales in a large sample of youths (n = 693, ages 8-23 years) to identify multi-scale patterns of network re-organization related to neurocognitive development. We found that developmental shifts in inter-network coupling reflect and strengthen a functional hierarchy of cortical organization. Furthermore, we observed that scale-dependent effects were present in lower-order, unimodal networks, but not higher-order, transmodal networks. Finally, we found that network maturation had clear behavioral relevance: the development of coupling in unimodal and transmodal networks are dissociably related to the emergence of executive function. These results suggest that the development of functional brain networks align with and refine a hierarchy linked to cognition.

Genetic and phylogenetic uncoupling of structure and function in human transmodal cortex

Brain structure scaffolds intrinsic function, supporting cognition and ultimately behavioral flexibility. However, it remains unclear how a static, genetically controlled architecture supports flexible cognition and behavior. Here, we synthesize genetic, phylogenetic and cognitive analyses to understand how the macroscale organization of structure-function coupling across the cortex can inform its role in cognition. In humans, structure-function coupling was highest in regions of unimodal cortex and lowest in transmodal cortex, a pattern that was mirrored by a reduced alignment with heritable connectivity profiles. Structure-function uncoupling in macaques had a similar spatial distribution, but we observed an increased coupling between structure and function in association cortices relative to humans. Meta-analysis suggested regions with the least genetic control (low heritable correspondence and different across primates) are linked to social-cognition and autobiographical memory. Our findings suggest that genetic and evolutionary uncoupling of structure and function in different transmodal systems may support the emergence of complex forms of cognition.

Focal neural perturbations reshape low-dimensional trajectories of brain activity supporting cognitive performance

The emergence of distributed patterns of neural activity supporting brain functions and behavior can be understood by study of the brain's low-dimensional topology. Functional neuroimaging demonstrates that brain activity linked to adaptive behavior is constrained to low-dimensional manifolds. In human participants, we tested whether these low-dimensional constraints preserve working memory performance following local neuronal perturbations. We combined multi-session functional magnetic resonance imaging, non-invasive transcranial magnetic stimulation (TMS), and methods translated from the fields of complex systems and computational biology to assess the functional link between changes in local neural activity and the reshaping of task-related low dimensional trajectories of brain activity. We show that specific reconfigurations of low-dimensional trajectories of brain activity sustain effective working memory performance following TMS manipulation of local activity on, but not off, the space traversed by these trajectories. We highlight an association between the multi-scale changes in brain activity underpinning cognitive function.

Pubertal timing and functional neurodevelopmental alterations independently mediate the effect of family conflict on adolescent psychopathology

This study tested the hypothesis that early life adversity (ELA) heightens psychopathology risk by concurrently altering pubertal and neurodevelopmental timing, and associated gene transcription signatures. Analyses focused on threat- (family conflict/neighbourhood crime) and deprivation-related ELAs (parental inattentiveness/unmet material needs), using longitudinal data from 1514 biologically unrelated youths in the Adolescent Brain and Cognitive Development (ABCD) study. Typical developmental changes in white matter microstructure corresponded to widespread BOLD signal variability (BOLDsv) increases (linked to cell communication and biosynthesis genes) and region-specific task-related BOLDsv increases/decreases (linked to signal transduction, immune and external environmental response genes). Increasing resting-state (RS), but decreasing task-related BOLDsv predicted normative functional network segregation. Family conflict was the strongest concurrent and prospective contributor to psychopathology, while material deprivation constituted an additive risk factor. ELA-linked psychopathology was predicted by higher Time 1 threat-evoked BOLDSV (associated with axonal development, myelination, cell differentiation and signal transduction genes), reduced Time 2 RS BOLDsv (associated with cell metabolism and attention genes) and greater Time 1 to Time 2 control/attention network segregation. Earlier pubertal timing and neurodevelopmental alterations independently mediated ELA effects on psychopathology. Our results underscore the differential roles of the immediate and wider external environment(s) in concurrent and longer-term ELA consequences.

Standardizing workflows in imaging transcriptomics with the abagen toolbox

Gene expression fundamentally shapes the structural and functional architecture of the human brain. Open-access transcriptomic datasets like the Allen Human Brain Atlas provide an unprecedented ability to examine these mechanisms in vivo; however, a lack of standardization across research groups has given rise to myriad processing pipelines for using these data. Here, we develop the abagen toolbox, an open-access software package for working with transcriptomic data, and use it to examine how methodological variability influences the outcomes of research using the Allen Human Brain Atlas. Applying three prototypical analyses to the outputs of 750,000 unique processing pipelines, we find that choice of pipeline has a large impact on research findings, with parameters commonly varied in the literature influencing correlations between derived gene expression and other imaging phenotypes by as much as ρ ≥ 1.0. Our results further reveal an ordering of parameter importance, with processing steps that influence gene normalization yielding the greatest impact on downstream statistical inferences and conclusions. The presented work and the development of the abagen toolbox lay the foundation for more standardized and systematic research in imaging transcriptomics, and will help to advance future understanding of the influence of gene expression in the human brain.

Early and late development of hub connectivity in the human brain

Human brain networks undergo pronounced changes during development. The emergence of highly connected hub regions that can support integrated brain function is central to this maturational process, with these areas undergoing a particularly protracted period of development that extends into adulthood. The location of cortical network hubs emerges early but connections to and from hubs continue to strengthen throughout childhood and adolescence. Patterns of functional coupling in cortical association hubs are immature and incomplete at birth, but gradually strengthen during development. Early establishment of hub connectivity may provide a stable substrate that is refined by changes in tissue organization and microstructure, resulting in the emergence of complex functional dynamics by adulthood.

Neurodevelopment of the association cortices: Patterns, mechanisms, and implications for psychopathology

The human brain undergoes a prolonged period of cortical development that spans multiple decades. During childhood and adolescence, cortical development progresses from lower-order, primary and unimodal cortices with sensory and motor functions to higher-order, transmodal association cortices subserving executive, socioemotional, and mentalizing functions. The spatiotemporal patterning of cortical maturation thus proceeds in a hierarchical manner, conforming to an evolutionarily rooted, sensorimotor-to-association axis of cortical organization. This developmental program has been characterized by data derived from multimodal human neuroimaging and is linked to the hierarchical unfolding of plasticity-related neurobiological events. Critically, this developmental program serves to enhance feature variation between lower-order and higher-order regions, thus endowing the brain's association cortices with unique functional properties. However, accumulating evidence suggests that protracted plasticity within late-maturing association cortices, which represents a defining feature of the human developmental program, also confers risk for diverse developmental psychopathologies.

Inter-individual body mass variations relate to fractionated functional brain hierarchies

Variations in body mass index (BMI) have been suggested to relate to atypical brain organization, yet connectome-level substrates of BMI and their neurobiological underpinnings remain unclear. Studying 325 healthy young adults, we examined associations between functional connectivity and inter-individual BMI variations. We utilized non-linear connectome manifold learning techniques to represent macroscale functional organization along continuous hierarchical axes that dissociate low level and higher order brain systems. We observed an increased differentiation between unimodal and heteromodal association networks in individuals with higher BMI, indicative of a disrupted modular architecture and hierarchy of the brain. Transcriptomic decoding and gene enrichment analyses identified genes previously implicated in genome-wide associations to BMI and specific cortical, striatal, and cerebellar cell types. These findings illustrate functional connectome substrates of BMI variations in healthy young adults and point to potential molecular associations.