Gene expression of Caenorhabditis elegans neurons carries information on their synaptic connectivity

Deciphering the genetic code of neuronal type connectivity through bilinear modeling

Understanding how different neuronal types connect and communicate is critical to interpreting brain function and behavior. However, it has remained a formidable challenge to decipher the genetic underpinnings that dictate the specific connections formed between neuronal types. To address this, we propose a novel bilinear modeling approach that leverages the architecture similar to that of recommendation systems. Our model transforms the gene expressions of presynaptic and postsynaptic neuronal types, obtained from single-cell transcriptomics, into a covariance matrix. The objective is to construct this covariance matrix that closely mirrors a connectivity matrix, derived from connectomic data, reflecting the known anatomical connections between these neuronal types. When tested on a dataset of Caenorhabditis elegans, our model achieved a performance comparable to, if slightly better than, the previously proposed spatial connectome model (SCM) in reconstructing electrical synaptic connectivity based on gene expressions. Through a comparative analysis, our model not only captured all genetic interactions identified by the SCM but also inferred additional ones. Applied to a mouse retinal neuronal dataset, the bilinear model successfully recapitulated recognized connectivity motifs between bipolar cells and retinal ganglion cells, and provided interpretable insights into genetic interactions shaping the connectivity. Specifically, it identified unique genetic signatures associated with different connectivity motifs, including genes important to cell-cell adhesion and synapse formation, highlighting their role in orchestrating specific synaptic connections between these neurons. Our work establishes an innovative computational strategy for decoding the genetic programming of neuronal type connectivity. It not only sets a new benchmark for single-cell transcriptomic analysis of synaptic connections but also paves the way for mechanistic studies of neural circuit assembly and genetic manipulation of circuit wiring.

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.

Connectomes: from a sparsity of networks to large-scale databases

The analysis of whole brain networks started in the 1980s when only a handful of connectomes were available. In these early days, information about the human connectome was absent and one could only dream about having information about connectivity in a single human subject. Thanks to non-invasive methods such as diffusion imaging, we now know about connectivity in many species and, for some species, in many individuals. To illustrate the rapid change in availability of connectome data, the UK Biobank is on track to record structural and functional connectivity in 100,000 human subjects. Moreover, connectome data from a range of species is now available: from Caenorhabditis elegans and the fruit fly to pigeons, rodents, cats, non-human primates, and humans. This review will give a brief overview of what structural connectivity data is now available, how connectomes are organized, and how their organization shows common features across species. Finally, I will outline some of the current challenges and potential future work in making use of connectome information.

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.

Molecular topography of an entire nervous system

We have produced gene expression profiles of all 302 neurons of the C. elegans nervous system that match the single-cell resolution of its anatomy and wiring diagram. Our results suggest that individual neuron classes can be solely identified by combinatorial expression of specific gene families. For example, each neuron class expresses distinct codes of ∼23 neuropeptide genes and ∼36 neuropeptide receptors, delineating a complex and expansive "wireless" signaling network. To demonstrate the utility of this comprehensive gene expression catalog, we used computational approaches to (1) identify cis-regulatory elements for neuron-specific gene expression and (2) reveal adhesion proteins with potential roles in process placement and synaptic specificity. Our expression data are available at https://cengen.org and can be interrogated at the web application CengenApp. We expect that this neuron-specific directory of gene expression will spur investigations of underlying mechanisms that define anatomy, connectivity, and function throughout the C. elegans nervous system.

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.

Uncovering the genetic blueprint of the C. elegans nervous system

Despite rapid advances in connectome mapping and neuronal genetics, we lack theoretical and computational tools to unveil, in an experimentally testable fashion, the genetic mechanisms that govern neuronal wiring. Here we introduce a computational framework to link the adjacency matrix of a connectome to the expression patterns of its neurons, helping us uncover a set of genetic rules that govern the interactions between neurons in contact. The method incorporates the biological realities of the system, accounting for noise from data collection limitations, as well as spatial restrictions. The resulting methodology allows us to infer a network of 19 innexin interactions that govern the formation of gap junctions in Caenorhabditis elegans, five of which are already supported by experimental data. As advances in single-cell gene expression profiling increase the accuracy and the coverage of the data, the developed framework will allow researchers to systematically infer experimentally testable connection rules, offering mechanistic predictions for synapse and gap junction formation.

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.

