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

Brain change trajectories in healthy adults correlate with Alzheimer's related genetic variation and memory decline across life

Throughout adulthood and ageing our brains undergo structural loss in an average pattern resembling faster atrophy in Alzheimer's disease (AD). Using a longitudinal adult lifespan sample (aged 30-89; 2-7 timepoints) and four polygenic scores for AD, we show that change in AD-sensitive brain features correlates with genetic AD-risk and memory decline in healthy adults. We first show genetic risk links with more brain loss than expected for age in early Braak regions, and find this extends beyond APOE genotype. Next, we run machine learning on AD-control data from the Alzheimer's Disease Neuroimaging Initiative using brain change trajectories conditioned on age, to identify AD-sensitive features and model their change in healthy adults. Genetic AD-risk linked with multivariate change across many AD-sensitive features, and we show most individuals over age ~50 are on an accelerated trajectory of brain loss in AD-sensitive regions. Finally, high genetic risk adults with elevated brain change showed more memory decline through adulthood, compared to high genetic risk adults with less brain change. Our findings suggest quantitative AD risk factors are detectable in healthy individuals, via a shared pattern of ageing- and AD-related neurodegeneration that occurs along a continuum and tracks memory decline through adulthood.

Transcriptional signatures of cortical structural changes in chronic insomnia disorder

Chronic insomnia disorder (CID) is a multidimensional disease that may influence various levels of brain organization, spanning the macroscopic structural connectome to microscopic gene expression. However, the connection between genomic variations and morphological alterations in CID remains unclear. Here, we investigated brain structural changes in CID patients at the whole-brain level and whether these link to transcriptional characteristics. Brain structural data from 104 CID patients and 102 matched healthy controls (HC) were acquired to examine cortical structural alterations using morphometric similarity (MS) analysis. Partial least squares (PLS) regression and transcriptome data from the Allen Human Brain Atlas were used to extract genomes related to MS changes. Gene-category enrichment analysis (GCEA) was used to identify potential molecular mechanisms behind the observed structural changes. We found that CID patients exhibited MS reductions in the parietal and limbic regions, along with enhancements in the temporal and frontal regions compared to HCs (pFDR < .05). Subsequently, PLS and GCEA revealed that these MS alterations were spatially correlated with a set of genes, especially those significantly correlated with excitatory and inhibitory neurons and chronic neuroinflammation. This neuroimaging-transcriptomic study bridges the gap between cortical structural changes and the molecular mechanisms in CID patients, providing novel insight into the pathophysiology of insomnia and targeted treatments.

Integrating amyloid and tau imaging with proteomics and genomics in Alzheimer's disease

Alzheimer's disease (AD) is the most common neurodegenerative disease and is characterized by the aggregation of β-amyloid (Aβ) and tau in the brain. Breakthroughs in disease-modifying treatments targeting Aβ bring new hope for the management of AD. But to effectively modify and someday even prevent AD, a better understanding is needed of the biological mechanisms that underlie and link Aβ and tau in AD. Developments of high-throughput omics, including genomics, proteomics, and transcriptomics, together with molecular imaging of Aβ and tau with positron emission tomography (PET), allow us to discover and understand the biological pathways that regulate the aggregation and spread of Aβ and tau in living humans. The field of integrated omics and PET studies of Aβ and tau in AD is growing rapidly. We here provide an update of this field, both in terms of biological insights and in terms of future clinical implications of integrated omics-molecular imaging studies.

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.

Deciphering the functional specialization of whole-brain spatiomolecular gradients in the adult brain

