T-REX: reconstructing and visualizing phylogenetic trees and reticulation networks.

T-REX (tree and reticulogram reconstruction) is an application to reconstruct phylogenetic trees and reticulation networks from distance matrices. The application includes a number of tree fitting methods like NJ, UNJ or ADDTREE which have been very popular in phylogenetic analysis. At the same time, the software comprises several new methods of phylogenetic analysis such as: tree reconstruction using weights, tree inference from incomplete distance matrices or modeling a reticulation network for a collection of objects or species. T-REX also allows the user to visualize obtained tree or network structures using Hierarchical, Radial or Axial types of tree drawing and manipulate them interactively. AVAILABILITY: T-REX is a freeware package available online at: http://www.fas.umontreal.ca/biol/casgrain/en/labo/t-rex     

Phenotypic plasticity and optimal timing of metamorphosis under uncertain time constraints

Life-history theory suggests that optimal timing of metamorphosis should depend on growth conditions and time constraints under which individuals develop. Current models cannot make reliable predictions for species in ephemeral habitats where individuals often face an increasing mortality risk over time because these models assume time-invariant mortality rates (i.e., daily mortality rates remain constant) and fixed seasons. We examined the plasticity of growth, development, and body mass at metamorphosis in tadpoles of the tree-hole breeding frog Phrynobatrachus guineensis in relation to an unpredictable time constraint in the field and in controlled experiments along a fixed density and food gradient. Mean mass and age at metamorphosis of sibships were positively correlated with per capita food level. Based on our results, we developed a simple model of the optimal timing of metamorphosis under time-dependent mortality rates showing that development rates are not only adjusted to growth conditions but also to time-variant mortality rates. The increasing mortality rate represents a time constraint that favors a reduced larval period, but because it is based on probabilities of survival it allows a trade-off between development time and mass. We extend this model to different types of time constraints and show that it can predict the range of documented reaction norms. Differences between species in␣the correlation of age and mass at metamorphosis may have evolved due to differences in their time-variant mortality rates.     

Evolution of biological regulation networks under complex environmental constraints

This work investigates the influence of environmental inducers on the organization of cell regulation networks, using a connectionist approach. Protein interactions are modeled by an asymmetrical recurrent network, the units of which take continuous values. In contrast to classical models, we explicitly introduce a genome to encode the architecture of the system. This feature enables us to introduce an evolution model, in which a genetic algorithm that mimics the effects of evolution on proteins mutual interactions is used. We assume an efficient system to respond to persistent environmental stimuli, independently of their amplitude. Results are presented that show a structuration of the network with the emergence of specialized hierarchical structures. These structures seem to drive the system at the edge of chaos, so that it can present adapted responses to significant environmental changes.     

Simultaneously reconstructing viral cross-species transmission history and identifying the underlying constraints

The factors that determine the origin and fate of cross-species transmission events remain unclear for the majority of human pathogens, despite being central for the development of predictive models and assessing the efficacy of prevention strategies. Here, we describe a flexible Bayesian statistical framework to reconstruct virus transmission between different host species based on viral gene sequences, while simultaneously testing and estimating the contribution of several potential predictors of cross-species transmission. Specifically, we use a generalized linear model extension of phylogenetic diffusion to perform Bayesian model averaging over candidate predictors. By further extending this model with branch partitioning, we allow for distinct host transition processes on external and internal branches, thus discriminating between recent cross-species transmissions, many of which are likely to result in dead-end infections, and host shifts that reflect successful onwards transmission in the new host species. Our approach corroborates genetic distance between hosts as a key determinant of both host shifts and cross-species transmissions of rabies virus in North American bats. Furthermore, our results indicate that geographical range overlap is a modest predictor for cross-species transmission, but not for host shifts. Although our evolutionary framework focused on the multi-host reservoir dynamics of bat rabies virus, it is applicable to other pathogens and to other discrete state transition processes.     

Analyzing protein-protein interaction networks