A Genetic Model of the Connectome

The connectomes of organisms of the same species show remarkable architectural and often local wiring similarity, raising the question: where and how is neuronal connectivity encoded? Here, we start from the hypothesis that the genetic identity of neurons guides synapse and gap-junction formation and show that such genetically driven wiring predicts the existence of specific biclique motifs in the connectome. We identify a family of large, statistically significant biclique subgraphs in the connectomes of three species and show that within many of the observed bicliques the neurons share statistically significant expression patterns and morphological characteristics, supporting our expectation of common genetic factors that drive the synapse formation within these subgraphs. The proposed connectome model offers a self-consistent framework to link the genetics of an organism to the reproducible architecture of its connectome, offering experimentally falsifiable predictions on the genetic factors that drive the formation of individual neuronal circuits.

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.

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.

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.

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.

A transcriptional signature of hub connectivity in the mouse connectome

Connectivity is not distributed evenly throughout the brain. Instead, it is concentrated on a small number of highly connected neural elements that act as network hubs. Across different species and measurement scales, these hubs show dense interconnectivity, forming a core or "rich club" that integrates information across anatomically distributed neural systems. Here, we show that projections between connectivity hubs of the mouse brain are both central (i.e., they play an important role in neural communication) and costly (i.e., they extend over long anatomical distances) aspects of network organization that carry a distinctive genetic signature. Analyzing the neuronal connectivity of 213 brain regions and the transcriptional coupling, across 17,642 genes, between each pair of regions, we find that coupling is highest for pairs of connected hubs, intermediate for links between hubs and nonhubs, and lowest for connected pairs of nonhubs. The high transcriptional coupling associated with hub connectivity is driven by genes regulating the oxidative synthesis and metabolism of ATP--the primary energetic currency of neuronal communication. This genetic signature contrasts that identified for neuronal connectivity in general, which is driven by genes regulating neuronal, synaptic, and axonal structure and function. Our findings establish a direct link between molecular function and the large-scale topology of neuronal connectivity, showing that brain hubs display a tight coordination of gene expression, often over long anatomical distances, that is intimately related to the metabolic requirements of these highly active network elements.

From Caenorhabditis elegans to the human connectome: a specific modular organization increases metabolic, functional and developmental efficiency

The connectome, or the entire connectivity of a neural system represented by a network, ranges across various scales from synaptic connections between individual neurons to fibre tract connections between brain regions. Although the modularity they commonly show has been extensively studied, it is unclear whether the connection specificity of such networks can already be fully explained by the modularity alone. To answer this question, we study two networks, the neuronal network of Caenorhabditis elegans and the fibre tract network of human brains obtained through diffusion spectrum imaging. We compare them to their respective benchmark networks with varying modularities, which are generated by link swapping to have desired modularity values. We find several network properties that are specific to the neural networks and cannot be fully explained by the modularity alone. First, the clustering coefficient and the characteristic path length of both C. elegans and human connectomes are higher than those of the benchmark networks with similar modularity. High clustering coefficient indicates efficient local information distribution, and high characteristic path length suggests reduced global integration. Second, the total wiring length is smaller than for the alternative configurations with similar modularity. This is due to lower dispersion of connections, which means each neuron in the C. elegans connectome or each region of interest in the human connectome reaches fewer ganglia or cortical areas, respectively. Third, both neural networks show lower algorithmic entropy compared with the alternative arrangements. This implies that fewer genes are needed to encode for the organization of neural systems. While the first two findings show that the neural topologies are efficient in information processing, this suggests that they are also efficient from a developmental point of view. Together, these results show that neural systems are organized in such a way as to yield efficient features beyond those given by their modularity alone.

Developmental time windows for axon growth influence neuronal network topology

