Imaging Transcriptomics of Brain Disorders

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.

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

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

Association of Gene Expression and Tremor Network Structure

BACKGROUND

Transcriptomic changes in the essential tremor (ET)-associated cerebello-thalamo-cortical "tremor network" and their association to brain structure have not been investigated.

OBJECTIVE

The aim was to characterize molecular changes associated with network-level imaging-derived phenotypes (IDP) found in ET.

METHODS

We performed an imaging-transcriptomic study in British adults using imaging-genome-wide association study summary statistics (UK Biobank "BIG40" cohort; n = 33,224, aged 40-69 years). We imputed imaging-transcriptomic associations for 184 IDPs and analyzed functional enrichment of gene modules and aggregate network-level phenotypes. Validation was performed in cerebellar-tissue RNA-sequencing data from ET patients and controls (n = 55).

RESULTS

Among 237,896 individual predicted gene expression levels for 6063 unique genes/transcripts, we detected 2269 genome-wide significant associations (Bonferroni P < 2.102e-7, 0.95%). These were concentrated in intracellular volume fraction measures of white matter pathways and in genes with putative links to tremor (MAPT, ARL17A, KANSL1, SPPL2C, LRRC37A4P, PLEKHM1, and FMNL1). Whole-tremor-network cortical thickness was associated with a gene module linked to mitochondrial organization and protein quality control (r = 0.91, P = 2e-70), whereas white-gray T1-weighted magnetic resonance imaging (MRI) contrast in the tremor network was associated with a gene module linked to sphingolipid synthesis and ethanolamine metabolism (r = -0.90, P = 2e-68). Imputed association effect sizes and RNA-sequencing log-fold change in the validation dataset were significantly correlated for cerebellar peduncular diffusion MRI phenotypes, and there was a close overlap of significant associations between both datasets for gray matter phenotypes (χ2 = 6.40, P = 0.006).

CONCLUSIONS

The identified genes and processes are potential treatment targets for ET, and our results help characterize molecular changes that could in future be used for patient treatment selection or prognosis prediction. © 2024 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Cortical gradient perturbation in attention deficit hyperactivity disorder correlates with neurotransmitter-, cell type-specific and chromosome- transcriptomic signatures

AIMS

This study aimed to illuminate the neuropathological landscape of attention deficit hyperactivity disorder (ADHD) by a multiscale macro-micro-molecular perspective from in vivo neuroimaging data.

METHODS

The "ADHD-200 initiative" repository provided multi-site high-quality resting-state functional connectivity (rsfc-) neuroimaging for ADHD children and matched typically developing (TD) cohort. Diffusion mapping embedding model to derive the functional connectome gradient detecting biologically plausible neural pattern was built, and the multivariate partial least square method to uncover the enrichment of neurotransmitomic, cellular and chromosomal gradient-transcriptional signatures of AHBA enrichment and meta-analytic decoding.

RESULTS

Compared to TD, ADHD children presented connectopic cortical gradient perturbations in almost all the cognition-involved brain macroscale networks (all pBH <0.001), but not in the brain global topology. As an intermediate phenotypic variant, such gradient perturbation was spatially enriched into distributions of GABAA/BZ and 5-HT2A receptors (all pBH <0.01) and co-varied with genetic transcriptional expressions (e.g. DYDC2, ATOH7, all pBH <0.01), associated with phenotypic variants in episodic memory and emotional regulations. Enrichment models demonstrated such gradient-transcriptional variants indicated the risk of both cell-specific and chromosome- dysfunctions, especially in enriched expression of oligodendrocyte precursors and endothelial cells (all pperm <0.05) as well enrichment into chromosome 18, 19 and X (pperm <0.05).

CONCLUSIONS

Our findings bridged brain macroscale neuropathological patterns to microscale/cellular biological architectures for ADHD children, demonstrating the neurobiologically pathological mechanism of ADHD into the genetic and molecular variants in GABA and 5-HT systems as well brain-derived enrichment of specific cellular/chromosomal expressions.

Homotopic functional connectivity disruptions in schizophrenia and their associated gene expression

