Multi-Fiber Tractography Visualizations for Diffusion MRI Data

In recent years, several new diffusion MRI approaches have been proposed to explore microstructural properties of the white matter, such as Q-ball imaging and spherical deconvolution-based techniques to estimate the orientation distribution function. These methods can describe the estimated diffusion profile with a higher accuracy than the more conventional second-rank diffusion tensor imaging technique. Despite many important advances, there are still inconsistent findings between different models that investigate the “crossing fibers” issue. Due to the high information content and the complex nature of the data, it becomes virtually impossible to interpret and compare results in a consistent manner. In this work, we present novel fiber tractography visualization approaches that provide a more complete picture of the microstructural architecture of fiber pathways: multi-fiber hyperstreamlines and streamribbons. By visualizing, for instance, the estimated fiber orientation distribution along the reconstructed tract in a continuous way, information of the local fiber architecture is combined with the global anatomical information derived from tractography. Facilitating the interpretation of diffusion MRI data, this approach can be useful for comparing different diffusion reconstruction techniques and may improve our understanding of the intricate white matter network.     

Explicating the Face Perception Network with White Matter Connectivity

A network of multiple brain regions is recruited in face perception. Our understanding of the functional properties of this network can be facilitated by explicating the structural white matter connections that exist between its functional nodes. We accomplished this using functional MRI (fMRI) in combination with fiber tractography on high angular resolution diffusion weighted imaging data. We identified the three nodes of the core face network: the “occipital face area” (OFA), the “fusiform face area” (mid-fusiform gyrus or mFus), and the superior temporal sulcus (STS). Additionally, a region of the anterior temporal lobe (aIT), implicated as being important for face perception was identified. Our data suggest that we can further divide the OFA into multiple anatomically distinct clusters – a partitioning consistent with several recent neuroimaging results. More generally, structural white matter connectivity within this network revealed: 1) Connectivity between aIT and mFus, and between aIT and occipital regions, consistent with studies implicating this posterior to anterior pathway as critical to normal face processing; 2) Strong connectivity between mFus and each of the occipital face-selective regions, suggesting that these three areas may subserve different functional roles; 3) Almost no connectivity between STS and mFus, or between STS and the other face-selective regions. Overall, our findings suggest a re-evaluation of the “core” face network with respect to what functional areas are or are not included in this network.     

Principal Networks

Graph representations of brain connectivity have attracted a lot of recent interest, but existing methods for dividing such graphs into connected subnetworks have a number of limitations in the context of neuroimaging. This is an important problem because most cognitive functions would be expected to involve some but not all brain regions. In this paper we outline a simple approach for decomposing graphs, which may be based on any measure of interregional association, into coherent “principal networks”. The technique is based on an eigendecomposition of the association matrix, and is closely related to principal components analysis. We demonstrate the technique using cortical thickness and diffusion tractography data, showing that the subnetworks which emerge are stable, meaningful and reproducible. Graph-theoretic measures of network cost and efficiency may be calculated separately for each principal network. Unlike some other approaches, all available connectivity information is taken into account, and vertices may appear in none or several of the subnetworks. Subject-by-subject “scores” for each principal network may also be obtained, under certain circumstances, and related to demographic or cognitive variables of interest.     

Contact Trees: Network Visualization beyond Nodes and Edges

Node-Link diagrams make it possible to take a quick glance at how nodes (or actors) in a network are connected by edges (or ties). A conventional network diagram of a “contact tree” maps out a root and branches that represent the structure of nodes and edges, often without further specifying leaves or fruits that would have grown from small branches. By furnishing such a network structure with leaves and fruits, we reveal details about “contacts” in our ContactTrees upon which ties and relationships are constructed. Our elegant design employs a bottom-up approach that resembles a recent attempt to understand subjective well-being by means of a series of emotions. Such a bottom-up approach to social-network studies decomposes each tie into a series of interactions or contacts, which can help deepen our understanding of the complexity embedded in a network structure. Unlike previous network visualizations, ContactTrees highlight how relationships form and change based upon interactions among actors, as well as how relationships and networks vary by contact attributes. Based on a botanical tree metaphor, the design is easy to construct and the resulting tree-like visualization can display many properties at both tie and contact levels, thus recapturing a key ingredient missing from conventional techniques of network visualization. We demonstrate ContactTrees using data sets consisting of up to three waves of 3-month contact diaries over the 2004-2012 period, and discuss how this design can be applied to other types of datasets.     