Early brain connectivity development consists of multiple stages: birth of neurons, their migration and the subsequent growth of axons and dendrites. Each stage occurs within a certain period of time depending on types of neurons and cortical layers. Forming synapses between neurons either by growing axons starting at similar times for all neurons (much-overlapped time windows) or at different time points (less-overlapped) may affect the topological and spatial properties of neuronal networks. Here, we explore the extreme cases of axon formation during early development, either starting at the same time for all neurons (parallel, i.e., maximally overlapped time windows) or occurring for each neuron separately one neuron after another (serial, i.e., no overlaps in time windows). For both cases, the number of potential and established synapses remained comparable. Topological and spatial properties, however, differed: Neurons that started axon growth early on in serial growth achieved higher out-degrees, higher local efficiency and longer axon lengths while neurons demonstrated more homogeneous connectivity patterns for parallel growth. Second, connection probability decreased more rapidly with distance between neurons for parallel growth than for serial growth. Third, bidirectional connections were more numerous for parallel growth. Finally, we tested our predictions with C. elegans data. Together, this indicates that time windows for axon growth influence the topological and spatial properties of neuronal networks opening up the possibility to a posteriori estimate developmental mechanisms based on network properties of a developed network.

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.

High-resolution prediction of mouse brain connectivity using gene expression patterns

The brain is a multi-level system in which the high-level functions are generated by low-level genetic mechanisms. Thus, elucidating the relationship among multiple brain levels via correlative and predictive analytics is an important area in brain research. Currently, studies in multiple species have indicated that the spatiotemporal gene expression patterns are predictive of brain wiring. Specifically, results on the worm Caenorhabditis elegans have shown that the prediction of neuronal connectivity using gene expression signatures yielded statistically significant results. Recent studies on the mammalian brain produced similar results at the coarse regional level. In this study, we provide the first high-resolution, large-scale integrative analysis of the transcriptome and connectome in a single mammalian brain at a fine voxel level. By using the Allen Brain Atlas data, we predict voxel-level brain connectivity based on the gene expressions in the adult mouse brain. We employ regularized models to show that gene expression is predictive of connectivity at the voxel-level with an accuracy of 93%. We also identify a set of genes playing the most important role in connectivity prediction. We use only this small number of genes to predict the brain wiring with an accuracy over 80%. We discover that these important genes are enriched in neurons as compared to glia, and they perform connectivity-related functions. We perform several interesting correlative studies to further elucidate the transcriptome-connectome relationship.

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

Unraveling the relationship between molecular signatures in the brain and their functional, architectonic, and anatomic correlates is an important neuroscientific goal. It is still not well understood whether the diversity demonstrated by histological studies in the human brain is reflected in the spatial patterning of whole brain transcriptional profiles. Using genome-wide maps of transcriptional distribution of the human brain by the Allen Brain Institute, we test the hypothesis that gene expression profiles are specific to anatomically described brain regions. In this work, we demonstrate that this is indeed the case by showing that gene similarity clusters appear to respect conventional basal-cortical and caudal-rostral gradients. To fully investigate the causes of this observed spatial clustering, we test a connectionist hypothesis that states that the spatial patterning of gene expression in the brain is simply reflective of the fiber tract connectivity between brain regions. We find that although gene expression and structural connectivity are not determined by each other, they do influence each other with a high statistical significance. This implies that spatial diversity of gene expressions is a result of mainly location-specific features but is influenced by neuronal connectivity, such that like cellular species preferentially connects with like cells.

Coupling of gene expression in medial prefrontal cortex and nucleus accumbens after neonatal ventral hippocampal lesions accompanies deficits in sensorimotor gating and auditory processing in rats

BACKGROUND

After neonatal ventral hippocampal lesions (NVHLs), adult rats exhibit evidence of neural processing deficits relevant to schizophrenia, including reduced prepulse inhibition (PPI) of acoustic startle and impaired sensory processing. In intact rats, the regulation of PPI by the ventral hippocampus (VH) is mediated via interactions with medial prefrontal cortex (mPFC) and nucleus accumbens (NAC). We assessed PPI, auditory-evoked responses and expression of 7 schizophrenia-related genes in mPFC and NAC, in adult rats after sham- or real NVHLs.

METHODS

Male inbred Buffalo (BUF) rat pups (d7; n=36) received either vehicle or ibotenic acid infusion into the VH. PPI and auditory-evoked dentate gyrus local field potentials (LFPs) were measured on d56 and d66, respectively. Brains were processed for RT-PCR measures of mPFC and NAC Comt, Erbb4, Grid2, Ncam1, Slc1a2, Nrg1 and Reln.

RESULTS

