No Arabic abstract
Eigens quasi-species model describes viruses as ensembles of different mutants of a high fitness master genotype. Mutants are assumed to have lower fitness than the master type, yet they coexist with it forming the quasi-species. When the mutation rate is sufficiently high, the master type no longer survives and gets replaced by a wide range of mutant types, thus destroying the quasi-species. It is the so-called error catastrophe. But natural selection acts on phenotypes, not genotypes, and huge amounts of genotypes yield the same phenotype. An important consequence of this is the appearance of beneficial mutations which increase the fitness of mutants. A model has been recently proposed to describe quasi-species in the presence of beneficial mutations. This model lacks the error catastrophe of Eigens model and predicts a steady state in which the viral population grows exponentially. Extinction can only occur if the infectivity of the quasi-species is so low that this exponential is negative. In this work I investigate the transient of this model when infection is started from a small amount of low fitness virions. I prove that, beyond an initial regime where viral population decreases (and can go extinct), the growth of the population is super-exponential. Hence this population quickly becomes so huge that selection due to lack of host cells to be infected begins to act before the steady state is reached. This result suggests that viral infection may widespread before the virus has developed its optimal form.
We study the early stages of viral infection, and the distribution of times to obtain a persistent infection. The virus population proliferates by entering and reproducing inside a target cell until a sufficient number of new virus particles are released via a burst, with a given burst size distribution, which results in the death of the infected cell. Starting with a 2D model describing the joint dynamics of the virus and infected cell populations, we analyze the corresponding master equation using the probability generating function formalism. Exploiting time-scale separation between the virus and infected cell dynamics, the 2D model can be cast into an effective 1D model. To this end, we solve the 1D model analytically for a particular choice of burst size distribution. In the general case, we solve the model numerically by performing extensive Monte-Carlo simulations, and demonstrate the equivalence between the 2D and 1D models by measuring the Kullback-Leibler divergence between the corresponding distributions. Importantly, we find that the distribution of infection times is highly skewed with a fat exponential right tail. This indicates that there is non-negligible portion of individuals with an infection time, significantly longer than the mean, which may have implications on when HIV tests should be performed.
We formulate a compartmental model for the propagation of a respiratory disease in a patchy environment. The patches are connected through the mobility of individuals, and we assume that disease transmission and recovery are possible during travel. Moreover, the migration terms are assumed to depend on the distance between patches and the perceived severity of the disease. The positivity and boundedness of the model solutions are discussed. We analytically show the existence and global asymptotic stability of the disease-free equilibrium. Without human movement, the global stability of the endemic equilibrium point is also established using Lyapunov functions. We study three different network topologies numerically and find that underlying network structure is crucial for disease transmission. Further numerical simulations reveal that infection during travel has the potential to change the stability of disease-free equilibrium from stable to unstable. The coupling strength and transmission coefficients are also very crucial in disease propagation. Different exit screening scenarios indicate that the patch with the highest prevalence may have adverse effects but other patches will be benefited from exit screening. Furthermore, while studying the multi-strain dynamics, it is observed that two co-circulating strains will not persist simultaneously in the community but only one of the strains may persist in the long run. Transmission coefficients corresponding to the second strain are very crucial and show threshold like behavior with respect to the equilibrium density of the second strain.
A question in evolutionary biology is why the number of males is approximately equal to that of females in many species, and Fishers theory of equal investment answers that it is the evolutionarily stable state. The Fisherian mechanism can be given a concrete form by a genetic model based on the following assumptions: (1) Males and females mate at random. (2) An allele acts on the father to determine the expected progeny sex ratio. (3) The offspring inherits the allele from either side of the parents with equal probability. The model is known to achieve the 1:1 sex ratio due to the invasion of mutant alleles with different progeny sex ratios. In this study, however, we argue that mutation plays a more subtle role in that fluctuations caused by mutation renormalize the sex ratio and thereby keep it away from 1:1 in general. This finding shows how the sex ratio is affected by mutation in a systematic way, whereby the effective mutation rate can be estimated from an observed sex ratio.
In this paper we present a discrete dynamical population modeling of invasive species, with reference to the swamp crayfish Procambarus clarkii. Since this species can cause environmental damage of various kinds, it is necessary to evaluate its expected in not yet infested areas. A structured discrete model is built, taking into account all biological information we were able to find, including the environmental variability implemented by means of stochastic parameters (coefficients of fertility, death, etc.). This model is based on a structure with 7 age classes, i.e. a Leslie mathematical population modeling type and it is calibrated with laboratory data provided by the Department of Evolutionary Biology (DEB) of Florence (Italy). The model presents many interesting aspects: the population has a high initial growth, then it stabilizes similarly to the logistic growth, but then it exhibits oscillations (a kind of limit-cycle attractor in the phase plane). The sensitivity analysis shows a good resilience of the model and, for low values of reproductive female fraction, the fluctuations may eventually lead to the extinction of the species: this fact might be exploited as a controlling factor. Moreover, the probability of extinction is valuated with an inverse Gaussian that indicates a high resilience of the species, confirmed by experimental data and field observation: this species has diffused in Italy since 1989 and it has shown a natural tendency to grow. Finally, the spatial mobility is introduced in the model, simulating the movement of the crayfishes in a virtual lake of elliptical form by means of simple cinematic rules encouraging the movement towards the banks of the catchment (as it happens in reality) while a random walk is imposed when the banks are reached.
Among the many aspects that characterize the COVID-19 pandemic, two seem particularly challenging to understand: (i) the great geographical differences in the degree of virus contagiousness and lethality which were found in the different phases of the epidemic progression, and (ii) the potential role of the infected peoples blood type in both the virus infectivity and the progression of the disease. A recent hypothesis could shed some light on both aspects. Specifically, it has been proposed that in the subject-to-subject transfer SARS-CoV-2 conserves on its capsid the erythrocytes antigens of the source subject. Thus these conserved antigens can potentially cause an immune reaction in a receiving subject that has previously acquired specific antibodies for the source subject antigens. This hypothesis implies a blood type-dependent infection rate. The strong geographical dependence of the blood type distribution could be, therefore, one of the factors at the origin of the observed heterogeneity in the epidemics spread. Here, we present an epidemiological deterministic model where the infection rules based on blood types are taken into account and compare our model outcomes with the exiting worldwide infection progression data. We found an overall good agreement, which strengthens the hypothesis that blood types do play a role in the COVID-19 infection.