Cortical arealization arises during neurodevelopment from the confluence of molecular gradients representing patterned expression of morphogens and transcription factors. However, whether similar gradients are maintained in the adult brain remains unknown. Here, we uncover three axes of topographic variation in gene expression in the adult human brain that specifically capture previously identified rostral-caudal, dorsal-ventral, and medial-lateral axes of early developmental patterning. The interaction of these spatiomolecular gradients i) accurately reconstructs the position of brain tissue samples, ii) delineates known functional territories, and iii) can model the topographical variation of diverse cortical features. The spatiomolecular gradients are distinct from canonical cortical axes differentiating the primary sensory cortex from the association cortex, but radiate in parallel with the axes traversed by local field potentials along the cortex. We replicate all three molecular gradients in three independent human datasets as well as two nonhuman primate datasets and find that each gradient shows a distinct developmental trajectory across the lifespan. The gradients are composed of several well-known transcription factors (e.g., PAX6 and SIX3), and a small set of genes shared across gradients are strongly enriched for multiple diseases. Together, these results provide insight into the developmental sculpting of functionally distinct brain regions, governed by three robust transcriptomic axes embedded within brain parenchyma.

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.

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.

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.

The alteration of cortical microstructure similarity in drug-resistant epilepsy correlated with mTOR pathway genes

BACKGROUND

Drug-resistant epilepsy (DRE) is associated with distributed laminar disruptions due to cytoarchitectonic pathologies, which may be characterized by multimodal MRI approaches such as morphometric similarity networks (MSNs). However, the genetic and histological underpinning of MSN alterations in DRE remains poorly understood, hampering its clinical application.

METHODS

We enrolled 60 patients with DRE and 23 controls, acquiring T1 and diffusion spectrum imaging data with a 3.0T GE SIGNA Premier scanner. Morphometric similarity networks (MSNs) were constructed and analyzed to identify microstructure similarity differences between patients and controls. Subsequently, patient-specific MSN alteration patterns were associated with gene expression using the GAMBA tool, and layer-specific neuronal signature mapping were also applied. During these analyses, sex and age were adjusted as covariates and multiple comparisons corrections were applied when appropriate.

FINDINGS

We observed widespread MSN changes in patients with DRE and identified five distinct MSN alteration patterns. Major patterns presented pattern-specific associations with expressions of epilepsy-related genes, particularly involving the mTOR pathway. Histological analysis confirmed the presence of cortical microstructure changes in areas with MSN alterations and revealed cellular abnormalities matching the aforementioned genetic risks.

INTERPRETATION

Our findings highlight the potential of quantifying laminar-related microstructure integrity using MSN to uncover the cytoarchitectonic changes in the pathophysiology of DRE. This approach may facilitate the identification of genetic and histological underpinnings of MSN alterations in DRE, ultimately aiding in the development of targeted therapeutic strategies.

FUNDINGS

The National Natural Science Foundation of China, the Ministry of Science and Technology of the People's Republic of China, and the Beijing Municipal Health Commission.

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.

Case-control virtual histology elucidates cell types associated with cortical thickness differences in Alzheimer's disease

Many neuropsychiatric disorders are characterised by altered cortical thickness, but the cell types underlying these changes remain largely unknown. Virtual histology (VH) approaches map regional patterns of gene expression with regional patterns of MRI-derived phenotypes, such as cortical thickness, to identify cell types associated with case-control differences in those MRI measures. However, this method does not incorporate valuable information of case-control differences in cell type abundance. We developed a novel method, termed case-control virtual histology (CCVH), and applied it to Alzheimer's disease (AD) and dementia cohorts. Leveraging a multi-region gene expression dataset of AD cases (n = 40) and controls (n = 20), we quantified AD case-control differential expression of cell type-specific markers across 13 brain regions. We then correlated these expression effects with MRI-derived AD case-control cortical thickness differences across the same regions. Cell types with spatially concordant AD-related effects were identified through resampling marker correlation coefficients. Among regions thinner in AD, gene expression patterns identified by CCVH suggested fewer excitatory and inhibitory neurons, and greater proportions of astrocytes, microglia, oligodendrocytes, oligodendrocyte precursor cells, and endothelial cells in AD cases vs. controls. In contrast, original VH identified expression patterns suggesting that excitatory but not inhibitory neuron abundance was associated with thinner cortex in AD, despite the fact that both types of neurons are known to be lost in the disorder. Compared to original VH, cell types identified through CCVH are more likely to directly underlie cortical thickness differences in AD. Sensitivity analyses suggest our results are largely robust to specific analysis choices, like numbers of cell type-specific marker genes used and background gene sets used to construct null models. As more multi-region brain expression datasets become available, CCVH will be useful for identifying the cellular correlates of cortical thickness across neuropsychiatric illnesses.

