Spatial patterns of genome-wide expression profiles reflect anatomic and fiber connectivity architecture of healthy human brain

Cerebral blood flow alterations and host genetic association in individuals with long COVID: A transcriptomic-neuroimaging study

Neuroimaging studies have indicated that altered cerebral blood flow (CBF) was associated with the long-term symptoms of postacute sequelae of SARS-CoV-2 infection (PASC), also known as "long COVID". COVID-19 and long COVID were found to be strongly associated with host gene expression. Nevertheless, the relationships between altered CBF, clinical symptoms, and gene expression in the central nervous system (CNS) remain unclear in individuals with long COVID. This study aimed to explore the genetic mechanisms of CBF abnormalities in individuals with long COVID by transcriptomic-neuroimaging spatial association. Lower CBF in the left frontal-temporal gyrus was associated with higher fatigue and worse cognition in individuals with long COVID. This CBF pattern was spatially associated with the expression of 2,178 genes, which were enriched in the molecular functions and biological pathways of COVID-19. Our study suggested that lower CBF is associated with persistent clinical symptoms in long COVID individuals, possibly as a consequence of the complex interactions among multiple COVID-19-related genes, which contributes to our understanding of the impact of adverse CNS outcomes and the trajectory of development to long COVID.

DSM: Deep sequential model for complete neuronal morphology representation and feature extraction

The full morphology of single neurons is indispensable for understanding cell types, the basic building blocks in brains. Projecting trajectories are critical to extracting biologically relevant information from neuron morphologies, as they provide valuable information for both connectivity and cell identity. We developed an artificial intelligence method, deep sequential model (DSM), to extract concise, cell-type-defining features from projections across brain regions. DSM achieves more than 90% accuracy in classifying 12 major neuron projection types without compromising performance when spatial noise is present. Such remarkable robustness enabled us to efficiently manage and analyze several major full-morphology data sources, showcasing how characteristic long projections can define cell identities. We also succeeded in applying our model to both discovering previously unknown neuron subtypes and analyzing exceptional co-expressed genes involved in neuron projection circuits.

Spatial cell-type enrichment predicts mouse brain connectivity

A fundamental neuroscience topic is the link between the brain's molecular, cellular, and cytoarchitectonic properties and structural connectivity. Recent studies relate inter-regional connectivity to gene expression, but the relationship to regional cell-type distributions remains understudied. Here, we utilize whole-brain mapping of neuronal and non-neuronal subtypes via the matrix inversion and subset selection algorithm to model inter-regional connectivity as a function of regional cell-type composition with machine learning. We deployed random forest algorithms for predicting connectivity from cell-type densities, demonstrating surprisingly strong prediction accuracy of cell types in general, and particular non-neuronal cells such as oligodendrocytes. We found evidence of a strong distance dependency in the cell connectivity relationship, with layer-specific excitatory neurons contributing the most for long-range connectivity, while vascular and astroglia were salient for short-range connections. Our results demonstrate a link between cell types and connectivity, providing a roadmap for examining this relationship in other species, including humans.

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.

BrainStat: A toolbox for brain-wide statistics and multimodal feature associations

Analysis and interpretation of neuroimaging datasets has become a multidisciplinary endeavor, relying not only on statistical methods, but increasingly on associations with respect to other brain-derived features such as gene expression, histological data, and functional as well as cognitive architectures. Here, we introduce BrainStat - a toolbox for (i) univariate and multivariate linear models in volumetric and surface-based brain imaging datasets, and (ii) multidomain feature association of results with respect to spatial maps of post-mortem gene expression and histology, task-based fMRI meta-analysis, as well as resting-state fMRI motifs across several common surface templates. The combination of statistics and feature associations into a turnkey toolbox streamlines analytical processes and accelerates cross-modal research. The toolbox is implemented in both Python and MATLAB, two widely used programming languages in the neuroimaging and neuroinformatics communities. BrainStat is openly available and complemented by an expandable documentation.

Alzheimer's genetic risk effects on cerebral blood flow across the lifespan are proximal to gene expression

Cerebrovascular dysregulation such as altered cerebral blood flow (CBF) can be observed in Alzheimer's disease (AD) and may precede symptom onset. Genome wide association studies show that AD has a polygenic aetiology, providing a tool for studying AD susceptibility across the lifespan. Here, we ascertain whether the AD genetic risk effects on CBF previously observed (Chandler et al., 2019) are also present in later life. Consistent with our prior observations, AD genetic risk score (AD-GRS) was associated with reduced CBF in the ADNI sample. The regional association between AD-GRS and CBF were also spatially similar. Furthermore, CBF was related to the regional mRNA transcript expression of AD risk genes proximal to AD-GRS risk loci. These observations suggest that AD risk alleles may reduce neurovascular process such as CBF, potentially via mechanisms such as regional expression of proximal AD risk genes as an antecedent AD pathophysiology. Our observations help establish processes that underpin AD genetic risk-related reductions in CBF as a therapeutic target prior to the onset of neurodegeneration.

