Parameterizable consensus connectomes from the Human Connectome Project: the Budapest Reference Connectome Server v3.0

Connections of the living human brain, on a macroscopic scale, can be mapped by a diffusion MR imaging based workflow. Since the same anatomic regions can be corresponded between distinct brains, one can compare the presence or the absence of the edges, connecting the very same two anatomic regions, among multiple cortices. Previously, we have constructed the consensus braingraphs on 1015 vertices first in five, then in 96 subjects in the Budapest Reference Connectome Server v1.0 and v2.0, respectively. Here we report the construction of the version 3.0 of the server, generating the common edges of the connectomes of variously parameterizable subsets of the 1015-vertex connectomes of 477 subjects of the Human Connectome Project’s 500-subject release. The consensus connectomes are downloadable in CSV and GraphML formats, and they are also visualized on the server’s page. The consensus connectomes of the server can be considered as the “average, healthy” human connectome since all of their connections are present in at least k subjects, where the default value of \(k=209\), but it can also be modified freely at the web server. The webserver is available at http://connectome.pitgroup.org.     

Information Transfer and Criticality in the Ising Model on the Human Connectome

We implement the Ising model on a structural connectivity matrix describing the brain at two different resolutions. Tuning the model temperature to its critical value, i.e. at the susceptibility peak, we find a maximal amount of total information transfer between the spin variables. At this point the amount of information that can be redistributed by some nodes reaches a limit and the net dynamics exhibits signature of the law of diminishing marginal returns, a fundamental principle connected to saturated levels of production. Our results extend the recent analysis of dynamical oscillators models on the connectome structure, taking into account lagged and directional influences, focusing only on the nodes that are more prone to became bottlenecks of information. The ratio between the outgoing and the incoming information at each node is related to the the sum of the weights to that node and to the average time between consecutive time flips of spins. The results for the connectome of 66 nodes and for that of 998 nodes are similar, thus suggesting that these properties are scale-independent. Finally, we also find that the brain dynamics at criticality is organized maximally to a rich-club w.r.t. the network of information flows.     

Probabilistic Models to Describe the Dynamics of Migrating Microbial Communities

In all but the most sterile environments bacteria will reside in fluid being transported through conduits and some of these will attach and grow as biofilms on the conduit walls. The concentration and diversity of bacteria in the fluid at the point of delivery will be a mix of those when it entered the conduit and those that have become entrained into the flow due to seeding from biofilms. Examples include fluids through conduits such as drinking water pipe networks, endotracheal tubes, catheters and ventilation systems. Here we present two probabilistic models to describe changes in the composition of bulk fluid microbial communities as they are transported through a conduit whilst exposed to biofilm communities. The first (discrete) model simulates absolute numbers of individual cells, whereas the other (continuous) model simulates the relative abundance of taxa in the bulk fluid. The discrete model is founded on a birth-death process whereby the community changes one individual at a time and the numbers of cells in the system can vary. The continuous model is a stochastic differential equation derived from the discrete model and can also accommodate changes in the carrying capacity of the bulk fluid. These models provide a novel Lagrangian framework to investigate and predict the dynamics of migrating microbial communities. In this paper we compare the two models, discuss their merits, possible applications and present simulation results in the context of drinking water distribution systems. Our results provide novel insight into the effects of stochastic dynamics on the composition of non-stationary microbial communities that are exposed to biofilms and provides a new avenue for modelling microbial dynamics in systems where fluids are being transported.     

Semi-Metric Topology of the Human Connectome: Sensitivity and Specificity to Autism and Major Depressive Disorder

Introduction

The human functional connectome is a graphical representation, consisting of nodes connected by edges, of the inter-relationships of blood oxygenation-level dependent (BOLD) time-series measured by MRI from regions encompassing the cerebral cortices and, often, the cerebellum. Semi-metric analysis of the weighted, undirected connectome distinguishes an edge as either direct (metric), such that there is no alternative path that is accumulatively stronger, or indirect (semi-metric), where one or more alternative paths exist that have greater strength than the direct edge. The sensitivity and specificity of this method of analysis is illustrated by two case-control analyses with independent, matched groups of adolescents with autism spectrum conditions (ASC) and major depressive disorder (MDD).

