Structural covariance networks are coupled to expression of genes enriched in supragranular layers of the human cortex

Accelerated Cortical Thinning in Schizophrenia Is Associated With Rare and Common Predisposing Variation to Schizophrenia and Neurodevelopmental Disorders

BACKGROUND

Schizophrenia is a highly heritable disorder characterized by increased cortical thinning throughout the life span. Studies have reported a shared genetic basis between schizophrenia and cortical thickness. However, no genes whose expression is related to abnormal cortical thinning in schizophrenia have been identified.

METHODS

We conducted linear mixed models to estimate the rates of accelerated cortical thinning across 68 regions from the Desikan-Killiany atlas in individuals with schizophrenia compared with healthy control participants from a large longitudinal sample (ncases = 169 and ncontrols = 298, ages 16-70 years). We studied the correlation between gene expression data from the Allen Human Brain Atlas and accelerated thinning estimates across cortical regions. Finally, we explored the functional and genetic underpinnings of the genes that contribute most to accelerated thinning.

RESULTS

We found a global pattern of accelerated cortical thinning in individuals with schizophrenia compared with healthy control participants. Genes underexpressed in cortical regions that exhibit this accelerated thinning were downregulated in several psychiatric disorders and were enriched for both common and rare disrupting variation for schizophrenia and neurodevelopmental disorders. In contrast, none of these enrichments were observed for baseline cross-sectional cortical thickness differences.

CONCLUSIONS

Our findings suggest that accelerated cortical thinning, rather than cortical thickness alone, serves as an informative phenotype for neurodevelopmental disruptions in schizophrenia. We highlight the genetic and transcriptomic correlates of this accelerated cortical thinning, emphasizing the need for future longitudinal studies to elucidate the role of genetic variation and the temporal-spatial dynamics of gene expression in brain development and aging in schizophrenia.

Voxel-based texture similarity networks reveal individual variability and correlate with biological ontologies

The human brain is organized as a complex, hierarchical network. However, the structural covariance patterns among brain regions and the underlying biological substrates of such covariance networks remain to be clarified. The present study proposed a novel individualized structural covariance network termed voxel-based texture similarity networks (vTSNs) based on 76 refined voxel-based textural features derived from structural magnetic resonance images. Validated in three independent longitudinal healthy cohorts (40, 23, and 60 healthy participants, respectively) with two common brain atlases, we found that the vTSN could robustly resolve inter-subject variability with high test-retest reliability. In contrast to the regional-based texture similarity networks (rTSNs) that calculate radiomic features based on region-of-interest information, vTSNs had higher inter- and intra-subject variability ratios and test-retest reliability in connectivity strength and network topological properties. Moreover, the Spearman correlation indicated a stronger association of the gene expression similarity network (GESN) with vTSNs than with rTSNs (vTSN: r = 0.600, rTSN: r = 0.433, z = 39.784, P < 0.001). Hierarchical clustering identified 3 vTSN subnets with differential association patterns with 13 coexpression modules, 16 neurotransmitters, 7 electrophysiology, 4 metabolism, and 2 large-scale structural and 4 functional organization maps. Moreover, these subnets had unique biological hierarchical organization from the subcortex-limbic system to the ventral neocortex and then to the dorsal neocortex. Based on 424 unrelated, qualified healthy subjects from the Human Connectome Project, we found that vTSNs could sensitively represent sex differences, especially for connections in the subcortex-limbic system and between the subcortex-limbic system and the ventral neocortex. Moreover, a multivariate variance component model revealed that vTSNs could explain a significant proportion of inter-subject behavioral variance in cognition (80.0 %) and motor functions (63.4 %). Finally, using 494 healthy adults (aged 19-80 years old) from the Southwest University Adult Lifespan Dataset, the Spearman correlation identified a significant association between aging and vTSN strength, especially within the subcortex-limbic system and between the subcortex-limbic system and the dorsal neocortex. In summary, our proposed vTSN is robust in uncovering individual variability and neurobiological brain processes, which can serve as biologically plausible measures for linking biological processes and human behavior.

Neural basis underlying the association between thought control ability and happiness: The moderating role of the amygdala

Thought control ability (TCA) plays an important role in individuals' health and happiness. Previous studies demonstrated that TCA was closely conceptually associated with happiness. However, empirical research supporting this relationship was limited. In addition, the neural basis underlying TCA and how this neural basis influences the relationship between TCA and happiness remain unexplored. In the present study, the voxel-based morphometry (VBM) method was adopted to investigate the neuroanatomical basis of TCA in 314 healthy subjects. The behavioral results revealed a significant positive association between TCA and happiness. On the neural level, there was a significant negative correlation between TCA and the gray matter density (GMD) of the bilateral amygdala. Split-half validation analysis revealed similar results, further confirming the stability of the VBM analysis findings. Furthermore, gray matter covariance network and graph theoretical analyses showed positive association between TCA and both the node degree and node strength of the amygdala. Moderation analysis revealed that the GMD of the amygdala moderated the relationship between TCA and happiness. Specifically, the positive association between TCA and self-perceived happiness was stronger in subjects with a lower GMD of the amygdala. The present study indicated the neural basis underlying the association between TCA and happiness and offered a method of improving individual well-being.

Cortical gene expression architecture links healthy neurodevelopment to the imaging, transcriptomics and genetics of autism and schizophrenia

Human brain organization involves the coordinated expression of thousands of genes. For example, the first principal component (C1) of cortical transcription identifies a hierarchy from sensorimotor to association regions. In this study, optimized processing of the Allen Human Brain Atlas revealed two new components of cortical gene expression architecture, C2 and C3, which are distinctively enriched for neuronal, metabolic and immune processes, specific cell types and cytoarchitectonics, and genetic variants associated with intelligence. Using additional datasets (PsychENCODE, Allen Cell Atlas and BrainSpan), we found that C1-C3 represent generalizable transcriptional programs that are coordinated within cells and differentially phased during fetal and postnatal development. Autism spectrum disorder and schizophrenia were specifically associated with C1/C2 and C3, respectively, across neuroimaging, differential expression and genome-wide association studies. Evidence converged especially in support of C3 as a normative transcriptional program for adolescent brain development, which can lead to atypical supragranular cortical connectivity in people at high genetic risk for schizophrenia.

Structural network topologies are associated with deep brain stimulation outcomes in Meige syndrome

Deep brain stimulation (DBS) is an effective therapy for Meige syndrome (MS). However, the DBS efficacy varies across MS patients and the factors contributing to the variable responses remain enigmatic. We aim to explain the difference in DBS efficacy from a network perspective. We collected preoperative T1-weighted MRI images of 76 MS patients who received DBS in our center. According to the symptomatic improvement rates, all MS patients were divided into two groups: the high improvement group (HIG) and the low improvement group (LIG). We constructed group-level structural covariance networks in each group and compared the graph-based topological properties and interregional connections between groups. Subsequent functional annotation and correlation analyses were also conducted. The results indicated that HIG showed a higher clustering coefficient, longer characteristic path length, lower small-world index, and lower global efficiency compared with LIG. Different nodal betweennesses and degrees between groups were mainly identified in the precuneus, sensorimotor cortex, and subcortical nuclei, among which the gray matter volume of the left precentral gyrus and left thalamus were positively correlated with the symptomatic improvement rates. Moreover, HIG had enhanced interregional connections within the somatomotor network and between the somatomotor network and default-mode network relative to LIG. We concluded that the high and low DBS responders have notable differences in large-scale network architectures. Our study sheds light on the structural network underpinnings of varying DBS responses in MS patients.

Exploring morphological similarity and randomness in Alzheimer's disease using adjacent grey matter voxel-based structural analysis

BACKGROUND

Alzheimer's disease is characterized by large-scale structural changes in a specific pattern. Recent studies developed morphological similarity networks constructed by brain regions similar in structural features to represent brain structural organization. However, few studies have used local morphological properties to explore inter-regional structural similarity in Alzheimer's disease.

METHODS

