Relationships between gene expression and brain wiring in the adult rodent brain

The neuroanatomical organization of the hypothalamus is driven by spatial and topological efficiency

The hypothalamus in the mammalian brain is responsible for regulating functions associated with survival and reproduction representing a complex set of highly interconnected, yet anatomically and functionally distinct, sub-regions. It remains unclear what factors drive the spatial organization of sub-regions within the hypothalamus. One potential factor may be structural connectivity of the network that promotes efficient function with well-connected sub-regions placed closer together geometrically, i.e., the strongest axonal signal transferred through the shortest geometrical distance. To empirically test for such efficiency, we use hypothalamic data derived from the Allen Mouse Brain Connectivity Atlas, which provides a structural connectivity map of mouse brain regions derived from a series of viral tracing experiments. Using both cost function minimization and comparison with a weighted, sphere-packing ensemble, we demonstrate that the sum of the distances between hypothalamic sub-regions are not close to the minimum possible distance, consistent with prior whole brain studies. However, if such distances are weighted by the inverse of the magnitude of the connectivity, their sum is among the lowest possible values. Specifically, the hypothalamus appears within the top 94th percentile of neural efficiencies of randomly packed configurations and within one standard deviation of the median efficiency when packings are optimized for maximal neural efficiency. Our results, therefore, indicate that a combination of geometrical and topological constraints help govern the structure of the hypothalamus.

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

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

Structural connectome architecture shapes the maturation of cortical morphology from childhood to adolescence

Cortical thinning is an important hallmark of the maturation of brain morphology during childhood and adolescence. However, the connectome-based wiring mechanism that underlies cortical maturation remains unclear. Here, we show cortical thinning patterns primarily located in the lateral frontal and parietal heteromodal nodes during childhood and adolescence, which are structurally constrained by white matter network architecture and are particularly represented using a network-based diffusion model. Furthermore, connectome-based constraints are regionally heterogeneous, with the largest constraints residing in frontoparietal nodes, and are associated with gene expression signatures of microstructural neurodevelopmental events. These results are highly reproducible in another independent dataset. These findings advance our understanding of network-level mechanisms and the associated genetic basis that underlies the maturational process of cortical morphology during childhood and adolescence.

Predicting evolutionary targets and parameters of gene deletion from expression data

Motivation

Gene deletion is traditionally thought of as a nonadaptive process that removes functional redundancy from genomes, such that it generally receives less attention than duplication in evolutionary turnover studies. Yet, mounting evidence suggests that deletion may promote adaptation via the "less-is-more" evolutionary hypothesis, as it often targets genes harboring unique sequences, expression profiles, and molecular functions. Hence, predicting the relative prevalence of redundant and unique functions among genes targeted by deletion, as well as the parameters underlying their evolution, can shed light on the role of gene deletion in adaptation.

Results

Here, we present CLOUDe, a suite of machine learning methods for predicting evolutionary targets of gene deletion events from expression data. Specifically, CLOUDe models expression evolution as an Ornstein-Uhlenbeck process, and uses multi-layer neural network, extreme gradient boosting, random forest, and support vector machine architectures to predict whether deleted genes are "redundant" or "unique", as well as several parameters underlying their evolution. We show that CLOUDe boasts high power and accuracy in differentiating between classes, and high accuracy and precision in estimating evolutionary parameters, with optimal performance achieved by its neural network architecture. Application of CLOUDe to empirical data from Drosophila suggests that deletion primarily targets genes with unique functions, with further analysis showing these functions to be enriched for protein deubiquitination. Thus, CLOUDe represents a key advance in learning about the role of gene deletion in functional evolution and adaptation.

Availability and implementation

