The probability of epidemic burnout in the stochastic SIR model with vital dynamics.

We present an approach to computing the probability of epidemic "burnout," i.e., the probability that a newly emergent pathogen will go extinct after a major epidemic. Our analysis is based on the standard stochastic formulation of the Susceptible-Infectious-Removed (SIR) epidemic model including host demography (births and deaths) and corresponds to the standard SIR ordinary differential equations (ODEs) in the infinite population limit. Exploiting a boundary layer approximation to the ODEs and a birth-death process approximation to the stochastic dynamics within the boundary layer, we derive convenient, fully analytical approximations for the burnout probability. We demonstrate-by comparing with computationally demanding individual-based stochastic simulations and with semi-analytical approximations derived previously-that our fully analytical approximations are highly accurate for biologically plausible parameters. We show that the probability of burnout always decreases with increased mean infectious period. However, for typical biological parameters, there is a relevant local minimum in the probability of persistence as a function of the basic reproduction number [Formula: see text]. For the shortest infectious periods, persistence is least likely if [Formula: see text]; for longer infectious periods, the minimum point decreases to [Formula: see text]. For typical acute immunizing infections in human populations of realistic size, our analysis of the SIR model shows that burnout is almost certain in a well-mixed population, implying that susceptible recruitment through births is insufficient on its own to explain disease persistence.

Inferring the differences in incubation-period and generation-interval distributions of the Delta and Omicron variants of SARS-CoV-2.

Estimating the differences in the incubation-period, serial-interval, and generation-interval distributions of SARS-CoV-2 variants is critical to understanding their transmission. However, the impact of epidemic dynamics is often neglected in estimating the timing of infection-for example, when an epidemic is growing exponentially, a cohort of infected individuals who developed symptoms at the same time are more likely to have been infected recently. Here, we reanalyze incubation-period and serial-interval data describing transmissions of the Delta and Omicron variants from the Netherlands at the end of December 2021. Previous analysis of the same dataset reported shorter mean observed incubation period (3.2 d vs. 4.4 d) and serial interval (3.5 d vs. 4.1 d) for the Omicron variant, but the number of infections caused by the Delta variant decreased during this period as the number of Omicron infections increased. When we account for growth-rate differences of two variants during the study period, we estimate similar mean incubation periods (3.8 to 4.5 d) for both variants but a shorter mean generation interval for the Omicron variant (3.0 d; 95% CI: 2.7 to 3.2 d) than for the Delta variant (3.8 d; 95% CI: 3.7 to 4.0 d). The differences in estimated generation intervals may be driven by the "network effect"-higher effective transmissibility of the Omicron variant can cause faster susceptible depletion among contact networks, which in turn prevents late transmission (therefore shortening realized generation intervals). Using up-to-date generation-interval distributions is critical to accurately estimating the reproduction advantage of the Omicron variant.

Incubation Period and Serial Interval of Mpox in 2022 Global Outbreak Compared with Historical Estimates.

Understanding changes in the transmission dynamics of mpox requires comparing recent estimates of key epidemiologic parameters with historical data. We derived historical estimates for the incubation period and serial interval for mpox and contrasted them with pooled estimates from the 2022 outbreak. Our findings show the pooled mean infection-to-onset incubation period was 8.1 days for the 2022 outbreak and 8.2 days historically, indicating the incubation periods remained relatively consistent over time, despite a shift in the major mode of transmission. However, we estimated the onset-to-onset serial interval at 8.7 days using 2022 data, compared with 14.2 days using historical data. Although the reason for this shortening of the serial interval is unclear, it may be because of increased public health interventions or a shift in the mode of transmission. Recognizing such temporal shifts is essential for informed response strategies, and public health measures remain crucial for controlling mpox and similar future outbreaks.

Forward-looking serial intervals correctly link epidemic growth to reproduction numbers.

The reproduction number R and the growth rate r are critical epidemiological quantities. They are linked by generation intervals, the time between infection and onward transmission. Because generation intervals are difficult to observe, epidemiologists often substitute serial intervals, the time between symptom onset in successive links in a transmission chain. Recent studies suggest that such substitution biases estimates of R based on r. Here we explore how these intervals vary over the course of an epidemic, and the implications for R estimation. Forward-looking serial intervals, measuring time forward from symptom onset of an infector, correctly describe the renewal process of symptomatic cases and therefore reliably link R with r. In contrast, backward-looking intervals, which measure time backward, and intrinsic intervals, which neglect population-level dynamics, give incorrect R estimates. Forward-looking intervals are affected both by epidemic dynamics and by censoring, changing in complex ways over the course of an epidemic. We present a heuristic method for addressing biases that arise from neglecting changes in serial intervals. We apply the method to early (21 January to February 8, 2020) serial interval-based estimates of R for the COVID-19 outbreak in China outside Hubei province; using improperly defined serial intervals in this context biases estimates of initial R by up to a factor of 2.6. This study demonstrates the importance of early contact tracing efforts and provides a framework for reassessing generation intervals, serial intervals, and R estimates for COVID-19.