Results

Significance differences in the global percentage of semi-metric edges was observed in both groups, with increases in ASC and decreases in MDD relative to controls. Furthermore, MDD was associated with regional differences in left frontal and temporal lobes, the right limbic system and cerebellum. In contrast, ASC had a broadly increased percentage of semi-metric edges with a more generalised distribution of effects and some areas of reduction. In summary, MDD was characterised by localised, large reductions in the percentage of semi-metric edges, whilst ASC is characterised by more generalised, subtle increases. These differences were corroborated in greater detail by inspection of the semi-metric backbone for each group; that is, the sub-graph of semi-metric edges present in >90% of participants, and by nodal degree differences in the semi-metric connectome.

Conclusion

These encouraging results, in what we believe is the first application of semi-metric analysis to neuroimaging data, raise confidence in the methodology as potentially capable of detection and characterisation of a range of neurodevelopmental and psychiatric disorders.     

Probabilistic Inference of Biochemical Reactions in Microbial Communities from Metagenomic Sequences

Shotgun metagenomics has been applied to the studies of the functionality of various microbial communities. As a critical analysis step in these studies, biological pathways are reconstructed based on the genes predicted from metagenomic shotgun sequences. Pathway reconstruction provides insights into the functionality of a microbial community and can be used for comparing multiple microbial communities. The utilization of pathway reconstruction, however, can be jeopardized because of imperfect functional annotation of genes, and ambiguity in the assignment of predicted enzymes to biochemical reactions (e.g., some enzymes are involved in multiple biochemical reactions). Considering that metabolic functions in a microbial community are carried out by many enzymes in a collaborative manner, we present a probabilistic sampling approach to profiling functional content in a metagenomic dataset, by sampling functions of catalytically promiscuous enzymes within the context of the entire metabolic network defined by the annotated metagenome. We test our approach on metagenomic datasets from environmental and human-associated microbial communities. The results show that our approach provides a more accurate representation of the metabolic activities encoded in a metagenome, and thus improves the comparative analysis of multiple microbial communities. In addition, our approach reports likelihood scores of putative reactions, which can be used to identify important reactions and metabolic pathways that reflect the environmental adaptation of the microbial communities. Source code for sampling metabolic networks is available online at http://omics.informatics.indiana.edu/mg/MetaNetSam/.     

Dynamic Hubs Show Competitive and Static Hubs Non-Competitive Regulation of Their Interaction Partners

Date hub proteins have 1 or 2 interaction interfaces but many interaction partners. This raises the question of whether all partner proteins compete for the interaction interface of the hub or if the cell carefully regulates aspects of this process? Here, we have used real-time rendering of protein interaction networks to analyse the interactions of all the 1 or 2 interface hubs of Saccharomyces cerevisiae during the cell cycle. By integrating previously determined structural and gene expression data, and visually hiding the nodes (proteins) and their edges (interactions) during their troughs of expression, we predict when interactions of hubs and their partners are likely to exist. This revealed that 20 out of all 36 one- or two- interface hubs in the yeast interactome fell within two main groups. The first was dynamic hubs with static partners, which can be considered as ‘competitive hubs’. Their interaction partners will compete for the interaction interface of the hub and the success of any interaction will be dictated by the kinetics of interaction (abundance and affinity) and subcellular localisation. The second was static hubs with dynamic partners, which we term ‘non-competitive hubs’. Regulatory mechanisms are finely tuned to lessen the presence and/or effects of competition between the interaction partners of the hub. It is possible that these regulatory processes may also be used by the cell for the regulation of other, non-cell cycle processes.     

On the Origins of Suboptimality in Human Probabilistic Inference

