Markus R. Owen

Root gravitropism is regulated by a transient lateral auxin gradient controlled by a tipping-point mechanism

Gravity profoundly influences plant growth and development. Plants respond to changes in orientation by using gravitropic responses to modify their growth. Cholodny and Went hypothesized over 80 years ago that plants bend in response to a gravity stimulus by generating a lateral gradient of a growth regulator at an organs apex, later found to be auxin. Auxin regulates root growth by targeting Aux/IAA repressor proteins for degradation. We used an Aux/IAA-based reporter, domain II (DII)-VENUS, in conjunction with a mathematical model to quantify auxin redistribution following a gravity stimulus. Our multidisciplinary approach revealed that auxin is rapidly redistributed to the lower side of the root within minutes of a 90° gravity stimulus. Unexpectedly, auxin asymmetry was rapidly lost as bending root tips reached an angle of 40° to the horizontal. We hypothesize roots use a “tipping point” mechanism that operates to reverse the asymmetric auxin flow at the midpoint of root bending. These mechanistic insights illustrate the scientific value of developing quantitative reporters such as DII-VENUS in conjunction with parameterized mathematical models to provide high-resolution kinetics of hormone redistribution.

Modelling aspects of cancer dynamics : a review

Cancer is a complex disease in which a variety of factors interact over a wide range of spatial and temporal scales with huge datasets relating to the different scales available. However, these data do not always reveal the mechanisms underpinning the observed phenomena. In this paper, we explain why mathematics is a powerful tool for interpreting such data by presenting case studies that illustrate the types of insight that realistic theoretical models of solid tumour growth may yield. These range from discriminating between competing hypotheses for the formation of collagenous capsules associated with benign tumours to predicting the most likely stimulus for protease production in early breast cancer. We will also illustrate the benefits that may result when experimentalists and theoreticians collaborate by considering a novel anti-cancer therapy.

Waves and bumps in neuronal networks with axo-dendritic synaptic interactions

We consider a firing rate model of a neuronal network continuum that incorporates axo-dendritic synaptic processing and the finite conduction velocities of action potentials. The model equation is an integral one defined on a spatially extended domain. Apart from a spatial integral mixing the network connectivity function with space-dependent delays, arising from non-instantaneous axonal communication, the integral model also includes a temporal integration over some appropriately identified distributed delay kernel. These distributed delay kernels are biologically motivated and represent the response of biological synapses to spiking inputs. They are interpreted as Green’s functions of some linear dierential operator. Exploiting this Green’s function description we discuss formal reductions of this non-local system to equivalent partial dierential equation (PDE) models. We distinguish between those spatial connectivity functions that give rise to local PDE models and those that give rise to PDE models with delayed non-local terms. For cases in which local PDEs are derived, we investigate traveling wave solutions in a comoving frame by numerically computing global heteroclinic connections for sigmoidal firing rate functions. We also calculate exact solutions, parameterized by axonal conduction velocity, for the Heaviside firing rate function (the sigmoidal firing rate function in the limit of infinte gain). The inclusion of synaptic adaptation is shown to alter traveling wave fronts to traveling pulses, which we study analytically and numerically in terms of a global homoclinic orbit. Finally, we consider the impact of dendritic interactions on waves and on static spatially localized solutions. Exact analysis for

Angiogenesis and vascular remodelling in normal and cancerous tissues

Vascular development and homeostasis are underpinned by two fundamental features: the generation of new vessels to meet the metabolic demands of under-perfused regions and the elimination of vessels that do not sustain flow. In this paper we develop the first multiscale model of vascular tissue growth that combines blood flow, angiogenesis, vascular remodelling and the subcellular and tissue scale dynamics of multiple cell populations. Simulations show that vessel pruning, due to low wall shear stress, is highly sensitive to the pressure drop across a vascular network, the degree of pruning increasing as the pressure drop increases. In the model, low tissue oxygen levels alter the internal dynamics of normal cells, causing them to release vascular endothelial growth factor (VEGF), which stimulates angiogenic sprouting. Consequently, the level of blood oxygenation regulates the extent of angiogenesis, with higher oxygenation leading to fewer vessels. Simulations show that network remodelling (and de novo network formation) is best achieved via an appropriate balance between pruning and angiogenesis. An important factor is the strength of endothelial tip cell chemotaxis in response to VEGF. When a cluster of tumour cells is introduced into normal tissue, as the tumour grows hypoxic regions form, producing high levels of VEGF that stimulate angiogenesis and cause the vascular density to exceed that for normal tissue. If the original vessel network is sufficiently sparse then the tumour may remain localised near its parent vessel until new vessels bridge the gap to an adjacent vessel. This can lead to metastable periods, during which the tumour burden is approximately constant, followed by periods of rapid growth.

Evans Functions for Integral Neural Field Equations with Heaviside Firing Rate Function

In this paper we show how to construct the Evans function for traveling wave solutions of integral neural field equations when the firing rate function is a Heaviside. This allows a discussion of wave stability and bifurcation as a function of system parameters, including the speed and strength of synaptic coupling and the speed of axonal signals. The theory is illustrated with the construction and stability analysis of front solutions to a scalar neural field model, and a limiting case is shown to recover recent results of Zhang [Differential Integral Equations, 16 (2003), pp. 513--536]. Traveling fronts and pulses are considered in more general models possessing either a linear or piecewise constant recovery variable. We establish the stability of coexisting traveling fronts beyond a front bifurcation and consider parameter regimes that support two stable traveling fronts of different speeds. Such fronts may be connected, and depending on their relative speed the resulting region of activity can widen o...