Acceleration of plague outbreaks in the second pandemic.

Historical records reveal the temporal patterns of a sequence of plague epidemics in London, United Kingdom, from the 14th to 17th centuries. Analysis of these records shows that later epidemics spread significantly faster ("accelerated"). Between the Black Death of 1348 and the later epidemics that culminated with the Great Plague of 1665, we estimate that the epidemic growth rate increased fourfold. Currently available data do not provide enough information to infer the mode of plague transmission in any given epidemic; nevertheless, order-of-magnitude estimates of epidemic parameters suggest that the observed slow growth rates in the 14th century are inconsistent with direct (pneumonic) transmission. We discuss the potential roles of demographic and ecological factors, such as climate change or human or rat population density, in driving the observed acceleration.

Awareness-driven behavior changes can shift the shape of epidemics away from peaks and toward plateaus, shoulders, and oscillations.

The COVID-19 pandemic has caused more than 1,000,000 reported deaths globally, of which more than 200,000 have been reported in the United States as of October 1, 2020. Public health interventions have had significant impacts in reducing transmission and in averting even more deaths. Nonetheless, in many jurisdictions, the decline of cases and fatalities after apparent epidemic peaks has not been rapid. Instead, the asymmetric decline in cases appears, in most cases, to be consistent with plateau- or shoulder-like phenomena-a qualitative observation reinforced by a symmetry analysis of US state-level fatality data. Here we explore a model of fatality-driven awareness in which individual protective measures increase with death rates. In this model, fast increases to the peak are often followed by plateaus, shoulders, and lag-driven oscillations. The asymmetric shape of model-predicted incidence and fatality curves is consistent with observations from many jurisdictions. Yet, in contrast to model predictions, we find that population-level mobility metrics usually increased from low levels before fatalities reached an initial peak. We show that incorporating fatigue and long-term behavior change can reconcile the apparent premature relaxation of mobility reductions and help understand when post-peak dynamics are likely to lead to a resurgence of cases.

Testing and Isolation Efficacy: Insights from a Simple Epidemic Model.

Testing individuals for pathogens can affect the spread of epidemics. Understanding how individual-level processes of sampling and reporting test results can affect community- or population-level spread is a dynamical modeling question. The effect of testing processes on epidemic dynamics depends on factors underlying implementation, particularly testing intensity and on whom testing is focused. Here, we use a simple model to explore how the individual-level effects of testing might directly impact population-level spread. Our model development was motivated by the COVID-19 epidemic, but has generic epidemiological and testing structures. To the classic SIR framework we have added a per capita testing intensity, and compartment-specific testing weights, which can be adjusted to reflect different testing emphases-surveillance, diagnosis, or control. We derive an analytic expression for the relative reduction in the basic reproductive number due to testing, test-reporting and related isolation behaviours. Intensive testing and fast test reporting are expected to be beneficial at the community level because they can provide a rapid assessment of the situation, identify hot spots, and may enable rapid contact-tracing. Direct effects of fast testing at the individual level are less clear, and may depend on how individuals' behaviour is affected by testing information. Our simple model shows that under some circumstances both increased testing intensity and faster test reporting can reduce the effectiveness of control, and allows us to explore the conditions under which this occurs. Conversely, we find that focusing testing on infected individuals always acts to increase effectiveness of control.

From science to politics: COVID-19 information fatigue on YouTube.