Humans have been shown to combine noisy sensory information with previous experience (priors), in qualitative and sometimes quantitative agreement with the statistically-optimal predictions of Bayesian integration. However, when the prior distribution becomes more complex than a simple Gaussian, such as skewed or bimodal, training takes much longer and performance appears suboptimal. It is unclear whether such suboptimality arises from an imprecise internal representation of the complex prior, or from additional constraints in performing probabilistic computations on complex distributions, even when accurately represented. Here we probe the sources of suboptimality in probabilistic inference using a novel estimation task in which subjects are exposed to an explicitly provided distribution, thereby removing the need to remember the prior. Subjects had to estimate the location of a target given a noisy cue and a visual representation of the prior probability density over locations, which changed on each trial. Different classes of priors were examined (Gaussian, unimodal, bimodal). Subjects'' performance was in qualitative agreement with the predictions of Bayesian Decision Theory although generally suboptimal. The degree of suboptimality was modulated by statistical features of the priors but was largely independent of the class of the prior and level of noise in the cue, suggesting that suboptimality in dealing with complex statistical features, such as bimodality, may be due to a problem of acquiring the priors rather than computing with them. We performed a factorial model comparison across a large set of Bayesian observer models to identify additional sources of noise and suboptimality. Our analysis rejects several models of stochastic behavior, including probability matching and sample-averaging strategies. Instead we show that subjects'' response variability was mainly driven by a combination of a noisy estimation of the parameters of the priors, and by variability in the decision process, which we represent as a noisy or stochastic posterior.     

Functional Cortical Hubs in the Eyes-Closed Resting Human Brain from an Electrophysiological Perspective Using Magnetoencephalography

It is not clear whether specific brain areas act as hubs in the eyes-closed (EC) resting state, which is an unconstrained state free from any passive or active tasks. Here, we used electrophysiological magnetoencephalography (MEG) signals to study functional cortical hubs in 88 participants. We identified several multispectral cortical hubs. Although cortical hubs vary slightly with different applied measures and frequency bands, the most consistent hubs were observed in the medial and posterior cingulate cortex, the left dorsolateral superior frontal cortex, and the left pole of the middle temporal cortex. Hubs were characterized as connector nodes integrating EC resting state functional networks. Hubs in the gamma band were more likely to include midline structures. Our results confirm the existence of multispectral cortical cores in EC resting state functional networks based on MEG and imply the existence of optimized functional networks in the resting brain.     

Connectotyping: Model Based Fingerprinting of the Functional Connectome

A better characterization of how an individual’s brain is functionally organized will likely bring dramatic advances to many fields of study. Here we show a model-based approach toward characterizing resting state functional connectivity MRI (rs-fcMRI) that is capable of identifying a so-called “connectotype”, or functional fingerprint in individual participants. The approach rests on a simple linear model that proposes the activity of a given brain region can be described by the weighted sum of its functional neighboring regions. The resulting coefficients correspond to a personalized model-based connectivity matrix that is capable of predicting the timeseries of each subject. Importantly, the model itself is subject specific and has the ability to predict an individual at a later date using a limited number of non-sequential frames. While we show that there is a significant amount of shared variance between models across subjects, the model’s ability to discriminate an individual is driven by unique connections in higher order control regions in frontal and parietal cortices. Furthermore, we show that the connectotype is present in non-human primates as well, highlighting the translational potential of the approach.     

An Integrated Mass-Spectrometry Pipeline Identifies Novel Protein Coding-Regions in the Human Genome

Background

Most protein mass spectrometry (MS) experiments rely on searches against a database of known or predicted proteins, limiting their ability as a gene discovery tool.

Results

Using a search against an in silico translation of the entire human genome, combined with a series of annotation filters, we identified 346 putative novel peptides [False Discovery Rate (FDR)<5%] in a MS dataset derived from two human breast epithelial cell lines. A subset of these were then successfully validated by a different MS technique. Two of these correspond to novel isoforms of Heterogeneous Ribonuclear Proteins, while the rest correspond to novel loci.

Conclusions

MS technology can be used for ab initio gene discovery in human data, which, since it is based on different underlying assumptions, identifies protein-coding genes not found by other techniques. As MS technology continues to evolve, such approaches will become increasingly powerful.     

High-throughput Functional Genomics Identifies Regulators of Primary Human Beta Cell Proliferation

The expansion of cells for regenerative therapy will require the genetic dissection of complex regulatory mechanisms governing the proliferation of non-transformed human cells. Here, we report the development of a high-throughput RNAi screening strategy specifically for use in primary cells and demonstrate that silencing the cell cycle-dependent kinase inhibitors CDKN2C/p18 or CDKN1A/p21 facilitates cell cycle entry of quiescent adult human pancreatic beta cells. This work identifies p18 and p21 as novel targets for promoting proliferation of human beta cells and demonstrates the promise of functional genetic screens for dissecting therapeutically relevant state changes in primary human cells.     

