Toward Best Practices for Imaging Transcriptomics of the Human Brain

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.

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.

Omics approaches open new horizons in major depressive disorder: from biomarkers to precision medicine

Major depressive disorder (MDD) is a recurrent episodic mood disorder that represents the third leading cause of disability worldwide. In MDD, several factors can simultaneously contribute to its development, which complicates its diagnosis. According to practical guidelines, antidepressants are the first-line treatment for moderate to severe major depressive episodes. Traditional treatment strategies often follow a one-size-fits-all approach, resulting in suboptimal outcomes for many patients who fail to experience a response or recovery and develop the so-called "therapy-resistant depression". The high biological and clinical inter-variability within patients and the lack of robust biomarkers hinder the finding of specific therapeutic targets, contributing to the high treatment failure rates. In this frame, precision medicine, a paradigm that tailors medical interventions to individual characteristics, would help allocate the most adequate and effective treatment for each patient while minimizing its side effects. In particular, multi-omic studies may unveil the intricate interplays between genetic predispositions and exposure to environmental factors through the study of epigenomics, transcriptomics, proteomics, metabolomics, gut microbiomics, and immunomics. The integration of the flow of multi-omic information into molecular pathways may produce better outcomes than the current psychopharmacological approach, which targets singular molecular factors mainly related to the monoamine systems, disregarding the complex network of our organism. The concept of system biomedicine involves the integration and analysis of enormous datasets generated with different technologies, creating a "patient fingerprint", which defines the underlying biological mechanisms of every patient. This review, centered on precision medicine, explores the integration of multi-omic approaches as clinical tools for prediction in MDD at a single-patient level. It investigates how combining the existing technologies used for diagnostic, stratification, prognostic, and treatment-response biomarkers discovery with artificial intelligence can improve the assessment and treatment of MDD.

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

AIMS

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Homotopic functional connectivity disruptions in schizophrenia and their associated gene expression

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

Exploring functional dysconnectivity in schizophrenia: alterations in eigenvector centrality mapping and insights into related genes from transcriptional profiles

Schizophrenia is a mental health disorder characterized by functional dysconnectivity. Eigenvector centrality mapping (ECM) has been employed to investigate alterations in functional connectivity in schizophrenia, yet the results lack consistency, and the genetic mechanisms underlying these changes remain unclear. In this study, whole-brain voxel-wise ECM analyses were conducted on resting-state functional magnetic resonance imaging data. A cohort of 91 patients with schizophrenia and 91 matched healthy controls were included during the discovery stage. Additionally, in the replication stage, 153 individuals with schizophrenia and 182 healthy individuals participated. Subsequently, a comprehensive analysis was performed using an independent transcriptional database derived from six postmortem healthy adult brains to explore potential genetic factors influencing the observed functional dysconnectivity, and to investigate the roles of identified genes in neural processes and pathways. The results revealed significant and reliable alterations in the ECM across multiple brain regions in schizophrenia. Specifically, there was a significant decrease in ECM in the bilateral superior and middle temporal gyrus, and an increase in the bilateral thalamus in both the discovery and replication stages. Furthermore, transcriptional analysis revealed 420 genes whose expression patterns were related to changes in ECM, and these genes were enriched mainly in biological processes associated with synaptic signaling and transmission. Together, this study enhances our knowledge of the neural processes and pathways involved in schizophrenia, shedding light on the genetic factors that may be linked to functional dysconnectivity in this disorder.

Open datasets and code for multi-scale relations on structure, function and neuro-genetics in the human brain

The human brain is an extremely complex network of structural and functional connections that operate at multiple spatial and temporal scales. Investigating the relationship between these multi-scale connections is critical to advancing our comprehension of brain function and disorders. However, accurately predicting structural connectivity from its functional counterpart remains a challenging pursuit. One of the major impediments is the lack of public repositories that integrate structural and functional networks at diverse resolutions, in conjunction with modular transcriptomic profiles, which are essential for comprehensive biological interpretation. To mitigate this limitation, our contribution encompasses the provision of an open-access dataset consisting of derivative matrices of functional and structural connectivity across multiple scales, accompanied by code that facilitates the investigation of their interrelations. We also provide additional resources focused on neuro-genetic associations of module-level network metrics, which present promising opportunities to further advance research in the field of network neuroscience, particularly concerning brain disorders.