Novel effects of Ras-MAPK pathogenic variants on the developing human brain and their link to gene expression and inhibition abilities

The RASopathies are genetic syndromes associated with pathogenic variants causing dysregulation of the Ras/mitogen-activated protein kinase (Ras-MAPK) pathway, essential for brain development, and increased risk for neurodevelopmental disorders. Yet, the effects of most pathogenic variants on the human brain are unknown. We examined: (1) How Ras-MAPK activating variants of PTPN11/SOS1 protein-coding genes affect brain anatomy. (2) The relationship between PTPN11 gene expression levels and brain anatomy, and (3) The relevance of subcortical anatomy to attention and memory skills affected in the RASopathies. We collected structural brain MRI and cognitive-behavioral data from 40 pre-pubertal children with Noonan syndrome (NS), caused by PTPN11 (n = 30) or SOS1 (n = 10) variants (age 8.53 ± 2.15, 25 females), and compared them to 40 age- and sex-matched typically developing controls (9.24 ± 1.62, 27 females). We identified widespread effects of NS on cortical and subcortical volumes and on determinants of cortical gray matter volume, surface area (SA), and cortical thickness (CT). In NS, we observed smaller volumes of bilateral striatum, precentral gyri, and primary visual area (d's < -0.8), and extensive effects on SA (d's > |0.8|) and CT (d's > |0.5|) relative to controls. Further, SA effects were associated with increasing PTPN11 gene expression, most prominently in the temporal lobe. Lastly, PTPN11 variants disrupted normative relationships between the striatum and inhibition functioning. We provide evidence for the effects of Ras-MAPK pathogenic variants on striatal and cortical anatomy as well as links between PTPN11 gene expression and cortical SA increases, and striatal volume and inhibition skills. These findings provide essential translational information on the Ras-MAPK pathway's effect on human brain development and function.

Cortical anatomical variations, gene expression profiles, and clinical phenotypes in patients with schizophrenia

BACKGROUND AND HYPOTHESIS

Schizophrenia (SZ) patients display significant structural brain abnormalities; nevertheless, the genetic mechanisms regulating cortical anatomical variations and their correlation with the disease phenotype are still ambiguous.

STUDY DESIGN

We characterized anatomical variation using a surface-based method derived from structural magnetic resonance imaging of patients with SZ and age- and sex-matched healthy controls (HCs). Partial least-squares regression was performed across cortex regions between anatomical variation and average transcriptional profiles of SZ risk genes and all qualified genes from the Allen Human Brain Atlas. The morphological features of each brain region were correlated to symptomology variables in patients with SZ using partial correlation analysis.

STUDY RESULTS

A total of 203 SZ and 201 HCs were included in the final analysis. We observed significant variation of 55 regions of cortical thickness, 23 regions of volume, 7 regions of area, and 55 regions of local gyrification index (LGI) between SZ and HC groups. Expression profiles of 4 SZ risk genes and 96 genes from all qualified genes showed a correlation to anatomical variability, however, after multiple comparisons, the correlations were no longer significant. LGI variability in multiple frontal subregions was associated with specific symptoms of SZ, whereas cognitive function involving attention/vigilance was linked to LGI variability across nine brain regions.

CONCLUSIONS

Cortical anatomical variation of patients with schizophrenia is associated with gene transcriptome profiles as well as clinical phenotypes.

Toward Best Practices for Imaging Transcriptomics of the Human Brain

Modern brainwide transcriptional atlases provide unprecedented opportunities for investigating the molecular correlates of brain organization, as quantified using noninvasive neuroimaging. However, integrating neuroimaging data with transcriptomic measures is not straightforward, and careful consideration is required to make valid inferences. In this article, we review recent work exploring how various methodological choices affect 3 main phases of imaging transcriptomic analyses, including 1) processing of transcriptional atlas data; 2) relating transcriptional measures to independently derived neuroimaging phenotypes; and 3) evaluating the functional implications of identified associations through gene enrichment analyses. Our aim is to facilitate the development of standardized and reproducible approaches for this rapidly growing field. We identify sources of methodological variability, key choices that can affect findings, and considerations for mitigating false positive and/or spurious results. Finally, we provide an overview of freely available open-source toolboxes implementing current best-practice procedures across all 3 analysis phases.