Here, we sourced T1-weighted MRI images of 342 cognitively normal participants and 276 individuals with Alzheimer's disease from the Alzheimer's Disease Neuroimaging Initiative database. The relationships of grey matter intensity between adjacent voxels were defined and converted to the structural pattern indices. We conducted the information-based similarity method to evaluate the structural similarity of structural pattern organization between brain regions. Besides, we examined the structural randomness on brain regions. Finally, the relationship between the structural randomness and cognitive performance of individuals with Alzheimer's disease was assessed by stepwise regression.

RESULTS

Compared to cognitively normal participants, individuals with Alzheimer's disease showed significant structural pattern changes in the bilateral posterior cingulate gyrus, hippocampus, and olfactory cortex. Additionally, individuals with Alzheimer's disease showed that the bilateral insula had decreased inter-regional structural similarity with frontal regions, while the bilateral hippocampus had increased inter-regional structural similarity with temporal and subcortical regions. For the structural randomness, we found significant decreases in the temporal and subcortical areas and significant increases in the occipital and frontal regions. The regression analysis showed that the structural randomness of five brain regions was correlated with the Mini-Mental State Examination scores of individuals with Alzheimer's disease.

CONCLUSIONS

Our study suggested that individuals with Alzheimer's disease alter micro-structural patterns and morphological similarity with the insula and hippocampus. Structural randomness of individuals with Alzheimer's disease changed in temporal, frontal, and occipital brain regions. Morphological similarity and randomness provide valuable insight into brain structural organization in Alzheimer's disease.

Transcriptomic Similarity Informs Neuromorphic Deviations in Depression Biotypes

BACKGROUND

Major depressive disorder (MDD) is complicated by population heterogeneity, motivating the investigation of biotypes through imaging-derived phenotypes. However, neuromorphic heterogeneity in MDD remains unclear, and how the correlated gene expression (CGE) connectome constrains these neuromorphic anomalies in MDD biotypes has not yet been studied.

METHODS

Here, we related cortical thickness deviations in MDD biotypes to a pattern of CGE connectome. Cortical thickness was estimated from 3-dimensional T1-weighted magnetic resonance images in 2 independent cohorts (discovery cohort: N = 425; replication cohort: N = 217). The transcriptional activity was measured according to Allen Human Brain Atlas. A density peak-based clustering algorithm was used to identify MDD biotypes.

RESULTS

We found that patients with MDD were clustered into 2 replicated biotypes based on single-patient regional deviations from healthy control participants across 2 datasets. Biotype 1 mainly exhibited cortical thinning across the brain, whereas biotype 2 mainly showed cortical thickening in the brain. Using brainwide gene expression data, we found that deviations of transcriptionally connected neighbors predicted regional deviation for both biotypes. Furthermore, putative CGE-informed epicenters of biotype 1 were concentrated on the cognitive control circuit, whereas biotype 2 epicenters were located in the social perception circuit. The patterns of epicenter likelihood were separately associated with depression- and anxiety-response maps, suggesting that epicenters of MDD biotypes may be associated with clinical efficacies.

CONCLUSIONS

Our findings linked the CGE connectome and neuromorphic deviations to identify distinct epicenters in MDD biotypes, providing insight into how microscale gene expressions informed MDD biotypes.

Transcriptional cartography integrates multiscale biology of the human cortex

The cerebral cortex underlies many of our unique strengths and vulnerabilities, but efforts to understand human cortical organization are challenged by reliance on incompatible measurement methods at different spatial scales. Macroscale features such as cortical folding and functional activation are accessed through spatially dense neuroimaging maps, whereas microscale cellular and molecular features are typically measured with sparse postmortem sampling. Here, we integrate these distinct windows on brain organization by building upon existing postmortem data to impute, validate, and analyze a library of spatially dense neuroimaging-like maps of human cortical gene expression. These maps allow spatially unbiased discovery of cortical zones with extreme transcriptional profiles or unusually rapid transcriptional change which index distinct microstructure and predict neuroimaging measures of cortical folding and functional activation. Modules of spatially coexpressed genes define a family of canonical expression maps that integrate diverse spatial scales and temporal epochs of human brain organization - ranging from protein-protein interactions to large-scale systems for cognitive processing. These module maps also parse neuropsychiatric risk genes into subsets which tag distinct cyto-laminar features and differentially predict the location of altered cortical anatomy and gene expression in patients. Taken together, the methods, resources, and findings described here advance our understanding of human cortical organization and offer flexible bridges to connect scientific fields operating at different spatial scales of human brain research.

Combining multimodal connectivity information improves modelling of pathology spread in Alzheimer's disease

Cortical atrophy and aggregates of misfolded tau proteins are key hallmarks of Alzheimer's disease. Computational models that simulate the propagation of pathogens between connected brain regions have been used to elucidate mechanistic information about the spread of these disease biomarkers, such as disease epicentres and spreading rates. However, the connectomes that are used as substrates for these models are known to contain modality-specific false positive and false negative connections, influenced by the biases inherent to the different methods for estimating connections in the brain. In this work, we compare five types of connectomes for modelling both tau and atrophy patterns with the network diffusion model, which are validated against tau PET and structural MRI data from individuals with either mild cognitive impairment or dementia. We then test the hypothesis that a joint connectome, with combined information from different modalities, provides an improved substrate for the model. We find that a combination of multimodal information helps the model to capture observed patterns of tau deposition and atrophy better than any single modality. This is validated with data from independent datasets. Overall, our findings suggest that combining connectivity measures into a single connectome can mitigate some of the biases inherent to each modality and facilitate more accurate models of pathology spread, thus aiding our ability to understand disease mechanisms, and providing insight into the complementary information contained in different measures of brain connectivity.

Similarity between structural and proxy estimates of brain connectivity

Functional magnetic resonance and diffusion weighted imaging have so far made a major contribution to delineation of the brain connectome at the macroscale. While functional connectivity (FC) was shown to be related to structural connectivity (SC) to a certain degree, their spatial overlap is unknown. Even less clear are relations of SC with estimates of connectivity from inter-subject covariance of regional F18-fluorodeoxyglucose uptake (FDGcov) and grey matter volume (GMVcov). Here, we asked to what extent SC underlies three proxy estimates of brain connectivity: FC, FDGcov and GMVcov. Simultaneous PET/MR acquisitions were performed in 56 healthy middle-aged individuals. Similarity between four networks was assessed using Spearman correlation and convergence ratio (CR), a measure of spatial overlap. Spearman correlation coefficient was 0.27 for SC-FC, 0.40 for SC-FDGcov, and 0.15 for SC-GMVcov. Mean CRs were 51% for SC-FC, 48% for SC-FDGcov, and 37% for SC-GMVcov. These results proved to be reproducible and robust against image processing steps. In sum, we found a relevant similarity of SC with FC and FDGcov, while GMVcov consistently showed the weakest similarity. These findings indicate that white matter tracts underlie FDGcov to a similar degree as FC, supporting FDGcov as estimate of functional brain connectivity.

Mapping individual structural covariance network in development brain with dynamic time warping

A conspicuous property of brain development or maturity is coupled with coordinated or synchronized brain structural co-variation. However, there is still a lack of effective approach to map individual structural covariance network. Here, we developed a novel individual structural covariance network method using dynamic time warping algorithm and applied it to delineate developmental trajectories of topological organizations of structural covariance network from childhood to early adulthood with a large sample of 655 individuals from Human Connectome Project-Development dataset. We found that the individual structural covariance network exhibited small-worldness property and the network global topological characteristics including small-worldness, global efficiency, local efficiency, and modularity linearly increase with age while the shortest path length linearly decreases with age. The nodal topological properties including betweenness and degree increased with age in language and emotion regulation related brain areas, while it decreased with age mainly in visual cortex, sensorimotor area, and hippocampus. Moreover, the topological attributes of structural covariance network as features could predict the age of each individual. Taken together, our results demonstrate that dynamic time warping can effectively map individual structural covariance network to uncover the developmental trajectories of network topology, which may facilitate future investigations to establish the links of structural co-variations with respect to cognition and disease vulnerability.

