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

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

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

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

Enhanced cerebral blood flow similarity of the somatomotor network in chronic insomnia: Transcriptomic decoding, gut microbial signatures and phenotypic roles

Chronic insomnia (CI) is a complex disease involving multiple factors including genetics, gut microbiota, and brain structure and function. However, there lacks a unified framework to elucidate how these factors interact in CI. By combining data of clinical assessment, sleep behavior recording, cognitive test, multimodal MRI (structural, functional, and perfusion), gene, and gut microbiota, this study demonstrated that enhanced cerebral blood flow (CBF) similarities of the somatomotor network (SMN) acted as a key mediator to link multiple factors in CI. Specifically, we first demonstrated that only CBF but not morphological or functional networks exhibited alterations in patients with CI, characterized by increases within the SMN and between the SMN and higher-order associative networks. Moreover, these findings were highly reproducible and the CBF similarity method was test-retest reliable. Further, we showed that transcriptional profiles explained 60.4 % variance of the pattern of the increased CBF similarities with the most correlated genes enriched in regulation of cellular and protein localization and material transport, and gut microbiota explained 69.7 % inter-individual variance in the increased CBF similarities with the most contributions from Negativicutes and Lactobacillales. Finally, we found that the increased CBF similarities were correlated with clinical variables, accounted for sleep behaviors and cognitive deficits, and contributed the most to the patient-control classification (accuracy = 84.4 %). Altogether, our findings have important implications for understanding the neuropathology of CI and may inform ways of developing new therapeutic strategies for the disease.

High-resolution diffusion magnetic resonance imaging and spatial-transcriptomic in developing mouse brain

Brain development is a highly complex process regulated by numerous genes at the molecular and cellular levels. Brain tissue exhibits serial microstructural changes during the development process. High-resolution diffusion magnetic resonance imaging (dMRI) affords a unique opportunity to probe these changes in the developing brain non-destructively. In this study, we acquired multi-shell dMRI datasets at 32 µm isotropic resolution to investigate the tissue microstructure alterations, which we believe to be the highest spatial resolution dMRI datasets obtained for postnatal mouse brains. We adapted the Allen Developing Mouse Brain Atlas (ADMBA) to integrate quantitative MRI metrics and spatial transcriptomics. Diffusion tensor imaging (DTI), diffusion kurtosis imaging (DKI), and neurite orientation dispersion and density imaging (NODDI) metrics were used to quantify brain development at different postnatal days. We demonstrated that the differential evolutions of fiber orientation distributions contribute to the distinct development patterns in white matter (WM) and gray matter (GM). Furthermore, the genes enriched in the nervous system that regulate brain structure and function were expressed in spatial correlation with age-matched dMRI. This study is the first one providing high-resolution dMRI, including DTI, DKI, and NODDI models, to trace mouse brain microstructural changes in WM and GM during postnatal development. This study also highlighted the genotype-phenotype correlation of spatial transcriptomics and dMRI, which may improve our understanding of brain microstructure changes at the molecular level.

Age-dependent functional development pattern in neonatal brain: An fMRI-based brain entropy study

The relationship between brain entropy (BEN) and early brain development has been established through animal studies. However, it remains unclear whether the BEN can be used to identify age-dependent functional changes in human neonatal brains and the genetic underpinning of the new neuroimaging marker remains to be elucidated. In this study, we analyzed resting-state fMRI data from the Developing Human Connectome Project, including 280 infants who were scanned at 37.5-43.5 weeks postmenstrual age. The BEN maps were calculated for each subject, and a voxel-wise analysis was conducted using a general linear model to examine the effects of age, sex, and preterm birth on BEN. Additionally, we evaluated the correlation between regional BEN and gene expression levels. Our results demonstrated that the BEN in the sensorimotor-auditory and association cortices, along the 'S-A' axis, was significantly positively correlated with postnatal age (PNA), and negatively correlated with gestational age (GA), respectively. Meanwhile, the BEN in the right rolandic operculum correlated significantly with both GA and PNA. Preterm-born infants exhibited increased BEN values in widespread cortical areas, particularly in the visual-motor cortex, when compared to term-born infants. Moreover, we identified five BEN-related genes (DNAJC12, FIG4, STX12, CETN2, and IRF2BP2), which were involved in protein folding, synaptic vesicle transportation and cell division. These findings suggest that the fMRI-based BEN can serve as an indicator of age-dependent brain functional development in human neonates, which may be influenced by specific genes.

Gray matter volume and corresponding covariance connectivity are biomarkers for major depressive disorder

The major depressive disorder (MDD) is a common and severe mental disorder. To identify a reliable biomarker for MDD is important for early diagnosis and prevention. Given easy access and high reproducibility, the structural magnetic resonance imaging (sMRI) is an ideal method to identify the biomarker for depression. In this study, sMRI data of first episode, treatment-naïve 66 MDD patients and 54 sex-, age-, and education-matched healthy controls (HC) were used to identify the differences in gray matter volume (GMV), group-level, individual-level covariance connections. Finally, the abnormal GMV and individual covariance connections were applied to classify MDD from HC. MDD patients showed higher GMV in middle occipital gyrus (MOG) and precuneus (PCun), and higher structural covariance connections between MOG and PCun. In addition, the Allen Human Brain Atlas (AHBA) was applied and revealed the genetic basis for the changes of gray matter volume. Importantly, we reported that GMV in MOG, PCun and structural covariance connectivity between MOG and PCun are able to discriminate MDD from HC. Our results revealed structural underpinnings for MDD, which may contribute towards early discriminating for depression.

Chronic hypercortisolism disrupts the principal functional gradient in Cushing's disease: A multi-scale connectomics and transcriptomics study

Cushing's disease (CD) represents a state of cortisol excess, serving as a model to investigate the effects of prolonged hypercortisolism on functional brain. Potential alterations in the functional connectome of the brain may explain frequently reported cognitive deficits and affective disorders in CD patients. This study aims to elucidate the effects of chronic hypercortisolism on the principal functional gradient, which represents a hierarchical architecture with gradual transitions across cognitive processes, by integrating connectomics and transcriptomics approaches. Utilizing resting-state functional magnetic resonance imaging data from 140 participants (86 CD patients, 54 healthy controls) recruited at a single center, we explored the alterations in the principal gradient in CD patients. Further, we thoroughly explored the underlying associative mechanisms of the observed characteristic alterations with cognitive function domains, biological attributes, and neuropsychiatric representations, as well as gene expression profiles. Compared to healthy controls, CD patients demonstrated changes in connectome patterns in both primary and higher-order networks, exhibiting an overall converged trend along the principal gradient axis. The gradient values in CD patients' right prefrontal cortex and bilateral sensorimotor cortices exhibited a significant correlation with cortisol levels. Moreover, the cortical regions showing gradient alterations were principally associated with sensory information processing and higher-cognitive functions, as well as correlated with the gene expression patterns which involved synaptic components and function. The findings suggest that converged alterations in the principal gradient in CD patients may mediate the relationship between hypercortisolism and cognitive impairments, potentially involving genes regulating synaptic components and function.

Transcriptomic decoding of regional cortical vulnerability to major depressive disorder

Previous studies in small samples have identified inconsistent cortical abnormalities in major depressive disorder (MDD). Despite genetic influences on MDD and the brain, it is unclear how genetic risk for MDD is translated into spatially patterned cortical vulnerability. Here, we initially examined voxel-wise differences in cortical function and structure using the largest multi-modal MRI data from 1660 MDD patients and 1341 controls. Combined with the Allen Human Brain Atlas, we then adopted transcription-neuroimaging spatial correlation and the newly developed ensemble-based gene category enrichment analysis to identify gene categories with expression related to cortical changes in MDD. Results showed that patients had relatively circumscribed impairments in local functional properties and broadly distributed disruptions in global functional connectivity, consistently characterized by hyper-function in associative areas and hypo-function in primary regions. Moreover, the local functional alterations were correlated with genes enriched for biological functions related to MDD in general (e.g., endoplasmic reticulum stress, mitogen-activated protein kinase, histone acetylation, and DNA methylation); and the global functional connectivity changes were associated with not only MDD-general, but also brain-relevant genes (e.g., neuron, synapse, axon, glial cell, and neurotransmitters). Our findings may provide important insights into the transcriptomic signatures of regional cortical vulnerability to MDD.

