Improving abdominal image segmentation with overcomplete shape priors.

The extraction of abdominal structures using deep learning has recently experienced a widespread interest in medical image analysis. Automatic abdominal organ and vessel segmentation is highly desirable to guide clinicians in computer-assisted diagnosis, therapy, or surgical planning. Despite a good ability to extract large organs, the capacity of U-Net inspired architectures to automatically delineate smaller structures remains a major issue, especially given the increase in receptive field size as we go deeper into the network. To deal with various abdominal structure sizes while exploiting efficient geometric constraints, we present a novel approach that integrates into deep segmentation shape priors from a semi-overcomplete convolutional auto-encoder (S-OCAE) embedding. Compared to standard convolutional auto-encoders (CAE), it exploits an over-complete branch that projects data onto higher dimensions to better characterize anatomical structures with a small spatial extent. Experiments on abdominal organs and vessel delineation performed on various publicly available datasets highlight the effectiveness of our method compared to state-of-the-art, including U-Net trained without and with shape priors from a traditional CAE. Exploiting a semi-overcomplete convolutional auto-encoder embedding as shape priors improves the ability of deep segmentation models to provide realistic and accurate abdominal structure contours.

A GPU Algorithm for Agent-Based Models to Simulate the Integration of Cell Membrane Signals.

Simulation of complex biological systems with agent-based models is becoming more relevant with the increase in Graphics Processing Unit (GPU) power. In those simulations, up to millions of virtual cells are individually computed, involving daunting processing times. An important part of computational models is the algorithm that manages how agents perceive their surroundings. This can be particularly problematic in three-dimensional environments where agents have deformable virtual membranes. This article presents a GPU algorithm that gives the possibility for agents to integrate the signals scattered on their virtual membrane. It is detailed to be coded in languages like OpenCL or Cuda. Its performances are tested to show its speed with current GPU devices. Finally, it was implemented inside an existing software to test and illustrate the possibilities it offers.

Probabilistic assessment method of the non-monotonic dose-responses-Part II: Robustness assessment.

In a previous study, we presented a new method that uses a large-scale sampling system to probabilistically assess non-monotonic dose-response curves. The statistical plausibility of the characterization was governed by the probability of the dominant category, but gave no information about the specific robustness of the curve. In this paper we propose an improvement to the method by integrating a scoring system based on 4 criteria which can be used to assess the slope robustness of each of the 10,000 sampled curves. The distribution criterion which assesses the number of doses forming a slope, the intensity criterion which assesses the amplitude of the response, and the minimum and maximum confirmation criteria which increase the certainty that the response is present. The probabilistic method was tested on 294 dose-response curves taken from 2 databases and 2 other methodologies currently proposed. A total of 544 dose-response curves have been processed. The developed method offers a concrete and probabilistic characterization of the type of curve analyzed. It evaluates its statistical plausibility and its robustness according to its sampling curves. This method is applicable to all types of data (continuous and discrete) and all experimental curves starting from theoretically 3 doses at least.

Diversity and survival of artificial lifeforms under sedimentation and random motion.

Cellular automata are often used to explore the numerous possible scenarios of what could have occurred at the origins of life and before, during the prebiotic ages, when very simple molecules started to assemble and organise into larger catalytic or informative structures, or to simulate ecosystems. Artificial self-maintained spatial structures emerge in cellular automata and are often used to represent molecules or living organisms. They converge generally towards homogeneous stationary soups of still-life creatures. It is hard for an observer to believe they are similar to living systems, in particular because nothing is moving anymore within such simulated environments after few computation steps, because they present isotropic spatial organisation, because the diversity of self-maintained morphologies is poor, and because when stationary states are reached the creatures are immortal. Natural living systems, on the contrary, are composed of a high diversity of creatures in interaction having limited lifetimes and generally present a certain anisotropy of their spatial organisation, in particular frontiers and interfaces. In the present work, we propose that the presence of directional weak fields such as gravity may counter-balance the excess of mixing and disorder caused by Brownian motion and favour the appearance of specific regions, i.e. different strata or environmental layers, in which physical-chemical conditions favour the emergence and the survival of self-maintained spatial structures including living systems. We test this hypothesis by way of numerical simulations of a very simplified ecosystem model. We use the well-known Game of Life to which we add rules simulating both sedimentation forces and thermal agitation. We show that this leads to more active (vitality and biodiversity) and robust (survival) dynamics. This effectively suggests that coupling such physical processes to reactive systems allows the separation of environments into different milieux and could constitute a simple mechanism to form ecosystem frontiers or elementary interfaces that would protect and favour the development of fragile auto-poietic systems.

Probabilistic assessment method of the non-monotonic dose-responses-Part I: Methodological approach.

More and more studies aim to characterize non-monotonic dose response curves (NMDRCs). The greatest difficulty is to assess the statistical plausibility of NMDRCs from previously conducted dose response studies. This difficulty is linked to the fact that these studies present (i) few doses tested, (ii) a low sample size per dose, and (iii) the absence of any raw data. In this study, we propose a new methodological approach to probabilistically characterize NMDRCs. The methodology is composed of three main steps: (i) sampling from summary data to cover all the possibilities that may be presented by the responses measured by dose and to obtain a new raw database, (ii) statistical analysis of each sampled dose-response curve to characterize the slopes and their signs, and (iii) characterization of these dose-response curves according to the variation of the sign in the slope. This method allows characterizing all types of dose-response curves and can be applied both to continuous data and to discrete data. The aim of this study is to present the general principle of this probabilistic method which allows to assess the non-monotonic dose responses curves, and to present some results.

Large Scale Tissue Morphogenesis Simulation on Heterogenous Systems Based on a Flexible Biomechanical Cell Model.