Spatial transcriptomic patterns underlying amyloid-β and tau pathology are associated with cognitive dysfunction in Alzheimer's disease

Amyloid-β (Aβ) and tau proteins accumulate within distinct neuronal systems in Alzheimer's disease (AD). Although it is not clear why certain brain regions are more vulnerable to Aβ and tau pathologies than others, gene expression may play a role. We study the association between brain-wide gene expression profiles and regional vulnerability to Aβ (gene-to-Aβ associations) and tau (gene-to-tau associations) pathologies by leveraging two large independent AD cohorts. We identify AD susceptibility genes and gene modules in a gene co-expression network with expression profiles specifically related to regional vulnerability to Aβ and tau pathologies in AD. In addition, we identify distinct biochemical pathways associated with the gene-to-Aβ and the gene-to-tau associations. These findings may explain the discordance between regional Aβ and tau pathologies. Finally, we propose an analytic framework, linking the identified gene-to-pathology associations to cognitive dysfunction in AD at the individual level, suggesting potential clinical implication of the gene-to-pathology associations.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Regional homogeneity alterations in patients with functional constipation and their associations with gene expression profiles

Functional constipation, a highly prevalent functional gastrointestinal disorder, often accompanies by mental and psychological disorders. Previous neuroimaging studies have demonstrated brain functional and structural alterations in patients with functional constipation. However, little is known about whether and how regional homogeneity is altered in these patients. Moreover, the potential genetic mechanisms associated with these alterations remain largely unknown. The study included 73 patients with functional constipation and 68 healthy controls, and regional homogeneity comparison was conducted to identify the abnormal spontaneous brain activities in patients with functional constipation. Using Allen Human Brain Atlas, we further investigated gene expression profiles associated with regional homogeneity alterations in functional constipation patients with partial least squares regression analysis applied. Compared with healthy controls, functional constipation patients demonstrated significantly decreased regional homogeneity in both bilateral caudate nucleus, putamen, anterior insula, thalamus and right middle cingulate cortex, supplementary motor area, and increased regional homogeneity in the bilateral orbitofrontal cortex. Genes related to synaptic signaling, central nervous system development, fatty acid metabolism, and immunity were spatially correlated with abnormal regional homogeneity patterns. Our findings showed significant regional homogeneity alterations in functional constipation patients, and the changes may be caused by complex polygenetic and poly-pathway mechanisms, which provides a new perspective on functional constipation's pathophysiology.

Longitudinal development of the human white matter structural connectome and its association with brain transcriptomic and cellular architecture

From childhood to adolescence, the spatiotemporal development pattern of the human brain white matter connectome and its underlying transcriptomic and cellular mechanisms remain largely unknown. With a longitudinal diffusion MRI cohort of 604 participants, we map the developmental trajectory of the white matter connectome from global to regional levels and identify that most brain network properties followed a linear developmental trajectory. Importantly, connectome-transcriptomic analysis reveals that the spatial development pattern of white matter connectome is potentially regulated by the transcriptomic architecture, with positively correlated genes involve in ion transport- and development-related pathways expressed in excitatory and inhibitory neurons, and negatively correlated genes enriches in synapse- and development-related pathways expressed in astrocytes, inhibitory neurons and microglia. Additionally, the macroscale developmental pattern is also associated with myelin content and thicknesses of specific laminas. These findings offer insights into the underlying genetics and neural mechanisms of macroscale white matter connectome development from childhood to adolescence.

Towards a neuroimaging biomarker for predicting cognitive behavioural therapy outcomes in treatment-naive depression: Preliminary findings

Clear prognostic indicators of cognitive behavioural therapy (CBT) are lacking for depression. This study aims to identify a biomarker that predicts CBT outcomes in depression. We developed a machine learning algorithm to predict post-CBT Hamilton Depression Rating Scale (HAMD) using pre-CBT regional homogeneity (ReHo). We examined transcriptomic signatures of regions with CBT-related ReHo changes. Twenty-five patients completed CBT and had increased ReHo in the dorsolateral prefrontal cortex (DLPFC) following CBT. Pre-CBT ReHo in left DLPFC was shown to be a predictor of post-HAMD scores. We identified left DLPFC ReHo as a neuroimaging biomarker for therapeutic effects of CBT in depression.