“Topological Significance” Analysis of Gene Expression and Proteomic Profiles from Prostate Cancer Cells Reveals Key Mechanisms of Androgen Response

Background

The problem of prostate cancer progression to androgen independence has been extensively studied. Several studies systematically analyzed gene expression profiles in the context of biological networks and pathways, uncovering novel aspects of prostate cancer. Despite significant research efforts, the mechanisms underlying tumor progression are poorly understood. We applied a novel approach to reconstruct system-wide molecular events following stimulation of LNCaP prostate cancer cells with synthetic androgen and to identify potential mechanisms of androgen-independent progression of prostate cancer.

Methodology/Principal Findings

We have performed concurrent measurements of gene expression and protein levels following the treatment using microarrays and iTRAQ proteomics. Sets of up-regulated genes and proteins were analyzed using our novel concept of “topological significance”. This method combines high-throughput molecular data with the global network of protein interactions to identify nodes which occupy significant network positions with respect to differentially expressed genes or proteins. Our analysis identified the network of growth factor regulation of cell cycle as the main response module for androgen treatment in LNCap cells. We show that the majority of signaling nodes in this network occupy significant positions with respect to the observed gene expression and proteomic profiles elicited by androgen stimulus. Our results further indicate that growth factor signaling probably represents a “second phase” response, not directly dependent on the initial androgen stimulus.

Conclusions/Significance

We conclude that in prostate cancer cells the proliferative signals are likely to be transmitted from multiple growth factor receptors by a multitude of signaling pathways converging on several key regulators of cell proliferation such as c-Myc, Cyclin D and CREB1. Moreover, these pathways are not isolated but constitute an interconnected network module containing many alternative routes from inputs to outputs. If the whole network is involved, a precisely formulated combination therapy may be required to fight the tumor growth effectively.     

A Conserved Rod Domain Phosphotyrosine That Is Targeted by the Phosphatase PTP1B Promotes Keratin 8 Protein Insolubility and Filament Organization

Post-translational modifications are important functional determinants for intermediate filament (IF) proteins. Phosphorylation of IF proteins regulates filament organization, solubility, and cell-protective functions. Most known IF protein phosphorylation sites are serines localized in the variable “head” and “tail” domain regions. By contrast, little is known about site-specific tyrosine phosphorylation or its implications on IF protein function. We used available proteomic data from large scale studies to narrow down potential phospho-tyrosine sites on the simple epithelial IF protein keratin 8 (K8). Validation of the predicted sites using a pan-phosphotyrosine and a site-specific antibody, which we generated, revealed that the highly conserved Tyr-267 in the K8 “rod” domain was basally phosphorylated. The charge at this site was critically important, as demonstrated by altered filament organization of site-directed mutants, Y267F and Y267D, the latter exhibiting significantly diminished solubility. Pharmacological inhibition of the protein-tyrosine phosphatase PTP1B increased K8 Tyr-267 phosphorylation, decreased solubility, and increased K8 filament bundling, whereas PTP1B overexpression had the opposite effects. Furthermore, there was significant co-localization between K8 and a “substrate-trapping” mutant of PTP1B (D181A). Because K8 Tyr-267 is conserved in many IFs (QYE motif), we tested the effect of the paralogous Tyr in glial fibrillary acidic protein (GFAP), which is mutated in Alexander disease (Y242D). Similar to K8, Y242D GFAP exhibited highly irregular filament organization and diminished solubility. Our results implicate the rod domain QYE motif tyrosine as an important determinant of IF assembly and solubility properties that can be dynamically modulated by phosphorylation.     

Automated Quantification and Integrative Analysis of 2D and 3D Mitochondrial Shape and Network Properties

