Data-driven contact structures: From homogeneous mixing to multilayer networks

Generalized contact matrices allow integrating socioeconomic variables into epidemic models

Variables related to socioeconomic status (SES), including income, ethnicity, and education, shape contact structures and affect the spread of infectious diseases. However, these factors are often overlooked in epidemic models, which typically stratify social contacts by age and interaction contexts. Here, we introduce and study generalized contact matrices that stratify contacts across multiple dimensions. We demonstrate a lower-bound theorem proving that disregarding additional dimensions, besides age and context, might lead to an underestimation of the basic reproductive number. By using SES variables in both synthetic and empirical data, we illustrate how generalized contact matrices enhance epidemic models, capturing variations in behaviors such as heterogeneous levels of adherence to nonpharmaceutical interventions among demographic groups. Moreover, we highlight the importance of integrating SES traits into epidemic models, as neglecting them might lead to substantial misrepresentation of epidemic outcomes and dynamics. Our research contributes to the efforts aiming at incorporating socioeconomic and other dimensions into epidemic modeling.

Reconciling heterogeneous dengue virus infection risk estimates from different study designs

Uncovering rates at which susceptible individuals become infected with a pathogen, i.e. the force of infection (FOI), is essential for assessing transmission risk and reconstructing distribution of immunity in a population. For dengue, reconstructing exposure and susceptibility statuses from the measured FOI is of particular significance as prior exposure is a strong risk factor for severe disease. FOI can be measured via many study designs. Longitudinal serology are considered gold standard measurements, as they directly track the transition of seronegative individuals to seropositive due to incident infections (seroincidence). Cross-sectional serology can provide estimates of FOI by contrasting seroprevalence across ages. Age of reported cases can also be used to infer FOI. Agreement of these measurements, however, have not been assessed. Using 26 years of data from cohort studies and hospital-attended cases from Kamphaeng Phet province, Thailand, we found FOI estimates from the three sources to be highly inconsistent. Annual FOI estimates from seroincidence was 2.46 to 4.33-times higher than case-derived FOI. Correlation between seroprevalence-derived and case-derived FOI was moderate (correlation coefficient=0.46) and no systematic bias. Through extensive simulations and theoretical analysis, we show that incongruences between methods can result from failing to account for dengue antibody kinetics, assay noise, and heterogeneity in FOI across ages. Extending standard inference models to include these processes reconciled the FOI and susceptibility estimates. Our results highlight the importance of comparing inferences across multiple data types to uncover additional insights not attainable through a single data type/analysis.

Building a Simple Model to Assess the Impact of Case Investigation and Contact Tracing for Sexually Transmitted Diseases: Lessons From COVID-19

Introduction

During the COVID-19 pandemic, the U.S. Centers for Disease Control and Prevention developed a simple spreadsheet-based tool to help state and local public health officials assess the performance and impact of COVID-19 case investigation and contact tracing in their jurisdiction. The applicability and feasibility of building such a tool for sexually transmitted diseases were assessed.

Methods

The key epidemiologic differences between sexually transmitted diseases and respiratory diseases (e.g., mixing patterns, incubation period, duration of infection, and the availability of treatment) were identified, and their implications for modeling case investigation and contact tracing impact with a simple spreadsheet tool were remarked on. Existing features of the COVID-19 tool that are applicable for evaluating the impact of case investigation and contact tracing for sexually transmitted diseases were also identified.

Results

Our findings offer recommendations for the future development of a spreadsheet-based modeling tool for evaluating the impact of sexually transmitted disease case investigation and contact tracing efforts. Generally, we advocate for simplifying sexually transmitted disease-specific complexities and performing sensitivity analyses to assess uncertainty. The authors also acknowledge that more complex modeling approaches might be required but note that it is possible that a sexually transmitted disease case investigation and contact tracing tool could incorporate features from more complex models while maintaining a user-friendly interface.

Conclusions

A sexually transmitted disease case investigation and contact tracing tool could benefit from the incorporation of key features of the COVID-19 model, namely its user-friendly interface. The inherent differences between sexually transmitted diseases and respiratory viruses should not be seen as a limitation to the development of such tool.