Linking Central Gene Expression Patterns and Mental States Using Transcriptomics and Large-Scale Meta-Analysis of fMRI Data: A Tutorial and Example Using the Oxytocin Signaling Pathway

The measurement of gene expression levels in the human brain can help accelerate our understanding of complex mental states and psychiatric illnesses. Mental states are typically associated with whole-brain networks; however, gene expression levels from postmortem brain samples have traditionally been measured in a limited number of brain regions due to resource limitations. The recent availability of whole-brain gene expression data from the Allen Human Brain Atlas (AHBA) provides the opportunity to generate gene expression patterns for over 20,000 genes. By linking these expression patterns with brain activity patterns that are associated with specific mental states, researchers can better understand which genes may support given mental states, via forward inference. Conversely, reverse inference can also be used to determine which mental state activation patterns are most strongly associated with a given gene expression map. This chapter provides a step-by-step guide on how to use the AHBA in conjunction with the NeuroSynth fMRI meta-analysis tool to identify the mental state correlates of specific gene expression patterns, using genes from oxytocin signaling pathway as an example. We also demonstrate how to perform an out-of-sample validation and assess the specificity of results for genes of interest.

Where the genome meets the connectome: Understanding how genes shape human brain connectivity

The integration of modern neuroimaging methods with genetically informative designs and data can shed light on the molecular mechanisms underlying the structural and functional organization of the human connectome. Here, we review studies that have investigated the genetic basis of human brain network structure and function through three complementary frameworks: (1) the quantification of phenotypic heritability through classical twin designs; (2) the identification of specific DNA variants linked to phenotypic variation through association and related studies; and (3) the analysis of correlations between spatial variations in imaging phenotypes and gene expression profiles through the integration of neuroimaging and transcriptional atlas data. We consider the basic foundations, strengths, limitations, and discoveries associated with each approach. We present converging evidence to indicate that anatomical connectivity is under stronger genetic influence than functional connectivity and that genetic influences are not uniformly distributed throughout the brain, with phenotypic variation in certain regions and connections being under stronger genetic control than others. We also consider how the combination of imaging and genetics can be used to understand the ways in which genes may drive brain dysfunction in different clinical disorders.

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.

A single cell transcriptomics map of paracrine networks in the intrinsic cardiac nervous system

We developed a spatially-tracked single neuron transcriptomics map of an intrinsic cardiac ganglion, the right atrial ganglionic plexus (RAGP) that is a critical mediator of sinoatrial node (SAN) activity. This 3D representation of RAGP used neuronal tracing to extensively map the spatial distribution of the subset of neurons that project to the SAN. RNA-seq of laser capture microdissected neurons revealed a distinct composition of RAGP neurons compared to the central nervous system and a surprising finding that cholinergic and catecholaminergic markers are coexpressed, suggesting multipotential phenotypes that can drive neuroplasticity within RAGP. High-throughput qPCR of hundreds of laser capture microdissected single neurons confirmed these findings and revealed a high dimensionality of neuromodulatory factors that contribute to dynamic control of the heart. Neuropeptide-receptor coexpression analysis revealed a combinatorial paracrine neuromodulatory network within RAGP informing follow-on studies on the vagal control of RAGP to regulate cardiac function in health and disease.

Brain Gene Expression Pattern Correlated with the Differential Brain Activation by Pain and Touch in Humans

Genes involved in pain and touch sensations have been studied extensively, but very few studies have tried to link them with neural activities in the brain. Here, we aimed to identify genes preferentially correlated to painful activation patterns by linking the spatial patterns of gene expression of Allen Human Brain Atlas with the pain-elicited neural responses in the human brain, with a parallel, control analysis for identification of genes preferentially correlated to tactile activation patterns. We identified 1828 genes whose expression patterns preferentially correlated to painful activation patterns and 411 genes whose expression patterns preferentially correlated to tactile activation pattern at the cortical level. In contrast to the enrichment for astrocyte and inhibitory synaptic transmission of genes preferentially correlated to tactile activation, the genes preferentially correlated to painful activation were mainly enriched for neuron and opioid- and addiction-related pathways and showed significant overlap with pain-related genes identified in previous studies. These findings not only provide important evidence for the differential genetic architectures of specific brain activation patterns elicited by painful and tactile stimuli but also validate a new approach to studying pain- and touch-related genes more directly from the perspective of neural responses in the human brain.

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

