Spatial analysis and high resolution mapping of the human whole-brain transcriptome for integrative analysis in neuroimaging

Challenges and Frontiers in Computational Metabolic Psychiatry

One of the primary challenges in metabolic psychiatry is that the disrupted brain functions that underlie psychiatric conditions arise from a complex set of downstream and feedback processes spanning across multiple spatiotemporal scales. Importantly, the same circuit can have multiple points of failure, each of which results in a different type of dysregulation, and thus elicits distinct cascades downstream that produce divergent signs and symptoms. Here, we illustrate this challenge by examining how subtle differences in circuit perturbations can lead to divergent clinical outcomes. We also discuss how computational models can perform the spatially heterogenous integration and bridge in vitro and in vivo paradigms. By leveraging recent methodological advances and tools, computational models can integrate relevant processes across scales (e.g., TCA-cycle, ion channel, neural microassembly, whole-brain macro-circuit) and across physiological systems (e.g., neural, endocrine, immune, vascular), providing a framework that can unite these mechanistic processes in a manner that goes beyond the conceptual and descriptive, to the quantitative and generative. These hold the potential to sharpen our intuitions towards circuit-based models for personalized diagnostics and treatment.

Synaptic signaling modeled by functional connectivity predicts metabolic demands of the human brain

PURPOSE

The human brain is characterized by interacting large-scale functional networks fueled by glucose metabolism. Since former studies could not sufficiently clarify how these functional connections shape glucose metabolism, we aimed to provide a neurophysiologically-based approach.

METHODS

51 healthy volunteers underwent simultaneous PET/MRI to obtain BOLD functional connectivity and [18F]FDG glucose metabolism. These multimodal imaging proxies of fMRI and PET were combined in a whole-brain extension of metabolic connectivity mapping. Specifically, functional connectivity of all brain regions were used as input to explain glucose metabolism of a given target region. This enabled the modeling of postsynaptic energy demands by incoming signals from distinct brain regions.

RESULTS

Functional connectivity input explained a substantial part of metabolic demands but with pronounced regional variations (34 - 76%). During cognitive task performance this multimodal association revealed a shift to higher network integration compared to resting state. In healthy aging, a dedifferentiation (decreased segregated/modular structure of the brain) of brain networks during rest was observed. Furthermore, by including data from mRNA maps, [11C]UCB-J synaptic density and aerobic glycolysis (oxygen-to-glucose index from PET data), we show that whole-brain functional input reflects non-oxidative, on-demand metabolism of synaptic signaling. The metabolically-derived directionality of functional inputs further marked them as top-down predictions. In addition, the approach uncovered formerly hidden networks with superior efficiency through metabolically informed network partitioning.

CONCLUSIONS

Applying multimodal imaging, we decipher a crucial part of the metabolic and neurophysiological basis of functional connections in the brain as interregional on-demand synaptic signaling fueled by anaerobic metabolism. The observed task- and age-related effects indicate promising future applications to characterize human brain function and clinical alterations.

In-vivo neuronal dysfunction by Aβ and tau overlaps with brain-wide inflammatory mechanisms in Alzheimer's disease

The molecular mechanisms underlying neuronal dysfunction in Alzheimer's disease (AD) remain uncharacterized. Here, we identify genes, molecular pathways and cellular components associated with whole-brain dysregulation caused by amyloid-beta (Aβ) and tau deposits in the living human brain. We obtained in-vivo resting-state functional MRI (rs-fMRI), Aβ- and tau-PET for 47 cognitively unimpaired and 16 AD participants from the Translational Biomarkers in Aging and Dementia cohort. Adverse neuronal activity impacts by Aβ and tau were quantified with personalized dynamical models by fitting pathology-mediated computational signals to the participant's real rs-fMRIs. Then, we detected robust brain-wide associations between the spatial profiles of Aβ-tau impacts and gene expression in the neurotypical transcriptome (Allen Human Brain Atlas). Within the obtained distinctive signature of in-vivo neuronal dysfunction, several genes have prominent roles in microglial activation and in interactions with Aβ and tau. Moreover, cellular vulnerability estimations revealed strong association of microglial expression patterns with Aβ and tau's synergistic impact on neuronal activity (q < 0.001). These results further support the central role of the immune system and neuroinflammatory pathways in AD pathogenesis. Neuronal dysregulation by AD pathologies also associated with neurotypical synaptic and developmental processes. In addition, we identified drug candidates from the vast LINCS library to halt or reduce the observed Aβ-tau effects on neuronal activity. Top-ranked pharmacological interventions target inflammatory, cancer and cardiovascular pathways, including specific medications undergoing clinical evaluation in AD. Our findings, based on the examination of molecular-pathological-functional interactions in humans, may accelerate the process of bringing effective therapies into clinical practice.

Functionality of arousal-regulating brain circuitry at rest predicts human cognitive abilities

Arousal state is regulated by subcortical neuromodulatory nuclei, such as locus coeruleus, which send wide-reaching projections to cortex. Whether higher-order cortical regions have the capacity to recruit neuromodulatory systems to aid cognition is unclear. Here, we hypothesized that select cortical regions activate the arousal system, which, in turn, modulates large-scale brain activity, creating a functional circuit predicting cognitive ability. We utilized the Human Connectome Project 7T functional magnetic resonance imaging dataset (n = 149), acquired at rest with simultaneous eye tracking, along with extensive cognitive assessment for each subject. First, we discovered select frontoparietal cortical regions that drive large-scale spontaneous brain activity specifically via engaging the arousal system. Second, we show that the functionality of the arousal circuit driven by bilateral posterior cingulate cortex (associated with the default mode network) predicts subjects' cognitive abilities. This suggests that a cortical region that is typically associated with self-referential processing supports cognition by regulating the arousal system.

Transcriptomic contributions to a modern cytoarchitectonic parcellation of the human cerebral cortex

Transcriptomic contributions to the anatomical, functional, and network layout of the human cerebral cortex (HCC) have become a major interest in cognitive and systems neuroscience. Here, we tested if transcriptomic differences support a modern, algorithmic cytoarchitectonic parcellation of HCC. Using a data-driven approach, we identified a sparse subset of genes that differentially contributed to the cytoarchitectonic parcellation of HCC. A combined metric of cortical thickness and myelination (CT/M ratio), as well as cell density, correlated with gene expression. Enrichment analyses showed that genes specific to the cytoarchitectonic parcellation of the HCC were related to molecular functions such as transmembrane transport and ion channel activity. Together, the relationship between transcriptomics and cytoarchitecture bridges the gap among (i) gradients at the macro-scale (including thickness and myelination), (ii) areas at the meso-scale, and (iii) cell density at the microscale, as well as supports the recently proposed cortical spectrum theory and structural model.

InSpectro-Gadget: A Tool for Estimating Neurotransmitter and Neuromodulator Receptor Distributions for MRS Voxels

