Stochastic social behavior coupled to COVID-19 dynamics leads to waves, plateaus, and an endemic state

Evidence for Behavioral Autorepression in Covid-19 Epidemiological Dynamics

It has long been hypothesized that behavioral reactions to epidemic severity autoregulate infection dynamics, for example when susceptible individuals self-sequester based on perceived levels of circulating disease. However, evidence for such 'behavioral autorepression' has remained elusive, and its presence could significantly affect epidemic forecasting and interventions. Here, we analyzed early COVID-19 dynamics at 708 locations over three epidemiological scales (96 countries, 50 US states, and 562 US counties). Signatures of behavioral autorepression were identified through: (i) a counterintuitive mobility-death correlation, (ii) fluctuation-magnitude analysis, and (iii) dynamics of SARS-CoV-2 infection waves. These data enabled calculation of the average behavioral-autorepression strength (i.e., negative feedback 'gain') across different populations. Surprisingly, incorporating behavioral autorepression into conventional models was required to accurately forecast COVID-19 mortality. Models also predicted that the strength of behavioral autorepression has the potential to alter the efficacy of non-pharmaceutical interventions. Overall, these results provide evidence for the long-hypothesized existence of behavioral autorepression, which could improve epidemic forecasting and enable more effective application of non-pharmaceutical interventions during future epidemics.

Significance

Challenges with epidemiological forecasting during the COVID-19 pandemic suggested gaps in underlying model architecture. One long-held hypothesis, typically omitted from conventional models due to lack of empirical evidence, is that human behaviors lead to intrinsic negative autoregulation of epidemics (termed 'behavioral autorepression'). This omission substantially alters model forecasts. Here, we provide independent lines of evidence for behavioral autorepression during the COVID-19 pandemic, demonstrate that it is sufficient to explain counterintuitive data on 'shutdowns', and provides a mechanistic explanation of why early shutdowns were more effective than delayed, high-intensity shutdowns. We empirically measure autorepression strength, and show that incorporating autorepression dramatically improves epidemiological forecasting. The autorepression phenomenon suggests that tailoring interventions to specific populations may be warranted.

Heterogeneous risk tolerance, in-groups, and epidemic waves

There is a growing interest in the joint modeling of the dynamics of disease and health-related beliefs and attitudes, but coupling mechanisms are yet to be understood. We introduce a model where risk information, which can be delayed, comes in two flavors, including historical risk derived from perceived incidence data and predicted risk information. Our model also includes an interpretation domain where the behavioral response to risk information is subject to in-group pressure. We then simulate how the strength of behavioral reaction impacts epidemic severity as measured by epidemic peak size, number of waves, and final size. Simulated behavioral response is not effective when the level of protection that prophylactic behavior provides is as small as 50% or lower. At a higher level of 75% or more, we see the emergence of multiple epidemic waves. In addition, simulations show that different behavioral response profiles can lead to various epidemic outcomes that are non-monotonic with the strength of reaction to risk information. We also modeled heterogeneity in the response profile of a population and find they can lead to less severe epidemic outcome in terms of peak size.

Estimating hidden relationships in dynamical systems: Discovering drivers of infection rates of COVID-19

Discovering causal influences among internal variables is a fundamental goal of complex systems research. This paper presents a framework for uncovering hidden relationships from limited time-series data by combining methods from nonlinear estimation and information theory. The approach is based on two sequential steps: first, we reconstruct a more complete state of the underlying dynamical system, and second, we calculate mutual information between pairs of internal state variables to detail causal dependencies. Equipped with time-series data related to the spread of COVID-19 from the past three years, we apply this approach to identify the drivers of falling and rising infections during the three main waves of infection in the Chicago metropolitan region. The unscented Kalman filter nonlinear estimation algorithm is implemented on an established epidemiological model of COVID-19, which we refine to include isolation, masking, loss of immunity, and stochastic transition rates. Through the systematic study of mutual information between infection rate and various stochastic parameters, we find that increased mobility, decreased mask use, and loss of immunity post sickness played a key role in rising infections, while falling infections were controlled by masking and isolation.

Stabilizing the return to normal behavior in an epidemic