Novel effects of Ras-MAPK pathogenic variants on the developing human brain and their link to gene expression and inhibition abilities

The RASopathies are genetic syndromes associated with pathogenic variants causing dysregulation of the Ras/mitogen-activated protein kinase (Ras-MAPK) pathway, essential for brain development, and increased risk for neurodevelopmental disorders. Yet, the effects of most pathogenic variants on the human brain are unknown. We examined: 1. How Ras-MAPK activating variants of PTPN11 / SOS1 protein-coding genes affect brain anatomy. 2. The relationship between PTPN11 gene expression levels and brain anatomy, and 3. The relevance of subcortical anatomy to attention and memory skills affected in the RASopathies. We collected structural brain MRI and cognitive-behavioral data from 40 pre-pubertal children with Noonan syndrome (NS), caused by PTPN11 ( n  = 30) or SOS1 ( n  = 10) variants (age 8.53 ± 2.15, 25 females), and compared them to 40 age- and sex-matched typically developing controls (9.24 ± 1.62, 27 females). We identified widespread effects of NS on cortical and subcortical volumes and on determinants of cortical gray matter volume, surface area (SA) and cortical thickness (CT). In NS, we observed smaller volumes of bilateral striatum, precentral gyri, and primary visual area ( d 's<-0.8), and extensive effects on SA ( d 's>|0.8|) and CT ( d 's>|0.5|) relative to controls. Further, SA effects were associated with increasing PTPN11 gene expression, most prominently in the temporal lobe. Lastly, PTPN11 variants disrupted normative relationships between the striatum and inhibition functioning. We provide evidence for effects of Ras-MAPK pathogenic variants on striatal and cortical anatomy as well as links between PTPN11 gene expression and cortical SA increases, and striatal volume and inhibition skills. These findings provide essential translational information on the Ras-MAPK pathway's effect on human brain development and function.

Long-range functional connections mirror and link microarchitectural and cognitive hierarchies in the human brain

BACKGROUND

Higher-order cognition is hypothesized to be implemented via distributed cortical networks that are linked via long-range connections. However, it is unknown how computational advantages of long-range connections reflect cortical microstructure and microcircuitry.

METHODS

We investigated this question by (i) profiling long-range cortical connectivity using resting-state functional magnetic resonance imaging (MRI) and cortico-cortical geodesic distance mapping, (ii) assessing how long-range connections reflect local brain microarchitecture, and (iii) examining the microarchitectural similarity of regions connected through long-range connections.

RESULTS

Analysis of 2 independent datasets indicated that sensory/motor areas had more clustered short-range connections, while transmodal association systems hosted distributed, long-range connections. Meta-analytical decoding suggested that this topographical difference mirrored shifts in cognitive function, from perception/action towards emotional/social processing. Analysis of myelin-sensitive in vivo MRI as well as postmortem histology and transcriptomics datasets established that gradients in functional connectivity distance are paralleled by those present in cortical microarchitecture. Notably, long-range connections were found to link spatially remote regions of association cortex with an unexpectedly similar microarchitecture.

CONCLUSIONS

By mapping covarying topographies of long-range functional connections and cortical microcircuits, the current work provides insights into structure-function relations in human neocortex.

Coordinated cortical thickness alterations across six neurodevelopmental and psychiatric disorders

Neuropsychiatric disorders are increasingly conceptualized as overlapping spectra sharing multi-level neurobiological alterations. However, whether transdiagnostic cortical alterations covary in a biologically meaningful way is currently unknown. Here, we studied co-alteration networks across six neurodevelopmental and psychiatric disorders, reflecting pathological structural covariance. In 12,024 patients and 18,969 controls from the ENIGMA consortium, we observed that co-alteration patterns followed normative connectome organization and were anchored to prefrontal and temporal disease epicenters. Manifold learning revealed frontal-to-temporal and sensory/limbic-to-occipitoparietal transdiagnostic gradients, differentiating shared illness effects on cortical thickness along these axes. The principal gradient aligned with a normative cortical thickness covariance gradient and established a transcriptomic link to cortico-cerebello-thalamic circuits. Moreover, transdiagnostic gradients segregated functional networks involved in basic sensory, attentional/perceptual, and domain-general cognitive processes, and distinguished between regional cytoarchitectonic profiles. Together, our findings indicate that shared illness effects occur in a synchronized fashion and along multiple levels of hierarchical cortical organization.

Meta-connectomic analysis maps consistent, reproducible, and transcriptionally relevant functional connectome hubs in the human brain

Human brain connectomes include sets of densely connected hub regions. However, the consistency and reproducibility of functional connectome hubs have not been established to date and the genetic signatures underlying robust hubs remain unknown. Here, we conduct a worldwide harmonized meta-connectomic analysis by pooling resting-state functional MRI data of 5212 healthy young adults across 61 independent cohorts. We identify highly consistent and reproducible connectome hubs in heteromodal and unimodal regions both across cohorts and across individuals, with the greatest effects observed in lateral parietal cortex. These hubs show heterogeneous connectivity profiles and are critical for both intra- and inter-network communications. Using post-mortem transcriptome datasets, we show that as compared to non-hubs, connectome hubs have a spatiotemporally distinctive transcriptomic pattern dominated by genes involved in the neuropeptide signaling pathway, neurodevelopmental processes, and metabolic processes. These results highlight the robustness of macroscopic connectome hubs and their potential cellular and molecular underpinnings, which markedly furthers our understanding of how connectome hubs emerge in development, support complex cognition in health, and are involved in disease.

Structural network alterations in focal and generalized epilepsy assessed in a worldwide ENIGMA study follow axes of epilepsy risk gene expression

Epilepsy is associated with genetic risk factors and cortico-subcortical network alterations, but associations between neurobiological mechanisms and macroscale connectomics remain unclear. This multisite ENIGMA-Epilepsy study examined whole-brain structural covariance networks in patients with epilepsy and related findings to postmortem epilepsy risk gene expression patterns. Brain network analysis included 578 adults with temporal lobe epilepsy (TLE), 288 adults with idiopathic generalized epilepsy (IGE), and 1328 healthy controls from 18 centres worldwide. Graph theoretical analysis of structural covariance networks revealed increased clustering and path length in orbitofrontal and temporal regions in TLE, suggesting a shift towards network regularization. Conversely, people with IGE showed decreased clustering and path length in fronto-temporo-parietal cortices, indicating a random network configuration. Syndrome-specific topological alterations reflected expression patterns of risk genes for hippocampal sclerosis in TLE and for generalized epilepsy in IGE. These imaging-transcriptomic signatures could potentially guide diagnosis or tailor therapeutic approaches to specific epilepsy syndromes.

Null models in network neuroscience

Recent advances in imaging and tracing technology provide increasingly detailed reconstructions of brain connectomes. Concomitant analytic advances enable rigorous identification and quantification of functionally important features of brain network architecture. Null models are a flexible tool to statistically benchmark the presence or magnitude of features of interest, by selectively preserving specific architectural properties of brain networks while systematically randomizing others. Here we describe the logic, implementation and interpretation of null models of connectomes. We introduce randomization and generative approaches to constructing null networks, and outline a taxonomy of network methods for statistical inference. We highlight the spectrum of null models - from liberal models that control few network properties, to conservative models that recapitulate multiple properties of empirical networks - that allow us to operationalize and test detailed hypotheses about the structure and function of brain networks. We review emerging scenarios for the application of null models in network neuroscience, including for spatially embedded networks, annotated networks and correlation-derived networks. Finally, we consider the limits of null models, as well as outstanding questions for the field.

Imaging Transcriptomics of Brain Disorders

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