Extracting gene networks for low-dose radiation using graph theoretical algorithms

Genes with common functions often exhibit correlated expression levels, which can be used to identify sets of interacting genes from microarray data. Microarrays typically measure expression across genomic space, creating a massive matrix of co-expression that must be mined to extract only the most relevant gene interactions. We describe a graph theoretical approach to extracting co-expressed sets of genes, based on the computation of cliques. Unlike the results of traditional clustering algorithms, cliques are not disjoint and allow genes to be assigned to multiple sets of interacting partners, consistent with biological reality. A graph is created by thresholding the correlation matrix to include only the correlations most likely to signify functional relationships. Cliques computed from the graph correspond to sets of genes for which significant edges are present between all members of the set, representing potential members of common or interacting pathways. Clique membership can be used to infer function about poorly annotated genes, based on the known functions of better-annotated genes with which they share clique membership (i.e., “guilt-by-association”). We illustrate our method by applying it to microarray data collected from the spleens of mice exposed to low-dose ionizing radiation. Differential analysis is used to identify sets of genes whose interactions are impacted by radiation exposure. The correlation graph is also queried independently of clique to extract edges that are impacted by radiation. We present several examples of multiple gene interactions that are altered by radiation exposure and thus represent potential molecular pathways that mediate the radiation response.     

Identification of regulatory properties of metabolic networks by graph theoretical modeling

An earlier graph theoretical model of metabolic and gene-expression networks has been modified and extended to include the effect of electrical potentials on binding constants, representation of uncatalyzed processes, and treatment of parallel reactions catalyzed by a single enzyme. Formal operations on the graph, which are facilitated by a set of standardized guidelines, identify the feedback signals in the network and rank them according to their influence. The technique was applied to a model of glycolysis in ascites tumor cells in the absence and presence of 12.5 mM exogenous glucose. Feedback regulation was widely distributed and mostly due to binding of adenine nucleotide cofactors to the enzymes of the network. The major changes in feedback regulation on adding glucose is the relief of inhibition of hexokinase and phosphofructokinase and the activation of pyruvate kinase. We conclude that regulation of tumor cell glycolysis is not restricted to hexokinase or to (Na+,K+)-ATPase as was previously suggested by others.     

Graph theoretical analysis of complex networks in the brain

Since the discovery of small-world and scale-free networks the study of complex systems from a network perspective has taken an enormous flight. In recent years many important properties of complex networks have been delineated. In particular, significant progress has been made in understanding the relationship between the structural properties of networks and the nature of dynamics taking place on these networks. For instance, the 'synchronizability' of complex networks of coupled oscillators can be determined by graph spectral analysis. These developments in the theory of complex networks have inspired new applications in the field of neuroscience. Graph analysis has been used in the study of models of neural networks, anatomical connectivity, and functional connectivity based upon fMRI, EEG and MEG. These studies suggest that the human brain can be modelled as a complex network, and may have a small-world structure both at the level of anatomical as well as functional connectivity. This small-world structure is hypothesized to reflect an optimal situation associated with rapid synchronization and information transfer, minimal wiring costs, as well as a balance between local processing and global integration. The topological structure of functional networks is probably restrained by genetic and anatomical factors, but can be modified during tasks. There is also increasing evidence that various types of brain disease such as Alzheimer's disease, schizophrenia, brain tumours and epilepsy may be associated with deviations of the functional network topology from the optimal small-world pattern.     

Disturbed functional brain networks and neurocognitive function in low-grade glioma patients: a graph theoretical analysis of resting-state MEG

Background  

To understand neurophysiological mechanisms underlying cognitive dysfunction in low-grade glioma (LGG) patients by evaluating the spatial structure of 'resting-state' brain networks with graph theory.     

The interaction graph structure of mass-action reaction networks