Abnormalities of structural covariance of insular subregions in drug-naïve OCD patients

The insula plays a significant role in the neural mechanisms of obsessive-compulsive disorder. Previous studies have identified functional and structural abnormalities in insula in obsessive-compulsive disorder patients. The predictive coding model in the context of interoception can explain the psychological and neuropathological manifestations observed in obsessive-compulsive disorder. The model is based on the degree of laminar differentiation of cerebral cortex. The interindividual differences in a local measure of brain structure often covary with interindividual differences in other brain regions. We investigated the anatomical network involving the insula in a drug-naïve obsessive-compulsive disorder sample. We recruited 58 obsessive-compulsive disorder patients and 84 matched health controls. The cortical thickness covariance maps between groups were compared at each vertex. We also evaluated the modulation of Yale-Brown Obsessive-Compulsive Scale scores and obsessive-compulsive disorder duration on thickness covariance. Our findings indicated that the thickness covariance seeded from granular and dysgranular insula are different compared with controls. The duration and severity of obsessive-compulsive disorder can modulate the thickness covariance of granular and dysgranular insula with posterior cingulate cortex and rostral anterior cingulate cortex. Our results revealed aberrant insular structural characteristics and cortical thickness covariance in obsessive-compulsive disorder patients, contributing to a better understanding of the involvement of insula in the pathological mechanisms underlying obsessive-compulsive disorder.

Organization of thalamocortical structural covariance and a corresponding 3D atlas of the mouse thalamus

For information from sensory organs to be processed by the brain, it is usually passed to appropriate areas of the cerebral cortex. Almost all of this information passes through the thalamus, a relay structure that reciprocally connects to the vast majority of the cortex. The thalamus facilitates this information transfer through a set of thalamocortical connections that vary in cellular structure, molecular profiles, innervation patterns, and firing rates. Additionally, corticothalamic connections allow for intracortical information transfer through the thalamus. These efferent and afferent connections between the thalamus and cortex have been the focus of many studies, and the importance of cortical connectivity in defining thalamus anatomy is demonstrated by multiple studies that parcellate the thalamus based on cortical connectivity profiles. Here, we examine correlated morphological variation between the thalamus and cortex, or thalamocortical structural covariance. For each voxel in the thalamus as a seed, we construct a cortical structural covariance map that represents correlated cortical volume variation, and examine whether high structural covariance is observed in cortical areas that are functionally relevant to the seed. Then, using these cortical structural covariance maps as features, we subdivide the thalamus into six non-overlapping regions (clusters of voxels), and assess whether cortical structural covariance is associated with cortical connectivity that specifically originates from these regions. We show that cortical structural covariance is high in areas of the cortex that are functionally related to the seed voxel, cortical structural covariance varies along cortical depth, and sharp transitions in cortical structural covariance profiles are observed when varying seed locations in the thalamus. Subdividing the thalamus based on structural covariance, we additionally demonstrate that the six thalamic clusters of voxels stratify cortical structural covariance along the dorsal-ventral, medial-lateral, and anterior-posterior axes. These cluster-associated structural covariance patterns are prominently detected in cortical regions innervated by fibers projecting out of their related thalamic subdivisions. Together, these results advance our understanding of how the thalamus and the cortex couple in their volumes. Our results indicate that these volume correlations reflect functional organization and structural connectivity, and further provides a novel segmentation of the mouse thalamus that can be used to examine thalamic structural variation and thalamocortical structural covariation in disease models.

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.

The Neurogenetics of Functional Connectivity Alterations in Autism: Insights From Subtyping in 657 Individuals

BACKGROUND

There is little consensus and controversial evidence on anatomical alterations in the brains of people with autism spectrum disorder (ASD), due in part to the large heterogeneity present in ASD, which in turn is a major drawback for developing therapies. One strategy to characterize this heterogeneity in ASD is to cluster large-scale functional brain connectivity profiles.

METHODS

A subtyping approach based on consensus clustering of functional brain connectivity patterns was applied to a population of 657 autistic individuals with quality-assured neuroimaging data. We then used high-resolution gene transcriptomic data to characterize the molecular mechanism behind each subtype by performing enrichment analysis of the set of genes showing a high spatial similarity with the profiles of functional connectivity alterations between each subtype and a group of typically developing control participants.

RESULTS

Two major stable subtypes were found: subtype 1 exhibited hypoconnectivity (less average connectivity than typically developing control participants) and subtype 2, hyperconnectivity. The 2 subtypes did not differ in structural imaging metrics in any of the analyzed regions (68 cortical and 14 subcortical) or in any of the behavioral scores (including IQ, Autism Diagnostic Interview, and Autism Diagnostic Observation Schedule). Finally, only subtype 2, comprising about 43% of ASD participants, led to significant enrichments after multiple testing corrections. Notably, the dominant enrichment corresponded to excitation/inhibition imbalance, a leading well-known primary mechanism in the pathophysiology of ASD.

CONCLUSIONS

Our results support a link between excitation/inhibition imbalance and functional connectivity alterations, but only in one ASD subtype, overall characterized by brain hyperconnectivity and major alterations in somatomotor and default mode networks.

The biological role of local and global fMRI BOLD signal variability in human brain organization

Variability drives the organization and behavior of complex systems, including the human brain. Understanding the variability of brain signals is thus necessary to broaden our window into brain function and behavior. Few empirical investigations of macroscale brain signal variability have yet been undertaken, given the difficulty in separating biological sources of variance from artefactual noise. Here, we characterize the temporal variability of the most predominant macroscale brain signal, the fMRI BOLD signal, and systematically investigate its statistical, topographical and neurobiological properties. We contrast fMRI acquisition protocols, and integrate across histology, microstructure, transcriptomics, neurotransmitter receptor and metabolic data, fMRI static connectivity, and empirical and simulated magnetoencephalography data. We show that BOLD signal variability represents a spatially heterogeneous, central property of multi-scale multi-modal brain organization, distinct from noise. Our work establishes the biological relevance of BOLD signal variability and provides a lens on brain stochasticity across spatial and temporal scales.

Network connectivity and local transcriptomic vulnerability underpin cortical atrophy progression in Parkinson's disease

Parkinson's disease pathology is hypothesized to spread through the brain via axonal connections between regions and is further modulated by local vulnerabilities within those regions. The resulting changes to brain morphology have previously been demonstrated in both prodromal and de novo Parkinson's disease patients. However, it remains unclear whether the pattern of atrophy progression in Parkinson's disease over time is similarly explained by network-based spreading and local vulnerability. We address this gap by mapping the trajectory of cortical atrophy rates in a large, multi-centre cohort of Parkinson's disease patients and relate this atrophy progression pattern to network architecture and gene expression profiles. Across 4-year follow-up visits, increased atrophy rates were observed in posterior, temporal, and superior frontal cortices. We demonstrated that this progression pattern was shaped by network connectivity. Regional atrophy rates were strongly related to atrophy rates across structurally and functionally connected regions. We also found that atrophy progression was associated with specific gene expression profiles. The genes whose spatial distribution in the brain was most related to atrophy rate were those enriched for mitochondrial and metabolic function. Taken together, our findings demonstrate that both global and local brain features influence vulnerability to neurodegeneration in Parkinson's disease.

Cortical structural differences following repeated ayahuasca use hold molecular signatures

Introduction

Serotonergic psychedelics such as ayahuasca are reported to promote both structural and functional neural plasticity via partial 5-HT2A agonism. However, little is known about how these molecular mechanisms may extend to repeated psychedelic administration in humans, let alone neuroanatomy. While early evidence suggests localised changes to cortical thickness in long-term ayahuasca users, it is unknown how such findings may be reflected by large-scale anatomical brain networks comprising cytoarchitecturally complex regions.

Methods

Here, we examined the relationship between cortical gene expression markers of psychedelic action and brain morphometric change following repeated ayahuasca usage, using high-field 7 Tesla neuroimaging data derived from 24 members of an ayahuasca-using church (Santo Daime) and case-matched controls.

