Overcoming false-positive gene-category enrichment in the analysis of spatially resolved transcriptomic brain atlas data

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.

A worldwide harmonized meta-connectomic analysis of 5212 healthy young adults across 61 independent cohorts provides consistent, reproducible, and transcriptionally relevant functional connectome hubs in the resting human brain.

Deciphering cell-cell interactions and communication in the tumor microenvironment and unraveling intratumoral genetic heterogeneity via single-cell genomic sequencing

A tumor's heterogeneity has important implications in terms of its clonal origin, progression, stemness, and drug resistance. Therefore, because of its significance in treatment, it is important to understand the gene expression pattern of a single cell, track gene expression or mutation in heterogeneous cells, evaluate the clonal origin of cancer cells, and determine the selective evolution of different subpopulations of cancer cells. Researchers are able to trace a cell's mutation and identify different types of tumor cells by measuring the whole transcriptome with single-cell sequencing (scRNA-seq). This technology provides a better understanding of the molecular mechanisms driving tumor growth than that offered by traditional RNA sequencing methods. In addition, it has revealed changes in the mutations and functions of somatic cells as a tumor evolves; it has also clarified immune cell infiltration and activation. Research on scRNA-seq technology has recently advanced significantly, suggesting new strategies for the treatment of cancer. In short, cancer researchers have become increasingly dependent on scRNA-seq. This paper reviews the development, detection principles, and processes of scRNA-seq technology and their application in tumor research. It also considers potential clinical applications.

Integrating neuroimaging and gene expression data using the imaging transcriptomics toolbox

The integration of neuroimaging and transcriptomics data, Imaging Transcriptomics, is becoming increasingly popular but standardized workflows for its implementation are still lacking. We describe the Imaging Transcriptomics toolbox, a new package that implements a full imaging transcriptomics pipeline using a user-friendly, command line interface. This toolbox allows the user to identify patterns of gene expression which correlates with a specific neuroimaging phenotype and perform gene set enrichment analyses to inform the biological interpretation of the findings using up-to-date methods.

For complete details on the use and execution of this protocol, please refer to Martins et al. (2021).

•

Protocol for using the Imaging Transcriptomics toolbox

•

Identifies transcriptomic correlates of neuroimaging data

•

Performs gene enrichment analysis to contextualize findings

•

Standardized and user-friendly pipeline using robust statistics

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 Genetics in Epilepsy: Current Knowledge and New Perspectives

Epilepsy is a neurological network disease with genetics playing a much greater role than was previously appreciated. Unfortunately, the relationship between genetic basis and imaging phenotype is by no means simple. Imaging genetics integrates multidimensional datasets within a unified framework, providing a unique opportunity to pursue a global vision for epilepsy. This review delineates the current knowledge of underlying genetic mechanisms for brain networks in different epilepsy syndromes, particularly from a neural developmental perspective. Further, endophenotypes and their potential value are discussed. Finally, we highlight current challenges and provide perspectives for the future development of imaging genetics in epilepsy.

Sexually divergent development of depression-related brain networks during healthy human adolescence

Sexual differences in human brain development could be relevant to sex differences in the incidence of depression during adolescence. We tested for sex differences in parameters of normative brain network development using fMRI data on N = 298 healthy adolescents, aged 14 to 26 years, each scanned one to three times. Sexually divergent development of functional connectivity was located in the default mode network, limbic cortex, and subcortical nuclei. Females had a more "disruptive" pattern of development, where weak functional connectivity at age 14 became stronger during adolescence. This fMRI-derived map of sexually divergent brain network development was robustly colocated with i prior loci of reward-related brain activation ii a map of functional dysconnectivity in major depressive disorder (MDD), and iii an adult brain gene transcriptional pattern enriched for genes on the X chromosome, neurodevelopmental genes, and risk genes for MDD. We found normative sexual divergence in adolescent development of a cortico-subcortical brain functional network that is relevant to depression.

A synergistic core for human brain evolution and cognition

How does the organization of neural information processing enable humans' sophisticated cognition? Here we decompose functional interactions between brain regions into synergistic and redundant components, revealing their distinct information-processing roles. Combining functional and structural neuroimaging with meta-analytic results, we demonstrate that redundant interactions are predominantly associated with structurally coupled, modular sensorimotor processing. Synergistic interactions instead support integrative processes and complex cognition across higher-order brain networks. The human brain leverages synergistic information to a greater extent than nonhuman primates, with high-synergy association cortices exhibiting the highest degree of evolutionary cortical expansion. Synaptic density mapping from positron emission tomography and convergent molecular and metabolic evidence demonstrate that synergistic interactions are supported by receptor diversity and human-accelerated genes underpinning synaptic function. This information-resolved approach provides analytic tools to disentangle information integration from coupling, enabling richer, more accurate interpretations of functional connectivity, and illuminating how the human neurocognitive architecture navigates the trade-off between robustness and integration.

Spatially resolved transcriptomics provide a new method for cancer research