NVHL rats exhibited significant deficits in PPI (p=0.005) and LFPs (p<0.015) proportional to lesion size. Sham vs. NVHL rats did not differ in gene expression levels in mPFC or NAC. As we previously reported, multiple gene expression levels were highly correlated within- (mean r's≈0.5), but not across-brain regions (mean r's≈0). However, for three genes--Comt, Slc1a2 and Ncam1--after NVHLs, expression levels became significantly correlated, or "coupled," across the mPFC and NAC (p's<0.03, 0.002 and 0.05, respectively), and the degree of "coupling" increased with VH lesion size.

CONCLUSIONS

After NVHLs that disrupt PPI and auditory processing, specific gene expression levels suggest an abnormal functional coupling of the mPFC and NAC. This model of VH-mPFC-NAC network dysfunction after NVHLs may have implications for understanding the neural basis for PPI- and related sensory processing deficits in schizophrenia patients.

Large-Scale Analysis of Gene Expression and Connectivity in the Rodent Brain: Insights through Data Integration

Recent research in C. elegans and the rodent has identified correlations between gene expression and connectivity. Here we extend this type of approach to examine complex patterns of gene expression in the rodent brain in the context of regional brain connectivity and differences in cellular populations. Using multiple large-scale data sets obtained from public sources, we identified two novel patterns of mouse brain gene expression showing a strong degree of anti-correlation, and relate this to multiple data modalities including macroscale connectivity. We found that these signatures are associated with differences in expression of neuronal and oligodendrocyte markers, suggesting they reflect regional differences in cellular populations. We also find that the expression level of these genes is correlated with connectivity degree, with regions expressing the neuron-enriched pattern having more incoming and outgoing connections with other regions. Our results exemplify what is possible when increasingly detailed large-scale cell- and gene-level data sets are integrated with connectivity data.

Gene expression in the rodent brain is associated with its regional connectivity

The putative link between gene expression of brain regions and their neural connectivity patterns is a fundamental question in neuroscience. Here this question is addressed in the first large scale study of a prototypical mammalian rodent brain, using a combination of rat brain regional connectivity data with gene expression of the mouse brain. Remarkably, even though this study uses data from two different rodent species (due to the data limitations), we still find that the connectivity of the majority of brain regions is highly predictable from their gene expression levels-the outgoing (incoming) connectivity is successfully predicted for 73% (56%) of brain regions, with an overall fairly marked accuracy level of 0.79 (0.83). Many genes are found to play a part in predicting both the incoming and outgoing connectivity (241 out of the 500 top selected genes, p-value<1e-5). Reassuringly, the genes previously known from the literature to be involved in axon guidance do carry significant information about regional brain connectivity. Surveying the genes known to be associated with the pathogenesis of several brain disorders, we find that those associated with schizophrenia, autism and attention deficit disorder are the most highly enriched in the connectivity-related genes identified here. Finally, we find that the profile of functional annotation groups that are associated with regional connectivity in the rodent is significantly correlated with the annotation profile of genes previously found to determine neural connectivity in C. elegans (Pearson correlation of 0.24, p<1e-6 for the outgoing connections and 0.27, p<1e-5 for the incoming). Overall, the association between connectivity and gene expression in a specific extant rodent species' brain is likely to be even stronger than found here, given the limitations of current data.

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

We studied the global relationship between gene expression and neuroanatomical connectivity in the adult rodent brain. We utilized a large data set of the rat brain “connectome” from the Brain Architecture Management System (942 brain regions and over 5000 connections) and used statistical approaches to relate the data to the gene expression signatures of 17,530 genes in 142 anatomical regions from the Allen Brain Atlas. Our analysis shows that adult gene expression signatures have a statistically significant relationship to connectivity. In particular, brain regions that have similar expression profiles tend to have similar connectivity profiles, and this effect is not entirely attributable to spatial correlations. In addition, brain regions which are connected have more similar expression patterns. Using a simple optimization approach, we identified a set of genes most correlated with neuroanatomical connectivity, and find that this set is enriched for genes involved in neuronal development and axon guidance. A number of the genes have been implicated in neurodevelopmental disorders such as autistic spectrum disorder. Our results have the potential to shed light on the role of gene expression patterns in influencing neuronal activity and connectivity, with potential applications to our understanding of brain disorders. Supplementary data are available at http://www.chibi.ubc.ca/ABAMS.

We tested the idea that the “wiring diagram” of the adult brain has a relationship with where genes are expressed. We were inspired by similar work carried out by groups examining the nematode worm Caenorhabditis elegans. By using large-scale databases of brain connectivity and gene expression in rodents, we found that many genes involved in the development of the brain show correlations with anatomical connectivity patterns. Some of the genes we found have been implicated in disorders such as autism, which is suspected to affect brain wiring. While the biological causes of the patterns we found are not yet known, we believe they provide new insight into the patterns of gene expression in the brain and will spur further study of this problem.

Random distance dependent attachment as a model for neural network generation in the Caenorhabditis elegans

MOTIVATION

The topology of the network induced by the neurons connectivity's in the Caenorhabditis elegans differs from most common random networks. The neurons positions of the C.elegans have been previously explained as being optimal to induce the required network wiring. We here propose a complementary explanation that the network wiring is the direct result of a local stochastic synapse formation process.

RESULTS

We show that a model based on the physical distance between neurons can explain the C.elegans neural network structure, specifically, we demonstrate that a simple model based on a geometrical synapse formation probability and the inhibition of short coherent cycles can explain the properties of the C.elegans' neural network. We suggest this model as an initial framework to discuss neural network generation and as a first step toward the development of models for more advanced creatures. In order to measure the circle frequency in the network, a novel graph-theory circle length measurement algorithm is proposed.

Navigating the cell: UNC-53 and the navigators, a family of cytoskeletal regulators with multiple roles in cell migration, outgrowth and trafficking

Changes in cell shape are associated with a variety of processes including cell migration, axon outgrowth, cell division and vesicle trafficking. C. elegans UNC-53 and its vertebrate homologs, the Navigators, are required for the migration of cells and the outgrowth of neuronal processes. The identification of novel molecular interactions and live imaging studies have revealed that UNC-53/NAVs are signal transducers associated with actin filaments, microtubules and intermediate filaments. In addition to modulating cytoskeletal dynamics at the leading edge of migrating or outgrowing cells, both UNC-53 and the navigators are expressed in adult cells, conspicuously those with specialized roles in endocytosis or secretion. Collectively, these results suggest that UNC-53/NAVs may be a central regulator of cytoskeletal dynamics, responsible for integrating signaling cues to multiple components of the cytoskeleton to coordinate rearrangement during cell outgrowth or trafficking.

Large-scale reconstruction and phylogenetic analysis of metabolic environments

The topology of metabolic networks may provide important insights not only into the metabolic capacity of species, but also into the habitats in which they evolved. Here we introduce the concept of a metabolic network's "seed set"--the set of compounds that, based on the network topology, are exogenously acquired--and provide a methodological framework to computationally infer the seed set of a given network. Such seed sets form ecological "interfaces" between metabolic networks and their surroundings, approximating the effective biochemical environment of each species. Analyzing the metabolic networks of 478 species and identifying the seed set of each species, we present a comprehensive large-scale reconstruction of such predicted metabolic environments. The seed sets' composition significantly correlates with several basic properties characterizing the species' environments and agrees with biological observations concerning major adaptations. Species whose environments are highly predictable (e.g., obligate parasites) tend to have smaller seed sets than species living in variable environments. Phylogenetic analysis of the seed sets reveals the complex dynamics governing gain and loss of seeds across the phylogenetic tree and the process of transition between seed and non-seed compounds. Our findings suggest that the seed state is transient and that seeds tend either to be dropped completely from the network or to become non-seed compounds relatively fast. The seed sets also permit a successful reconstruction of a phylogenetic tree of life. The "reverse ecology" approach presented lays the foundations for studying the evolutionary interplay between organisms and their habitats on a large scale.

Using expression profiles of Caenorhabditis elegans neurons to identify genes that mediate synaptic connectivity

Synaptic wiring of neurons in Caenorhabditis elegans is largely invariable between animals. It has been suggested that this feature stems from genetically encoded molecular markers that guide the neurons in the final stage of synaptic formation. Identifying these markers and unraveling the logic by which they direct synapse formation is a key challenge. Here, we address this task by constructing a probabilistic model that attempts to explain the neuronal connectivity diagram of C. elegans as a function of the expression patterns of its neurons. By only considering neuron pairs that are known to be connected by chemical or electrical synapses, we focus on the final stage of synapse formation, in which neurons identify their designated partners. Our results show that for many neurons the neuronal expression map of C. elegans can be used to accurately predict the subset of adjacent neurons that will be chosen as its postsynaptic partners. Notably, these predictions can be achieved using the expression patterns of only a small number of specific genes that interact in a combinatorial fashion.

Geometric constraints on neuronal connectivity facilitate a concise synaptic adhesive code

The nervous system contains trillions of neurons, each forming thousands of synaptic connections. It has been suggested that this complex connectivity is determined by a synaptic "adhesive code," where connections are dictated by a variable set of cell surface proteins, combinations of which form neuronal addresses. The estimated number of neuronal addresses is orders of magnitude smaller than the number of neurons. Here, we show that the limited number of addresses dictates constraints on the possible neuronal network topologies. We show that to encode arbitrary networks, in which each neuron can potentially connect to any other neuron, the number of neuronal addresses needed scales linearly with network size. In contrast, the number of addresses needed to encode the wiring of geometric networks grows only as the square root of network size. The more efficient encoding in geometric networks is achieved through the reutilization of the same addresses in physically independent portions of the network. We also find that ordered geometric networks, in which the same connectivity patterns are iterated throughout the network, further reduce the required number of addresses. We demonstrate our findings using simulated networks and the C. elegans neuronal network. Geometric neuronal connectivity with recurring connectivity patterns have been suggested to confer an evolutionary advantage by saving biochemical resources on the one hand and reutilizing functionally efficient neuronal circuits. Our study suggests an additional advantage of these prominent topological features--the facilitation of the ability to genetically encode neuronal networks given constraints on the number of addresses.

The extracellular leucine-rich repeat superfamily; a comparative survey and analysis of evolutionary relationships and expression patterns

BACKGROUND

Leucine-rich repeats (LRRs) are highly versatile and evolvable protein-ligand interaction motifs found in a large number of proteins with diverse functions, including innate immunity and nervous system development. Here we catalogue all of the extracellular LRR (eLRR) proteins in worms, flies, mice and humans. We use convergent evidence from several transmembrane-prediction and motif-detection programs, including a customised algorithm, LRRscan, to identify eLRR proteins, and a hierarchical clustering method based on TribeMCL to establish their evolutionary relationships.

RESULTS

This yields a total of 369 proteins (29 in worm, 66 in fly, 135 in mouse and 139 in human), many of them of unknown function. We group eLRR proteins into several classes: those with only LRRs, those that cluster with Toll-like receptors (Tlrs), those with immunoglobulin or fibronectin-type 3 (FN3) domains and those with some other domain. These groups show differential patterns of expansion and diversification across species. Our analyses reveal several clusters of novel genes, including two Elfn genes, encoding transmembrane proteins with eLRRs and an FN3 domain, and six genes encoding transmembrane proteins with eLRRs only (the Elron cluster). Many of these are expressed in discrete patterns in the developing mouse brain, notably in the thalamus and cortex. We have also identified a number of novel fly eLRR proteins with discrete expression in the embryonic nervous system.

CONCLUSION

This study provides the necessary foundation for a systematic analysis of the functions of this class of genes, which are likely to include prominently innate immunity, inflammation and neural development, especially the specification of neuronal connectivity.

Distributed representations accelerate evolution of adaptive behaviours

Animals with rudimentary innate abilities require substantial learning to transform those abilities into useful skills, where a skill can be considered as a set of sensory–motor associations. Using linear neural network models, it is proved that if skills are stored as distributed representations, then within-lifetime learning of part of a skill can induce automatic learning of the remaining parts of that skill. More importantly, it is shown that this “free-lunch” learning (FLL) is responsible for accelerated evolution of skills, when compared with networks which either 1) cannot benefit from FLL or 2) cannot learn. Specifically, it is shown that FLL accelerates the appearance of adaptive behaviour, both in its innate form and as FLL-induced behaviour, and that FLL can accelerate the rate at which learned behaviours become innate.