It has been revealed that abnormal voxel-mirrored homotopic connectivity (VMHC) is present in patients with schizophrenia, yet there are inconsistencies in the relevant findings. Moreover, little is known about their association with brain gene expression profiles. In this study, transcription-neuroimaging association analyses using gene expression data from Allen Human Brain Atlas and case-control VMHC differences from both the discovery (meta-analysis, including 9 studies with a total of 386 patients and 357 controls) and replication (separate group-level comparisons within two datasets, including a total of 258 patients and 287 controls) phases were performed to identify genes associated with VMHC alterations. Enrichment analyses were conducted to characterize the biological functions and specific expression of identified genes, and Neurosynth decoding analysis was performed to examine the correlation between cognitive-related processes and VMHC alterations in schizophrenia. In the discovery and replication phases, patients with schizophrenia exhibited consistent VMHC changes compared to controls, which were correlated with a series of cognitive-related processes; meta-regression analysis revealed that illness duration was negatively correlated with VMHC abnormalities in the cerebellum and postcentral/precentral gyrus. The abnormal VMHC patterns were stably correlated with 1287 genes enriched for fundamental biological processes like regulation of cell communication, nervous system development, and cell communication. In addition, these genes were overexpressed in astrocytes and immune cells, enriched in extensive cortical regions and wide developmental time windows. The present findings may contribute to a more comprehensive understanding of the molecular mechanisms underlying VMHC alterations in patients with schizophrenia.

Strength of spatial correlation between gray matter connectivity and patterns of proto-oncogene and neural network construction gene expression is associated with diffuse glioma survival

Introduction

Like other forms of neuropathology, gliomas appear to spread along neural pathways. Accordingly, our group and others have previously shown that brain network connectivity is highly predictive of glioma survival. In this study, we aimed to examine the molecular mechanisms of this relationship via imaging transcriptomics.

Methods

We retrospectively obtained presurgical, T1-weighted MRI datasets from 669 adult patients, newly diagnosed with diffuse glioma. We measured brain connectivity using gray matter networks and coregistered these data with a transcriptomic brain atlas to determine the spatial co-localization between brain connectivity and expression patterns for 14 proto-oncogenes and 3 neural network construction genes.

Results

We found that all 17 genes were significantly co-localized with brain connectivity (p < 0.03, corrected). The strength of co-localization was highly predictive of overall survival in a cross-validated Cox Proportional Hazards model (mean area under the curve, AUC = 0.68 +/- 0.01) and significantly (p < 0.001) more so for a random forest survival model (mean AUC = 0.97 +/- 0.06). Bayesian network analysis demonstrated direct and indirect causal relationships among gene-brain co-localizations and survival. Gene ontology analysis showed that metabolic processes were overexpressed when spatial co-localization between brain connectivity and gene transcription was highest (p < 0.001). Drug-gene interaction analysis identified 84 potential candidate therapies based on our findings.

Discussion

Our findings provide novel insights regarding how gene-brain connectivity interactions may affect glioma survival.

Distinctive whole-brain cell types predict tissue damage patterns in thirteen neurodegenerative conditions

For over a century, brain research narrative has mainly centered on neuron cells. Accordingly, most neurodegenerative studies focus on neuronal dysfunction and their selective vulnerability, while we lack comprehensive analyses of other major cell types' contribution. By unifying spatial gene expression, structural MRI, and cell deconvolution, here we describe how the human brain distribution of canonical cell types extensively predicts tissue damage in 13 neurodegenerative conditions, including early- and late-onset Alzheimer's disease, Parkinson's disease, dementia with Lewy bodies, amyotrophic lateral sclerosis, mutations in presenilin-1, and 3 clinical variants of frontotemporal lobar degeneration (behavioral variant, semantic and non-fluent primary progressive aphasia) along with associated three-repeat and four-repeat tauopathies and TDP43 proteinopathies types A and C. We reconstructed comprehensive whole-brain reference maps of cellular abundance for six major cell types and identified characteristic axes of spatial overlapping with atrophy. Our results support the strong mediating role of non-neuronal cells, primarily microglia and astrocytes, in spatial vulnerability to tissue loss in neurodegeneration, with distinct and shared across-disorder pathomechanisms. These observations provide critical insights into the multicellular pathophysiology underlying spatiotemporal advance in neurodegeneration. Notably, they also emphasize the need to exceed the current neuro-centric view of brain diseases, supporting the imperative for cell-specific therapeutic targets in neurodegeneration.

C9orf72 gene networks in the human brain correlate with cortical thickness in C9-FTD and implicate vulnerable cell types

Introduction