Results

Using a morphometric similarity network (MSN) analysis, repeated ayahuasca use was associated with a spatially distributed cortical patterning of both structural differentiation in sensorimotor areas and de-differentiation in transmodal areas. Cortical MSN remodelling was found to be spatially correlated with dysregulation of 5-HT2A gene expression as well as a broader set of genes encoding target receptors pertinent to ayahuasca's effects. Furthermore, these associations were similarly interrelated with altered gene expression of specific transcriptional factors and immediate early genes previously identified in preclinical assays as relevant to psychedelic-induced neuroplasticity.

Conclusion

Taken together, these findings provide preliminary evidence that the molecular mechanisms of psychedelic action may scale up to a macroscale level of brain organisation in vivo. Closer attention to the role of cortical transcriptomics in structural-functional coupling may help account for the behavioural differences observed in experienced psychedelic users.

Genomic Structural Equation Modeling Reveals Latent Phenotypes in the Human Cortex with Distinct Genetic Architecture

Genetic contributions to human cortical structure manifest pervasive pleiotropy. This pleiotropy may be harnessed to identify unique genetically-informed parcellations of the cortex that are neurobiologically distinct from functional, cytoarchitectural, or other cortical parcellation schemes. We investigated genetic pleiotropy by applying genomic structural equation modeling (SEM) to map the genetic architecture of cortical surface area (SA) and cortical thickness (CT) for the 34 brain regions recently reported in the ENIGMA cortical GWAS. Genomic SEM uses the empirical genetic covariance estimated from GWAS summary statistics with LD score regression (LDSC) to discover factors underlying genetic covariance, which we are denoting genetically informed brain networks (GIBNs). Genomic SEM can fit a multivariate GWAS from summary statistics for each of the GIBNs, which can subsequently be used for LD score regression (LDSC). We found the best-fitting model of cortical SA identified 6 GIBNs and CT identified 4 GIBNs. The multivariate GWASs of these GIBNs identified 74 genome-wide significant (GWS) loci (p<5×10-8), including many previously implicated in neuroimaging phenotypes, behavioral traits, and psychiatric conditions. LDSC of GIBN GWASs found that SA-derived GIBNs had a positive genetic correlation with bipolar disorder (BPD), and cannabis use disorder, indicating genetic predisposition to a larger SA in the specific GIBN is associated with greater genetic risk of these disorders. A negative genetic correlation was observed with attention deficit hyperactivity disorder (ADHD), major depressive disorder (MDD), and insomnia, indicating genetic predisposition to a larger SA in the specific GIBN is associated with lower genetic risk of these disorders. CT GIBNs displayed a negative genetic correlation with alcohol dependence. Jointly modeling the genetic architecture of complex traits and investigating multivariate genetic links across phenotypes offers a new vantage point for mapping the cortex into genetically informed networks.

Neurophysiological signatures of cortical micro-architecture

Systematic spatial variation in micro-architecture is observed across the cortex. These micro-architectural gradients are reflected in neural activity, which can be captured by neurophysiological time-series. How spontaneous neurophysiological dynamics are organized across the cortex and how they arise from heterogeneous cortical micro-architecture remains unknown. Here we extensively profile regional neurophysiological dynamics across the human brain by estimating over 6800 time-series features from the resting state magnetoencephalography (MEG) signal. We then map regional time-series profiles to a comprehensive multi-modal, multi-scale atlas of cortical micro-architecture, including microstructure, metabolism, neurotransmitter receptors, cell types and laminar differentiation. We find that the dominant axis of neurophysiological dynamics reflects characteristics of power spectrum density and linear correlation structure of the signal, emphasizing the importance of conventional features of electromagnetic dynamics while identifying additional informative features that have traditionally received less attention. Moreover, spatial variation in neurophysiological dynamics is co-localized with multiple micro-architectural features, including gene expression gradients, intracortical myelin, neurotransmitter receptors and transporters, and oxygen and glucose metabolism. Collectively, this work opens new avenues for studying the anatomical basis of neural activity.

Integrating multimodal and multiscale connectivity blueprints of the human cerebral cortex in health and disease

The brain is composed of disparate neural populations that communicate and interact with one another. Although fiber bundles, similarities in molecular architecture, and synchronized neural activity all reflect how brain regions potentially interact with one another, a comprehensive study of how all these interregional relationships jointly reflect brain structure and function remains missing. Here, we systematically integrate 7 multimodal, multiscale types of interregional similarity ("connectivity modes") derived from gene expression, neurotransmitter receptor density, cellular morphology, glucose metabolism, haemodynamic activity, and electrophysiology in humans. We first show that for all connectivity modes, feature similarity decreases with distance and increases when regions are structurally connected. Next, we show that connectivity modes exhibit unique and diverse connection patterns, hub profiles, spatial gradients, and modular organization. Throughout, we observe a consistent primacy of molecular connectivity modes-namely correlated gene expression and receptor similarity-that map onto multiple phenomena, including the rich club and patterns of abnormal cortical thickness across 13 neurological, psychiatric, and neurodevelopmental disorders. Finally, to construct a single multimodal wiring map of the human cortex, we fuse all 7 connectivity modes and show that the fused network maps onto major organizational features of the cortex including structural connectivity, intrinsic functional networks, and cytoarchitectonic classes. Altogether, this work contributes to the integrative study of interregional relationships in the human cerebral cortex.

Multimodal and multiscale evidence for network-based cortical thinning in major depressive disorder

BACKGROUND

Major depressive disorder (MDD) is associated with widespread, irregular cortical thickness (CT) reductions across the brain. However, little is known regarding mechanisms that govern spatial distribution of the reductions.

METHODS

We combined multimodal MRI and genetic, cytoarchitectonic and chemoarchitectonic data to examine structural covariance, functional synchronization, gene co-expression, cytoarchitectonic similarity and chemoarchitectonic covariance between regions atrophied in MDD.

RESULTS

Regions atrophied in MDD were associated with significantly higher structural covariance, functional synchronization, gene co-expression and chemoarchitectonic covariance. These results were robust against methodological variations in brain parcellation and null model, reproducible in patients and controls, and independent of age at onset of MDD. Despite no significant differences in the cytoarchitectonic similarity, MDD-related CT reductions were susceptible to specific cytoarchitectonic class of association cortex. Further, we found that nodal shortest path lengths to disease epicenters derived from structural (right supramarginal gyrus) and chemoarchitectonic covariance (right sulcus intermedius primus) networks of healthy brains were correlated with the extent to which a region was atrophied in MDD, supporting the transneuronal spread hypothesis that regions closer to the epicenters are more susceptible to MDD. Finally, we showed that structural covariance and functional synchronization among regions atrophied in MDD were mainly related to genes enriched in metabolic and membrane-related processes, driven by genes in excitatory neurons, and associated with specific neurotransmitter transporters and receptors.

CONCLUSIONS

Altogether, our findings provide empirical evidence for and genetic and molecular insights into connectivity-constrained CT thinning in MDD.

The Cellular Underpinnings of the Human Cortical Connectome

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

Hippocampal anterior- posterior shift in childhood and adolescence

Hippocampal-cortical networks play an important role in neurocognitive development. Applying the method of Connectivity-Based Parcellation (CBP) on hippocampal-cortical structural covariance (SC) networks computed from T1-weighted magnetic resonance images, we examined how the hippocampus differentiates into subregions during childhood and adolescence (N = 1105, 6-18 years). In late childhood, the hippocampus mainly differentiated along the anterior-posterior axis similar to previous reported functional differentiation patterns of the hippocampus. In contrast, in adolescence a differentiation along the medial-lateral axis was evident, reminiscent of the cytoarchitectonic division into cornu ammonis and subiculum. Further meta-analytical characterization of hippocampal subregions in terms of related structural co-maturation networks, behavioural and gene profiling suggested that the hippocampal head is related to higher order functions (e.g. language, theory of mind, autobiographical memory) in late childhood morphologically co-varying with almost the whole brain. In early adolescence but not in childhood, posterior subicular SC networks were associated with action-oriented and reward systems. The findings point to late childhood as an important developmental period for hippocampal head morphology and to early adolescence as a crucial period for hippocampal integration into action- and reward-oriented cognition. The latter may constitute a developmental feature that conveys increased propensity for addictive disorders.