Sociodemographic Variables Can Guide Prioritized Testing Strategies for Epidemic Control in Resource-Limited Contexts

BACKGROUND

Targeted surveillance allows public health authorities to implement testing and isolation strategies when diagnostic resources are limited, and can be implemented via the consideration of social network topologies. However, it remains unclear how to implement such surveillance and control when network data are unavailable.

METHODS

We evaluated the ability of sociodemographic proxies of degree centrality to guide prioritized testing of infected individuals compared to known degree centrality. Proxies were estimated via readily available sociodemographic variables (age, gender, marital status, educational attainment, household size). We simulated severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) epidemics via a susceptible-exposed-infected-recovered individual-based model on 2 contact networks from rural Madagascar to test applicability of these findings to low-resource contexts.

RESULTS

Targeted testing using sociodemographic proxies performed similarly to targeted testing using known degree centralities. At low testing capacity, using proxies reduced infection burden by 22%-33% while using 20% fewer tests, compared to random testing. By comparison, using known degree centrality reduced the infection burden by 31%-44% while using 26%-29% fewer tests.

CONCLUSIONS

We demonstrate that incorporating social network information into epidemic control strategies is an effective countermeasure to low testing capacity and can be implemented via sociodemographic proxies when social network data are unavailable.

Evaluating the impact of test-trace-isolate for COVID-19 management and alternative strategies

There are many contrasting results concerning the effectiveness of Test-Trace-Isolate (TTI) strategies in mitigating SARS-CoV-2 spread. To shed light on this debate, we developed a novel static-temporal multiplex network characterizing both the regular (static) and random (temporal) contact patterns of individuals and a SARS-CoV-2 transmission model calibrated with historical COVID-19 epidemiological data. We estimated that the TTI strategy alone could not control the disease spread: assuming R0 = 2.5, the infection attack rate would be reduced by 24.5%. Increased test capacity and improved contact trace efficiency only slightly improved the effectiveness of the TTI. We thus investigated the effectiveness of the TTI strategy when coupled with reactive social distancing policies. Limiting contacts on the temporal contact layer would be insufficient to control an epidemic and contacts on both layers would need to be limited simultaneously. For example, the infection attack rate would be reduced by 68.1% when the reactive distancing policy disconnects 30% and 50% of contacts on static and temporal layers, respectively. Our findings highlight that, to reduce the overall transmission, it is important to limit contacts regardless of their types in addition to identifying infected individuals through contact tracing, given the substantial proportion of asymptomatic and pre-symptomatic SARS-CoV-2 transmission.

Impact of contact data resolution on the evaluation of interventions in mathematical models of infectious diseases

Computational models offer a unique setting to test strategies to mitigate the spread of infectious diseases, providing useful insights to applied public health. To be actionable, models need to be informed by data, which can be available at different levels of detail. While high-resolution data describing contacts between individuals are increasingly available, data gathering remains challenging, especially during a health emergency. Many models thus use synthetic data or coarse information to evaluate intervention protocols. Here, we evaluate how the representation of contact data might affect the impact of various strategies in models, in the realm of COVID-19 transmission in educational and work contexts. Starting from high-resolution contact data, we use detailed to coarse data representations to inform a model of SARS-CoV-2 transmission and simulate different mitigation strategies. We find that coarse data representations estimate a lower risk of superspreading events. However, the rankings of protocols according to their efficiency or cost remain coherent across representations, ensuring the consistency of model findings to inform public health advice. Caution should be taken, however, on the quantitative estimations of those benefits and costs triggering the adoption of protocols, as these may depend on data representation.

Quantifying human mixing patterns in Chinese provinces outside Hubei after the 2020 lockdown was lifted

Background

Contact patterns play a key role in the spread of respiratory infectious diseases in human populations. During the COVID-19 pandemic, the regular contact patterns of the population have been disrupted due to social distancing both imposed by the authorities and individual choices. Many studies have focused on age-mixing patterns before the COVID-19 pandemic, but they provide very little information about the mixing patterns in the COVID-19 era. In this study, we aim at quantifying human heterogeneous mixing patterns immediately after lockdowns implemented to contain COVID-19 spread in China were lifted. We also provide an illustrative example of how the collected mixing patterns can be used in a simulation study of SARS-CoV-2 transmission.

