Local signal time-series during rest used for areal boundary mapping in individual human brains

It is widely thought that resting state functional connectivity likely reflects functional interaction among brain areas and that different functional areas interact with different sets of brain areas. A method for mapping areal boundaries has been formulated based on the large-scale spatial characteristics of regional interaction revealed by resting state functional connectivity. In the present study, we present a novel analysis for areal boundary mapping that requires only the signal timecourses within a region of interest, without reference to the information from outside the region. The areal boundaries were generated by the novel analysis and were compared with those generated by the previously-established standard analysis. The boundaries were robust and reproducible across the two analyses, in two regions of interest tested. These results suggest that the information for areal boundaries is readily available inside the region of interest.     

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.     

On variability in the density of corticocortical and thalamocortical connections

Variability is an important but neglected aspect of connectional neuroanatomy. The quantitative density of the 'same' corticocortical or thalamocortical connection may vary by over two orders of magnitude between different injections of the same tracer. At present, however, the frequency distribution of connection densities is unknown. Therefore, it is unclear what kind of sampling strategies or statistical methods are appropriate for quantitative studies of connectivity. Nor is it clear if the measured variability represents differences between subjects, or if it is simply a consequence of intra-individual differences resulting from experimental technique and the exact placement of tracers relative to local spatial and laminar variation in connectivity. We used quantitative measurements of the density of a large number of corticocortical and thalamocortical connections from our own laboratories and from the literature. Variability in the density of given corticocortical and thalamocortical connections is high, with the standard deviation of density proportional to the mean. The frequency distribution is close to exponential. Therefore, analysis methods relying on the normal distribution are not appropriate. We provide an appendix that gives simple statistical guidance for samples drawn from exponentially distributed data. For a given corticocortical or thalamocortical connection density, between-individual standard deviation is 0.85 to 1.25 times the within-individual standard deviation. Therefore, much of the variability reported in conventional neuroanatomical studies (with one tracer deposited per animal) is due to within-individual factors. We also find that strong, but not weak, corticocortical connections are substantially more variable than thalamocortical connections. We propose that the near exponential distribution of connection densities is a simple consequence of 'patchy' connectivity. We anticipate that connection data will be well described by the negative binomial, a class of distribution that applies to events occurring in clumped or patchy substrates. Local patchiness may be a feature of all corticocortical connections and could explain why strong corticocortical connections are more variable than strong thalamocortical connections. This idea is supported by the columnar patterns of many corticocortical but few thalamocortical connections in the literature.     

Small-world connectivity, motif composition, and complexity of fractal neuronal connections

Connection patterns of the cerebral cortex consist of pathways linking neuronal populations across multiple levels of scale, from whole brain regions to local minicolumns. This nested interconnectivity suggests the hypothesis that cortical connections are arranged in fractal or self-similar patterns. We describe a simple procedure to generate fractal connection patterns that aim at capturing the potential self-similarity and hierarchical ordering of neuronal connections. We examine these connection patterns by calculating a broad range of structural measures, including small-world attributes and motif composition, as well as some global measures of functional connectivity, including complexity. As we vary fractal patterns by changing a critical control parameter, we find strongly correlated changes in several structural and functional measures, suggesting that they emerge together and are mutually linked. Measures obtained from some modeled fractal patterns closely resemble those of real neuroanatomical data sets, supporting the original hypothesis.     

Combining structural connectivity and response latencies to model the structure of the visual system

Several approaches exist to ascertain the connectivity of the brain, and these approaches lead to markedly different topologies, often incompatible with each other. Specifically, recent single-cell recording results seem incompatible with current structural connectivity models. We present a novel method that combines anatomical and temporal constraints to generate biologically plausible connectivity patterns of the visual system of the macaque monkey. Our method takes structural connectivity data from the CoCoMac database and recent single-cell recording data as input and employs an optimization technique to arrive at a new connectivity pattern of the visual system that is in agreement with both types of experimental data. The new connectivity pattern yields a revised model that has fewer levels than current models. In addition, it introduces subcortical-cortical connections. We show that these connections are essential for explaining latency data, are consistent with our current knowledge of the structural connectivity of the visual system, and might explain recent functional imaging results in humans. Furthermore we show that the revised model is not underconstrained like previous models and can be extended to include newer data and other kinds of data. We conclude that the revised model of the connectivity of the visual system reflects current knowledge on the structure and function of the visual system and addresses some of the limitations of previous models.     

Human connectomics