A hexanucleotide repeat expansion (HRE) intronic to chromosome 9 open reading frame 72 (C9orf72) is recognized as the most common genetic cause of amyotrophic lateral sclerosis (ALS), frontotemporal dementia (FTD), and ALS-FTD. Identifying genes that show similar regional co-expression patterns to C9orf72 may help identify novel gene targets and biological mechanisms that mediate selective vulnerability to ALS and FTD pathogenesis.

Methods

We leveraged mRNA expression data in healthy brain from the Allen Human Brain Atlas to evaluate C9orf72 co-expression patterns. To do this, we correlated average C9orf72 expression values in 51 regions across different anatomical divisions (cortex, subcortex, and cerebellum) with average gene expression values for 15,633 protein-coding genes, including 54 genes known to be associated with ALS, FTD, or ALS-FTD. We then performed imaging transcriptomic analyses to evaluate whether the identified C9orf72 co-expressed genes correlated with patterns of cortical thickness in symptomatic C9orf72 pathogenic HRE carriers (n = 19) compared to controls (n = 23). Lastly, we explored whether genes with significant C9orf72 imaging transcriptomic correlations (i.e., "C9orf72 imaging transcriptomic network") were enriched in specific cell populations in the brain and enriched for specific biological and molecular pathways.

Results

A total of 2,120 genes showed an anatomical distribution of gene expression in the brain similar to C9orf72 and significantly correlated with patterns of cortical thickness in C9orf72 HRE carriers. This C9orf72 imaging transcriptomic network was differentially expressed in cell populations previously implicated in ALS and FTD, including layer 5b cells, cholinergic neurons in the spinal cord and brainstem and medium spiny neurons of the striatum, and was enriched for biological and molecular pathways associated with protein ubiquitination, autophagy, cellular response to DNA damage, endoplasmic reticulum to Golgi vesicle-mediated transport, among others.

Conclusion

Considered together, we identified a network of C9orf72 associated genes that may influence selective regional and cell-type-specific vulnerabilities in ALS/FTD.

Acute Stress Increases Striatal Connectivity With Cortical Regions Enriched for μ and κ Opioid Receptors

BACKGROUND

Understanding the neurobiological effects of stress is critical for addressing the etiology of major depressive disorder (MDD). Using a dimensional approach involving individuals with differing degree of MDD risk, we investigated 1) the effects of acute stress on cortico-cortical and subcortical-cortical functional connectivity (FC) and 2) how such effects are related to gene expression and receptor maps.

METHODS

Across 115 participants (37 control, 39 remitted MDD, 39 current MDD), we evaluated the effects of stress on FC during the Montreal Imaging Stress Task. Using partial least squares regression, we investigated genes whose expression in the Allen Human Brain Atlas was associated with anatomical patterns of stress-related FC change. Finally, we correlated stress-related FC change maps with opioid and GABAA (gamma-aminobutyric acid A) receptor distribution maps derived from positron emission tomography.

RESULTS

Results revealed robust effects of stress on global cortical connectivity, with increased global FC in frontoparietal and attentional networks and decreased global FC in the medial default mode network. Moreover, robust increases emerged in FC of the caudate, putamen, and amygdala with regions from the ventral attention/salience network, frontoparietal network, and motor networks. Such regions showed preferential expression of genes involved in cell-to-cell signaling (OPRM1, OPRK1, SST, GABRA3, GABRA5), similar to previous genetic MDD studies.

CONCLUSIONS

Acute stress altered global cortical connectivity and increased striatal connectivity with cortical regions that express genes that have previously been associated with imaging abnormalities in MDD and are rich in μ and κ opioid receptors. These findings point to overlapping circuitry underlying stress response, reward, and MDD.

Integrative omics analysis reveals epigenomic and transcriptomic signatures underlying brain structural deficits in major depressive disorder

Several lines of evidence support the involvement of transcriptomic and epigenetic mechanisms in the brain structural deficits of major depressive disorder (MDD) separately. However, research in these two areas has remained isolated. In this study, we proposed an integrative strategy that combined neuroimaging, brain-wide gene expression, and peripheral DNA methylation data to investigate the genetic basis of gray matter abnormalities in MDD. The MRI T1-weighted images and Illumina 850 K DNA methylation microarrays were obtained from 269 patients and 416 healthy controls, and brain-wide transcriptomic data were collected from Allen Human Brain Atlas. The between-group differences in gray matter volume (GMV) and differentially methylated CpG positions (DMPs) were examined. The genes with their expression patterns spatially related to GMV changes and genes with DMPs were overlapped and selected. Using principal component regression, the associations between DMPs in overlapped genes and GMV across individual patients were investigated, and the region-specific correlations between methylation status and gene expression were examined. We found significant associations between the decreased GMV and DMPs methylation status in the anterior cingulate cortex, inferior frontal cortex, and fusiform face cortex regions. These DMPs genes were primarily enriched in the neurodevelopmental and synaptic transmission process. There was a significant negative correlation between DNA methylation and gene expression in genes associated with GMV changes of the frontal cortex in MDD. Our findings suggest that GMV abnormalities in MDD may have a transcriptomic and epigenetic basis. This imaging-transcriptomic-epigenetic integrative analysis provides spatial and biological links between cortical morphological deficits and peripheral epigenetic signatures in MDD.