A structural magnetic resonance imaging review of clinical motor outcomes from deep brain stimulation in movement disorders

Patients with movement disorders treated by deep brain stimulation do not always achieve successful therapeutic alleviation of motor symptoms, even in cases where surgery is without complications. Magnetic resonance imaging (MRI) offers methods to investigate structural brain-related factors that may be predictive of clinical motor outcomes. This review aimed to identify features which have been associated with variability in clinical post-operative motor outcomes in patients with Parkinson's disease, dystonia, and essential tremor from structural MRI modalities. We performed a literature search for articles published between 1 January 2000 and 1 April 2022 and identified 5197 articles. Following screening through our inclusion criteria, we identified 60 total studies (39 = Parkinson's disease, 11 = dystonia syndromes and 10 = essential tremor). The review captured a range of structural MRI methods and analysis techniques used to identify factors related to clinical post-operative motor outcomes from deep brain stimulation. Morphometric markers, including volume and cortical thickness were commonly identified in studies focused on patients with Parkinson's disease and dystonia syndromes. Reduced metrics in basal ganglia, sensorimotor and frontal regions showed frequent associations with reduced motor outcomes. Increased structural connectivity to subcortical nuclei, sensorimotor and frontal regions was also associated with greater motor outcomes. In patients with tremor, increased structural connectivity to the cerebellum and cortical motor regions showed high prevalence across studies for greater clinical motor outcomes. In addition, we highlight conceptual issues for studies assessing clinical response with structural MRI and discuss future approaches towards optimizing individualized therapeutic benefits. Although quantitative MRI markers are in their infancy for clinical purposes in movement disorder treatments, structural features obtained from MRI offer the powerful potential to identify candidates who are more likely to benefit from deep brain stimulation and provide insight into the complexity of disorder pathophysiology.

Mapping Human Brain Pathways: Challenges and Opportunities in the Integration of Scales

The human brain is composed of a complex web of pathways. Diffusion magnetic resonance (MR) tractography is a neuroimaging technique that relies on the principle of diffusion to reconstruct brain pathways. Its tractography is broadly applicable to a range of problems as it is amenable for study in individuals of any age and from any species. However, it is well known that this technique can generate biologically implausible pathways, especially in regions of the brain where multiple fibers cross. This review highlights potential misconnections in two cortico-cortical association pathways with a focus on the aslant tract and inferior frontal occipital fasciculus. The lack of alternative methods to validate observations from diffusion MR tractography means there is a need to develop new integrative approaches to trace human brain pathways. This review discusses integrative approaches in neuroimaging, anatomical, and transcriptional variation as having much potential to trace the evolution of human brain pathways.

Opposite effects of positive and negative symptoms on resting-state brain networks in schizophrenia

Schizophrenia is a severe psychotic disorder characterized by positive and negative symptoms, but their neural bases remain poorly understood. Here, we utilized a nested-spectral partition (NSP) approach to detect hierarchical modules in resting-state brain functional networks in schizophrenia patients and healthy controls, and we studied dynamic transitions of segregation and integration as well as their relationships with clinical symptoms. Schizophrenia brains showed a more stable integrating process and a more variable segregating process, thus maintaining higher segregation, especially in the limbic system. Hallucinations were associated with higher integration in attention systems, and avolition was related to a more variable segregating process in default-mode network (DMN) and control systems. In a machine-learning model, NSP-based features outperformed graph measures at predicting positive and negative symptoms. Multivariate analysis confirmed that positive and negative symptoms had opposite effects on dynamic segregation and integration of brain networks. Gene ontology analysis revealed that the effect of negative symptoms was related to autistic, aggressive and violent behavior; the effect of positive symptoms was associated with hyperammonemia and acidosis; and the interaction effect was correlated with abnormal motor function. Our findings could contribute to the development of more accurate diagnostic criteria for positive and negative symptoms in schizophrenia.

Connectional Hierarchy in Human Brain Revealed by Individual Variability of Functional Network Edges

The human cerebral cortex is connected by intricate inter-areal wiring at the macroscale. The cortical hierarchy from primary sensorimotor to higher-order association areas is a unifying organizational principle across various neurobiological properties; however, previous studies have not clarified whether the connections between cortical regions exhibit a similar hierarchical pattern. Here, we identify a connectional hierarchy indexed by inter-individual variability of functional connectivity edges, which continuously progresses along a hierarchical gradient from within-network connections to between-network edges connecting sensorimotor and association networks. We found that this connectional hierarchy of variability aligns with both hemodynamic and electromagnetic connectivity strength and is constrained by structural connectivity strength. Moreover, the patterning of connectional hierarchy is related to inter-regional similarity in transcriptional and neurotransmitter receptor profiles. Using the Neurosynth cognitive atlas and cortical vulnerability maps in 13 brain disorders, we found that the connectional hierarchy of variability is associated with similarity networks of cognitive relevance and that of disorder vulnerability. Finally, we found that the prominence of this hierarchical gradient of connectivity variability declines during youth. Together, our results reveal a novel hierarchal organizational principle at the connectional level that links multimodal and multiscale human connectomes to individual variability in functional connectivity.

Transcriptomic and connectomic correlates of differential spatial patterning among gliomas

Unravelling the complex events driving grade-specific spatial distribution of brain tumour occurrence requires rich datasets from both healthy individuals and patients. Here, we combined open-access data from The Cancer Genome Atlas, the UK Biobank and the Allen Brain Human Atlas to disentangle how the different spatial occurrences of glioblastoma multiforme and low-grade gliomas are linked to brain network features and the normative transcriptional profiles of brain regions. From MRI of brain tumour patients, we first constructed a grade-related frequency map of the regional occurrence of low-grade gliomas and the more aggressive glioblastoma multiforme. Using associated mRNA transcription data, we derived a set of differential gene expressions from glioblastoma multiforme and low-grade gliomas tissues of the same patients. By combining the resulting values with normative gene expressions from post-mortem brain tissue, we constructed a grade-related expression map indicating which brain regions express genes dysregulated in aggressive gliomas. Additionally, we derived an expression map of genes previously associated with tumour subtypes in a genome-wide association study (tumour-related genes). There were significant associations between grade-related frequency, grade-related expression and tumour-related expression maps, as well as functional brain network features (specifically, nodal strength and participation coefficient) that are implicated in neurological and psychiatric disorders. These findings identify brain network dynamics and transcriptomic signatures as key factors in regional vulnerability for glioblastoma multiforme and low-grade glioma occurrence, placing primary brain tumours within a well established framework of neurological and psychiatric cortical alterations.

Morphometric dis-similarity between cortical and subcortical areas underlies cognitive function and psychiatric symptomatology: a preadolescence study from ABCD

Preadolescence is a critical period characterized by dramatic morphological changes and accelerated cortico-subcortical development. Moreover, the coordinated development of cortical and subcortical regions underlies the emerging cognitive functions during this period. Deviations in this maturational coordination may underlie various psychiatric disorders that begin during preadolescence, but to date these deviations remain largely uncharted. We constructed a comprehensive whole-brain morphometric similarity network (MSN) from 17 neuroimaging modalities in a large preadolescence sample (N = 8908) from Adolescent Brain Cognitive Development (ABCD) study and investigated its association with 10 cognitive subscales and 27 psychiatric subscales or diagnoses. Based on the MSNs, each brain was clustered into five modules with distinct cytoarchitecture and evolutionary relevance. While morphometric correlation was positive within modules, it was negative between modules, especially between isocortical and paralimbic/subcortical modules; this developmental dissimilarity was genetically linked to synapse and neurogenesis. The cortico-subcortical dissimilarity becomes more pronounced longitudinally in healthy children, reflecting developmental differentiation of segregated cytoarchitectonic areas. Higher cortico-subcortical dissimilarity (between the isocortical and paralimbic/subcortical modules) were related to better cognitive performance. In comparison, children with poor modular differentiation between cortex and subcortex displayed higher burden of externalizing and internalizing symptoms. These results highlighted cortical-subcortical morphometric dissimilarity as a dynamic maturational marker of cognitive and psychiatric status during the preadolescent stage and provided insights into brain development.