Recent advances in non-invasive neuroimaging have enabled the measurement of connections between distant regions in the living human brain, thus opening up a new field of research: Human connectomics. Different imaging modalities allow the mapping of structural connections (axonal fibre tracts) as well as functional connections (correlations in time series), and individual variations in these connections may be related to individual variations in behaviour and cognition. Connectivity analysis has already led to a number of new insights about brain organization. For example, segregated brain regions may be identified by their unique patterns of connectivity, structural and functional connectivity may be compared to elucidate how dynamic interactions arise from the anatomical substrate, and the architecture of large-scale networks connecting sets of brain regions may be analysed in detail. The combined analysis of structural and functional networks has begun to reveal components or modules with distinct patterns of connections that become engaged in different cognitive tasks. Collectively, advances in human connectomics open up the possibility of studying how brain connections mediate regional brain function and thence behaviour.     

Estimation of Directed Effective Connectivity from fMRI Functional Connectivity Hints at Asymmetries of Cortical Connectome

The brain exhibits complex spatio-temporal patterns of activity. This phenomenon is governed by an interplay between the internal neural dynamics of cortical areas and their connectivity. Uncovering this complex relationship has raised much interest, both for theory and the interpretation of experimental data (e.g., fMRI recordings) using dynamical models. Here we focus on the so-called inverse problem: the inference of network parameters in a cortical model to reproduce empirically observed activity. Although it has received a lot of interest, recovering directed connectivity for large networks has been rather unsuccessful so far. The present study specifically addresses this point for a noise-diffusion network model. We develop a Lyapunov optimization that iteratively tunes the network connectivity in order to reproduce second-order moments of the node activity, or functional connectivity. We show theoretically and numerically that the use of covariances with both zero and non-zero time shifts is the key to infer directed connectivity. The first main theoretical finding is that an accurate estimation of the underlying network connectivity requires that the time shift for covariances is matched with the time constant of the dynamical system. In addition to the network connectivity, we also adjust the intrinsic noise received by each network node. The framework is applied to experimental fMRI data recorded for subjects at rest. Diffusion-weighted MRI data provide an estimate of anatomical connections, which is incorporated to constrain the cortical model. The empirical covariance structure is reproduced faithfully, especially its temporal component (i.e., time-shifted covariances) in addition to the spatial component that is usually the focus of studies. We find that the cortical interactions, referred to as effective connectivity, in the tuned model are not reciprocal. In particular, hubs are either receptors or feeders: they do not exhibit both strong incoming and outgoing connections. Our results sets a quantitative ground to explore the propagation of activity in the cortex.     

A Statistical Method of Identifying Interactions in Neuron–Glia Systems Based on Functional Multicell Ca2+ Imaging

Crosstalk between neurons and glia may constitute a significant part of information processing in the brain. We present a novel method of statistically identifying interactions in a neuron–glia network. We attempted to identify neuron–glia interactions from neuronal and glial activities via maximum-a-posteriori (MAP)-based parameter estimation by developing a generalized linear model (GLM) of a neuron–glia network. The interactions in our interest included functional connectivity and response functions. We evaluated the cross-validated likelihood of GLMs that resulted from the addition or removal of connections to confirm the existence of specific neuron-to-glia or glia-to-neuron connections. We only accepted addition or removal when the modification improved the cross-validated likelihood. We applied the method to a high-throughput, multicellular in vitro Ca2+ imaging dataset obtained from the CA3 region of a rat hippocampus, and then evaluated the reliability of connectivity estimates using a statistical test based on a surrogate method. Our findings based on the estimated connectivity were in good agreement with currently available physiological knowledge, suggesting our method can elucidate undiscovered functions of neuron–glia systems.     

Functional Brain Networks: Random, “Small World” or Deterministic?

Lately the problem of connectivity in brain networks is being approached frequently by graph theoretical analysis. In several publications based on bivariate estimators of relations between EEG channels authors reported random or “small world” structure of networks. The results of these works often have no relation to other evidence based on imaging, inverse solutions methods, physiological and anatomical data. Herein we try to find reasons for this discrepancy. We point out that EEG signals are very much interdependent, thus bivariate measures applied to them may produce many spurious connections. In fact, they may outnumber the true connections. Giving all connections equal weights, as it is usual in the framework of graph theoretical analysis, further enhances these spurious links. In effect, close to random and disorganized patterns of connections emerge. On the other hand, multivariate connectivity estimators, which are free of the artificial links, show specific, well determined patterns, which are in a very good agreement with other evidence. The modular structure of brain networks may be identified by multivariate estimators based on Granger causality and formalism of assortative mixing. In this way, the strength of coupling may be evaluated quantitatively. During working memory task, by means of multivariate Directed Transfer Function, it was demonstrated that the modules characterized by strong internal bonds exchange the information by weaker connections.     

Role of structural inhomogeneities in resting-state brain dynamics

