Demography and diffusion in epidemics: malaria and black death spread

Data-driven mathematical modeling approaches for COVID-19: A survey

In this review, we successively present the methods for phenomenological modeling of the evolution of reported and unreported cases of COVID-19, both in the exponential phase of growth and then in a complete epidemic wave. After the case of an isolated wave, we present the modeling of several successive waves separated by endemic stationary periods. Then, we treat the case of multi-compartmental models without or with age structure. Eventually, we review the literature, based on 260 articles selected in 11 sections, ranging from the medical survey of hospital cases to forecasting the dynamics of new cases in the general population. This review favors the phenomenological approach over the mechanistic approach in the choice of references and provides simulations of the evolution of the number of observed cases of COVID-19 for 10 states (California, China, France, India, Israel, Japan, New York, Peru, Spain and United Kingdom).

Mathematical Modelling of the Spatial Distribution of a COVID-19 Outbreak with Vaccination Using Diffusion Equation

The formulation of mathematical models using differential equations has become crucial in predicting the evolution of viral diseases in a population in order to take preventive and curative measures. In December 2019, a novel variety of Coronavirus (SARS-CoV-2) was identified in Wuhan, Hubei Province, China, which causes a severe and potentially fatal respiratory syndrome. Since then, it has been declared a pandemic by the World Health Organization and has spread around the globe. A reaction−diffusion system is a mathematical model that describes the evolution of a phenomenon subjected to two processes: a reaction process, in which different substances are transformed, and a diffusion process, which causes their distribution in space. This article provides a mathematical study of the Susceptible, Exposed, Infected, Recovered, and Vaccinated population model of the COVID-19 pandemic using the bias of reaction−diffusion equations. Both local and global asymptotic stability conditions for the equilibria were determined using a Lyapunov function, and the nature of the stability was determined using the Routh−Hurwitz criterion. Furthermore, we consider the conditions for the existence and uniqueness of the model solution and show the spatial distribution of the model compartments when the basic reproduction rate R0<1 and R0>1. Thereafter, we conducted a sensitivity analysis to determine the most sensitive parameters in the proposed model. We demonstrate the model’s effectiveness by performing numerical simulations and investigating the impact of vaccination, together with the significance of spatial distribution parameters in the spread of COVID-19. The findings indicate that reducing contact with an infected person and increasing the proportion of susceptible people who receive high-efficacy vaccination will lessen the burden of COVID-19 in the population. Therefore, we offer to the public health policymakers a better understanding of COVID-19 management.

Timing the race of vaccination, new variants, and relaxing restrictions during COVID-19 pandemic

Late in 2019, China identified a new type of coronavirus, SARS-CoV-2, and due to its fast spread, the World Health Organisation (WHO) declared a pandemic named COVID-19. Some variants of this virus were detected, including the Delta, which caused new waves of infections. This work uses an extended version of a SIRD model that includes vaccination effects to measure the impact of the Delta variant in three countries: Germany, Israel and Brazil. The calibrated models were able to reproduce the dynamics of the above countries. In addition, hypothetical scenarios were simulated to quantify the impact of vaccination and mitigation policies during the Delta wave. The results showed that the model could reproduce the complex dynamics observed in the different countries. The estimated increase of transmission rate due to the Delta variant was highest in Israel (7.9), followed by Germany (2.7) and Brazil (1.5). These values may support the hypothesis that people immunised against COVID-19 may lose their defensive antibodies with time since Israel, Germany, and Brazil fully vaccinated half of the population in March, July, and October. The scenario to study the impact of vaccination revealed relative reductions in the total number of deaths between 30% and 250%; an absolute reduction of 300 thousand deaths in Brazil due to vaccination during the Delta wave. The second hypothetical scenario revealed that mitigation policies saved up to 300 thousand Brazilians; relative reductions in the total number of deaths between 24% and 120% in the three analysed countries. Therefore, the results suggest that both vaccination and mitigation policies were crucial in decreasing the spread and the number of deaths during the Delta wave.

Characterisation of Omicron Variant during COVID-19 Pandemic and the Impact of Vaccination, Transmission Rate, Mortality, and Reinfection in South Africa, Germany, and Brazil