Mitochondrial morphology and function are coupled in healthy cells, during pathological conditions and (adaptation to) endogenous and exogenous stress. In this sense mitochondrial shape can range from small globular compartments to complex filamentous networks, even within the same cell. Understanding how mitochondrial morphological changes (i.e. “mitochondrial dynamics”) are linked to cellular (patho) physiology is currently the subject of intense study and requires detailed quantitative information. During the last decade, various computational approaches have been developed for automated 2-dimensional (2D) analysis of mitochondrial morphology and number in microscopy images. Although these strategies are well suited for analysis of adhering cells with a flat morphology they are not applicable for thicker cells, which require a three-dimensional (3D) image acquisition and analysis procedure. Here we developed and validated an automated image analysis algorithm allowing simultaneous 3D quantification of mitochondrial morphology and network properties in human endothelial cells (HUVECs). Cells expressing a mitochondria-targeted green fluorescence protein (mitoGFP) were visualized by 3D confocal microscopy and mitochondrial morphology was quantified using both the established 2D method and the new 3D strategy. We demonstrate that both analyses can be used to characterize and discriminate between various mitochondrial morphologies and network properties. However, the results from 2D and 3D analysis were not equivalent when filamentous mitochondria in normal HUVECs were compared with circular/spherical mitochondria in metabolically stressed HUVECs treated with rotenone (ROT). 2D quantification suggested that metabolic stress induced mitochondrial fragmentation and loss of biomass. In contrast, 3D analysis revealed that the mitochondrial network structure was dissolved without affecting the amount and size of the organelles. Thus, our results demonstrate that 3D imaging and quantification are crucial for proper understanding of mitochondrial shape and topology in non-flat cells. In summary, we here present an integrative method for unbiased 3D quantification of mitochondrial shape and network properties in mammalian cells.     

Analysing Local Sparseness in the Macaque Brain Network

Understanding the network structure of long distance pathways in the brain is a necessary step towards developing an insight into the brain’s function, organization and evolution. Dense global subnetworks of these pathways have often been studied, primarily due to their functional implications. Instead we study sparse local subnetworks of the pathways to establish the role of a brain area in enabling shortest path communication between its non-adjacent topological neighbours. We propose a novel metric to measure the topological communication load on a vertex due to its immediate neighbourhood, and show that in terms of distribution of this local communication load, a network of Macaque long distance pathways is substantially different from other real world networks and random graph models. Macaque network contains the entire range of local subnetworks, from star-like networks to clique-like networks, while other networks tend to contain a relatively small range of subnetworks. Further, sparse local subnetworks in the Macaque network are not only found across topographical super-areas, e.g., lobes, but also within a super-area, arguing that there is conservation of even relatively short-distance pathways. To establish the communication role of a vertex we borrow the concept of brokerage from social science, and present the different types of brokerage roles that brain areas play, highlighting that not only the thalamus, but also cingulate gyrus and insula often act as “relays” for areas in the neocortex. These and other analysis of communication load and roles of the sparse subnetworks of the Macaque brain provide new insights into the organisation of its pathways.     

Molecular Architecture of Spinal Cord Injury Protein Interaction Network

Spinal cord injury (SCI) is associated with complex pathophysiological processes that follow the primary traumatic event and determine the extent of secondary damage and functional recovery. Numerous reports have used global and hypothesis-driven approaches to identify protein changes that contribute to the overall pathology of SCI in an effort to identify potential therapeutic interventions. In this study, we use a semi-automatic annotation approach to detect terms referring to genes or proteins dysregulated in the SCI literature and develop a curated SCI interactome. Network analysis of the SCI interactome revealed the presence of a rich-club organization corresponding to a “powerhouse” of highly interacting hub-proteins. Studying the modular organization of the network have shown that rich-club proteins cluster into modules that are specifically enriched for biological processes that fall under the categories of cell death, inflammation, injury recognition and systems development. Pathway analysis of the interactome and the rich-club revealed high similarity indicating the role of the rich-club proteins as hubs of the most prominent pathways in disease pathophysiology and illustrating the centrality of pro-and anti-survival signal competition in the pathology of SCI. In addition, evaluation of centrality measures of single nodes within the rich-club have revealed that neuronal growth factor (NGF), caspase 3, and H-Ras are the most central nodes and potentially an interesting targets for therapy. Our integrative approach uncovers the molecular architecture of SCI interactome, and provide an essential resource for evaluating significant therapeutic candidates.     