Predicting the interplay between infectious disease and behavior has been an intractable problem because behavioral response is so varied. We introduce a general framework for feedback between incidence and behavior for an infectious disease. By identifying stable equilibria, we provide policy end-states that are self-managing and self-maintaining. We prove mathematically the existence of two new endemic equilibria depending on the vaccination rate: one in the presence of low vaccination but with reduced societal activity (the "new normal"), and one with return to normal activity but with vaccination rate below that required for disease elimination. This framework allows us to anticipate the long-term consequence of an emerging disease and design a vaccination response that optimizes public health and limits societal consequences.

Communicating COVID-19 exposure risk with an interactive website counteracts risk misestimation

During the COVID-19 pandemic, individuals depended on risk information to make decisions about everyday behaviors and public policy. Here, we assessed whether an interactive website influenced individuals' risk tolerance to support public health goals. We collected data from 11,169 unique users who engaged with the online COVID-19 Event Risk Tool (https://covid19risk.biosci.gatech.edu/) between 9/22/21 and 1/22/22. The website featured interactive elements, including a dynamic risk map, survey questions, and a risk quiz with accuracy feedback. After learning about the risk of COVID-19 exposure, participants reported being less willing to participate in events that could spread COVID-19, especially for high-risk large events. We also uncovered a bias in risk estimation: Participants tended to overestimate the risk of small events but underestimate the risk of large events. Importantly, even participants who voluntarily sought information about COVID risks tended to misestimate exposure risk, demonstrating the need for intervention. Participants from liberal-leaning counties were more likely to use the website tools and more responsive to feedback about risk misestimation, indicating that political partisanship influences how individuals seek and engage with COVID-19 information. Lastly, we explored temporal dynamics and found that user engagement and risk estimation fluctuated over the course of the Omicron variant outbreak. Overall, we report an effective large-scale method for communicating viral exposure risk; our findings are relevant to broader research on risk communication, epidemiological modeling, and risky decision-making.

Physical distancing versus testing with self-isolation for controlling an emerging epidemic

Two distinct strategies for controlling an emerging epidemic are physical distancing and regular testing with self-isolation. These strategies are especially important before effective vaccines or treatments become widely available. The testing strategy has been promoted frequently but used less often than physical distancing to mitigate COVID-19. We compared the performance of these strategies in an integrated epidemiological and economic model that includes a simple representation of transmission by "superspreading," wherein a relatively small fraction of infected individuals cause a large share of infections. We examined the economic benefits of distancing and testing over a wide range of conditions, including variations in the transmissibility and lethality of the disease meant to encompass the most prominent variants of COVID-19 encountered so far. In a head-to-head comparison using our primary parameter values, both with and without superspreading and a declining marginal value of mortality risk reductions, an optimized testing strategy outperformed an optimized distancing strategy. In a Monte Carlo uncertainty analysis, an optimized policy that combined the two strategies performed better than either one alone in more than 25% of random parameter draws. Insofar as diagnostic tests are sensitive to viral loads, and individuals with high viral loads are more likely to contribute to superspreading events, superspreading enhances the relative performance of testing over distancing in our model. Both strategies performed best at moderate levels of transmissibility, somewhat lower than the transmissibility of the ancestral strain of SARS-CoV-2.

A model of behavioural response to risk accurately predicts the statistical distribution of COVID-19 infection and reproduction numbers

One clear aspect of behaviour in the COVID-19 pandemic has been people's focus on, and response to, reported or observed infection numbers in their community. We describe a simple model of infectious disease spread in a pandemic situation where people's behaviour is influenced by the current risk of infection and where this behavioural response acts homeostatically to return infection risk to a certain preferred level. This homeostatic response is active until approximate herd immunity is reached: in this domain the model predicts that the reproduction rate R will be centred around a median of 1, that proportional change in infection numbers will follow the standard Cauchy distribution with location and scale parameters 0 and 1, and that high infection numbers will follow a power-law frequency distribution with exponent 2. To test these predictions we used worldwide COVID-19 data from 1st February 2020 to 30th June 2022 to calculate [Formula: see text] confidence interval estimates across countries for these R, location, scale and exponent parameters. The resulting median R estimate was [Formula: see text] (predicted value 1) the proportional change location estimate was [Formula: see text] (predicted value 0), the proportional change scale estimate was [Formula: see text] (predicted value 1), and the frequency distribution exponent estimate was [Formula: see text] (predicted value 2); in each case the observed estimate agreed with model predictions.

Waves of infection emerging from coupled social and epidemiological dynamics