Magnetic resonance spectroscopy (MRS) is widely used to estimate concentrations of glutamate and γ -aminobutyric acid (GABA) in specific regions of the living human brain. As cytoarchitectural properties differ across the brain, interpreting these measurements can be assisted by having knowledge of such properties for the MRS region(s) studied. In particular, some knowledge of likely local neurotransmitter receptor patterns can potentially give insights into the mechanistic environment GABA- and glutamatergic neurons are functioning in. This may be of particular utility when comparing two or more regions, given that the receptor populations may differ substantially across them. At the same time, when studying MRS data from multiple participants or timepoints, the homogeneity of the sample becomes relevant, as measurements taken from areas with different cytoarchitecture may be difficult to compare. To provide insights into the likely cytoarchitectural environment of user-defined regions-of-interest, we produced an easy to use tool - InSpectro-Gadget - that interfaces with receptor mRNA expression information from the Allen Human Brain Atlas. This Python tool allows users to input masks and automatically obtain a graphical overview of the receptor population likely to be found within. This includes comparison between multiple masks or participants where relevant. The receptors and receptor subunit genes featured include GABA- and glutamatergic classes, along with a wide range of neuromodulators. The functionality of the tool is explained here and its use is demonstrated through a set of example analyses. The tool is available at https://github.com/lizmcmanus/Inspectro-Gadget .

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.

Transcriptional cartography integrates multiscale biology of the human cortex

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

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

OBJECTIVE

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

METHODS

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

RESULTS

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

INTERPRETATION

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

Propofol Disrupts the Functional Core-Matrix Architecture of the Thalamus in Humans

Research into the role of thalamocortical circuits in anesthesia-induced unconsciousness is difficult due to anatomical and functional complexity. Prior neuroimaging studies have examined either the thalamus as a whole or focused on specific subregions, overlooking the distinct neuronal subtypes like core and matrix cells. We conducted a study of heathy volunteers and functional magnetic resonance imaging during conscious baseline, deep sedation, and recovery. We advanced the functional gradient mapping technique to delineate the functional geometry of thalamocortical circuits, within a framework of the unimodal-transmodal functional axis of the cortex. We observed a significant shift in this geometry during unconsciousness, marked by the dominance of unimodal over transmodal geometry. This alteration was closely linked to the spatial variations in the density of matrix cells within the thalamus. This research bridges cellular and systems-level understanding, highlighting the crucial role of thalamic core-matrix functional architecture in understanding the neural mechanisms of states of consciousness.

Association of protein distribution and gene expression revealed by positron emission tomography and postmortem gene expression in the dopaminergic system of the human brain

PURPOSE

The topological distribution of dopamine-related proteins is determined by gene transcription and subsequent regulations. Recent research strategies integrating positron emission tomography with a transcriptome atlas have opened new opportunities to understand the influence of regulation after transcription on protein distribution. Previous studies have reported that messenger (m)-RNA expression levels spatially correlate with the density maps of serotonin receptors but not with those of transporters. This discrepancy may be due to differences in regulation after transcription between presynaptic and postsynaptic proteins, which have not been studied in the dopaminergic system. Here, we focused on dopamine D1 and D2/D3 receptors and dopamine transporters and investigated their region-wise relationship between mRNA expression and protein distribution.

METHODS

We examined the region-wise correlation between regional binding potentials of the target region relative to that of non-displaceable tissue (BPND) values of 11C-SCH-23390 and mRNA expression levels of dopamine D1 receptors (D1R); regional BPND values of 11C-FLB-457 and mRNA expression levels of dopamine D2/D3 receptors (D2/D3R); and regional total distribution volume (VT) values of 18F-FE-PE2I and mRNA expression levels of dopamine transporters (DAT) using Spearman's rank correlation.

RESULTS

We found significant positive correlations between regional BPND values of 11C-SCH-23390 and the mRNA expression levels of D1R (r = 0.769, p = 0.0021). Similar to D1R, regional BPND values of 11C-FLB-457 positively correlated with the mRNA expression levels of D2R (r = 0.809, p = 0.0151) but not with those of D3R (r = 0.413, p = 0.3095). In contrast to D1R and D2R, no significant correlation between VT values of 18F-FE-PE2I and mRNA expression levels of DAT was observed (r = -0.5934, p = 0.140).

CONCLUSION

We found a region-wise correlation between the mRNA expression levels of dopamine D1 and D2 receptors and their respective protein distributions. However, we found no region-wise correlation between the mRNA expression levels of dopamine transporters and their protein distributions, indicating different regulatory mechanisms for the localization of pre- and postsynaptic proteins. These results provide a broader understanding of the application of the transcriptome atlas to neuroimaging studies of the dopaminergic nervous system.

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

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

From neurotransmitters to networks: Transcending organisational hierarchies with molecular-informed functional imaging

The human brain exhibits complex interactions across micro, meso-, and macro-scale organisational principles. Recent synergistic multi-modal approaches have begun to link micro-scale information to systems level dynamics, transcending organisational hierarchies and offering novel perspectives into the brain's function and dysfunction. Specifically, the distribution of micro-scale properties (such as receptor density or gene expression) can be mapped onto macro-scale measures from functional MRI to provide novel neurobiological insights. Methodological approaches to enrich functional imaging analyses with molecular information are rapidly evolving, with several streams of research having developed relatively independently, each offering unique potential to explore the trans-hierarchical functioning of the brain. Here, we address the three principal streams of research - spatial correlation, molecular-enriched network, and in-silico whole brain modelling analyses - to provide a critical overview of the different sources of molecular information, how this information can be utilised within analyses of fMRI data, the merits and pitfalls of each methodology, and, through the use of key examples, highlight their promise to shed new light on key domains of neuroscientific inquiry.

APOEε4 associates with microglial activation independently of Aβ plaques and tau tangles

Animal studies suggest that the apolipoprotein E ε4 (APOEε4) allele is a culprit of early microglial activation in Alzheimer's disease (AD). Here, we tested the association between APOEε4 status and microglial activation in living individuals across the aging and AD spectrum. We studied 118 individuals with positron emission tomography for amyloid-β (Aβ; [18F]AZD4694), tau ([18F]MK6240), and microglial activation ([11C]PBR28). We found that APOEε4 carriers presented increased microglial activation relative to noncarriers in early Braak stage regions within the medial temporal cortex accounting for Aβ and tau deposition. Furthermore, microglial activation mediated the Aβ-independent effects of APOEε4 on tau accumulation, which was further associated with neurodegeneration and clinical impairment. The physiological distribution of APOE mRNA expression predicted the patterns of APOEε4-related microglial activation in our population, suggesting that APOE gene expression may regulate the local vulnerability to neuroinflammation. Our results support that the APOEε4 genotype exerts Aβ-independent effects on AD pathogenesis by activating microglia in brain regions associated with early tau deposition.

Distribution of cholinergic nerve terminals in the aged human brain measured with [18F]FEOBV PET and its correlation with histological data

INTRODUCTION