Toward Best Practices for Imaging Transcriptomics of the Human Brain

Modern brainwide transcriptional atlases provide unprecedented opportunities for investigating the molecular correlates of brain organization, as quantified using noninvasive neuroimaging. However, integrating neuroimaging data with transcriptomic measures is not straightforward, and careful consideration is required to make valid inferences. In this article, we review recent work exploring how various methodological choices affect 3 main phases of imaging transcriptomic analyses, including 1) processing of transcriptional atlas data; 2) relating transcriptional measures to independently derived neuroimaging phenotypes; and 3) evaluating the functional implications of identified associations through gene enrichment analyses. Our aim is to facilitate the development of standardized and reproducible approaches for this rapidly growing field. We identify sources of methodological variability, key choices that can affect findings, and considerations for mitigating false positive and/or spurious results. Finally, we provide an overview of freely available open-source toolboxes implementing current best-practice procedures across all 3 analysis phases.

Neurophysiological signatures of cortical micro-architecture

Systematic spatial variation in micro-architecture is observed across the cortex. These micro-architectural gradients are reflected in neural activity, which can be captured by neurophysiological time-series. How spontaneous neurophysiological dynamics are organized across the cortex and how they arise from heterogeneous cortical micro-architecture remains unknown. Here we extensively profile regional neurophysiological dynamics across the human brain by estimating over 6 800 timeseries features from the resting state magnetoencephalography (MEG) signal. We then map regional time-series profiles to a comprehensive multi-modal, multi-scale atlas of cortical micro-architecture, including microstructure, metabolism, neurotransmitter receptors, cell types and laminar differentiation. We find that the dominant axis of neurophysiological dynamics reflects characteristics of power spectrum density and linear correlation structure of the signal, emphasizing the importance of conventional features of electromagnetic dynamics while identifying additional informative features that have traditionally received less attention. Moreover, spatial variation in neurophysiological dynamics is colocalized with multiple micro-architectural features, including genomic gradients, intracortical myelin, neurotransmitter receptors and transporters, and oxygen and glucose metabolism. Collectively, this work opens new avenues for studying the anatomical basis of neural activity.

Asymptomatic carotid stenosis is associated with both edge and network reconfigurations identified by single-subject cortical thickness networks

Background and purpose

Patients with asymptomatic carotid stenosis, even without stroke, are at high risk for cognitive impairment, and the neuroanatomical basis remains unclear. Using a novel edge-centric structural connectivity (eSC) analysis from individualized single-subject cortical thickness networks, we aimed to examine eSC and network measures in severe (> 70%) asymptomatic carotid stenosis (SACS).

Methods

Twenty-four SACS patients and 24 demographically- and comorbidities-matched controls were included, and structural MRI and multidomain cognitive data were acquired. Individual eSC was estimated via the Manhattan distances of pairwise cortical thickness histograms.

Results

In the eSC analysis, SACS patients showed longer interhemispheric but shorter intrahemispheric Manhattan distances seeding from left lateral temporal regions; in network analysis the SACS patients had a decreased system segregation paralleling with white matter hyperintensity burden and recall memory. Further network-based statistic analysis identified several eSC and subgraph features centred around the Perisylvian regions that predicted silent lesion load and cognitive tests.

Conclusion

We conclude that SACS exhibits abnormal eSC and a less-optimized trade-off between physical cost and network segregation, providing a reference and perspective for identifying high-risk individuals.

Correspondence between gene expression and neurotransmitter receptor and transporter density in the human brain

Neurotransmitter receptors modulate signaling between neurons. Thus, neurotransmitter receptors and transporters play a key role in shaping brain function. Due to the lack of comprehensive neurotransmitter receptor/transporter density datasets, microarray gene expression measuring mRNA transcripts is often used as a proxy for receptor densities. In the present report, we comprehensively test the spatial correlation between gene expression and protein density for a total of 27 neurotransmitter receptors, receptor binding-sites, and transporters across 9 different neurotransmitter systems, using both PET and autoradiography radioligand-based imaging modalities. We find poor spatial correspondences between gene expression and density for all neurotransmitter receptors and transporters except four single-protein metabotropic receptors (5-HT_1A, CB₁, D₂, and MOR). These expression-density associations are related to gene differential stability and can vary between cortical and subcortical structures. Altogether, we recommend using direct measures of receptor and transporter density when relating neurotransmitter systems to brain structure and function.

Neurofilament light-associated connectivity in young-adult Huntington's disease is related to neuronal genes

Upregulation of functional network connectivity in the presence of structural degeneration is seen in the premanifest stages of Huntington's disease (preHD) 10-15 years from clinical diagnosis. However, whether widespread network connectivity changes are seen in gene carriers much further from onset has yet to be explored. We characterized functional network connectivity throughout the brain and related it to a measure of disease pathology burden (CSF neurofilament light, NfL) and measures of structural connectivity in asymptomatic gene carriers, on average 24 years from onset. We related these measurements to estimates of cortical and subcortical gene expression. We found no overall differences in functional (or structural) connectivity anywhere in the brain comparing control and preHD participants. However, increased functional connectivity, particularly between posterior cortical areas, correlated with increasing CSF NfL level in preHD participants. Using the Allen Human Brain Atlas and expression-weighted cell-type enrichment analysis, we demonstrated that this functional connectivity upregulation occurred in cortical regions associated with regional expression of genes specific to neuronal cells. This relationship was validated using single-nucleus RNAseq data from post-mortem Huntington's disease and control brains showing enrichment of neuronal-specific genes that are differentially expressed in Huntington's disease. Functional brain networks in asymptomatic preHD gene carriers very far from disease onset show evidence of upregulated connectivity correlating with increased disease burden. These changes occur among brain areas that show regional expression of genes specific to neuronal GABAergic and glutamatergic cells.

Brain Structural Covariance Network in Asperger Syndrome Differs From Those in Autism Spectrum Disorder and Healthy Controls

Introduction

Autism is a heterogeneous neurodevelopmental disorder associated with social, cognitive and behavioral impairments. These impairments are often reported along with alteration of the brain structure such as abnormal changes in the grey matter (GM) density. However, it is not yet clear whether these changes could be used to differentiate various subtypes of autism spectrum disorder (ASD).

Method

We compared the regional changes of GM density in ASD, Asperger's Syndrome (AS) individuals and a group of healthy controls (HC). In addition to regional changes itself, the amount of GM density changes in one region as compared to other brain regions was also calculated. We hypothesized that this structural covariance network could differentiate the AS individuals from the ASD and HC groups. Therefore, statistical analysis was performed on the MRI data of 70 male subjects including 26 ASD (age=14-50, IQ=92-132), 16 AS (age=7-58, IQ=93-133) and 28 HC (age=9-39, IQ=95-144).

Result

The one-way ANOVA on the GM density of 116 anatomically separated regions showed significant differences among the groups. The pattern of structural covariance network indicated that covariation of GM density between the brain regions is altered in ASD.

Conclusion

This changed structural covariance could be considered as a reason for less efficient segregation and integration of information in the brain that could lead to cognitive dysfunctions in autism. We hope these findings could improve our understanding about the pathobiology of autism and may pave the way towards a more effective intervention paradigm.

Remodeling of the Cortical Structural Connectome in Posttraumatic Stress Disorder: Results From the ENIGMA-PGC Posttraumatic Stress Disorder Consortium

BACKGROUND

