Search

Living species, ranging from bacteria to animals, exist in environmental conditions that exhibit spatial and temporal heterogeneity which requires them to adapt. Risk-spreading through spontaneous phenotypic variations is a known concept in ecology, which is used to explain how species may survive when faced with the evolutionary risks associated with temporally varying environments. In order to support a deeper understanding of the adaptive role of spontaneous phenotypic variations in fluctuating environments, we consider a system of non-local partial differential equations modelling the evolutionary dynamics of two competing phenotype-structured populations in the presence of periodically oscillating nutrient levels. The two populations undergo heritable, spontaneous phenotypic variations at different rates. The phenotypic state of each individual is represented by a continuous variable, and the phenotypic landscape of the populations evolves in time due to variations in the nutrient level. Exploiting the analytical tractability of our model, we study the long-time behaviour of the solutions to obtain a detailed mathematical depiction of the evolutionary dynamics. The results suggest that when nutrient levels undergo small and slow oscillations, it is evolutionarily more convenient to rarely undergo spontaneous phenotypic variations. Conversely, under relatively large and fast periodic oscillations in the nutrient levels, which bring about alternating cycles of starvation and nutrient abundance, higher rates of spontaneous phenotypic variations confer a competitive advantage. We discuss the implications of our results in the context of cancer metabolism. [ABSTRACT FROM AUTHOR]

We consider mechanically-induced pattern formation within the framework of a growing, planar, elastic rod attached to an elastic foundation. Through a combination of weakly nonlinear analysis and numerical methods, we identify how the shape and type of buckling (super- or subcritical) depend on material parameters, and a complex phase-space of transition from super- to subcritical is uncovered. We then examine the effect of heterogeneity on buckling and post-buckling behaviour, in the context of a heterogeneous substrate adhesion, elastic stiffness, or growth. We show how the same functional form of heterogeneity in different properties is manifest in a vastly differing post-buckled shape. Finally, a fourth form of heterogeneity, an imperfect foundation, is incorporated and shown to have a more dramatic impact on the buckling instability, a difference that can be qualitatively understood via the weakly nonlinear analysis. [ABSTRACT FROM AUTHOR]

Highlights • A continuum model of macrophages releasing an oncolytic virus within a tumour spheroid. • Predictive modelling of this treatment given in combination with radiotherapy. • Investigation into how radiotherapy and oncolytic virotherapy should be scheduled. Abstract Tumour hypoxia has long presented a challenge for cancer therapy: Poor vascularisation in hypoxic regions hinders both the delivery of chemotherapeutic agents and the response to radiotherapy, and hypoxic cancer cells that survive treatment can trigger tumour regrowth after treatment has ended. Tumour-associated macrophages are attractive vehicles for drug delivery because they localise in hypoxic areas of the tumour. In this paper, we derive a mathematical model for the infiltration of an in vitro tumour spheroid by macrophages that have been engineered to release an oncolytic adenovirus under hypoxic conditions. We use this model to predict the efficacy of treatment schedules in which radiotherapy and the engineered macrophages are given in combination. Our work suggests that engineered macrophages should be introduced immediately after radiotherapy for maximum treatment efficacy. Our model provides a framework that may guide future experiments to determine how multiple rounds of radiotherapy and macrophage virotherapy should be coordinated to maximise therapeutic responses. [ABSTRACT FROM AUTHOR]

Angiogenesis is the process by which new blood vessels form from existing vessels. During angiogenesis, tip cells migrate via diffusion and chemotaxis, new tip cells are introduced through branching, loops form via tip-to-tip and tip-to-sprout anastomosis, and a vessel network forms as endothelial cells, known as stalk cells, follow the paths of tip cells (a process known as the snail-trail). Using a mean-field approximation, we systematically derive one-dimensional non-linear continuum models from a lattice-based cellular automaton model of angiogenesis in the corneal assay, explicitly accounting for cell volume. We compare our continuum models and a well-known phenomenological snail-trail model that is linear in the diffusive, chemotactic and branching terms, with averaged cellular automaton simulation results to distinguish macroscale volume exclusion effects and determine whether linear models can capture them. We conclude that, in general, both linear and non-linear models can be used at low cell densities when single or multi-species exclusion effects are negligible at the macroscale. When cell densities increase, our non-linear model should be used to capture non-linear tip cell behavior that occurs when single-species exclusion effects are pronounced, and alternative models should be derived for non-negligible multi-species exclusion effects. [ABSTRACT FROM AUTHOR]

