ATria: a novel centrality algorithm applied to biological networks

Background The notion of centrality is used to identify ?important? nodes in social networks. Importance of nodes is not well-defined, and many different notions exist in the literature. The challenge of defining centrality in meaningful ways when network edges can be positively or negatively weighted has not been adequately addressed in the literature. Existing centrality algorithms also have a second shortcoming, i.e., the list of the most central nodes are often clustered in a specific region of the network and are not well represented across the network. Methods We address both by proposing Ablatio Triadum (ATria), an iterative centrality algorithm that uses the concept of ?payoffs? from economic theory. Results We compare our algorithm with other known centrality algorithms and demonstrate how ATria overcomes several of their shortcomings. We demonstrate the applicability of our algorithm to synthetic networks as well as biological networks including bacterial co-occurrence networks, sometimes referred to as microbial social networks. Conclusions We show evidence that ATria identifies three different kinds of ?important? nodes in microbial social networks with different potential roles in the community

Inferring directional relationships in microbial communities using signed Bayesian networks

Background: Microbe-microbe and host-microbe interactions in a microbiome play a vital role in both health and disease. However, the structure of the microbial community and the colonization patterns are highly complex to infer even under controlled wet laboratory conditions. In this study, we investigate what information, if any, can be provided by a Bayesian Network (BN) about a microbial community. Unlike the previously proposed Co-occurrence Networks (CoNs), BNs are based on conditional dependencies and can help in revealing complex associations. Results: In this paper, we propose a way of combining a BN and a CoN to construct a signed Bayesian Network (sBN). We report a surprising association between directed edges in signed BNs and known colonization orders. Conclusions: BNs are powerful tools for community analysis and extracting influences and colonization patterns, even though the analysis only uses an abundance matrix with no temporal information. We conclude that directed edges in sBNs when combined with negative correlations are consistent with and strongly suggestive of colonization order

Metagenomics, Metatranscriptomics, and Metabolomics Approaches for Microbiome Analysis

Microbiomes are ubiquitous and are found in the ocean, the soil, and in/on other living organisms. Changes in the microbiome can impact the health of the environmental niche in which they reside. In order to learn more about these communities, different approaches based on data from mul-tiple omics have been pursued. Metagenomics produces a taxonomical profile of the sample, metatranscriptomics helps us to obtain a functional profile, and metabolomics completes the picture by determining which byproducts are being released into the environment. Although each approach provides valuable information separately, we show that, when combined, they paint a more comprehensive picture. We conclude with a review of network-based approaches as applied to integrative studies, which we believe holds the key to in-depth understanding of microbiomes

The Multiscale Systems Immunology project: software for cell-based immunological simulation

<p>Abstract</p> <p>Background</p> <p>Computer simulations are of increasing importance in modeling biological phenomena. Their purpose is to predict behavior and guide future experiments. The aim of this project is to model the early immune response to vaccination by an agent based immune response simulation that incorporates realistic biophysics and intracellular dynamics, and which is sufficiently flexible to accurately model the multi-scale nature and complexity of the immune system, while maintaining the high performance critical to scientific computing.</p> <p>Results</p> <p>The Multiscale Systems Immunology (MSI) simulation framework is an object-oriented, modular simulation framework written in C++ and Python. The software implements a modular design that allows for flexible configuration of components and initialization of parameters, thus allowing simulations to be run that model processes occurring over different temporal and spatial scales.</p> <p>Conclusion</p> <p>MSI addresses the need for a flexible and high-performing agent based model of the immune system.</p

Adhesion between cells, diffusion of growth factors, and elasticity of the AER produce the paddle shape of the chick limb

This paper has been withdrawnComment: This paper has been withdraw

A modelling approach towards Epidermal homoeostasis control