Behaviour of chemical networks that are described by systems of ordinary differential equations can be analysed via the associated graph structures. This paper deals with observations based on the interaction graph which is defined by the signs of the Jacobian matrix entries. Some of the important graph structures linked to network dynamics are signed circuits and the nucleus (or Hamiltonian hooping). We use mass-action chemical reaction networks as an example to showcase interesting observations about the aforementioned interaction graph structures. We show that positive circuits and specific nucleus structures (associated to multistationarity) are always present in a great generic class of mass-action chemical and biological networks. The theory of negative circuits remains poorly understood, but there is some evidence that they are indicators of stable periodicity. Here we introduce the concept of non-isolated circuits which indicate the presence of a negative circuit.     

The control of ventilation: a theoretical analysis of the response to transient disturbances

Serial measurements of respiration in infants from birth to 7 months have revealed a changing pattern of response with age to spontaneous transient disturbances (e.g. deep sighs). These responses may be interpreted in terms of changing stability of the respiratory control system. The simplest possible model of respiratory control is analysed here. The control system is described by a feedback loop comprising a lung compartment, central and peripheral CO2 detectors and a circulatory delay. The differential equations representing this model are analysed in terms of the phase portrait. First we define the bounds of the solutions in the phase space and the equilibrium points. Next we discuss the local stability and damping of the system under small displacements from these equilibrium points. This local stability can be described in terms of the solution of a linear approximation about the equilibrium points. The trajectories of the solution vector under larger displacements from the equilibrium points are discussed in terms of Taylor expansions. It is shown by this analyses that this simple model can account qualitatively for the observed patterns of respiration following a deep sigh and the changes in that response from birth to 7 months.     

Modular response analysis of cellular regulatory networks

The sheer complexity of intracellular regulatory networks, which involve signal transducing, metabolic, and genetic circuits, hampers our ability to carry out a quantitative analysis of their functions. Here, we describe an approach that greatly simplifies this type of analysis by capitalizing on the modular organization of such networks. Steady-state responses of the network as a whole are accounted for in terms of intermodular interactions between the modules alone; processes operating solely within modules need not be considered when analysing signal transfer through the entire network. The intermodular interactions are quantified through (local) response coefficients which populate an interaction map (matrix). This matrix can be derived from a biochemical or molecular biological analysis of (macro) molecular interactions that constitute the regulatory network. The approach is illustrated by two examples: (i) mitogenic signalling through the mitogen-activated protein kinase cascade in the epidermal growth factor receptor network and (ii) regulation of ammonium assimilation in Escherichia coli.     

From knockouts to networks: establishing direct cause-effect relationships through graph analysis

Background

Reverse-engineering gene networks from expression profiles is a difficult problem for which a multitude of techniques have been developed over the last decade. The yearly organized DREAM challenges allow for a fair evaluation and unbiased comparison of these methods.

Results

We propose an inference algorithm that combines confidence matrices, computed as the standard scores from single-gene knockout data, with the down-ranking of feed-forward edges. Substantial improvements on the predictions can be obtained after the execution of this second step.

Conclusions

Our algorithm was awarded the best overall performance at the DREAM4 In Silico 100-gene network sub-challenge, proving to be effective in inferring medium-size gene regulatory networks. This success demonstrates once again the decisive importance of gene expression data obtained after systematic gene perturbations and highlights the usefulness of graph analysis to increase the reliability of inference.     

Study of biological networks using graph theory

As an effective modeling, analysis and computational tool, graph theory is widely used in biological mathematics to deal with various biology problems. In the field of microbiology, graph can express the molecular structure, where cell, gene or protein can be denoted as a vertex, and the connect element can be regarded as an edge. In this way, the biological activity characteristic can be measured via topological index computing in the corresponding graphs. In our article, we mainly study the biology features of biological networks in terms of eccentric topological indices computation. By means of graph structure analysis and distance calculating, the exact expression of several important eccentric related indices of hypertree network and X-tree are determined. The conclusions we get in this paper illustrate that the bioengineering has the promising application prospects.     

Graph theoretical analysis of the phytosociological structure of plant communities: The theoretical basis