Some behaviours are purely innate (e.g., blinking), whereas other, “apparently innate,” behaviours require a degree of learning to refine them into a useful skill (e.g., nest building). In terms of biological fitness, it matters how quickly such learning occurs, because time spent learning is time spent not eating, or time spent being eaten, both of which reduce fitness. Using artificial neural networks as model organisms, it is proven that it is possible for an organism to be born with a set of “primed” connections which guarantee that learning part of a skill induces automatic learning of other skill components, an effect known as free-lunch learning (FLL). Critically, this effect depends on the assumption that associations are stored as distributed representations. Using a genetic algorithm, it is shown that primed organisms can evolve within 30 generations. This has three important consequences. First, primed organisms learn quickly, which increases their fitness. Second, the presence of FLL effectively accelerates the rate of evolution, for both learned and innate skill components. Third, FLL can accelerate the rate at which learned behaviours become innate. These findings suggest that species may depend on the presence of distributed representations to ensure rapid evolution of adaptive behaviours.

Predicting the connectivity of primate cortical networks from topological and spatial node properties

Background

The organization of the connectivity between mammalian cortical areas has become a major subject of study, because of its important role in scaffolding the macroscopic aspects of animal behavior and intelligence. In this study we present a computational reconstruction approach to the problem of network organization, by considering the topological and spatial features of each area in the primate cerebral cortex as subsidy for the reconstruction of the global cortical network connectivity. Starting with all areas being disconnected, pairs of areas with similar sets of features are linked together, in an attempt to recover the original network structure.