Transcriptomic atlases have improved our understanding of the correlations between gene-expression patterns and spatially varying properties of brain structure and function. Gene-category enrichment analysis (GCEA) is a common method to identify functional gene categories that drive these associations, using gene-to-category annotation systems like the Gene Ontology (GO). Here, we show that applying standard GCEA methodology to spatial transcriptomic data is affected by substantial false-positive bias, with GO categories displaying an over 500-fold average inflation of false-positive associations with random neural phenotypes in mouse and human. The estimated false-positive rate of a GO category is associated with its rate of being reported as significantly enriched in the literature, suggesting that published reports are affected by this false-positive bias. We show that within-category gene-gene coexpression and spatial autocorrelation are key drivers of the false-positive bias and introduce flexible ensemble-based null models that can account for these effects, made available as a software toolbox.

Imaging Transcriptomics in Neurodegenerative Diseases

Imaging transcriptomics investigates the relationship between neuroanatomical/neuroimaging features and gene expression. The spatial and temporal distribution of the expressed genes and their pattern of spreading over time can contribute to elucidating cellular and molecular processes involved in neurodegeneration. In this study, we review recent findings regarding the correlation between neuroimaging and expression data in neurodegenerative diseases with a focus on Alzheimer's disease and Parkinson's disease. An association between gene expression data and different neuroimaging neurodegeneration features, such as R2 relaxation time and volumetric cortical atrophy, was established. Several positive and negative expression correlations were identified, and they confirmed the focal nature of neurodegeneration. Positively correlated genes were associated with cell motility, immune system activity, neuroinflammation, and microglia. Data from connectome studies support the hypothesis of selective network vulnerability and a temporal spreading pattern in neurodegenerative pathologies. Genes related to cellular mobility and transport are overexpressed in the neuroimaging-defined delineated areas of degeneration. In addition, expression enrichment of genes involved in immunological processes in vulnerable regions-such as the Toll-like receptor, a receptor involved in innate immunity-plays a major role in neuroinflammation in neurodegenerative diseases. However, substantial limitations must be overcome in future studies: the lack of high-quality resolution expression data, the lack of standardized study protocols, and insufficient sensitive early stage neuroimaging markers of degeneration. Identifying neuroimaging and expression prodromal biomarkers and investigating their causal relation in the preclinical disease stage may enable early targeted therapy before the onset of irreversible brain changes.

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

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

Gene Expression Correlates of the Cortical Network Underlying Sentence Processing

A pivotal question in modern neuroscience is which genes regulate brain circuits that underlie cognitive functions. However, the field is still in its infancy. Here we report an integrated investigation of the high-level language network (i.e., sentence-processing network) in the human cerebral cortex, combining regional gene expression profiles, task fMRI, large-scale neuroimaging meta-analysis, and resting-state functional network approaches. We revealed reliable gene expression-functional network correlations using three different network definition strategies, and identified a consensus set of genes related to connectivity within the sentence-processing network. The genes involved showed enrichment for neural development and actin-related functions, as well as association signals with autism, which can involve disrupted language functioning. Our findings help elucidate the molecular basis of the brain's infrastructure for language. The integrative approach described here will be useful for studying other complex cognitive traits.

Neural Transcription Correlates of Multimodal Cortical Phenotypes during Development

During development, cellular events such as cell proliferation, migration, and synaptogenesis determine the structural organization of the brain. These processes are driven in part by spatiotemporally regulated gene expression. We investigated how the genetic signatures of specific neural cell types shape cortical organization of the human brain throughout infancy and childhood. Using a transcriptional atlas and in vivo magnetic resonance imaging (MRI) data, we demonstrated time-dependent associations between the expression levels of neuronal and glial genes and cortical macro- and microstructure. Neonatal cortical phenotypes were associated with prenatal glial but not neuronal gene expression. These associations reflect cell migration and proliferation during fetal development. Childhood cortical phenotypes were associated with neuronal and astrocyte gene expression related to synaptic signaling processes, reflecting the refinement of cortical connections. These findings indicate that sequential developmental stages contribute to distinct MRI measures at different time points. This helps to bridge the gap between the genetic mechanisms driving cellular changes and widely used neuroimaging techniques.

Application of Computational Biology to Decode Brain Transcriptomes

The rapid development of high-throughput sequencing technologies has generated massive valuable brain transcriptome atlases, providing great opportunities for systematically investigating gene expression characteristics across various brain regions throughout a series of developmental stages. Recent studies have revealed that the transcriptional architecture is the key to interpreting the molecular mechanisms of brain complexity. However, our knowledge of brain transcriptional characteristics remains very limited. With the immense efforts to generate high-quality brain transcriptome atlases, new computational approaches to analyze these high-dimensional multivariate data are greatly needed. In this review, we summarize some public resources for brain transcriptome atlases and discuss the general computational pipelines that are commonly used in this field, which would aid in making new discoveries in brain development and disorders.