Posttraumatic stress disorder (PTSD) is accompanied by disrupted cortical neuroanatomy. We investigated alteration in covariance of structural networks associated with PTSD in regions that demonstrate the case-control differences in cortical thickness (CT) and surface area (SA).

METHODS

Neuroimaging and clinical data were aggregated from 29 research sites in >1300 PTSD cases and >2000 trauma-exposed control subjects (ages 6.2-85.2 years) by the ENIGMA-PGC (Enhancing Neuro Imaging Genetics through Meta Analysis-Psychiatric Genomics Consortium) PTSD working group. Cortical regions in the network were rank ordered by the effect size of PTSD-related cortical differences in CT and SA. The top-n (n = 2-148) regions with the largest effect size for PTSD > non-PTSD formed hypertrophic networks, the largest effect size for PTSD < non-PTSD formed atrophic networks, and the smallest effect size of between-group differences formed stable networks. The mean structural covariance (SC) of a given n-region network was the average of all positive pairwise correlations and was compared with the mean SC of 5000 randomly generated n-region networks.

RESULTS

Patients with PTSD, relative to non-PTSD control subjects, exhibited lower mean SC in CT-based and SA-based atrophic networks. Comorbid depression, sex, and age modulated covariance differences of PTSD-related structural networks.

CONCLUSIONS

Covariance of structural networks based on CT and cortical SA are affected by PTSD and further modulated by comorbid depression, sex, and age. The SC networks that are perturbed in PTSD comport with converging evidence from resting-state functional connectivity networks and networks affected by inflammatory processes and stress hormones in PTSD.

Sex differences in the functional topography of association networks in youth

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

Local molecular and global connectomic contributions to cross-disorder cortical abnormalities

Numerous brain disorders demonstrate structural brain abnormalities, which are thought to arise from molecular perturbations or connectome miswiring. The unique and shared contributions of these molecular and connectomic vulnerabilities to brain disorders remain unknown, and has yet to be studied in a single multi-disorder framework. Using MRI morphometry from the ENIGMA consortium, we construct maps of cortical abnormalities for thirteen neurodevelopmental, neurological, and psychiatric disorders from N = 21,000 participants and N = 26,000 controls, collected using a harmonised processing protocol. We systematically compare cortical maps to multiple micro-architectural measures, including gene expression, neurotransmitter density, metabolism, and myelination (molecular vulnerability), as well as global connectomic measures including number of connections, centrality, and connection diversity (connectomic vulnerability). We find a relationship between molecular vulnerability and white-matter architecture that drives cortical disorder profiles. Local attributes, particularly neurotransmitter receptor profiles, constitute the best predictors of both disorder-specific cortical morphology and cross-disorder similarity. Finally, we find that cross-disorder abnormalities are consistently subtended by a small subset of network epicentres in bilateral sensory-motor, inferior temporal lobe, precuneus, and superior parietal cortex. Collectively, our results highlight how local molecular attributes and global connectivity jointly shape cross-disorder cortical abnormalities.

Human electromagnetic and haemodynamic networks systematically converge in unimodal cortex and diverge in transmodal cortex

Whole-brain neural communication is typically estimated from statistical associations among electromagnetic or haemodynamic time-series. The relationship between functional network architectures recovered from these 2 types of neural activity remains unknown. Here, we map electromagnetic networks (measured using magnetoencephalography (MEG)) to haemodynamic networks (measured using functional magnetic resonance imaging (fMRI)). We find that the relationship between the 2 modalities is regionally heterogeneous and systematically follows the cortical hierarchy, with close correspondence in unimodal cortex and poor correspondence in transmodal cortex. Comparison with the BigBrain histological atlas reveals that electromagnetic-haemodynamic coupling is driven by laminar differentiation and neuron density, suggesting that the mapping between the 2 modalities can be explained by cytoarchitectural variation. Importantly, haemodynamic connectivity cannot be explained by electromagnetic activity in a single frequency band, but rather arises from the mixing of multiple neurophysiological rhythms. Correspondence between the two is largely driven by MEG functional connectivity at the beta (15 to 29 Hz) frequency band. Collectively, these findings demonstrate highly organized but only partly overlapping patterns of connectivity in MEG and fMRI functional networks, opening fundamentally new avenues for studying the relationship between cortical microarchitecture and multimodal connectivity patterns.

Morphometric features of drug-resistant essential tremor and recovery after stereotactic radiosurgical thalamotomy

Essential tremor (ET) is the most common movement disorder. Its neural underpinnings remain unclear. Here, we quantified structural covariance between cortical thickness (CT), surface area (SA), and mean curvature (MC) estimates in patients with ET before and 1 year after ventro-intermediate nucleus stereotactic radiosurgical thalamotomy, and contrasted the observed patterns with those from matched healthy controls. For SA, complex rearrangements within a network of motion-related brain areas characterized patients with ET. This was complemented by MC alterations revolving around the left middle temporal cortex and the disappearance of positive-valued covariance across both modalities in the right fusiform gyrus. Recovery following thalamotomy involved MC readjustments in frontal brain centers, the amygdala, and the insula, capturing nonmotor characteristics of the disease. The appearance of negative-valued CT covariance between the left parahippocampal gyrus and hippocampus was another recovery mechanism involving high-level visual areas. This was complemented by the appearance of negative-valued CT/MC covariance, and positive-valued SA/MC covariance, in the right inferior temporal cortex and bilateral fusiform gyrus. Our results demonstrate that different morphometric properties provide complementary information to understand ET, and that their statistical cross-dependences are also valuable. They pinpoint several anatomical features of the disease and highlight routes of recovery following thalamotomy.

Multivariate GWAS of psychiatric disorders and their cardinal symptoms reveal two dimensions of cross-cutting genetic liabilities

Understanding which biological pathways are specific versus general across diagnostic categories and levels of symptom severity is critical to improving nosology and treatment of psychopathology. Here, we combine transdiagnostic and dimensional approaches to genetic discovery for the first time, conducting a novel multivariate genome-wide association study of eight psychiatric symptoms and disorders broadly related to mood disturbance and psychosis. We identify two transdiagnostic genetic liabilities that distinguish between common forms of psychopathology versus rarer forms of serious mental illness. Biological annotation revealed divergent genetic architectures that differentially implicated prenatal neurodevelopment and neuronal function and regulation. These findings inform psychiatric nosology and biological models of psychopathology, as they suggest that the severity of mood and psychotic symptoms present in serious mental illness may reflect a difference in kind rather than merely in degree.

A synergistic core for human brain evolution and cognition

How does the organization of neural information processing enable humans' sophisticated cognition? Here we decompose functional interactions between brain regions into synergistic and redundant components, revealing their distinct information-processing roles. Combining functional and structural neuroimaging with meta-analytic results, we demonstrate that redundant interactions are predominantly associated with structurally coupled, modular sensorimotor processing. Synergistic interactions instead support integrative processes and complex cognition across higher-order brain networks. The human brain leverages synergistic information to a greater extent than nonhuman primates, with high-synergy association cortices exhibiting the highest degree of evolutionary cortical expansion. Synaptic density mapping from positron emission tomography and convergent molecular and metabolic evidence demonstrate that synergistic interactions are supported by receptor diversity and human-accelerated genes underpinning synaptic function. This information-resolved approach provides analytic tools to disentangle information integration from coupling, enabling richer, more accurate interpretations of functional connectivity, and illuminating how the human neurocognitive architecture navigates the trade-off between robustness and integration.

Graph Theoretical Analysis of Structural Covariance Reveals the Relevance of Visuospatial and Attentional Areas in Essential Tremor Recovery After Stereotactic Radiosurgical Thalamotomy