Several variants of SARS-CoV-2 have been identified in different parts of the world, including Gamma, detected in Brazil, Delta, detected in India, and the recent Omicron variant, detected in South Africa. The emergence of a new variant is a cause of great concern. This work considers an extended version of an SIRD model capable of incorporating the effects of vaccination, time-dependent transmissibility rates, mortality, and even potential reinfections during the pandemic. We use this model to characterise the Omicron wave in Brazil, South Africa, and Germany. During Omicron, the transmissibility increased by five for Brazil and Germany and eight for South Africa, whereas the estimated mortality was reduced by three-fold. We estimated that the reported cases accounted for less than 25% of the actual cases during Omicron. The mortality among the nonvaccinated population in these countries is, on average, three to four times higher than the mortality among the fully vaccinated. Finally, we could only reproduce the observed dynamics after introducing a new parameter that accounts for the percentage of the population that can be reinfected. Reinfection was as high as 40% in South Africa, which has only 29% of its population fully vaccinated and as low as 13% in Brazil, which has over 70% and 80% of its population fully vaccinated and with at least one dose, respectively. The calibrated models were able to estimate essential features of the complex virus and vaccination dynamics and stand as valuable tools for quantifying the impact of protocols and decisions in different populations.

Age Dependent Epidemic Modeling of COVID-19 Outbreak in Kuwait, France, and Cameroon

Revisiting the classical model by Ross and Kermack-McKendrick, the Susceptible−Infectious−Recovered (SIR) model used to formalize the COVID-19 epidemic, requires improvements which will be the subject of this article. The heterogeneity in the age of the populations concerned leads to considering models in age groups with specific susceptibilities, which makes the prediction problem more difficult. Basically, there are three age groups of interest which are, respectively, 0−19 years, 20−64 years, and >64 years, but in this article, we only consider two (20−64 years and >64 years) age groups because the group 0−19 years is widely seen as being less infected by the virus since this age group had a low infection rate throughout the pandemic era of this study, especially the countries under consideration. In this article, we proposed a new mathematical age-dependent (Susceptible−Infectious−Goneanewsusceptible−Recovered (SIGR)) model for the COVID-19 outbreak and performed some mathematical analyses by showing the positivity, boundedness, stability, existence, and uniqueness of the solution. We performed numerical simulations of the model with parameters from Kuwait, France, and Cameroon. We discuss the role of these different parameters used in the model; namely, vaccination on the epidemic dynamics. We open a new perspective of improving an age-dependent model and its application to observed data and parameters.

Estimation of Daily Reproduction Numbers during the COVID-19 Outbreak

(1) Background: The estimation of daily reproduction numbers throughout the contagiousness period is rarely considered, and only their sum R0 is calculated to quantify the contagiousness level of an infectious disease. (2) Methods: We provide the equation of the discrete dynamics of the epidemic’s growth and obtain an estimation of the daily reproduction numbers by using a deconvolution technique on a series of new COVID-19 cases. (3) Results: We provide both simulation results and estimations for several countries and waves of the COVID-19 outbreak. (4) Discussion: We discuss the role of noise on the stability of the epidemic’s dynamics. (5) Conclusions: We consider the possibility of improving the estimation of the distribution of daily reproduction numbers during the contagiousness period by taking into account the heterogeneity due to several host age classes.

Reaction-diffusion spatial modeling of COVID-19: Greece and Andalusia as case examples

We examine the spatial modeling of the outbreak of COVID-19 in two regions: the autonomous community of Andalusia in Spain and the mainland of Greece. We start with a zero-dimensional (0D; ordinary-differential-equation-level) compartmental epidemiological model consisting of Susceptible, Exposed, Asymptomatic, (symptomatically) Infected, Hospitalized, Recovered, and deceased populations (SEAIHR model). We emphasize the importance of the viral latent period (reflected in the exposed population) and the key role of an asymptomatic population. We optimize model parameters for both regions by comparing predictions to the cumulative number of infected and total number of deaths, the reported data we found to be most reliable, via minimizing the ℓ^{2} norm of the difference between predictions and observed data. We consider the sensitivity of model predictions on reasonable variations of model parameters and initial conditions, and we address issues of parameter identifiability. We model both the prequarantine and postquarantine evolution of the epidemic by a time-dependent change of the viral transmission rates that arises in response to containment measures. Subsequently, a spatially distributed version of the 0D model in the form of reaction-diffusion equations is developed. We consider that, after an initial localized seeding of the infection, its spread is governed by the diffusion (and 0D model "reactions") of the asymptomatic and symptomatically infected populations, which decrease with the imposed restrictive measures. We inserted the maps of the two regions, and we imported population-density data into the finite-element software package COMSOL Multiphysics®, which was subsequently used to numerically solve the model partial differential equations. Upon discussing how to adapt the 0D model to this spatial setting, we show that these models bear significant potential towards capturing both the well-mixed, zero-dimensional description and the spatial expansion of the pandemic in the two regions. Veins of potential refinement of the model assumptions towards future work are also explored.