[¹⁸F]fluoroetoxybenzovesamicol ([¹⁸F]FEOBV) is a positron emission topography (PET) tracer for the vesicular acetylcholine transporter (VAChT), a protein located predominantly in synaptic vesicles in cholinergic nerve terminals. We aimed to use [¹⁸F]FEOBV PET to study the cholinergic topography of the healthy human brain.

MATERIALS AND METHODS

[¹⁸F]FEOBV PET brain data volumes of healthy elderly humans were normalized to standard space and intensity-normalized to the white matter. Stereotactic atlases of regions of interest were superimposed to describe and quantify tracer distribution. The spatial distribution of [¹⁸F]FEOBV PET uptake was compared with histological and gene expression data.

RESULTS

Twenty participants of both sexes and a mean age of 73.9 ± 6.0 years, age-range [64; 86], were recruited. Highest tracer binding was present in the striatum, some thalamic nuclei, and the basal forebrain. Intermediate binding was found in most nuclei of the brainstem, thalamus, and hypothalamus; the vermis and flocculonodular lobe; and the hippocampus, amygdala, insula, cingulate, olfactory cortex, and Heschl's gyrus. Lowest binding was present in most areas of the cerebral cortex, and in the cerebellar nuclei and hemispheres. The spatial distribution of tracer correlated with immunohistochemical post-mortem data, as well as with regional expression levels of SLC18A3, the VAChT coding gene.

DISCUSSION

Our in vivo findings confirm the regional cholinergic distribution in specific brain structures as described post-mortem. A positive spatial correlation between tracer distribution and regional gene expression levels further corroborates [¹⁸F]FEOBV PET as a validated tool for in vivo cholinergic imaging. The study represents an advancement in the continued efforts to delineate the spatial topography of the human cholinergic system in vivo.

Grey matter volume loss in Parkinson's disease psychosis and its relationship with serotonergic gene expression: A meta-analysis

BACKGROUND

Neuroanatomical alterations underlying psychosis in Parkinson's Disease (PDP) remain unclear. We carried out a meta-analysis of MRI studies investigating the neural correlates of PDP and examined its relation with dopaminergic and serotonergic receptor gene expression.

METHODS

PubMed, Web of Science and Embase were searched for MRI studies (k studies = 10) of PDP compared to PD patients without psychosis (PDnP). Seed-based d Mapping with Permutation of Subject Images and multiple linear regression analyses was used to examine the relationship between pooled estimates of grey matter volume (GMV) loss in PDP and D1/D2 and 5-HT1a/5-HT2a receptor gene expression estimates from Allen Human Brain Atlas.

RESULTS

We observed lower grey matter volume in parietal-temporo-occipital regions (PDP n = 211, PDnP, n = 298). GMV loss in PDP was associated with local expression of 5-HT1a (b = 0.109, p = 0.012) and 5-HT2a receptors (b= -0.106, p = 0.002) but not dopaminergic receptors.

CONCLUSION

Widespread GMV loss in the parieto-temporo-occipital regions may underlie PDP. Association between grey matter volume and local expression of serotonergic receptor genes may suggest a role for serotonergic receptors in PDP.

Brain-wide and cell-specific transcriptomic insights into MRI-derived cortical morphology in macaque monkeys

Integrative analyses of transcriptomic and neuroimaging data have generated a wealth of information about biological pathways underlying regional variability in imaging-derived brain phenotypes in humans, but rarely in nonhuman primates due to the lack of a comprehensive anatomically-defined atlas of brain transcriptomics. Here we generate complementary bulk RNA-sequencing dataset of 819 samples from 110 brain regions and single-nucleus RNA-sequencing dataset, and neuroimaging data from 162 cynomolgus macaques, to examine the link between brain-wide gene expression and regional variation in morphometry. We not only observe global/regional expression profiles of macaque brain comparable to human but unravel a dorsolateral-ventromedial gradient of gene assemblies within the primate frontal lobe. Furthermore, we identify a set of 971 protein-coding and 34 non-coding genes consistently associated with cortical thickness, specially enriched for neurons and oligodendrocytes. These data provide a unique resource to investigate nonhuman primate models of human diseases and probe cross-species evolutionary mechanisms.

Association of CNR1 gene and cannabinoid 1 receptor protein in the human brain

We aimed to integrate genomic mapping from brain mRNA atlas with the protein expression from positron emission tomography (PET) scans of type 1 cannabinoid (CB1) receptor and to compare the predictive power of CB1 receptor with those of other neuroreceptor/transporters using a meta-analysis. Volume of distribution (V_T ) from F18-FMPEP-d2 PET scans, CNR1 gene (Cannabinoid receptor 1) expression, and H3-CP55940 binding were calculated and correlation analysis was performed. Between V_T of F18-FMPEP-d2 PET scans and CNR1 mRNA expression, moderate strength of correlation was observed (rho = .5067, p = .0337). Strong positive correlation was also found between CNR1 mRNA expression and H3-CP55940 binding (r = .6336, p = .0364), validating the finding between F18-FMPEP-d2 PET scans and CNR1 mRNA. The correlation between V_T of F18-FMPEP-d2 PET scans and H3-CP55940 binding was marginally significant (r = .5025, p = .0563). From the meta-analysis, the correlation coefficient between mRNA expression and protein expressions ranged from -.10 to .99, with a pooled effect of .76. In conclusion, we observed the moderate to strong associations between gene and protein expression for CB1 receptor in the human brain, which was validated by autoradiography. We combined the autoradiographic finding with PET of CB1 receptor, producing the density atlas map of CB1 receptor. From the meta-analysis, the moderate to strong correlation was observed between mRNA expression and protein expressions across multiple genes. Further study is needed to investigate the relationship between multiple genes and in vivo proteins to improve and accelerate drug development.

Regional gene expression patterns are associated with task-specific brain activation during reward and emotion processing measured with functional MRI

The exploration of the spatial relationship between gene expression profiles and task-evoked response patterns known to be altered in neuropsychiatric disorders, for example depression, can guide the development of more targeted therapies. Here, we estimated the correlation between human transcriptome data and two different brain activation maps measured with functional magnetic resonance imaging (fMRI) in healthy subjects. Whole-brain activation patterns evoked during an emotional face recognition task were associated with topological mRNA expression of genes involved in cellular transport. In contrast, fMRI activation patterns related to the acceptance of monetary rewards were associated with genes implicated in cellular localization processes, metabolism, translation, and synapse regulation. An overlap of these genes with risk genes from major depressive disorder genome-wide association studies revealed the involvement of the master regulators TCF4 and PAX6 in emotion and reward processing. Overall, the identification of stable relationships between spatial gene expression profiles and fMRI data may reshape the prospects for imaging transcriptomics studies.

APOE ε4 genotype, amyloid-β, and sex interact to predict tau in regions of high APOE mRNA expression