Strength of spatial correlation between structural brain network connectivity and brain-wide patterns of proto-oncogene and neural network construction gene expression is associated with diffuse glioma survival

Like other forms of neuropathology, gliomas appear to spread along neural pathways. Accordingly, our group and others have previously shown that brain network connectivity is highly predictive of glioma survival. In this study, we aimed to examine the molecular mechanisms of this relationship via imaging transcriptomics. We retrospectively obtained presurgical, T1-weighted MRI datasets from 669 adult patients, newly diagnosed with diffuse glioma. We measured brain connectivity using gray matter networks and coregistered these data with a transcriptomic brain atlas to determine the spatial co-localization between brain connectivity and expression patterns for 14 proto-oncogenes and 3 neural network construction genes. We found that all 17 genes were significantly co-localized with brain connectivity (p < 0.03, corrected). The strength of co-localization was highly predictive of overall survival in a cross-validated Cox Proportional Hazards model (mean area under the curve, AUC = 0.68 +/- 0.01) and significantly (p < 0.001) more so for a random forest survival model (mean AUC = 0.97 +/- 0.06). Bayesian network analysis demonstrated direct and indirect causal relationships among gene-brain co-localizations and survival. Gene ontology analysis showed that metabolic processes were overexpressed when spatial co-localization between brain connectivity and gene transcription was highest (p < 0.001). Drug-gene interaction analysis identified 84 potential candidate therapies based on our findings. Our findings provide novel insights regarding how gene-brain connectivity interactions may affect glioma survival.

Genetic topography and cortical cell loss in Huntington's disease link development and neurodegeneration

Cortical cell loss is a core feature of Huntington's disease (HD), beginning many years before clinical motor diagnosis, during the premanifest stage. However, it is unclear how genetic topography relates to cortical cell loss. Here, we explore the biological processes and cell types underlying this relationship and validate these using cell-specific post-mortem data. Eighty premanifest participants on average 15 years from disease onset and 71 controls were included. Using volumetric and diffusion MRI we extracted HD-specific whole brain maps where lower grey matter volume and higher grey matter mean diffusivity, relative to controls, were used as proxies of cortical cell loss. These maps were combined with gene expression data from the Allen Human Brain Atlas (AHBA) to investigate the biological processes relating genetic topography and cortical cell loss. Cortical cell loss was positively correlated with the expression of developmental genes (i.e. higher expression correlated with greater atrophy and increased diffusivity) and negatively correlated with the expression of synaptic and metabolic genes that have been implicated in neurodegeneration. These findings were consistent for diffusion MRI and volumetric HD-specific brain maps. As wild-type huntingtin is known to play a role in neurodevelopment, we explored the association between wild-type huntingtin (HTT) expression and developmental gene expression across the AHBA. Co-expression network analyses in 134 human brains free of neurodegenerative disorders were also performed. HTT expression was correlated with the expression of genes involved in neurodevelopment while co-expression network analyses also revealed that HTT expression was associated with developmental biological processes. Expression weighted cell-type enrichment (EWCE) analyses were used to explore which specific cell types were associated with HD cortical cell loss and these associations were validated using cell specific single nucleus RNAseq (snRNAseq) data from post-mortem HD brains. The developmental transcriptomic profile of cortical cell loss in preHD was enriched in astrocytes and endothelial cells, while the neurodegenerative transcriptomic profile was enriched for neuronal and microglial cells. Astrocyte-specific genes differentially expressed in HD post-mortem brains relative to controls using snRNAseq were enriched in the developmental transcriptomic profile, while neuronal and microglial-specific genes were enriched in the neurodegenerative transcriptomic profile. Our findings suggest that cortical cell loss in preHD may arise from dual pathological processes, emerging as a consequence of neurodevelopmental changes, at the beginning of life, followed by neurodegeneration in adulthood, targeting areas with reduced expression of synaptic and metabolic genes. These events result in age-related cell death across multiple brain cell types.