Methods and results

In this work, a contact survey was conducted in Chinese provinces outside Hubei in March 2020, right after lockdowns were lifted. We then leveraged the estimated mixing patterns to calibrate a mathematical model of SARS-CoV-2 transmission. Study participants reported 2.3 contacts per day (IQR: 1.0–3.0) and the mean per-contact duration was 7.0 h (IQR: 1.0–10.0). No significant differences in average contact number and contact duration were observed between provinces, the number of recorded contacts did not show a clear trend by age, and most of the recorded contacts occurred with family members (about 78%). The simulation study highlights the importance of considering age-specific contact patterns to estimate the COVID-19 burden.

Conclusions

Our findings suggest that, despite lockdowns were no longer in place at the time of the survey, people were still heavily limiting their contacts as compared to the pre-pandemic situation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12879-022-07455-7.

Assessing the role of multiple mechanisms increasing the age of dengue cases in Thailand

The age of reported dengue hemorrhagic fever (DHF) cases, the severe form of dengue infections, has been increasing in Thailand for four decades. Factors underlying this shift remain poorly understood, challenging public health planning. Here, we found aging of the population and its effect on the hazard of transmission to be the dominant contributors, with temporal changes in surveillance practices playing a lesser role. With ongoing population aging, we expect a continuing shift of DHF toward older individuals, heightening the chance of clinical complications with comorbidities. With most other highly endemic countries facing similar shifts in age structure, the pattern is expected to appear elsewhere. Awareness is needed to improve diagnosis and treatment.

The mean age of dengue hemorrhagic fever (DHF) cases increased considerably in Thailand from 8.1 to 24.3 y between 1981 and 2017 (mean annual increase of 0.45 y). Alternative proposed explanations for this trend, such as changes in surveillance practices, reduced mosquito–human contact, and shifts in population demographics, have different implications for global dengue epidemiology. To evaluate the contribution of each of these hypothesized mechanisms to the observed data, we developed 20 nested epidemiological models of dengue virus infection, allowing for variation over time in population demographics, infection hazards, and reporting rates. We also quantified the effect of removing or retaining each source of variation in simulations of the age trajectory. Shifts in the age structure of susceptibility explained 58% of the observed change in age. Adding heterogeneous reporting by age and reductions in per-serotype infection hazard to models with shifts in susceptibility explained an additional 42%. Reductions in infection hazards were mostly driven by changes in the number of infectious individuals at any time (another consequence of shifting age demographics) rather than changes in the transmissibility of individual infections. We conclude that the demographic transition drives the overwhelming majority of the observed change as it changes both the age structure of susceptibility and the number of infectious individuals. With the projected Thai population age structure, our results suggest a continuing increase in age of DHF cases, shifting the burden toward individuals with more comorbidity. These insights into dengue epidemiology may be relevant to many regions of the globe currently undergoing comparable changes in population demographics.

Why Plasmodium vivax and Plasmodium falciparum are so different? A tale of two clades and their species diversities

The global malaria burden sometimes obscures that the genus Plasmodium comprises diverse clades with lineages that independently gave origin to the extant human parasites. Indeed, the differences between the human malaria parasites were highlighted in the classical taxonomy by dividing them into two subgenera, the subgenus Plasmodium, which included all the human parasites but Plasmodium falciparum that was placed in its separate subgenus, Laverania. Here, the evolution of Plasmodium in primates will be discussed in terms of their species diversity and some of their distinct phenotypes, putative molecular adaptations, and host–parasite biocenosis. Thus, in addition to a current phylogeny using genome-level data, some specific molecular features will be discussed as examples of how these parasites have diverged. The two subgenera of malaria parasites found in primates, Plasmodium and Laverania, reflect extant monophyletic groups that originated in Africa. However, the subgenus Plasmodium involves species in Southeast Asia that were likely the result of adaptive radiation. Such events led to the Plasmodium vivax lineage. Although the Laverania species, including P. falciparum, has been considered to share “avian characteristics,” molecular traits that were likely in the common ancestor of primate and avian parasites are sometimes kept in the Plasmodium subgenus while being lost in Laverania. Assessing how molecular traits in the primate malaria clades originated is a fundamental science problem that will likely provide new targets for interventions. However, given that the genus Plasmodium is paraphyletic (some descendant groups are in other genera), understanding the evolution of malaria parasites will benefit from studying “non-Plasmodium” Haemosporida.