Multidimensional analysis and detection of informative features in human brain white matter

The white matter contains long-range connections between different brain regions and the organization of these connections holds important implications for brain function in health and disease. Tractometry uses diffusion-weighted magnetic resonance imaging (dMRI) to quantify tissue properties along the trajectories of these connections. Statistical inference from tractometry usually either averages these quantities along the length of each fiber bundle or computes regression models separately for each point along every one of the bundles. These approaches are limited in their sensitivity, in the former case, or in their statistical power, in the latter. We developed a method based on the sparse group lasso (SGL) that takes into account tissue properties along all of the bundles and selects informative features by enforcing both global and bundle-level sparsity. We demonstrate the performance of the method in two settings: i) in a classification setting, patients with amyotrophic lateral sclerosis (ALS) are accurately distinguished from matched controls. Furthermore, SGL identifies the corticospinal tract as important for this classification, correctly finding the parts of the white matter known to be affected by the disease. ii) In a regression setting, SGL accurately predicts “brain age.” In this case, the weights are distributed throughout the white matter indicating that many different regions of the white matter change over the lifespan. Thus, SGL leverages the multivariate relationships between diffusion properties in multiple bundles to make accurate phenotypic predictions while simultaneously discovering the most relevant features of the white matter.     

Exact Hybrid Particle/Population Simulation of Rule-Based Models of Biochemical Systems

Detailed modeling and simulation of biochemical systems is complicated by the problem of combinatorial complexity, an explosion in the number of species and reactions due to myriad protein-protein interactions and post-translational modifications. Rule-based modeling overcomes this problem by representing molecules as structured objects and encoding their interactions as pattern-based rules. This greatly simplifies the process of model specification, avoiding the tedious and error prone task of manually enumerating all species and reactions that can potentially exist in a system. From a simulation perspective, rule-based models can be expanded algorithmically into fully-enumerated reaction networks and simulated using a variety of network-based simulation methods, such as ordinary differential equations or Gillespie''s algorithm, provided that the network is not exceedingly large. Alternatively, rule-based models can be simulated directly using particle-based kinetic Monte Carlo methods. This “network-free” approach produces exact stochastic trajectories with a computational cost that is independent of network size. However, memory and run time costs increase with the number of particles, limiting the size of system that can be feasibly simulated. Here, we present a hybrid particle/population simulation method that combines the best attributes of both the network-based and network-free approaches. The method takes as input a rule-based model and a user-specified subset of species to treat as population variables rather than as particles. The model is then transformed by a process of “partial network expansion” into a dynamically equivalent form that can be simulated using a population-adapted network-free simulator. The transformation method has been implemented within the open-source rule-based modeling platform BioNetGen, and resulting hybrid models can be simulated using the particle-based simulator NFsim. Performance tests show that significant memory savings can be achieved using the new approach and a monetary cost analysis provides a practical measure of its utility.     

Multi-Shell Hybrid Diffusion Imaging (HYDI) at 7 Tesla in TgF344-AD Transgenic Alzheimer Rats

Diffusion weighted imaging (DWI) is widely used to study microstructural characteristics of the brain. Diffusion tensor imaging (DTI) and high-angular resolution imaging (HARDI) are frequently used in radiology and neuroscience research but can be limited in describing the signal behavior in composite nerve fiber structures. Here, we developed and assessed the benefit of a comprehensive diffusion encoding scheme, known as hybrid diffusion imaging (HYDI), composed of 300 DWI volumes acquired at 7-Tesla with diffusion weightings at b = 1000, 3000, 4000, 8000 and 12000 s/mm² and applied it in transgenic Alzheimer rats (line TgF344-AD) that model the full clinico-pathological spectrum of the human disease. We studied and visualized the effects of the multiple concentric “shells” when computing three distinct anisotropy maps–fractional anisotropy (FA), generalized fractional anisotropy (GFA) and normalized quantitative anisotropy (NQA). We tested the added value of the multi-shell q-space sampling scheme, when reconstructing neural pathways using mathematical frameworks from DTI and q-ball imaging (QBI). We show a range of properties of HYDI, including lower apparent anisotropy when using high b-value shells in DTI-based reconstructions, and increases in apparent anisotropy in QBI-based reconstructions. Regardless of the reconstruction scheme, HYDI improves FA-, GFA- and NQA-aided tractography. HYDI may be valuable in human connectome projects and clinical research, as well as magnetic resonance research in experimental animals.     