Making Sense of Pharmacovigilance and Drug Adverse Event Reporting: Comparative Similarity Association Analysis Using AI Machine Learning Algorithms in Dogs and Cats

Drug-associated adverse events cause approximately 30 billion dollars a year of added health care expense, along with negative health outcomes including patient death. This constitutes a major public health concern. The US Food and Drug Administration (FDA) requires drug labeling to include potential adverse effects for each newly developed drug product. With the advancement in incidence of adverse drug events (ADEs) and potential adverse drug events, published studies have mainly concluded potential ADEs from labeling documents obtained from the FDA's preapproval clinical trials, and very few analyzed their research work based on reported ADEs after widespread use of a drug to animal subjects. The aforesaid procedure of deriving practice based on information from preapproval labeling may misrepresent or deprecate the incidence and prevalence of specific ADEs. In this study, we make the most of the recently disseminated ADE data by the FDA for animal drugs and devices used in animals to address this public and welfare concern. For this purpose, we implemented 5 different methods (Pearson distance, Spearman distance, cosine distance, Yule distance, and Euclidean distance) to determine the most efficient and robust approach to properly discover highly associated ADEs from the reported data and accurately exclude noise-induced reported events, while maintaining a high level of correlation precision. Our comparative analysis of ADEs based on an artificial intelligence (AI) approach for the 5 robust similarity methods revealed high ADE associations for 2 drugs used in dogs and cats. In addition, the described distance methods systematically analyzed and compared ADEs from the drug labeling sections with a specific emphasis on analyzing serious ADEs. Our finding showed that the cosine method significantly outperformed all the other methods by correctly detecting and validating ADEs based on the comparative similarity association analysis compared with ADEs reported by preapproval clinical trials, premarket testing, or postapproval complication experience of FDA-approved animal drugs.

Uncovering the Transcriptional Correlates of Hub Connectivity in Neural Networks

Connections in nervous systems are disproportionately concentrated on a small subset of neural elements that act as network hubs. Hubs have been found across different species and scales ranging from C. elegans to mouse, rat, cat, macaque, and human, suggesting a role for genetic influences. The recent availability of brain-wide gene expression atlases provides new opportunities for mapping the transcriptional correlates of large-scale network-level phenotypes. Here we review studies that use these atlases to investigate gene expression patterns associated with hub connectivity in neural networks and present evidence that some of these patterns are conserved across species and scales.

Brain structural alterations are distributed following functional, anatomic and genetic connectivity

The pathological brain is characterized by distributed morphological or structural alterations in the grey matter, which tend to follow identifiable network-like patterns. We analysed the patterns formed by these alterations (increased and decreased grey matter values detected with the voxel-based morphometry technique) conducting an extensive transdiagnostic search of voxel-based morphometry studies in a large variety of brain disorders. We devised an innovative method to construct the networks formed by the structurally co-altered brain areas, which can be considered as pathological structural co-alteration patterns, and to compare these patterns with three associated types of connectivity profiles (functional, anatomical, and genetic). Our study provides transdiagnostical evidence that structural co-alterations are influenced by connectivity constraints rather than being randomly distributed. Analyses show that although all the three types of connectivity taken together can account for and predict with good statistical accuracy, the shape and temporal development of the co-alteration patterns, functional connectivity offers the better account of the structural co-alteration, followed by anatomic and genetic connectivity. These results shed new light on the possible mechanisms at the root of neuropathological processes and open exciting prospects in the quest for a better understanding of brain disorders.

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

The recent availability of comprehensive, brain-wide gene expression atlases such as the Allen Human Brain Atlas (AHBA) has opened new opportunities for understanding how spatial variations on molecular scale relate to the macroscopic neuroimaging phenotypes. A rapidly growing body of literature is demonstrating relationships between gene expression and diverse properties of brain structure and function, but approaches for combining expression atlas data with neuroimaging are highly inconsistent, with substantial variations in how the expression data are processed. The degree to which these methodological variations affect findings is unclear. Here, we outline a seven-step analysis pipeline for relating brain-wide transcriptomic and neuroimaging data and compare how different processing choices influence the resulting data. We suggest that studies using the AHBA should work towards a unified data processing pipeline to ensure consistent and reproducible results in this burgeoning field.

Molecular properties underlying regional vulnerability to Alzheimer's disease pathology