Essential tremor (ET) is the most common movement disorder. Its pathophysiology is only partially understood. Here, we leveraged graph theoretical analysis on structural covariance patterns quantified from morphometric estimates for cortical thickness, surface area, and mean curvature in patients with ET before and one year after (to account for delayed clinical effect) ventro-intermediate nucleus (Vim) stereotactic radiosurgical thalamotomy. We further contrasted the observed patterns with those from matched healthy controls (HCs). Significant group differences at the level of individual morphometric properties were specific to mean curvature and the post-/pre-thalamotomy contrast, evidencing brain plasticity at the level of the targeted left thalamus, and of low-level visual, high-level visuospatial and attentional areas implicated in the dorsal visual stream. The introduction of cross-correlational analysis across pairs of morphometric properties strengthened the presence of dorsal visual stream readjustments following thalamotomy, as cortical thickness in the right lingual gyrus, bilateral rostral middle frontal gyrus, and left pre-central gyrus was interrelated with mean curvature in the rest of the brain. Overall, our results position mean curvature as the most relevant morphometric feature to understand brain plasticity in drug-resistant ET patients following Vim thalamotomy. They also highlight the importance of examining not only individual features, but also their interactions, to gain insight into the routes of recovery following intervention.

In-scanner head motion and structural covariance networks

In-scanner head motion systematically reduces estimated regional gray matter volumes obtained from structural brain MRI. Here, we investigate how head motion affects structural covariance networks that are derived from regional gray matter volumetric estimates. We acquired motion-affected and low-motion whole brain T1-weighted MRI in 29 healthy adult subjects and estimated relative regional gray matter volumes using a voxel-based morphometry approach. Structural covariance network analyses were undertaken while systematically increasing the number of included motion-affected scans. We demonstrate that the standard deviation in regional gray matter estimates increases as the number of motion-affected scans increases. This increases pairwise correlations between regions, a key determinant for construction of structural covariance networks. We further demonstrate that head motion systematically alters graph theoretic metrics derived from these networks. Finally, we present evidence that weighting correlations using image quality metrics can mitigate the effects of head motion. Our findings suggest that in-scanner head motion is a source of error that violates the assumption that structural covariance networks reflect neuroanatomical connectivity between brain regions. Results of structural covariance studies should be interpreted with caution, particularly when subject groups are likely to move their heads in the scanner.

Predictive Biomarkers for Antipsychotic Treatment Response in Early Phase of Schizophrenia: Multi-Omic Measures Linking Subcortical Covariant Network, Transcriptomic Signatures, and Peripheral Epigenetics

Background

Volumetric alterations of subcortical structures as predictors of antipsychotic treatment response have been previously corroborated, but less is known about whether their morphological covariance relates to treatment outcome and is driven by gene expression and epigenetic modifications.

Methods

Subcortical volumetric covariance was analyzed by using baseline T1-weighted magnetic resonance imaging (MRI) in 38 healthy controls and 38 drug-naïve first-episode schizophrenia patients. Patients were treated with 8-week risperidone monotherapy and divided into responder and non-responder groups according to the Remission in Schizophrenia Working Group (RSWG). We utilized partial least squares (PLS) regression to examine the spatial associations between gene expression of subcortical structures from a publicly available transcriptomic dataset and between-group variances of structural covariance. The peripheral DNA methylation (DNAm) status of a gene of interest (GOI), overlapping between genes detected in the PLS and 108 schizophrenia candidate gene loci previously reported, was examined in parallel with MRI scanning.

Results

In the psychotic symptom dimension, non-responders had a higher baseline structural covariance in the putamen-hippocampus-pallidum-accumbens pathway compared with responders. For disorganized symptoms, significant differences in baseline structural covariant connections were found in the putamen-hippocampus-pallidum-thalamus circuit between the two subgroups. The imaging variances related to psychotic symptom response were spatially related to the expression of genes enriched in neurobiological processes and dopaminergic pathways. The DNAm of GOI demonstrated significant associations with patients' improvement of psychotic symptoms.

Conclusion

Baseline subcortical structural covariance and peripheral DNAm may relate to antipsychotic treatment response. Phenotypic variations in subcortical connectome related to psychotic symptom response may be transcriptomically and epigenetically underlaid. This study defines a roadmap for future studies investigating multimodal imaging epigenetic biomarkers for treatment response in schizophrenia.

Transcriptomic and cellular decoding of functional brain connectivity changes reveal regional brain vulnerability to pro- and anti-inflammatory therapies

BACKGROUND

Systemic inflammation induces acute changes in mood, motivation and cognition that closely resemble those observed in depressed individuals. However, the mechanistic pathways linking peripheral inflammation to depression-like psychopathology via intermediate effects on brain function remain incompletely understood.

METHODS

We combined data from 30 patients initiating interferon-α treatment for Hepatitis-C and 20 anti-tumour necrosis factor (TNF) therapy for inflammatory arthritis and used resting-state functional magnetic resonance imaging to investigate acute effects of each treatment on regional global brain connectivity (GBC). We leveraged transcriptomic data from the Allen Human Brain Atlas to uncover potential biological and cellular pathways underpinning regional vulnerability to GBC changes induced by each treatment.

RESULTS

Interferon-α and anti-TNF therapies both produced differential small-to-medium sized decreases in regional GBC. However, these were observed within distinct brain regions and the regional patterns of GBC changes induced by each treatment did not correlate suggesting independent underlying processes. Further, the spatial distribution of these differential GBC decreases could be captured by multivariate patterns of constitutive regional expression of genes respectively related to: i) neuroinflammation and glial cells; and ii) glutamatergic neurotransmission and neurons. The extent to which each participant expressed patterns of GBC changes aligning with these patterns of transcriptomic vulnerability also correlated with both acute treatment-induced changes in interleukin-6 (IL-6) and, for Interferon-α, longer-term treatment-associated changes in depressive symptoms.

CONCLUSIONS

Together, we present two transcriptomic models separately linking regional vulnerability to the acute effects of interferon-α and anti-TNF treatments on brain function to glial neuroinflammation and glutamatergic neurotransmission. These findings generate hypotheses about two potential brain mechanisms through which bidirectional changes in peripheral inflammation may contribute to the development/resolution of psychopathology.

Disrupted default mode network and executive control network are associated with depression severity on structural network

OBJECTIVE

Major depressive disorder (MDD) is a psychiatric disorder with a relatively limited response to treatment. It is necessary to better understand the neuroanatomical mechanisms of structural networks.

METHODS

The current study recruited 181 first-onset, untreated adult MDD patients: slight MDD (SD, N = 23), moderate MDD (MD, N = 77), Heavy MDD (HD, N = 81) groups; along with a healthy control group (HC, N = 81) with matched general clinical data. FreeSurfer was used to preprocess T1 images for gray matter volume (GMV), and the default mode network (DMN) and the execution control network (ECN) were analyzed by structural covariance network (SCN).

RESULTS

Present study found that the GMV of brain regions reduced with the severity of the disease. Specifically, the GMV of the left anterior cingulate gyrus (ACC.L) is negatively correlated with MDD severity. In addition, the SCN connectivity of the whole-brain network increases with the increase of severity in MDD. ACC.L is a key brain region with increased connectivity between the left orbitofrontal in DMN and between the right orbitofrontal in ECN, which leads to damage to the balance of neural circuits.

CONCLUSIONS

Patients with smaller GMV of ACC.L are more likely to develop severe MDD, and as a key region in both networks which have distinct structural network models in DMN and ECN. MDD patients with different severity have different neuroimaging changes in DMN and ECN.

Edges in brain networks: Contributions to models of structure and function

Network models describe the brain as sets of nodes and edges that represent its distributed organization. So far, most discoveries in network neuroscience have prioritized insights that highlight distinct groupings and specialized functional contributions of network nodes. Importantly, these functional contributions are determined and expressed by the web of their interrelationships, formed by network edges. Here, we underscore the important contributions made by brain network edges for understanding distributed brain organization. Different types of edges represent different types of relationships, including connectivity and similarity among nodes. Adopting a specific definition of edges can fundamentally alter how we analyze and interpret a brain network. Furthermore, edges can associate into collectives and higher order arrangements, describe time series, and form edge communities that provide insights into brain network topology complementary to the traditional node-centric perspective. Focusing on the edges, and the higher order or dynamic information they can provide, discloses previously underappreciated aspects of structural and functional network organization.