The coronavirus (SARS-CoV-2) exhibited waves of infection in 2020 and 2021 in Japan. The number of infected had multiple distinct peaks at intervals of several months. One possible process causing these waves of infection is people switching their activities in response to the prevalence of infection. In this paper, we present a simple model for the coupling of social and epidemiological dynamics. The assumptions are as follows. Each person switches between active and restrained states. Active people move more often to crowded areas, interact with each other, and suffer a higher rate of infection than people in the restrained state. The rate of transition from restrained to active states is enhanced by the fraction of currently active people (conformity), whereas the rate of backward transition is enhanced by the abundance of infected people (risk avoidance). The model may show transient or sustained oscillations, initial-condition dependence, and various bifurcations. The infection is maintained at a low level if the recovery rate is between the maximum and minimum levels of the force of infection. In addition, waves of infection may emerge instead of converging to the stationary abundance of infected people if both conformity and risk avoidance of people are strong.

A generalized distributed delay model of COVID-19: An endemic model with immunity waning

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been spreading worldwide for over two years, with millions of reported cases and deaths. The deployment of mathematical modeling in the fight against COVID-19 has recorded tremendous success. However, most of these models target the epidemic phase of the disease. The development of safe and effective vaccines against SARS-CoV-2 brought hope of safe reopening of schools and businesses and return to pre-COVID normalcy, until mutant strains like the Delta and Omicron variants, which are more infectious, emerged. A few months into the pandemic, reports of the possibility of both vaccine- and infection-induced immunity waning emerged, thereby indicating that COVID-19 may be with us for longer than earlier thought. As a result, to better understand the dynamics of COVID-19, it is essential to study the disease with an endemic model. In this regard, we developed and analyzed an endemic model of COVID-19 that incorporates the waning of both vaccine- and infection-induced immunities using distributed delay equations. Our modeling framework assumes that the waning of both immunities occurs gradually over time at the population level. We derived a nonlinear ODE system from the distributed delay model and showed that the model could exhibit either a forward or backward bifurcation depending on the immunity waning rates. Having a backward bifurcation implies that $ R_c < 1 $ is not sufficient to guarantee disease eradication, and that the immunity waning rates are critical factors in eradicating COVID-19. Our numerical simulations show that vaccinating a high percentage of the population with a safe and moderately effective vaccine could help in eradicating COVID-19.

A systematic procedure for incorporating separable static heterogeneity into compartmental epidemic models

In this paper, we show how to modify a compartmental epidemic model, without changing the dimension, such that separable static heterogeneity is taken into account. The derivation is based on the Kermack-McKendrick renewal equation.

Host age structure reshapes parasite symbiosis: collaboration begets pathogens, competition begets virulent mutualists

Background

Symbiotic relationships are ubiquitous in the biosphere. Inter-species symbiosis is impacted by intra-specific distinctions, in particular, those defined by the age structure of a population. Older individuals compete with younger individuals for resources despite being less likely to reproduce, diminishing the fitness of the population. Conversely, however, older individuals can support the reproduction of younger individuals, increasing the population fitness. Parasitic relationships are commonly age structured, typically, more adversely affecting older hosts.

Results

We employ mathematical modeling to explore the differential effects of collaborative or competitive host age structures on host-parasite relationships. A classical epidemiological compartment model is constructed with three disease states: susceptible, infected, and recovered. Each of these three states is partitioned into two compartments representing young, potentially reproductive, and old, post-reproductive, hosts, yielding 6 compartments in total. In order to describe competition and collaboration between old and young compartments, we model the reproductive success to depend on the fraction of young individuals in the population. Collaborative populations with relatively greater numbers of post-reproductive hosts enjoy greater reproductive success whereas in purely competitive populations, increasing the post-reproductive subpopulation reduces reproductive success.

Conclusions

We demonstrate that, in collaborative host populations, pathogens strictly impacting older, post-reproductive individuals can reduce population fitness even more than pathogens that directly impact younger, potentially reproductive individuals. In purely competitive populations, the reverse is observed, and we demonstrate that endemic, virulent pathogens can oxymoronically form a mutualistic relationship with the host, increasing the fitness of the host population. Applications to endangered species conservation and invasive species containment are discussed.

The effect of competition between health opinions on epidemic dynamics