Unraveling multi-scale neuroimaging biomarkers and molecular foundations for schizophrenia: A combined multivariate pattern analysis and transcriptome-neuroimaging association study

AIMS

Schizophrenia is characterized by alterations in resting-state spontaneous brain activity; however, it remains uncertain whether variations at diverse spatial scales are capable of effectively distinguishing patients from healthy controls. Additionally, the genetic underpinnings of these alterations remain poorly elucidated. We aimed to address these questions in this study to gain better understanding of brain alterations and their underlying genetic factors in schizophrenia.

METHODS

A cohort of 103 individuals with diagnosed schizophrenia and 110 healthy controls underwent resting-state functional MRI scans. Spontaneous brain activity was assessed using the regional homogeneity (ReHo) metric at four spatial scales: voxel-level (Scale 1) and regional-level (Scales 2-4: 272, 53, 17 regions, respectively). For each spatial scale, multivariate pattern analysis was performed to classify schizophrenia patients from healthy controls, and a transcriptome-neuroimaging association analysis was performed to establish connections between gene expression data and ReHo alterations in schizophrenia.

RESULTS

The ReHo metrics at all spatial scales effectively discriminated schizophrenia from healthy controls. Scale 2 showed the highest classification accuracy at 84.6%, followed by Scale 1 (83.1%) and Scale 3 (78.5%), while Scale 4 exhibited the lowest accuracy (74.2%). Furthermore, the transcriptome-neuroimaging association analysis showed that there were not only shared but also unique enriched biological processes across the four spatial scales. These related biological processes were mainly linked to immune responses, inflammation, synaptic signaling, ion channels, cellular development, myelination, and transporter activity.

CONCLUSIONS

This study highlights the potential of multi-scale ReHo as a valuable neuroimaging biomarker in the diagnosis of schizophrenia. By elucidating the complex molecular basis underlying the ReHo alterations of this disorder, this study not only enhances our understanding of its pathophysiology, but also pave the way for future advancements in genetic diagnosis and treatment of schizophrenia.

Compressed representation of brain genetic transcription

The architecture of the brain is too complex to be intuitively surveyable without the use of compressed representations that project its variation into a compact, navigable space. The task is especially challenging with high-dimensional data, such as gene expression, where the joint complexity of anatomical and transcriptional patterns demands maximum compression. The established practice is to use standard principal component analysis (PCA), whose computational felicity is offset by limited expressivity, especially at great compression ratios. Employing whole-brain, voxel-wise Allen Brain Atlas transcription data, here we systematically compare compressed representations based on the most widely supported linear and non-linear methods-PCA, kernel PCA, non-negative matrix factorisation (NMF), t-stochastic neighbour embedding (t-SNE), uniform manifold approximation and projection (UMAP), and deep auto-encoding-quantifying reconstruction fidelity, anatomical coherence, and predictive utility across signalling, microstructural, and metabolic targets, drawn from large-scale open-source MRI and PET data. We show that deep auto-encoders yield superior representations across all metrics of performance and target domains, supporting their use as the reference standard for representing transcription patterns in the human brain.

Processing Language Partly Shares Neural Genetic Basis with Processing Tools and Body Parts

Language is an evolutionarily salient faculty for humans that relies on a distributed brain network spanning across frontal, temporal, parietal, and subcortical regions. To understand whether the complex language network shares common or distinct genetic mechanisms, we examined the relationships between the genetic effects underlying the brain responses to language and a set of object domains that have been suggested to coevolve with language: tools, faces (indicating social), and body parts (indicating social and gesturing). Analyzing the twin datasets released by the Human Connectome Project that had functional magnetic resonance imaging data from human twin subjects (monozygotic and dizygotic) undergoing language and working memory tasks contrasting multiple object domains (198 females and 144 males for the language task; 192 females and 142 males for the working memory task), we identified a set of cortical regions in the frontal and temporal cortices and subcortical regions whose activity to language was significantly genetically influenced. The heterogeneity of the genetic effects among these language clusters was corroborated by significant differences of the human gene expression profiles (Allen Human Brain Atlas dataset). Among them, the bilateral basal ganglia (mainly dorsal caudate) exhibited a common genetic basis for language, tool, and body part processing, and the right superior temporal gyrus exhibited a common genetic basis for language and tool processing across multiple types of analyses. These results uncovered the heterogeneous genetic patterns of language neural processes, shedding light on the evolution of language and its shared origins with tools and bodily functions.

Neural mechanism of non-adaptive cognitive emotion regulation in patients with non-suicidal self-injury

BACKGROUND

The incidence of non-suicidal self-injury (NSSI) has been on the rise in recent years. Studies have shown that people with NSSI have difficulties in emotion regulation and cognitive control. In addition, some studies have investigated the cognitive emotion regulation of people with NSSI which found that they have difficulties in cognitive emotion regulation, but there was a lack of research on cognitive emotion regulation strategies and related neural mechanisms.

METHODS

This study included 117 people with NSSI (age = 19.47 ± 5.13, male = 17) and 84 non-NSSI participants (age = 19.86 ± 4.14, male = 16). People with NSSI met the DSM-5 diagnostic criteria, and non-NSSI participants had no mental or physical disorders. The study collected all participants' data of Cognitive Emotion Regulation Questionnaire (CERQ) and functional magnetic resonance imaging (fMRI) to explore the differences in psychological performance and brain between two groups. Afterwards, Machine learning was used to select the found differential brain regions to obtain the highest correlation regions with NSSI. Then, Allen's Human Brain Atlas database was used to compare with the information on the abnormal brain regions of people with NSSI to find the genetic information related to NSSI. In addition, gene enrichment analysis was carried out to find the related pathways and specific cells that may have differences.

RESULTS

The differences between NSSI participants and non-NSSI participants were as follows: positive refocusing (t = -4.74, p < 0.01); refocusing on plans (t = -4.11, p < 0.01); positive reappraisal (t = -9.22, p < 0.01); self-blame (t = 6.30, p < 0.01); rumination (t = 3.64, p < 0.01); catastrophizing (t = 9.10, p < 0.01), and blaming others (t = 2.52, p < 0.01), the precentral gyrus (t = 6.04, pFDR < 0.05) and the rolandic operculum (t = -4.57, pFDR < 0.05). Rolandic operculum activity was negatively correlated with blaming others (r = -0.20, p < 0.05). Epigenetic results showed that excitatory neurons (p < 0.01) and inhibitory neurons (p < 0.01) were significant differences in two pathways, "trans-synaptic signaling" (p < -log108) and "modulation of chemical synaptic transmission" (p < -log108) in both cells.

CONCLUSIONS

People with NSSI are more inclined to adopt non-adaptive cognitive emotion regulation strategies. Rolandic operculum is also abnormally active. Abnormal changes in the rolandic operculum of them are associated with non-adaptive cognitive emotion regulation strategies. Changes in the excitatory and inhibitory neurons provide hints to explore the abnormalities of the neurological mechanisms at the cellular level of them. Trial registration number NCT04094623.

Relationship between brain white matter damage and grey matter atrophy in hereditary spastic paraplegia types 4 and 5

BACKGROUND AND PURPOSE

