Cancer is a disease of cellular regulation, often initiated by genetic mutation within cells, and leading to a heterogeneous cell population within tissues. In the competition for nutrients and growth space within the tumors the phenotype of each cell determines its success. Selection in this process is imposed by both the microenvironment (neighboring cells, extracellular matrix, and diffusing substances), and the whole of the organism through for example the blood supply. In this view, the development of tumor cells is in close interaction with their increasingly changing environment: the more cells can change, the more their environment will change. Furthermore, instabilities are also introduced on the organism level: blood supply can be blocked by increased tissue pressure or the tortuousity of the tumor-neovascular vessels. This coupling between cell, microenvironment, and organism results in behavior that is hard to predict. Here we introduce a cell-based computational model to study the effect of blood flow obstruction on the micro-evolution of cells within a tumorous tissue. We demonstrate that stages of tumor development emerge naturally, without the need for sequential mutation of specific genes. Secondly, we show that instabilities in blood supply can impact the overall development of tumors and lead to the extinction of the dominant aggressive phenotype, showing a clear distinction between the fitness at the cell level and fitness at the population level. This provides new insights into potential side effects of recent tumor vasculature renormalization approaches.
Whole genome duplication (WGD) is one of the most important events in the molecular evolution of organisms. In fish species, a WGD is considered to have occurred in the ancestral lineage of teleosts. Recent comprehensive ortholog comparisons among teleost genomes have provided useful data and insights into the fate of redundant genes generated by WGD. Based on these data, a mathematical model is proposed to explain the evolutionary scenario of genes after WGD. The model is parameterized taking into account an equilibrium between i) rapid loss of either of the duplicate genes and ii) moderate functional differentiation of each of duplicate genes, both of which are followed by slow gene loss under purifying selection. This model predicts that, in the teleost lineage, a maximum of about 3000 gene pairs may have differentiated functionally during 90 million years after WGD. Thus, the present study provides a possibility that the whole impact of WGD can be quantitatively assessed according to the model parameters, before details of genomic structural changes or functional differentiation are investigated. If the equilibrium model is valid not only for teleosts but also for other lineages that have undergone WGDs, correlations between the assessment indices and evolutionarily significant events, such as the diversification of species or the occurrence of novel phenotypes, could be tested and compared among those lineages.
Zika virus (ZIKV) exhibits unique transmission dynamics in that it is concurrently spread by a mosquito vector and through sexual contact. We show that this sexual component of ZIKV transmission induces novel processes on networks through the highly asymmetric durations of infectiousness between males and females -- it is estimated that males are infectious for periods up to ten times longer than females -- leading to an asymmetric percolation process on the network of sexual contacts. We exactly solve the properties of this asymmetric percolation on random sexual contact networks and show that this process exhibits two epidemic transitions corresponding to a core-periphery structure. This structure is not present in the underlying contact networks, which are not distinguishable from random networks, and emerges \textitbecause of the asymmetric percolation. We provide an exact analytical description of this double transition and discuss the implications of our results in the context of ZIKV epidemics. Most importantly, our study suggests a bias in our current ZIKV surveillance as the community most at risk is also one of the least likely to get tested.
We analyse the asymptotic behaviour of integro-differential equations modelling N populations in interaction, where interactions are modelled by non-local terms involving linear combinations of the total number of individuals in each population. This model generalises the usual Lotka-Volterra ordinary differential equations. Our aim is to give conditions under which there is global asymptotical stability of coexistence steady-states at the level of the total number of individuals in each species. Through the analysis of a Lyapunov function, our first main result gives a simple and general condition on the matrix of interactions, together with a convergence rate. The second main result establishes another type of condition in the specific case of mutualistic interactions. These conditions are compared to the well-known condition given by Goh for classical Lotka-Volterra ordinary differential equations.