A practical guide to linking brain-wide gene expression and neuroimaging data

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.

Gene co-expression changes underlying the functional connectomic alterations in Alzheimer's disease

BACKGROUND

There is growing evidence indicating that a number of functional connectivity networks are disrupted at each stage of the full clinical Alzheimer's disease spectrum. Such differences are also detectable in cognitive normal (CN) carrying mutations of AD risk genes, suggesting a substantial relationship between genetics and AD-altered functional brain networks. However, direct genetic effect on functional connectivity networks has not been measured.

METHODS

Leveraging existing AD functional connectivity studies collected in NeuroSynth, we performed a meta-analysis to identify two sets of brain regions: ones with altered functional connectivity in resting state network and ones without. Then with the brain-wide gene expression data in the Allen Human Brain Atlas, we applied a new biclustering method to identify a set of genes with differential co-expression patterns between these two set of brain regions.

RESULTS

Differential co-expression analysis using biclustering method led to a subset of 38 genes which showed distinctive co-expression patterns between AD-related and non AD-related brain regions in default mode network. More specifically, we observed 4 sub-clusters with noticeable co-expression difference, where the difference in correlations is above 0.5 on average.

CONCLUSIONS

This work applies a new biclustering method to search for a subset of genes with altered co-expression patterns in AD-related default mode network regions. Compared with traditional differential expression analysis, differential co-expression analysis yielded many more significant hits with extra insights into the wiring mechanism between genes. Particularly, the differential co-expression pattern was observed between two sets of genes, suggesting potential upstream genetic regulators in AD development.

Cellular and molecular signatures of in vivo imaging measures of GABAergic neurotransmission in the human brain

Diverse GABAergic interneuron networks orchestrate information processing in the brain. Understanding the principles underlying the organisation of this system in the human brain, and whether these principles are reflected by available non-invasive in vivo neuroimaging methods, is crucial for the study of GABAergic neurotransmission. Here, we use human gene expression data and state-of-the-art imaging transcriptomics to uncover co-expression patterns between genes encoding GABAA receptor subunits and inhibitory interneuron subtype-specific markers, and their association with binding patterns of the gold-standard GABA PET radiotracers [11C]Ro15-4513 and [11C]flumazenil. We found that the inhibitory interneuron marker somatostatin covaries with GABAA receptor-subunit genes GABRA5 and GABRA2, and that their distribution followed [11C]Ro15-4513 binding. In contrast, the inhibitory interneuron marker parvalbumin covaried with GABAA receptor-subunit genes GABRA1, GABRB2 and GABRG2, and their distribution tracked [11C]flumazenil binding. Our findings indicate that existing PET radiotracers may provide complementary information about key components of the GABAergic system.

Psychological Resilience and Neurodegenerative Risk: A Connectomics-Transcriptomics Investigation in Healthy Adolescent and Middle-Aged Females

Adverse life events can inflict substantial long-term damage, which, paradoxically, has been posited to stem from initially adaptative responses to the challenges encountered in one's environment. Thus, identification of the mechanisms linking resilience against recent stressors to longer-term psychological vulnerability is key to understanding optimal functioning across multiple timescales. To address this issue, our study tested the relevance of neuro-reproductive maturation and senescence, respectively, to both resilience and longer-term risk for pathologies characterised by accelerated brain aging, specifically, Alzheimer's Disease (AD). Graph theoretical and partial least squares analyses were conducted on multimodal imaging, reported biological aging and recent adverse experience data from the Lifespan Human Connectome Project (HCP). Availability of reproductive maturation/senescence measures restricted our investigation to adolescent (N =178) and middle-aged (N=146) females. Psychological resilience was linked to age-specific brain senescence patterns suggestive of precocious functional development of somatomotor and control-relevant networks (adolescence) and earlier aging of default mode and salience/ventral attention systems (middle adulthood). Biological aging showed complementary associations with the neural patterns relevant to resilience in adolescence (positive relationship) versus middle-age (negative relationship). Transcriptomic and expression quantitative trait locus data analyses linked the neural aging patterns correlated with psychological resilience in middle adulthood to gene expression patterns suggestive of increased AD risk. Our results imply a partially antagonistic relationship between resilience against proximal stressors and longer-term psychological adjustment in later life. They thus underscore the importance of fine-tuning extant views on successful coping by considering the multiple timescales across which age-specific processes may unfold.

Connectome gradient dysfunction in major depression and its association with gene expression profiles and treatment outcomes