A major feature of cancer is the heterogeneity, both intratumoral and intertumoral. Traditional single-cell techniques have given us a comprehensive understanding of the biological characteristics of individual tumor cells, but the lack of spatial context of the transcriptome has limited the study of cell-to-cell interaction patterns and hindered further exploration of tumor heterogeneity. In recent years, the advent of spatially resolved transcriptomics (SRT) technology has made possible the multidimensional analysis of the tumor microenvironment in the context of intact tissues. Different SRT methods are applicable to different working ranges due to different working principles. In this paper, we review the advantages and disadvantages of various current SRT methods and the overall idea of applying these techniques to oncology studies, hoping to help researchers find breakthroughs. Finally, we discussed the future direction of SRT technology, and deeper investigation into the complex mechanisms of tumor development from different perspectives through multi-omics fusion, paving the way for precisely targeted tumor therapy.

Uncovering the genetic profiles underlying the intrinsic organization of the human cerebellum

The functional diversity of the human cerebellum is largely believed to be derived more from its extensive connections rather than being limited to its mostly invariant architecture. However, whether and how the determination of cerebellar connections in its intrinsic organization interact with microscale gene expression is still unknown. Here we decode the genetic profiles of the cerebellar functional organization by investigating the genetic substrates simultaneously linking cerebellar functional heterogeneity and its drivers, i.e., the connections. We not only identified 443 network-specific genes but also discovered that their co-expression pattern correlated strongly with intra-cerebellar functional connectivity (FC). Ninety of these genes were also linked to the FC of cortico-cerebellar cognitive-limbic networks. To further discover the biological functions of these genes, we performed a "virtual gene knock-out" by observing the change in the coupling between gene co-expression and FC and divided the genes into two subsets, i.e., a positive gene contribution indicator (GCI⁺) involved in cerebellar neurodevelopment and a negative gene set (GCI^-) related to neurotransmission. A more interesting finding is that GCI^- is significantly linked with the cerebellar connectivity-behavior association and many recognized brain diseases that are closely linked with the cerebellar functional abnormalities. Our results could collectively help to rethink the genetic substrates underlying the cerebellar functional organization and offer possible micro-macro interacted mechanistic interpretations of the cerebellum-involved high order functions and dysfunctions in neuropsychiatric disorders.

Transcriptomic and cellular decoding of functional brain connectivity changes reveal regional brain vulnerability to pro- and anti-inflammatory therapies

BACKGROUND

Systemic inflammation induces acute changes in mood, motivation and cognition that closely resemble those observed in depressed individuals. However, the mechanistic pathways linking peripheral inflammation to depression-like psychopathology via intermediate effects on brain function remain incompletely understood.

METHODS

We combined data from 30 patients initiating interferon-α treatment for Hepatitis-C and 20 anti-tumour necrosis factor (TNF) therapy for inflammatory arthritis and used resting-state functional magnetic resonance imaging to investigate acute effects of each treatment on regional global brain connectivity (GBC). We leveraged transcriptomic data from the Allen Human Brain Atlas to uncover potential biological and cellular pathways underpinning regional vulnerability to GBC changes induced by each treatment.

RESULTS

Interferon-α and anti-TNF therapies both produced differential small-to-medium sized decreases in regional GBC. However, these were observed within distinct brain regions and the regional patterns of GBC changes induced by each treatment did not correlate suggesting independent underlying processes. Further, the spatial distribution of these differential GBC decreases could be captured by multivariate patterns of constitutive regional expression of genes respectively related to: i) neuroinflammation and glial cells; and ii) glutamatergic neurotransmission and neurons. The extent to which each participant expressed patterns of GBC changes aligning with these patterns of transcriptomic vulnerability also correlated with both acute treatment-induced changes in interleukin-6 (IL-6) and, for Interferon-α, longer-term treatment-associated changes in depressive symptoms.

CONCLUSIONS

Together, we present two transcriptomic models separately linking regional vulnerability to the acute effects of interferon-α and anti-TNF treatments on brain function to glial neuroinflammation and glutamatergic neurotransmission. These findings generate hypotheses about two potential brain mechanisms through which bidirectional changes in peripheral inflammation may contribute to the development/resolution of psychopathology.

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

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

Neurogenetic traits outline vulnerability to cortical disruption in Parkinson's disease

The genetic traits that underlie vulnerability to neuronal damage across specific brain circuits in Parkinson's disease (PD) remain to be elucidated. In this study, we characterized the brain topological intersection between propagating connectivity networks in controls and PD participants and gene expression patterns across the human cortex - such as the SNCA gene. We observed that brain connectivity originated from PD-related pathology epicenters in the brainstem recapitulated the anatomical distribution of alpha-synuclein histopathology in postmortem data. We also discovered that the gene set most related to cortical propagation patterns of PD-related pathology was primarily involved in microtubule cellular components. Thus, this study sheds light on new avenues for enhancing detection of PD neuronal vulnerability via an evaluation of in vivo connectivity trajectories across the human brain and successful integration of neuroimaging-genetic strategies.