Directed percolation in random temporal network models with heterogeneities

The event graph representation of temporal networks suggests that the connectivity of temporal structures can be mapped to a directed percolation problem. However, similarly to percolation theory on static networks, this mapping is valid under the approximation that the structure and interaction dynamics of the temporal network are determined by its local properties, and, otherwise, it is maximally random. We challenge these conditions and demonstrate the robustness of this mapping in case of more complicated systems. We systematically analyze random and regular network topologies and heterogeneous link-activation processes driven by bursty renewal or self-exciting processes using numerical simulation and finite-size scaling methods. We find that the critical percolation exponents characterizing the temporal network are not sensitive to many structural and dynamical network heterogeneities, while they recover known scaling exponents characterizing directed percolation on low-dimensional lattices. While it is not possible to demonstrate the validity of this mapping for all temporal network models, our results establish the first batch of evidence supporting the robustness of the scaling relationships in the limited-time reachability of temporal networks.

Rapid Review of Social Contact Patterns During the COVID-19 Pandemic

BACKGROUND

Physical distancing measures aim to reduce person-to-person contact, a key driver of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) transmission. In response to unprecedented restrictions on human contact during the coronavirus disease 2019 (COVID-19) pandemic, studies measured social contact patterns under the implementation of physical distancing measures. This rapid review synthesizes empirical data on the changing social contact patterns during the COVID-19 pandemic.

METHOD

We conducted a systematic review using PubMed, Medline, Embase, and Google Scholar following the Preferred Reporting Items for Systematic Reviews and Meta-Analyses guidelines. We descriptively compared the distribution of contacts observed during the pandemic to pre-COVID data across countries to explore changes in contact patterns during physical distancing measures.

RESULTS

We identified 12 studies reporting social contact patterns during the COVID-19 pandemic. Eight studies were conducted in European countries and eleven collected data during the initial mitigation period in the spring of 2020 marked by government-declared lockdowns. Some studies collected additional data after relaxation of initial mitigation. Most study settings reported a mean of between 2 and 5 contacts per person per day, a substantial reduction compared to pre-COVID rates, which ranged from 7 to 26 contacts per day. This reduction was pronounced for contacts outside of the home. Consequently, levels of assortative mixing by age substantially declined. After relaxation of initial mitigation, mean contact rates increased but did not return to pre-COVID levels. Increases in contacts post-relaxation were driven by working-age adults.

CONCLUSION

Information on changes in contact patterns during physical distancing measures can guide more realistic representations of contact patterns in mathematical models for SARS-CoV-2 transmission.

Predicting the effects of COVID-19 related interventions in urban settings by combining activity-based modelling, agent-based simulation, and mobile phone data

Epidemiological simulations as a method are used to better understand and predict the spreading of infectious diseases, for example of COVID-19. This paper presents an approach that combines a well-established approach from transportation modelling that uses person-centric data-driven human mobility modelling with a mechanistic infection model and a person-centric disease progression model. The model includes the consequences of different room sizes, air exchange rates, disease import, changed activity participation rates over time (coming from mobility data), masks, indoors vs. outdoors leisure activities, and of contact tracing. It is validated against the infection dynamics in Berlin (Germany). The model can be used to understand the contributions of different activity types to the infection dynamics over time. It predicts the effects of contact reductions, school closures/vacations, masks, or the effect of moving leisure activities from outdoors to indoors in fall, and is thus able to quantitatively predict the consequences of interventions. It is shown that these effects are best given as additive changes of the reproduction number R. The model also explains why contact reductions have decreasing marginal returns, i.e. the first 50% of contact reductions have considerably more effect than the second 50%. Our work shows that is is possible to build detailed epidemiological simulations from microscopic mobility models relatively quickly. They can be used to investigate mechanical aspects of the dynamics, such as the transmission from political decisions via human behavior to infections, consequences of different lockdown measures, or consequences of wearing masks in certain situations. The results can be used to inform political decisions.