Brain imaging methods allow a non-invasive assessment of both structural and functional connectivity. However, the mechanism of how functional connectivity arises in a structured network of interacting neural populations is as yet poorly understood. Here we use a modeling approach to explore the way in which functional correlations arise from underlying structural connections taking into account inhomogeneities in the interactions between the brain regions of interest. The local dynamics of a neural population is assumed to be of phase-oscillator type. The considered structural connectivity patterns describe long-range anatomical connections between interacting neural elements. We find a dependence of the simulated functional connectivity patterns on the parameters governing the dynamics. We calculate graph-theoretic measures of the functional network topology obtained from numerical simulations. The effect of structural inhomogeneities in the coupling term on the observed network state is quantified by examining the relation between simulated and empirical functional connectivity. Importantly, we show that simulated and empirical functional connectivity agree for a narrow range of coupling strengths. We conclude that identification of functional connectivity during rest requires an analysis of the network dynamics.     

Time-variant estimation of directed influences during Parkinsonian tremor

The inference of interaction structures in multidimensional time series is a major challenge not only in neuroscience but in many fields of research. To gather information about the connectivity in a network from measured data, several parametric as well as non-parametric approaches have been proposed and widely examined. Today a lot of interest is focused on the evolution of the network connectivity in time which might contain information about ongoing tasks in the brain or possible dynamic dysfunctions. Therefore an extension of the current approaches towards time-resolved analysis techniques is desired. We present a parametric approach for time variant analysis, test its performance for simulated data, and apply it to real-world data.     

The wiring economy principle: connectivity determines anatomy in the human brain

Minimization of the wiring cost of white matter fibers in the human brain appears to be an organizational principle. We investigate this aspect in the human brain using whole brain connectivity networks extracted from high resolution diffusion MRI data of 14 normal volunteers. We specifically address the question of whether brain anatomy determines its connectivity or vice versa. Unlike previous studies we use weighted networks, where connections between cortical nodes are real-valued rather than binary off-on connections. In one set of analyses we found that the connectivity structure of the brain has near optimal wiring cost compared to random networks with the same number of edges, degree distribution and edge weight distribution. A specifically designed minimization routine could not find cheaper wiring without significantly degrading network performance. In another set of analyses we kept the observed brain network topology and connectivity but allowed nodes to freely move on a 3D manifold topologically identical to the brain. An efficient minimization routine was written to find the lowest wiring cost configuration. We found that beginning from any random configuration, the nodes invariably arrange themselves in a configuration with a striking resemblance to the brain. This confirms the widely held but poorly tested claim that wiring economy is a driving principle of the brain. Intriguingly, our results also suggest that the brain mainly optimizes for the most desirable network connectivity, and the observed brain anatomy is merely a result of this optimization.     

Maximizing Negative Correlations in Resting-State Functional Connectivity MRI by Time-Lag

This paper aims to better understand the physiological meaning of negative correlations in resting state functional connectivity MRI (r-fcMRI). The correlations between anatomy-based brain regions of 18 healthy humans were calculated and analyzed with and without a correction for global signal and with and without spatial smoothing. In addition, correlations between anatomy-based brain regions of 18 naïve anesthetized rats were calculated and compared to the human data. T-statistics were used to differentiate between positive and negative connections. The application of spatial smoothing and global signal correction increased the number of significant positive connections but their effect on negative connections was complex. Positive connections were mainly observed between cortical structures while most negative connections were observed between cortical and non-cortical structures with almost no negative connections between non-cortical structures. In both human and rats, negative connections were never observed between bilateral homologous regions. The main difference between positive and negative connections in both the human and rat data was that positive connections became less significant with time-lags, while negative connections became more significant with time-lag. This effect was evident in all four types of analyses (with and without global signal correction and spatial smoothing) but was most significant in the analysis with no correction for the global signal. We hypothesize that the valence of r-fcMRI connectivity reflects the relative contributions of cerebral blood volume (CBV) and flow (CBF) to the BOLD signal and that these relative contributions are location-specific. If cerebral circulation is primarily regulated by CBF in one region and by CBV in another, a functional connection between these regions can manifest as an r-fcMRI negative and time-delayed correlation. Similarly, negative correlations could result from spatially inhomogeneous responses of rCBV or rCBF alone. Consequently, neuronal regulation of brain circulation may be deduced from the valence of r-fcMRI connectivity.     

Disrupted functional brain connectivity in partial epilepsy: a resting-state fMRI study