Stochastic Blockmodeling of the Modules and Core of the Caenorhabditis elegans Connectome

Recently, there has been much interest in the community structure or mesoscale organization of complex networks. This structure is characterised either as a set of sparsely inter-connected modules or as a highly connected core with a sparsely connected periphery. However, it is often difficult to disambiguate these two types of mesoscale structure or, indeed, to summarise the full network in terms of the relationships between its mesoscale constituents. Here, we estimate a community structure with a stochastic blockmodel approach, the Erdős-Rényi Mixture Model, and compare it to the much more widely used deterministic methods, such as the Louvain and Spectral algorithms. We used the Caenorhabditis elegans (C. elegans) nervous system (connectome) as a model system in which biological knowledge about each node or neuron can be used to validate the functional relevance of the communities obtained. The deterministic algorithms derived communities with 4–5 modules, defined by sparse inter-connectivity between all modules. In contrast, the stochastic Erdős-Rényi Mixture Model estimated a community with 9 blocks or groups which comprised a similar set of modules but also included a clearly defined core, made of 2 small groups. We show that the “core-in-modules” decomposition of the worm brain network, estimated by the Erdős-Rényi Mixture Model, is more compatible with prior biological knowledge about the C. elegans nervous system than the purely modular decomposition defined deterministically. We also show that the blockmodel can be used both to generate stochastic realisations (simulations) of the biological connectome, and to compress network into a small number of super-nodes and their connectivity. We expect that the Erdős-Rényi Mixture Model may be useful for investigating the complex community structures in other (nervous) systems.     

Human and Environmental Impacts on River Sediment Microbial Communities

Sediment microbial communities are responsible for a majority of the metabolic activity in river and stream ecosystems. Understanding the dynamics in community structure and function across freshwater environments will help us to predict how these ecosystems will change in response to human land-use practices. Here we present a spatiotemporal study of sediments in the Tongue River (Montana, USA), comprising six sites along 134 km of river sampled in both spring and fall for two years. Sequencing of 16S rRNA amplicons and shotgun metagenomes revealed that these sediments are the richest (∼65,000 microbial ‘species’ identified) and most novel (93% of OTUs do not match known microbial diversity) ecosystems analyzed by the Earth Microbiome Project to date, and display more functional diversity than was detected in a recent review of global soil metagenomes. Community structure and functional potential have been significantly altered by anthropogenic drivers, including increased pathogenicity and antibiotic metabolism markers near towns and metabolic signatures of coal and coalbed methane extraction byproducts. The core (OTUs shared across all samples) and the overall microbial community exhibited highly similar structure, and phylogeny was weakly coupled with functional potential. Together, these results suggest that microbial community structure is shaped by environmental drivers and niche filtering, though stochastic assembly processes likely play a role as well. These results indicate that sediment microbial communities are highly complex and sensitive to changes in land use practices.     

High-Density Genomewide Linkage Analysis of Exceptional Human Longevity Identifies Multiple Novel Loci

Background

Human lifespan is approximately 25% heritable, and genetic factors may be particularly important for achieving exceptional longevity. Accordingly, siblings of centenarians have a dramatically higher probability of reaching extreme old age than the general population.

Methodology/Principal Findings

To map the loci conferring a survival advantage, we performed the second genomewide linkage scan on human longevity and the first using a high-density marker panel of single nucleotide polymorphisms. By systematically testing a range of minimum age cutoffs in 279 families with multiple long-lived siblings, we identified a locus on chromosome 3p24-22 with a genomewide significant allele-sharing LOD score of 4.02 (empirical P = 0.037) and a locus on chromosome 9q31-34 with a highly suggestive LOD score of 3.89 (empirical P = 0.054). The empirical P value for the combined result was 0.002. A third novel locus with a LOD score of 4.05 on chromosome 12q24 was detected in a subset of the data, and we also obtained modest evidence for a previously reported interval on chromosome 4q22-25.

Conclusions/Significance

Our linkage data should facilitate the discovery of both common and rare variants that determine genetic variability in lifespan.