White matter (WM) damage is the main target of hereditary spastic paraplegia (HSP), but mounting evidence indicates that genotype-specific grey matter (GM) damage is not uncommon. Our aim was to identify and compare brain GM and WM damage patterns in HSP subtypes and investigate how gene expression contributes to these patterns, and explore the relationship between GM and WM damage.

METHODS

In this prospective single-centre cohort study from 2019 to 2022, HSP patients and controls underwent magnetic resonance imaging evaluations. The alterations of GM and WM patterns were compared between groups by applying a source-based morphometry approach. Spearman rank correlation was used to explore the associations between gene expression and GM atrophy patterns in HSP subtypes. Mediation analysis was conducted to investigate the interplay between GM and WM damage.

RESULTS

Twenty-one spastic paraplegia type 4 (SPG4) patients (mean age 50.7 years ± 12.0 SD, 15 men), 21 spastic paraplegia type 5 (SPG5) patients (mean age 29.1 years ± 12.8 SD, 14 men) and 42 controls (sex- and age-matched) were evaluated. Compared to controls, SPG4 and SPG5 showed similar WM damage but different GM atrophy patterns. GM atrophy patterns in SPG4 and SPG5 were correlated with corresponding gene expression (ρ = 0.30, p = 0.008, ρ = 0.40, p < 0.001, respectively). Mediation analysis indicated that GM atrophy patterns were mediated by WM damage in HSP.

CONCLUSIONS

Grey matter atrophy patterns were distinct between SPG4 and SPG5 and were not only secondary to WM damage but also associated with disease-related gene expression.

CLINICAL TRIAL REGISTRATION NO

NCT04006418.

A shifting role of thalamocortical connectivity in the emergence of cortical functional organization

The cortical patterning principle has been a long-standing question in neuroscience, yet how this translates to macroscale functional specialization in the human brain remains largely unknown. Here we examine age-dependent differences in resting-state thalamocortical connectivity to investigate its role in the emergence of large-scale functional networks during early life, using a primarily cross-sectional but also longitudinal approach. We show that thalamocortical connectivity during infancy reflects an early differentiation of sensorimotor networks and genetically influenced axonal projection. This pattern changes in childhood, when connectivity is established with the salience network, while decoupling externally and internally oriented functional systems. A developmental simulation using generative network models corroborated these findings, demonstrating that thalamic connectivity contributes to developing key features of the mature brain, such as functional segregation and the sensory-association axis, especially across 12-18 years of age. Our study suggests that the thalamus plays an important role in functional specialization during development, with potential implications for studying conditions with compromised internal and external processing.

Whole brain alignment of spatial transcriptomics between humans and mice with BrainAlign

The increasing utilization of mouse models in human neuroscience research places higher demands on computational methods to translate findings from the mouse brain to the human one. In this study, we develop BrainAlign, a self-supervised learning approach, for the whole brain alignment of spatial transcriptomics (ST) between humans and mice. BrainAlign encodes spots and genes simultaneously in two separated shared embedding spaces by a heterogeneous graph neural network. We demonstrate that BrainAlign could integrate cross-species spots into the embedding space and reveal the conserved brain regions supported by ST information, which facilitates the detection of homologous regions between humans and mice. Genomic analysis further presents gene expression connections between humans and mice and reveals similar expression patterns for marker genes. Moreover, BrainAlign can accurately map spatially similar homologous regions or clusters onto a unified spatial structural domain while preserving their relative positions.

Macaque Brainnetome Atlas: A multifaceted brain map with parcellation, connection, and histology

The rhesus macaque (Macaca mulatta) is a crucial experimental animal that shares many genetic, brain organizational, and behavioral characteristics with humans. A macaque brain atlas is fundamental to biomedical and evolutionary research. However, even though connectivity is vital for understanding brain functions, a connectivity-based whole-brain atlas of the macaque has not previously been made. In this study, we created a new whole-brain map, the Macaque Brainnetome Atlas (MacBNA), based on the anatomical connectivity profiles provided by high angular and spatial resolution ex vivo diffusion MRI data. The new atlas consists of 248 cortical and 56 subcortical regions as well as their structural and functional connections. The parcellation and the diffusion-based tractography were evaluated with invasive neuronal-tracing and Nissl-stained images. As a demonstrative application, the structural connectivity divergence between macaque and human brains was mapped using the Brainnetome atlases of those two species to uncover the genetic underpinnings of the evolutionary changes in brain structure. The resulting resource includes: (1) the thoroughly delineated Macaque Brainnetome Atlas (MacBNA), (2) regional connectivity profiles, (3) the postmortem high-resolution macaque diffusion and T2-weighted MRI dataset (Brainnetome-8), and (4) multi-contrast MRI, neuronal-tracing, and histological images collected from a single macaque. MacBNA can serve as a common reference frame for mapping multifaceted features across modalities and spatial scales and for integrative investigation and characterization of brain organization and function. Therefore, it will enrich the collaborative resource platform for nonhuman primates and facilitate translational and comparative neuroscience research.

A network control theory pipeline for studying the dynamics of the structural connectome

Network control theory (NCT) is a simple and powerful tool for studying how network topology informs and constrains the dynamics of a system. Compared to other structure-function coupling approaches, the strength of NCT lies in its capacity to predict the patterns of external control signals that may alter the dynamics of a system in a desired way. An interesting development for NCT in the neuroscience field is its application to study behavior and mental health symptoms. To date, NCT has been validated to study different aspects of the human structural connectome. NCT outputs can be monitored throughout developmental stages to study the effects of connectome topology on neural dynamics and, separately, to test the coherence of empirical datasets with brain function and stimulation. Here, we provide a comprehensive pipeline for applying NCT to structural connectomes by following two procedures. The main procedure focuses on computing the control energy associated with the transitions between specific neural activity states. The second procedure focuses on computing average controllability, which indexes nodes' general capacity to control the dynamics of the system. We provide recommendations for comparing NCT outputs against null network models, and we further support this approach with a Python-based software package called 'network control theory for python'. The procedures in this protocol are appropriate for users with a background in network neuroscience and experience in dynamical systems theory.

The myelin water imaging transcriptome: myelin water fraction regionally varies with oligodendrocyte-specific gene expression

Identifying sensitive and specific measures that can quantify myelin are instrumental in characterizing microstructural changes in neurological conditions. Neuroimaging transcriptomics is emerging as a valuable technique in this regard, offering insights into the molecular basis of promising candidates for myelin quantification, such as myelin water fraction (MWF). We aimed to demonstrate the utility of neuroimaging transcriptomics by validating MWF as a myelin measure. We utilized data from a normative MWF brain atlas, comprised of 50 healthy subjects (mean age = 25 years, range = 17-42 years) scanned at 3 Tesla. Magnetic resonance imaging data included myelin water imaging to extract MWF and T1 anatomical scans for image registration and segmentation. We investigated the inter-regional distributions of gene expression data from the Allen Human Brain Atlas in conjunction with inter-regional MWF distribution patterns. Pearson correlations were used to identify genes with expression profiles mirroring MWF. The Single Cell Type Atlas from the Human Protein Atlas was leveraged to classify genes into gene sets with high cell type specificity, and a control gene set with low cell type specificity. Then, we compared the Pearson correlation coefficients for each gene set to determine if cell type-specific gene expression signatures correlate with MWF. Pearson correlation coefficients between MWF and gene expression for oligodendrocytes and adipocytes were significantly higher than for the control gene set, whereas correlations between MWF and inhibitory/excitatory neurons were significantly lower. Our approach in integrating transcriptomics with neuroimaging measures supports an emerging technique for understanding and validating MRI-derived markers such as MWF.

Functional magnetic resonance imaging alternations in suicide attempts individuals and their association with gene expression

BACKGROUND