Summary A method employing samples of variable size is assumed and shown to be useful in the detection of unidirectional associations of one species with another, whether caused by direct influence or by coincidence of ecological requirements. The pattern of these associations from a signed digraph (directed graph), the graph theoretical properties of which give insight into the phytosociological structure of the plant community. The properties of transitivity, balance, homogeneity, reachable sets and types of subgraphs are examined.From a thesis submitted to the University of Toronto in partial fulfilment of the requirements for the degree of Master of Science.     

Loss of 'small-world' networks in Alzheimer's disease: graph analysis of FMRI resting-state functional connectivity

Background

Local network connectivity disruptions in Alzheimer''s disease patients have been found using graph analysis in BOLD fMRI. Other studies using MEG and cortical thickness measures, however, show more global long distance connectivity changes, both in functional and structural imaging data. The form and role of functional connectivity changes thus remains ambiguous. The current study shows more conclusive data on connectivity changes in early AD using graph analysis on resting-state condition fMRI data.

Methodology/Principal Findings

18 mild AD patients and 21 healthy age-matched control subjects without memory complaints were investigated in resting-state condition with MRI at 1.5 Tesla. Functional coupling between brain regions was calculated on the basis of pair-wise synchronizations between regional time-series. Local (cluster coefficient) and global (path length) network measures were quantitatively defined. Compared to controls, the characteristic path length of AD functional networks is closer to the theoretical values of random networks, while no significant differences were found in cluster coefficient. The whole-brain average synchronization does not differ between Alzheimer and healthy control groups. Post-hoc analysis of the regional synchronization reveals increased AD synchronization involving the frontal cortices and generalized decreases located at the parietal and occipital regions. This effectively translates in a global reduction of functional long-distance links between frontal and caudal brain regions.

Conclusions/Significance

We present evidence of AD-induced changes in global brain functional connectivity specifically affecting long-distance connectivity. This finding is highly relevant for it supports the anterior-posterior disconnection theory and its role in AD. Our results can be interpreted as reflecting the randomization of the brain functional networks in AD, further suggesting a loss of global information integration in disease.     

Spatio-temporal dependence of the signaling response in immune-receptor trafficking networks regulated by cell density: a theoretical model

Cell signaling processes involve receptor trafficking through highly connected networks of interacting components. The binding of surface receptors to their specific ligands is a key factor for the control and triggering of signaling pathways. In most experimental systems, ligand concentration and cell density vary within a wide range of values. Dependence of the signal response on cell density is related with the extracellular volume available per cell. This dependence has previously been studied using non-spatial models which assume that signaling components are well mixed and uniformly distributed in a single compartment. In this paper, a mathematical model that shows the influence exerted by cell density on the spatio-temporal evolution of ligands, cell surface receptors, and intracellular signaling molecules is developed. To this end, partial differential equations were used to model ligand and receptor trafficking dynamics through the different domains of the whole system. This enabled us to analyze several interesting features involved with these systems, namely: a) how the perturbation caused by the signaling response propagates through the system; b) receptor internalization dynamics and how cell density affects the robustness of dose-response curves upon variation of the binding affinity; and c) that enhanced correlations between ligand input and system response are obtained under conditions that result in larger perturbations of the equilibrium ligand + surface receptor [Please see text] ligand - receptor complex. Finally, the results are compared with those obtained by considering that the above components are well mixed in a single compartment.     

Dynamic generation and qualitative analysis of metabolic pathways by a joint database/graph theoretical approach.

The dynamic generation and qualitative analysis of metabolic networks relying on continuously growing qualified metabolic data by a joint database/graph theoretical approach is described. The procedure is applied to analyze the connectivity of a metabolic network after enzyme removal and to subsequently perform shortest path analyses. The focus lies on the analysis of the connectivity of the metabolic network depending on model assumptions. Here we analyze the influence of the number of strongly connected components on the assignment of reversibility or irreversibility of the biochemical reactions.