Heterogeneity matters: Contact structure and individual variation shape epidemic dynamics

In the recent COVID-19 pandemic, mathematical modeling constitutes an important tool to evaluate the prospective effectiveness of non-pharmaceutical interventions (NPIs) and to guide policy-making. Most research is, however, centered around characterizing the epidemic based on point estimates like the average infectiousness or the average number of contacts. In this work, we use stochastic simulations to investigate the consequences of a population's heterogeneity regarding connectivity and individual viral load levels. Therefore, we translate a COVID-19 ODE model to a stochastic multi-agent system. We use contact networks to model complex interaction structures and a probabilistic infection rate to model individual viral load variation. We observe a large dependency of the dispersion and dynamical evolution on the population's heterogeneity that is not adequately captured by point estimates, for instance, used in ODE models. In particular, models that assume the same clinical and transmission parameters may lead to different conclusions, depending on different types of heterogeneity in the population. For instance, the existence of hubs in the contact network leads to an initial increase of dispersion and the effective reproduction number, but to a lower herd immunity threshold (HIT) compared to homogeneous populations or a population where the heterogeneity stems solely from individual infectivity variations.

Measles dynamics on network models with optimal control strategies

To investigate the influences of heterogeneity and waning immunity on measles transmission, we formulate a network model with periodic transmission rate, and theoretically examine the threshold dynamics. We numerically find that the waning of immunity can lead to an increase in the basic reproduction number R 0 and the density of infected individuals. Moreover, there exists a critical level for average degree above which R 0 increases quicker in the scale-free network than in the random network. To design the effective control strategies for the subpopulations with different activities, we examine the optimal control problem of the heterogeneous model. Numerical studies suggest us no matter what the network is, we should implement control measures as soon as possible once the outbreak takes off, and particularly, the subpopulation with high connectivity should require high intensity of interventions. However, with delayed initiation of controls, relatively strong control measures should be given to groups with medium degrees. Furthermore, the allocation of costs (or resources) should coincide with their contact patterns.

Modeling partial lockdowns in multiplex networks using partition strategies

National stay-at-home orders, or lockdowns, were imposed in several countries to drastically reduce the social interactions mainly responsible for the transmission of the SARS-CoV-2 virus. Despite being essential to slow down the COVID-19 pandemic, these containment measures are associated with an economic burden. In this work, we propose a network approach to model the implementation of a partial lockdown, breaking the society into disconnected components, or partitions. Our model is composed by two main ingredients: a multiplex network representing human contacts within different contexts, formed by a Household layer, a Work layer, and a Social layer including generic social interactions, and a Susceptible-Infected-Recovered process that mimics the epidemic spreading. We compare different partition strategies, with a twofold aim: reducing the epidemic outbreak and minimizing the economic cost associated to the partial lockdown. We also show that the inclusion of unconstrained social interactions dramatically increases the epidemic spreading, while different kinds of restrictions on social interactions help in keeping the benefices of the network partition.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s41109-021-00366-7.

Oscillatory dynamics in the dilemma of social distancing

Social distancing as one of the main non-pharmaceutical interventions can help slow down the spread of diseases, like in the COVID-19 pandemic. Effective social distancing, unless enforced as drastic lockdowns and mandatory cordon sanitaire, requires consistent strict collective adherence. However, it remains unknown what the determinants for the resultant compliance of social distancing and their impact on disease mitigation are. Here, we incorporate into the epidemiological process with an evolutionary game theory model that governs the evolution of social distancing behaviour. In our model, we assume an individual acts in their best interest and their decisions are driven by adaptive social learning of the real-time risk of infection in comparison with the cost of social distancing. We find interesting oscillatory dynamics of social distancing accompanied with waves of infection. Moreover, the oscillatory dynamics are dampened with a non-trivial dependence on model parameters governing decision-makings and gradually cease when the cumulative infections exceed the herd immunity. Compared to the scenario without social distancing, we quantify the degree to which social distancing mitigates the epidemic and its dependence on individuals’ responsiveness and rationality in their behaviour changes. Our work offers new insights into leveraging human behaviour in support of pandemic response.