Examining the spontaneous activity to understand the neural mechanism of brain disorder is a focus in recent resting-state fMRI. In the current study, to investigate the alteration of brain functional connectivity in partial epilepsy in a systematical way, two levels of analyses (functional connectivity analysis within resting state networks (RSNs) and functional network connectivity (FNC) analysis) were carried out on resting-state fMRI data acquired from the 30 participants including 14 healthy controls(HC) and 16 partial epilepsy patients. According to the etiology, all patients are subdivided into temporal lobe epilepsy group (TLE, included 7 patients) and mixed partial epilepsy group (MPE, 9 patients). Using group independent component analysis, eight RSNs were identified, and selected to evaluate functional connectivity and FNC between groups. Compared with the controls, decreased functional connectivity within all RSNs was found in both TLE and MPE. However, dissociating patterns were observed within the 8 RSNs between two patient groups, i.e, compared with TLE, we found decreased functional connectivity in 5 RSNs increased functional connectivity in 1 RSN, and no difference in the other 2 RSNs in MPE. Furthermore, the hierarchical disconnections of FNC was found in two patient groups, in which the intra-system connections were preserved for all three subsystems while the lost connections were confined to intersystem connections in patients with partial epilepsy. These findings may suggest that decreased resting state functional connectivity and disconnection of FNC are two remarkable characteristics of partial epilepsy. The selective impairment of FNC implicated that it is unsuitable to understand the partial epilepsy only from global or local perspective. We presumed that studying epilepsy in the multi-perspective based on RSNs may be a valuable means to assess the functional changes corresponding to specific RSN and may contribute to the understanding of the neuro-pathophysiological mechanism of epilepsy.     

Neurotree: A Collaborative,Graphical Database of the Academic Genealogy of Neuroscience

Neurotree is an online database that documents the lineage of academic mentorship in neuroscience. Modeled on the tree format typically used to describe biological genealogies, the Neurotree web site provides a concise summary of the intellectual history of neuroscience and relationships between individuals in the current neuroscience community. The contents of the database are entirely crowd-sourced: any internet user can add information about researchers and the connections between them. As of July 2012, Neurotree has collected information from 10,000 users about 35,000 researchers and 50,000 mentor relationships, and continues to grow. The present report serves to highlight the utility of Neurotree as a resource for academic research and to summarize some basic analysis of its data. The tree structure of the database permits a variety of graphical analyses. We find that the connectivity and graphical distance between researchers entered into Neurotree early has stabilized and thus appears to be mostly complete. The connectivity of more recent entries continues to mature. A ranking of researcher fecundity based on their mentorship reveals a sustained period of influential researchers from 1850–1950, with the most influential individuals active at the later end of that period. Finally, a clustering analysis reveals that some subfields of neuroscience are reflected in tightly interconnected mentor-trainee groups.     

How Does Connectivity Between Cortical Areas Depend on Brain Size? Implications for Efficient Computation

A formula for an average connectivity between cortical areas in mammals is derived. Based on comparative neuroanatomical data, it is found, surprisingly, that this connectivity is either only weakly dependent or independent of brain size. It is discussed how this formula can be used to estimate the average length of axons in white matter. Other allometric relations, such as cortical patches and area sizes vs. brain size, are also provided. Finally, some functional implications, with an emphasis on efficient cortical computation, are discussed as well.     

The Human PAX6 Mutation Database.

The Human PAX6 Mutation Database contains details of 94 mutations of the PAX6 gene. A Microsoft Access program is used by the Curator to store, update and search the database entries. Mutations can be entered directly by the Curator, or imported from submissions made via the World Wide Web. The PAX6 Mutation Database web page at URL http://www.hgu.mrc.ac.uk/Softdata/PAX6/ provides information about PAX6, as well as a fill-in form through which new mutations can be submitted to the Curator. A search facility allows remote users to query the database. A plain text format file of the data can be downloaded via the World Wide Web. The Curation program contains prior knowledge of the genetic code and of the PAX6 gene including cDNA sequence, location of intron/exon boundaries, and protein domains, so that the minimum of information need be provided by the submitter or Curator.     

Statistics of weighted brain networks reveal hierarchical organization and Gaussian degree distribution

Whole brain weighted connectivity networks were extracted from high resolution diffusion MRI data of 14 healthy volunteers. A statistically robust technique was proposed for the removal of questionable connections. Unlike most previous studies our methods are completely adapted for networks with arbitrary weights. Conventional statistics of these weighted networks were computed and found to be comparable to existing reports. After a robust fitting procedure using multiple parametric distributions it was found that the weighted node degree of our networks is best described by the normal distribution, in contrast to previous reports which have proposed heavy tailed distributions. We show that post-processing of the connectivity weights, such as thresholding, can influence the weighted degree asymptotics. The clustering coefficients were found to be distributed either as gamma or power-law distribution, depending on the formula used. We proposed a new hierarchical graph clustering approach, which revealed that the brain network is divided into a regular base-2 hierarchical tree. Connections within and across this hierarchy were found to be uncommonly ordered. The combined weight of our results supports a hierarchically ordered view of the brain, whose connections have heavy tails, but whose weighted node degrees are comparable.