Two human brain systems micro-structurally associated with obesity

The relationship between obesity and human brain structure is incompletely understood. Using diffusion-weighted MRI from ∼30,000 UK Biobank participants, we test the hypothesis that obesity (waist-to-hip ratio, WHR) is associated with regional differences in two micro-structural MRI metrics: isotropic volume fraction (ISOVF), an index of free water, and intra-cellular volume fraction (ICVF), an index of neurite density. We observed significant associations with obesity in two coupled but distinct brain systems: a prefrontal/temporal/striatal system associated with ISOVF and a medial temporal/occipital/striatal system associated with ICVF. The ISOVF~WHR system colocated with expression of genes enriched for innate immune functions, decreased glial density, and high mu opioid (MOR) and other neurotransmitter receptor density. Conversely, the ICVF~WHR system co-located with expression of genes enriched for G-protein coupled receptors and decreased density of MOR and other receptors. To test whether these distinct brain phenotypes might differ in terms of their underlying shared genetics or relationship to maps of the inflammatory marker C-reactive Protein (CRP), we estimated the genetic correlations between WHR and ISOVF (rg = 0.026, P = 0.36) and ICVF (rg = 0.112, P < 9×10-4) as well as comparing correlations between WHR maps and equivalent CRP maps for ISOVF and ICVF (P<0.05). These correlational results are consistent with a two-way mechanistic model whereby genetically determined differences in neurite density in the medial temporal system may contribute to obesity, whereas water content in the prefrontal system could reflect a consequence of obesity mediated by innate immune system activation.

Cortical structural differences following repeated ayahuasca use hold molecular signatures

Introduction

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

Methods

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

Results

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

Conclusion

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

Spatial transcriptomic patterns underlying regional vulnerability to amyloid-β and tau pathologies and their relationships to cognitive dysfunction in Alzheimer's disease

Amyloid-β (Aβ) and tau proteins accumulate within distinct neuronal systems in Alzheimer's disease (AD). Although it is not clear why certain brain regions are more vulnerable to Aβ and tau pathologies than others, gene expression may play a role. We studied the association between brain-wide gene expression profiles and regional vulnerability to Aβ (gene-to-Aβ associations) and tau (gene-to-tau associations) pathologies leveraging two large independent cohorts (n = 715) of participants along the AD continuum. We identified several AD susceptibility genes and gene modules in a gene co-expression network with expression profiles related to regional vulnerability to Aβ and tau pathologies in AD. In particular, we found that the positive APOE -to-tau association was only seen in the AD cohort, whereas patients with AD and frontotemporal dementia shared similar positive MAPT -to-tau association. Some AD candidate genes showed sex-dependent negative gene-to-Aβ and gene-to-tau associations. In addition, we identified distinct biochemical pathways associated with the gene-to-Aβ and the gene-to-tau associations. Finally, we proposed a novel analytic framework, linking the identified gene-to-pathology associations to cognitive dysfunction in AD at the individual level, suggesting potential clinical implication of the gene-to-pathology associations. Taken together, our study identified distinct gene expression profiles and biochemical pathways that may explain the discordance between regional Aβ and tau pathologies, and filled the gap between gene-to-pathology associations and cognitive dysfunction in individual AD patients that may ultimately help identify novel personalized pathogenetic biomarkers and therapeutic targets.

One Sentence Summary

We identified replicable cognition-related associations between regional gene expression profiles and selectively regional vulnerability to amyloid-β and tau pathologies in AD.

Mitochondrial function-associated genes underlie cortical atrophy in prodromal synucleinopathies