The apolipoprotein E (APOE) ε4 allele is strongly linked with cerebral β-amyloidosis, but its relationship with tauopathy is less established. We investigated the relationship between APOE ε4 carrier status, regional amyloid-β (Aβ), magnetic resonance imaging (MRI) volumetrics, tau positron emission tomography (PET), APOE messenger RNA (mRNA) expression maps, and cerebrospinal fluid phosphorylated tau (CSF ptau₁₈₁). Three hundred fifty participants underwent imaging, and 270 had ptau₁₈₁. We used computational models to evaluate the main effect of APOE ε4 carrier status on regional neuroimaging values and then the interaction of ε4 status and global Aβ on regional tau PET and brain volumes as well as CSF ptau₁₈₁. Separately, we also examined the additional interactive influence of sex. We found that, for the same degree of Aβ burden, APOE ε4 carriers showed greater tau PET signal relative to noncarriers in temporal regions, but no interaction was present for MRI volumes or CSF ptau₁₈₁. This potentiation of tau aggregation irrespective of sex occurred in brain regions with high APOE mRNA expression, suggesting local vulnerabilities to tauopathy. There were greater effects of APOE genotype in females, although the interactive sex effects did not strongly mirror mRNA expression. Pathology is not homogeneously expressed throughout the brain but mirrors underlying biological patterns such as gene expression.

Task-independent acute effects of delta-9-tetrahydrocannabinol on human brain function and its relationship with cannabinoid receptor gene expression: A neuroimaging meta-regression analysis

The neurobiological mechanisms underlying the effects of delta-9-tetrahydrocannabinol (THC) remain unclear. Here, we examined the spatial acute effect of THC on human regional brain activation or blood flow (hereafter called 'activation signal') in a 'core' network of brain regions from 372 participants, tested using a within-subject repeated measures design under experimental conditions. We also investigated whether the neuromodulatory effects of THC are related to the local expression of the cannabinoid-type-1 (CB1R) and type-2 (CB2R) receptors. Finally, we investigated the dose-response relationship between THC and key brain substrates. These meta-analytic findings shed new light on the localisation of the effects of THC in the human brain, suggesting that THC has neuromodulatory effects in regions central to many cognitive tasks and processes, related to dose, with greater effects in regions with higher levels of CB1R expression.

Correlation of receptor density and mRNA expression patterns in the human cerebral cortex

Changes in distribution of associated molecular targets have been reported across several neuropsychiatric disorders. However, the high-resolution topology of most proteins is unknown and simultaneous in vivo measurement in multi-receptor systems is complicated. To account for the missing proteomic information, messenger ribonucleic acid (mRNA) transcripts are typically used as a surrogate. Nonetheless, post-transcriptional and post-translational processes might cause the discrepancy between the final distribution of proteins and gene expression patterns. Therefore, this study aims to investigate ex vivo links between mRNA expression and corresponding receptor density in the human cerebral cortex. To this end, autoradiography data on the density of 15 different receptors in 38 brain regions were correlated with the expression patterns of 50 associated genes derived from microarray data (mA), RNA sequencing data (RNA-Seq) provided by the Allen Human Brain Atlas and predicted mRNA expression patterns (pred-mRNA). Spearman's rank correlation was used to evaluate the possible links between proteomic data and mRNA expression patterns. Correlations between mRNA and protein density varied greatly between targets: Positive associations were found for e.g. the serotonin 1A (pred-mRNA: rs = 0.708; mA: rs = 0.601) or kainate receptor (pred-mRNA: rs = 0.655; mA: rs = 0.601; RNA-Seq: rs = 0.575) as well as a few negative associations e.g. γ-Aminobutyric acid (GABA) A receptor subunit α3 (pred-mRNA: rs = -0.638; mA: rs = -0.619) or subunit α5 (pred-mRNA: rs = -0.565; mA: rs = -0.563), while most of the other investigated target receptors showed low correlations. The high variability in the correspondence of mRNA expression and receptor spatial distribution warrants caution when inferring the topology of molecular targets in the brain from transcriptome data. This not only highlights the longstanding value of molecular imaging but also indicates a need for comprehensive proteomic studies.

[11C]Martinostat PET analysis reveals reduced HDAC I availability in Alzheimer's disease

Alzheimer’s disease (AD) is characterized by the brain accumulation of amyloid-β and tau proteins. A growing body of literature suggests that epigenetic dysregulations play a role in the interplay of hallmark proteinopathies with neurodegeneration and cognitive impairment. Here, we aim to characterize an epigenetic dysregulation associated with the brain deposition of amyloid-β and tau proteins. Using positron emission tomography (PET) tracers selective for amyloid-β, tau, and class I histone deacetylase (HDAC I isoforms 1–3), we find that HDAC I levels are reduced in patients with AD. HDAC I PET reduction is associated with elevated amyloid-β PET and tau PET concentrations. Notably, HDAC I reduction mediates the deleterious effects of amyloid-β and tau on brain atrophy and cognitive impairment. HDAC I PET reduction is associated with 2-year longitudinal neurodegeneration and cognitive decline. We also find HDAC I reduction in the postmortem brain tissue of patients with AD and in a transgenic rat model expressing human amyloid-β plus tau pathology in the same brain regions identified in vivo using PET. These observations highlight HDAC I reduction as an element associated with AD pathophysiology.

The link between amyloid and tau proteins with Alzheimer’s disease progression remains unclear. Here, the authors propose HDACs I downregulation as an element linking the deleterious effects of brain proteinopathies with disease progression.

Effects of APOE4 allelic dosage on lipidomic signatures in the entorhinal cortex of aged mice

Apolipoprotein E ε4 (APOE4) is the primary genetic risk factor for the late-onset form of Alzheimer's disease (AD). Although the reason for this association is not completely understood, researchers have uncovered numerous effects of APOE4 expression on AD-relevant brain processes, including amyloid beta (Aβ) accumulation, lipid metabolism, endosomal-lysosomal trafficking, and bioenergetics. In this study, we aimed to determine the effect of APOE4 allelic dosage on regional brain lipid composition in aged mice, as well as in cultured neurons. We performed a targeted lipidomic analysis on an AD-vulnerable brain region (entorhinal cortex; EC) and an AD-resistant brain region (primary visual cortex; PVC) from 14-15 month-old APOE3/3, APOE3/4, and APOE4/4 targeted replacement mice, as well as on neurons cultured with conditioned media from APOE3/3 or APOE4/4 astrocytes. Our results reveal that the EC possesses increased susceptibility to APOE4-associated lipid alterations compared to the PVC. In the EC, APOE4 expression showed a dominant effect in decreasing diacylglycerol (DAG) levels, and a semi-dominant, additive effect in the upregulation of multiple ceramide, glycosylated sphingolipid, and bis(monoacylglycerol)phosphate (BMP) species, lipids known to accumulate as a result of endosomal-lysosomal dysfunction. Neurons treated with conditioned media from APOE4/4 vs. APOE3/3 astrocytes showed similar alterations of DAG and BMP species to those observed in the mouse EC. Our results suggest that APOE4 expression differentially modulates regional neuronal lipid signatures, which may underlie the increased susceptibility of EC-localized neurons to AD pathology.

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.