Multiscale Modelling of Vascular Tumour Growth in 3D: The Roles of Domain Size and Boundary Conditions

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 tumours 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 tumours 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.

Mathematical modeling elucidates the role of transcriptional feedback in gibberellin signaling

The hormone gibberellin (GA) is a key regulator of plant growth. Many of the components of the gibberellin signal transduction [e.g., GIBBERELLIN INSENSITIVE DWARF 1 (GID1) and DELLA], biosynthesis [e.g., GA 20-oxidase (GA20ox) and GA3ox], and deactivation pathways have been identified. Gibberellin binds its receptor, GID1, to form a complex that mediates the degradation of DELLA proteins. In this way, gibberellin relieves DELLA-dependent growth repression. However, gibberellin regulates expression of GID1, GA20ox, and GA3ox, and there is also evidence that it regulates DELLA expression. In this paper, we use integrated mathematical modeling and experiments to understand how these feedback loops interact to control gibberellin signaling. Model simulations are in good agreement with in vitro data on the signal transduction and biosynthesis pathways and in vivo data on the expression levels of gibberellin-responsive genes. We find that GA–GID1 interactions are characterized by two timescales (because of a lid on GID1 that can open and close slowly relative to GA–GID1 binding and dissociation). Furthermore, the model accurately predicts the response to exogenous gibberellin after a number of chemical and genetic perturbations. Finally, we investigate the role of the various feedback loops in gibberellin signaling. We find that regulation of GA20ox transcription plays a significant role in both modulating the level of endogenous gibberellin and generating overshoots after the removal of exogenous gibberellin. Moreover, although the contribution of other individual feedback loops seems relatively small, GID1 and DELLA transcriptional regulation acts synergistically with GA20ox feedback.

Bumps and rings in a two-dimensional neural field: splitting and rotational instabilities

In this paper we consider instabilities of localised solutions in planar neural field firing rate models of Wilson-Cowan or Amari type. Importantly we show that angular perturbations can destabilise spatially localised solutions. For a scalar model with Heaviside firing rate function we calculate symmetric one-bump and ring solutions explicitly and use an Evans function approach to predict the point of instability and the shapes of the dominant growing modes. Our predictions are shown to be in excellent agreement with direct numerical simulations. Moreover, beyond the instability our simulations demonstrate the emergence of multi-bump and labyrinthine patterns. With the addition of spike-frequency adaptation, numerical simulations of the resulting vector model show that it is possible for structures without rotational symmetry, and in particular multi-bumps, to undergo an instability to a rotating wave. We use a general argument, valid for smooth firing rate functions, to establish the conditions necessary to generate such a rotational instability. Numerical continuation of the rotating wave is used to quantify the emergent angular velocity as a bifurcation parameter is varied. Wave stability is found via the numerical evaluation of an associated eigenvalue problem.

Growth-induced hormone dilution can explain the dynamics of plant root cell elongation

In the elongation zone of the Arabidopsis thaliana plant root, cells undergo rapid elongation, increasing their length by ∼10-fold over 5 h while maintaining a constant radius. Although progress is being made in understanding how this growth is regulated, little consideration has been given as to how cell elongation affects the distribution of the key regulating hormones. Using a multiscale mathematical model and measurements of growth dynamics, we investigate the distribution of the hormone gibberellin in the root elongation zone. The model quantifies how rapid cell expansion causes gibberellin to dilute, creating a significant gradient in gibberellin levels. By incorporating the gibberellin signaling network, we simulate how gibberellin dilution affects the downstream components, including the growth-repressing DELLA proteins. We predict a gradient in DELLA that provides an explanation of the reduction in growth exhibited as cells move toward the end of the elongation zone. These results are validated at the molecular level by comparing predicted mRNA levels with transcriptomic data. To explore the dynamics further, we simulate perturbed systems in which gibberellin levels are reduced, considering both genetically modified and chemically treated roots. By modeling these cases, we predict how these perturbations affect gibberellin and DELLA levels and thereby provide insight into their altered growth dynamics.

Sequential induction of auxin efflux and influx carriers regulates lateral root emergence

In Arabidopsis, lateral roots originate from pericycle cells deep within the primary root. New lateral root primordia (LRP) have to emerge through several overlaying tissues. Here, we report that auxin produced in new LRP is transported towards the outer tissues where it triggers cell separation by inducing both the auxin influx carrier LAX3 and cell‐wall enzymes. LAX3 is expressed in just two cell files overlaying new LRP. To understand how this striking pattern of LAX3 expression is regulated, we developed a mathematical model that captures the network regulating its expression and auxin transport within realistic three‐dimensional cell and tissue geometries. Our model revealed that, for the LAX3 spatial expression to be robust to natural variations in root tissue geometry, an efflux carrier is required—later identified to be PIN3. To prevent LAX3 from being transiently expressed in multiple cell files, PIN3 and LAX3 must be induced consecutively, which we later demonstrated to be the case. Our study exemplifies how mathematical models can be used to direct experiments to elucidate complex developmental processes.