OBJECTIVE: The COVID-19 pandemic is the first pandemic where social media platforms relayed information on a large scale, enabling an "infodemic" of conflicting information which undermined the global response to the pandemic. Understanding how the information circulated and evolved on social media platforms is essential for planning future public health campaigns. This study investigated what types of themes about COVID-19 were most viewed on YouTube during the first 8 months of the pandemic, and how COVID-19 themes progressed over this period. METHODS: We analyzed top-viewed YouTube COVID-19-related videos in English from December 1, 2019 to August 16, 2020 with an open inductive content analysis. We coded 536 videos associated with 1.1 billion views across the study period. East Asian countries were the first to report the virus, while most of the top-viewed videos in English were from the US. Videos from straight news outlets dominated the top-viewed videos throughout the outbreak, and public health authorities contributed the fewest. Although straight news was the dominant COVID-19 video source with various types of themes, its viewership per video was similar to that for entertainment news and YouTubers after March. RESULTS: We found, first, that collective public attention to the COVID-19 pandemic on YouTube peaked around March 2020, before the outbreak peaked, and flattened afterwards despite a spike in worldwide cases. Second, more videos focused on prevention early on, but videos with political themes increased through time. Third, regarding prevention and control measures, masking received much less attention than lockdown and social distancing in the study period. CONCLUSION: Our study suggests that a transition of focus from science to politics on social media intensified the COVID-19 infodemic and may have weakened mitigation measures during the first waves of the COVID-19 pandemic. It is recommended that authorities should consider co-operating with reputable social media influencers to promote health campaigns and improve health literacy. In addition, given high levels of globalization of social platforms and polarization of users, tailoring communication towards different digital communities is likely to be essential.

Speed and strength of an epidemic intervention.

An epidemic can be characterized by its strength (i.e., the reproductive number [Formula: see text]) and speed (i.e., the exponential growth rate r). Disease modellers have historically placed much more emphasis on strength, in part because the effectiveness of an intervention strategy is typically evaluated on this scale. Here, we develop a mathematical framework for the classic, strength-based paradigm and show that there is a dual speed-based paradigm which can provide complementary insights. In particular, we note that r = 0 is a threshold for disease spread, just like [Formula: see text] [ 1], and show that we can measure the strength and speed of an intervention on the same scale as the strength and speed of an epidemic, respectively. We argue that, while the strength-based paradigm provides the clearest insight into certain questions, the speed-based paradigm provides the clearest view in other cases. As an example, we show that evaluating the prospects of 'test-and-treat' interventions against the human immunodeficiency virus (HIV) can be done more clearly on the speed than strength scale, given uncertainty in the proportion of HIV spread that happens early in the course of infection. We also discuss evaluating the effects of the importance of pre-symptomatic transmission of the SARS-CoV-2 virus. We suggest that disease modellers should avoid over-emphasizing the reproductive number at the expense of the exponential growth rate, but instead look at these as complementary measures.

Calibration of individual-based models to epidemiological data: A systematic review.

Individual-based models (IBMs) informing public health policy should be calibrated to data and provide estimates of uncertainty. Two main components of model-calibration methods are the parameter-search strategy and the goodness-of-fit (GOF) measure; many options exist for each of these. This review provides an overview of calibration methods used in IBMs modelling infectious disease spread. We identified articles on PubMed employing simulation-based methods to calibrate IBMs informing public health policy in HIV, tuberculosis, and malaria epidemiology published between 1 January 2013 and 31 December 2018. Articles were included if models stored individual-specific information, and calibration involved comparing model output to population-level targets. We extracted information on parameter-search strategies, GOF measures, and model validation. The PubMed search identified 653 candidate articles, of which 84 met the review criteria. Of the included articles, 40 (48%) combined a quantitative GOF measure with an algorithmic parameter-search strategy-either an optimisation algorithm (14/40) or a sampling algorithm (26/40). These 40 articles varied widely in their choices of parameter-search strategies and GOF measures. For the remaining 44 (52%) articles, the parameter-search strategy could either not be identified (32/44) or was described as an informal, non-reproducible method (12/44). Of these 44 articles, the majority (25/44) were unclear about the GOF measure used; of the rest, only five quantitatively evaluated GOF. Only a minority of the included articles, 14 (17%) provided a rationale for their choice of model-calibration method. Model validation was reported in 31 (37%) articles. Reporting on calibration methods is far from optimal in epidemiological modelling studies of HIV, malaria and TB transmission dynamics. The adoption of better documented, algorithmic calibration methods could improve both reproducibility and the quality of inference in model-based epidemiology. There is a need for research comparing the performance of calibration methods to inform decisions about the parameter-search strategies and GOF measures.

Does Counting Different Life Stages Impact Estimates for Extinction Probabilities for Tsetse (Glossina spp)?