Patients with major depressive disorder (MDD) exhibit concurrent deficits in both sensory and higher-order cognitive processing. Connectome studies have suggested a principal primary-to-transmodal gradient in functional brain networks, supporting the spectrum from sensation to cognition. However, whether this gradient structure is disrupted in patients with MDD and how this disruption associates with gene expression profiles and treatment outcome remain unknown. Using a large cohort of resting-state fMRI data from 2227 participants (1148 MDD patients and 1079 healthy controls) recruited at nine sites, we investigated MDD-related alterations in the principal connectome gradient. We further used Neurosynth, postmortem gene expression, and an 8-week antidepressant treatment (20 MDD patients) data to assess the meta-analytic cognitive functions, transcriptional profiles, and treatment outcomes related to MDD gradient alterations, respectively. Relative to the controls, MDD patients exhibited global topographic alterations in the principal primary-to-transmodal gradient, including reduced explanation ratio, gradient range, and gradient variation (Cohen's d = 0.16-0.21), and focal alterations mainly in the primary and transmodal systems (d = 0.18-0.25). These gradient alterations were significantly correlated with meta-analytic terms involving sensory processing and higher-order cognition. The transcriptional profiles explained 53.9% variance of the altered gradient pattern, with the most correlated genes enriched in transsynaptic signaling and calcium ion binding. The baseline gradient maps of patients significantly predicted symptomatic improvement after treatment. These results highlight the connectome gradient dysfunction in MDD and its linkage with gene expression profiles and clinical management, providing insight into the neurobiological underpinnings and potential biomarkers for treatment evaluation in this disorder.

Gene expression associated with human brain activations in facial expression recognition

Previous studies identified some genetic loci of emotion, but few focused on human emotion-related gene expression. In this study, the facial expression recognition (FER) task-based high-resolution fMRI data of 203 subjects in the Human Connectome Project (HCP) and expression data of the six healthy human postmortem brain tissues in the Allen Human Brain Atlas (AHBA) were used to conduct a transcriptome-neuroimaging spatial association analysis. Finally, 371 genes were identified to be significantly associated with FER-related brain activations. Enrichment analyses revealed that FER-related genes were mainly expressed in the brain, especially neurons, and might be related to cell junction organization, synaptic functions, and nervous system development regulation, indicating that FER was a complex polygenetic biological process involving multiple pathways. Moreover, these genes exhibited higher enrichment for psychiatric diseases with heavy emotion impairments. This study provided new insight into understanding the FER-related biological mechanisms and might be helpful to explore treatment methods for emotion-related psychiatric disorders.

Network structure and transcriptomic vulnerability shape atrophy in frontotemporal dementia

Connections among brain regions allow pathological perturbations to spread from a single source region to multiple regions. Patterns of neurodegeneration in multiple diseases, including behavioural variant of frontotemporal dementia (bvFTD), resemble the large-scale functional systems, but how bvFTD-related atrophy patterns relate to structural network organization remains unknown. Here we investigate whether neurodegeneration patterns in sporadic and genetic bvFTD are conditioned by connectome architecture. Regional atrophy patterns were estimated in both genetic bvFTD (75 patients, 247 controls) and sporadic bvFTD (70 patients, 123 controls). First, we identified distributed atrophy patterns in bvFTD, mainly targeting areas associated with the limbic intrinsic network and insular cytoarchitectonic class. Regional atrophy was significantly correlated with atrophy of structurally- and functionally-connected neighbours, demonstrating that network structure shapes atrophy patterns. The anterior insula was identified as the predominant group epicentre of brain atrophy using data-driven and simulation-based methods, with some secondary regions in frontal ventromedial and antero-medial temporal areas. We found that FTD-related genes, namely C9orf72 and TARDBP, confer local transcriptomic vulnerability to the disease, modulating the propagation of pathology through the connectome. Collectively, our results demonstrate that atrophy patterns in sporadic and genetic bvFTD are jointly shaped by global connectome architecture and local transcriptomic vulnerability, providing an explanation as to how heterogenous pathological entities can lead to the same clinical syndrome.

Integrative Functional, Molecular, and Transcriptomic Analyses of Altered Intrinsic Timescale Gradient in Depression

The pathophysiology and pharmacology of depression are hypothesized to be related to the imbalance of excitation–inhibition that gives rise to hierarchical dynamics (or intrinsic timescale gradient), further supporting a hierarchy of cortical functions. On this assumption, intrinsic timescale gradient is theoretically altered in depression. However, it remains unknown. We investigated altered intrinsic timescale gradient recently developed to measure hierarchical brain dynamics gradient and its underlying molecular architecture and brain-wide gene expression in depression. We first presented replicable intrinsic timescale gradient in two independent Chinese Han datasets and then investigated altered intrinsic timescale gradient and its possible underlying molecular and transcriptional bases in patients with depression. As a result, patients with depression showed stage-specifically shorter timescales compared with healthy controls according to illness duration. The shorter timescales were spatially correlated with monoamine receptor/transporter densities, suggesting the underlying molecular basis of timescale aberrance and providing clues to treatment. In addition, we identified that timescale aberrance-related genes ontologically enriched for synapse-related and neurotransmitter (receptor) terms, elaborating the underlying transcriptional basis of timescale aberrance. These findings revealed atypical timescale gradient in depression and built a link between neuroimaging, transcriptome, and neurotransmitter information, facilitating an integrative understanding of depression.

Dynamic network impairments underlie cognitive fluctuations in Lewy body dementia

Cognitive fluctuations are a characteristic and distressing disturbance of attention and consciousness seen in patients with Dementia with Lewy bodies and Parkinson's disease dementia. It has been proposed that fluctuations result from disruption of key neuromodulatory systems supporting states of attention and wakefulness which are normally characterised by temporally variable and highly integrated functional network architectures. In this study, patients with DLB (n = 25) and age-matched controls (n = 49) were assessed using dynamic resting state fMRI. A dynamic network signature of reduced temporal variability and integration was identified in DLB patients compared to controls. Reduced temporal variability correlated significantly with fluctuation-related measures using a sustained attention task. A less integrated (more segregated) functional network architecture was seen in DLB patients compared to the control group, with regions of reduced integration observed across dorsal and ventral attention, sensorimotor, visual, cingulo-opercular and cingulo-parietal networks. Reduced network integration correlated positively with subjective and objective measures of fluctuations. Regions of reduced integration and unstable regional assignments significantly matched areas of expression of specific classes of noradrenergic and cholinergic receptors across the cerebral cortex. Correlating topological measures with maps of neurotransmitter/neuromodulator receptor gene expression, we found that regions of reduced integration and unstable modular assignments correlated significantly with the pattern of expression of subclasses of noradrenergic and cholinergic receptors across the cerebral cortex. Altogether, these findings demonstrate that cognitive fluctuations are associated with an imaging signature of dynamic network impairment linked to specific neurotransmitters/neuromodulators within the ascending arousal system, highlighting novel potential diagnostic and therapeutic approaches for this troubling symptom.

"Less is more": a dose-response account of intranasal oxytocin pharmacodynamics in the human brain

Intranasal oxytocin is attracting attention as a potential treatment for several brain disorders due to promising preclinical results. However, translating findings to humans has been hampered by remaining uncertainties about its pharmacodynamics and the methods used to probe its effects in the human brain. Using a dose-response design (9, 18 and 36 IU), we demonstrate that intranasal oxytocin-induced changes in local regional cerebral blood flow (rCBF) in the amygdala at rest, and in the covariance between rCBF in the amygdala and other key hubs of the brain oxytocin system, follow a dose-response curve with maximal effects for lower doses. Yet, the effects on local rCBF might vary by amygdala subdivision, highlighting the need to qualify dose-response curves within subregion. We further link physiological changes with the density of the oxytocin receptor gene mRNA across brain regions, strengthening our confidence in intranasal oxytocin as a valid approach to engage central targets. Finally, we demonstrate that intranasal oxytocin does not disrupt cerebrovascular reactivity, which corroborates the validity of haemodynamic neuroimaging to probe the effects of intranasal oxytocin in the human brain. DATA AVAILABILITY: Participants did not consent for open sharing of the data. Therefore, data can only be accessed from the corresponding author upon reasonable request.

Multi-task learning based structured sparse canonical correlation analysis for brain imaging genetics

The advances in technologies for acquiring brain imaging and high-throughput genetic data allow the researcher to access a large amount of multi-modal data. Although the sparse canonical correlation analysis is a powerful bi-multivariate association analysis technique for feature selection, we are still facing major challenges in integrating multi-modal imaging genetic data and yielding biologically meaningful interpretation of imaging genetic findings. In this study, we propose a novel multi-task learning based structured sparse canonical correlation analysis (MTS2CCA) to deliver interpretable results and improve integration in imaging genetics studies. We perform comparative studies with state-of-the-art competing methods on both simulation and real imaging genetic data. On the simulation data, our proposed model has achieved the best performance in terms of canonical correlation coefficients, estimation accuracy, and feature selection accuracy. On the real imaging genetic data, our proposed model has revealed promising features of single-nucleotide polymorphisms and brain regions related to sleep. The identified features can be used to improve clinical score prediction using promising imaging genetic biomarkers. An interesting future direction is to apply our model to additional neurological or psychiatric cohorts such as patients with Alzheimer's or Parkinson's disease to demonstrate the generalizability of our method.

Dynamic Functional Connectivity Alterations and Their Associated Gene Expression Pattern in Autism Spectrum Disorders

Autism spectrum disorders (ASDs) are a group of heterogeneous neurodevelopmental disorders that are highly heritable and are associated with impaired dynamic functional connectivity (DFC). However, the molecular mechanisms behind DFC alterations remain largely unknown. Eighty-eight patients with ASDs and 87 demographically matched typical controls (TCs) from the Autism Brain Imaging Data Exchange II database were included in this study. A seed-based sliding window approach was then performed to investigate the DFC changes in each of the 29 seeds in 10 classic resting-state functional networks and the whole brain. Subsequently, the relationships between DFC alterations in patients with ASDs and their symptom severity were assessed. Finally, transcription-neuroimaging association analyses were conducted to explore the molecular mechanisms of DFC disruptions in patients with ASDs. Compared with TCs, patients with ASDs showed significantly increased DFC between the right dorsolateral prefrontal cortex (DLPFC) and left fusiform/lingual gyrus, between the DLPFC and the superior temporal gyrus, between the right frontal eye field (FEF) and left middle frontal gyrus, between the FEF and the right angular gyrus, and between the left intraparietal sulcus and the right middle temporal gyrus. Moreover, significant relationships between DFC alterations and symptom severity were observed. Furthermore, the genes associated with DFC changes in ASDs were identified by performing gene-wise across-sample spatial correlation analysis between gene expression extracted from six donors’ brain of the Allen Human Brain Atlas and case-control DFC difference. In enrichment analysis, these genes were enriched for processes associated with synaptic signaling and voltage-gated ion channels and calcium pathways; also, these genes were highly expressed in autistic disorder, chronic alcoholic intoxication and several disorders related to depression. These results not only demonstrated higher DFC in patients with ASDs but also provided novel insight into the molecular mechanisms underlying these alterations.

Transcriptomic Signatures Associated With Gray Matter Volume Changes in Patients With Functional Constipation

Functional constipation, which belongs to the functional gastrointestinal disorder (FGID), is a common disease and significantly impacts daily life. FGID patients have been progressively proven with functional and structural alterations in various brain regions, but whether and how functional constipation affects the brain gray matter volume (GMV) remains unclear; besides, which genes are associated with the GMV changes in functional constipation is largely unknown. On account of the structural MRI image from the 30 functional constipation patients and 30 healthy controls (HCs), GMV analysis showed that functional constipation patients had significantly decreased GMV in the right orbital prefrontal cortex (OFC), left precentral gyrus (PreG), and bilateral thalamus (THA). Correlation analysis showed that the self-rating depressive scale, patient assessment of constipation quality of life (PAC-QOL), and Wexner constipation scores were negatively correlated with GMV of the OFC and negative correlations between PAC-QOL score and GMV of the bilateral THA. Based on the Allen Human Brain Atlas, a cross-sample spatial correlation was conducted and found that 18 genes’ expression values showed robust correlations with GMV changes in functional constipation patients. These outcomes highlight our recognition of the transcriptional features related to GMV changes in functional constipation and could be regarded as candidates to detect biological mechanisms of abnormality in functional constipation patients.

Structural brain and spinal cord damage in symptomatic and pre-symptomatic VAPB-related ALS

INTRODUCTION

The VAPB gene is associated with fALS (fALS 8). This disease presents a variable phenotype and no study sought to characterize its neuroanatomical abnormalities until now. This study aims to evaluate structural brain and spinal cord abnormalities in symptomatic and pre-symptomatic VAPB-related ALS.

METHODS

This cohort included 10 presymptomatic and 20 symptomatic carriers of the Pro56Ser VAPB variant as well as 30 matched controls and 20 individuals with sporadic ALS. They underwent detailed clinical evaluation and MRI in a 3 T scanner. Using volumetric T1 sequence, we computed cerebral cortical thickness (FreeSurfer), basal ganglia volumetry (T1 Multi-atlas) and SC morphometry (SpineSeg). DTI was used to assess white matter integrity (DTI Multi-atlas). Groups were compared using a generalized linear model with Bonferroni-corrected p values<0.05. We also plotted VAPB brain expression map using Allen Human Brain Atlas to compare with imaging findings.

RESULTS

Mean age of presymptomatic and symptomatic subjects were 43.2 and 51.9 years, respectively. Most patients had a predominant lower motor neuron phenotype (16/20). Sleep complaints and tremor were the most frequent additional manifestations. Compared to controls, symptomatic subjects had pallidal, brainstem and SC atrophy, whereas presymptomatic only had SC atrophy. This pattern also contrasted with the sALS group that presented motor cortex and corticospinal abnormalities. Brain structural damage and VAPB expression maps were highly overlapping.

CONCLUSION

VAPB-related ALS has a distinctive structural signature that targets the basal ganglia, brainstem and SC, which are regions with high VAPB expression. Neuroanatomical SC changes are evident before clinical onset of the disease.

Imaging transcriptomics: Convergent cellular, transcriptomic, and molecular neuroimaging signatures in the healthy adult human brain

The integration of transcriptomic and neuroimaging data, "imaging transcriptomics," has recently emerged to generate hypotheses about potential biological pathways underlying regional variability in neuroimaging features. However, the validity of this approach is yet to be examined in depth. Here, we sought to bridge this gap by performing transcriptomic decoding of the regional distribution of well-known molecular markers spanning different elements of the biology of the healthy human brain. Imaging transcriptomics identifies biological and cell pathways that are consistent with the known biology of a wide range of molecular neuroimaging markers. The extent to which it can capture patterns of gene expression that align well with elements of the biology of the neuroinflammatory axis, at least in healthy controls without a proinflammatory challenge, is inconclusive. Imaging transcriptomics might constitute an interesting approach to improve our understanding of the biological pathways underlying regional variability in a wide range of neuroimaging phenotypes.

Dissect Relationships Between Gene Co-expression and Functional Connectivity in Human Brain

Although recent evidence indicates an association between gene co-expression and functional connectivity in human brain, specific association patterns remain largely unknown. Here, using neuroimaging-based functional connectivity data of living brains and brain-wide gene expression data of postmortem brains, we performed comprehensive analyses to dissect relationships between gene co-expression and functional connectivity. We identified 125 connectivity-related genes (20 novel genes) enriched for dendrite extension, signaling pathway and schizophrenia, and 179 gene-related functional connections mainly connecting intra-network regions, especially homologous cortical regions. In addition, 51 genes were associated with connectivity in all brain functional networks and enriched for action potential and schizophrenia; in contrast, 51 genes showed network-specific modulatory effects and enriched for ion transportation. These results indicate that functional connectivity is unequally affected by gene expression, and connectivity-related genes with different biological functions are involved in connectivity modulation of different networks.

Choroid plexus enlargement is associated with neuroinflammation and reduction of blood brain barrier permeability in depression

BACKGROUND

Recent studies have shown that choroid plexuses (CP) may be involved in the neuro-immune axes, playing a role in the interaction between the central and peripheral inflammation. Here we aimed to investigate CP volume alterations in depression and their associations with inflammation.

METHODS

51 depressed participants (HDRS score > 13) and 25 age- and sex-matched healthy controls (HCs) from the Wellcome Trust NIMA consortium were re-analysed for the study. All the participants underwent full peripheral cytokine profiling and simultaneous [11C]PK11195 PET/structural MRI imaging for measuring neuroinflammation and CP volume respectively.

RESULTS

We found a significantly greater CP volume in depressed subjects compared to HCs (t(76) = +2.17) that was positively correlated with [11C]PK11195 PET binding in the anterior cingulate cortex (r = 0.28, p = 0.02), prefrontal cortex (r = 0.24, p = 0.04), and insular cortex (r = 0.24, p = 0.04), but not with the peripheral inflammatory markers: CRP levels (r = 0.07, p = 0.53), IL-6 (r = -0.08, p = 0.61), and TNF-α (r = -0.06, p = 0.70). The CP volume correlated with the [11C]PK11195 PET binding in CP (r = 0.34, p = 0.005). Integration of transcriptomic data from the Allen Human Brain Atlas with the brain map depicting the correlations between CP volume and PET imaging found significant gene enrichment for several pathways involved in neuroinflammatory response.

CONCLUSION

This result supports the hypothesis that changes in brain barriers may cause reduction in solute exchanges between blood and CSF, disturbing the brain homeostasis and ultimately contributing to inflammation in depression. Given that CP anomalies have been recently detected in other brain disorders, these results may not be specific to depression and might extend to other conditions with a peripheral inflammatory component.

Coexpression of Gene Transcripts with Monoamine Oxidase A Quantified by Human In Vivo Positron Emission Tomography

The monoamine oxidase A (MAO-A) is integral to monoamine metabolism and is thus relevant to the pathophysiology of various neuropsychiatric disorders; however, associated gene-enzyme relations are not well understood. This study aimed to unveil genes coexpressed with MAO-A. Therefore, 18 179 mRNA expression maps (based on the Allen Human Brain Atlas) were correlated with the cerebral distribution volume (VT) of MAO-A assessed in 36 healthy subjects (mean age ± standard deviation: 32.9 ± 8.8 years, 18 female) using [11C]harmine positron emission tomography scans. Coexpression analysis was based on Spearman's ρ, over-representation tests on Fisher's exact test with false discovery rate (FDR) correction. The analysis revealed 35 genes in cortex (including B-cell translocation gene family, member 3, implicated in neuroinflammation) and 247 genes in subcortex (including kallikrein-related peptidase 10, implicated in Alzheimer's disease). Significantly over-represented Gene Ontology terms included "neuron development", "neuron differentiation", and "cell-cell signaling" as well as "axon" and "neuron projection". In vivo MAO-A enzyme distribution and MAOA expression did not correlate in cortical areas (ρ = 0.08) while correlation was found in subcortical areas (ρ = 0.52), suggesting influences of region-specific post-transcriptional and -translational modifications. The herein reported information could contribute to guide future genetic studies, deepen the understanding of associated pathomechanisms and assist in the pursuit of novel therapeutic targets.

Transcriptomic and macroscopic architectures of intersubject functional variability in human brain white-matter

Intersubject variability is a fundamental characteristic of brain organizations, and not just "noise". Although intrinsic functional connectivity (FC) is unique to each individual and varies across brain gray-matter, the underlying mechanisms of intersubject functional variability in white-matter (WM) remain unknown. This study identified WMFC variabilities and determined the genetic basis and macroscale imaging in 45 healthy subjects. The functional localization pattern of intersubject variability across WM is heterogeneous, with most variability observed in the heteromodal cortex. The variabilities of heteromodal regions in expression profiles of genes are related to neuronal cells, involved in synapse-related and glutamic pathways, and associated with psychiatric disorders. In contrast, genes overexpressed in unimodal regions are mostly expressed in glial cells and were related to neurological diseases. Macroscopic variability recapitulates the functional and structural specializations and behavioral phenotypes. Together, our results provide clues to intersubject variabilities of the WMFC with convergent transcriptomic and cellular signatures, which relate to macroscale brain specialization.

Statistical testing in transcriptomic-neuroimaging studies: A how-to and evaluation of methods assessing spatial and gene specificity

Multiscale integration of gene transcriptomic and neuroimaging data is becoming a widely used approach for exploring the molecular underpinnings of large‐scale brain organization in health and disease. Proper statistical evaluation of determined associations between imaging‐based phenotypic and transcriptomic data is key in these explorations, in particular to establish whether observed associations exceed “chance level” of random, nonspecific effects. Recent approaches have shown the importance of statistical models that can correct for spatial autocorrelation effects in the data to avoid inflation of reported statistics. Here, we discuss the need for examination of a second category of statistical models in transcriptomic‐neuroimaging analyses, namely those that can provide “gene specificity.” By means of a couple of simple examples of commonly performed transcriptomic‐neuroimaging analyses, we illustrate some of the potentials and challenges of transcriptomic‐imaging analyses, showing that providing gene specificity on observed transcriptomic‐neuroimaging effects is of high importance to avoid reports of nonspecific effects. Through means of simulations we show that the rate of reported nonspecific effects (i.e., effects that cannot be specifically linked to a specific gene or gene‐set) can run as high as 60%, with only less than 5% of transcriptomic‐neuroimaging associations observed through ordinary linear regression analyses showing both spatial and gene specificity. We provide a discussion, a tutorial, and an easy‐to‐use toolbox for the different options of null models in transcriptomic‐neuroimaging analyses.

Where the genome meets the connectome: Understanding how genes shape human brain connectivity

The integration of modern neuroimaging methods with genetically informative designs and data can shed light on the molecular mechanisms underlying the structural and functional organization of the human connectome. Here, we review studies that have investigated the genetic basis of human brain network structure and function through three complementary frameworks: (1) the quantification of phenotypic heritability through classical twin designs; (2) the identification of specific DNA variants linked to phenotypic variation through association and related studies; and (3) the analysis of correlations between spatial variations in imaging phenotypes and gene expression profiles through the integration of neuroimaging and transcriptional atlas data. We consider the basic foundations, strengths, limitations, and discoveries associated with each approach. We present converging evidence to indicate that anatomical connectivity is under stronger genetic influence than functional connectivity and that genetic influences are not uniformly distributed throughout the brain, with phenotypic variation in certain regions and connections being under stronger genetic control than others. We also consider how the combination of imaging and genetics can be used to understand the ways in which genes may drive brain dysfunction in different clinical disorders.

Transcriptional signatures of synaptic vesicle genes define myotonic dystrophy type I neurodegeneration

AIM

To delineate the neurogenetic profiles of brain degeneration patterns in myotonic dystrophy type I (DM1).

METHODS

In two cohorts of DM1 patients, brain maps of volume loss (VL) and neuropsychological deficits (NDs) were intersected to large-scale transcriptome maps provided by the Allen Human Brain Atlas (AHBA). For validation, neuropathological and RNA analyses were performed in a small series of DM1 brain samples.

RESULTS

Twofold: (1) From a list of preselected hypothesis-driven genes, confirmatory analyses found that three genes play a major role in brain degeneration: dystrophin (DMD), alpha-synuclein (SNCA) and the microtubule-associated protein tau (MAPT). Neuropathological analyses confirmed a highly heterogeneous Tau-pathology in DM1, different to the one in Alzheimer's disease. (2) Exploratory analyses revealed gene clusters enriched for key biological processes in the central nervous system, such as synaptic vesicle recycling, localization, endocytosis and exocytosis, and the serotonin and dopamine neurotransmitter pathways. RNA analyses confirmed synaptic vesicle dysfunction.

CONCLUSIONS

The combination of large-scale transcriptome interactions with brain imaging and cognitive function sheds light on the neurobiological mechanisms of brain degeneration in DM1 that might help define future therapeutic strategies and research into this condition.

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

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

Brain atrophy progression in Parkinson's disease is shaped by connectivity and local vulnerability

Brain atrophy has been reported in the early stages of Parkinson's disease, but there have been few longitudinal studies. How intrinsic properties of the brain, such as anatomical connectivity, local cell-type distribution and gene expression combine to determine the pattern of disease progression also remains unknown. One hypothesis proposes that the disease stems from prion-like propagation of misfolded alpha-synuclein via the connectome that might cause varying degrees of tissue damage based on local properties. Here, we used MRI data from the Parkinson Progression Markers Initiative to map the progression of brain atrophy over 1, 2 and 4 years compared with baseline. We derived atrophy maps for four time points using deformation-based morphometry applied to T1-weighted MRI from 120 de novo Parkinson's disease patients, 74 of whom had imaging at all four time points (50 Men: 24 Women) and 157 healthy control participants (115 Men: 42 Women). In order to determine factors that may influence neurodegeneration, we related atrophy progression to brain structural and functional connectivity, cell-type expression and gene ontology enrichment analyses. After regressing out the expected age and sex effects associated with normal ageing, we found that atrophy significantly progressed over 2 and 4 years in the caudate, nucleus accumbens, hippocampus and posterior cortical regions. This progression was shaped by both structural and functional brain connectivity. Also, the progression of atrophy was more pronounced in regions with a higher expression of genes related to synapses and was inversely related to the prevalence of oligodendrocytes and endothelial cells. In sum, we demonstrate that the progression of atrophy in Parkinson's disease is in line with the prion-like propagation hypothesis of alpha-synuclein and provide evidence that synapses may be especially vulnerable to synucleinopathy. In addition to identifying vulnerable brain regions, this study reveals different factors that may be implicated in the neurotoxic mechanisms leading to progression in Parkinson's disease. All brain maps generated here are available on request.

Standardizing workflows in imaging transcriptomics with the abagen toolbox

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

A dopamine gradient controls access to distributed working memory in the large-scale monkey cortex

Dopamine is required for working memory, but how it modulates the large-scale cortex is unknown. Here, we report that dopamine receptor density per neuron, measured by autoradiography, displays a macroscopic gradient along the macaque cortical hierarchy. This gradient is incorporated in a connectome-based large-scale cortex model endowed with multiple neuron types. The model captures an inverted U-shaped dependence of working memory on dopamine and spatial patterns of persistent activity observed in over 90 experimental studies. Moreover, we show that dopamine is crucial for filtering out irrelevant stimuli by enhancing inhibition from dendrite-targeting interneurons. Our model revealed that an activity-silent memory trace can be realized by facilitation of inter-areal connections and that adjusting cortical dopamine induces a switch from this internal memory state to distributed persistent activity. Our work represents a cross-level understanding from molecules and cell types to recurrent circuit dynamics underlying a core cognitive function distributed across the primate cortex.

Transcriptomics of the depressed and PTSD brain

Stress is the response of an organism to demands for change, yet excessive or chronic stress contributes to nearly all psychiatric disorders. The advent of high-throughput transcriptomic methods such as single cell RNA sequencing poses new opportunities to understand the neurobiology of stress, yet substantial barriers to understanding stress remain. Stress adaptation is an organismal survival mechanism conserved across all organisms, yet there is an infinity of potential stressful experiences. Unraveling shared and separate transcriptional programs for adapting to stressful experience remains a challenge, despite methodological and analytic advances. Here we review the state of the field focusing on the technologies used to study the transcriptome for the stress neurobiologist, and also attempt to identify central questions about the heterogeneity of stress for those applying transcriptomic approaches. We further explore how postmortem transcriptome studies aided by preclinical animal models are converging on common molecular pathways for adaptation to aversive experience. Finally, we discuss approaches to integrate large genomic datasets with human neuroimaging and other datasets.

Cross-Sectional Characterization of Local Brain Network Connectivity Pre and Post Breast Cancer Treatment and Distinct Association With Subjective Cognitive and Psychological Function

Objective: We aimed to characterize local brain network connectivity in long-term breast cancer survivors compared to newly diagnosed patients. Methods: Functional magnetic resonance imaging (fMRI) and subjective cognitive and psychological function data were obtained from a group of 76 newly diagnosed, pre-treatment female patients with breast cancer (mean age 57 ± 7 years) and a separate group of 80, post-treatment, female breast cancer survivors (mean age 58 ± 8; mean time since treatment 44 ± 43 months). The network-based statistic (NBS) was used to compare connectivity of local brain edges between groups. Hubs were defined as nodes with connectivity indices one standard deviation or more above network mean and were further classified as provincial (higher intra-subnetwork connectivity) or connector (higher inter-subnetwork connectivity) using the participation coefficient. We determined the hub status of nodes encompassing significantly different edges and correlated the centralities of edges with behavioral measures. Results: The post-treatment group demonstrated significantly lower subjective cognitive function (W = 3,856, p = 0.004) but there were no group differences in psychological distress (W = 2,866, p = 0.627). NBS indicated significantly altered connectivity (p < 0.042, corrected) in the post-treatment group compared to the pre-treatment group largely in temporal, frontal-temporal and temporal-parietal areas. The majority of the regions projecting these connections (78%) met criteria for hub status and significantly less of these hubs were connectors in the post-treatment group (z = 1.85, p = 0.031). Subjective cognitive function and psychological distress were correlated with largely non-overlapping edges in the post-treatment group (p < 0.05). Conclusion: Widespread functional network alterations are evident in long-term survivors of breast cancer compared to newly diagnosed patients. We also demonstrated that there are both overlapping and unique brain network signatures for subjective cognitive function vs. psychological distress.

Relating quantitative 7T MRI across cortical depths to cytoarchitectonics, gene expression and connectomics

Ultra-high field MRI across the depth of the cortex has the potential to provide anatomically precise biomarkers and mechanistic insights into neurodegenerative disease like Huntington's disease that show layer-selective vulnerability. Here we compare multi-parametric mapping (MPM) measures across cortical depths for a 7T 500 μm whole brain acquisition to (a) layer-specific cell measures from the von Economo histology atlas, (b) layer-specific gene expression, using the Allen Human Brain atlas and (c) white matter connections using high-fidelity diffusion tractography, at a 1.3 mm isotropic voxel resolution, from a 300mT/m Connectom MRI system. We show that R2*, but not R1, across cortical depths is highly correlated with layer-specific cell number and layer-specific gene expression. R1- and R2*-weighted connectivity strength of cortico-striatal and intra-hemispheric cortical white matter connections was highly correlated with grey matter R1 and R2* across cortical depths. Limitations of the layer-specific relationships demonstrated are at least in part related to the high cross-correlations of von Economo atlas cell counts and layer-specific gene expression across cortical layers. These findings demonstrate the potential and limitations of combining 7T MPMs, gene expression and white matter connections to provide an anatomically precise framework for tracking neurodegenerative disease.

The relationship between transcription and eccentricity in human V1

Gene expression gradients radiating from regions of primary sensory cortices have recently been described and are thought to underlie the large-scale organization of the human cerebral cortex. However, the role of transcription in the functional layout of a single region within the adult brain has yet to be clarified, likely owing to the difficulty of identifying a brain region anatomically consistent enough across individuals with dense enough tissue sampling. Overcoming these hurdles in human primary visual cortex (V1), we show a relationship between differential gene expression and the cortical layout of eccentricity in human V1. Interestingly, these genes are unique from those previously identified that contribute to the positioning of cortical areas in the visual processing hierarchy. Enrichment analyses show that a subset of the identified genes encode for structures related to inhibitory interneurons, ion channels, as well as cellular projections, and are expressed more in foveal compared to peripheral portions of human V1. These findings predict that tissue density should be higher in foveal compared to peripheral V1. Using a histological pipeline, we validate this prediction using Nissl-stained sections of postmortem occipital cortex. We discuss these findings relative to previous studies in non-human primates, as well as in the context of an organizational pattern in which the adult human brain employs transcription gradients at multiple spatial scales: across the cerebral cortex, across areas within processing hierarchies, and within single cortical areas.

Genetic Architecture Underlying Differential Resting-state Functional Connectivity of Subregions Within the Human Visual Cortex

The human visual cortex is a heterogeneous entity that has multiple subregions showing substantial variability in their functions and connections. We aimed to identify genes associated with resting-state functional connectivity (rsFC) of visual subregions using transcriptome-neuroimaging spatial correlations in discovery and validation datasets. Results showed that rsFC of eight visual subregions were associated with expression measures of eight gene sets, which were specifically expressed in brain tissue and showed the strongest correlations with visual behavioral processes. Moreover, there was a significant divergence in these gene sets and their functional features between medial and lateral visual subregions. Relative to those associated with lateral subregions, more genes associated with medial subregions were found to be enriched for neuropsychiatric diseases and more diverse biological functions and pathways, and to be specifically expressed in multiple types of neurons and immune cells and during the middle and late stages of cortical development. In addition to shared behavioral processes, lateral subregion associated genes were uniquely correlated with high-order cognition. These findings of commonalities and differences in the identified rsFC-related genes and their functional features across visual subregions may improve our understanding of the functional heterogeneity of the visual cortex from the perspective of underlying genetic architecture.

Transcriptomic Signatures Associated With Regional Cortical Thickness Changes in Parkinson's Disease

Cortical atrophy is a common manifestation in Parkinson's disease (PD), particularly in advanced stages of the disease. To elucidate the molecular underpinnings of cortical thickness changes in PD, we performed an integrated analysis of brain-wide healthy transcriptomic data from the Allen Human Brain Atlas and patterns of cortical thickness based on T1-weighted anatomical MRI data of 149 PD patients and 369 controls. For this purpose, we used partial least squares regression to identify gene expression patterns correlated with cortical thickness changes. In addition, we identified gene expression patterns underlying the relationship between cortical thickness and clinical domains of PD. Our results show that genes whose expression in the healthy brain is associated with cortical thickness changes in PD are enriched in biological pathways related to sumoylation, regulation of mitotic cell cycle, mitochondrial translation, DNA damage responses, and ER-Golgi traffic. The associated pathways were highly related to each other and all belong to cellular maintenance mechanisms. The expression of genes within most pathways was negatively correlated with cortical thickness changes, showing higher expression in regions associated with decreased cortical thickness (atrophy). On the other hand, sumoylation pathways were positively correlated with cortical thickness changes, showing higher expression in regions with increased cortical thickness (hypertrophy). Our findings suggest that alterations in the balanced interplay of these mechanisms play a role in changes of cortical thickness in PD and possibly influence motor and cognitive functions.

Imaging Transcriptomics of Brain Disorders

Noninvasive neuroimaging is a powerful tool for quantifying diverse aspects of brain structure and function in vivo, and it has been used extensively to map the neural changes associated with various brain disorders. However, most neuroimaging techniques offer only indirect measures of underlying pathological mechanisms. The recent development of anatomically comprehensive gene expression atlases has opened new opportunities for studying the transcriptional correlates of noninvasively measured neural phenotypes, offering a rich framework for evaluating pathophysiological hypotheses and putative mechanisms. Here, we provide an overview of some fundamental methods in imaging transcriptomics and outline their application to understanding brain disorders of neurodevelopment, adulthood, and neurodegeneration. Converging evidence indicates that spatial variations in gene expression are linked to normative changes in brain structure during age-related maturation and neurodegeneration that are in part associated with cell-specific gene expression markers of gene expression. Transcriptional correlates of disorder-related neuroimaging phenotypes are also linked to transcriptionally dysregulated genes identified in ex vivo analyses of patient brains. Modeling studies demonstrate that spatial patterns of gene expression are involved in regional vulnerability to neurodegeneration and the spread of disease across the brain. This growing body of work supports the utility of transcriptional atlases in testing hypotheses about the molecular mechanism driving disease-related changes in macroscopic neuroimaging phenotypes.

Connectivity alterations in autism reflect functional idiosyncrasy

Autism spectrum disorder (ASD) is commonly understood as an alteration of brain networks, yet case-control analyses against typically-developing controls (TD) have yielded inconsistent results. Here, we devised a novel approach to profile the inter-individual variability in functional network organization and tested whether such idiosyncrasy contributes to connectivity alterations in ASD. Studying a multi-centric dataset with 157 ASD and 172 TD, we obtained robust evidence for increased idiosyncrasy in ASD relative to TD in default mode, somatomotor and attention networks, but also reduced idiosyncrasy in lateral temporal cortices. Idiosyncrasy increased with age and significantly correlated with symptom severity in ASD. Furthermore, while patterns of functional idiosyncrasy were not correlated with ASD-related cortical thickness alterations, they co-localized with the expression patterns of ASD risk genes. Notably, we could demonstrate that patterns of atypical idiosyncrasy in ASD closely overlapped with connectivity alterations that are measurable with conventional case-control designs and may, thus, be a principal driver of inconsistency in the autism connectomics literature. These findings support important interactions between inter-individual heterogeneity in autism and functional signatures. Our findings provide novel biomarkers to study atypical brain development and may consolidate prior research findings on the variable nature of connectome level anomalies in autism.

Brain areas affected by intranasal oxytocin show higher oxytocin receptor expression

Neuroimaging studies suggest that intranasal oxytocin (IN‐OXT) may modulate emotional and social processes by altering neural activity patterns. The extent of brain penetration after IN‐OXT is unclear, and it is currently unknown whether IN‐OXT can directly bind central oxytocin receptors (OXTRs). We investigated oxytocin pathway gene expression in regions affected by IN‐OXT on task‐based fMRI. We found that OXTR is more highly expressed in affected than unaffected subcortical regions; this effect did not vary by task type or sex. Cortical results revealed higher OXTR expression in regions affected by IN‐OXT in emotional processing tasks and in male‐only data. No significant differences were found in expression of the closely related vasopressin receptors. Our findings suggest that the mechanism by which IN‐OXT may alter brain functionality involves direct activation of central OXTRs.

Convergent structural network and gene signatures for MRgFUS thalamotomy in patients with Parkinson's disease

MRgFUS has just been made available for the 1.7 million Parkinson's disease patients in China. Despite its non-invasive and rapid therapeutic advantages for involuntary tremor, some concerns have emerged about outcomes variability, non-specificity, and side-effects, as little is known about its impact on the long-term plasticity of brain structure. We sought to dissect the characteristics of long-term changes in brain structure caused by MRgFUS lesion and explored potential biological mechanisms. One-year multimodal imaging follow-ups were conducted for nine tremor-dominant Parkinson's disease patients undergoing unilateral MRgFUS thalamotomy. A structural connectivity map was generated for each patient to analyze dynamic changes in brain structure. The human brain transcriptome was extracted and spatially registered for connectivity vulnerability. Genetic functional enrichment analysis was performed and further clarified using in vivo emission computed tomography data. MRgFUS not only abolished tremors but also significantly disrupted the brain network topology. Network-based statistics identified a U-shape MRgFUS-sensitive subnetwork reflective of hand tremor recovery and surgical process, accompanied by relevant cerebral blood flow and gray matter alteration. Using human brain gene expression data, we observed that dopaminergic signatures were responsible for the preferential vulnerability associated with these architectural alterations. Additional PET/SPECT data not only validated these gene signatures, but also suggested that structural alteration was significantly correlated with D1 and D2 receptors, DAT, and F-DOPA measures. There was a long-term dynamic loop between structural alteration and dopaminergic signature for MRgFUS thalamotomy, which may be closely related to the long-term improvements in clinical tremor.

Mapping gene transcription and neurocognition across human neocortex

Regulation of gene expression drives protein interactions that govern synaptic wiring and neuronal activity. The resulting coordinated activity among neuronal populations supports complex psychological processes, yet how gene expression shapes cognition and emotion remains unknown. Here, we directly bridge the microscale and macroscale by mapping gene expression patterns to functional activation patterns across the cortical sheet. Applying unsupervised learning to the Allen Human Brain Atlas and Neurosynth databases, we identify a ventromedial-dorsolateral gradient of gene assemblies that separate affective and perceptual domains. This topographic molecular-psychological signature reflects the hierarchical organization of the neocortex, including systematic variations in cell type, myeloarchitecture, laminar differentiation and intrinsic network affiliation. In addition, this molecular-psychological signature strengthens over neurodevelopment and can be replicated in two independent repositories. Collectively, our results reveal spatially covarying transcriptomic and cognitive architectures, highlighting the influence that molecular mechanisms exert on psychological processes.

Regional radiomics similarity networks (R2SNs) in the human brain: Reproducibility, small-world properties and a biological basis

A structural covariance network (SCN) has been used successfully in structural magnetic resonance imaging (sMRI) studies. However, most SCNs have been constructed by a unitary marker that is insensitive for discriminating different disease phases. The aim of this study was to devise a novel regional radiomics similarity network (R2SN) that could provide more comprehensive information in morphological network analysis. R2SNs were constructed by computing the Pearson correlations between the radiomics features extracted from any pair of regions for each subject (AAL atlas). We further assessed the small-world property of R2SNs, and we evaluated the reproducibility in different datasets and through test-retest analysis. The relationships between the R2SNs and general intelligence/interregional coexpression of genes were also explored. R2SNs could be replicated in different datasets, regardless of the use of different feature subsets. R2SNs showed high reproducibility in the test-retest analysis (intraclass correlation coefficient > 0.7). In addition, the small-word property (σ > 2) and the high correlation between gene expression (R = 0.29, p < 0.001) and general intelligence were determined for R2SNs. Furthermore, the results have also been repeated in the Brainnetome atlas. R2SNs provide a novel, reliable, and biologically plausible method to understand human morphological covariance based on sMRI.

Genes associated with grey matter volume reduction in multiple sclerosis

There is extensive grey matter volume (GMV) reduction in multiple sclerosis (MS), which may account for cognitive impairment in this disabling disorder. Although genome-wide association studies (GWASs) have identified hundreds of genes associated with MS, we know little about which genes associated with GMV reduction and cognitive decline in MS. In the present study, we aimed to uncover genes associated with GMV reduction in MS by performing cross-sample (1473 brain tissue samples) partial least squares regression between gene expression from 6 postmortem brains and case-control GMV difference of MS from a meta-analysis of 1391 patients and 1189 controls (discovery phase) and from the intergroup comparison between 69 patients and 70 controls (replication phase). We identified 623 genes whose brain spatial expression profiles were significantly associated with GMV reduction in MS. These genes showed significant enrichment for MS-related genes identified by GWAS; were functionally associated with ion channel, synaptic transmission, axon and neuron projection; and showed more significant cell type-specific expression in neurons than other cell types. More importantly, the identified genes showed significant enrichment for those genes with downregulated rather than upregulated expression in MS. The spatial distribution patterns of the expression of the identified genes showed more significant correlations with brain activation patterns of memory and language tasks. These findings indicate that grey matter atrophy in MS may be resulted from the joint effects of multiple genes that are associated with this disorder, especially genes with downregulated expression in MS.

Cortical patterning of morphometric similarity gradient reveals diverged hierarchical organization in sensory-motor cortices

The topological organization of the cerebral cortex provides hierarchical axes, namely gradients, which reveal systematic variations of brain structure and function. However, the hierarchical organization of macroscopic brain morphology and how it constrains cortical function along the organizing axes remains unclear. We map the gradient of cortical morphometric similarity (MS) connectome, combining multiple features conceptualized as a "fingerprint" of an individual's brain. The principal MS gradient is anchored by motor and sensory cortices at two extreme ends, which are reliable and reproducible. Notably, divergences between motor and sensory hierarchies are consistent with the laminar histological thickness gradient but contrary to the canonical functional connectivity (FC) gradient. Moreover, the MS dissociates with FC gradients in the higher-order association cortices. The MS gradient recapitulates fundamental properties of cortical organization, from gene expression and cyto- and myeloarchitecture to evolutionary expansion. Collectively, our findings provide a heuristic hierarchical organization of cortical morphological neuromarkers.

A Transcriptome Community-and-Module Approach of the Human Mesoconnectome

Graph analysis allows exploring transcriptome compartments such as communities and modules for brain mesostructures. In this work, we proposed a bottom-up model of a gene regulatory network to brain-wise connectome workflow. We estimated the gene communities across all brain regions from the Allen Brain Atlas transcriptome database. We selected the communities method to yield the highest number of functional mesostructures in the network hierarchy organization, which allowed us to identify specific brain cell functions (e.g., neuroplasticity, axonogenesis and dendritogenesis communities). With these communities, we built brain-wise region modules that represent the connectome. Our findings match with previously described anatomical and functional brain circuits, such the default mode network and the default visual network, supporting the notion that the brain dynamics that carry out low- and higher-order functions originate from the modular composition of a GRN complex network.

Progressive Stabilization of Brain Network Dynamics during Childhood and Adolescence

Functional brain networks require dynamic reconfiguration to support flexible cognitive function. However, the developmental principles shaping brain network dynamics remain poorly understood. Here, we report the longitudinal development of large-scale brain network dynamics during childhood and adolescence, and its connection with gene expression profiles. Using a multilayer network model, we show the temporally varying modular architecture of child brain networks, with higher network switching primarily in the association cortex and lower switching in the primary regions. This topographical profile exhibits progressive maturation, which manifests as reduced modular dynamics, particularly in the transmodal (e.g., default-mode and frontoparietal) and sensorimotor regions. These developmental refinements mediate age-related enhancements of global network segregation and are linked with the expression profiles of genes associated with the enrichment of ion transport and nucleobase-containing compound transport. These results highlight a progressive stabilization of brain dynamics, which expand our understanding of the neural mechanisms that underlie cognitive development.

Unveiling the neuroimaging-genetic intersections in the human brain

PURPOSE OF REVIEW

The prevalence of new public datasets of brain-wide and single-cell transcriptome data has created new opportunities to link neuroimaging findings with genetic data. The aim of this study is to present the different methodological approaches that have been used to combine this data.

RECENT FINDINGS

Drawing from various sources of open access data, several studies have been able to correlate neuroimaging maps with spatial distribution of brain expression. These efforts have enabled researchers to identify functional annotations of related genes, identify specific cell types related to brain phenotypes, study the expression of genes across life span and highlight the importance of selected brain genes in disease genetic networks.

SUMMARY

New transcriptome datasets and methodological approaches complement current neuroimaging work and will be crucial to improve our understanding of the biological mechanism that underlies many neurological conditions.

Early-life trauma endophenotypes and brain circuit-gene expression relationships in functional neurological (conversion) disorder

Functional neurological (conversion) disorder (FND) is a neuropsychiatric condition whereby individuals present with sensorimotor symptoms incompatible with other neurological disorders. Early-life maltreatment (ELM) is a risk factor for developing FND, yet few studies have investigated brain network-trauma relationships in this population. In this neuroimaging-gene expression study, we used two graph theory approaches to elucidate ELM subtype effects on resting-state functional connectivity architecture in 30 patients with motor FND. Twenty-one individuals with comparable depression, anxiety, and ELM scores were used as psychiatric controls. Thereafter, we compared trauma endophenotypes in FND with regional differences in transcriptional gene expression as measured by the Allen Human Brain Atlas (AHBA). In FND patients only, we found that early-life physical abuse severity, and to a lesser extent physical neglect, correlated with corticolimbic weighted-degree functional connectivity. Connectivity profiles influenced by physical abuse occurred in limbic (amygdalar-hippocampal), paralimbic (cingulo-insular and ventromedial prefrontal), and cognitive control (ventrolateral prefrontal) areas, as well as in sensorimotor and visual cortices. These findings held adjusting for individual differences in depression/anxiety, PTSD, and motor phenotypes. In FND, physical abuse also correlated with amygdala and insula coupling to motor cortices. In exploratory analyses, physical abuse correlated connectivity maps overlapped with the AHBA spatial expression of three gene clusters: (i) neuronal morphogenesis and synaptic transmission genes in limbic/paralimbic areas; (ii) locomotory behavior and neuronal generation genes in left-lateralized structures; and (iii) nervous system development and cell motility genes in right-lateralized structures. These circuit-specific architectural profiles related to individual differences in childhood physical abuse burden advance our understanding of the pathophysiology of FND.

Regional brain iron and gene expression provide insights into neurodegeneration in Parkinson's disease

The mechanisms responsible for the selective vulnerability of specific neuronal populations in Parkinson's disease are poorly understood. Oxidative stress secondary to brain iron accumulation is one postulated mechanism. We measured iron deposition in 180 cortical regions of 96 patients with Parkinson's disease and 35 control subjects using quantitative susceptibility mapping. We estimated the expression of 15 745 genes in the same regions using transcriptomic data from the Allen Human Brain Atlas. Using partial least squares regression, we then identified the profile of gene transcription in the healthy brain that underlies increased cortical iron in patients with Parkinson's disease relative to controls. Applying gene ontological tools, we investigated the biological processes and cell types associated with this transcriptomic profile and identified the sets of genes with spatial expression profiles in control brains that correlated significantly with the spatial pattern of cortical iron deposition in Parkinson's disease. Gene ontological analyses revealed that these genes were enriched for biological processes relating to heavy metal detoxification, synaptic function and nervous system development and were predominantly expressed in astrocytes and glutamatergic neurons. Furthermore, we demonstrated that the genes differentially expressed in Parkinson's disease are associated with the pattern of cortical expression identified in this study. Our findings provide mechanistic insights into regional selective vulnerabilities in Parkinson's disease, particularly the processes involving iron accumulation.

Structural Features of the Human Connectome That Facilitate the Switching of Brain Dynamics via Noradrenergic Neuromodulation

The structural connectivity of human brain allows the coexistence of segregated and integrated states of activity. Neuromodulatory systems facilitate the transition between these functional states and recent computational studies have shown how an interplay between the noradrenergic and cholinergic systems define these transitions. However, there is still much to be known about the interaction between the structural connectivity and the effect of neuromodulation, and to what extent the connectome facilitates dynamic transitions. In this work, we use a whole brain model, based on the Jasen and Rit equations plus a human structural connectivity matrix, to find out which structural features of the human connectome network define the optimal neuromodulatory effects. We simulated the effect of the noradrenergic system as changes in filter gain, and studied its effects related to the global-, local-, and meso-scale features of the connectome. At the global-scale, we found that the ability of the network of transiting through a variety of dynamical states is disrupted by randomization of the connection weights. By simulating neuromodulation of partial subsets of nodes, we found that transitions between integrated and segregated states are more easily achieved when targeting nodes with greater connection strengths-local feature-or belonging to the rich club-meso-scale feature. Overall, our findings clarify how the network spatial features, at different levels, interact with neuromodulation to facilitate the switching between segregated and integrated brain states and to sustain a richer brain dynamics.

Genetic influences on hub connectivity of the human connectome

Brain network hubs are both highly connected and highly inter-connected, forming a critical communication backbone for coherent neural dynamics. The mechanisms driving this organization are poorly understood. Using diffusion-weighted magnetic resonance imaging in twins, we identify a major role for genes, showing that they preferentially influence connectivity strength between network hubs of the human connectome. Using transcriptomic atlas data, we show that connected hubs demonstrate tight coupling of transcriptional activity related to metabolic and cytoarchitectonic similarity. Finally, comparing over thirteen generative models of network growth, we show that purely stochastic processes cannot explain the precise wiring patterns of hubs, and that model performance can be improved by incorporating genetic constraints. Our findings indicate that genes play a strong and preferential role in shaping the functionally valuable, metabolically costly connections between connectome hubs.

A generative network model of neurodevelopmental diversity in structural brain organization

The formation of large-scale brain networks, and their continual refinement, represent crucial developmental processes that can drive individual differences in cognition and which are associated with multiple neurodevelopmental conditions. But how does this organization arise, and what mechanisms drive diversity in organization? We use generative network modeling to provide a computational framework for understanding neurodevelopmental diversity. Within this framework macroscopic brain organization, complete with spatial embedding of its organization, is an emergent property of a generative wiring equation that optimizes its connectivity by renegotiating its biological costs and topological values continuously over time. The rules that govern these iterative wiring properties are controlled by a set of tightly framed parameters, with subtle differences in these parameters steering network growth towards different neurodiverse outcomes. Regional expression of genes associated with the simulations converge on biological processes and cellular components predominantly involved in synaptic signaling, neuronal projection, catabolic intracellular processes and protein transport. Together, this provides a unifying computational framework for conceptualizing the mechanisms and diversity in neurodevelopment, capable of integrating different levels of analysis-from genes to cognition.