Results

Inferring primate cortical connectivity from the properties of the nodes, remarkably good reconstructions of the global network organization could be obtained, with the topological features allowing slightly superior accuracy to the spatial ones. Analogous reconstruction attempts for the C. elegans neuronal network resulted in substantially poorer recovery, indicating that cortical area interconnections are relatively stronger related to the considered topological and spatial properties than neuronal projections in the nematode.

Conclusion

The close relationship between area-based features and global connectivity may hint on developmental rules and constraints for cortical networks. Particularly, differences between the predictions from topological and spatial properties, together with the poorer recovery resulting from spatial properties, indicate that the organization of cortical networks is not entirely determined by spatial constraints.

Cell-specific microarray profiling experiments reveal a comprehensive picture of gene expression in the C. elegans nervous system

BACKGROUND

With its fully sequenced genome and simple, well-defined nervous system, the nematode Caenorhabditis elegans offers a unique opportunity to correlate gene expression with neuronal differentiation. The lineal origin, cellular morphology and synaptic connectivity of each of the 302 neurons are known. In many instances, specific behaviors can be attributed to particular neurons or circuits. Here we describe microarray-based methods that monitor gene expression in C. elegans neurons and, thereby, link comprehensive profiles of neuronal transcription to key developmental and functional properties of the nervous system.