Functional Magnetic Resonance Imaging (fMRI) has shown brain activity alterations in individuals with a history of attempted suicide (SA) who are diagnosed with depression disorder (DD) or bipolar disorder (BD). However, patterns of spontaneous brain activity and their genetic correlations need further investigation.

METHODS

A voxel-based meta-analysis of 19 studies including 26 datasets, involving 742 patients with a history of SA and 978 controls (both nonsuicidal patients and healthy controls) was conducted. We examined fMRI changes in SA patients and analyzed the association between these changes and gene expression profiles using data from the Allen Human Brain Atlas by partial least squares regression analysis.

RESULTS

SA patients demonstrated increased spontaneous brain activity in several brain regions including the bilateral inferior temporal gyrus, hippocampus, fusiform gyrus, and right insula, and decreased activity in areas like the bilateral paracentral lobule and inferior frontal gyrus. Additionally, 5,077 genes were identified, exhibiting expression patterns associated with SA-related fMRI alterations. Functional enrichment analyses demonstrated that these SA-related genes were enriched for biological functions including glutamatergic synapse and mitochondrial structure. Concurrently, specific expression analyses showed that these genes were specifically expressed in the brain tissue, in neurons cells, and during early developmental periods.

CONCLUSION

Our findings suggest a neurobiological basis for fMRI abnormalities in SA patients with DD or BD, potentially guiding future genetic and therapeutic research.

Morphological Brain Networks of White Matter: Mapping, Evaluation, Characterization, and Application

Although white matter (WM) accounts for nearly half of adult brain, its wiring diagram is largely unknown. Here, an approach is developed to construct WM networks by estimating interregional morphological similarity based on structural magnetic resonance imaging. It is found that morphological WM networks showed nontrivial topology, presented good-to-excellent test-retest reliability, accounted for phenotypic interindividual differences in cognition, and are under genetic control. Through integration with multimodal and multiscale data, it is further showed that morphological WM networks are able to predict the patterns of hamodynamic coherence, metabolic synchronization, gene co-expression, and chemoarchitectonic covariance, and associated with structural connectivity. Moreover, the prediction followed WM functional connectomic hierarchy for the hamodynamic coherence, is related to genes enriched in the forebrain neuron development and differentiation for the gene co-expression, and is associated with serotonergic system-related receptors and transporters for the chemoarchitectonic covariance. Finally, applying this approach to multiple sclerosis and neuromyelitis optica spectrum disorders, it is found that both diseases exhibited morphological dysconnectivity, which are correlated with clinical variables of patients and are able to diagnose and differentiate the diseases. Altogether, these findings indicate that morphological WM networks provide a reliable and biologically meaningful means to explore WM architecture in health and disease.

Transcriptomics : Approaches to Quantifying Gene Expression and Their Application to Studying the Human Brain

To date, the field of transcriptomics has been characterized by rapid methods development and technological advancement, with new technologies continuously rendering older ones obsolete.This chapter traces the evolution of approaches to quantifying gene expression and provides an overall view of the current state of the field of transcriptomics, its applications to the study of the human brain, and its place in the broader emerging multiomics landscape.

Excitation/Inhibition balance relates to cognitive function and gene expression in temporal lobe epilepsy: a high density EEG assessment with aperiodic exponent

Patients with epilepsy are characterized by a dysregulation of excitation/inhibition balance (E/I). The assessment of E/I may inform clinicians during the diagnosis and therapy management, even though it is rarely performed. An accessible measure of the E/I of the brain represents a clinically relevant feature. Here, we exploited the exponent of the aperiodic component of the power spectrum of the electroencephalography (EEG) signal, as a non-invasive and cost-effective proxy of the E/I balance. We recorded resting-state activity with high-density EEG from 67 patients with temporal lobe epilepsy and 35 controls. We extracted the exponent of the aperiodic fit of the power spectrum from source-reconstructed EEG and tested differences between patients with epilepsy and controls. Spearman's correlation was performed between the exponent and clinical variables (age of onset, epilepsy duration and neuropsychology) and cortical expression of epilepsy-related genes derived from the Allen Human Brain Atlas. Patients with temporal lobe epilepsy showed a significantly larger exponent, corresponding to inhibition-directed E/I balance, in bilateral frontal and temporal regions. Lower E/I in the left entorhinal and bilateral dorsolateral prefrontal cortices corresponded to a lower performance of short-term verbal memory. Limited to patients with temporal lobe epilepsy, we detected a significant correlation between the exponent and the cortical expression of GABRA1, GRIN2A, GABRD, GABRG2, KCNA2 and PDYN genes. EEG aperiodic exponent maps the E/I balance non-invasively in patients with epilepsy and reveals a close relationship between altered E/I patterns, cognition and genetics.

MR Imaging Findings in Anti-Leucine-Rich Glioma Inactivated Protein 1 Encephalitis: A Systematic Review and Meta-analysis

BACKGROUND

Antibodies against leucine-rich glioma inactivated protein 1 (LGI1) constitute a common form of autoimmune encephalitis. On MR imaging, it may show T2 FLAIR hyperintensities of the medial temporal lobe (T2 FLAIR-MTL), involve the basal ganglia, or be unremarkable.

PURPOSE

We performed a systematic review and meta-analysis to obtain prevalence estimates of abnormal findings on MR imaging in anti-LGI1 encephalitis. A human brain map of the LGI1 microarray gene expression was derived from the Allen Human Brain Atlas.

DATA SOURCES

PubMed and Web of Science were searched with the terms "LGI1" and "encephalitis" from inception to April 7, 2022.

STUDY SELECTION

Thirty-one research publications, encompassing case series and retrospective cohort and case-control studies, with >10 patients with anti-LGI1 encephalitis and MR imaging data were included.

DATA ANALYSIS

Pooled prevalence estimates were calculated using Freeman-Tukey double-arcsine transformation. Meta-analysis used DerSimonian and Laird random effects models.

DATA SYNTHESIS

Of 1318 patients in 30 studies, T2 FLAIR-MTL hyperintensities were present in 54% (95% CI, 0.48-0.60; I2 = 76%). Of 394 patients in 13 studies, 27% showed bilateral (95% CI, 0.19-0.36; I2 = 71%) and 24% unilateral T2 FLAIR-MTL abnormalities (95% CI, 0.17-0.32; I2 = 61%). Of 612 patients in 15 studies, basal ganglia abnormalities were present in 10% (95% CI, 0.06-0.15; I2 = 67%). LGI1 expression was highest in the amygdala, hippocampus, and caudate nucleus.

LIMITATIONS

Only part of the spectrum of MR imaging abnormalities in anti-LGI1 encephalitis could be included in a meta-analysis. MR imaging findings were not the main outcomes in most studies, limiting available information. I2 values ranged from 62% to 76%, representing moderate-to-large heterogeneity.

CONCLUSIONS

T2 FLAIR-MTL hyperintensities were present in around one-half of patients with anti-LGI1. The prevalence of unilateral and bilateral presentations was similar, suggesting unilaterality should raise the suspicion of this disease in the appropriate clinical context. Around 10% of patients showed basal ganglia abnormalities, indicating that special attention should be given to this region. LGI1 regional expression coincided with the most frequently reported abnormal findings on MR imaging. Regional specificity might be partially determined by expression levels of the target protein.

Genes associated with cortical thickness alterations in behavioral addiction