Exploring Bacterial Organelle Interactomes: A Model of the Protein-Protein Interaction Network in the Pdu Microcompartment

Bacterial microcompartments (MCPs) are protein-bound organelles that carry out diverse metabolic pathways in a wide range of bacteria. These supramolecular assemblies consist of a thin outer protein shell, reminiscent of a viral capsid, which encapsulates sequentially acting enzymes. The most complex MCP elucidated so far is the propanediol utilizing (Pdu) microcompartment. It contains the reactions for degrading 1,2-propanediol. While several experimental studies on the Pdu system have provided hints about its organization, a clear picture of how all the individual components interact has not emerged yet. Here we use co-evolution-based methods, involving pairwise comparisons of protein phylogenetic trees, to predict the protein-protein interaction (PPI) network governing the assembly of the Pdu MCP. We propose a model of the Pdu interactome, from which selected PPIs are further inspected via computational docking simulations. We find that shell protein PduA is able to serve as a “universal hub” for targeting an array of enzymes presenting special N-terminal extensions, namely PduC, D, E, L and P. The varied N-terminal peptides are predicted to bind in the same cleft on the presumptive luminal face of the PduA hexamer. We also propose that PduV, a protein of unknown function with remote homology to the Ras-like GTPase superfamily, is likely to localize outside the MCP, interacting with the protruding β-barrel of the hexameric PduU shell protein. Preliminary experiments involving a bacterial two-hybrid assay are presented that corroborate the existence of a PduU-PduV interaction. This first systematic computational study aimed at characterizing the interactome of a bacterial microcompartment provides fresh insight into the organization of the Pdu MCP.     

An Integrative Approach to the Study of Filamentous Oligomeric Assemblies,with Application to RecA

Oligomeric macromolecules in the cell self-organize into a wide variety of geometrical motifs such as helices, rings or linear filaments. The recombinase proteins involved in homologous recombination present many such assembly motifs. Here, we examine in particular the polymorphic characteristics of RecA, the most studied member of the recombinase family, using an integrative approach that relates local modes of monomer/monomer association to the global architecture of their screw-type organization. In our approach, local modes of association are sampled via docking or Monte Carlo simulations. This enables shedding new light on fiber morphologies that may be adopted by the RecA protein. Two distinct RecA helical morphologies, the so-called “extended” and “compressed” forms, are known to play a role in homologous recombination. We investigate the variability within each form in terms of helical parameters and steric accessibility. We also address possible helical discontinuities in RecA filaments due to multiple monomer-monomer association modes. By relating local interface organization to global filament morphology, the strategies developed here to study RecA self-assembly are particularly well suited to other DNA-binding proteins and to filamentous protein assemblies in general.     

Quilt Plots: A Simple Tool for the Visualisation of Large Epidemiological Data

Background

Graphical representation of data is one of the most easily comprehended forms of explanation. The current study describes a simple visualization tool which may allow greater understanding of medical and epidemiological data.

Method

We propose a simple tool for visualization of data, known as a “quilt plot”, that provides an alternative to presenting large volumes of data as frequency tables. Data from the Australian Needle and Syringe Program survey are used to illustrate “quilt plots”.

Conclusion

Visualization of large volumes of data using “quilt plots” enhances interpretation of medical and epidemiological data. Such intuitive presentations are particularly useful for the rapid assessment of problems in the data which cannot be readily identified by manual review. We recommend that, where possible, “quilt plots” be used along with traditional quantitative assessments of the data as an explanatory data analysis tool.     

Development of Large-Scale Functional Brain Networks in Children