A parsimonious approach for spatial transmission and heterogeneity in the COVID-19 propagation

Raw data on the number of deaths at a country level generally indicate a spatially variable distribution of COVID-19 incidence. An important issue is whether this pattern is a consequence of environmental heterogeneities, such as the climatic conditions, during the course of the outbreak. Another fundamental issue is to understand the spatial spreading of COVID-19. To address these questions, we consider four candidate epidemiological models with varying complexity in terms of initial conditions, contact rates and non-local transmissions, and we fit them to French mortality data with a mixed probabilistic-ODE approach. Using statistical criteria, we select the model with non-local transmission corresponding to a diffusion on the graph of counties that depends on the geographic proximity, with time-dependent contact rate and spatially constant parameters. This suggests that in a geographically middle size centralized country such as France, once the epidemic is established, the effect of global processes such as restriction policies and sanitary measures overwhelms the effect of local factors. Additionally, this approach reveals the latent epidemiological dynamics including the local level of immunity, and allows us to evaluate the role of non-local interactions on the future spread of the disease.

A Brief Theory of Epidemic Kinetics

In the context of the COVID-19 epidemic, and on the basis of the Theory of Dynamical Systems, we propose a simple theoretical approach for the expansion of contagious diseases, with a particular focus on viral respiratory tracts. The infection develops through contacts between contagious and exposed people, with a rate proportional to the number of contagious and of non-immune individuals, to contact duration and turnover, inversely proportional to the efficiency of protection measures, and balanced by the average individual recovery response. The obvious initial exponential increase is readily hindered by the growing recovery rate, and also by the size reduction of the exposed population. The system converges towards a stable attractor whose value is expressed in terms of the “reproductive rate” R0, depending on contamination and recovery factors. Various properties of the attractor are examined, and particularly its relations with R0. Decreasing this ratio below a critical value leads to a tipping threshold beyond which the epidemic is over. By contrast, significant values of the above ratio may bring the system through a bifurcating hierarchy of stable cycles up to a chaotic behaviour.

Temperature Decreases Spread Parameters of the New Covid-19 Case Dynamics

(1) Background: The virulence of coronavirus diseases due to viruses like SARS-CoV or MERS-CoV decreases in humid and hot weather. The putative temperature dependence of infectivity by the new coronavirus SARS-CoV-2 or covid-19 has a high predictive medical interest. (2) Methods: External temperature and new covid-19 cases in 21 countries and in the French administrative regions were collected from public data. Associations between epidemiological parameters of the new case dynamics and temperature were examined using an ARIMA model. (3) Results: We show that, in the first stages of the epidemic, the velocity of contagion decreases with country- or region-wise temperature. (4) Conclusions: Results indicate that high temperatures diminish initial contagion rates, but seasonal temperature effects at later stages of the epidemy remain questionable. Confinement policies and other eviction rules should account for climatological heterogeneities, in order to adapt the public health decisions to possible geographic or seasonal gradients.

CRISPR-based herd immunity can limit phage epidemics in bacterial populations

Herd immunity, a process in which resistant individuals limit the spread of a pathogen among susceptible hosts has been extensively studied in eukaryotes. Even though bacteria have evolved multiple immune systems against their phage pathogens, herd immunity in bacteria remains unexplored. Here we experimentally demonstrate that herd immunity arises during phage epidemics in structured and unstructured Escherichia coli populations consisting of differing frequencies of susceptible and resistant cells harboring CRISPR immunity. In addition, we develop a mathematical model that quantifies how herd immunity is affected by spatial population structure, bacterial growth rate, and phage replication rate. Using our model we infer a general epidemiological rule describing the relative speed of an epidemic in partially resistant spatially structured populations. Our experimental and theoretical findings indicate that herd immunity may be important in bacterial communities, allowing for stable coexistence of bacteria and their phages and the maintenance of polymorphism in bacterial immunity.

