Abstract:Pathological angiogenesis has been extensively explored by the mathematical modelling community over the past few decades, specifically in the contexts of tumour-induced vascularisation and wound healing. However, there have been relatively few attempts to model angiogenesis associated with normal development, despite the availability of animal models with experimentally accessible and highly ordered vascular topologies: for example, growth and development of the vascular plexus layers in the murine retina. The current study aims to address this issue through the development of a hybrid discrete-continuum mathematical model of the developing retinal vasculature in neonatal mice that is closely coupled with an ongoing experimental programme. The model of the functional vasculature is informed by a range of morphological and molecular data obtained over a period of several days, from 6 days prior to birth to approximately 8 days after birth.The spatio-temporal formation of the superficial retinal vascular plexus (RVP) in wild-type mice occurs in a well-defined sequence. Prior to birth, astrocytes migrate from the optic nerve over the surface of the inner retina in response to a chemotactic gradient of PDGF-A, formed at an earlier stage by migrating retinal ganglion cells (RGCs). Astrocytes express a variety of chemotactic and haptotactic proteins, including VEGF and fibronectin (respectively), which subsequently induce endothelial cell sprouting and modulate growth of the RVP. The developing RVP is not an inert structure; however, the vascular bed adapts and remodels in response to a wide variety of metabolic and biomolecular stimuli. The main focus of this investigation is to understand how these interacting cellular, molecular, and metabolic cues regulate RVP growth and formation.In an earlier one-dimensional continuum model of astrocyte and endothelial migration, we showed that the measured frontal velocities of the two cell types could be accurately reproduced by means of a system of five coupled partial differential equations (Aubert et al. in Bull. Math. Biol. 73:2430–2451, 2011). However, this approach was unable to generate spatial information and structural detail for the entire retinal surface. Building upon this earlier work, a more realistic two-dimensional hybrid PDE-discrete model is derived here that tracks the migration of individual astrocytes and endothelial tip cells towards the outer retinal boundary. Blood perfusion is included throughout plexus development and the emergent retinal architectures adapt and remodel in response to various biological factors. The resulting in silico RVP structures are compared with whole-mounted retinal vasculatures at various stages of development, and the agreement is found to be excellent. Having successfully benchmarked the model against wild-type data, the effect of transgenic over-expression of various genes is predicted, based on the ocular-specific expression of VEGF-A during murine development. These results can be used to help inform future experimental investigations of signalling pathways in ocular conditions characterised by aberrant angiogenesis.

Integrative modeling of sprout formation in angiogenesis: coupling the VEGFA-Notch signaling in a dynamic stalk-tip cell selection

Tip cell overtaking occurs as a side effect of sprouting in computational models of angiogenesis

A hybrid continuum–discrete modelling approach to predict and control angiogenesis: analysis of combinatorial growth factor and matrix effects on vessel-sprouting morphology

Computational Screening of Tip and Stalk Cell Behavior Proposes a Role for Apelin Signaling in Sprout Progression

Angiogenesis Dynamics: A Computational Model of Intravascular Flow Within a Structural Adaptive Vascular Network

Computational modelling of angiogenesis: The importance of cell rearrangements during vascular growth

A Hybrid Discrete-Continuum Mathematical Model of Pattern Prediction in the Developing Retinal Vasculature

A three-dimensional model with two-body interactions for endothelial cells in angiogenesis

The endothelial tip-stalk cell selection and shuffling during angiogenesis

Dynamic Endothelial Stalk Cell–Matrix Interactions Regulate Angiogenic Sprout Diameter

Spatial–temporal order–disorder transition in angiogenic NOTCH signaling controls cell fate specification

Molecular control of endothelial cell behaviour during blood vessel morphogenesis

Dll4 and Notch signalling couples sprouting angiogenesis and artery formation

Parenchymal cues define Vegfa-driven venous angiogenesis by activating a sprouting competent venous endothelial subtype

Falsifying computational models of angiogenesis through quantitative comparison with in vitro models

Mesoscopic and continuum modelling of angiogenesis

Identification of New Markers of Angiogenic Sprouting Using Transcriptomics: New Role for RND3

Role of cell rearrangement and related signaling pathways in the dynamic process of tip cell selection

Matrix deformations around angiogenic sprouts correlate to sprout dynamics and suggest pulling activity

A nanoporous hydrogel-based model to study chemokine gradient-driven angiogenesis under luminal flow

Functional Angiogenesis Requires Microenvironmental Cues Balancing Endothelial Cell Migration and Proliferation.