Amyloid deposition and neurofibrillary degeneration in Alzheimer's disease specifically affect discrete neuronal systems, but the underlying mechanisms that render some brain regions more vulnerable to Alzheimer's disease pathology than others remain largely unknown. Here we studied molecular properties underlying these distinct regional vulnerabilities by analysing Alzheimer's disease-typical neuroimaging patterns of amyloid deposition and neurodegeneration in relation to regional gene expression profiles of the human brain. Graded patterns of brain-wide vulnerability to amyloid deposition and neurodegeneration in Alzheimer's disease were estimated by contrasting multimodal amyloid-sensitive PET and structural MRI data between patients with Alzheimer's disease dementia (n = 76) and healthy controls (n = 126) enrolled in the Alzheimer's Disease Neuroimaging Initiative (ADNI). Regional gene expression profiles were derived from brain-wide microarray measurements provided by the Allen brain atlas of the adult human brain transcriptome. In a hypothesis-driven analysis focusing on the genes coding for the amyloid precursor (APP) and tau proteins (MAPT), regional expression levels of APP were positively correlated with the severity of regional amyloid deposition (r = 0.44, P = 0.009), but not neurodegeneration (r = 0.01, P = 0.96), whereas the opposite pattern was observed for MAPT (neurodegeneration: r = 0.46, P = 0.006; amyloid: r = 0.08, P = 0.65). Using explorative gene set enrichment analysis, amyloid-vulnerable regions were found to be characterized by relatively low expression levels of gene sets implicated in protein synthesis and mitochondrial respiration. By contrast, neurodegeneration-vulnerable regions were characterized by relatively high expression levels of gene sets broadly implicated in neural plasticity, with biological functions ranging from neurite outgrowth and synaptic contact over intracellular signalling cascades to proteoglycan metabolism. At the individual gene level this data-driven analysis further corroborated the association between neurodegeneration and MAPT expression, and additionally identified associations with known tau kinases (CDK5, MAPK1/ERK2) alongside components of their intracellular (Ras-ERK) activation pathways. Sensitivity analyses showed that these pathology-specific imaging-genetic associations were largely robust against changes in some of the methodological parameters, including variation in the brain donor sample used for estimating regional gene expression profiles, and local variations in the Alzheimer's disease-typical imaging patterns when these were derived from an independent patient cohort (BioFINDER study). These findings highlight that the regionally selective vulnerability to Alzheimer's disease pathology relates to specific molecular-functional properties of the affected neural systems, and that the implicated biochemical pathways largely differ for amyloid accumulation versus neurodegeneration. The data provide novel insights into the complex pathophysiological mechanisms of Alzheimer's disease and point to pathology-specific treatment targets that warrant further exploration in independent studies.

Regional vulnerability in Alzheimer's disease: The role of cell-autonomous and transneuronal processes

INTRODUCTION

The stereotypical progression of Alzheimer's disease (AD) pathology is not fully understood. The selective impact of AD on distinct regions has led the field to question if innate vulnerability exists. This study aims to determine if the causative factors of regional vulnerability are dependent on cell-autonomous or transneuronal (non-cell autonomous) processes.

METHODS

Using mathematical and statistical models, we analyzed the contribution of cell-autonomous and non-cell autonomous factors to predictive linear models of AD pathology.

RESULTS

Results indicate gene expression as a weak contributor to predictive linear models of AD. Instead, the network diffusion model acts as a strong predictor of observed AD atrophy and hypometabolism.

DISCUSSION

We propose a convenient methodology for identifying genes and their role in determining AD topography, in comparison with network spread. Results reinforce the role of transneuronal network spread on disease progression and suggest that innate gene expression plays a secondary role in seeding and subsequent disease progression.

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

The quantification of big pools of diverse molecules provides important insights on brain function, but is often restricted to a limited number of observations, which impairs integration with other modalities. To resolve this issue, a method allowing for the prediction of mRNA expression in the entire brain based on microarray data provided in the Allen Human Brain Atlas was developed. Microarray data of 3702 samples from 6 brain donors was registered to MNI and cortical surface space using FreeSurfer. For each of 18,686 genes, spatial dependence of transcription was assessed using variogram modelling. Variogram models were employed in Gaussian process regression to calculate best linear unbiased predictions for gene expression at all locations represented in well-established imaging atlases for cortex, subcortical structures and cerebellum. For validation, predicted whole-brain transcription of the HTR1A gene was correlated with [carbonyl-¹¹C]WAY-100635 positron emission tomography data collected from 30 healthy subjects. Prediction results showed minimal bias ranging within ±0.016 (cortical surface), ±0.12 (subcortical regions) and ±0.14 (cerebellum) in units of log2 expression intensity for all genes. Across genes, the correlation of predicted and observed mRNA expression in leave-one-out cross-validation correlated with the strength of spatial dependence (cortical surface: r = 0.91, subcortical regions: r = 0.85, cerebellum: r = 0.84). 816 out of 18,686 genes exhibited a high spatial dependence accounting for more than 50% of variance in the difference of gene expression on the cortical surface. In subcortical regions and cerebellum, different sets of genes were implicated by high spatially structured variability. For the serotonin 1A receptor, correlation between PET binding potentials and predicted comprehensive mRNA expression was markedly higher (Spearman ρ = 0.72 for cortical surface, ρ = 0.84 for subcortical regions) than correlation of PET and discrete samples only (ρ = 0.55 and ρ = 0.63, respectively). Prediction of mRNA expression in the entire human brain allows for intuitive visualization of gene transcription and seamless integration in multimodal analysis without bias arising from non-uniform distribution of available samples. Extension of this methodology promises to facilitate translation of omics research and enable investigation of human brain function at a systems level.

