No Arabic abstract
Dengue viral infections show unique infection patterns arising from its four serot- ypes, (DENV-1,2,3,4). Its effects range from simple fever in primary infections to potentially fatal secondary infections. We analytically and numerically analyse virus dynamics and humoral response in a host during primary and secondary dengue infection for long periods using micro-epidemic models. The models presented here incorporate time delays, antibody dependent enhancement (ADE), a dynamic switch and a correlation factor between different DENV serotypes. We find that the viral load goes down to undetectable levels within 7-14 days as is observed for dengue infection, in both cases. For primary infection, the stability analysis of steady states shows interesting dependence on the time delay involved in the production of antibodies from plasma cells. We demonstrate the existence of a critical value for the immune response parameter, beyond which the infection gets completely cured. For secondary infections with a different serotype, the homologous antibody production is enhanced due to the influence of heterologous antibodies. The antibody production is also controlled by the correlation factor, which is a measure of similarities between the different DENV serotypes involved. Our results agree with clinically observed humoral responses for primary and secondary infections.
We discuss several models of the dynamics of interacting populations. The models are constructed by nonlinear differential equations and have two sets of parameters: growth rates and coefficients of interaction between populations. We assume that the parameters depend on the densities of the populations. In addition the parameters can be influenced by different factors of the environment. This influence is modelled by noise terms in the equations for the growth rates and interaction coefficients. Thus the model differential equations become stochastic. In some particular cases these equations can be reduced to a Foker-Plancnk equation for the probability density function of the densities of the interacting populations.
In this work we have investigated the evolutionary dynamics of a generalist pathogen, e.g. a virus population, that evolves towards specialisation in an environment with multiple host types. We have particularly explored under which conditions generalist viral strains may rise in frequency and coexist with specialist strains or even dominate the population. By means of a nonlinear mathematical model and bifurcation analysis, we have determined the theoretical conditions for stability of nine identified equilibria and provided biological interpretation in terms of the infection rates for the viral specialist and generalist strains. By means of a stability diagram we identified stable fixed points and stable periodic orbits, as well as regions of bistability. For arbitrary biologically feasible initial population sizes, the probability of evolving towards stable solutions is obtained for each point of the analyzed parameter space. This probability map shows combinations of infection rates of the generalist and specialist strains that might lead to equal chances for each type becoming the dominant strategy. Furthermore, we have identified infection rates for which the model predicts the onset of chaotic dynamics. Several degenerate Bogdanov-Takens and zero-Hopf bifurcations are detected along with generalized Hopf and zero-Hopf bifurcations. This manuscript provides additional insights into the dynamical complexity of host-pathogen evolution towards different infection strategies.
Near real-time monitoring of outbreak transmission dynamics and evaluation of public health interventions are critical for interrupting the spread of the novel coronavirus (SARS-CoV-2) and mitigating morbidity and mortality caused by coronavirus disease (COVID-19). Formulating a regional mechanistic model of SARS-CoV-2 transmission dynamics and frequently estimating parameters of this model using streaming surveillance data offers one way to accomplish data-driven decision making. For example, to detect an increase in new SARS-CoV-2 infections due to relaxation of previously implemented mitigation measures one can monitor estimates of the basic and effective reproductive numbers. However, parameter estimation can be imprecise, and sometimes even impossible, because surveillance data are noisy and not informative about all aspects of the mechanistic model, even for reasonably parsimonious epidemic models. To overcome this obstacle, at least partially, we propose a Bayesian modeling framework that integrates multiple surveillance data streams. Our model uses both COVID-19 incidence and mortality time series to estimate our model parameters. Importantly, our data generating model for incidence data takes into account changes in the total number of tests performed. We apply our Bayesian data integration method to COVID-19 surveillance data collected in Orange County, California. Our results suggest that California Department of Public Health stay-at-home order, issued on March 19, 2020, lowered the SARS-CoV-2 effective reproductive number $R_{e}$ in Orange County below 1.0, which means that the order was successful in suppressing SARS-CoV-2 infections. However, subsequent re-opening steps took place when thousands of infectious individuals remained in Orange County, so $R_{e}$ increased to approximately 1.0 by mid-June and above 1.0 by mid-July.
We propose and study a new mathematical model of the human immunodeficiency virus (HIV). The main novelty is to consider that the antibody growth depends not only on the virus and on the antibodies concentration but also on the uninfected cells concentration. The model consists of five nonlinear differential equations describing the evolution of the uninfected cells, the infected ones, the free viruses, and the adaptive immunity. The adaptive immune response is represented by the cytotoxic T-lymphocytes (CTL) cells and the antibodies with the growth function supposed to be trilinear. The model includes two kinds of treatments. The objective of the first one is to reduce the number of infected cells, while the aim of the second is to block free viruses. Firstly, the positivity and the boundedness of solutions are established. After that, the local stability of the disease free steady state and the infection steady states are characterized. Next, an optimal control problem is posed and investigated. Finally, numerical simulations are performed in order to show the behavior of solutions and the effectiveness of the two incorporated treatments via an efficient optimal control strategy.
The abundance of a species population in an ecosystem is rarely stationary, often exhibiting large fluctuations over time. Using historical data on marine species, we show that the year-to-year fluctuations of population growth rate obey a well-defined double-exponential (Laplace) distribution. This striking regularity allows us to devise a stochastic model despite seemingly irregular variations in population abundances. The model identifies the effect of reduced growth at low population density as a key factor missed in current approaches of population variability analysis and without which extinction risks are severely underestimated. The model also allows us to separate the effect of demographic stochasticity and show that single-species growth rates are dominantly determined by stochasticity common to all species. This dominance---and the implications it has for interspecies correlations, including co-extinctions---emphasizes the need of ecosystem-level management approaches to reduce the extinction risk of the individual species themselves.