YOLO2U-Net: Detection-Guided 3D Instance Segmentation for Microscopy

Microscopy imaging techniques are instrumental for characterization and analysis of biological structures. As these techniques typically render 3D visualization of cells by stacking 2D projections, issues such as out-of-plane excitation and low resolution in the $z$-axis may pose challenges (even for human experts) to detect individual cells in 3D volumes as these non-overlapping cells may appear as overlapping. In this work, we introduce a comprehensive method for accurate 3D instance segmentation of cells in the brain tissue. The proposed method combines the 2D YOLO detection method with a multi-view fusion algorithm to construct a 3D localization of the cells. Next, the 3D bounding boxes along with the data volume are input to a 3D U-Net network that is designed to segment the primary cell in each 3D bounding box, and in turn, to carry out instance segmentation of cells in the entire volume. The promising performance of the proposed method is shown in comparison with some current deep learning-based 3D instance segmentation methods

Computer Simulations of Cell Sorting Due to Differential Adhesion

The actions of cell adhesion molecules, in particular, cadherins during embryonic development and morphogenesis more generally, regulate many aspects of cellular interactions, regulation and signaling. Often, a gradient of cadherin expression levels drives collective and relative cell motions generating macroscopic cell sorting. Computer simulations of cell sorting have focused on the interactions of cells with only a few discrete adhesion levels between cells, ignoring biologically observed continuous variations in expression levels and possible nonlinearities in molecular binding. In this paper, we present three models relating the surface density of cadherins to the net intercellular adhesion and interfacial tension for both discrete and continuous levels of cadherin expression. We then use then the Glazier-Graner-Hogeweg (GGH) model to investigate how variations in the distribution of the number of cadherins per cell and in the choice of binding model affect cell sorting. We find that an aggregate with a continuous variation in the level of a single type of cadherin molecule sorts more slowly than one with two levels. The rate of sorting increases strongly with the interfacial tension, which depends both on the maximum difference in number of cadherins per cell and on the binding model. Our approach helps connect signaling at the molecular level to tissue-level morphogenesis

Oxygen Tension and the VHL-Hif1α Pathway Determine Onset of Neuronal Polarization and Cerebellar Germinal Zone Exit.

Postnatal brain circuit assembly is driven by temporally regulated intrinsic and cell-extrinsic cues that organize neurogenesis, migration, and axo-dendritic specification in post-mitotic neurons. While cell polarity is an intrinsic organizer of morphogenic events, environmental cues in the germinal zone (GZ) instructing neuron polarization and their coupling during postnatal development are unclear. We report that oxygen tension, which rises at birth, and the von Hippel-Lindau (VHL)-hypoxia-inducible factor 1α (Hif1α) pathway regulate polarization and maturation of post-mitotic cerebellar granule neurons (CGNs). At early postnatal stages with low GZ vascularization, Hif1α restrains CGN-progenitor cell-cycle exit. Unexpectedly, cell-intrinsic VHL-Hif1α pathway activation also delays the timing of CGN differentiation, germinal zone exit, and migration initiation through transcriptional repression of the partitioning-defective (Pard) complex. As vascularization proceeds, these inhibitory mechanisms are downregulated, implicating increasing oxygen tension as a critical switch for neuronal polarization and cerebellar GZ exit

Development of a coupled simulation toolkit for computational radiation biology based on Geant4 and CompuCell3D