Hub connectivity, neuronal diversity, and gene expression in the Caenorhabditis elegans connectome

Studies of nervous system connectivity, in a wide variety of species and at different scales of resolution, have identified several highly conserved motifs of network organization. One such motif is a heterogeneous distribution of connectivity across neural elements, such that some elements act as highly connected and functionally important network hubs. These brain network hubs are also densely interconnected, forming a so-called rich club. Recent work in mouse has identified a distinctive transcriptional signature of neural hubs, characterized by tightly coupled expression of oxidative metabolism genes, with similar genes characterizing macroscale inter-modular hub regions of the human cortex. Here, we sought to determine whether hubs of the neuronal C. elegans connectome also show tightly coupled gene expression. Using open data on the chemical and electrical connectivity of 279 C. elegans neurons, and binary gene expression data for each neuron across 948 genes, we computed a correlated gene expression score for each pair of neurons, providing a measure of their gene expression similarity. We demonstrate that connections between hub neurons are the most similar in their gene expression while connections between nonhubs are the least similar. Genes with the greatest contribution to this effect are involved in glutamatergic and cholinergic signaling, and other communication processes. We further show that coupled expression between hub neurons cannot be explained by their neuronal subtype (i.e., sensory, motor, or interneuron), separation distance, chemically secreted neurotransmitter, birth time, pairwise lineage distance, or their topological module affiliation. Instead, this coupling is intrinsically linked to the identity of most hubs as command interneurons, a specific class of interneurons that regulates locomotion. Our results suggest that neural hubs may possess a distinctive transcriptional signature, preserved across scales and species, that is related to the involvement of hubs in regulating the higher-order behaviors of a given organism.

Regional expression of genes mediating trans-synaptic alpha-synuclein transfer predicts regional atrophy in Parkinson disease

Multiple genes have been implicated in Parkinson disease pathogenesis, but the relationship between regional expression of these genes and regional dysfunction across the brain is unknown. We address this question by joint analysis of high resolution magnetic resonance imaging data from the Parkinson's Progression Markers Initiative and regional genetic microarray expression data from the Allen Brain Atlas. Regional brain atrophy and genetic expression was co-registered to a common 86 region brain atlas and robust multivariable regression analysis was performed to identify genetic predictors of regional brain atrophy. Top candidate genes from GWAS analysis, as well as genes implicated in trans-synaptic alpha-synuclein transfer and autosomal recessive PD were included in our analysis. We identify three genes with expression patterns that are highly significant predictors of regional brain atrophy. The two most significant predictors are LAG3 and RAB5A, genes implicated in trans-synaptic synuclein transfer. Other well-validated PD-related genes do not have expression patterns that predict regional atrophy, suggesting that they may serve other roles such as disease initiation factors.

•

Joint volumetric and microarray analysis identifies gene expression patterns that predict the PD atrophy pattern.

•

The most highly predictive genes, LAG3 and RAB5A, are implicated in trans-synaptic alpha-synuclein transfer.

•

The expression patterns of alpha-synuclein and otherPD-related genes do not predict atrophy.

Gene networks show associations with seed region connectivity

Primary patterns in adult brain connectivity are established during development by coordinated networks of transiently expressed genes; however, neural networks remain malleable throughout life. The present study hypothesizes that structural connectivity from key seed regions may induce effects on their connected targets, which are reflected in gene expression at those targeted regions. To test this hypothesis, analyses were performed on data from two brains from the Allen Human Brain Atlas, for which both gene expression and DW-MRI were available. Structural connectivity was estimated from the DW-MRI data and an approach motivated by network topology, that is, weighted gene coexpression network analysis (WGCNA), was used to cluster genes with similar patterns of expression across the brain. Group exponential lasso models were then used to predict gene cluster expression summaries as a function of seed region structural connectivity patterns. In several gene clusters, brain regions located in the brain stem, diencephalon, and hippocampal formation were identified that have significant predictive power for these expression summaries. These connectivity-associated clusters are enriched in genes associated with synaptic signaling and brain plasticity. Furthermore, using seed region based connectivity provides a novel perspective in understanding relationships between gene expression and connectivity. Hum Brain Mapp 38:3126-3140, 2017. © 2017 Wiley Periodicals, Inc.

Brain transcriptome atlases: a computational perspective