The ontogeny of large-scale functional organization of the human brain is not well understood. Here we use network analysis of intrinsic functional connectivity to characterize the organization of brain networks in 23 children (ages 7–9 y) and 22 young-adults (ages 19–22 y). Comparison of network properties, including path-length, clustering-coefficient, hierarchy, and regional connectivity, revealed that although children and young-adults' brains have similar “small-world” organization at the global level, they differ significantly in hierarchical organization and interregional connectivity. We found that subcortical areas were more strongly connected with primary sensory, association, and paralimbic areas in children, whereas young-adults showed stronger cortico-cortical connectivity between paralimbic, limbic, and association areas. Further, combined analysis of functional connectivity with wiring distance measures derived from white-matter fiber tracking revealed that the development of large-scale brain networks is characterized by weakening of short-range functional connectivity and strengthening of long-range functional connectivity. Importantly, our findings show that the dynamic process of over-connectivity followed by pruning, which rewires connectivity at the neuronal level, also operates at the systems level, helping to reconfigure and rebalance subcortical and paralimbic connectivity in the developing brain. Our study demonstrates the usefulness of network analysis of brain connectivity to elucidate key principles underlying functional brain maturation, paving the way for novel studies of disrupted brain connectivity in neurodevelopmental disorders such as autism.     

Distance-based assessment of the localization of functional annotations in 3D genome reconstructions

Background

Recent studies used the contact data or three-dimensional (3D) genome reconstructions from Hi-C (chromosome conformation capture with next-generation sequencing) to assess the co-localization of functional genomic annotations in the nucleus. These analyses dichotomized data point pairs belonging to a functional annotation as “close” or “far” based on some threshold and then tested for enrichment of “close” pairs. We propose an alternative approach that avoids dichotomization of the data and instead directly estimates the significance of distances within the 3D reconstruction.

Results

We applied this approach to 3D genome reconstructions for Plasmodium falciparum, the causative agent of malaria, and Saccharomyces cerevisiae and compared the results to previous approaches. We found significant 3D co-localization of centromeres, telomeres, virulence genes, and several sets of genes with developmentally regulated expression in P. falciparum; and significant 3D co-localization of centromeres and long terminal repeats in S. cerevisiae. Additionally, we tested the experimental observation that telomeres form three to seven clusters in P. falciparum and S. cerevisiae. Applying affinity propagation clustering to telomere coordinates in the 3D reconstructions yielded six telomere clusters for both organisms.

Conclusions

Distance-based assessment replicated key findings, while avoiding dichotomization of the data (which previously yielded threshold-sensitive results).

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-992) contains supplementary material, which is available to authorized users.     

Cerebellar Cortex Granular Layer Interneurons in the Macaque Monkey Are Functionally Driven by Mossy Fiber Pathways through Net Excitation or Inhibition

The granular layer is the input layer of the cerebellar cortex. It receives information through mossy fibers, which contact local granular layer interneurons (GLIs) and granular layer output neurons (granule cells). GLIs provide one of the first signal processing stages in the cerebellar cortex by exciting or inhibiting granule cells. Despite the importance of this early processing stage for later cerebellar computations, the responses of GLIs and the functional connections of mossy fibers with GLIs in awake animals are poorly understood. Here, we recorded GLIs and mossy fibers in the macaque ventral-paraflocculus (VPFL) during oculomotor tasks, providing the first full inventory of GLI responses in the VPFL of awake primates. We found that while mossy fiber responses are characterized by a linear monotonic relationship between firing rate and eye position, GLIs show complex response profiles characterized by “eye position fields” and single or double directional tunings. For the majority of GLIs, prominent features of their responses can be explained by assuming that a single GLI receives inputs from mossy fibers with similar or opposite directional preferences, and that these mossy fiber inputs influence GLI discharge through net excitatory or inhibitory pathways. Importantly, GLIs receiving mossy fiber inputs through these putative excitatory and inhibitory pathways show different firing properties, suggesting that they indeed correspond to two distinct classes of interneurons. We propose a new interpretation of the information flow through the cerebellar cortex granular layer, in which mossy fiber input patterns drive the responses of GLIs not only through excitatory but also through net inhibitory pathways, and that excited and inhibited GLIs can be identified based on their responses and their intrinsic properties.