Hospitalisation for COVID-19 predicts long lasting cerebrovascular impairment: A prospective observational cohort study

Human coronavirus disease 2019 (COVID-19) due to severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has multiple neurological consequences, but its long-term effect on brain health is still uncertain. The cerebrovascular consequences of COVID-19 may also affect brain health. We studied the chronic effect of COVID-19 on cerebrovascular health, in relation to acute severity, adverse clinical outcomes and in contrast to control group data. Here we assess cerebrovascular health in 45 patients six months after hospitalisation for acute COVID-19 using the resting state fluctuation amplitudes (RSFA) from functional magnetic resonance imaging, in relation to disease severity and in contrast with 42 controls. Acute COVID-19 severity was indexed by COVID-19 WHO Progression Scale, inflammatory and coagulatory biomarkers. Chronic widespread changes in frontoparietal RSFA were related to the severity of the acute COVID-19 episode. This relationship was not explained by chronic cardiorespiratory dysfunction, age, or sex. The level of cerebrovascular dysfunction was associated with cognitive, mental, and physical health at follow-up. The principal findings were consistent across univariate and multivariate approaches. The results indicate chronic cerebrovascular impairment following severe acute COVID-19, with the potential for long-term consequences on cognitive function and mental wellbeing.

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

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

Gene expression associated with individual variability in intrinsic functional connectivity

It has been revealed that intersubject variability (ISV) in intrinsic functional connectivity (FC) is associated with a wide variety of cognitive and behavioral performances. However, the underlying organizational principle of ISV in FC and its related gene transcriptional profiles remain unclear. Using resting-state fMRI data from the Human Connectome Project (299 adult participants) and microarray gene expression data from the Allen Human Brain Atlas, we conducted a transcription-neuroimaging association study to investigate the spatial configurations of ISV in intrinsic FC and their associations with spatial gene transcriptional profiles. We found that the multimodal association cortices showed the greatest ISV in FC, while the unimodal cortices and subcortical areas showed the least ISV. Importantly, partial least squares regression analysis revealed that the transcriptional profiles of genes associated with human accelerated regions (HARs) could explain 31.29% of the variation in the spatial distribution of ISV in FC. The top-related genes in the transcriptional profiles were enriched for the development of the central nervous system, neurogenesis and the cellular components of synapse. Moreover, we observed that the effect of gene expression profile on the heterogeneous distribution of ISV in FC was significantly mediated by the cerebral blood flow configuration. These findings highlighted the spatial arrangement of ISV in FC and their coupling with variations in transcriptional profiles and cerebral blood flow supply.

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

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

NuMorph: Tools for cortical cellular phenotyping in tissue-cleared whole-brain images

Tissue-clearing methods allow every cell in the mouse brain to be imaged without physical sectioning. However, the computational tools currently available for cell quantification in cleared tissue images have been limited to counting sparse cell populations in stereotypical mice. Here, we introduce NuMorph, a group of analysis tools to quantify all nuclei and nuclear markers within the mouse cortex after clearing and imaging by light-sheet microscopy. We apply NuMorph to investigate two distinct mouse models: a Topoisomerase 1 (Top1) model with severe neurodegenerative deficits and a Neurofibromin 1 (Nf1) model with a more subtle brain overgrowth phenotype. In each case, we identify differential effects of gene deletion on individual cell-type counts and distribution across cortical regions that manifest as alterations of gross brain morphology. These results underline the value of whole-brain imaging approaches, and the tools are widely applicable for studying brain structure phenotypes at cellular resolution.

Divergent connectomic organization delineates genetic evolutionary traits in the human brain

The relationship between human brain connectomics and genetic evolutionary traits remains elusive due to the inherent challenges in combining complex associations within cerebral tissue. In this study, insights are provided about the relationship between connectomics, gene expression and divergent evolutionary pathways from non-human primates to humans. Using in vivo human brain resting-state data, we detected two co-existing idiosyncratic functional systems: the segregation network, in charge of module specialization, and the integration network, responsible for information flow. Their topology was approximated to whole-brain genetic expression (Allen Human Brain Atlas) and the co-localization patterns yielded that neuron communication functionalities-linked to Neuron Projection-were overrepresented cell traits. Homologue-orthologue comparisons using dN/dS-ratios bridged the gap between neurogenetic outcomes and biological data, summarizing the known evolutionary divergent pathways within the Homo Sapiens lineage. Evidence suggests that a crosstalk between functional specialization and information flow reflects putative biological qualities of brain architecture, such as neurite cellular functions like axonal or dendrite processes, hypothesized to have been selectively conserved in the species through positive selection. These findings expand our understanding of human brain function and unveil aspects of our cognitive trajectory in relation to our simian ancestors previously left unexplored.

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.

Genetic influences on hub connectivity of the human connectome

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

A generative network model of neurodevelopmental diversity in structural brain organization

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