In order to grasp the features arising from cellular discreteness and individuality, in large parts of cell tissue modelling agent-based models are favoured. The subclass of off-lattice models allows for a physical motivation of the intercellular interaction rules. We apply an improved version of a previously introduced off-lattice agent-based model to the steady-state flow equilibrium of skin. The dynamics of cells is determined by conservative and drag forces,supplemented with delta-correlated random forces. Cellular adjacency is detected by a weighted Delaunay triangulation. The cell cycle time of keratinocytes is controlled by a diffusible substance provided by the dermis. Its concentration is calculated from a diffusion equation with time-dependent boundary conditions and varying diffusion coefficients. The dynamics of a nutrient is also taken into account by a reaction-diffusion equation. It turns out that the analysed control mechanism suffices to explain several characteristics of epidermal homoeostasis formation. In addition, we examine the question of how {\em in silico} melanoma with decreased basal adhesion manage to persist within the steady-state flow-equilibrium of the skin.Interestingly, even for melanocyte cell cycle times being substantially shorter than for keratinocytes, tiny stochastic effects can lead to completely different outcomes. The results demonstrate that the understanding of initial states of tumour growth can profit significantly from the application of off-lattice agent-based models in computer simulations.Comment: 23 pages, 7 figures, 1 table; version that is to appear in Journal of Theoretical Biolog

Chaste: an open source C++ library for computational physiology and biology

Chaste - Cancer, Heart And Soft Tissue Environment - is an open source C++ library for the computational simulation of mathematical models developed for physiology and biology. Code development has been driven by two initial applications: cardiac electrophysiology and cancer development. A large number of cardiac electrophysiology studies have been enabled and performed, including high performance computational investigations of defibrillation on realistic human cardiac geometries. New models for the initiation and growth of tumours have been developed. In particular, cell-based simulations have provided novel insight into the role of stem cells in the colorectal crypt. Chaste is constantly evolving and is now being applied to a far wider range of problems. The code provides modules for handling common scientific computing components, such as meshes and solvers for ordinary and partial differential equations (ODEs/PDEs). Re-use of these components avoids the need for researchers to "re-invent the wheel" with each new project, accelerating the rate of progress in new applications. Chaste is developed using industrially-derived techniques, in particular test-driven development, to ensure code quality, re-use and reliability. In this article we provide examples that illustrate the types of problems Chaste can be used to solve, which can be run on a desktop computer. We highlight some scientific studies that have used or are using Chaste, and the insights they have provided. The source code, both for specific releases and the development version, is available to download under an open source Berkeley Software Distribution (BSD) licence at http://www.cs.ox.ac.uk/chaste, together with details of a mailing list and links to documentation and tutorials

MDCK Cystogenesis Driven by Cell Stabilization within Computational Analogues

The study of epithelial morphogenesis is fundamental to increasing our understanding of organ function and disease. Great progress has been made through study of culture systems such as Madin-Darby canine kidney (MDCK) cells, but many aspects of even simple morphogenesis remain unclear. For example, are specific cell actions tightly coupled to the characteristics of the cell's environment or are they more often cell state dependent? How does the single lumen, single cell layer cyst consistently emerge from a variety of cell actions? To improve insight, we instantiated in silico analogues that used hypothesized cell behavior mechanisms to mimic MDCK cystogenesis. We tested them through in vitro experimentation and quantitative validation. We observed novel growth patterns, including a cell behavior shift that began around day five of growth. We created agent-oriented analogues that used the cellular Potts model along with an Iterative Refinement protocol. Following several refinements, we achieved a degree of validation for two separate mechanisms. Both survived falsification and achieved prespecified measures of similarity to cell culture properties. In silico components and mechanisms mapped to in vitro counterparts. In silico, the axis of cell division significantly affects lumen number without changing cell number or cyst size. Reducing the amount of in silico luminal cell death had limited effect on cystogenesis. Simulations provide an observable theory for cystogenesis based on hypothesized, cell-level operating principles

Growth Based Morphogenesis of Vertebrate Limb Bud

Many genes and their regulatory relationships are involved in developmental phenomena. However, by chemical information alone, we cannot fully understand changing organ morphologies through tissue growth because deformation and growth of the organ are essentially mechanical processes. Here, we develop a mathematical model to describe the change of organ morphologies through cell proliferation. Our basic idea is that the proper specification of localized volume source (e.g., cell proliferation) is able to guide organ morphogenesis, and that the specification is given by chemical gradients. We call this idea “growth-based morphogenesis.” We find that this morphogenetic mechanism works if the tissue is elastic for small deformation and plastic for large deformation. To illustrate our concept, we study the development of vertebrate limb buds, in which a limb bud protrudes from a flat lateral plate and extends distally in a self-organized manner. We show how the proportion of limb bud shape depends on different parameters and also show the conditions needed for normal morphogenesis, which can explain abnormal morphology of some mutants. We believe that the ideas shown in the present paper are useful for the morphogenesis of other organs