Behavioral addiction (BA) is a conceptually new addictive phenotype characterized by compulsive reward-seeking behaviors despite adverse consequences. Currently, its underlying neurogenetic mechanism remains unclear. Here, this study aimed to investigate the association between cortical thickness (CTh) and genetic phenotypes in BA. We conducted a systematic search in five databases and extracted gene expression data from the Allen Human Brain Atlas. Meta-analysis of 10 studies (343 addicted individuals and 355 controls) revealed that the BA group showed thinner CTh in the precuneus, postcentral gyrus, orbital-frontal cortex, and dorsolateral prefrontal cortex (P < 0.005). Meta-regression showed that the CTh in the precuneus and postcentral gyrus were negatively associated with the addiction severity (P < 0.0005). More importantly, the CTh phenotype of BA was spatially correlated with the expression of 12 genes (false discovery rate [FDR] < 0.05), and the dopamine D2 receptor had the highest correlation (rho = 0.55). Gene enrichment analysis further revealed that the 12 genes were involved in the biological processes of behavior regulation and response to stimulus (FDR < 0.05). In conclusion, our findings demonstrated the thinner CTh in cognitive control-related brain areas in BA, which could be associated with the expression of genes involving dopamine metabolism and behavior regulation.

Association of Gene Expression and Tremor Network Structure

BACKGROUND

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

OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Associations between accelerometer-measured circadian rest-activity rhythm, brain structural and genetic mechanisms, and dementia

AIM

Knowledge of how circadian rhythm influences brain health remains limited. We aimed to investigate the associations of accelerometer-measured circadian rest-activity rhythm (CRAR) with incident dementia, cognitive dysfunction, and structural brain abnormalities in the general population and underlying biological mechanisms.

METHODS

Fifty-seven thousand five hundred and two participants aged over 60 years with accelerometer data were included to investigate the association of CRAR with incidental dementia. Non-parametric CRAR parameters were utilized, including activity level during active periods of the day (M10), activity level during rest periods of the day (L5), and the relative difference between the M10 and L5 (relative amplitude, RA). Associations of CRAR with cognitive dysfunction and brain structure were studied in a subset of participants. Neuroimaging-transcriptomics analysis was utilized to identify the underlying molecular mechanisms.

RESULTS

Over 6.86 (4.94-8.78) years of follow-up, 494 participants developed dementia. The risk of incident dementia was associated with decreasing M10 (hazard ratio [HR] 1.45; 95% conference interval [CI], 1.28-1.64) and RA (HR 1.37; 95% CI, 1.28-1.64), increasing L5 (HR 1.14, 95% CI 1.07-1.21) and advanced L5 onset time (HR 1.12; 95% CI, 1.02-1.23). The detrimental associations were exacerbated by APOE ε4 status and age (>65 years). Decreased RA was associated with lower processing speed (Beta -0.04; SE 0.011), predominantly mediated by abnormalities in subcortical regions and white matter microstructure. The genes underlying CRAR-related brain regional structure variation were enriched for synaptic function.

CONCLUSIONS

Our study underscores the potential of intervention targeting at maintaining a healthy CRAR pattern to prevent dementia risk.

State-specific Regulation of Electrical Stimulation in the Intralaminar Thalamus of Macaque Monkeys: Network and Transcriptional Insights into Arousal

Long-range thalamocortical communication is central to anesthesia-induced loss of consciousness and its reversal. However, isolating the specific neural networks connecting thalamic nuclei with various cortical regions for state-specific anesthesia regulation is challenging, with the biological underpinnings still largely unknown. Here, simultaneous electroencephalogram-fuctional magnetic resonance imaging (EEG-fMRI) and deep brain stimulation are applied to the intralaminar thalamus in macaques under finely-tuned propofol anesthesia. This approach led to the identification of an intralaminar-driven network responsible for rapid arousal during slow-wave oscillations. A network-based RNA-sequencing analysis is conducted of region-, layer-, and cell-specific gene expression data from independent transcriptomic atlases and identifies 2489 genes preferentially expressed within this arousal network, notably enriched in potassium channels and excitatory, parvalbumin-expressing neurons, and oligodendrocytes. Comparison with human RNA-sequencing data highlights conserved molecular and cellular architectures that enable the matching of homologous genes, protein interactions, and cell types across primates, providing novel insight into network-focused transcriptional signatures of arousal.

Spatial and temporal pattern of structure-function coupling of human brain connectome with development

Brain structural circuitry shapes a richly patterned functional synchronization, supporting for complex cognitive and behavioural abilities. However, how coupling of structural connectome (SC) and functional connectome (FC) develops and its relationships with cognitive functions and transcriptomic architecture remain unclear. We used multimodal magnetic resonance imaging data from 439 participants aged 5.7-21.9 years to predict functional connectivity by incorporating intracortical and extracortical structural connectivity, characterizing SC-FC coupling. Our findings revealed that SC-FC coupling was strongest in the visual and somatomotor networks, consistent with evolutionary expansion, myelin content, and functional principal gradient. As development progressed, SC-FC coupling exhibited heterogeneous alterations dominated by an increase in cortical regions, broadly distributed across the somatomotor, frontoparietal, dorsal attention, and default mode networks. Moreover, we discovered that SC-FC coupling significantly predicted individual variability in general intelligence, mainly influencing frontoparietal and default mode networks. Finally, our results demonstrated that the heterogeneous development of SC-FC coupling is positively associated with genes in oligodendrocyte-related pathways and negatively associated with astrocyte-related genes. This study offers insight into the maturational principles of SC-FC coupling in typical development.

Mapping Brain Synergy Dysfunction in Schizophrenia: Understanding Individual Differences and Underlying Molecular Mechanisms

To elucidate the brain-wide information interactions that vary and contribute to individual differences in schizophrenia (SCZ), an information-resolved method is employed to construct individual synergistic and redundant interaction matrices based on regional pairwise BOLD time-series from 538 SCZ and 540 normal controls (NC). This analysis reveals a stable pattern of regionally-specific synergy dysfunction in SCZ. Furthermore, a hierarchical Bayesian model is applied to deconstruct the patterns of whole-brain synergy dysfunction into three latent factors that explain symptom heterogeneity in SCZ. Factor 1 exhibits a significant positive correlation with Positive and Negative Syndrome Scale (PANSS) positive scores, while factor 3 demonstrates significant negative correlations with PANSS negative and general scores. By integrating the neuroimaging data with normative gene expression information, this study identifies that each of these three factors corresponded to a subset of the SCZ risk gene set. Finally, by combining data from NeuroSynth and open molecular imaging sources, along with a spatially heterogeneous mean-field model, this study delineates three SCZ synergy factors corresponding to distinct symptom profiles and implicating unique cognitive, neurodynamic, and neurobiological mechanisms.

Comparative Analysis of Human-Chimpanzee Divergence in Brain Connectivity and its Genetic Correlates

Chimpanzees (Pan troglodytes) are humans' closest living relatives, making them the most directly relevant comparison point for understanding human brain evolution. Zeroing in on the differences in brain connectivity between humans and chimpanzees can provide key insights into the specific evolutionary changes that might have occured along the human lineage. However, conducting comparisons of brain connectivity between humans and chimpanzees remains challenging, as cross-species brain atlases established within the same framework are currently lacking. Without the availability of cross-species brain atlases, the region-wise connectivity patterns between humans and chimpanzees cannot be directly compared. To address this gap, we built the first Chimpanzee Brainnetome Atlas (ChimpBNA) by following a well-established connectivity-based parcellation framework. Leveraging this new resource, we found substantial divergence in connectivity patterns across most association cortices, notably in the lateral temporal and dorsolateral prefrontal cortex between the two species. Intriguingly, these patterns significantly deviate from the patterns of cortical expansion observed in humans compared to chimpanzees. Additionally, we identified regions displaying connectional asymmetries that differed between species, likely resulting from evolutionary divergence. Genes associated with these divergent connectivities were found to be enriched in cell types crucial for cortical projection circuits and synapse formation. These genes exhibited more pronounced differences in expression patterns in regions with higher connectivity divergence, suggesting a potential foundation for brain connectivity evolution. Therefore, our study not only provides a fine-scale brain atlas of chimpanzees but also highlights the connectivity divergence between humans and chimpanzees in a more rigorous and comparative manner and suggests potential genetic correlates for the observed divergence in brain connectivity patterns between the two species. This can help us better understand the origins and development of uniquely human cognitive capabilities.