The immense complexity of the mammalian brain is largely reflected in the underlying molecular signatures of its billions of cells. Brain transcriptome atlases provide valuable insights into gene expression patterns across different brain areas throughout the course of development. Such atlases allow researchers to probe the molecular mechanisms which define neuronal identities, neuroanatomy, and patterns of connectivity. Despite the immense effort put into generating such atlases, to answer fundamental questions in neuroscience, an even greater effort is needed to develop methods to probe the resulting high-dimensional multivariate data. We provide a comprehensive overview of the various computational methods used to analyze brain transcriptome atlases.

Cognitive genomics: Linking genes to behavior in the human brain

Correlations of genetic variation in DNA with functional brain activity have already provided a starting point for delving into human cognitive mechanisms. However, these analyses do not provide the specific genes driving the associations, which are complicated by intergenic localization as well as tissue-specific epigenetics and expression. The use of brain-derived expression datasets could build upon the foundation of these initial genetic insights and yield genes and molecular pathways for testing new hypotheses regarding the molecular bases of human brain development, cognition, and disease. Thus, coupling these human brain gene expression data with measurements of brain activity may provide genes with critical roles in brain function. However, these brain gene expression datasets have their own set of caveats, most notably a reliance on postmortem tissue. In this perspective, I summarize and examine the progress that has been made in this realm to date, and discuss the various frontiers remaining, such as the inclusion of cell-type-specific information, additional physiological measurements, and genomic data from patient cohorts.

HASE: Framework for efficient high-dimensional association analyses

High-throughput technology can now provide rich information on a person’s biological makeup and environmental surroundings. Important discoveries have been made by relating these data to various health outcomes in fields such as genomics, proteomics, and medical imaging. However, cross-investigations between several high-throughput technologies remain impractical due to demanding computational requirements (hundreds of years of computing resources) and unsuitability for collaborative settings (terabytes of data to share). Here we introduce the HASE framework that overcomes both of these issues. Our approach dramatically reduces computational time from years to only hours and also requires several gigabytes to be exchanged between collaborators. We implemented a novel meta-analytical method that yields identical power as pooled analyses without the need of sharing individual participant data. The efficiency of the framework is illustrated by associating 9 million genetic variants with 1.5 million brain imaging voxels in three cohorts (total N = 4,034) followed by meta-analysis, on a standard computational infrastructure. These experiments indicate that HASE facilitates high-dimensional association studies enabling large multicenter association studies for future discoveries.

Automatic Segmentation of Drosophila Neural Compartments Using GAL4 Expression Data Reveals Novel Visual Pathways