The advent of the "omics" era in biology research has brought new challenges and requires the development of novel strategies to answer previously intractable questions. Molecular interaction networks provide a framework to visualize cellular processes, but their complexity often makes their interpretation an overwhelming task. The inherently artificial nature of interaction detection methods and the incompleteness of currently available interaction maps call for a careful and well-informed utilization of this valuable data. In this tutorial, we aim to give an overview of the key aspects that any researcher needs to consider when working with molecular interaction data sets and we outline an example for interactome analysis. Using the molecular interaction database IntAct, the software platform Cytoscape, and its plugins BiNGO and clusterMaker, and taking as a starting point a list of proteins identified in a mass spectrometry-based proteomics experiment, we show how to build, visualize, and analyze a protein-protein interaction network.     

Analyzing molecular reaction networks

Pathways are typically the central concept in the analysis of biochemical reaction networks. A pathway can be interpreted as a chain of enzymatical reactions performing a specific biological function. A common way to study metabolic networks are minimal pathways that can operate at steady state called elementary modes. The theory of chemical organizations has recently been used to decompose biochemical networks into algebraically closed and self-maintaining subnetworks termed organizations. The aim of this paper is to elucidate the relation between these two concepts. Whereas elementary modes represent the boundaries of the potential behavior of the network, organizations define metabolite compositions that are likely to be present in biological feasible situations. Hence, steady state organizations consist of combinations of elementary modes. On the other hand, it is possible to assign a unique (and possibly empty) set of organizations to each elementary mode, indicating the metabolites accompanying the active pathway in a feasible steady state.     

A practical algorithm for reconstructing level-1 phylogenetic networks

Recently, much attention has been devoted to the construction of phylogenetic networks which generalize phylogenetic trees in order to accommodate complex evolutionary processes. Here, we present an efficient, practical algorithm for reconstructing level-1 phylogenetic networks--a type of network slightly more general than a phylogenetic tree--from triplets. Our algorithm has been made publicly available as the program LEV1ATHAN. It combines ideas from several known theoretical algorithms for phylogenetic tree and network reconstruction with two novel subroutines. Namely, an exponential-time exact and a greedy algorithm both of which are of independent theoretical interest. Most importantly, LEV1ATHAN runs in polynomial time and always constructs a level-1 network. If the data are consistent with a phylogenetic tree, then the algorithm constructs such a tree. Moreover, if the input triplet set is dense and, in addition, is fully consistent with some level-1 network, it will find such a network. The potential of LEV1ATHAN is explored by means of an extensive simulation study and a biological data set. One of our conclusions is that LEV1ATHAN is able to construct networks consistent with a high percentage of input triplets, even when these input triplets are affected by a low to moderate level of noise.     

Evolutionary constraints permeate large metabolic networks

Background  

Metabolic networks show great evolutionary plasticity, because they can differ substantially even among closely related prokaryotes. Any one metabolic network can also effectively compensate for the blockage of individual reactions by rerouting metabolic flux through other pathways. These observations, together with the continual discovery of new microbial metabolic pathways and enzymes, raise the possibility that metabolic networks are only weakly constrained in changing their complement of enzymatic reactions.     

EXAMINE: a computational approach to reconstructing gene regulatory networks

Reverse-engineering of gene networks using linear models often results in an underdetermined system because of excessive unknown parameters. In addition, the practical utility of linear models has remained unclear. We address these problems by developing an improved method, EXpression Array MINing Engine (EXAMINE), to infer gene regulatory networks from time-series gene expression data sets. EXAMINE takes advantage of sparse graph theory to overcome the excessive-parameter problem with an adaptive-connectivity model and fitting algorithm. EXAMINE also guarantees that the most parsimonious network structure will be found with its incremental adaptive fitting process. Compared to previous linear models, where a fully connected model is used, EXAMINE reduces the number of parameters by O(N), thereby increasing the chance of recovering the underlying regulatory network. The fitting algorithm increments the connectivity during the fitting process until a satisfactory fit is obtained. We performed a systematic study to explore the data mining ability of linear models. A guideline for using linear models is provided: If the system is small (3-20 elements), more than 90% of the regulation pathways can be determined correctly. For a large-scale system, either clustering is needed or it is necessary to integrate information in addition to expression profile. Coupled with the clustering method, we applied EXAMINE to rat central nervous system development (CNS) data with 112 genes. We were able to efficiently generate regulatory networks with statistically significant pathways that have been predicted previously.     

Analyzing and visualizing residue networks of protein structures