Quantitative susceptibility mapping in the brain reflects spatial expression of genes involved in iron homeostasis and myelination

Quantitative susceptibility mapping (QSM) is an MRI modality used to non-invasively measure iron content in the brain. Iron exhibits a specific anatomically varying pattern of accumulation in the brain across individuals. The highest regions of accumulation are the deep grey nuclei, where iron is stored in paramagnetic molecule ferritin. This form of iron is considered to be what largely contributes to the signal measured by QSM in the deep grey nuclei. It is also known that QSM is affected by diamagnetic myelin contents. Here, we investigate spatial gene expression of iron and myelin related genes, as measured by the Allen Human Brain Atlas, in relation to QSM images of age-matched subjects. We performed multiple linear regressions between gene expression and the average QSM signal within 34 distinct deep grey nuclei regions. Our results show a positive correlation (p < .05, corrected) between expression of ferritin and the QSM signal in deep grey nuclei regions. We repeated the analysis for other genes that encode proteins thought to be involved in the transport and storage of iron in the brain, as well as myelination. In addition to ferritin, our findings demonstrate a positive correlation (p < .05, corrected) between the expression of ferroportin, transferrin, divalent metal transporter 1, several gene markers of myelinating oligodendrocytes, and the QSM signal in deep grey nuclei regions. Our results suggest that the QSM signal reflects both the storage and active transport of iron in the deep grey nuclei regions of the brain.

BrainTACO: an explorable multi-scale multi-modal brain transcriptomic and connectivity data resource

Exploring the relationships between genes and brain circuitry can be accelerated by joint analysis of heterogeneous datasets from 3D imaging data, anatomical data, as well as brain networks at varying scales, resolutions, and modalities. Generating an integrated view, beyond the individual resources' original purpose, requires the fusion of these data to a common space, and a visualization that bridges the gap across scales. However, despite ever expanding datasets, few platforms for integration and exploration of this heterogeneous data exist. To this end, we present the BrainTACO (Brain Transcriptomic And Connectivity Data) resource, a selection of heterogeneous, and multi-scale neurobiological data spatially mapped onto a common, hierarchical reference space, combined via a holistic data integration scheme. To access BrainTACO, we extended BrainTrawler, a web-based visual analytics framework for spatial neurobiological data, with comparative visualizations of multiple resources. This enables gene expression dissection of brain networks with, to the best of our knowledge, an unprecedented coverage and allows for the identification of potential genetic drivers of connectivity in both mice and humans that may contribute to the discovery of dysconnectivity phenotypes. Hence, BrainTACO reduces the need for time-consuming manual data aggregation often required for computational analyses in script-based toolboxes, and supports neuroscientists by directly leveraging the data instead of preparing it.

Neuroimaging epicenters as potential sites of onset of the neuroanatomical pathology in schizophrenia

Schizophrenia lacks a clear definition at the neuroanatomical level, capturing the sites of origin and progress of this disorder. Using a network-theory approach called epicenter mapping on cross-sectional magnetic resonance imaging from 1124 individuals with schizophrenia, we identified the most likely "source of origin" of the structural pathology. Our results suggest that the Broca's area and adjacent frontoinsular cortex may be the epicenters of neuroanatomical pathophysiology in schizophrenia. These epicenters can predict an individual's response to treatment for psychosis. In addition, cross-diagnostic similarities based on epicenter mapping over of 4000 individuals diagnosed with neurological, neurodevelopmental, or psychiatric disorders appear to be limited. When present, these similarities are restricted to bipolar disorder, major depressive disorder, and obsessive-compulsive disorder. We provide a comprehensive framework linking schizophrenia-specific epicenters to multiple levels of neurobiology, including cognitive processes, neurotransmitter receptors and transporters, and human brain gene expression. Epicenter mapping may be a reliable tool for identifying the potential onset sites of neural pathophysiology in schizophrenia.

Disrupted topological properties of structural brain networks present a glutamatergic neuropathophysiology in people with narcolepsy

STUDY OBJECTIVES

Growing evidences have documented various abnormalities of the white matter bundles in people with narcolepsy. We sought to evaluate topological properties of brain structural networks, and their association with symptoms and neuropathophysiological features in people with narcolepsy.

METHODS

Diffusion tensor imaging was conducted for people with narcolepsy (n = 30) and matched healthy controls as well as symptoms assessment. Structural connectivity for each participant was generated to analyze global and regional topological properties and their correlations with narcoleptic features. Further 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.

RESULTS

A wide and dramatic decrease in structural connectivities was observed in people with narcolepsy, with descending network degree and global efficiency. These metrics were not only correlated with sleep latency and awakening features, but also reflected alterations of sleep macrostructure in people with narcolepsy. Network-based statistics identified a small hyperenhanced subnetwork of cingulate gyrus that was closely related to rapid eye movement sleep behavior disorder (RBD) in narcolepsy. Further imaging genetics analysis suggested glutamatergic signatures were responsible for the preferential vulnerability of connectivity alterations in people with narcolepsy, while additional PET/SPECT data verified that structural alteration was significantly correlated with metabotropic glutamate receptor 5 (mGlutR5) and N-methyl-D-aspartate receptor (NMDA).

CONCLUSIONS

People with narcolepsy endured a remarkable decrease in the structural architecture, which was not only closely related to narcolepsy symptoms but also glutamatergic signatures.

Atypical dynamic neural configuration in autism spectrum disorder and its relationship to gene expression profiles

Although it is well recognized that autism spectrum disorder (ASD) is associated with atypical dynamic functional connectivity patterns, the dynamic changes in brain intrinsic activity over each time point and the potential molecular mechanisms associated with atypical dynamic temporal characteristics in ASD remain unclear. Here, we employed the Hidden Markov Model (HMM) to explore the atypical neural configuration at every scanning time point in ASD, based on resting-state functional magnetic resonance imaging (rs-fMRI) data from the Autism Brain Imaging Data Exchange. Subsequently, partial least squares regression and pathway enrichment analysis were employed to explore the potential molecular mechanism associated with atypical neural dynamics in ASD. 8 HMM states were inferred from rs-fMRI data. Compared to typically developing, individuals on the autism spectrum showed atypical state-specific temporal characteristics, including number of states and occurrences, mean life time and transition probability between states. Moreover, these atypical temporal characteristics could predict communication difficulties of ASD, and states assoicated with negative activation in default mode network and frontoparietal network, and positive activation in somatomotor network, ventral attention network, and limbic network, had higher predictive contribution. Furthermore, a total of 321 genes was revealed to be significantly associated with atypical dynamic brain states of ASD, and these genes are mainly enriched in neurodevelopmental pathways. Our study provides new insights into characterizing the atypical neural dynamics from a moment-to-moment perspective, and indicates a linkage between atypical neural configuration and gene expression in ASD.

Repairing the in situ hybridization missing data in the hippocampus region by using a 3D residual U-Net model

The hippocampus is a critical brain region. Transcriptome data provides valuable insights into the structure and function of the hippocampus at the gene level. However, transcriptome data is often incomplete. To address this issue, we use the convolutional neural network model to repair the missing voxels in the hippocampus region, based on Allen institute coronal slices in situ hybridization (ISH) dataset. Moreover, we analyze the gene expression correlation between coronal and sagittal dataset in the hippocampus region. The results demonstrated that the trend of gene expression correlation between the coronal and sagittal datasets remained consistent following the repair of missing data in the coronal ISH dataset. In the last, we use repaired ISH dataset to identify novel genes specific to hippocampal subregions. Our findings demonstrate the accuracy and effectiveness of using deep learning method to repair ISH missing data. After being repaired, ISH has the potential to improve our comprehension of the hippocampus's structure and function.