Past major epidemic events showed that when an infectious disease is perceived to cause severe health outcomes, individuals modify health behavior affecting epidemic dynamics. To investigate the effect of this feedback relationship on epidemic dynamics, we developed a compartmental model that couples a disease spread framework with competition of two mutually exclusive health opinions (health-positive and health-neutral) associated with different health behaviors. The model is based on the assumption that individuals switch health opinions as a result of exposure to opinions of others through interpersonal communications. To model opinion switch rates, we considered a family of functions and identified the ones that allow health opinions to coexist. Finally, the model includes assortative mixing by opinions. In the disease-free population, either the opinions cannot coexist and one of them is always dominating (mono-opinion equilibrium) or there is at least one stable coexistence of opinions equilibrium. In the latter case, there is multistability between the coexistence equilibrium and the two mono-opinion equilibria. When two opinions coexist, it depends on their distribution whether the infection can invade. If presence of the infection leads to increased switching to a health-positive opinion, the epidemic burden becomes smaller than indicated by the basic reproduction number. Additionally, a feedback between epidemic dynamics and health opinion dynamics may result in (sustained) oscillatory dynamics and a switch to a different stable opinion distribution. Our model captures feedback between spread of awareness through social interactions and infection dynamics and can serve as a basis for more elaborate individual-based models.

Sars-Cov2 world pandemic recurrent waves controlled by variants evolution and vaccination campaign

While understanding the time evolution of Covid-19 pandemic is needed to plan economics and tune sanitary policies, a quantitative information of the recurrent epidemic waves is elusive. This work describes a statistical physics study of the subsequent waves in the epidemic spreading of Covid-19 and disclose the frequency components of the epidemic waves pattern over two years in United States, United Kingdom and Japan. These countries have been taken as representative cases of different containment policies such as "Mitigation" (USA and UK) and "Zero Covid" (Japan) policies. The supercritical phases in spreading have been identified by intervals with RIC-index > 0. We have used the wavelet transform of infection and fatality waves to get the spectral analysis showing a dominant component around 130 days. Data of the world dynamic clearly indicates also the crossover to a different phase due to the enforcement of vaccination campaign. In Japan and United Kingdom, we observed the emergence in the infection waves of a long period component (~ 170 days) during vaccination campaign. These results indicate slowing down of the epidemic spreading dynamics due to the vaccination campaign. Finally, we find an intrinsic difference between infection and fatality waves pointing to a non-trivial variation of the lethality due to different gene variants.

Controlling Multiple COVID-19 Epidemic Waves: An Insight from a Multi-scale Model Linking the Behaviour Change Dynamics to the Disease Transmission Dynamics

COVID-19 epidemics exhibited multiple waves regionally and globally since 2020. It is important to understand the insight and underlying mechanisms of the multiple waves of COVID-19 epidemics in order to design more efficient non-pharmaceutical interventions (NPIs) and vaccination strategies to prevent future waves. We propose a multi-scale model by linking the behaviour change dynamics to the disease transmission dynamics to investigate the effect of behaviour dynamics on COVID-19 epidemics using game theory. The proposed multi-scale models are calibrated and key parameters related to disease transmission dynamics and behavioural dynamics with/without vaccination are estimated based on COVID-19 epidemic data (daily reported cases and cumulative deaths) and vaccination data. Our modeling results demonstrate that the feedback loop between behaviour changes and COVID-19 transmission dynamics plays an essential role in inducing multiple epidemic waves. We find that the long period of high-prevalence or persistent deterioration of COVID-19 epidemics could drive almost all of the population to change their behaviours and maintain the altered behaviours. However, the effect of behaviour changes fades out gradually along the progress of epidemics. This suggests that it is essential to have not only persistent, but also effective behaviour changes in order to avoid subsequent epidemic waves. In addition, our model also suggests the importance to maintain the effective altered behaviours during the initial stage of vaccination, and to counteract relaxation of NPIs, it requires quick and massive vaccination to avoid future epidemic waves.

Comparison of the 2021 COVID-19 roadmap projections against public health data in England

Control and mitigation of the COVID-19 pandemic in England has relied on a combination of vaccination and non-pharmaceutical interventions (NPIs). Some of these NPIs are extremely costly (economically and socially), so it was important to relax these promptly without overwhelming already burdened health services. The eventual policy was a Roadmap of four relaxation steps throughout 2021, taking England from lock-down to the cessation of all restrictions on social interaction. In a series of six Roadmap documents generated throughout 2021, models assessed the potential risk of each relaxation step. Here we show that the model projections generated a reliable estimation of medium-term hospital admission trends, with the data points up to September 2021 generally lying within our 95% prediction intervals. The greatest uncertainties in the modelled scenarios came from vaccine efficacy estimates against novel variants, and from assumptions about human behaviour in the face of changing restrictions and risk.