The study of individual amino acid residues and their molecular interactions in protein structures is crucial for understanding structure-function relationships. Recent work has indicated that residue networks derived from 3D protein structures provide additional insights into the structural and functional roles of interacting residues. Here, we present the new software tools RINerator and RINalyzer for the automatized generation, 2D visualization, and interactive analysis of residue interaction networks, and highlight their use in different application scenarios.     

Genotype networks shed light on evolutionary constraints

An evolutionary constraint is a bias or limitation in phenotypic variation that a biological system produces. One can distinguish physicochemical, selective, genetic and developmental causes of such constraints. Here, I discuss these causes in three classes of system that bring forth many phenotypic traits and evolutionary innovations: regulatory circuits, macromolecules and metabolic networks. In these systems, genotypes with the same phenotype form large genotype networks that extend throughout a vast genotype space. Such genotype networks can help unify different causes of evolutionary constraints. They can show that these causes ultimately emerge from the process of development; that is, how phenotypes form from genotypes. Furthermore, they can explain important consequences of constraints, such as punctuated stasis and canalization.     

Investigating meta-approaches for reconstructing gene networks in a mammalian cellular context

The output of state-of-the-art reverse-engineering methods for biological networks is often based on the fitting of a mathematical model to the data. Typically, different datasets do not give single consistent network predictions but rather an ensemble of inconsistent networks inferred under the same reverse-engineering method that are only consistent with the specific experimentally measured data. Here, we focus on an alternative approach for combining the information contained within such an ensemble of inconsistent gene networks called meta-analysis, to make more accurate predictions and to estimate the reliability of these predictions. We review two existing meta-analysis approaches; the Fisher transformation combined coefficient test (FTCCT) and Fisher's inverse combined probability test (FICPT); and compare their performance with five well-known methods, ARACNe, Context Likelihood or Relatedness network (CLR), Maximum Relevance Minimum Redundancy (MRNET), Relevance Network (RN) and Bayesian Network (BN). We conducted in-depth numerical ensemble simulations and demonstrated for biological expression data that the meta-analysis approaches consistently outperformed the best gene regulatory network inference (GRNI) methods in the literature. Furthermore, the meta-analysis approaches have a low computational complexity. We conclude that the meta-analysis approaches are a powerful tool for integrating different datasets to give more accurate and reliable predictions for biological networks.     

Physiological regulatory networks: ecological roles and evolutionary constraints

Ecological and evolutionary physiology has traditionally focused on one aspect of physiology at a time. Here, we discuss the implications of considering physiological regulatory networks (PRNs) as integrated wholes, a perspective that reveals novel roles for physiology in organismal ecology and evolution. For example, evolutionary response to changes in resource abundance might be constrained by the role of dietary micronutrients in immune response regulation, given a particular pathogen environment. Because many physiological components impact more than one process, organismal homeostasis is maintained, individual fitness is determined and evolutionary change is constrained (or facilitated) by interactions within PRNs. We discuss how PRN structure and its system-level properties could determine both individual performance and patterns of physiological evolution.     

Deterministic Effects Propagation Networks for reconstructing protein signaling networks from multiple interventions

Background  

Modern gene perturbation techniques, like RNA interference (RNAi), enable us to study effects of targeted interventions in cells efficiently. In combination with mRNA or protein expression data this allows to gain insights into the behavior of complex biological systems.     

Analyzing large biological datasets with association networks

Due to advances in high-throughput biotechnologies biological information is being collected in databases at an amazing rate, requiring novel computational approaches that process collected data into new knowledge in a timely manner. In this study, we propose a computational framework for discovering modular structure, relationships and regularities in complex data. The framework utilizes a semantic-preserving vocabulary to convert records of biological annotations of an object, such as an organism, gene, chemical or sequence, into networks (Anets) of the associated annotations. An association between a pair of annotations in an Anet is determined by the similarity of their co-occurrence pattern with all other annotations in the data. This feature captures associations between annotations that do not necessarily co-occur with each other and facilitates discovery of the most significant relationships in the collected data through clustering and visualization of the Anet. To demonstrate this approach, we applied the framework to the analysis of metadata from the Genomes OnLine Database and produced a biological map of sequenced prokaryotic organisms with three major clusters of metadata that represent pathogens, environmental isolates and plant symbionts.