Association between hearing ability and cortical morphology in the elderly: multiparametric mapping, cognitive relevance, and neurobiological underpinnings

BACKGROUND

Hearing impairment is a common condition in the elderly. However, a comprehensive understanding of its neural correlates is still lacking.

METHODS

We recruited 284 elderly adults who underwent structural MRI, magnetic resonance spectroscopy, audiometry, and cognitive assessments. Individual hearing abilities indexed by pure tone average (PTA) were correlated with multiple structural MRI-derived cortical morphological indices. For regions showing significant correlations, mediation analyses were performed to examine their role in the relationship between hearing ability and cognitive function. Finally, the correlation maps between hearing ability and cortical morphology were linked with publicly available connectomic gradient, transcriptomic, and neurotransmitter maps.

FINDINGS

Poorer hearing was related to cortical thickness (CT) reductions in widespread regions and gyrification index (GI) reductions in the right Area 52 and Insular Granular Complex. The GI in the right Area 52 mediated the relationship between hearing ability and executive function. This mediating effect was further modulated by glutamate and N-acetylaspartate levels in the right auditory region. The PTA-CT correlation map followed microstructural connectomic hierarchy, were related to genes involved in certain biological processes (e.g., glutamate metabolic process), cell types (e.g., excitatory neurons and astrocytes), and developmental stages (i.e., childhood to young adulthood), and covaried with dopamine receptor 1, dopamine transporter, and fluorodopa. The PTA-GI correlation map was related to 5-hydroxytryptamine receptor 2a.

INTERPRETATION

Poorer hearing is associated with cortical thinning and folding reductions, which may be engaged in the relationship between hearing impairment and cognitive decline in the elderly and have different neurobiological substrates.

FUNDING

See the Acknowledgements section.

Data-driven decomposition and staging of flortaucipir uptake in Alzheimer's disease

INTRODUCTION

Previous approaches pursuing in vivo staging of tau pathology in Alzheimer's disease (AD) have typically relied on neuropathologically defined criteria. In using predefined systems, these studies may miss spatial deposition patterns which are informative of disease progression.

METHODS

We selected discovery (n = 418) and replication (n = 132) cohorts with flortaucipir imaging. Non-negative matrix factorization (NMF) was applied to learn tau covariance patterns and develop a tau staging system. Flortaucipir components were also validated by comparison with amyloid burden, gray matter loss, and the expression of AD-related genes.

RESULTS

We found eight flortaucipir covariance patterns which were reproducible and overlapped with relevant gene expression maps. Tau stages were associated with AD severity as indexed by dementia status and neuropsychological performance. Comparisons of flortaucipir uptake with amyloid and atrophy also supported our model of tau progression.

DISCUSSION

Data-driven decomposition of flortaucipir uptake provides a novel framework for tau staging which complements existing systems.

HIGHLIGHTS

NMF reveals patterns of tau deposition in AD. Data-driven staging of flortaucipir tracks AD severity. Learned flortaucipir patterns overlap with AD-related gene expression.

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

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

Morphometric brain organization across the human lifespan reveals increased dispersion linked to cognitive performance

The human brain is organized as segregation and integration units and follows complex developmental trajectories throughout life. The cortical manifold provides a new means of studying the brain's organization in a multidimensional connectivity gradient space. However, how the brain's morphometric organization changes across the human lifespan remains unclear. Here, leveraging structural magnetic resonance imaging scans from 1,790 healthy individuals aged 8 to 89 years, we investigated age-related global, within- and between-network dispersions to reveal the segregation and integration of brain networks from 3D manifolds based on morphometric similarity network (MSN), combining multiple features conceptualized as a "fingerprint" of an individual's brain. Developmental trajectories of global dispersion unfolded along patterns of molecular brain organization, such as acetylcholine receptor. Communities were increasingly dispersed with age, reflecting more disassortative morphometric similarity profiles within a community. Increasing within-network dispersion of primary motor and association cortices mediated the influence of age on the cognitive flexibility of executive functions. We also found that the secondary sensory cortices were decreasingly dispersed with the rest of the cortices during aging, possibly indicating a shift of secondary sensory cortices across the human lifespan from an extreme to a more central position in 3D manifolds. Together, our results reveal the age-related segregation and integration of MSN from the perspective of a multidimensional gradient space, providing new insights into lifespan changes in multiple morphometric features of the brain, as well as the influence of such changes on cognitive performance.

Functional Connectome Hierarchy Distortions in Female Nurses With Occupational Burnout and Its Gene Expression Signatures

BACKGROUND

Burnout has become a serious public health issue worldwide, particularly during the COVID-19 pandemic. Functional connectome impairments associated with occupational burnout were widely distributed, involving both low-level sensorimotor cortices and high-level association cortices.

PURPOSE

To investigate whether there are hierarchical perturbations in the functional connectomes and if these perturbations are potentially influenced by genetic factors in nurses who feel "burned out."

STUDY TYPE

Prospective, case control.

POPULATION

Thirty-three female nurses with occupational burnout (aged 27-40, 32.42 ± 3.37) and 32 matched nurses who were not feeling burned out (aged 27-42, 32.50 ± 4.21).

FIELD STRENGTH/SEQUENCE

3.0 T, gradient-echo echo-planar imaging sequence (GE-EPI).

ASSESSMENT

Gradient-based techniques were used to depict the perturbations in the multi-dimensional hierarchical structure of the macroscale connectome. Gene expression data were acquired from the Allen Human Brain Atlas.

STATISTICAL TESTS

Cortex-wide multivariate analyses were used for between-group differences in gradients as well as association analyses between the hierarchy distortions and the MBI score (FDR corrected). Partial least squares, spin test and bootstrapping were utilized together to select the gene sets (FDR corrected). Gene enrichment analyses (GO, KEGG and cell-type) were further performed. Significance level: P < 0.05.

RESULTS

There were significant gradient distortions, with strong between-group effects in the somatosensory network and moderate effects in the higher-order default-mode network, which were significantly correlated with the gene expression profiles (r = 0.3171). The most related genes were broadly involved in the cellular response to minerals, neuronal plasticity, and the circadian rhythm pathway (q value < 0.01). Significant enrichments were found in excitatory (r = 0.2588), inhibitory neurons (r = 0.2610), and astrocytes cells (r = 0.2633). Regions affected by burnout severity were mainly distributed in the association and visual cortices.

DATA CONCLUSION

By connecting in vivo imaging to genes, cell classes, and clinical data, this study provides a framework to understand functional impairments in occupational burnout and how the microscopic genetic architecture drive macroscopic distortions.

EVIDENCE LEVEL

1 TECHNICAL EFFICACY: Stage 2.

Unraveling the molecular relevance of brain phenotypes: A comparative analysis of null models and test statistics