Mathematical models for predicting human mobility in the context of infectious disease spread: introducing the impedance model

Background

Mathematical models of human mobility have demonstrated a great potential for infectious disease epidemiology in contexts of data scarcity. While the commonly used gravity model involves parameter tuning and is thus difficult to implement without reference data, the more recent radiation model based on population densities is parameter-free, but biased. In this study we introduce the new impedance model, by analogy with electricity. Previous research has compared models on the basis of a few specific available spatial patterns. In this study, we use a systematic simulation-based approach to assess the performances.

Methods

Five hundred spatial patterns were generated using various area sizes and location coordinates. Model performances were evaluated based on these patterns. For simulated data, comparison measures were average root mean square error (aRMSE) and bias criteria. Modeling of the 2010 Haiti cholera epidemic with a basic susceptible–infected–recovered (SIR) framework allowed an empirical evaluation through assessing the goodness-of-fit of the observed epidemic curve.

Results

The new, parameter-free impedance model outperformed previous models on simulated data according to average aRMSE and bias criteria. The impedance model achieved better performances with heterogeneous population densities and small destination populations. As a proof of concept, the basic compartmental SIR framework was used to confirm the results obtained with the impedance model in predicting the spread of cholera in Haiti in 2010.

Conclusions

The proposed new impedance model provides accurate estimations of human mobility, especially when the population distribution is highly heterogeneous. This model can therefore help to achieve more accurate predictions of disease spread in the context of an epidemic.

Electronic supplementary material

The online version of this article (10.1186/s12942-017-0115-7) contains supplementary material, which is available to authorized users.

Discrete dynamics of contagious social diseases: Example of obesity

Modeling contagious diseases needs to incorporate information about social networks through which the disease spreads as well as data about demographic and genetic changes in the susceptible population. In this paper, we propose a theoretical framework (conceptualization and formalization) which seeks to model obesity as a process of transformation of one's own body determined by individual (physical and psychological), inter-individual (relational, i.e., relative to the relationship between the individual and others) and socio-cultural (environmental, i.e., relative to the relationship between the individual and his milieu) factors. Individual and inter-individual factors are tied to each other in a socio-cultural context whose impact is notably related to the visibility of anybody being exposed on the public stage in a non-contingent way. The question we are dealing with in this article is whether such kind of social diseases, i.e., depending upon socio-environmental exposure, can be considered as "contagious". In other words, can obesity be propagated from individual to individual or from environmental sources throughout an entire population?

Mathematical modeling of infectious disease dynamics

Over the last years, an intensive worldwide effort is speeding up the developments in the establishment of a global surveillance network for combating pandemics of emergent and re-emergent infectious diseases. Scientists from different fields extending from medicine and molecular biology to computer science and applied mathematics have teamed up for rapid assessment of potentially urgent situations. Toward this aim mathematical modeling plays an important role in efforts that focus on predicting, assessing, and controlling potential outbreaks. To better understand and model the contagious dynamics the impact of numerous variables ranging from the micro host–pathogen level to host-to-host interactions, as well as prevailing ecological, social, economic, and demographic factors across the globe have to be analyzed and thoroughly studied. Here, we present and discuss the main approaches that are used for the surveillance and modeling of infectious disease dynamics. We present the basic concepts underpinning their implementation and practice and for each category we give an annotated list of representative works.

Random modelling of contagious diseases

Modelling contagious diseases needs to include a mechanistic knowledge about contacts between hosts and pathogens as specific as possible, e.g., by incorporating in the model information about social networks through which the disease spreads. The unknown part concerning the contact mechanism can be modelled using a stochastic approach. For that purpose, we revisit SIR models by introducing first a microscopic stochastic version of the contacts between individuals of different populations (namely Susceptible, Infective and Recovering), then by adding a random perturbation in the vicinity of the endemic fixed point of the SIR model and eventually by introducing the definition of various types of random social networks. We propose as example of application to contagious diseases the HIV, and we show that a micro-simulation of individual based modelling (IBM) type can reproduce the current stable incidence of the HIV epidemic in a population of HIV-positive men having sex with men (MSM).