The complexity of biological tissue morphogenesis makes in silico simulations of such system very interesting in order to gain a better understanding of the underlying mechanisms ruling the development of multicellular tissues. This complexity is mainly due to two elements: firstly, biological tissues comprise a large amount of cells; secondly, these cells exhibit complex interactions and behaviors. To address these two issues, we propose two tools: the first one is a virtual cell model that comprise two main elements: firstly, a mechanical structure (membrane, cytoskeleton, and cortex) and secondly, the main behaviors exhibited by biological cells, i.e., mitosis, growth, differentiation, molecule consumption, and production as well as the consideration of the physical constraints issued from the environment. An artificial chemistry is also included in the model. This virtual cell model is coupled to an agent-based formalism. The second tool is a simulator that relies on the OpenCL framework. It allows efficient parallel simulations on heterogenous devices such as micro-processors or graphics processors. We present two case studies validating the implementation of our model in our simulator: cellular proliferation controlled by cell signalling and limb growth in a virtual organism.

A software architecture for multi-cellular system simulations on graphics processing units.

The first aim of simulation in virtual environment is to help biologists to have a better understanding of the simulated system. The cost of such simulation is significantly reduced compared to that of in vivo simulation. However, the inherent complexity of biological system makes it hard to simulate these systems on non-parallel architectures: models might be made of sub-models and take several scales into account; the number of simulated entities may be quite large. Today, graphics cards are used for general purpose computing which has been made easier thanks to frameworks like CUDA or OpenCL. Parallelization of models may however not be easy: parallel computer programing skills are often required; several hardware architectures may be used to execute models. In this paper, we present the software architecture we built in order to implement various models able to simulate multi-cellular system. This architecture is modular and it implements data structures adapted for graphics processing units architectures. It allows efficient simulation of biological mechanisms.

Computing with bacterial constituents, cells and populations: from bioputing to bactoputing.

The relevance of biological materials and processes to computing-alias bioputing-has been explored for decades. These materials include DNA, RNA and proteins, while the processes include transcription, translation, signal transduction and regulation. Recently, the use of bacteria themselves as living computers has been explored but this use generally falls within the classical paradigm of computing. Computer scientists, however, have a variety of problems to which they seek solutions, while microbiologists are having new insights into the problems bacteria are solving and how they are solving them. Here, we envisage that bacteria might be used for new sorts of computing. These could be based on the capacity of bacteria to grow, move and adapt to a myriad different fickle environments both as individuals and as populations of bacteria plus bacteriophage. New principles might be based on the way that bacteria explore phenotype space via hyperstructure dynamics and the fundamental nature of the cell cycle. This computing might even extend to developing a high level language appropriate to using populations of bacteria and bacteriophage. Here, we offer a speculative tour of what we term bactoputing, namely the use of the natural behaviour of bacteria for calculating.

Computational energetic model of morphogenesis based on multi-agent Cellular Potts Model.

The Cellular Potts Model (CPM) is a cellular automaton (CA), developed by Glazier and Graner in 1992, to model the morphogenesis. In this model, the entities are the cells. It has already been improved in many ways; however, a key point in biological systems, not defined in CPM, is energetic exchange between entities. We integrate this energetic concept inside the CPM. We simulate a cell differentiation inside a growing cell tissue. The results are the emergence of dynamic patterns coming from the consumption and production of energy. A model described by CA is less scalable than one described by a multi-agent system (MAS). We have developed a MAS based on the CPM, where a cell agent is implemented from the cell of CPM together with several behaviours, in particular the consumption and production of energy from the consumption of molecules.

Model driven quantification of individual and collective cell migration.

While the control of cell migration by biochemical and biophysical factors is largely documented, a precise quantification of cell migration parameters in different experimental contexts is still questionable. Indeed, these phenomenological parameters can be evaluated from data obtained either at the cell population level or at the individual cell level. However, the range within which both characterizations of cell migration are equivalent remains unclear. We analyse here to which extent both sources of data could be integrated within a unified description of cell migration by considering the motility of the endothelial cell line EAhy926. Using time-lapse video-microscopy and associated analysis of digital image time series, we quantified EAhy926 random motility coefficient, migration speed and trajectory persistence time in two different migration assays: the in vitro wound healing assay, and the cell-populated agarose drop assay. In order to analyse the agreement between independent quantifications of cell motility based either on individual cell analysis or cell population dynamic analysis, a theoretical multi-agents cellular model was developed and discussed as a possible theoretical framework able to unify these multi-scale data. Model simulations especially reveal the potential bias induced by cell proliferation and cell-cell adhesion when cell migration parameters are estimated from the extensively used in vitro wound healing assay.

Hyperstructures, genome analysis and I-cells.

New concepts may prove necessary to profit from the avalanche of sequence data on the genome, transcriptome, proteome and interactome and to relate this information to cell physiology. Here, we focus on the concept of large activity-based structures, or hyperstructures, in which a variety of types of molecules are brought together to perform a function. We review the evidence for the existence of hyperstructures responsible for the initiation of DNA replication, the sequestration of newly replicated origins of replication, cell division and for metabolism. The processes responsible for hyperstructure formation include changes in enzyme affinities due to metabolite-induction, lipid-protein affinities, elevated local concentrations of proteins and their binding sites on DNA and RNA, and transertion. Experimental techniques exist that can be used to study hyperstructures and we review some of the ones less familiar to biologists. Finally, we speculate on how a variety of in silico approaches involving cellular automata and multi-agent systems could be combined to develop new concepts in the form of an Integrated cell (I-cell) which would undergo selection for growth and survival in a world of artificial microbiology.