Correlating transcriptional profiles with imaging-derived phenotypes has the potential to reveal possible molecular architectures associated with cognitive functions, brain development and disorders. Competitive null models built by resampling genes and self-contained null models built by spinning brain regions, along with varying test statistics, have been used to determine the significance of transcriptional associations. However, there has been no systematic evaluation of their performance in imaging transcriptomics analyses. Here, we evaluated the performance of eight different test statistics (mean, mean absolute value, mean squared value, max mean, median, Kolmogorov-Smirnov (KS), Weighted KS and the number of significant correlations) in both competitive null models and self-contained null models. Simulated brain maps (n = 1,000) and gene sets (n = 500) were used to calculate the probability of significance (Psig) for each statistical test. Our results suggested that competitive null models may result in false positive results driven by co-expression within gene sets. Furthermore, we demonstrated that the self-contained null models may fail to account for distribution characteristics (e.g., bimodality) of correlations between all available genes and brain phenotypes, leading to false positives. These two confounding factors interacted differently with test statistics, resulting in varying outcomes. Specifically, the sign-sensitive test statistics (i.e., mean, median, KS, Weighted KS) were influenced by co-expression bias in the competitive null models, while median and sign-insensitive test statistics were sensitive to the bimodality bias in the self-contained null models. Additionally, KS-based statistics produced conservative results in the self-contained null models, which increased the risk of false negatives. Comprehensive supplementary analyses with various configurations, including realistic scenarios, supported the results. These findings suggest utilizing sign-insensitive test statistics such as mean absolute value, max mean in the competitive null models and the mean as the test statistic for the self-contained null models. Additionally, adopting the confounder-matched (e.g., coexpression-matched) null models as an alternative to standard null models can be a viable strategy. Overall, the present study offers insights into the selection of statistical tests for imaging transcriptomics studies, highlighting areas for further investigation and refinement in the evaluation of novel and commonly used tests.

Convergent Neuroimaging and Molecular Signatures in Mild Cognitive Impairment and Alzheimer's Disease: A Data-Driven Meta-Analysis with N = 3,118

The current study aimed to evaluate the susceptibility to regional brain atrophy and its biological mechanism in Alzheimer's disease (AD). We conducted data-driven meta-analyses to combine 3,118 structural magnetic resonance images from three datasets to obtain robust atrophy patterns. Then we introduced a set of radiogenomic analyses to investigate the biological basis of the atrophy patterns in AD. Our results showed that the hippocampus and amygdala exhibit the most severe atrophy, followed by the temporal, frontal, and occipital lobes in mild cognitive impairment (MCI) and AD. The extent of atrophy in MCI was less severe than that in AD. A series of biological processes related to the glutamate signaling pathway, cellular stress response, and synapse structure and function were investigated through gene set enrichment analysis. Our study contributes to understanding the manifestations of atrophy and a deeper understanding of the pathophysiological processes that contribute to atrophy, providing new insight for further clinical research on AD.

Neural correlates of obesity across the lifespan

Associations between brain and obesity are bidirectional: changes in brain structure and function underpin over-eating, while chronic adiposity leads to brain atrophy. Investigating brain-obesity interactions across the lifespan can help better understand these relationships. This study explores the interaction between obesity and cortical morphometry in children, young adults, adults, and older adults. We also investigate the genetic, neurochemical, and cognitive correlates of the brain-obesity associations. Our findings reveal a pattern of lower cortical thickness in fronto-temporal brain regions associated with obesity across all age cohorts and varying age-dependent patterns in the remaining brain regions. In adults and older adults, obesity correlates with neurochemical changes and expression of inflammatory and mitochondrial genes. In children and older adults, adiposity is associated with modifications in brain regions involved in emotional and attentional processes. Thus, obesity might originate from cognitive changes during early adolescence, leading to neurodegeneration in later life through mitochondrial and inflammatory mechanisms.

Characterizing the distinct imaging phenotypes, clinical behavior, and genetic vulnerability of brain maturational subtypes in mood disorders

BACKGROUND

Mood disorders are characterized by great heterogeneity in clinical manifestation. Uncovering such heterogeneity using neuroimaging-based individual biomarkers, clinical behaviors, and genetic risks, might contribute to elucidating the etiology of these diseases and support precision medicine.

METHODS

We recruited 174 drug-naïve and drug-free patients with major depressive disorder and bipolar disorder, as well as 404 healthy controls. T1 MRI imaging data, clinical symptoms, and neurocognitive assessments, and genetics were obtained and analyzed. We applied regional gray matter volumes (GMV) and quantile normative modeling to create maturation curves, and then calculated individual deviations to identify subtypes within the patients using hierarchical clustering. We compared the between-subtype differences in GMV deviations, clinical behaviors, cell-specific transcriptomic associations, and polygenic risk scores. We also validated the GMV deviations based subtyping analysis in a replication cohort.

RESULTS

Two subtypes emerged: subtype 1, characterized by increased GMV deviations in the frontal cortex, cognitive impairment, a higher genetic risk for Alzheimer's disease, and transcriptionally associated with Alzheimer's disease pathways, oligodendrocytes, and endothelial cells; and subtype 2, displaying globally decreased GMV deviations, more severe depressive symptoms, increased genetic vulnerability to major depressive disorder and transcriptionally related to microglia and inhibitory neurons. The distinct patterns of GMV deviations in the frontal, cingulate, and primary motor cortices between subtypes were shown to be replicable.

CONCLUSIONS

Our current results provide vital links between MRI-derived phenotypes, spatial transcriptome, genetic vulnerability, and clinical manifestation, and uncover the heterogeneity of mood disorders in biological and behavioral terms.

Exploring Schizophrenia Classification Through Multimodal MRI and Deep Graph Neural Networks: Unveiling Brain Region-Specific Weight Discrepancies and Their Association With Cell-Type Specific Transcriptomic Features

BACKGROUND AND HYPOTHESIS

Schizophrenia (SZ) is a prevalent mental disorder that imposes significant health burdens. Diagnostic accuracy remains challenging due to clinical subjectivity. To address this issue, we explore magnetic resonance imaging (MRI) as a tool to enhance SZ diagnosis and provide objective references and biomarkers. Using deep learning with graph convolution, we represent MRI data as graphs, aligning with brain structure, and improving feature extraction, and classification. Integration of multiple modalities is expected to enhance classification.

STUDY DESIGN

Our study enrolled 683 SZ patients and 606 healthy controls from 7 hospitals, collecting structural MRI and functional MRI data. Both data types were represented as graphs, processed by 2 graph attention networks, and fused for classification. Grad-CAM with graph convolution ensured interpretability, and partial least squares analyzed gene expression in brain regions.

STUDY RESULTS

Our method excelled in the classification task, achieving 83.32% accuracy, 83.41% sensitivity, and 83.20% specificity in 10-fold cross-validation, surpassing traditional methods. And our multimodal approach outperformed unimodal methods. Grad-CAM identified potential brain biomarkers consistent with gene analysis and prior research.

CONCLUSIONS

Our study demonstrates the effectiveness of deep learning with graph attention networks, surpassing previous SZ diagnostic methods. Multimodal MRI's superiority over unimodal MRI confirms our initial hypothesis. Identifying potential brain biomarkers alongside gene biomarkers holds promise for advancing objective SZ diagnosis and research in SZ.

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.

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

AIMS

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Similar imaging changes and their relations to genetic profiles in bipolar disorder across different clinical stages

Bipolar disorder (BD) across different clinical stages may present shared and distinct changes in brain activity. We aimed to reveal the neuroimaging homogeneity and heterogeneity of BD and its relationship with clinical variables and genetic variations. In present study, we conducted fractional amplitude of low-frequency fluctuations (fALFF), functional connectivity (FC) and genetic neuroimaging association analyses with 32 depressed, 26 manic, 35 euthymic BD patients and 87 healthy controls (HCs). Significant differences were found in the bilateral pre/subgenual anterior cingulate cortex (ACC) across the four groups, and all bipolar patients exhibited decreased fALFF values in the ACC when compared to HCs. Furthermore, positive associations were significantly observed between fALFF values in the pre/subgenual ACC and participants' cognitive functioning. No significant changes were found in ACC-based FC. We identified fALFF-alteration-related genes in BD, with enrichment in biological progress including synaptic and ion transmission. Taken together, abnormal activity in ACC is a characteristic change associated with BD, regardless of specific mood stages, serving as a potential neuroimaging feature in BD patients. Our genetic neuroimaging association analysis highlights possible heterogeneity in biological processes that could be responsible for different clinical stages in BD.

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

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