CLOUDe is freely available on GitHub (https://github.com/anddssan/CLOUDe).

Organization of thalamocortical structural covariance and a corresponding 3D atlas of the mouse thalamus

For information from sensory organs to be processed by the brain, it is usually passed to appropriate areas of the cerebral cortex. Almost all of this information passes through the thalamus, a relay structure that reciprocally connects to the vast majority of the cortex. The thalamus facilitates this information transfer through a set of thalamocortical connections that vary in cellular structure, molecular profiles, innervation patterns, and firing rates. Additionally, corticothalamic connections allow for intracortical information transfer through the thalamus. These efferent and afferent connections between the thalamus and cortex have been the focus of many studies, and the importance of cortical connectivity in defining thalamus anatomy is demonstrated by multiple studies that parcellate the thalamus based on cortical connectivity profiles. Here, we examine correlated morphological variation between the thalamus and cortex, or thalamocortical structural covariance. For each voxel in the thalamus as a seed, we construct a cortical structural covariance map that represents correlated cortical volume variation, and examine whether high structural covariance is observed in cortical areas that are functionally relevant to the seed. Then, using these cortical structural covariance maps as features, we subdivide the thalamus into six non-overlapping regions (clusters of voxels), and assess whether cortical structural covariance is associated with cortical connectivity that specifically originates from these regions. We show that cortical structural covariance is high in areas of the cortex that are functionally related to the seed voxel, cortical structural covariance varies along cortical depth, and sharp transitions in cortical structural covariance profiles are observed when varying seed locations in the thalamus. Subdividing the thalamus based on structural covariance, we additionally demonstrate that the six thalamic clusters of voxels stratify cortical structural covariance along the dorsal-ventral, medial-lateral, and anterior-posterior axes. These cluster-associated structural covariance patterns are prominently detected in cortical regions innervated by fibers projecting out of their related thalamic subdivisions. Together, these results advance our understanding of how the thalamus and the cortex couple in their volumes. Our results indicate that these volume correlations reflect functional organization and structural connectivity, and further provides a novel segmentation of the mouse thalamus that can be used to examine thalamic structural variation and thalamocortical structural covariation in disease models.

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

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

Towards a biologically annotated brain connectome

The brain is a network of interleaved neural circuits. In modern connectomics, brain connectivity is typically encoded as a network of nodes and edges, abstracting away the rich biological detail of local neuronal populations. Yet biological annotations for network nodes - such as gene expression, cytoarchitecture, neurotransmitter receptors or intrinsic dynamics - can be readily measured and overlaid on network models. Here we review how connectomes can be represented and analysed as annotated networks. Annotated connectomes allow us to reconceptualize architectural features of networks and to relate the connection patterns of brain regions to their underlying biology. Emerging work demonstrates that annotated connectomes help to make more veridical models of brain network formation, neural dynamics and disease propagation. Finally, annotations can be used to infer entirely new inter-regional relationships and to construct new types of network that complement existing connectome representations. In summary, biologically annotated connectomes offer a compelling way to study neural wiring in concert with local biological features.

The regional variation of laminar thickness in the human isocortex is related to cortical hierarchy and interregional connectivity

The human isocortex consists of tangentially organized layers with unique cytoarchitectural properties. These layers show spatial variations in thickness and cytoarchitecture across the neocortex, which is thought to support function through enabling targeted corticocortical connections. Here, leveraging maps of the 6 cortical layers based on 3D human brain histology, we aimed to quantitatively characterize the systematic covariation of laminar structure in the cortex and its functional consequences. After correcting for the effect of cortical curvature, we identified a spatial pattern of changes in laminar thickness covariance from lateral frontal to posterior occipital regions, which differentiated the dominance of infra- versus supragranular layer thickness. Corresponding to the laminar regularities of cortical connections along cortical hierarchy, the infragranular-dominant pattern of laminar thickness was associated with higher hierarchical positions of regions, mapped based on resting-state effective connectivity in humans and tract-tracing of structural connections in macaques. Moreover, we show that regions with similar laminar thickness patterns have a higher likelihood of structural connections and strength of functional connections. In sum, here we characterize the organization of laminar thickness in the human isocortex and its association with cortico-cortical connectivity, illustrating how laminar organization may provide a foundational principle of cortical function.

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.

NeST: nested hierarchical structure identification in spatial transcriptomic data

Spatial gene expression in tissue is characterized by regions in which particular genes are enriched or depleted. Frequently, these regions contain nested inside them subregions with distinct expression patterns. Segmentation methods in spatial transcriptomic (ST) data extract disjoint regions maximizing similarity over the greatest number of genes, typically on a particular spatial scale, thus lacking the ability to find region-within-region structure. We present NeST, which extracts spatial structure through coexpression hotspots-regions exhibiting localized spatial coexpression of some set of genes. Coexpression hotspots identify structure on any spatial scale, over any possible subset of genes, and are highly explainable. NeST also performs spatial analysis of cell-cell interactions via ligand-receptor, identifying active areas de novo without restriction of cell type or other groupings, in both two and three dimensions. Through application on ST datasets of varying type and resolution, we demonstrate the ability of NeST to reveal a new level of biological structure.

Robust estimation of cortical similarity networks from brain MRI

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

Assortative mixing in micro-architecturally annotated brain connectomes

The wiring of the brain connects micro-architecturally diverse neuronal populations, but the conventional graph model, which encodes macroscale brain connectivity as a network of nodes and edges, abstracts away the rich biological detail of each regional node. Here, we annotate connectomes with multiple biological attributes and formally study assortative mixing in annotated connectomes. Namely, we quantify the tendency for regions to be connected based on the similarity of their micro-architectural attributes. We perform all experiments using four cortico-cortical connectome datasets from three different species, and consider a range of molecular, cellular, and laminar annotations. We show that mixing between micro-architecturally diverse neuronal populations is supported by long-distance connections and find that the arrangement of connections with respect to biological annotations is associated to patterns of regional functional specialization. By bridging scales of cortical organization, from microscale attributes to macroscale connectivity, this work lays the foundation for next-generation annotated connectomics.

Adult expression of Semaphorins and Plexins is essential for motor neuron survival

Axon guidance cues direct the growth and steering of neuronal growth cones, thus guiding the axons to their targets during development. Nonetheless, after axons have reached their targets and established functional circuits, many mature neurons continue to express these developmental cues. The role of axon guidance cues in the adult nervous system has not been fully elucidated. Using the expression pattern data available on FlyBase, we found that more than 96% of the guidance genes that are expressed in the Drosophila melanogaster embryo continue to be expressed in adults. We utilized the GeneSwitch and TARGET systems to spatiotemporally knockdown the expression of these guidance genes selectively in the adult neurons, once the development was completed. We performed an RNA interference (RNAi) screen against 44 guidance genes in the adult Drosophila nervous system and identified 14 genes that are required for adult survival and normal motility. Additionally, we show that adult expression of Semaphorins and Plexins in motor neurons is necessary for neuronal survival, indicating that guidance genes have critical functions in the mature nervous system.

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.

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.

Maximum Entropy Principle Underlies Wiring Length Distribution in Brain Networks

A brain network comprises a substantial amount of short-range connections with an admixture of long-range connections. The portion of long-range connections in brain networks is observed to be quantitatively dissimilar across species. It is hypothesized that the length of connections is constrained by the spatial embedding of brain networks, yet fundamental principles that underlie the wiring length distribution remain unclear. By quantifying the structural diversity of a brain network using Shannon's entropy, here we show that the wiring length distribution across multiple species-including Drosophila, mouse, macaque, human, and C. elegans-follows the maximum entropy principle (MAP) under the constraints of limited wiring material and the spatial locations of brain areas or neurons. In addition, by considering stochastic axonal growth, we propose a network formation process capable of reproducing wiring length distributions of the 5 species, thereby implementing MAP in a biologically plausible manner. We further develop a generative model incorporating MAP, and show that, for the 5 species, the generated network exhibits high similarity to the real network. Our work indicates that the brain connectivity evolves to be structurally diversified by maximizing entropy to support efficient interareal communication, providing a potential organizational principle of brain networks.

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.

Uncovering Statistical Links Between Gene Expression and Structural Connectivity Patterns in the Mouse Brain

Finding links between genes and structural connectivity is of the utmost importance for unravelling the underlying mechanism of the brain connectome. In this study we identify links between the gene expression and the axonal projection density in the mouse brain, by applying a modified version of the Linked ICA method to volumetric data from the Allen Institute for Brain Science for identifying independent sources of information that link both modalities at the voxel level. We performed separate analyses on sets of projections from the visual cortex, the caudoputamen and the midbrain reticular nucleus, and we determined those brain areas, injections and genes that were most involved in independent components that link both gene expression and projection density data, while we validated their biological context through enrichment analysis. We identified representative and literature-validated cortico-midbrain and cortico-striatal projections, whose gene subsets were enriched with annotations for neuronal and synaptic function and related developmental and metabolic processes. The results were highly reproducible when including all available projections, as well as consistent with factorisations obtained using the Dictionary Learning and Sparse Coding technique. Hence, Linked ICA yielded reproducible independent components that were preserved under increasing data variance. Taken together, we have developed and validated a novel paradigm for linking gene expression and structural projection patterns in the mouse mesoconnectome, which can power future studies aiming to relate genes to brain function.

Learning Retention Mechanisms and Evolutionary Parameters of Duplicate Genes from Their Expression Data

Learning about the roles that duplicate genes play in the origins of novel phenotypes requires an understanding of how their functions evolve. A previous method for achieving this goal, CDROM, employs gene expression distances as proxies for functional divergence and then classifies the evolutionary mechanisms retaining duplicate genes from comparisons of these distances in a decision tree framework. However, CDROM does not account for stochastic shifts in gene expression or leverage advances in contemporary statistical learning for performing classification, nor is it capable of predicting the parameters driving duplicate gene evolution. Thus, here we develop CLOUD, a multi-layer neural network built on a model of gene expression evolution that can both classify duplicate gene retention mechanisms and predict their underlying evolutionary parameters. We show that not only is the CLOUD classifier substantially more powerful and accurate than CDROM, but that it also yields accurate parameter predictions, enabling a better understanding of the specific forces driving the evolution and long-term retention of duplicate genes. Further, application of the CLOUD classifier and predictor to empirical data from Drosophila recapitulates many previous findings about gene duplication in this lineage, showing that new functions often emerge rapidly and asymmetrically in younger duplicate gene copies, and that functional divergence is driven by strong natural selection. Hence, CLOUD represents a major advancement in classifying retention mechanisms and predicting evolutionary parameters of duplicate genes, thereby highlighting the utility of incorporating sophisticated statistical learning techniques to address long-standing questions about evolution after gene duplication.

Scaling of gene transcriptional gradients with brain size across mouse development

The structure of the adult brain is the result of complex physical mechanisms acting in three-dimensional space through development. Consequently, the brain's spatial embedding plays a key role in its organization, including the gradient-like patterning of gene expression that encodes the molecular underpinning of functional specialization. However, we do not yet understand how changes in brain shape and size that occur across development influence the brain's transcriptional architecture. Here we investigate the spatial embedding of transcriptional patterns of over 1800 genes across seven time points through mouse-brain development using data from the Allen Developing Mouse Brain Atlas. We find that transcriptional similarity decreases exponentially with separation distance across all developmental time points, with a correlation length scale that follows a power-law scaling relationship with a linear dimension of brain size. This scaling suggests that the mouse brain achieves a characteristic balance between local molecular similarity (homogeneous gene expression within a specialized brain area) and longer-range diversity (between functionally specialized brain areas) throughout its development. Extrapolating this mouse developmental scaling relationship to the human cortex yields a prediction consistent with the value measured from microarray data. We introduce a simple model of brain growth as spatially autocorrelated gene-expression gradients that expand through development, which captures key features of the mouse developmental data. Complementing the well-known exponential distance rule for structural connectivity, our findings characterize an analogous exponential distance rule for transcriptional gradients that scales across mouse brain development, providing new understanding of spatial constraints on the brain's molecular patterning.

Prediction of a Cell-Class-Specific Mouse Mesoconnectome Using Gene Expression Data

Reconstructing brain connectivity at sufficient resolution for computational models designed to study the biophysical mechanisms underlying cognitive processes is extremely challenging. For such a purpose, a mesoconnectome that includes laminar and cell-class specificity would be a major step forward. We analyzed the ability of gene expression patterns to predict cell-class and layer-specific projection patterns and assessed the functional annotations of the most predictive groups of genes. To achieve our goal we used publicly available volumetric gene expression and connectivity data and we trained computational models to learn and predict cell-class and layer-specific axonal projections using gene expression data. Predictions were done in two ways, namely predicting projection strengths using the expression of individual genes and using the co-expression of genes organized in spatial modules, as well as predicting binary forms of projection. For predicting the strength of projections, we found that ridge (L2-regularized) regression had the highest cross-validated accuracy with a median r² score of 0.54 which corresponded for binarized predictions to a median area under the ROC value of 0.89. Next, we identified 200 spatial gene modules using a dictionary learning and sparse coding approach. We found that these modules yielded predictions of comparable accuracy, with a median r² score of 0.51. Finally, a gene ontology enrichment analysis of the most predictive gene groups resulted in significant annotations related to postsynaptic function. Taken together, we have demonstrated a prediction workflow that can be used to perform multimodal data integration to improve the accuracy of the predicted mesoconnectome and support other neuroscience use cases.

Deep Neural Networks for In Situ Hybridization Grid Completion and Clustering

Transcriptome in brain plays a crucial role in understanding the cortical organization and the development of brain structure and function. Two challenges, incomplete data and high dimensionality of transcriptome, remain unsolved. Here, we present a novel training scheme that successfully adapts the U-net architecture to the problem of volume recovery. By analogy to denoising autoencoder, we hide a portion of each training sample so that the network can learn to recover missing voxels from context. Then on the completed volumes, we show that Restricted Boltzmann Machines (RBMs) can be used to infer co-occurrences among voxels, providing foundations for dividing the cortex into discrete subregions. As we stack multiple RBMs to form a deep belief network (DBN), we progressively map the high-dimensional raw input into abstract representations and create a hierarchy of transcriptome architecture. A coarse to fine organization emerges from the network layers. This organization incidentally corresponds to the anatomical structures, suggesting a close link between structures and the genetic underpinnings. Thus, we demonstrate a new way of learning transcriptome-based hierarchical organization using RBM and DBN.

Comparison of brain connectomes by MRI and genomics and its implication in Alzheimer's disease

BACKGROUND

The human brain is complex and interconnected structurally. Brain connectome change is associated with Alzheimer's disease (AD) and other neurodegenerative diseases. Genetics and genomics studies have identified molecular changes in AD; however, the results are often limited to isolated brain regions and are difficult to interpret its findings in respect to brain connectome. The mechanisms of how one brain region impacts the molecular pathways in other regions have not been systematically studied. And how the brain regions susceptible to AD pathology interact with each other at the transcriptome level and how these interactions relate to brain connectome change are unclear.

METHODS

Here, we compared structural brain connectomes defined by probabilistic tracts using diffusion magnetic resonance imaging data in Alzheimer's Disease Neuroimaging Initiative database and a brain transcriptome dataset covering 17 brain regions.

RESULTS

We observed that the changes in diffusion measures associated with AD diagnosis status and the associations were replicated in an independent cohort. The result suggests that disease associated white matter changes are focal. Analysis of the brain connectome by genomic data, tissue-tissue transcriptional synchronization between 17 brain regions, indicates that the regions connected by AD-associated tracts were likely connected at the transcriptome level with high number of tissue-to-tissue correlated (TTC) gene pairs (P = 0.03). And genes involved in TTC gene pairs between white matter tract connected brain regions were enriched in signaling pathways (P = 6.08 × 10-9). Further pathway interaction analysis identified ionotropic glutamate receptor pathway and Toll receptor signaling pathways to be important for tissue-tissue synchronization at the transcriptome level. Transcript profile entailing Toll receptor signaling in the blood was significantly associated with diffusion properties of white matter tracts, notable association between fractional anisotropy and bilateral cingulum angular bundles (Ppermutation = 1.0 × 10-2 and 4.9 × 10-4 for left and right respectively).

CONCLUSIONS

In summary, our study suggests that brain connectomes defined by MRI and transcriptome data overlap with each other.

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.

Altered Neocortical Gene Expression, Brain Overgrowth and Functional Over-Connectivity in Chd8 Haploinsufficient Mice

Truncating CHD8 mutations are amongst the highest confidence risk factors for autism spectrum disorder (ASD) identified to date. Here, we report that Chd8 heterozygous mice display increased brain size, motor delay, hypertelorism, pronounced hypoactivity, and anomalous responses to social stimuli. Whereas gene expression in the neocortex is only mildly affected at midgestation, over 600 genes are differentially expressed in the early postnatal neocortex. Genes involved in cell adhesion and axon guidance are particularly prominent amongst the downregulated transcripts. Resting-state functional MRI identified increased synchronized activity in cortico-hippocampal and auditory-parietal networks in Chd8 heterozygous mutant mice, implicating altered connectivity as a potential mechanism underlying the behavioral phenotypes. Together, these data suggest that altered brain growth and diminished expression of important neurodevelopmental genes that regulate long-range brain wiring are followed by distinctive anomalies in functional brain connectivity in Chd8+/- mice. Human imaging studies have reported altered functional connectivity in ASD patients, with long-range under-connectivity seemingly more frequent. Our data suggest that CHD8 haploinsufficiency represents a specific subtype of ASD where neuropsychiatric symptoms are underpinned by long-range over-connectivity.

Structural Insights Into the Dynamic Evolution of Neuronal Networks as Synaptic Density Decreases

The human brain is thought to be an extremely complex but efficient computing engine, processing vast amounts of information from a changing world. The decline in the synaptic density of neuronal networks is one of the most important characteristics of brain development, which is closely related to synaptic pruning, synaptic growth, synaptic plasticity, and energy metabolism. However, because of technical limitations in observing large-scale neuronal networks dynamically connected through synapses, how neuronal networks are organized and evolve as their synaptic density declines remains unclear. Here, by establishing a biologically reasonable neuronal network model, we show that despite a decline in the synaptic density, the connectivity, and efficiency of neuronal networks can be improved. Importantly, by analyzing the degree distribution, we also find that both the scale-free characteristic of neuronal networks and the emergence of hub neurons rely on the spatial distance between neurons. These findings may promote our understanding of neuronal networks in the brain and have guiding significance for the design of neuronal network models.

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.

Development of Tetrahydroindazole-Based Potent and Selective Sigma-2 Receptor Ligands

The sigma-2 receptor has been shown to play important roles in a number of important diseases, including central nervous system (CNS) disorders and cancer. However, mechanisms by which sigma-2 contributes to these diseases remain unclear. The development of new sigma-2 ligands that can be used to probe the function of this protein and potentially as drug discovery leads is therefore of great importance. Herein we report the development of a series of tetrahydroindazole compounds that are highly potent and selective for sigma-2. Structure-activity relationship data were used to generate a pharmacophore model that summarizes the common features present in the potent ligands. Assays for solubility and microsomal stability showed that several members of this compound series possess promising characteristics for further development of useful chemical probes or drug discovery leads.

Human visual cortex is organized along two genetically opposed hierarchical gradients with unique developmental and evolutionary origins

Human visual cortex is organized with striking consistency across individuals. While recent findings demonstrate an unexpected coupling between functional and cytoarchitectonic regions relative to the folding of human visual cortex, a unifying principle linking these anatomical and functional features of the cortex remains elusive. To fill this gap in knowledge, we combined independent and ground truth measurements of cytoarchitectonic regions and genetic tissue characterization within human occipitotemporal cortex. Using a data-driven approach, we examined whether differential gene expression among cytoarchitectonic areas could contribute to the arealization of occipitotemporal cortex into a hierarchy based on transcriptomics. This approach revealed two opposing gene expression gradients: one that contains a series of genes with expression magnitudes that ascend from posterior (e.g., areas human occipital [hOc]1, hOc2, hOc3, etc.) to anterior cytoarchitectonic areas (e.g., areas fusiform gyrus [FG]1–FG4) and another that contains a separate series of genes that show a descending gradient from posterior to anterior areas. Using data from the living human brain, we show that each of these gradients correlates strongly with variations in measures related to either thickness or myelination of cortex, respectively. We further reveal that these genetic gradients emerge along unique trajectories in human development: the ascending gradient is present at 10–12 gestational weeks, while the descending gradient emerges later (19–24 gestational weeks). Interestingly, it is not until early childhood (before 5 years of age) that the two expression gradients achieve their adult-like mean expression values. Additional analyses in nonhuman primates (NHPs) reveal that homologous genes do not generate the same ascending and descending expression gradients as in humans. We discuss these findings relative to previously proposed hierarchies based on functional and cytoarchitectonic features of visual cortex. Altogether, these findings bridge macroscopic features of human cytoarchitectonic areas in visual cortex with microscopic features of cellular organization and genetic expression, which, despite the complexity of this multiscale correspondence, can be described by a sparse subset (approximately 200) of genes. These findings help pinpoint the genes contributing to healthy cortical development and explicate the cortical biology distinguishing humans from other primates, as well as establishing essential groundwork for understanding future work linking genetic mutations with the function and development of the human brain.

The expression of a sparse subset of human genes forms two opposed gradients that capture the processing hierarchy of visual cortex; these transcription gradients emerge at different points during human development and distinguish human from nonhuman primates.

Exploiting Gene Expression Profiles for the Automated Prediction of Connectivity between Brain Regions

The brain comprises a complex system of neurons interconnected by an intricate network of anatomical links. While recent studies demonstrated the correlation between anatomical connectivity patterns and gene expression of neurons, using transcriptomic information to automatically predict such patterns is still an open challenge. In this work, we present a completely data-driven approach relying on machine learning (i.e., neural networks) to learn the anatomical connection directly from a training set of gene expression data. To do so, we combined gene expression and connectivity data from the Allen Mouse Brain Atlas to generate thousands of gene expression profile pairs from different brain regions. To each pair, we assigned a label describing the physical connection between the corresponding brain regions. Then, we exploited these data to train neural networks, designed to predict brain area connectivity. We assessed our solution on two prediction problems (with three and two connectivity class categories) involving cortical and cerebellum regions. As demonstrated by our results, we distinguish between connected and unconnected regions with 85% prediction accuracy and good balance of precision and recall. In our future work we may extend the analysis to more complex brain structures and consider RNA-Seq data as additional input to our model.

Neuronal Activities in the Mouse Visual Cortex Predict Patterns of Sensory Stimuli

Visual cortex forms the basis of visual processing and plays important roles in visual encoding. By using the recently published Allen Brain Observatory dataset consisting of large-scale calcium imaging of mouse V1 activities under visual stimuli, we were able to obtain high-quality data capturing simultaneous neuronal activities at multiple sub-areas and cortical depths of V1. Using prediction models, we analyzed the activity profiles related to static and drifting grating stimuli. We conducted a comprehensive survey of the coding ability of multiple cortical locations toward different stimulus attributes. Specifically, we focused on orientations and spatial frequencies (for static stimuli), as well as moving directions and speed (for drifting stimuli). By using results produced from a prediction model, we quantified the decoding performance profile at different sub-areas and layers of V1. In addition, we analyzed the interactions and interference between different stimulus attributes. The insights obtained from these discoveries would contribute to more precise and quantitative understanding of V1 coding mechanisms.

Lineage-Specific Expression Divergence in Grasses Is Associated with Male Reproduction, Host-Pathogen Defense, and Domestication

Poaceae (grasses) is an agriculturally important and widely distributed family of plants with extraordinary phenotypic diversity, much of which was generated under recent lineage-specific evolution. Yet, little is known about the genes and functional modules involved in the lineage-specific divergence of grasses. Here, I address this question on a genome-wide scale by applying a novel branch-based statistic of lineage-specific expression divergence, LED, to RNA-seq data from nine tissues of the wild grass Brachypodium distachyon and its domesticated relatives Oryza sativa japonica (rice) and Sorghum bicolor (sorghum). I find that LED is generally smallest in B. distachyon and largest in O. sativa japonica, which underwent domestication earlier than S. bicolor, supporting the hypothesis that domestication may increase the rate of lineage-specific expression divergence in grasses. Moreover, in all three species, LED is positively correlated with protein-coding sequence divergence and tissue specificity, and negatively correlated with network connectivity. Further analysis reveals that genes with large LED are often primarily expressed in anther, implicating lineage-specific expression divergence in the evolution of male reproductive phenotypes. Gene ontology enrichment analysis also identifies an overrepresentation of terms related to male reproduction in the two domesticated grasses, as well as to those involved in host-pathogen defense in all three species. Last, examinations of genes with the largest LED reveal that their lineage-specific expression divergence may have contributed to antimicrobial functions in B. distachyon, to enhanced adaptation and yield during domestication in O. sativa japonica, and to defense against a widespread and devastating fungal pathogen in S. bicolor. Together, these findings suggest that lineage-specific expression divergence in grasses may increase under domestication and preferentially target rapidly evolving genes involved in male reproduction, host-pathogen defense, and the origin of domesticated phenotypes.

Spatial Embedding Imposes Constraints on Neuronal Network Architectures

Recent progress towards understanding circuit function has capitalized on tools from network science to parsimoniously describe the spatiotemporal architecture of neural systems. Such tools often address systems topology divorced from its physical instantiation. Nevertheless, for embedded systems such as the brain, physical laws directly constrain the processes of network growth, development, and function. We review here the rules imposed by the space and volume of the brain on the development of neuronal networks, and show that these rules give rise to a specific set of complex topologies. These rules also affect the repertoire of neural dynamics that can emerge from the system, and thereby inform our understanding of network dysfunction in disease. We close by discussing new tools and models to delineate the effects of spatial embedding.

Different Behavioral Experiences Produce Distinctive Parallel Changes in, and Correlate With, Frontal Cortex and Hippocampal Global Post-translational Histone Levels

While it is clear that behavioral experience modulates epigenetic profiles, it is less evident how the nature of that experience influences outcomes and whether epigenetic/genetic "biomarkers" could be extracted to classify different types of behavioral experience. To begin to address this question, male and female mice were subjected to either a Fixed Interval (FI) schedule of food reward, or a single episode of forced swim followed by restraint stress, or no explicit behavioral experience after which global expression levels of two activating (H3K9ac and H3K4me3) and two repressive (H3K9me2 and H3k27me3) post-translational histone modifications (PTHMs), were measured in hippocampus (HIPP) and frontal cortex (FC). The specific nature of the behavioral experience differentiated profiles of PTHMs in a sex- and brain region-dependent manner, with all 4 PTHMs changing in parallel in response to different behavioral experiences. These different behavioral experiences also modified the pattern of correlations of PTHMs both within and across FC and HIPP. Unexpectedly, highly robust correlations were found between global PTHM levels and behavioral performances, suggesting that global PTHMs may provide a higher-order pattern recognition function. Further efforts are needed to determine the generality of such findings and what characteristics of behavioral experience are critical for modulating PTHM responses.

Structural covariance of brain region volumes is associated with both structural connectivity and transcriptomic similarity

An organizational pattern seen in the brain, termed structural covariance, is the statistical association of pairs of brain regions in their anatomical properties. These associations, measured across a population as covariances or correlations usually in cortical thickness or volume, are thought to reflect genetic and environmental underpinnings. Here, we examine the biological basis of structural volume covariance in the mouse brain. We first examined large scale associations between brain region volumes using an atlas-based approach that parcellated the entire mouse brain into 318 regions over which correlations in volume were assessed, for volumes obtained from 153 mouse brain images via high-resolution MRI. We then used a seed-based approach and determined, for 108 different seed regions across the brain and using mouse gene expression and connectivity data from the Allen Institute for Brain Science, the variation in structural covariance data that could be explained by distance to seed, transcriptomic similarity to seed, and connectivity to seed. We found that overall, correlations in structure volumes hierarchically clustered into distinct anatomical systems, similar to findings from other studies and similar to other types of networks in the brain, including structural connectivity and transcriptomic similarity networks. Across seeds, this structural covariance was significantly explained by distance (17% of the variation, up to a maximum of 49% for structural covariance to the visceral area of the cortex), transcriptomic similarity (13% of the variation, up to maximum of 28% for structural covariance to the primary visual area) and connectivity (15% of the variation, up to a maximum of 36% for structural covariance to the intermediate reticular nucleus in the medulla) of covarying structures. Together, distance, connectivity, and transcriptomic similarity explained 37% of structural covariance, up to a maximum of 63% for structural covariance to the visceral area. Additionally, this pattern of explained variation differed spatially across the brain, with transcriptomic similarity playing a larger role in the cortex than subcortex, while connectivity explains structural covariance best in parts of the cortex, midbrain, and hindbrain. These results suggest that both gene expression and connectivity underlie structural volume covariance, albeit to different extents depending on brain region, and this relationship is modulated by distance.

Multimodal Connectomics in Psychiatry: Bridging Scales From Micro to Macro

The human brain is a highly complex system, with a large variety of microscale cellular morphologies and macroscale global properties. Working at multiple scales, it forms an efficient system for processing and integration of multimodal information. Studies have repeatedly demonstrated strong associations between modalities of both microscales and macroscales of brain organization. These consistent observations point toward potential common organization principles where regions with a microscale architecture supportive of a larger computational load have more and stronger connections in the brain network on the macroscale. Conversely, disruptions observed on one organizational scale could modulate the other. First neuropsychiatric micro-macro comparisons in, among other conditions, Alzheimer's disease and schizophrenia, have, for example, shown overlapping alterations across both scales. We give an overview of recent findings on associations between microscale and macroscale organization observed in the healthy brain, followed by a summary of microscale and macroscale findings reported in the context of brain disorders. We conclude with suggestions for future multiscale connectome comparisons linking multiple scales and modalities of organization and suggest how such comparisons could contribute to a more complete fundamental understanding of brain organization and associated disease-related alterations.

Gene expression links functional networks across cortex and striatum

The human brain is comprised of a complex web of functional networks that link anatomically distinct regions. However, the biological mechanisms supporting network organization remain elusive, particularly across cortical and subcortical territories with vastly divergent cellular and molecular properties. Here, using human and primate brain transcriptional atlases, we demonstrate that spatial patterns of gene expression show strong correspondence with limbic and somato/motor cortico-striatal functional networks. Network-associated expression is consistent across independent human datasets and evolutionarily conserved in non-human primates. Genes preferentially expressed within the limbic network (encompassing nucleus accumbens, orbital/ventromedial prefrontal cortex, and temporal pole) relate to risk for psychiatric illness, chloride channel complexes, and markers of somatostatin neurons. Somato/motor associated genes are enriched for oligodendrocytes and markers of parvalbumin neurons. These analyses indicate that parallel cortico-striatal processing channels possess dissociable genetic signatures that recapitulate distributed functional networks, and nominate molecular mechanisms supporting cortico-striatal circuitry in health and disease.

The functional connectivity of brain regions can be reflected in a shared molecular architecture. This cross-modal study demonstrates correspondence of spatial patterns of gene expression to limbic and somato/motor cortico-striatal networks in human and non-human primates.

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.

The Pathoconnectivity Profile of Alzheimer's Disease: A Morphometric Coalteration Network Analysis

Gray matter alterations are typical features of brain disorders. However, they do not impact on the brain randomly. Indeed, it has been suggested that neuropathological processes can selectively affect certain assemblies of neurons, which typically are at the center of crucial functional networks. Because of their topological centrality, these areas form a core set that is more likely to be affected by neuropathological processes. In order to identify and study the pattern formed by brain alterations in patients’ with Alzheimer’s disease (AD), we devised an innovative meta-analytic method for analyzing voxel-based morphometry data. This methodology enabled us to discover that in AD gray matter alterations do not occur randomly across the brain but, on the contrary, follow identifiable patterns of distribution. This alteration pattern exhibits a network-like structure composed of coaltered areas that can be defined as coatrophy network. Within the coatrophy network of AD, we were able to further identify a core subnetwork of coaltered areas that includes the left hippocampus, left and right amygdalae, right parahippocampal gyrus, and right temporal inferior gyrus. In virtue of their network centrality, these brain areas can be thought of as pathoconnectivity hubs.

Heritable genotype contrast mining reveals novel gene associations specific to autism subgroups

Though the genetic etiology of autism is complex, our understanding can be improved by identifying genes and gene-gene interactions that contribute to the development of specific autism subtypes. Identifying such gene groupings will allow individuals to be diagnosed and treated according to their precise characteristics. To this end, we developed a method to associate gene combinations with groups with shared autism traits, targeting genetic elements that distinguish patient populations with opposing phenotypes. Our computational method prioritizes genetic variants for genome-wide association, then utilizes Frequent Pattern Mining to highlight potential interactions between variants. We introduce a novel genotype assessment metric, the Unique Inherited Combination support, which accounts for inheritance patterns observed in the nuclear family while estimating the impact of genetic variation on phenotype manifestation at the individual level. High-contrast variant combinations are tested for significant subgroup associations. We apply this method by contrasting autism subgroups defined by severe or mild manifestations of a phenotype. Significant associations connected 286 genes to the subgroups, including 193 novel autism candidates. 71 pairs of genes have joint associations with subgroups, presenting opportunities to investigate interacting functions. This study analyzed 12 autism subgroups, but our informatics method can explore other meaningful divisions of autism patients, and can further be applied to reveal precise genetic associations within other phenotypically heterogeneous disorders, such as Alzheimer's disease.

Natural Selection Drives Rapid Functional Evolution of Young Drosophila Duplicate Genes

Gene duplication is thought to play a major role in phenotypic evolution. Yet the forces involved in the functional divergence of young duplicate genes remain unclear. Here, we use population-genetic inference to elucidate the role of natural selection in the functional evolution of young duplicate genes in Drosophila melanogaster. We find that negative selection acts on young duplicates with ancestral functions, and positive selection on those with novel functions, suggesting that natural selection may determine whether and how young duplicate genes are retained. Moreover, evidence of natural selection is strongest in protein-coding regions and 3' UTRs of young duplicates, indicating that selection may primarily target encoded proteins and regulatory sequences specific to 3' UTRs. Further analysis reveals that natural selection acts immediately after duplication and weakens over time, possibly explaining the observed bias toward the acquisition of new functions by young, rather than old, duplicate gene copies. Last, we find an enrichment of testis-related functions in young duplicates that underwent recent positive selection, but not in young duplicates that did not undergo recent positive selection, or in old duplicates that either did or did not undergo recent positive selection. Thus, our findings reveal that natural selection is a key player in the functional evolution of young duplicate genes, acts rapidly and in a region-specific manner, and may underlie the origin of novel testis-specific phenotypes in Drosophila.

Gene transcription profiles associated with inter-modular hubs and connection distance in human functional magnetic resonance imaging networks

Human functional magnetic resonance imaging (fMRI) brain networks have a complex topology comprising integrative components, e.g. long-distance inter-modular edges, that are theoretically associated with higher biological cost. Here, we estimated intra-modular degree, inter-modular degree and connection distance for each of 285 cortical nodes in multi-echo fMRI data from 38 healthy adults. We used the multivariate technique of partial least squares (PLS) to reduce the dimensionality of the relationships between these three nodal network parameters and prior microarray data on regional expression of 20 737 genes. The first PLS component defined a transcriptional profile associated with high intra-modular degree and short connection distance, whereas the second PLS component was associated with high inter-modular degree and long connection distance. Nodes in superior and lateral cortex with high inter-modular degree and long connection distance had local transcriptional profiles enriched for oxidative metabolism and mitochondria, and for genes specific to supragranular layers of human cortex. In contrast, primary and secondary sensory cortical nodes in posterior cortex with high intra-modular degree and short connection distance had transcriptional profiles enriched for RNA translation and nuclear components. We conclude that, as predicted, topologically integrative hubs, mediating long-distance connections between modules, are more costly in terms of mitochondrial glucose metabolism.This article is part of the themed issue 'Interpreting BOLD: a dialogue between cognitive and cellular neuroscience'.

Predicting functional neuroanatomical maps from fusing brain networks with genetic information

Functional neuroanatomical maps provide a mesoscale reference framework for studies from molecular to systems neuroscience and psychiatry. The underlying structure-function relationships are typically derived from functional manipulations or imaging approaches. Although highly informative, these are experimentally costly. The increasing amount of publicly available brain and genetic data offers a rich source that could be mined to address this problem computationally. Here, we developed an algorithm that fuses gene expression and connectivity data with functional genetic meta data and exploits cumulative effects to derive neuroanatomical maps related to multi-genic functions. We validated the approach by using public available mouse and human data. The generated neuroanatomical maps recapture known functional anatomical annotations from literature and functional MRI data. When applied to multi-genic meta data from mouse quantitative trait loci (QTL) studies and human neuropsychiatric databases, this method predicted known functional maps underlying behavioral or psychiatric traits. Taken together, genetically weighted connectivity analysis (GWCA) allows for high throughput functional exploration of brain anatomy in silico. It maps functional genetic associations onto brain circuitry for refining functional neuroanatomy, or identifying trait-associated brain circuitry, from genetic data.

Commentary: BRAIN NETWORKS. Correlated Gene Expression Supports Synchronous Activity in Brain Networks. Science 348, 1241-4

A recent report claims that functional brain networks defined with resting-state functional magnetic resonance imaging (fMRI) can be recapitulated with correlated gene expression (i.e., high within-network tissue-tissue "strength fraction," SF) (Richiardi et al., 2015). However, the authors do not adequately control for spatial proximity. We replicated their main analysis, performed a more effective adjustment for spatial proximity, and tested whether "null networks" (i.e., clusters with center coordinates randomly placed throughout cortex) also exhibit high SF. Removing proximal tissue-tissue correlations by Euclidean distance, as opposed to removing correlations within arbitrary tissue labels as in Richiardi et al. (2015), reduces within-network SF to no greater than null. Moreover, randomly placed clusters also have significantly high SF, indicating that high within-network SF is entirely attributable to proximity and is unrelated to functional brain networks defined by resting-state fMRI. We discuss why additional validations in the original article are invalid and/or misleading and suggest future directions.

Small molecule modulator of sigma 2 receptor is neuroprotective and reduces cognitive deficits and neuroinflammation in experimental models of Alzheimer's disease

Accumulating evidence suggests that modulating the sigma 2 receptor (Sig2R) can provide beneficial effects for neurodegenerative diseases. Herein, we report the identification of a novel class of Sig2R ligands and their cellular and in vivo activity in experimental models of Alzheimer's disease (AD). We report that SAS-0132 and DKR-1051, selective ligands of Sig2R, modulate intracellular Ca²⁺ levels in human SK-N-SH neuroblastoma cells. The Sig2R ligands SAS-0132 and JVW-1009 are neuroprotective in a C. elegans model of amyloid precursor protein-mediated neurodegeneration. Since this neuroprotective effect is replicated by genetic knockdown and knockout of vem-1, the ortholog of progesterone receptor membrane component-1 (PGRMC1), these results suggest that Sig2R ligands modulate a PGRMC1-related pathway. Last, we demonstrate that SAS-0132 improves cognitive performance both in the Thy-1 hAPP^Lond/Swe+ transgenic mouse model of AD and in healthy wild-type mice. These results demonstrate that Sig2R is a promising therapeutic target for neurocognitive disorders including AD.

Transcriptomic Perspectives on Neocortical Structure, Development, Evolution, and Disease

The cerebral cortex is the source of our most complex cognitive capabilities and a vulnerable target of many neurological and neuropsychiatric disorders. Transcriptomics offers a new approach to understanding the cortex at the level of its underlying genetic code, and rapid technological advances have propelled this field to the high-throughput study of the complete set of transcribed genes at increasingly fine resolution to the level of individual cells. These tools have revealed features of the genetic architecture of adult cortical areas, layers, and cell types, as well as spatiotemporal patterning during development. This has allowed a fresh look at comparative anatomy as well, illustrating surprisingly large differences between mammals while at the same time revealing conservation of some features from avians to mammals. Finally, transcriptomics is fueling progress in understanding the causes of neurodevelopmental diseases such as autism, linking genetic association studies to specific molecular pathways and affected brain regions.

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.

Network component analysis reveals developmental trajectories of structural connectivity and specific alterations in autism spectrum disorder

The structural organization of the brain can be characterized as a hierarchical ensemble of segregated modules linked by densely interconnected hub regions that facilitate distributed functional interactions. Disturbances to this network may be an important marker of abnormal development. Recently, several neurodevelopmental disorders, including autism spectrum disorder (ASD), have been framed as disorders of connectivity but the full nature and timing of these disturbances remain unclear. In this study, we use non-negative matrix factorization, a data-driven, multivariate approach, to model the structural network architecture of the brain as a set of superposed subnetworks, or network components. In an openly available dataset of 196 subjects scanned between 5 and 85 years we identify a set of robust and reliable subnetworks that develop in tandem with age and reflect both anatomically local and long-range, network hub connections. In a second experiment, we compare network components in a cohort of 51 high-functioning ASD adolescents to a group of age-matched controls. We identify a specific subnetwork representing an increase in local connection strength in the cingulate cortex in ASD (t = 3.44, P < 0.001). This work highlights possible long-term implications of alterations to the developmental trajectories of specific cortical subnetworks. Hum Brain Mapp 38:4169-4184, 2017. © 2017 Wiley Periodicals, Inc.