The adult mammalian kidney has a complex, highly-branched collecting duct epithelium that arises as a ureteric bud sidebranch from an epithelial tube known as the nephric duct. Subsequent branching of the ureteric bud to form the collecting duct tree is regulated by subcellular interactions between the epithelium and a population of mesenchymal cells that surround the tips of outgrowing branches. The mesenchymal cells produce glial cell-line derived neurotrophic factor (GDNF), that binds with RET receptors on the surface of the epithelial cells to stimulate several subcellular pathways in the epithelium. Such interactions are known to be a prerequisite for normal branching development, although competing theories exist for their role in morphogenesis. Here we introduce the first agent-based model of ex vivo kidney uretic branching. Through comparison with experimental data, we show that growth factor-regulated growth mechanisms can explain early epithelial cell branching, but only if epithelial cell division depends in a switch-like way on the local growth factor concentration; cell division occurring only if the driving growth factor level exceeds a threshold. We also show how a recently-developed method, “Approximate Approximate Bayesian Computation”, can be used to infer key model parameters, and reveal the dependency between the parameters controlling a growth factor-dependent growth switch. These results are consistent with a requirement for signals controlling proliferation and chemotaxis, both of which are previously identified roles for GDNF. [ABSTRACT FROM AUTHOR]

Current protocols for delivering radiotherapy are based primarily on tumour stage and nodal and metastases status, even though it is well known that tumours and their microenvironments are highly heterogeneous. It is well established that the local oxygen tension plays an important role in radiation-induced cell death, with hypoxic tumour regions responding poorly to irradiation. Therefore, to improve radiation response, it is important to understand more fully the spatiotemporal distribution of oxygen within a growing tumour before and during fractionated radiation. To this end, we have extended a spatially resolved mathematical model of tumour growth, first proposed by Greenspan (Stud Appl Math 51:317-340, 1972), to investigate the effects of oxygen heterogeneity on radiation-induced cell death. In more detail, cell death due to radiation at each location in the tumour, as determined by the well-known linear-quadratic model, is assumed also to depend on the local oxygen concentration. The oxygen concentration is governed by a reaction-diffusion equation that is coupled to an integro-differential equation that determines the size of the assumed spherically symmetric tumour. We combine numerical and analytical techniques to investigate radiation response of tumours with different intratumoral oxygen distribution profiles. Model simulations reveal a rapid transient increase in hypoxia upon regrowth of the tumour spheroid post-irradiation. We investigate the response to different radiation fractionation schedules and identify a tumour-specific relationship between inter-fraction time and dose per fraction to achieve cure. The rich dynamics exhibited by the model suggest that spatial heterogeneity may be important for predicting tumour response to radiotherapy for clinical applications. [ABSTRACT FROM AUTHOR]

Arboviruses can emerge rapidly and cause explosive epidemics of severe disease. Some of the most epidemiologically important arboviruses, including dengue virus (DENV), Zika virus (ZIKV), Chikungunya (CHIKV) and yellow fever virus (YFV), are transmitted by Aedes mosquitoes, most notably Aedes aegypti and Aedes albopictus. After a mosquito blood feeds on an infected host, virus enters the midgut and infects the midgut epithelium. The virus must then overcome a series of barriers before reaching the mosquito saliva and being transmitted to a new host. The virus must escape from the midgut (known as the midgut escape barrier; MEB), which is thought to be mediated by transient changes in the permeability of the midgut-surrounding basal lamina layer (BL) following blood feeding. Here, we present a mathematical model of the within-mosquito population dynamics of DENV (as a model system for mosquito-borne viruses more generally) that includes the interaction of the midgut and BL which can account for the MEB. Our results indicate a dose-dependency of midgut establishment of infection as well as rate of escape from the midgut: collectively, these suggest that the extrinsic incubation period (EIP)—the time taken for DENV virus to be transmissible after infection—is shortened when mosquitoes imbibe more virus. Additionally, our experimental data indicate that multiple blood feeding events, which more closely mimic mosquito-feeding behavior in the wild, can hasten the course of infections, and our model predicts that this effect is sensitive to the amount of virus imbibed. Our model indicates that mutations to the virus which impact its replication rate in the midgut could lead to even shorter EIPs when double-feeding occurs. Mechanistic models of within-vector viral infection dynamics provide a quantitative understanding of infection dynamics and could be used to evaluate novel interventions that target the mosquito stages of the infection. Author summary: Aedes mosquitoes are the main vectors of dengue virus (DENV), Zika virus (ZIKV) and yellow fever virus (YFV), all of which can cause severe disease in humans with dengue alone infecting an estimated 100–400 million people each year. Understanding the processes that affect whether, and at which rate, mosquitoes may transmit such viruses is, hence, paramount. Here, we present a mathematical model of virus dynamics within infected mosquitoes. By combining the model with novel experimental data, we show that the course of infection is sensitive to the initial dose of virus ingested by the mosquito. The data also indicate that mosquitoes which blood feed subsequent to becoming infected may be able to transmit infection earlier, which is reproduced in the model. This is important as many mosquito species feed multiple times during their lifespan, and any reduction in time to dissemination will increase the number of days that a mosquito is infectious and so enhance the risk of transmission. Our study highlights the key and complementary roles played by mathematical models and experimental data for understanding within-mosquito virus dynamics. [ABSTRACT FROM AUTHOR]