Isolated rapid eye movement sleep behaviour disorder (iRBD) is a sleep disorder characterized by the loss of rapid eye movement sleep muscle atonia and the appearance of abnormal movements and vocalizations during rapid eye movement sleep. It is a strong marker of incipient synucleinopathy such as dementia with Lewy bodies and Parkinson's disease. Patients with iRBD already show brain changes that are reminiscent of manifest synucleinopathies including brain atrophy. However, the mechanisms underlying the development of this atrophy remain poorly understood. In this study, we performed cutting-edge imaging transcriptomics and comprehensive spatial mapping analyses in a multicentric cohort of 171 polysomnography-confirmed iRBD patients [67.7 ± 6.6 (49-87) years; 83% men] and 238 healthy controls [66.6 ± 7.9 (41-88) years; 77% men] with T1-weighted MRI to investigate the gene expression and connectivity patterns associated with changes in cortical thickness and surface area in iRBD. Partial least squares regression was performed to identify the gene expression patterns underlying cortical changes in iRBD. Gene set enrichment analysis and virtual histology were then done to assess the biological processes, cellular components, human disease gene terms, and cell types enriched in these gene expression patterns. We then used structural and functional neighbourhood analyses to assess whether the atrophy patterns in iRBD were constrained by the brain's structural and functional connectome. Moreover, we used comprehensive spatial mapping analyses to assess the specific neurotransmitter systems, functional networks, cytoarchitectonic classes, and cognitive brain systems associated with cortical changes in iRBD. All comparisons were tested against null models that preserved spatial autocorrelation between brain regions and compared to Alzheimer's disease to assess the specificity of findings to synucleinopathies. We found that genes involved in mitochondrial function and macroautophagy were the strongest contributors to the cortical thinning occurring in iRBD. Moreover, we demonstrated that cortical thinning was constrained by the brain's structural and functional connectome and that it mapped onto specific networks involved in motor and planning functions. In contrast with cortical thickness, changes in cortical surface area were related to distinct genes, namely genes involved in the inflammatory response, and to different spatial mapping patterns. The gene expression and connectivity patterns associated with iRBD were all distinct from those observed in Alzheimer's disease. In summary, this study demonstrates that the development of brain atrophy in synucleinopathies is constrained by specific genes and networks.

Robust estimation of cortical similarity networks from brain MRI

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

The Cellular Underpinnings of the Human Cortical Connectome

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

A systematic review of neuroimaging epigenetic research: calling for an increased focus on development

Epigenetic mechanisms, such as DNA methylation (DNAm), have gained increasing attention as potential biomarkers and mechanisms underlying risk for neurodevelopmental, psychiatric and other brain-based disorders. Yet, surprisingly little is known about the extent to which DNAm is linked to individual differences in the brain itself, and how these associations may unfold across development - a time of life when many of these disorders emerge. Here, we systematically review evidence from the nascent field of Neuroimaging Epigenetics, combining structural or functional neuroimaging measures with DNAm, and the extent to which the developmental period (birth to adolescence) is represented in these studies. We identified 111 articles published between 2011-2021, out of which only a minority (21%) included samples under 18 years of age. Most studies were cross-sectional (85%), employed a candidate-gene approach (67%), and examined DNAm-brain associations in the context of health and behavioral outcomes (75%). Nearly half incorporated genetic data, and a fourth investigated environmental influences. Overall, studies support a link between peripheral DNAm and brain imaging measures, but there is little consistency in specific findings and it remains unclear whether DNAm markers present a cause, correlate or consequence of brain alterations. Overall, there is large heterogeneity in sample characteristics, peripheral tissue and brain outcome examined as well as the methods used. Sample sizes were generally low to moderate (median nall = 98, ndevelopmental = 80), and attempts at replication or meta-analysis were rare. Based on the strengths and weaknesses of existing studies, we propose three recommendations on how advance the field of Neuroimaging Epigenetics. We advocate for: (1) a greater focus on developmentally oriented research (i.e. pre-birth to adolescence); (2) the analysis of large, prospective, pediatric cohorts with repeated measures of DNAm and imaging to assess directionality; and (3) collaborative, interdisciplinary science to identify robust signals, triangulate findings and enhance translational potential.

Transcriptomics and magnetic resonance imaging in major psychiatric disorders

Major psychiatric disorders create a significant public health burden, and mental disorders such as major depressive disorder, bipolar disorder, and schizophrenia are major contributors to the national disease burden. The search for biomarkers has been a leading endeavor in the field of biological psychiatry in recent decades. And the application of cross-scale and multi-omics approaches combining genes and imaging in major psychiatric studies has facilitated the elucidation of gene-related pathogenesis and the exploration of potential biomarkers. In this article, we summarize the results of using combined transcriptomics and magnetic resonance imaging to understand structural and functional brain changes associated with major psychiatric disorders in the last decade, demonstrating the neurobiological mechanisms of genetically related structural and functional brain alterations in multiple directions, and providing new avenues for the development of quantifiable objective biomarkers, as well as clinical diagnostic and prognostic indicators.