As insect populations decline, due to climate change and other environmental disruptions, there has been an increased interest in understanding extinction probabilities. Generally, the life cycle of insects occurs in well-defined stages: when counting insects, questions naturally arise about which life stage to count. Using tsetse flies (vectors of trypanosomiasis) as a case study, we develop a model that works when different life stages are counted. Previous branching process models for tsetse populations only explicitly represent newly emerged adult female tsetse and use that subpopulation to keep track of population growth/decline. Here, we directly model other life stages. We analyse reproduction numbers and extinction probabilities and show that several previous models used for estimating extinction probabilities for tsetse populations are special cases of the current model. We confirm that the reproduction number is the same regardless of which life stage is counted, and show how the extinction probability depends on which life stage we start from. We demonstrate, and provide a biological explanation for, a simple relationship between extinction probabilities for the different life stages, based on the probability of recruitment between stages. These results offer insights into insect population dynamics and provide tools that will help with more detailed models of tsetse populations. Population dynamics studies of insects should be clear about life stages and counting points.

Age-dependence of healthcare interventions for COVID-19 in Ontario, Canada.

BACKGROUND: Patient age is one of the most salient clinical indicators of risk from COVID-19. Age-specific distributions of known SARS-CoV-2 infections and COVID-19-related deaths are available for many regions. Less attention has been given to the age distributions of serious medical interventions administered to COVID-19 patients, which could reveal sources of potential pressure on the healthcare system should SARS-CoV-2 prevalence increase, and could inform mass vaccination strategies. The aim of this study is to quantify the relationship between COVID-19 patient age and serious outcomes of the disease, beyond fatalities alone. METHODS: We analysed 277,555 known SARS-CoV-2 infection records for Ontario, Canada, from 23 January 2020 to 16 February 2021 and estimated the age distributions of hospitalizations, Intensive Care Unit admissions, intubations, and ventilations. We quantified the probability of hospitalization given known SARS-CoV-2 infection, and of survival given COVID-19-related hospitalization. RESULTS: The distribution of hospitalizations peaks with a wide plateau covering ages 60-90, whereas deaths are concentrated in ages 80+. The estimated probability of hospitalization given known infection reaches a maximum of 27.8% at age 80 (95% CI 26.0%-29.7%). The probability of survival given hospitalization is nearly 100% for adults younger than 40, but declines substantially after this age; for example, a hospitalized 54-year-old patient has a 91.7% chance of surviving COVID-19 (95% CI 88.3%-94.4%). CONCLUSIONS: Our study demonstrates a significant need for hospitalization in middle-aged individuals and young seniors. This need is not captured by the distribution of deaths, which is heavily concentrated in very old ages. The probability of survival given hospitalization for COVID-19 is lower than is generally perceived for patients over 40. If acute care capacity is exceeded due to an increase in COVID-19 prevalence, the distribution of deaths could expand toward younger ages. These results suggest that vaccine programs should aim to prevent infection not only in old seniors, but also in young seniors and middle-aged individuals, to protect them from serious illness and to limit stress on the healthcare system.

Potential Role of Social Distancing in Mitigating Spread of Coronavirus Disease, South Korea.

In South Korea, the coronavirus disease outbreak peaked at the end of February and subsided in mid-March. We analyzed the likely roles of social distancing in reducing transmission. Our analysis indicated that although transmission might persist in some regions, epidemics can be suppressed with less extreme measures than those taken by China.

Traditional Male Circumcision is Associated with Sexual Risk Behaviors in Sub-Saharan Countries Prioritized for Male Circumcision.

To understand the sexual risk behavior of men with traditional male circumcision and medical male circumcision in the context of the World Health Organization's (WHO) campaign for voluntary medical male circumcision (VMMC) scale-up, we investigated ten countries prioritized for the scale-up from the Demographic and Health Surveys. Male respondents aged 15-49 were selected. Ordinal regression was used to analyze the relationship between three sexual risk behaviors-condom use with non-cohabiting partners, number of non-cohabiting partners, and partner type-and circumcision status (traditionally circumcised before and after the VMMC scale-up, medically circumcised before and after the scale-up, and not circumcised), while controlling for social demographic covariates. We found evidence that some sexual risky behavior, specifically lower condom use and higher number of sexual partners, was associated with traditional circumcision. This finding suggests that messages about the protective effect of male circumcision may not have reached men with traditional circumcision. We suggest that WHO's VMMC campaign should include communities where traditional male circumcision is popular. We looked for, but did not find, evidence of differences between groups circumcised at different times, which could have indicated sexual risk compensation.

Patterns of seasonal and pandemic influenza-associated health care and mortality in Ontario, Canada.