RL acknowledges support from Consortium for Risk Evaluation and Stakeholder Participation (http://cresp.org). JAG acknowledges support from National Science Foundation grant NSF 1720625 and National Institutes of Health, National Institute of General Medical Sciences grants U01 GM111243 and R01 GM076692, JAG and MS acknowledge support from National Institutes of Health, National Institute of General Medical Sciences grant R01 GM122424.Understanding and designing clinical radiation therapy is one of the most important areas of state-of-the-art oncological treatment regimens. Decades of research have gone into developing sophisticated treatment devices and optimization protocols for schedules and dosages. In this paper, we presented a comprehensive computational platform that facilitates building of the sophisticated multi-cell-based model of how radiation affects the biology of living tissue. We designed and implemented a coupled simulation method, including a radiation transport model, and a cell biology model, to simulate the tumor response after irradiation. The radiation transport simulation was implemented through Geant4 which is an open-source Monte Carlo simulation platform that provides many flexibilities for users, as well as low energy DNA damage simulation physics, Geant4-DNA. The cell biology simulation was implemented using CompuCell3D (CC3D) which is a cell biology simulation platform. In order to couple Geant4 solver with CC3D, we developed a "bridging" module, RADCELL, that extracts tumor cellular geometry of the CC3D simulation (including specification of the individual cells) and ported it to the Geant4 for radiation transport simulation. The cell dose and cell DNA damage distribution in multicellular system were obtained using Geant4. The tumor response was simulated using cell-based tissue models based on CC3D, and the cell dose and cell DNA damage information were fed back through RADCELL to CC3D for updating the cell properties. By merging two powerful and widely used modeling platforms, CC3D and Geant4, we delivered a novel tool that can give us the ability to simulate the dynamics of biological tissue in the presence of ionizing radiation, which provides a framework for quantifying the biological consequences of radiation therapy. In this introductory methods paper, we described our modeling platform in detail and showed how it can be applied to study the application of radiotherapy to a vascularized tumor.PostprintPeer reviewe

Multiscale modelling of vascular tumour growth in 3D: the roles of domain size & boundary condition

We investigate a three-dimensional multiscale model of vascular tumour growth, which couples blood flow, angiogenesis, vascular remodelling, nutrient/growth factor transport, movement of, and interactions between, normal and tumour cells, and nutrient-dependent cell cycle dynamics within each cell. In particular, we determine how the domain size, aspect ratio and initial vascular network influence the tumour's growth dynamics and its long-time composition. We establish whether it is possible to extrapolate simulation results obtained for small domains to larger ones, by constructing a large simulation domain from a number of identical subdomains, each subsystem initially comprising two parallel parent vessels, with associated cells and diffusible substances. We find that the subsystem is not representative of the full domain and conclude that, for this initial vessel geometry, interactions between adjacent subsystems contribute to the overall growth dynamics. We then show that extrapolation of results from a small subdomain to a larger domain can only be made if the subdomain is sufficiently large and is initialised with a sufficiently complex vascular network. Motivated by these results, we perform simulations to investigate the tumour's response to therapy and show that the probability of tumour elimination in a larger domain can be extrapolated from simulation results on a smaller domain. Finally, we demonstrate how our model may be combined with experimental data, to predict the spatio-temporal evolution of a vascular tumour

Adhesion Failures Determine the Pattern of Choroidal Neovascularization in the Eye: A Computer Simulation Study

Choroidal neovascularization (CNV) of the macular area of the retina is the major cause of severe vision loss in adults. In CNV, after choriocapillaries initially penetrate Bruch's membrane (BrM), invading vessels may regress or expand (CNV initiation). Next, during Early and Late CNV, the expanding vasculature usually spreads in one of three distinct patterns: in a layer between BrM and the retinal pigment epithelium (sub-RPE or Type 1 CNV), in a layer between the RPE and the photoreceptors (sub-retinal or Type 2 CNV) or in both loci simultaneously (combined pattern or Type 3 CNV). While most studies hypothesize that CNV primarily results from growth-factor effects or holes in BrM, our three-dimensional simulations of multi-cell model of the normal and pathological maculae recapitulate the three growth patterns, under the hypothesis that CNV results from combinations of impairment of: 1) RPE-RPE epithelial junctional adhesion, 2) Adhesion of the RPE basement membrane complex to BrM (RPE-BrM adhesion), and 3) Adhesion of the RPE to the photoreceptor outer segments (RPE-POS adhesion). Our key findings are that when an endothelial tip cell penetrates BrM: 1) RPE with normal epithelial junctions, basal attachment to BrM and apical attachment to POS resists CNV. 2) Small holes in BrM do not, by themselves, initiate CNV. 3) RPE with normal epithelial junctions and normal apical RPE-POS adhesion, but weak adhesion to BrM (e.g. due to lipid accumulation in BrM) results in Early sub-RPE CNV. 4) Normal adhesion of RBaM to BrM, but reduced apical RPE-POS or epithelial RPE-RPE adhesion (e.g. due to inflammation) results in Early sub-retinal CNV. 5) Simultaneous reduction in RPE-RPE epithelial binding and RPE-BrM adhesion results in either sub-RPE or sub-retinal CNV which often progresses to combined pattern CNV. These findings suggest that defects in adhesion dominate CNV initiation and progression

Contact-inhibited chemotaxis in de novo and sprouting blood-vessel growth

Blood vessels form either when dispersed endothelial cells (the cells lining the inner walls of fully-formed blood vessels) organize into a vessel network (vasculogenesis), or by sprouting or splitting of existing blood vessels (angiogenesis). Although they are closely related biologically, no current model explains both phenomena with a single biophysical mechanism. Most computational models describe sprouting at the level of the blood vessel, ignoring how cell behavior drives branch splitting during sprouting. We present a cell-based, Glazier-Graner-Hogeweg-model simulation of the initial patterning before the vascular cords form lumens, based on plausible behaviors of endothelial cells. The endothelial cells secrete a chemoattractant, which attracts other endothelial cells. As in the classic Keller-Segel model, chemotaxis by itself causes cells to aggregate into isolated clusters. However, including experimentally-observed adhesion-driven contact inhibition of chemotaxis in the simulation causes randomly-distributed cells to organize into networks and cell aggregates to sprout, reproducing aspects of both de novo and sprouting blood-vessel growth. We discuss two branching instabilities responsible for our results. Cells at the surfaces of cell clusters attempting to migrate to the centers of the clusters produce a buckling instability. In a model variant that eliminates the surface-normal force, a dissipative mechanism drives sprouting, with the secreted chemical acting both as a chemoattractant and as an inhibitor of pseudopod extension. The branching instabilities responsible for our results, which result from contact inhibition of chemotaxis, are both generic developmental mechanisms and interesting examples of unusual patterning instabilities.Comment: Thoroughly revised version, now in press in PLoS Computational Biology. 53 pages, 13 figures, 2 supporting figures, 56 supporting movies, source code and parameters files for computer simulations provided. Supporting information: http://www.psb.ugent.be/~romer/ploscompbiol/ Source code: http://sourceforge.net/projects/tst

Comparing individual-based approaches to modelling the self-organization of multicellular tissues.

The coordinated behaviour of populations of cells plays a central role in tissue growth and renewal. Cells react to their microenvironment by modulating processes such as movement, growth and proliferation, and signalling. Alongside experimental studies, computational models offer a useful means by which to investigate these processes. To this end a variety of cell-based modelling approaches have been developed, ranging from lattice-based cellular automata to lattice-free models that treat cells as point-like particles or extended shapes. However, it remains unclear how these approaches compare when applied to the same biological problem, and what differences in behaviour are due to different model assumptions and abstractions. Here, we exploit the availability of an implementation of five popular cell-based modelling approaches within a consistent computational framework, Chaste (http://www.cs.ox.ac.uk/chaste). This framework allows one to easily change constitutive assumptions within these models. In each case we provide full details of all technical aspects of our model implementations. We compare model implementations using four case studies, chosen to reflect the key cellular processes of proliferation, adhesion, and short- and long-range signalling. These case studies demonstrate the applicability of each model and provide a guide for model usage

3D Multi-Cell Simulation of Tumor Growth and Angiogenesis

We present a 3D multi-cell simulation of a generic simplification of vascular tumor growth which can be easily extended and adapted to describe more specific vascular tumor types and host tissues. Initially, tumor cells proliferate as they take up the oxygen which the pre-existing vasculature supplies. The tumor grows exponentially. When the oxygen level drops below a threshold, the tumor cells become hypoxic and start secreting pro-angiogenic factors. At this stage, the tumor reaches a maximum diameter characteristic of an avascular tumor spheroid. The endothelial cells in the pre-existing vasculature respond to the pro-angiogenic factors both by chemotaxing towards higher concentrations of pro-angiogenic factors and by forming new blood vessels via angiogenesis. The tumor-induced vasculature increases the growth rate of the resulting vascularized solid tumor compared to an avascular tumor, allowing the tumor to grow beyond the spheroid in these linear-growth phases. First, in the linear-spherical phase of growth, the tumor remains spherical while its volume increases. Second, in the linear-cylindrical phase of growth the tumor elongates into a cylinder. Finally, in the linear-sheet phase of growth, tumor growth accelerates as the tumor changes from cylindrical to paddle-shaped. Substantial periods during which the tumor grows slowly or not at all separate the exponential from the linear-spherical and the linear-spherical from the linear-cylindrical growth phases. In contrast to other simulations in which avascular tumors remain spherical, our simulated avascular tumors form cylinders following the blood vessels, leading to a different distribution of hypoxic cells within the tumor. Our simulations cover time periods which are long enough to produce a range of biologically reasonable complex morphologies, allowing us to study how tumor-induced angiogenesis affects the growth rate, size and morphology of simulated tumors

Module-based multiscale simulation of angiogenesis in skeletal muscle

<p>Abstract</p> <p>Background</p> <p>Mathematical modeling of angiogenesis has been gaining momentum as a means to shed new light on the biological complexity underlying blood vessel growth. A variety of computational models have been developed, each focusing on different aspects of the angiogenesis process and occurring at different biological scales, ranging from the molecular to the tissue levels. Integration of models at different scales is a challenging and currently unsolved problem.</p> <p>Results</p> <p>We present an object-oriented module-based computational integration strategy to build a multiscale model of angiogenesis that links currently available models. As an example case, we use this approach to integrate modules representing microvascular blood flow, oxygen transport, vascular endothelial growth factor transport and endothelial cell behavior (sensing, migration and proliferation). Modeling methodologies in these modules include algebraic equations, partial differential equations and agent-based models with complex logical rules. We apply this integrated model to simulate exercise-induced angiogenesis in skeletal muscle. The simulation results compare capillary growth patterns between different exercise conditions for a single bout of exercise. Results demonstrate how the computational infrastructure can effectively integrate multiple modules by coordinating their connectivity and data exchange. Model parameterization offers simulation flexibility and a platform for performing sensitivity analysis.</p> <p>Conclusions</p> <p>This systems biology strategy can be applied to larger scale integration of computational models of angiogenesis in skeletal muscle, or other complex processes in other tissues under physiological and pathological conditions.</p