Identifying distinct anatomical structures within the brain and developing genetic tools to target them are fundamental steps for understanding brain function. We hypothesize that enhancer expression patterns can be used to automatically identify functional units such as neuropils and fiber tracts. We used two recent, genome-scale Drosophila GAL4 libraries and associated confocal image datasets to segment large brain regions into smaller subvolumes. Our results (available at https://strawlab.org/braincode) support this hypothesis because regions with well-known anatomy, namely the antennal lobes and central complex, were automatically segmented into familiar compartments. The basis for the structural assignment is clustering of voxels based on patterns of enhancer expression. These initial clusters are agglomerated to make hierarchical predictions of structure. We applied the algorithm to central brain regions receiving input from the optic lobes. Based on the automated segmentation and manual validation, we can identify and provide promising driver lines for 11 previously identified and 14 novel types of visual projection neurons and their associated optic glomeruli. The same strategy can be used in other brain regions and likely other species, including vertebrates.

•

Genome-scale enhancer expression patterns can be used to predict brain structure

•

Automated clustering of images finds known structures such as olfactory glomeruli

•

Results identify GAL4 lines with strong expression in the predicted structures

•

We validate novel predictions to reveal previously undescribed optic glomeruli

Global analysis of gene expression and projection target correlations in the mouse brain

Recent studies have shown that projection targets in the mouse neocortex are correlated with their gene expression patterns. However, a brain-wide quantitative analysis of the relationship between voxel genetic composition and their projection targets is lacking to date. Here we extended those studies to perform a global, integrative analysis of gene expression and projection target correlations in the mouse brain. By using the Allen Brain Atlas data, we analyzed the relationship between gene expression and projection targets. We first visualized and clustered the two data sets separately and showed that they both exhibit strong spatial autocorrelation. Building upon this initial analysis, we conducted an integrative correlation analysis of the two data sets while correcting for their spatial autocorrelation. This resulted in a correlation of 0.19 with significant p value. We further identified the top genes responsible for this correlation using two greedy gene ranking techniques. Using only the top genes identified by those techniques, we recomputed the correlation between these two data sets. This led to correlation values up to 0.49 with significant p values. Our results illustrated that although the target specificity of neurons is in fact complex and diverse, yet they are strongly affected by their genetic and molecular compositions.

Network Diffusion Model of Progression Predicts Longitudinal Patterns of Atrophy and Metabolism in Alzheimer's Disease

Alzheimer's disease pathology (AD) originates in the hippocampus and subsequently spreads to temporal, parietal, and prefrontal association cortices in a relatively stereotyped progression. Current evidence attributes this orderly progression to transneuronal transmission of misfolded proteins along the projection pathways of affected neurons. A network diffusion model was recently proposed to mathematically predict disease topography resulting from transneuronal transmission on the brain's connectivity network. Here, we use this model to predict future patterns of regional atrophy and metabolism from baseline regional patterns of 418 subjects. The model accurately predicts end-of-study regional atrophy and metabolism starting from baseline data, with significantly higher correlation strength than given by the baseline statistics directly. The model's rate parameter encapsulates overall atrophy progression rate; group analysis revealed this rate to depend on diagnosis as well as baseline cerebrospinal fluid (CSF) biomarker levels. This work helps validate the model as a prognostic tool for Alzheimer's disease assessment.

Differences in human cortical gene expression match the temporal properties of large-scale functional networks

We explore the relationships between the cortex functional organization and genetic expression (as provided by the Allen Human Brain Atlas). Previous work suggests that functional cortical networks (resting state and task based) are organized as two large networks (differentiated by their preferred information processing mode) shaped like two rings. The first ring–Visual-Sensorimotor-Auditory (VSA)–comprises visual, auditory, somatosensory, and motor cortices that process real time world interactions. The second ring–Parieto-Temporo-Frontal (PTF)–comprises parietal, temporal, and frontal regions with networks dedicated to cognitive functions, emotions, biological needs, and internally driven rhythms. We found–with correspondence analysis–that the patterns of expression of the 938 genes most differentially expressed across the cortex organized the cortex into two sets of regions that match the two rings. We confirmed this result using discriminant correspondence analysis by showing that the genetic profiles of cortical regions can reliably predict to what ring these regions belong. We found that several of the proteins–coded by genes that most differentiate the rings–were involved in neuronal information processing such as ionic channels and neurotransmitter release. The systematic study of families of genes revealed specific proteins within families preferentially expressed in each ring. The results showed strong congruence between the preferential expression of subsets of genes, temporal properties of the proteins they code, and the preferred processing modes of the rings. Ionic channels and release-related proteins more expressed in the VSA ring favor temporal precision of fast evoked neural transmission (Sodium channels SCNA1, SCNB1 potassium channel KCNA1, calcium channel CACNA2D2, Synaptotagmin SYT2, Complexin CPLX1, Synaptobrevin VAMP1). Conversely, genes expressed in the PTF ring favor slower, sustained, or rhythmic activation (Sodium channels SCNA3, SCNB3, SCN9A potassium channels KCNF1, KCNG1) and facilitate spontaneous transmitter release (calcium channel CACNA1H, Synaptotagmins SYT5, Complexin CPLX3, and synaptobrevin VAMP2).

Molecular networks and the evolution of human cognitive specializations

Inroads into elucidating the origins of human cognitive specializations have taken many forms, including genetic, genomic, anatomical, and behavioral assays that typically compare humans to non-human primates. While the integration of all of these approaches is essential for ultimately understanding human cognition, here, we review the usefulness of coexpression network analysis for specifically addressing this question. An increasing number of studies have incorporated coexpression networks into brain expression studies comparing species, disease versus control tissue, brain regions, or developmental time periods. A clearer picture has emerged of the key genes driving brain evolution, as well as the developmental and regional contributions of gene expression patterns important for normal brain development and those misregulated in cognitive diseases.

Structure-Function Relationships behind the Phenomenon of Cognitive Resilience in Neurology: Insights for Neuroscience and Medicine

The phenomenon of cognitive resilience, that is, the dynamical preservation of normal functions despite neurological disorders, demonstrates that cognition can be highly robust to devastating brain injury. Here, cognitive resilience is considered across a range of neurological conditions. Simple computational models of structure-function relationships are used to discuss hypotheses about the neural mechanisms of resilience. Resilience expresses functional redundancies in brain networks and suggests a process of dynamic rerouting of brain signals. This process is underlined by a global renormalization of effective connectivity, capable of restoring information transfer between spared brain structures via alternate pathways. Local mechanisms of synaptic plasticity mediate the renormalization at the lowest level of implementation, but it is also driven by top-down cognition, with a key role of self-awareness in fostering resilience. The presence of abstraction layers in brain computation and networking is hypothesized to account for the renormalization process. Future research directions and challenges are discussed regarding the understanding and control of resilience based on multimodal neuroimaging and computational neuroscience. The study of resilience will illuminate ways by which the brain can overcome adversity and help inform prevention and treatment strategies. It is relevant to combating the negative neuropsychological impact of aging and fostering cognitive enhancement.