RESULTS

We employed complementary microarray-based strategies to profile gene expression in the embryonic and larval nervous systems. In the MAPCeL (Microarray Profiling C. elegans cells) method, we used fluorescence activated cell sorting (FACS) to isolate GFP-tagged embryonic neurons for microarray analysis. To profile the larval nervous system, we used the mRNA-tagging technique in which an epitope-labeled mRNA binding protein (FLAG-PAB-1) was transgenically expressed in neurons for immunoprecipitation of cell-specific transcripts. These combined approaches identified approximately 2,500 mRNAs that are highly enriched in either the embryonic or larval C. elegans nervous system. These data are validated in part by the detection of gene classes (for example, transcription factors, ion channels, synaptic vesicle components) with established roles in neuronal development or function. Of particular interest are 19 conserved transcripts of unknown function that are also expressed in the mammalian brain. In addition to utilizing these profiling approaches to define stage-specific gene expression, we also applied the mRNA-tagging method to fingerprint a specific neuron type, the A-class group of cholinergic motor neurons, during early larval development. A comparison of these data to a MAPCeL profile of embryonic A-class motor neurons identified genes with common functions in both types of A-class motor neurons as well as transcripts with roles specific to each motor neuron type.

CONCLUSION

We describe microarray-based strategies for generating expression profiles of embryonic and larval C. elegans neurons. These methods can be applied to particular neurons at specific developmental stages and, therefore, provide an unprecedented opportunity to obtain spatially and temporally defined snapshots of gene expression in a simple model nervous system.