A developmental reduction of the excitation:inhibition ratio in association cortex during adolescence

Adolescence is hypothesized to be a critical period for the development of association cortex. A reduction of the excitation:inhibition (E:I) ratio is a hallmark of critical period development; however, it has been unclear how to assess the development of the E:I ratio using noninvasive neuroimaging techniques. Here, we used pharmacological fMRI with a GABAergic benzodiazepine challenge to empirically generate a model of E:I ratio based on multivariate patterns of functional connectivity. In an independent sample of 879 youth (ages 8 to 22 years), this model predicted reductions in the E:I ratio during adolescence, which were specific to association cortex and related to psychopathology. These findings support hypothesized shifts in E:I balance of association cortices during a neurodevelopmental critical period in adolescence.

Spatial transcriptomics using combinatorial fluorescence spectral and lifetime encoding, imaging and analysis

Multiplexed mRNA profiling in the spatial context provides new information enabling basic research and clinical applications. Unfortunately, existing spatial transcriptomics methods are limited due to either low multiplexing or complexity. Here, we introduce a spatialomics technology, termed Multi Omic Single-scan Assay with Integrated Combinatorial Analysis (MOSAICA), that integrates in situ labeling of mRNA and protein markers in cells or tissues with combinatorial fluorescence spectral and lifetime encoded probes, spectral and time-resolved fluorescence imaging, and machine learning-based decoding. We demonstrate MOSAICA's multiplexing scalability in detecting 10-plex targets in fixed colorectal cancer cells using combinatorial labeling of five fluorophores with facile error-detection and removal of autofluorescence. MOSAICA's analysis is strongly correlated with sequencing data (Pearson's r = 0.96) and was further benchmarked using RNAscopeTM and LGC StellarisTM. We further apply MOSAICA for multiplexed analysis of clinical melanoma Formalin-Fixed Paraffin-Embedded (FFPE) tissues. We finally demonstrate simultaneous co-detection of protein and mRNA in cancer cells.

APOE mediated neuroinflammation and neurodegeneration in Alzheimer's disease

Neuroinflammation is a central mechanism involved in neurodegeneration as observed in Alzheimer's disease (AD), the most prevalent form of neurodegenerative disease. Apolipoprotein E4 (APOE4), the strongest genetic risk factor for AD, directly influences disease onset and progression by interacting with the major pathological hallmarks of AD including amyloid-β plaques, neurofibrillary tau tangles, as well as neuroinflammation. Microglia and astrocytes, the two major immune cells in the brain, exist in an immune-vigilant state providing immunological defense as well as housekeeping functions that promote neuronal well-being. It is becoming increasingly evident that under disease conditions, these immune cells become progressively dysfunctional in regulating metabolic and immunoregulatory pathways, thereby promoting chronic inflammation-induced neurodegeneration. Here, we review and discuss how APOE and specifically APOE4 directly influences amyloid-β and tau pathology, and disrupts microglial as well as astroglial immunomodulating functions leading to chronic inflammation that contributes to neurodegeneration in AD.

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.

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.

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.

The ascending arousal system promotes optimal performance through mesoscale network integration in a visuospatial attentional task

Previous research has shown that the autonomic nervous system provides essential constraints over ongoing cognitive function. However, there is currently a relative lack of direct empirical evidence for how this interaction manifests in the brain at the macroscale level. Here, we examine the role of ascending arousal and attentional load on large-scale network dynamics by combining pupillometry, functional MRI, and graph theoretical analysis to analyze data from a visual motion-tracking task with a parametric load manipulation. We found that attentional load effects were observable in measures of pupil diameter and in a set of brain regions that parametrically modulated their BOLD activity and mesoscale network-level integration. In addition, the regional patterns of network reconfiguration were correlated with the spatial distribution of the α2a adrenergic receptor. Our results further solidify the relationship between ascending noradrenergic activity, large-scale network integration, and cognitive task performance.

VGF as a biomarker and therapeutic target in neurodegenerative and psychiatric diseases

Neurosecretory protein VGF (non-acronymic) belongs to the granin family of neuropeptides. VGF and VGF-derived peptides have been repeatedly identified in well-powered and well-designed multi-omic studies as dysregulated in neurodegenerative and psychiatric diseases. New therapeutics is urgently needed for these devastating and costly diseases, as are new biomarkers to improve disease diagnosis and mechanistic understanding. From a list of 537 genes involved in Alzheimer's disease pathogenesis, VGF was highlighted by the Accelerating Medicines Partnership in Alzheimer's disease as the potential therapeutic target of greatest interest. VGF levels are consistently decreased in brain tissue and CSF samples from patients with Alzheimer's disease compared to controls, and its levels correlate with disease severity and Alzheimer's disease pathology. In the brain, VGF exists as multiple functional VGF-derived peptides. Full-length human VGF1-615 undergoes proteolytic processing by prohormone convertases and other proteases in the regulated secretory pathway to produce at least 12 active VGF-derived peptides. In cell and animal models, these VGF-derived peptides have been linked to energy balance regulation, neurogenesis, synaptogenesis, learning and memory, and depression-related behaviours throughout development and adulthood. The C-terminal VGF-derived peptides, TLQP-62 (VGF554-615) and TLQP-21 (VGF554-574) have differential effects on Alzheimer's disease pathogenesis, neuronal and microglial activity, and learning and memory. TLQP-62 activates neuronal cell-surface receptors and regulates long-term hippocampal memory formation. TLQP-62 also prevents immune-mediated memory impairment, depression-like and anxiety-like behaviours in mice. TLQP-21 binds to microglial cell-surface receptors, triggering microglial chemotaxis and phagocytosis. These actions were reported to reduce amyloid-β plaques and decrease neuritic dystrophy in a transgenic mouse model of familial Alzheimer's disease. Expression differences of VGF-derived peptides have also been associated with frontotemporal lobar dementias, amyotrophic lateral sclerosis, Lewy body diseases, Huntington's disease, pain, schizophrenia, bipolar disorder, depression and antidepressant response. This review summarizes current knowledge and highlights questions for future investigation regarding the roles of VGF and its dysregulation in neurodegenerative and psychiatric disease. Finally, the potential of VGF and VGF-derived peptides as biomarkers and novel therapeutic targets for neurodegenerative and psychiatric diseases is highlighted.

Neuropsychological evidence of multi-domain network hubs in the human thalamus

Hubs in the human brain support behaviors that arise from brain network interactions. Previous studies have identified hub regions in the human thalamus that are connected with multiple functional networks. However, the behavioral significance of thalamic hubs has yet to be established. Our framework predicts that thalamic subregions with strong hub properties are broadly involved in functions across multiple cognitive domains. To test this prediction, we studied human patients with focal thalamic lesions in conjunction with network analyses of the human thalamocortical functional connectome. In support of our prediction, lesions to thalamic subregions with stronger hub properties were associated with widespread deficits in executive, language, and memory functions, whereas lesions to thalamic subregions with weaker hub properties were associated with more limited deficits. These results highlight how a large-scale network model can broaden our understanding of thalamic function for human cognition.

Microglial activation and tau propagate jointly across Braak stages

Compelling experimental evidence suggests that microglial activation is involved in the spread of tau tangles over the neocortex in Alzheimer's disease (AD). We tested the hypothesis that the spatial propagation of microglial activation and tau accumulation colocalize in a Braak-like pattern in the living human brain. We studied 130 individuals across the aging and AD clinical spectrum with positron emission tomography brain imaging for microglial activation ([¹¹C]PBR28), amyloid-β (Aβ) ([¹⁸F]AZD4694) and tau ([¹⁸F]MK-6240) pathologies. We further assessed microglial triggering receptor expressed on myeloid cells 2 (TREM2) cerebrospinal fluid (CSF) concentrations and brain gene expression patterns. We found that [¹¹C]PBR28 correlated with CSF soluble TREM2 and showed regional distribution resembling TREM2 gene expression. Network analysis revealed that microglial activation and tau correlated hierarchically with each other following Braak-like stages. Regression analysis revealed that the longitudinal tau propagation pathways depended on the baseline microglia network rather than the tau network circuits. The co-occurrence of Aβ, tau and microglia abnormalities was the strongest predictor of cognitive impairment in our study population. Our findings support a model where an interaction between Aβ and activated microglia sets the pace for tau spread across Braak stages.

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.

Computational models link cellular mechanisms of neuromodulation to large-scale neural dynamics

Decades of neurobiological research have disclosed the diverse manners in which the response properties of neurons are dynamically modulated to support adaptive cognitive functions. This neuromodulation is achieved through alterations in the biophysical properties of the neuron. However, changes in cognitive function do not arise directly from the modulation of individual neurons, but are mediated by population dynamics in mesoscopic neural ensembles. Understanding this multiscale mapping is an important but nontrivial issue. Here, we bridge these different levels of description by showing how computational models parametrically map classic neuromodulatory processes onto systems-level models of neural activity. The ensuing critical balance of systems-level activity supports perception and action, although our knowledge of this mapping remains incomplete. In this way, quantitative models that link microscale neuronal neuromodulation to systems-level brain function highlight gaps in knowledge and suggest new directions for integrating theoretical and experimental work.

Integrated transcriptomic and neuroimaging brain model decodes biological mechanisms in aging and Alzheimer's disease

Both healthy aging and Alzheimer's disease (AD) are characterized by concurrent alterations in several biological factors. However, generative brain models of aging and AD are limited in incorporating the measures of these biological factors at different spatial resolutions. Here, we propose a personalized bottom-up spatiotemporal brain model that accounts for the direct interplay between hundreds of RNA transcripts and multiple macroscopic neuroimaging modalities (PET, MRI). In normal elderly and AD participants, the model identifies top genes modulating tau and amyloid-β burdens, vascular flow, glucose metabolism, functional activity, and atrophy to drive cognitive decline. The results also revealed that AD and healthy aging share specific biological mechanisms, even though AD is a separate entity with considerably more altered pathways. Overall, this personalized model offers novel insights into the multiscale alterations in the elderly brain, with important implications for identifying effective genetic targets for extending healthy aging and treating AD progression.

Differential patterns of gray matter volumes and associated gene expression profiles in cognitively-defined Alzheimer's disease subgroups

The clinical presentation of Alzheimer's disease (AD) varies widely across individuals but the neurobiological mechanisms underlying this heterogeneity are largely unknown. Here, we compared regional gray matter (GM) volumes and associated gene expression profiles between cognitively-defined subgroups of amyloid-β positive individuals clinically diagnosed with AD dementia (age: 66 ± 7, 47% male, MMSE: 21 ± 5). All participants underwent neuropsychological assessment with tests covering memory, executive-functioning, language and visuospatial-functioning domains. Subgroup classification was achieved using a psychometric framework that assesses which cognitive domain shows substantial relative impairment compared to the intra-individual average across domains, which yielded the following subgroups in our sample; AD-Memory (n = 41), AD-Executive (n = 117), AD-Language (n = 33), AD-Visuospatial (n = 171). We performed voxel-wise contrasts of GM volumes derived from 3Tesla structural MRI between subgroups and controls (n = 127, age 58 ± 9, 42% male, MMSE 29 ± 1), and observed that differences in regional GM volumes compared to controls closely matched the respective cognitive profiles. Specifically, we detected lower medial temporal lobe GM volumes in AD-Memory, lower fronto-parietal GM volumes in AD-Executive, asymmetric GM volumes in the temporal lobe (left < right) in AD-Language, and lower GM volumes in posterior areas in AD-Visuospatial. In order to examine possible biological drivers of these differences in regional GM volumes, we correlated subgroup-specific regional GM volumes to brain-wide gene expression profiles based on a stereotactic characterization of the transcriptional architecture of the human brain as provided by the Allen human brain atlas. Gene-set enrichment analyses revealed that variations in regional expression of genes involved in processes like mitochondrial respiration and metabolism of proteins were associated with patterns of regional GM volume across multiple subgroups. Other gene expression vs GM volume-associations were only detected in particular subgroups, e.g., genes involved in the cell cycle for AD-Memory, specific sets of genes related to protein metabolism in AD-Language, and genes associated with modification of gene expression in AD-Visuospatial. We conclude that cognitively-defined AD subgroups show neurobiological differences, and distinct biological pathways may be involved in the emergence of these differences.

Integration of human whole-brain transcriptome and neuroimaging data: Practical considerations of current available methods

The Allen Human Brain Atlas (AHBA) is the first example of human brain transcriptomic mappings and detailed anatomical annotations which, for the first time, has allowed the integration of human brain transcriptomics with neuroimaging. This has been made possible because the AHBA offered an original dataset that contains mRNA expression measures for >20,000 genes covering the whole brain and, critically, in a standard stereotaxic space. In recent years many different methods have been used to integrate this data set with brain imaging data, although this endeavour has lacked harmony in terms of the workflow of data processing and subsequent analyses. In this work we discuss five main issues that experience has highlighted as in need of thorough consideration when integrating the AHBA with neuroimaging. These concerns are corroborated by comparing the performance of three different publicly available methods in correlating the same measures of serotonin receptors density with the correspondent AHBA mRNA maps. In this representative case, we were able to show how these methods can lead to discrepant results, suggesting that processing options are not neutral. We believe that the field should take into serious consideration these issues as they could undermine reproducibility and, in the end, the intrinsic value of the AHBA. We also advise on possible strategies to overcome these discrepancies. Finally, we encourage authors towards practices that will improve reproducibility such as transparency in reporting, algorithm and data sharing, collaboration.

Enrichment of Disease-Associated Genes in Cortical Areas Defined by Transcriptome-Based Parcellation

BACKGROUND

Parcellation of the cerebral cortex serves the investigation of the emergence of uniquely human brain functions and disorders. Transcriptome data enable the characterization of the molecular properties of cortical areas in unprecedented detail. Previously, we predicted the expression of 18,686 genes in the entire human brain based on microarray data. Here, we employed these data to parcellate the cortex and study the regional enrichment of disease-associated genes.

METHODS

We performed agglomerative hierarchical clustering based on normalized transcriptome data to delineate areas with distinct gene expression profiles. Subsequently, we tested these profiles for the enrichment of gene sets associated with brain disorders by genome-wide association studies and expert-curated databases using gene set enrichment analysis.

RESULTS

Transcriptome-based parcellation identified borders in line with major anatomical landmarks and the functional differentiation of primary motor, somatosensory, visual, and auditory areas. Gene set enrichment analysis based on curated databases suggested new roles of specific areas in psychiatric and neurological disorders while reproducing well-established links for movement and neurodegenerative disorders, for example, amyotrophic lateral sclerosis (motor cortex) and Alzheimer's disease (entorhinal cortex). Meanwhile, gene sets derived from genome-wide association studies on psychiatric disorders exhibited similar enrichment patterns driven by pleiotropic genes expressed in the posterior fusiform gyrus and inferior parietal lobule.

CONCLUSIONS

The identified enrichment patterns suggest the vulnerability of specific cortical areas to various influences that might alter the risk of developing one or several brain disorders. For several diseases, specific genes were highlighted, which could lead to the discovery of novel disease mechanisms and urgently needed treatments.

Association of dopamine D2/3 receptor binding potential measured using PET and [11C]-(+)-PHNO with post-mortem DRD2/3 gene expression in the human brain

Open access post-mortem transcriptome atlases such as the Allen Human Brain Atlas (AHBA) can inform us about mRNA expression of numerous proteins of interest across the whole brain, while in vivo protein binding in the human brain can be quantified by means of neuroreceptor positron emission tomography (PET). By combining both modalities, the association between regional gene expression and receptor distribution in the living brain can be approximated. Here, we compare the characteristics of D2 and D3 dopamine receptor distribution by applying the dopamine D2/3 receptor agonist radioligand [11C]-(+)-PHNO and human gene expression data. Since [11C]-(+)-PHNO has a higher affinity for D3 compared to D2 receptors, we hypothesized that there is a stronger relationship between D2/3 non-displaceable binding potentials (BPND) and D3 mRNA expression. To investigate the relationship between D2/3 BPND and mRNA expression of DRD2 and DRD3 we performed [11C]-(+)-PHNO PET scans in 27 healthy subjects (12 females) and extracted gene expression data from the AHBA. We also calculated D2/D3 mRNA expression ratios to imitate the mixed D2/3 signal of [11C]-(+)-PHNO. In accordance with our a priori hypothesis, a strong correlation between [11C]-(+)-PHNO and DRD3 expression was found. However, there was no significant correlation with DRD2 expression. Calculated D2/D3 mRNA expression ratios also showed a positive correlation with [11C]-(+)-PHNO binding, reflecting the mixed D2/3 signal of the radioligand. Our study supports the usefulness of combining gene expression data from open access brain atlases with in vivo imaging data in order to gain more detailed knowledge on neurotransmitter signaling.

Neuronal timescales are functionally dynamic and shaped by cortical microarchitecture

Complex cognitive functions such as working memory and decision-making require information maintenance over seconds to years, from transient sensory stimuli to long-term contextual cues. While theoretical accounts predict the emergence of a corresponding hierarchy of neuronal timescales, direct electrophysiological evidence across the human cortex is lacking. Here, we infer neuronal timescales from invasive intracranial recordings. Timescales increase along the principal sensorimotor-to-association axis across the entire human cortex, and scale with single-unit timescales within macaques. Cortex-wide transcriptomic analysis shows direct alignment between timescales and expression of excitation- and inhibition-related genes, as well as genes specific to voltage-gated transmembrane ion transporters. Finally, neuronal timescales are functionally dynamic: prefrontal cortex timescales expand during working memory maintenance and predict individual performance, while cortex-wide timescales compress with aging. Thus, neuronal timescales follow cytoarchitectonic gradients across the human cortex and are relevant for cognition in both short and long terms, bridging microcircuit physiology with macroscale dynamics and behavior.

APOE4 is associated with cognitive and pathological heterogeneity in patients with Alzheimer's disease: a systematic review

Possession of the ε4 allele of apolipoprotein E (APOE) is the primary genetic risk factor for the sporadic form of Alzheimer’s disease (AD). While researchers have extensively characterized the impact that APOE ε4 (APOE4) has on the susceptibility of AD, far fewer studies have investigated the phenotypic differences of patients with AD who are APOE4 carriers vs. those who are non-carriers. In order to understand these differences, we performed a qualitative systematic literature review of the reported cognitive and pathological differences between APOE4-positive (APOE4+) vs. APOE4-negative (APOE4−) AD patients. The studies performed on this topic to date suggest that APOE4 is not only an important mediator of AD susceptibility, but that it likely confers specific phenotypic heterogeneity in AD presentation, as well. Specifically, APOE4+ AD patients appear to possess more tau accumulation and brain atrophy in the medial temporal lobe, resulting in greater memory impairment, compared to APOE4− AD patients. On the other hand, APOE4− AD patients appear to possess more tau accumulation and brain atrophy in the frontal and parietal lobes, resulting in greater impairment in executive function, visuospatial abilities, and language, compared to APOE4+ AD patients. Although more work is necessary to validate and interrogate these findings, these initial observations of pathological and cognitive heterogeneity between APOE4+ vs. APOE4− AD patients suggest that there is a fundamental divergence in AD manifestation related to APOE genotype, which may have important implications in regard to the therapeutic treatment of these two patient populations.

Generative modeling of brain maps with spatial autocorrelation

Studies of large-scale brain organization have revealed interesting relationships between spatial gradients in brain maps across multiple modalities. Evaluating the significance of these findings requires establishing statistical expectations under a null hypothesis of interest. Through generative modeling of synthetic data that instantiate a specific null hypothesis, quantitative benchmarks can be derived for arbitrarily complex statistical measures. Here, we present a generative null model, provided as an open-access software platform, that generates surrogate maps with spatial autocorrelation (SA) matched to SA of a target brain map. SA is a prominent and ubiquitous property of brain maps that violates assumptions of independence in conventional statistical tests. Our method can simulate surrogate brain maps, constrained by empirical data, that preserve the SA of cortical, subcortical, parcellated, and dense brain maps. We characterize how SA impacts p-values in pairwise brain map comparisons. Furthermore, we demonstrate how SA-preserving surrogate maps can be used in gene set enrichment analyses to test hypotheses of interest related to brain map topography. Our findings demonstrate the utility of SA-preserving surrogate maps for hypothesis testing in complex statistical analyses, and underscore the need to disambiguate meaningful relationships from chance associations in studies of large-scale brain organization.

Dysfunction of the Blood-Brain Barrier-A Key Step in Neurodegeneration and Dementia

The vascular endothelium in the brain is an essential part of the blood-brain-barrier (BBB) because of its very tight structure to secure a functional and molecular separation of the brain from the rest of the body and to protect neurons from pathogens and toxins. Impaired transport of metabolites across the BBB due to its increasing dysfunction affects brain health and cognitive functioning, thus providing a starting point of neurodegenerative diseases. The term “cerebral metabolic syndrome” is proposed to highlight the importance of lifestyle factors in neurodegeneration and to describe the impact of increasing BBB dysfunction on neurodegeneration and dementia, especially in elderly patients. If untreated, the cerebral metabolic syndrome may evolve into dementia. Due to the high energy demand of the brain, impaired glucose transport across the BBB via glucose transporters as GLUT1 renders the brain increasingly susceptible to neurodegeneration. Apoptotic processes are further supported by the lack of essential metabolites of the phosphocholine synthesis. In Alzheimer’s disease (AD), inflammatory and infectious processes at the BBB increase the dysfunction and might be pace-making events. At this point, the potentially highly relevant role of the thrombocytic amyloid precursor protein (APP) in endothelial inflammation of the BBB is discussed. Chronic inflammatory processes of the BBB transmitted to an increasing number of brain areas might cause a lasting build-up of spreading, pore-forming β-amyloid fragments explaining the dramatic progression of the disease. In the view of the essential requirement of an early diagnosis to investigate and implement causal therapeutic strategies against dementia, brain imaging methods are of great importance. Therefore, status and opportunities in the field of diagnostic imaging of the living human brain will be portrayed, comprising diverse techniques such as positron emissions tomography (PET) and functional magnetic resonance imaging (fMRI) to uncover the patterns of atrophy, protein deposits, hypometabolism, and molecular as well as functional alterations in AD.

Association of Apolipoprotein E ε4 With Medial Temporal Tau Independent of Amyloid-β

Importance

Apolipoprotein E ε4 (APOEε4) is the single most important genetic risk factor for Alzheimer disease. While APOEε4 is associated with increased amyloid-β burden, its association with cerebral tau pathology has been controversial.

Objective

To determine whether APOEε4 is associated with medial temporal tau pathology independently of amyloid-β, sex, clinical status, and age.

Design, Setting, and Participants

This is a study of 2 cross-sectional cohorts of volunteers who were cognitively normal, had mild cognitive impairment (MCI), or had Alzheimer disease dementia: the Translational Biomarkers in Aging and Dementia (TRIAD) study (data collected between October 2017 and July 2019) and the Alzheimer's Disease Neuroimaging Initiative (ADNI) (collected between November 2015 and June 2019). The first cohort (TRIAD) comprised cognitively normal elderly participants (n = 124), participants with MCI (n = 50), and participants with Alzheimer disease (n = 50) who underwent tau positron emission tomography (PET) with fluorine 18-labeled MK6240 and amyloid-β PET with [18F]AZD4694. The second sample (ADNI) was composed of cognitively normal elderly participants (n = 157), participants with MCI (n = 83), and participants with Alzheimer disease (n = 25) who underwent tau PET with [18F]flortaucipir and amyloid-β PET with [18F]florbetapir. Exclusion criteria were a history of other neurological disorders, stroke, or head trauma. There were 489 eligible participants, selected based on availability of amyloid-PET, tau-PET, magnetic resonance imaging, and genotyping for APOEε4. Forty-five young adults (<30 years) from the TRIAD cohort were not selected for this study.

Main Outcomes and Measures

A main association between APOEε4 and tau-PET standardized uptake value ratio, correcting for age, sex, clinical status, and neocortical amyloid-PET standardized uptake value ratio.

Results

The mean (SD) age of the 489 participants was 70.5 (7.1) years; 171 were APOEε4 carriers (34.9%), and 230 of 489 were men. In both cohorts, APOEε4 was associated in increased tau-PET uptake in the entorhinal cortex and hippocampus independently of amyloid-β, sex, age, and clinical status after multiple comparisons correction (TRIAD: β = 0.33; 95% CI, 0.19-0.49; ADNI: β = 0.13; 95% CI, 0.08-0.19; P < .001).

Conclusions and Relevance

Our results indicate that the elevated risk of developing dementia conferred by APOEε4 genotype involves mechanisms associated with both amyloid-β and tau aggregation. These results contribute to an evolving framework in which APOEε4 has deleterious consequences in Alzheimer disease beyond its link with amyloid-β and suggest APOEε4 as a potential target for future disease-modifying therapeutic trials targeting tau pathology.

Effects of route of administration on oxytocin-induced changes in regional cerebral blood flow in humans

Could nose-to-brain pathways mediate the effects of peptides such as oxytocin (OT) on brain physiology when delivered intranasally? We address this question by contrasting two methods of intranasal administration (a standard nasal spray, and a nebulizer expected to improve OT deposition in nasal areas putatively involved in direct nose-to-brain transport) to intravenous administration in terms of effects on regional cerebral blood flow during two hours post-dosing. We demonstrate that OT-induced decreases in amygdala perfusion, a key hub of the OT central circuitry, are explained entirely by OT increases in systemic circulation following both intranasal and intravenous OT administration. Yet we also provide robust evidence confirming the validity of the intranasal route to target specific brain regions. Our work has important translational implications and demonstrates the need to carefully consider the method of administration in our efforts to engage specific central oxytocinergic targets for the treatment of neuropsychiatric disorders.

Topologically Guided Prioritization of Candidate Gene Transcripts Coexpressed with the 5-HT1A Receptor by Combining In Vivo PET and Allen Human Brain Atlas Data

The serotonin-1A receptor (5-HT_1AR) represents a viable target in the treatment of disorders of the brain. However, development of psychiatric drugs continues to be hindered by the relative inaccessibility of brain tissue. Although the efficacy of drugs selective for the 5-HT_1AR has not been proven, research continues to focus on drugs that influence this receptor subtype. To further knowledge on this topic, we investigated the topological coexpression patterns of the 5-HT_1AR. We calculated Spearman’s rho for the correlation of positron emission tomography-binding potentials (BP_ND) of the 5-HT_1AR assessed in 30 healthy subjects using the tracer [carbonyl-¹¹C]WAY-100635 and predicted whole-brain mRNA expression of 18 686 genes. After applying a threshold of r > 0.3 in a leave-one-out cross-validation of the prediction of mRNA expression, genes with ρ ≥ 0.7 were considered to be relevant. In cortical regions, 199 genes showed high correlation with the BP_ND of the 5-HT_1AR, in subcortical regions 194 genes. Using our approach, we could consolidate the role of BDNF and implicate new genes (AnxA8, NeuroD2) in serotonergic functioning. Despite its explorative nature, the analysis can be seen as a gene prioritization approach to reduce the number of genes potentially connected to 5-HT_1AR functioning and guide future in vitro studies.