BACKGROUND: Mathematical and statistical models are used to project the future time course of infectious disease epidemics and the expected future burden on health care systems and economies. Influenza is a particularly important disease in this context because it causes annual epidemics and occasional pandemics. In order to forecast health care utilization during epidemics-and the effects of hospitalizations and deaths on the contact network and, in turn, on transmission dynamics-modellers must make assumptions about the lengths of time between infection, visiting a physician, being admitted to hospital, leaving hospital, and death. More reliable forecasts could be be made if the distributions of times between these types of events ("delay distributions") were known. METHODS: We estimated delay distributions in the province of Ontario, Canada, between 2006 and 2010. To do so, we used encrypted health insurance numbers to link 1.34 billion health care billing records to 4.27 million hospital inpatient stays. Because the four year period we studied included three typical influenza seasons and the 2009 influenza pandemic, we were able to compare the delay distributions in non-pandemic and pandemic settings. We also estimated conditional probabilities such as the probability of hospitalization within the year if diagnosed with influenza. RESULTS: In non-pandemic [pandemic] years, delay distribution medians (inter-quartile ranges) were: Service to Admission 6.3 days (0.1-17.6 days) [2.4 days (-0.3-13.6 days)], Admission to Discharge 3 days (1.4-5.9 days) [2.6 days (1.2-5.1 days)], Admission to Death 5.3 days (2.1-11 days) [6 days (2.6-13.1 days)]. (Service date is defined as the date, within the year, of the first health care billing that included a diagnostic code for influenza-like-illness.) Among individuals diagnosed with either pneumonia or influenza in a given year, 19% [16%] were hospitalized within the year and 3% [2%] died in hospital. Among all individuals who were hospitalized, 10% [12%] were diagnosed with pneumonia or influenza during the year and 5% [5%] died in hospital. CONCLUSION: Our empirical delay distributions and conditional probabilities should help facilitate more accurate forecasts in the future, including improved predictions of hospital bed demands during influenza outbreaks, and the expected effects of hospitalizations on epidemic dynamics.

Stochastic simulation of multiscale complex systems with PISKaS: A rule-based approach.

Computational simulation is a widely employed methodology to study the dynamic behavior of complex systems. Although common approaches are based either on ordinary differential equations or stochastic differential equations, these techniques make several assumptions which, when it comes to biological processes, could often lead to unrealistic models. Among others, model approaches based on differential equations entangle kinetics and causality, failing when complexity increases, separating knowledge from models, and assuming that the average behavior of the population encompasses any individual deviation. To overcome these limitations, simulations based on the Stochastic Simulation Algorithm (SSA) appear as a suitable approach to model complex biological systems. In this work, we review three different models executed in PISKaS: a rule-based framework to produce multiscale stochastic simulations of complex systems. These models span multiple time and spatial scales ranging from gene regulation up to Game Theory. In the first example, we describe a model of the core regulatory network of gene expression in Escherichia coli highlighting the continuous model improvement capacities of PISKaS. The second example describes a hypothetical outbreak of the Ebola virus occurring in a compartmentalized environment resembling cities and highways. Finally, in the last example, we illustrate a stochastic model for the prisoner's dilemma; a common approach from social sciences describing complex interactions involving trust within human populations. As whole, these models demonstrate the capabilities of PISKaS providing fertile scenarios where to explore the dynamics of complex systems.

The Hayflick Limit May Determine the Effective Clonal Diversity of Naive T Cells.

Having a large number of sufficiently abundant T cell clones is important for adequate protection against diseases. However, as shown in this paper and elsewhere, between young adulthood and >70 y of age the effective clonal diversity of naive CD4/CD8 T cells found in human blood declines by a factor of >10. (Effective clonal diversity accounts for both the number and the abundance of T cell clones.) The causes of this observation are incompletely understood. A previous study proposed that it might result from the emergence of certain rare, replication-enhancing mutations in T cells. In this paper, we propose an even simpler explanation: that it results from the loss of T cells that have attained replicative senescence (i.e., the Hayflick limit). Stochastic numerical simulations of naive T cell population dynamics, based on experimental parameters, show that the rate of homeostatic T cell proliferation increases after the age of ∼60 y because naive T cells collectively approach replicative senescence. This leads to a sharp decline of effective clonal diversity after ∼70 y, in agreement with empirical data. A mathematical analysis predicts that, without an increase in the naive T cell proliferation rate, this decline will occur >50 yr later than empirically observed. These results are consistent with a model in which exhaustion of the proliferative capacity of naive T cells causes a sharp decline of their effective clonal diversity and imply that therapeutic potentiation of thymopoiesis might either prevent or reverse this outcome.