Connectome-based modelling of neurodegenerative diseases: towards precision medicine and mechanistic insight

Neurodegenerative diseases are the most common cause of dementia. Although their underlying molecular pathologies have been identified, there is substantial heterogeneity in the patterns of progressive brain alterations across and within these diseases. Recent advances in neuroimaging methods have revealed that pathological proteins accumulate along specific macroscale brain networks, implicating the network architecture of the brain in the system-level pathophysiology of neurodegenerative diseases. However, the extent to which 'network-based neurodegeneration' applies across the wide range of neurodegenerative disorders remains unclear. Here, we discuss the state-of-the-art of neuroimaging-based connectomics for the mapping and prediction of neurodegenerative processes. We review findings supporting brain networks as passive conduits through which pathological proteins spread. As an alternative view, we also discuss complementary work suggesting that network alterations actively modulate the spreading of pathological proteins between connected brain regions. We conclude this Perspective by proposing an integrative framework in which connectome-based models can be advanced along three dimensions of innovation: incorporating parameters that modulate propagation behaviour on the basis of measurable biological features; building patient-tailored models that use individual-level information and allowing model parameters to interact dynamically over time. We discuss promises and pitfalls of these strategies for improving disease insights and moving towards precision medicine.

New Challenges for Anatomists in the Era of Omics

Anatomic studies have traditionally relied on macroscopic, microscopic, and histological techniques to investigate the structure of tissues and organs. Anatomic studies are essential in many fields, including medicine, biology, and veterinary science. Advances in technology, such as imaging techniques and molecular biology, continue to provide new insights into the anatomy of living organisms. Therefore, anatomy remains an active and important area in the scientific field. The consolidation in recent years of some omics technologies such as genomics, transcriptomics, proteomics, and metabolomics allows for a more complete and detailed understanding of the structure and function of cells, tissues, and organs. These have been joined more recently by "omics" such as radiomics, pathomics, and connectomics, supported by computer-assisted technologies such as neural networks, 3D bioprinting, and artificial intelligence. All these new tools, although some are still in the early stages of development, have the potential to strongly contribute to the macroscopic and microscopic characterization in medicine. For anatomists, it is time to hitch a ride and get on board omics technologies to sail to new frontiers and to explore novel scenarios in anatomy.

Robust estimation of cortical similarity networks from brain MRI

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

Cortical profiles of numerous psychiatric disorders and normal development share a common pattern

The neurobiological bases of the association between development and psychopathology remain poorly understood. Here, we identify a shared spatial pattern of cortical thickness (CT) in normative development and several psychiatric and neurological disorders. Principal component analysis (PCA) was applied to CT of 68 regions in the Desikan-Killiany atlas derived from three large-scale datasets comprising a total of 41,075 neurotypical participants. PCA produced a spatially broad first principal component (PC1) that was reproducible across datasets. Then PC1 derived from healthy adult participants was compared to the pattern of CT differences associated with psychiatric and neurological disorders comprising a total of 14,886 cases and 20,962 controls from seven ENIGMA disease-related working groups, normative maturation and aging comprising a total of 17,697 scans from the ABCD Study® and the IMAGEN developmental study, and 17,075 participants from the ENIGMA Lifespan working group, as well as gene expression maps from the Allen Human Brain Atlas. Results revealed substantial spatial correspondences between PC1 and widespread lower CT observed in numerous psychiatric disorders. Moreover, the PC1 pattern was also correlated with the spatial pattern of normative maturation and aging. The transcriptional analysis identified a set of genes including KCNA2, KCNS1 and KCNS2 with expression patterns closely related to the spatial pattern of PC1. The gene category enrichment analysis indicated that the transcriptional correlations of PC1 were enriched to multiple gene ontology categories and were specifically over-represented starting at late childhood, coinciding with the onset of significant cortical maturation and emergence of psychopathology during the prepubertal-to-pubertal transition. Collectively, the present study reports a reproducible latent pattern of CT that captures interregional profiles of cortical changes in both normative brain maturation and a spectrum of psychiatric disorders. The pubertal timing of the expression of PC1-related genes implicates disrupted neurodevelopment in the pathogenesis of the spectrum of psychiatric diseases emerging during adolescence.

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.