Continuing inequalities in COVID-19 mortality in England and Wales, and the changing importance of regional, over local, deprivation

BACKGROUND

Observational studies have highlighted that where individuals live is far more important for risk of dying with COVID-19, than for dying of other causes. Deprivation is commonly proposed as explaining such differences. During the period of localised restrictions in late 2020, areas with higher restrictions tended to be more deprived. We explore how this impacted the relationship between deprivation and mortality and see whether local or regional deprivation matters more for inequalities in COVID-19 mortality.

METHODS

We use publicly available population data on deaths due to COVID-19 and all-cause mortality between March 2020 and April 2021 to investigate the scale of spatial inequalities. We use a multiscale approach to simultaneously consider three spatial scales through which processes driving inequalities may act. We go on to explore whether deprivation explains such inequalities.

RESULTS

Adjusting for population age structure and number of care homes, we find highest regional inequality in October 2020, with a COVID-19 mortality rate ratio of 5.86 (95% CI 3.31 to 19.00) for the median between-region comparison. We find spatial context is most important, and spatial inequalities higher, during periods of low mortality. Almost all unexplained spatial inequality in October 2020 is removed by adjusting for deprivation. During October 2020, one standard deviation increase in regional deprivation was associated with 20% higher local mortality (95% CI, 1.10 to 1.30).

CONCLUSIONS

Spatial inequalities are greatest in periods of lowest overall mortality, implying that as mortality declines it does not do so equally. During the prolonged period of low restrictions and low mortality in summer 2020, spatial inequalities strongly increased. Contrary to previous months, we show that the strong spatial patterning during autumn 2020 is almost entirely explained by deprivation. As overall mortality declines, policymakers must be proactive in detecting areas where this is not happening, or risk worsening already strong health inequalities.

Mitigation of SARS-CoV-2 transmission at a large public university

In Fall 2020, universities saw extensive transmission of SARS-CoV-2 among their populations, threatening health of the university and surrounding communities, and viability of in-person instruction. Here we report a case study at the University of Illinois at Urbana-Champaign, where a multimodal “SHIELD: Target, Test, and Tell” program, with other non-pharmaceutical interventions, was employed to keep classrooms and laboratories open. The program included epidemiological modeling and surveillance, fast/frequent testing using a novel low-cost and scalable saliva-based RT-qPCR assay for SARS-CoV-2 that bypasses RNA extraction, called covidSHIELD, and digital tools for communication and compliance. In Fall 2020, we performed >1,000,000 covidSHIELD tests, positivity rates remained low, we had zero COVID-19-related hospitalizations or deaths amongst our university community, and mortality in the surrounding Champaign County was reduced more than 4-fold relative to expected. This case study shows that fast/frequent testing and other interventions mitigated transmission of SARS-CoV-2 at a large public university.

Mitigation of SARS-CoV-2 transmission at a large public university

In Fall 2020, universities saw extensive transmission of SARS-CoV-2 among their populations, threatening health of the university and surrounding communities, and viability of in-person instruction. Here we report a case study at the University of Illinois at Urbana-Champaign, where a multimodal “SHIELD: Target, Test, and Tell” program, with other non-pharmaceutical interventions, was employed to keep classrooms and laboratories open. The program included epidemiological modeling and surveillance, fast/frequent testing using a novel low-cost and scalable saliva-based RT-qPCR assay for SARS-CoV-2 that bypasses RNA extraction, called covidSHIELD, and digital tools for communication and compliance. In Fall 2020, we performed >1,000,000 covidSHIELD tests, positivity rates remained low, we had zero COVID-19-related hospitalizations or deaths amongst our university community, and mortality in the surrounding Champaign County was reduced more than 4-fold relative to expected. This case study shows that fast/frequent testing and other interventions mitigated transmission of SARS-CoV-2 at a large public university.

Safely opening university campuses has been a major challenge during the COVID-19 pandemic. Here, the authors describe a program of public health measures employed at a university in the United States which, combined with other non-pharmaceutical interventions, allowed the university to stay open in fall 2020 with limited evidence of transmission.