Impact of reproduction number on the multiwave spreading dynamics of COVID-19 with temporary immunity: A mathematical model.
Ontology highlight
ABSTRACT: Objectives
The recent discoveries of phylogenetically confirmed COVID-19 reinfection cases worldwide, together with studies suggesting that antibody titres decrease over time, raise the question of what course the epidemic trajectories may take if immunity were really to be temporary in a significant fraction of the population. The objective of this study is to obtain an answer for this important question.Methods
We construct a ground-up delay differential equation model tailored to incorporate different types of immune response. We considered two immune responses: (a) short-lived immunity of all types, and (b) short-lived sterilizing immunity with durable severity-reducing immunity.Results
Multiple wave solutions to the model are manifest for intermediate values of the reproduction number R; interestingly, for sufficiently low as well as sufficiently high R, we find conventional single-wave solutions despite temporary immunity.Conclusions
The versatility of our model, and its very modest demands on computational resources, ensure that a set of disease trajectories can be computed virtually on the same day that a new and relevant immune response study is released. Our work can also be used to analyse the disease dynamics after a vaccine is certified for use and information regarding its immune response becomes available.
SUBMITTER: Shayak B
PROVIDER: S-EPMC7836705 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:The SIR ('susceptible-infectious-recovered') formulation is used to uncover the generic spread mechanisms observed by COVID-19 dynamics globally, especially in the early phases of infectious spread. During this early period, potential controls were not effectively put in place or enforced in many countries. Hence, the early phases of COVID-19 spread in countries where controls were weak offer a unique perspective on the ensemble-behavior of COVID-19 basic reproduction number Ro inferred from SIR formulation. The work here shows that there is global convergence (i.e., across many nations) to an uncontrolled Ro = 4.5 that describes the early time spread of COVID-19. This value is in agreement with independent estimates from other sources reviewed here and adds to the growing consensus that the early estimate of Ro = 2.2 adopted by the World Health Organization is low. A reconciliation between power-law and exponential growth predictions is also featured within the confines of the SIR formulation. The effects of testing ramp-up and the role of 'super-spreaders' on the inference of Ro are analyzed using idealized scenarios. Implications for evaluating potential control strategies from this uncontrolled Ro are briefly discussed in the context of the maximum possible infected fraction of the population (needed to assess health care capacity) and mortality (especially in the USA given diverging projections). Model results indicate that if intervention measures still result in Ro > 2.7 within 44 days after first infection, intervention is unlikely to be effective in general for COVID-19.
| S-EPMC7514051 | biostudies-literature
Project description:To monitor local and global COVID-19 outbreaks, and to plan containment measures, accessible and comprehensible decision-making tools need to be based on the growth rates of new confirmed infections, hospitalization or case fatality rates. Growth rates of new cases form the empirical basis for estimates of a variety of reproduction numbers, dimensionless numbers whose value, when larger than unity, describes surging infections and generally worsening epidemiological conditions. Typically, these determinations rely on noisy or incomplete data gained over limited periods of time, and on many parameters to estimate. This paper examines how estimates from data and models of time-evolving reproduction numbers of national COVID-19 infection spread change by using different techniques and assumptions. Given the importance acquired by reproduction numbers as diagnostic tools, assessing their range of possible variations obtainable from the same epidemiological data is relevant. We compute control reproduction numbers from Swiss and Italian COVID-19 time series adopting both data convolution (renewal equation) and a SEIR-type model. Within these two paradigms we run a comparative analysis of the possible inferences obtained through approximations of the distributions typically used to describe serial intervals, generation, latency and incubation times, and the delays between onset of symptoms and notification. Our results suggest that estimates of reproduction numbers under these different assumptions may show significant temporal differences, while the actual variability range of computed values is rather small.
| S-EPMC7723757 | biostudies-literature
Project description:Early estimates of the transmission potential of emerging and re-emerging infections are increasingly used to inform public health authorities on the level of risk posed by outbreaks. Existing methods to estimate the reproduction number generally assume exponential growth in case incidence in the first few disease generations, before susceptible depletion sets in. In reality, outbreaks can display subexponential (i.e. polynomial) growth in the first few disease generations, owing to clustering in contact patterns, spatial effects, inhomogeneous mixing, reactive behaviour changes or other mechanisms. Here, we introduce the generalized growth model to characterize the early growth profile of outbreaks and estimate the effective reproduction number, with no need for explicit assumptions about the shape of epidemic growth. We demonstrate this phenomenological approach using analytical results and simulations from mechanistic models, and provide validation against a range of empirical disease datasets. Our results suggest that subexponential growth in the early phase of an epidemic is the rule rather the exception. Mechanistic simulations show that slight modifications to the classical susceptible-infectious-removed model result in subexponential growth, and in turn a rapid decline in the reproduction number within three to five disease generations. For empirical outbreaks, the generalized-growth model consistently outperforms the exponential model for a variety of directly and indirectly transmitted diseases datasets (pandemic influenza, measles, smallpox, bubonic plague, cholera, foot-and-mouth disease, HIV/AIDS and Ebola) with model estimates supporting subexponential growth dynamics. The rapid decline in effective reproduction number predicted by analytical results and observed in real and synthetic datasets within three to five disease generations contrasts with the expectation of invariant reproduction number in epidemics obeying exponential growth. The generalized-growth concept also provides us a compelling argument for the unexpected extinction of certain emerging disease outbreaks during the early ascending phase. Overall, our approach promotes a more reliable and data-driven characterization of the early epidemic phase, which is important for accurate estimation of the reproduction number and prediction of disease impact.
| S-EPMC5095223 | biostudies-literature
Project description:Keratin intermediate filaments (IFs) form cross-linked arrays to fulfill their structural support function in epithelial cells and tissues subjected to external stress. How the cross-linking of keratin IFs impacts the morphology and differentiation of keratinocytes in the epidermis and related surface epithelia remains an open question. Experimental measurements have established that keratinocyte spreading area is inversely correlated to the extent of keratin IF bundling in two-dimensional culture. In an effort to quantitatively explain this relationship, we developed a mathematical model in which isotropic cell spreading is considered as a first approximation. Relevant physical properties such as actin protrusion, adhesion events, and the corresponding response of lamellum formation at the cell periphery are included in this model. Through optimization with experimental data that relate time-dependent changes in keratinocyte surface area during spreading, our simulation results confirm the notion that the organization and mechanical properties of cross-linked keratin filaments affect cell spreading; in addition, our results provide details of the kinetics of this effect. These in silico findings provide further support for the notion that differentiation-related changes in the density and intracellular organization of keratin IFs affect tissue architecture in epidermis and related stratified epithelia.
| S-EPMC3491728 | biostudies-literature
Project description:BACKGROUND:COVID-19 disease is becoming a global pandemic and more than 200 countries were affected because of this disease. Italy is one of the countries is largely suffered with this virus outbreak, and about 180,000 cases (as of 20 April 2020) were registered which explains the large transmissibility and reproduction case numbers. OBJECTIVE:In this study, we considered the Marche region of Italy to compute different daily transmission rates (Rt) including five provinces in it. We also present forecasting of daily and cumulative incidences associated after the next thirty days. The Marche region is the 8th in terms of number of people infected in Italy and the first in terms of diffusion of the infection among the 4 regions of the center of Italy. METHODS:Epidemic statistics were extracted from the national Italian Health Ministry website. We considered outbreak information where the first case registered in Marche with onset symptoms (26 February 2020) to the present date (20 April 2020). Adoption of incidence and projections with R statistics was done. RESULTS:The median values of Rt for the five provinces of Pesaro and Urbano, Ancona, Fermo, Ascoli Piceno, and Macerata, was 2.492 (1.1-4.5), 2.162 (1.0-4.0), 1.512 (0.75-2.75), 1.141 (1.0-1.6), and 1.792 (1.0-3.5) with 95% of CI achieved. The projections at end of 30th day of the cumulative incidences 323 (95% CI), and daily incidences 45 (95% CI) could be possible. CONCLUSIONS:This study highlights the knowledge of essential insights into the Marche region in particular to virus transmission dynamics, geographical characteristics of positive incidences, and the necessity of implementing mitigation procedures to fight against the COVID-19 outbreak.
| S-EPMC7211603 | biostudies-literature
Project description:ObjectiveTo estimate the basic reproduction number (R0) for COVID-19 in Western Europe.MethodsData (official statistics) on the cumulative incidence of COVID-19 at the start of the outbreak (before any confinement rules were declared) were retrieved in the 15 largest countries in Western Europe, allowing us to estimate the exponential growth rate of the disease. The rate was then combined with estimates of the distribution of the generation interval as reconstructed from the literature.ResultsDespite the possible unreliability of some official statistics about COVID-19, the spread of the disease appears to be remarkably similar in most European countries, allowing us to estimate an average R0 in Western Europe of 2.2 (95% CI: 1.9-2.6).ConclusionsThe value of R0 for COVID-19 in Western Europe appears to be significantly lower than that in China. The proportion of immune persons in the European population required to stop the outbreak could thus be closer to 50% than to 70%.
| S-EPMC7968714 | biostudies-literature
Project description:Background: The current propagation models of COVID-19 are poorly consistent with existing epidemiological data and with evidence that the SARS-CoV-2 genome is mutating, for potential aggressive evolution of the disease. Objectives: We looked for fundamental variables that were missing from current analyses. Among them were regional climate heterogeneity, viral evolution processes versus founder effects, and large-scale virus containment measures. Methods: We challenged regional versus genetic evolution models of COVID-19 at a whole-population level, over 168,089 laboratory-confirmed SARS-CoV-2 infection cases in Italy, Spain, and Scandinavia at early time-points of the pandemic. Diffusion data in Germany, France, and the United Kingdom provided a validation dataset of 210,239 additional cases. Results: Mean doubling time of COVID-19 cases was 6.63 days in Northern versus 5.38 days in Southern Italy. Spain extended this trend of faster diffusion in Southern Europe, with a doubling time of 4.2 days. Slower doubling times were observed in Sweden (9.4 days), Finland (10.8 days), and Norway (12.95 days). COVID-19 doubling time in Germany (7.0 days), France (7.5 days), and the United Kingdom (7.2 days) supported the North/South gradient model. Clusters of SARS-CoV-2 mutations upon sequential diffusion were not found to clearly correlate with regional distribution dynamics. Conclusion: Acquisition of mutations upon SARS-CoV-2 spreading failed to explain regional diffusion heterogeneity at early pandemic times. Our findings indicate that COVID-19 transmission rates are rather associated with a sharp North/South climate gradient, with faster spreading in Southern regions. Thus, warmer climate conditions may not limit SARS-CoV-2 infectivity. Very cold regions may be better spared by recurrent courses of SARS-CoV-2 infection.
| S-EPMC8740632 | biostudies-literature
Project description:UnlabelledUnderstanding the transmission dynamics of pathogens is essential to determine the epidemiology, ecology, and ways of controlling enterohemorrhagic Escherichia coli (EHEC) in animals and their environments. Our objective was to estimate the epidemiological fitness of common EHEC strains in cattle populations. For that purpose, we developed a Markov chain model to characterize the dynamics of 7 serogroups of enterohemorrhagic Escherichia coli (O26, O45, O103, O111, O121, O145, and O157) in cattle production environments based on a set of cross-sectional data on infection prevalence in 2 years in two U.S. states. The basic reproduction number (R0) was estimated using a Bayesian framework for each serogroup based on two criteria (using serogroup alone [the O-group data] and using O serogroup, Shiga toxin gene[s], and intimin [eae] gene together [the EHEC data]). In addition, correlations between external covariates (e.g., location, ambient temperature, dietary, and probiotic usage) and prevalence/R0 were quantified. R0 estimates varied substantially among different EHEC serogroups, with EHEC O157 having an R0 of >1 (?1.5) and all six other EHEC serogroups having an R0 of less than 1. Using the O-group data substantially increased R0 estimates for the O26, O45, and O103 serogroups (R0 > 1) but not for the others. Different covariates had distinct influences on different serogroups: the coefficients for each covariate were different among serogroups. Our modeling and analysis of this system can be readily expanded to other pathogen systems in order to estimate the pathogen and external factors that influence spread of infectious agents.ImportanceIn this paper we describe a Bayesian modeling framework to estimate basic reproduction numbers of multiple serotypes of Shiga toxin-producing Escherichia coli according to a cross-sectional study. We then coupled a compartmental model to reconstruct the infection dynamics of these serotypes and quantify their risk in the population. We incorporated different sensitivity levels of detecting different serotypes and evaluated their potential influence on the estimation of basic reproduction numbers.
| S-EPMC5007764 | biostudies-literature
Project description:Currently, investigations are intensively conducted on modeling, forecasting, and studying the dynamic spread of coronavirus (Covid-19) new pandemic. In the present work, the sigmoidal-Boltzmann mathematical model was applied to study the Covid-19 spread in 15 different countries. The cumulative number of infected persons I has been accurately fitted by the sigmoidal-Boltzmann equation (SBE), giving rise to different epidemiological parameters such as the pandemic peak tp, the maximum number of infected persons Imax, and the time of the epidemic stabilization t?. The time constant relative to the sigmoid ?t (called also the slope factor) was revealed to be the determining parameter which influences all the epidemiological parameters. Empirical laws between the different parameters allowed us to propose a modified sigmoidal-Boltzmann equation describing the spread of the pandemic. The expression of the spread speed Vp was further determined as a function of the sigmoid parameters. This made it possible to assess the maximum speed of spread of the virus Vpmax and to trace the speed profile in each country. In addition, for countries undergoing a second pandemic wave, the cumulative number of infected people I has been successfully adjusted by a double sigmoidal-Boltzmann equation (DSBE) allowing the comparison between the two waves. Finally, the comparison between the maximum virus spread of two waves Vp?max?1 and Vp?max?2 showed that the intensity of the second wave of Covid-19 is low compared to the first for all the countries studied.
| S-EPMC7557153 | biostudies-literature
Project description:BACKGROUND:The 2019 novel coronavirus (COVID-19) outbreak in Wuhan, China has attracted world-wide attention. As of March 31, 2020, a total of 82,631 cases of COVID-19 in China were confirmed by the National Health Commission (NHC) of China. METHODS:Three approaches, namely Poisson likelihood-based method (ML), exponential growth rate-based method (EGR) and stochastic Susceptible-Infected-Removed dynamic model-based method (SIR), were implemented to estimate the basic and controlled reproduction numbers. RESULTS:A total of 198 chains of transmission together with dates of symptoms onset and 139 dates of infections were identified among 14,829 confirmed cases outside Hubei Province as reported as of March 31, 2020. Based on this information, we found that the serial interval had an average of 4.60 days with a standard deviation of 5.55 days, the incubation period had an average of 8.00 days with a standard deviation of 4.75 days and the infectious period had an average of 13.96 days with a standard deviation of 5.20 days. The estimated controlled reproduction numbers, Rc, produced by all three methods in all analyzed regions of China are significantly smaller compared with the basic reproduction numbers R0. CONCLUSIONS:The controlled reproduction number in China is much lower than one in all regions of China by now. It fell below one within 30 days from the implementations of unprecedent containment measures, which indicates that the strong measures taken by China government was effective to contain the epidemic. Nonetheless, efforts are still needed in order to end the current epidemic as imported cases from overseas pose a high risk of a second outbreak.
| S-EPMC7211652 | biostudies-literature