Projecting social contact matrices in 152 countries using contact surveys and demographic data

Learning from the COVID-19 pandemic: A systematic review of mathematical vaccine prioritization models

As the world becomes ever more connected, the chance of pandemics increases as well. The recent COVID-19 pandemic and the concurrent global mass vaccine roll-out provides an ideal setting to learn from and refine our understanding of infectious disease models for better future preparedness. In this review, we systematically analyze and categorize mathematical models that have been developed to design optimal vaccine prioritization strategies of an initially limited vaccine. As older individuals are disproportionately affected by COVID-19, the focus is on models that take age explicitly into account. The lower mobility and activity level of older individuals gives rise to non-trivial trade-offs. Secondary research questions concern the optimal time interval between vaccine doses and spatial vaccine distribution. This review showcases the effect of various modeling assumptions on model outcomes. A solid understanding of these relationships yields better infectious disease models and thus public health decisions during the next pandemic.

Mathematical assessment of the roles of age heterogeneity and vaccination on the dynamics and control of SARS-CoV-2

The COVID-19 pandemic, caused by SARS-CoV-2, disproportionately affected certain segments of society, particularly the elderly population (which suffered the brunt of the burden of the pandemic in terms of severity of the disease, hospitalization, and death). This study presents a generalized multigroup model, with m heterogeneous sub-populations, to assess the population-level impact of age heterogeneity and vaccination on the transmission dynamics and control of the SARS-CoV-2 pandemic in the United States. Rigorous analysis of the model for the homogeneous case (i.e., the model with m = 1) reveal that its disease-free equilibrium is globally-asymptotically stable for two special cases (with perfect vaccine efficacy or negligible disease-induced mortality) whenever the associated reproduction number is less than one. The model has a unique and globally-asymptotically stable endemic equilibrium, for special a case, when the associated reproduction threshold exceeds one. The homogeneous model was fitted using the observed cumulative mortality data for the United States during three distinct waves (Waves A (October 17, 2020 to April 5, 2021), B (July 9, 2021 to November 7, 2021) and C (January 1, 2022 to May 7, 2022)) chosen to align with time periods when the Alpha, Delta and Omicron were, respectively, the predominant variants in the United States. The calibrated model was used to derive a theoretical expression for achieving vaccine-derived herd immunity (needed to eliminate the disease in the United States). It was shown that, using the one-group homogeneous model, vaccine-derived herd immunity is not attainable during Wave C of the pandemic in the United States, regardless of the coverage level of the fully-vaccinated individuals. Global sensitivity analysis was carried out to determine the parameters of the model that have the most influence on the disease dynamics and burden. These analyses reveal that control and mitigation strategies that may be very effective during one wave may not be so very effective during the other wave or waves. However, strategies that target asymptomatic and pre-symptomatic infectious individuals are shown to be consistently effective across all waves. To study the impact of the disproportionate effect of COVID-19 on the elderly population, we considered the heterogeneous model for the case where the total population is subdivided into the sub-populations of individuals under 65 years of age and those that are 65 and older. The resulting two-group heterogeneous model, which was also fitted using the cumulative mortality data for wave C, was also rigorously analysed. Unlike for the case of the one-group model, it was shown, for the two-group model, that vaccine-derived herd immunity can indeed be achieved during Wave C of the pandemic if at least 61% of the populace is fully vaccinated. Thus, this study shows that adding age heterogeneity into a SARS-CoV-2 vaccination model with homogeneous mixing significantly reduces the level of vaccination coverage needed to achieve vaccine-derived herd immunity (specifically, for the heterogeneous model, herd-immunity can be attained during Wave C if a moderate proportion of susceptible individuals are fully vaccinated). The consequence of this result is that vaccination models for SARS-CoV-2 that do not explicitly account for age heterogeneity may be overestimating the level of vaccine-derived herd immunity threshold needed to eliminate the SARS-CoV-2 pandemic.

High-resolution epidemiological landscape from ~290,000 SARS-CoV-2 genomes from Denmark

Vast amounts of pathogen genomic, demographic and spatial data are transforming our understanding of SARS-CoV-2 emergence and spread. We examined the drivers of molecular evolution and spread of 291,791 SARS-CoV-2 genomes from Denmark in 2021. With a sequencing rate consistently exceeding 60%, and up to 80% of PCR-positive samples between March and November, the viral genome set is broadly whole-epidemic representative. We identify a consistent rise in viral diversity over time, with notable spikes upon the importation of novel variants (e.g., Delta and Omicron). By linking genomic data with rich individual-level demographic data from national registers, we find that individuals aged  < 15 and  > 75 years had a lower contribution to molecular change (i.e., branch lengths) compared to other age groups, but similar molecular evolutionary rates, suggesting a lower likelihood of introducing novel variants. Similarly, we find greater molecular change among vaccinated individuals, suggestive of immune evasion. We also observe evidence of transmission in rural areas to follow predictable diffusion processes. Conversely, urban areas are expectedly more complex due to their high mobility, emphasising the role of population structure in driving virus spread. Our analyses highlight the added value of integrating genomic data with detailed demographic and spatial information, particularly in the absence of structured infection surveys.

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

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

Estimating the relative importance of epidemiological and behavioural parameters for epidemic mpox transmission: a modelling study

BACKGROUND

Many European countries experienced outbreaks of mpox in 2022, and there was an mpox outbreak in 2023 in the Democratic Republic of Congo. There were many apparent differences between these outbreaks and previous outbreaks of mpox; the recent outbreaks were observed in men who have sex with men after sexual encounters at common events, whereas earlier outbreaks were observed in a wider population with no identifiable link to sexual contacts. These apparent differences meant that data from previous outbreaks could not reliably be used to parametrise infectious disease models during the 2022 and 2023 mpox outbreaks, and modelling efforts were hampered by uncertainty around key transmission and immunity parameters.

METHODS

We developed a stochastic, discrete-time metapopulation model for mpox that allowed for sexual and non-sexual transmission and the implementation of non-pharmaceutical interventions, specifically contact tracing and pre- and post-exposure vaccinations. We calibrated the model to case data from Berlin and used Sobol sensitivity analysis to identify parameters that mpox transmission is especially sensitive to. We also briefly analysed the sensitivity of the effectiveness of non-pharmaceutical interventions to various efficacy parameters.

RESULTS

We found that variance in the transmission probabilities due to both sexual and non-sexual transmission had a large effect on mpox transmission in the model, as did the level of immunity to mpox conferred by a previous smallpox vaccination. Furthermore, variance in the number of pre-exposure vaccinations offered was the dominant contributor to variance in mpox dynamics in men who have sex with men. If pre-exposure vaccinations were not available, both the accuracy and timeliness of contact tracing had a large impact on mpox transmission in the model.

CONCLUSIONS

Our results are valuable for guiding epidemiological studies for parameter ascertainment and identifying key factors for success of non-pharmaceutical interventions.

Emergence of the reproduction matrix in epidemic forecasting

During the recent COVID-19 pandemic, the instantaneous reproduction number, R(t), has surged as a widely used measure to target public health interventions aiming at curbing the infection rate. In analogy with the basic reproduction number that arises from the linear stability analysis, R(t) is typically interpreted as a threshold parameter that separates exponential growth (R(t) > 1) from exponential decay (R(t) < 1). In real epidemics, however, the finite number of susceptibles, the stratification of the population (e.g. by age or vaccination state), and heterogeneous mixing lead to more complex epidemic courses. In the context of the multidimensional renewal equation, we generalize the scalar R(t) to a reproduction matrix, [Formula: see text], which details the epidemic state of the stratified population, and offers a concise epidemic forecasting scheme. First, the reproduction matrix is computed from the available incidence data (subject to some a priori assumptions), then it is projected into the future by a transfer functional to predict the epidemic course. We demonstrate that this simple scheme allows realistic and accurate epidemic trajectories both in synthetic test cases and with reported incidence data from the COVID-19 pandemic. Accounting for the full heterogeneity and nonlinearity of the infection process, the reproduction matrix improves the prediction of the infection peak. In contrast, the scalar reproduction number overestimates the possibility of sustaining the initial infection rate and leads to an overshoot in the incidence peak. Besides its simplicity, the devised forecasting scheme offers rich flexibility to be generalized to time-dependent mitigation measures, contact rate, infectivity and vaccine protection.

Social contact patterns following the COVID-19 pandemic: a snapshot of post-pandemic behaviour from the CoMix study

The COVID-19 pandemic led to unprecedented changes in behaviour. To estimate if these persisted, a final round of the CoMix social contact survey was conducted in four countries at a time when all societal restrictions had been lifted for several months. We conducted a survey on a nationally representative sample in the UK, Netherlands (NL), Belgium (BE), and Switzerland (CH). Participants were asked about their contacts and behaviours on the previous day. We calculated contact matrices and compared the contact levels to a pre-pandemic baseline to estimate R0. Data collection occurred from 17 November to 7 December 2022. 7477 participants were recruited. Some were asked to undertake the survey on behalf of their children. Only 14.4 % of all participants reported wearing a facemask on the previous day. Self-reported vaccination rates in adults were similar for each country at around 86 %. Trimmed mean recorded contacts were highest in NL with 9.9 (95 % confidence interval [CI] 9.0-10.8) contacts per person per day and lowest in CH at 6.0 (95 % CI 5.4-6.6). Contacts at work were lowest in the UK (1.4 contacts per person per day) and highest in NL at 2.8 contacts per person per day. Other contacts were also lower in the UK at 1.6 per person per day (95 % CI 1.4-1.9) and highest in NL at 3.4 recorded per person per day (95 % CI 43.0-4.0). The next-generation approach suggests that R0 for a close-contact disease would be roughly half pre-pandemic levels in the UK, 80 % in NL and intermediate in the other two countries. The pandemic appears to have resulted in lasting changes in contact patterns expected to have an impact on the epidemiology of many different pathogens. Further post-pandemic surveys are necessary to confirm this finding.

Role of heterogeneity: National scale data-driven agent-based modeling for the US COVID-19 Scenario Modeling Hub

UVA-EpiHiper is a national scale agent-based model to support the US COVID-19 Scenario Modeling Hub (SMH). UVA-EpiHiper uses a detailed representation of the underlying social contact network along with data measured during the course of the pandemic to initialize and calibrate the model. In this paper, we study the role of heterogeneity on model complexity and resulting epidemic dynamics using UVA-EpiHiper. We discuss various sources of heterogeneity that we encounter in the use of UVA-EpiHiper to support modeling and analysis of epidemic dynamics under various scenarios. We also discuss how this affects model complexity and computational complexity of the corresponding simulations. Using round 13 of the SMH as an example, we discuss how UVA-EpiHiper was initialized and calibrated. We then discuss how the detailed output produced by UVA-EpiHiper can be analyzed to obtain interesting insights. We find that despite the complexity in the model, the software, and the computation incurred to an agent-based model in scenario modeling, it is capable of capturing various heterogeneities of real-world systems, especially those in networks and behaviors, and enables analyzing heterogeneities in epidemiological outcomes between different demographic, geographic, and social cohorts. In applying UVA-EpiHiper to round 13 scenario modeling, we find that disease outcomes are different between and within states, and between demographic groups, which can be attributed to heterogeneities in population demographics, network structures, and initial immunity.

Impact of vaccine coverage and disruption to health services on COVID-19 in Ukraine

COVID-19 surveillance in Ukraine ceased after the Russian invasion of the country in 2022, on a background of low vaccination rates of 34.5% for two doses at this time. We conducted a modelling study to estimate the epidemic trajectory of SARS-COV-2 in Ukraine after the start of the war. We use a COVID-19 deterministic Susceptible-Exposed-Infected-Recovered (SEIR) model for Ukraine to estimate the impact of increased vaccination coverage and masking as public health interventions. We fit the model output to case notification data between 6 January and 25 February 2022, then we forecast the COVID-19 epidemic trajectory in different scenarios of mask use and vaccine coverage. In the best-case scenario, 69% of the Ukrainian population would have been infected in the first half of 2022. Increasing mask use from 50 to 80% reduces cases and deaths by 17% and 30% respectively, while increasing vaccination rates to 60% and 9.6% for two and three doses respectively results in a 3% reduction in cases and 28% in deaths. However, if vaccination is increased to a higher coverage of 80% with two doses and 12.8% with three, or mask effectiveness is reduced to 40%, increasing vaccination coverage is more effective. The loss of health services, displacement, and destruction of infrastructure will amplify the risk of COVID-19 in Ukraine and make vaccine programs less feasible. Masks do not need the health infrastructure or cold-chain logistics required for vaccines and are more feasible for rapid epidemic control during war. However, increasing vaccine coverage will save more lives. Vaccination of refugees who have fled to other countries can be more feasibly achieved.

Impact of immunity loss on the optimal vaccination strategy for an age-structured epidemiological model

The pursuit of effective vaccination strategies against COVID-19 remains a critical endeavour in global public health, particularly amidst challenges posed by immunity loss and evolving epidemiological dynamics. This study investigated optimal vaccination strategies by considering age structure, immunity dynamics, and varying maximal vaccination rates. To this end, we formulated an SEIR model stratified into $ n $ age classes, with the vaccination rate as an age-dependent control variable in an optimal control problem. We developed an objective function aimed at minimising critical infections while optimising vaccination efforts and then conducted rigorous mathematical analyses to ensure the existence and characterization of the optimal control. Using data from three countries with diverse age distributions, in expansive, constrictive, and stationary pyramids, we performed numerical simulations to evaluate the optimal age-dependent vaccination strategy, number of critical infections, and vaccination frequency. Our findings highlight the significant influence of maximal vaccination rates on shaping optimal vaccination strategies. Under constant maximal vaccination rates, prioritising age groups based on population demographics proves effective, with higher rates resulting in fewer critically infected individuals across all age distributions. Conversely, adopting age-dependent maximal vaccination rates, akin to the WHO strategy, may not always lead to the lowest critical infection peaks but offers a viable alternative in resource-constrained settings.

Comprehensive profiling of social mixing patterns in resource poor countries: A mixed methods research protocol

BACKGROUND

Low-and-middle-income countries (LMICs) bear a disproportionate burden of communicable diseases. Social interaction data inform infectious disease models and disease prevention strategies. The variations in demographics and contact patterns across ages, cultures, and locations significantly impact infectious disease dynamics and pathogen transmission. LMICs lack sufficient social interaction data for infectious disease modeling.

METHODS

To address this gap, we will collect qualitative and quantitative data from eight study sites (encompassing both rural and urban settings) across Guatemala, India, Pakistan, and Mozambique. We will conduct focus group discussions and cognitive interviews to assess the feasibility and acceptability of our data collection tools at each site. Thematic and rapid analyses will help to identify key themes and categories through coding, guiding the design of quantitative data collection tools (enrollment survey, contact diaries, exit survey, and wearable proximity sensors) and the implementation of study procedures. We will create three age-specific contact matrices (physical, nonphysical, and both) at each study site using data from standardized contact diaries to characterize the patterns of social mixing. Regression analysis will be conducted to identify key drivers of contacts. We will comprehensively profile the frequency, duration, and intensity of infants' interactions with household members using high resolution data from the proximity sensors and calculating infants' proximity score (fraction of time spent by each household member in proximity with the infant, over the total infant contact time) for each household member.

DISCUSSION

Our qualitative data yielded insights into the perceptions and acceptability of contact diaries and wearable proximity sensors for collecting social mixing data in LMICs. The quantitative data will allow a more accurate representation of human interactions that lead to the transmission of pathogens through close contact in LMICs. Our findings will provide more appropriate social mixing data for parameterizing mathematical models of LMIC populations. Our study tools could be adapted for other studies.

Impact assessment of self-medication on COVID-19 prevalence in Gauteng, South Africa, using an age-structured disease transmission modelling framework

OBJECTIVE

To assess the impact of self-medication on the transmission dynamics of COVID-19 across different age groups, examine the interplay of vaccination and self-medication in disease spread, and identify the age group most prone to self-medication.

METHODS

We developed an age-structured compartmentalized epidemiological model to track the early dynamics of COVID-19. Age-structured data from the Government of Gauteng, encompassing the reported cumulative number of cases and daily confirmed cases, were used to calibrate the model through a Markov Chain Monte Carlo (MCMC) framework. Subsequently, uncertainty and sensitivity analyses were conducted on the model parameters.

RESULTS

We found that self-medication is predominant among the age group 15-64 (74.52%), followed by the age group 0-14 (34.02%), and then the age group 65+ (11.41%). The mean values of the basic reproduction number, the size of the first epidemic peak (the highest magnitude of the disease), and the time of the first epidemic peak (when the first highest magnitude occurs) are 4.16499, 241,715 cases, and 190.376 days, respectively. Moreover, we observed that self-medication among individuals aged 15-64 results in the highest spreading rate of COVID-19 at the onset of the outbreak and has the greatest impact on the first epidemic peak and its timing.

CONCLUSION

Studies aiming to understand the dynamics of diseases in areas prone to self-medication should account for this practice. There is a need for a campaign against COVID-19-related self-medication, specifically targeting the active population (ages 15-64).

Uncovering COVID-19 transmission tree: identifying traced and untraced infections in an infection network

Introduction

This paper presents a comprehensive analysis of COVID-19 transmission dynamics using an infection network derived from epidemiological data in South Korea, covering the period from January 3, 2020, to July 11, 2021. The network illustrates infector-infectee relationships and provides invaluable insights for managing and mitigating the spread of the disease. However, significant missing data hinder conventional analysis of such networks from epidemiological surveillance.

Methods

To address this challenge, this article suggests a novel approach for categorizing individuals into four distinct groups, based on the classification of their infector or infectee status as either traced or untraced cases among all confirmed cases. The study analyzes the changes in the infection networks among untraced and traced cases across five distinct periods.

Results

The four types of cases emphasize the impact of various factors, such as the implementation of public health strategies and the emergence of novel COVID-19 variants, which contribute to the propagation of COVID-19 transmission. One of the key findings is the identification of notable transmission patterns in specific age groups, particularly in those aged 20-29, 40-69, and 0-9, based on the four type classifications. Furthermore, we develop a novel real-time indicator to assess the potential for infectious disease transmission more effectively. By analyzing the lengths of connected components, this indicator facilitates improved predictions and enables policymakers to proactively respond, thereby helping to mitigate the effects of the pandemic on global communities.

Conclusion

This study offers a novel approach to categorizing COVID-19 cases, provides insights into transmission patterns, and introduces a real-time indicator for better assessment and management of the disease transmission, thereby supporting more effective public health interventions.

A multiscale modeling framework for Scenario Modeling: Characterizing the heterogeneity of the COVID-19 epidemic in the US

The Scenario Modeling Hub (SMH) initiative provides projections of potential epidemic scenarios in the United States (US) by using a multi-model approach. Our contribution to the SMH is generated by a multiscale model that combines the global epidemic metapopulation modeling approach (GLEAM) with a local epidemic and mobility model of the US (LEAM-US), first introduced here. The LEAM-US model consists of 3142 subpopulations each representing a single county across the 50 US states and the District of Columbia, enabling us to project state and national trajectories of COVID-19 cases, hospitalizations, and deaths under different epidemic scenarios. The model is age-structured, and multi-strain. It integrates data on vaccine administration, human mobility, and non-pharmaceutical interventions. The model contributed to all 17 rounds of the SMH, and allows for the mechanistic characterization of the spatio-temporal heterogeneities observed during the COVID-19 pandemic. Here we describe the mathematical and computational structure of our model, and present the results concerning the emergence of the SARS-CoV-2 Alpha variant (lineage designation B.1.1.7) as a case study. Our findings show considerable spatial and temporal heterogeneity in the introduction and diffusion of the Alpha variant, both at the level of individual states and combined statistical areas, as it competes against the ancestral lineage. We discuss the key factors driving the time required for the Alpha variant to rise to dominance within a population, and quantify the impact that the emergence of the Alpha variant had on the effective reproduction number at the state level. Overall, we show that our multiscale modeling approach is able to capture the complexity and heterogeneity of the COVID-19 pandemic response in the US.

Burden of bacterial antimicrobial resistance in low-income and middle-income countries avertible by existing interventions: an evidence review and modelling analysis

National action plans enumerate many interventions as potential strategies to reduce the burden of bacterial antimicrobial resistance (AMR). However, knowledge of the benefits achievable by specific approaches is needed to inform policy making, especially in low-income and middle-income countries (LMICs) with substantial AMR burden and low health-care system capacity. In a modelling analysis, we estimated that improving infection prevention and control programmes in LMIC health-care settings could prevent at least 337 000 (95% CI 250 200-465 200) AMR-associated deaths annually. Ensuring universal access to high-quality water, sanitation, and hygiene services would prevent 247 800 (160 000-337 800) AMR-associated deaths and paediatric vaccines 181 500 (153 400-206 800) AMR-associated deaths, from both direct prevention of resistant infections and reductions in antibiotic consumption. These estimates translate to prevention of 7·8% (5·6-11·0) of all AMR-associated mortality in LMICs by infection prevention and control, 5·7% (3·7-8·0) by water, sanitation, and hygiene, and 4·2% (3·4-5·1) by vaccination interventions. Despite the continuing need for research and innovation to overcome limitations of existing approaches, our findings indicate that reducing global AMR burden by 10% by the year 2030 is achievable with existing interventions. Our results should guide investments in public health interventions with the greatest potential to reduce AMR burden.

Estimating measures to reduce the transmission of SARS-CoV-2 in Australia to guide a 'National Plan' to reopening

The availability of COVID-19 vaccines promised a reduction in the severity of disease and relief from the strict public health and social measures (PHSMs) imposed in many countries to limit spread and burden of COVID-19. We were asked to define vaccine coverage thresholds for Australia's transition to easing restrictions and reopening international borders. Using evidence of vaccine effectiveness against the then-circulating Delta variant, we used a mathematical model to determine coverage targets. The absence of any COVID-19 infections in many sub-national jurisdictions in Australia posed particular methodological challenges. We used a novel metric called Transmission Potential (TP) as a proxy measure of the population-level effective reproduction number. We estimated TP of the Delta variant under a range of PHSMs, test-trace-isolate-quarantine (TTIQ) efficiencies, vaccination coverage thresholds, and age-based vaccine allocation strategies. We found that high coverage across all ages (≥70%) combined with ongoing TTIQ and minimal PHSMs was sufficient to avoid lockdowns. At lesser coverage (≤60%) rapid case escalation risked overwhelming of the health sector or the need to reimpose stricter restrictions. Maintaining low case numbers was most beneficial for health and the economy, and at higher coverage levels (≥80%) further easing of restrictions was deemed possible. These results directly informed easing of COVID-19 restrictions in Australia.

Modelling respiratory syncytial virus age-specific risk of hospitalisation in term and preterm infants

BACKGROUND

Respiratory syncytial virus (RSV) is the most common cause of acute lower respiratory infections in children worldwide. The highest incidence of severe disease is in the first 6 months of life, with infants born preterm at greatest risk for severe RSV infections. The licensure of new RSV therapeutics (a long-acting monoclonal antibody and a maternal vaccine) in Europe, USA, UK and most recently in Australia, has driven the need for strategic decision making on the implementation of RSV immunisation programs. Data driven approaches, considering the local RSV epidemiology, are critical to advise on the optimal use of these therapeutics for effective RSV control.

METHODS

We developed a dynamic compartmental model of RSV transmission fitted to individually-linked population-based laboratory, perinatal and hospitalisation data for 2000-2012 from metropolitan Western Australia (WA), stratified by age and prior exposure. We account for the differential risk of RSV-hospitalisation in full-term and preterm infants (defined as < 37 weeks gestation). We formulated a function relating age, RSV exposure history, and preterm status to the risk of RSV-hospitalisation given infection.

RESULTS

The age-to-risk function shows that risk of hospitalisation, given RSV infection, declines quickly in the first 12 months of life for all infants and is 2.6 times higher in preterm compared with term infants. The hospitalisation risk, given infection, declines to < 10% of the risk at birth by age 7 months for term infants and by 9 months for preterm infants.

CONCLUSIONS

The dynamic model, using the age-to-risk function, characterises RSV epidemiology for metropolitan WA and can now be extended to predict the impact of prevention measures. The stratification of the model by preterm status will enable the comparative assessment of potential strategies in the extended model that target this RSV risk group relative to all-population approaches. Furthermore, the age-to-risk function developed in this work has wider relevance to the epidemiological characterisation of RSV.

Importance of social inequalities to contact patterns, vaccine uptake, and epidemic dynamics

Individuals' socio-demographic and economic characteristics crucially shape the spread of an epidemic by largely determining the exposure level to the virus and the severity of the disease for those who got infected. While the complex interplay between individual characteristics and epidemic dynamics is widely recognised, traditional mathematical models often overlook these factors. In this study, we examine two important aspects of human behaviour relevant to epidemics: contact patterns and vaccination uptake. Using data collected during the COVID-19 pandemic in Hungary, we first identify the dimensions along which individuals exhibit the greatest variation in their contact patterns and vaccination uptake. We find that generally higher socio-economic groups of the population have a higher number of contacts and a higher vaccination uptake with respect to disadvantaged groups. Subsequently, we propose a data-driven epidemiological model that incorporates these behavioural differences. Finally, we apply our model to analyse the fourth wave of COVID-19 in Hungary, providing valuable insights into real-world scenarios. By bridging the gap between individual characteristics and epidemic spread, our research contributes to a more comprehensive understanding of disease dynamics and informs effective public health strategies.

Vaccination strategies impact the probability of outbreak extinction: A case study of COVID-19 transmission

Mass vaccination has proven to be an effective control measure for mitigating the transmission of infectious diseases. Throughout history, various vaccination strategies have been employed to control infections and terminate outbreaks. In this study, we utilized the transmission of COVID-19 as a case study and constructed a stochastic age-structured compartmental model to investigate the effectiveness of different vaccination strategies. Our analysis focused on estimating the outbreak extinction probability under different vaccination scenarios in both homogeneous and heterogeneous populations. Notably, we found that population heterogeneity can enhance the likelihood of outbreak extinction at varying levels of vaccine coverage. Prioritizing vaccinations for individuals with higher infection risk was found to maximize outbreak extinction probability and reduce overall infections, while allocating vaccines to those with higher mortality risk has been proven more effective in reducing deaths. Moreover, our study highlighted the significance of booster doses as the vaccine effectiveness wanes over time, showing that they can significantly enhance the extinction probability and mitigate disease transmission.

Estimating household contact matrices structure from easily collectable metadata

Contact matrices are a commonly adopted data representation, used to develop compartmental models for epidemic spreading, accounting for the contact heterogeneities across age groups. Their estimation, however, is generally time and effort consuming and model-driven strategies to quantify the contacts are often needed. In this article we focus on household contact matrices, describing the contacts among the members of a family and develop a parametric model to describe them. This model combines demographic and easily quantifiable survey-based data and is tested on high resolution proximity data collected in two sites in South Africa. Given its simplicity and interpretability, we expect our method to be easily applied to other contexts as well and we identify relevant questions that need to be addressed during the data collection procedure.

Morocco's population contact matrices: A crowd dynamics-based approach using aggregated literature data

Estimation of contact patterns is often based on questionnaires and time-use data. The results obtained using these methods have been used extensively over the years and recently to predict the spread of the COVID-19 pandemic. They have also been used to test the effectiveness of non-pharmaceutical measures such as social distance. The latter is integrated into epidemiological models by multiplying contact matrices by control functions. We present a novel method that allows the integration of social distancing and other scenarios such as panic. Our method is based on a modified social force model. The model is calibrated using data relating to the movements of individuals and their interactions such as desired walking velocities and interpersonal distances as well as demographic data. We used the framework to assess contact patterns in different social contexts in Morocco. The estimated matrices are extremely assortative and exhibit patterns similar to those observed in other studies including the POLYMOD project. Our findings suggest social distancing would reduce the numbers of contacts by 95%. Further, we estimated the effect of panic on contact patterns, which indicated an increase in the number of contacts of 11%. This approach could be an alternative to questionnaire-based methods in the study of non-pharmaceutical measures and other specific scenarios such as rush hours. It also provides a substitute for estimating children's contact patterns which are typically assessed through parental proxy reporting in surveys.

Learning from the COVID-19 pandemic: a systematic review of mathematical vaccine prioritization models

As the world becomes ever more connected, the chance of pandemics increases as well. The recent COVID-19 pandemic and the concurrent global mass vaccine roll-out provides an ideal setting to learn from and refine our understanding of infectious disease models for better future preparedness. In this review, we systematically analyze and categorize mathematical models that have been developed to design optimal vaccine prioritization strategies of an initially limited vaccine. As older individuals are disproportionately affected by COVID-19, the focus is on models that take age explicitly into account. The lower mobility and activity level of older individuals gives rise to non-trivial trade-offs. Secondary research questions concern the optimal time interval between vaccine doses and spatial vaccine distribution. This review showcases the effect of various modeling assumptions on model outcomes. A solid understanding of these relationships yields better infectious disease models and thus public health decisions during the next pandemic.

Comparing frequency of booster vaccination to prevent severe COVID-19 by risk group in the United States

There is a public health need to understand how different frequencies of COVID-19 booster vaccines may mitigate the risk of severe COVID-19, while accounting for waning of protection and differential risk by age and immune status. By analyzing United States COVID-19 surveillance and seroprevalence data in a microsimulation model, here we show that more frequent COVID-19 booster vaccination (every 6-12 months) in older age groups and the immunocompromised population would effectively reduce the burden of severe COVID-19, while frequent boosters in the younger population may only provide modest benefit against severe disease. In persons 75+ years, the model estimated that annual boosters would reduce absolute annual risk of severe COVID-19 by 199 (uncertainty interval: 183-232) cases per 100,000 persons, compared to a one-time booster vaccination. In contrast, for persons 18-49 years, the model estimated that annual boosters would reduce this risk by 14 (10-19) cases per 100,000 persons. Those with prior infection had lower benefit of more frequent boosting, and immunocompromised persons had larger benefit. Scenarios with emerging variants with immune evasion increased the benefit of more frequent variant-targeted boosters. This study underscores the benefit of considering key risk factors to inform frequency of COVID-19 booster vaccines in public health guidance and ensuring at least annual boosters in high-risk populations.

Transmission matrices used in epidemiologic modelling

Mixing matrices are included in infectious disease models to reflect transmission opportunities between population strata. These matrices were originally constructed on the basis of theoretical considerations and most of the early work in this area originates from research on sexually transferred diseases in the 80s, in response to AIDS. Later work in the 90s populated these matrices on the basis of survey data gathered to capture transmission risks for respiratory diseases. We provide an overview of developments in the construction of matrices for capturing transmission opportunities in populations. Such transmission matrices are useful for epidemiologic modelling to capture within and between stratum transmission and can be informed from theoretical mixing assumptions, informed by empirical evidence gathered through investigation as well as generated on the basis of data. Links to summary measures and threshold conditions are also provided.

Modeling the impact of distancing measures on infectious disease spread: a case study of COVID-19 in the Moroccan population

This paper explores the impact of various distancing measures on the spread of infectious diseases, focusing on the spread of COVID-19 in the Moroccan population as a case study. Contact matrices, generated through a social force model, capture population interactions within distinct activity locations and age groups. These matrices, tailored for each distancing scenario, have been incorporated into an SEIR model. The study models the region as a network of interconnected activity locations, enabling flexible analysis of the effects of different distancing measures within social contexts and between age groups. Additionally, the method assesses the influence of measures targeting potential superspreaders (i.e., agents with a very high contact rate) and explores the impact of inter-activity location flows, providing insights beyond scalar contact rates or survey-based contact matrices. The results suggest that implementing intra-activity location distancing measures significantly reduces in the number of infected individuals relative to the act of imposing restrictions on individuals with a high contact rate in each activity location. The combination of both measures proves more advantageous. On a regional scale, characterized as a network of interconnected activity locations, restrictions on the movement of individuals with high contact rates was found to result in a $ 2 \% $ reduction, while intra-activity location-based distancing measures was found to achieve a $ 44 \% $ reduction. The combination of these two measures yielded a $ 48\% $ reduction.

Impact of alternative Non-Pharmaceutical Interventions strategies for controlling COVID-19 outbreak in Bangladesh: A modeling study

The COVID-19 pandemic has been a major health concern in Bangladesh until very recently. Although the Bangladesh government has employed various infection control strategies, more targeted Non-Pharmaceutical interventions (NPIs), including school closure, mask-wearing, hand washing, and social distancing have gained special attention. Despite significant long-term adverse effects of school closures, authorities have opted to keep schools closed to curb the spread of COVID-19 infection. However, there is limited knowledge about the impact of reopening schools alongside other NPI measures on the course of the epidemic. In this study, we implemented a mathematical modeling framework developed by the CoMo Consortium to explore the impact of NPIs on the dynamics of the COVID-19 outbreak and deaths for Bangladesh. For robustness, the results of prediction models are then validated through model calibration with incidence and mortality data and using external sources. Hypothetical projections are made under alternative NPIs where we compare the impact of current NPIs with school closures versus enhanced NPIs with school openings. Results suggest that enhanced NPIs with schools opened may have lower COVID-19 related prevalence and deaths. This finding indicates that enhanced NPIs and school openings may mitigate the long-term negative impacts of COVID-19 in low- and middle-income countries. Potential shortcomings and ways to improve the research are also discussed.

Seasonal variations in social contact patterns in a rural population in north India: Implications for pandemic control

Social contact mixing patterns are critical to model the transmission of communicable diseases, and have been employed to model disease outbreaks including COVID-19. Nonetheless, there is a paucity of studies on contact mixing in low and middle-income countries such as India. Furthermore, mathematical models of disease outbreaks do not account for the temporal nature of social contacts. We conducted a longitudinal study of social contacts in rural north India across three seasons and analysed the temporal differences in contact patterns. A contact diary survey was performed across three seasons from October 2015-16, in which participants were queried on the number, duration, and characteristics of contacts that occurred on the previous day. A total of 8,421 responses from 3,052 respondents (49% females) recorded characteristics of 180,073 contacts. Respondents reported a significantly higher number and duration of contacts in the winter, followed by the summer and the monsoon season (Nemenyi post-hoc, p<0.001). Participants aged 0-9 years and 10-19 years of age reported the highest median number of contacts (16 (IQR 12-21), 17 (IQR 13-24) respectively) and were found to have the highest node centrality in the social network of the region (pageranks = 0.20, 0.17). A large proportion (>80%) of contacts that were reported in schools or on public transport involved physical contact. To the best of our knowledge, our study is the first from India to show that contact mixing patterns vary by the time of the year and provides useful implications for pandemic control. We compared the differences in the number, duration and location of contacts by age-group and gender, and studied the impact of the season, age-group, employment and day of the week on the number and duration of contacts using multivariate negative binomial regression. We created a social network to further understand the age and gender-specific contact patterns, and used the contact matrices in each season to parameterise a nine-compartment agent-based model for simulating a COVID-19 epidemic in each season. Our results can be used to parameterize more accurate mathematical models for prediction of epidemiological trends of infections in rural India.

Examining the Influence of Imbalanced Social Contact Matrices in Epidemic Models

Transmissible infections such as those caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spread according to who contacts whom. Therefore, many epidemic models incorporate contact patterns through contact matrices. Contact matrices can be generated from social contact survey data. However, the resulting matrices are often imbalanced, such that the total number of contacts reported by group A with group B do not match those reported by group B with group A. We examined the theoretical influence of imbalanced contact matrices on the estimated basic reproduction number (R0). We then explored how imbalanced matrices may bias model-based epidemic projections using an illustrative simulation model of SARS-CoV-2 with 2 age groups (<15 and ≥15 years). Models with imbalanced matrices underestimated the initial spread of SARS-CoV-2, had later time to peak incidence, and had smaller peak incidence. Imbalanced matrices also influenced cumulative infections observed per age group, as well as the estimated impact of an age-specific vaccination strategy. Stratified transmission models that do not consider contact balancing may generate biased projections of epidemic trajectory and the impact of targeted public health interventions. Therefore, modeling studies should implement and report methods used to balance contact matrices for stratified transmission models.

Modelling the Cost-Effectiveness of Hepatitis A in South Africa

The World Health Organization (WHO) recommends the consideration of introducing routine hepatitis A vaccination into national immunization schedules for children ≥ 1 years old in countries with intermediate HAV endemicity. Recent data suggest that South Africa is transitioning from high to intermediate HAV endemicity, thus it is important to consider the impact and cost of potential routine hepatitis A vaccination strategies in the country. An age-structured compartmental model of hepatitis A transmission was calibrated with available data from South Africa, incorporating direct costs of hepatitis A treatment and vaccination. We used the calibrated model to evaluate the impact and costs of several childhood hepatitis A vaccination scenarios from 2023 to 2030. We assessed how each scenario impacted the burden of hepatitis A (symptomatic hepatitis A cases and mortality) as well as calculated the incremental cost per DALY averted as compared to the South African cost-effectiveness threshold. All costs and outcomes were discounted at 5%. For the modelled scenarios, the median estimated cost of the different vaccination strategies ranged from USD 1.71 billion to USD 2.85 billion over the period of 2023 to 2030, with the cost increasing for each successive scenario and approximately 39-52% of costs being due to vaccination. Scenario 1, which represented the administration of one dose of the hepatitis A vaccine in children < 2 years old, requires approximately 5.3 million vaccine doses over 2023-2030 and is projected to avert a total of 136,042 symptomatic cases [IQR: 88,842-221,483] and 31,106 [IQR: 22,975-36,742] deaths due to hepatitis A over the period of 2023 to 2030. The model projects that Scenario 1 would avert 8741 DALYs over the period of 2023 to 2030; however, it is not cost-effective against the South African cost-effectiveness threshold with an ICER per DALY averted of USD 21,006. While Scenario 3 and 4 included the administration of more vaccine doses and averted more symptomatic cases of hepatitis A, these scenarios were absolutely dominated owing to the population being infected before vaccination through the mass campaigns at older ages. The model was highly sensitive to variation of access to liver transplant in South Africa. When increasing the access to liver transplant to 100% for the baseline and Scenario 1, the ICER for Scenario 1 becomes cost-effective against the CET (ICER = USD 2425). Given these findings, we recommend further research is conducted to understand the access to liver transplants in South Africa and better estimate the cost of liver transplant care for hepatitis A patients. The modelling presented in this paper has been used to develop a user-friendly application for vaccine policy makers to further interrogate the model outcomes and consider the costs and benefits of introducing routine hepatitis A vaccination in South Africa.

Modelling the population-level benefits and cost-effectiveness of cell-based quadrivalent influenza vaccine for children and adolescents aged 6 months to 17 years in the US

BACKGROUND

Cell-based quadrivalent influenza vaccines (QIVc) can increase effectiveness against seasonal influenza by avoiding mismatch from egg adaption of vaccine viruses. This study evaluates the population-level cost-effectiveness and impacts on health outcomes of QIVc versus an egg-based vaccine (QIVe) in children aged 6 months to 17 years in the US.

RESEARCH DESIGN AND METHODS

A dynamic age-structured susceptible-exposed-infected-recovered model was used to simulate influenza transmission in low and high incidence seasons for two scenarios: 1. QIVe for 6 months-17 year-olds, QIVc for 18-64 year-olds, and adjuvanted QIV (aQIV) for ≥ 65 year-olds, and 2. QIVc for 6 months-64 year-olds, and aQIV for ≥ 65 year-olds. Probabilistic sensitivity analysis was performed to account for uncertainty in parameter estimates. Cost-effectiveness was evaluated as incremental cost-effectiveness ratios (ICERs).

RESULTS

Extension of QIVc to children resulted in 3-4% reductions in cases (1,656,271), hospitalizations (16,688), and deaths (2,126) at a population level in a high incidence season, and 65% reductions (cases: 2,856,384; hospitalizations: 31667; deaths: 4,163) in a low incidence season. Use of QIVc would be cost-saving, with ICERs of -$16,427/QALY and -$8,100/QALY from a payer perspective and -$22,669/QALY and -$15,015/QALY from a societal perspective, for low and high incidence seasons respectively. Cost savings were estimated at approximately $468 million and $1.366 billion for high and low incidence seasons, respectively.

CONCLUSION

Use of QIVc instead of QIVe in children > 6 months of age in the US would reduce the disease burden and be cost-saving from both a payer and societal perspective.

Planning for healthcare services during the COVID-19 pandemic in the Southeast of England: a system dynamics modelling approach

OBJECTIVES

To develop, test, validate and implement a system dynamics model to simulate the pandemic progress and the impact of various interventions on viral spread, healthcare utilisation and demand in secondary care.

DESIGN

We adopted the system dynamics model incorporating susceptible, exposed, infection and recovery framework to simulate the progress of the pandemic and how the interventions for the COVID-19 response influence the outcomes with a focus on secondary care.

SETTING

This study was carried out covering all the local health systems in Southeast of England with a catchment population of six million with a specific focus on Kent and Medway system.

PARTICIPANTS

Six local health systems in Southeast of England using Kent and Medway as a case study.

INTERVENTIONS

Short to medium 'what if' scenarios incorporating human behaviour, non-pharmaceutical interventions and medical interventions were tested using the model with regular and continuous feedback of the model results to the local health system leaders for monitoring, planning and rapid response as needed.

MAIN OUTCOME MEASURES

Daily output from the model which included number infected in the population, hospital admissions needing COVID-19 care, occupied general beds, continuous positive airway pressure beds, intensive care beds, hospital discharge pathways and deaths.

RESULTS

We successfully implemented a regional series of models based on the local population needs which were used in healthcare planning as part of the pandemic response.

CONCLUSIONS

In this study, we have demonstrated the utility of system dynamics modelling incorporating local intelligence and collaborative working during the pandemic to respond rapidly and take decisions to protect the population. This led to strengthened cooperation among partners and ensured that the local population healthcare needs were met.

Canadian Covid-19 pandemic public health mitigation measures at the province level

The Covid-19 pandemic has prompted governments across the world to enforce a range of public health interventions. We introduce the Covid-19 Policy Response Canadian tracker (CPRCT) database that tracks and records implemented public health measures in every province and territory in Canada. The implementations are recorded on a four-level ordinal scale (0-3) for three domains, (Schools, Work, and Other), capturing differences in degree of response. The data-set allows the exploration of the effects of public health mitigation on the spread of Covid-19, as well as provides a near-real-time record in an accessible format that is useful for a diverse range of modeling and research questions.

Data-driven assessment of immune evasion and dynamic Zero-COVID policy on fast-spreading Omicron in Changchun

Due to its immune evasion capability, the SARS-CoV-2 Omicron variant was declared a variant of concern by the World Health Organization. The spread of Omicron in Changchun (i.e., the capital of Jilin province in northeast of China) during the spring of 2022 was successfully curbed under the strategy of a dynamic Zero-COVID policy. To evaluate the impact of immune evasion on vaccination and other measures, and to understand how the dynamic Zero-COVID measure stopped the epidemics in Changchun, we establish a compartmental model over different stages and parameterized the model with actual reported data. The model simulation firstly shows a reasonably good fit between our model prediction and the data. Second, we estimate the testing rate in the early stage of the outbreak to reveal the real infection size. Third, numerical simulations show that the coverage of vaccine immunization in Changchun and the regular nucleic acid testing could not stop the epidemic, while the 'non-pharmaceutical' intervention measures utilized in the dynamic Zero-COVID policy could play significant roles in the containment of Omicron. Based on the parameterized model, numerical analysis demonstrates that if one wants to achieve epidemic control by fully utilizing the effect of 'dynamic Zero-COVID' measures, therefore social activities are restricted to the minimum level, and then the economic development may come to a halt. The insight analysis in this work could provide reference for infectious disease prevention and control measures in the future.

Comprehensive profiling of social mixing patterns in resource poor countries: a mixed methods research protocol

Background

Low-and-middle-income countries (LMICs) bear a disproportionate burden of communicable diseases. Social interaction data inform infectious disease models and disease prevention strategies. The variations in demographics and contact patterns across ages, cultures, and locations significantly impact infectious disease dynamics and pathogen transmission. LMICs lack sufficient social interaction data for infectious disease modeling.

Methods

To address this gap, we will collect qualitative and quantitative data from eight study sites (encompassing both rural and urban settings) across Guatemala, India, Pakistan, and Mozambique. We will conduct focus group discussions and cognitive interviews to assess the feasibility and acceptability of our data collection tools at each site. Thematic and rapid analyses will help to identify key themes and categories through coding, guiding the design of quantitative data collection tools (enrollment survey, contact diaries, exit survey, and wearable proximity sensors) and the implementation of study procedures.We will create three age-specific contact matrices (physical, nonphysical, and both) at each study site using data from standardized contact diaries to characterize the patterns of social mixing. Regression analysis will be conducted to identify key drivers of contacts. We will comprehensively profile the frequency, duration, and intensity of infants' interactions with household members using high resolution data from the proximity sensors and calculating infants' proximity score (fraction of time spent by each household member in proximity with the infant, over the total infant contact time) for each household member.

Discussion

Our qualitative data yielded insights into the perceptions and acceptability of contact diaries and wearable proximity sensors for collecting social mixing data in LMICs. The quantitative data will allow a more accurate representation of human interactions that lead to the transmission of pathogens through close contact in LMICs. Our findings will provide more appropriate social mixing data for parameterizing mathematical models of LMIC populations. Our study tools could be adapted for other studies.

Real-time estimation of the effective reproduction number of COVID-19 from behavioral data

Monitoring the effective reproduction number [Formula: see text] of a rapidly unfolding pandemic in real-time is key to successful mitigation and prevention strategies. However, existing methods based on case numbers, hospital admissions or fatalities suffer from multiple measurement biases and temporal lags due to high test positivity rates or delays in symptom development or administrative reporting. Alternative methods such as web search and social media tracking are less directly indicating epidemic prevalence over time. We instead record age-stratified anonymous contact matrices at a daily resolution using a longitudinal online-offline survey in Hungary during the first two waves of the COVID-19 pandemic. This approach is innovative, cheap, and provides information in near real-time for estimating [Formula: see text] at a daily resolution. Moreover, it allows to complement traditional surveillance systems by signaling periods when official monitoring infrastructures are unreliable due to observational biases.

Evaluation of vaccine rollout strategies for emerging infectious diseases: A model-based approach including protection attitudes

Vaccine allocation strategies become crucial during vaccine shortages, especially in the face of potential outbreaks of new infectious diseases, as witnessed during the COVID-19 pandemic. To address this, a specialized compartmental model is created, which simulates an emerging infectious disease similar to COVID-19. This model divides the population into different age groups and is used to compare various vaccine prioritisation approaches, aiming to minimize the total number of fatalities. The model is an improvement upon previous ones as it incorporates essential behavioural factors and is adapted to account for the protective effects of vaccination against both disease infection and transmission. It takes into account human behaviors such as mask-wearing and social distancing by utilizing specific parameters related to self-protection, awareness levels, and the frequency of daily person-to-person interactions within each age group. Furthermore, a novel method for dynamic vaccine prioritisation was introduced in this study. This approach is model-independent and relies on the dynamic R number. It is the first time such a method has been developed, offering a decision-making approach that is not tied to any specific model. This innovation provides a flexible and adaptable strategy for determining vaccine priorities based on real-time data and the current state of the outbreak. Our findings reveal crucial insights into vaccine allocation strategies. When the daily rollout rates are fast (0.75% or higher) and children are eligible for vaccination, prioritising groups with high daily person-to-person interactions can lead to substantial reductions in total fatalities (up to approximately 40% lower). On the other hand, if rollout rates are slower and overall vaccination coverage is high, focusing on vaccinating elders emerges as the most effective strategy, resulting in up to approximately 10% fewer fatalities. However, the scenario changes significantly when children are not eligible for vaccination, as they constitute a highly interactive population group. In this case, the differences between priority strategies become smaller. With fast daily rollout rates, prioritisation based on interactions achieves only a 7% reduction in total fatalities, while a slower rollout with vaccination of elders first leads to an approximately 11% reduction in fatalities compared to the scenario where children are eligible for vaccination. The impact of behavioural parameters is equally critical. When the self-protection levels exercised by the population are low, it significantly affects the optimal vaccine prioritisation strategy to be followed, making it essential to consider behavioural factors in decision-making.

Modeling vaccination strategies with limited early COVID-19 vaccine access in low- and middle-income countries: A case study of Thailand

Low- and middle-income countries faced significant challenges in accessing COVID-19 vaccines during the early stages of the pandemic. In this study, we utilized an age-structured modeling approach to examine the implications of various vaccination strategies, vaccine prioritization, and vaccine rollout speeds in Thailand, an upper-middle-income country experiencing vaccine shortages during the early stages of the pandemic. The model directly compares the effectiveness of several vaccination strategies, including the heterologous vaccination where CoronaVac (CV) vaccine was administered as the first dose, followed by ChAdOx1 nCoV-19 (AZ) vaccine as the second dose, under varying disease transmission dynamics. We found that the traditional AZ homologous vaccination was more effective than the CV homologous vaccination, regardless of disease transmission dynamics. However, combining CV and AZ vaccines via either parallel homologous or heterologous vaccinations was more effective than relying solely on AZ homologous vaccination. Additionally, prioritizing vaccination for the elderly aged 60 years and above was the most effective way to reduce mortality when community transmission is well-controlled. On the other hand, prioritizing workers aged 20-59 was most effective in lowering COVID-19 cases, irrespective of the transmission dynamics. Lastly, despite the vaccine prioritization strategy, rapid vaccine rollout speeds were crucial in reducing COVID-19 infections and deaths. These findings suggested that in low- and middle-income countries where early access to high-efficacy vaccines might be limited, obtaining any accessible vaccines as early as possible and using them in parallel with other higher-efficacy vaccines might be a better strategy than waiting for and relying solely on higher-efficacy vaccines.

Modeling the impact of COVID-19 nonpharmaceutical interventions on respiratory syncytial virus transmission in South Africa

Background

The South African government employed various nonpharmaceutical interventions (NPIs) to reduce the spread of SARS-CoV-2. Surveillance data from South Africa indicates reduced circulation of respiratory syncytial virus (RSV) throughout the 2020-2021 seasons. Here, we use a mechanistic transmission model to project the rebound of RSV in the two subsequent seasons.

Methods

We fit an age-structured epidemiological model to hospitalization data from national RSV surveillance in South Africa, allowing for time-varying reduction in RSV transmission during periods of COVID-19 circulation. We apply the model to project the rebound of RSV in the 2022 and 2023 seasons.

Results

We projected an early and intense outbreak of RSV in April 2022, with an age shift to older infants (6-23 months old) experiencing a larger portion of severe disease burden than typical. In March 2022, government alerts were issued to prepare the hospital system for this potentially intense outbreak. We then assess the 2022 predictions and project the 2023 season. Model predictions for 2023 indicate that RSV activity has not fully returned to normal, with a projected early and moderately intense wave. We estimate that NPIs reduced RSV transmission between 15% and 50% during periods of COVID-19 circulation.

Conclusions

A wide range of NPIs impacted the dynamics of the RSV outbreaks throughout 2020-2023 in regard to timing, magnitude, and age structure, with important implications in a low- and middle-income countries (LMICs) setting where RSV interventions remain limited. More efforts should focus on adapting RSV models to LMIC data to project the impact of upcoming medical interventions for this disease.

Seasonality as a driver of pH1N12009 influenza vaccination campaign impact

Although the most recent respiratory virus pandemic was triggered by a Coronavirus, sustained and elevated prevalence of highly pathogenic avian influenza viruses able to infect mammalian hosts highlight the continued threat of pandemics of influenza A virus (IAV) to global health. Retrospective analysis of pandemic outcomes, including comparative investigation of intervention efficacy in different regions, provide important contributions to the evidence base for future pandemic planning. The swine-origin IAV pandemic of 2009 exhibited regional variation in onset, infection dynamics and annual infection attack rates (IARs). For example, the UK experienced three severe peaks of infection over two influenza seasons, whilst Australia experienced a single severe wave. We adopt a seasonally forced 2-subtype model for the transmission of pH1N12009 and seasonal H3N2 to examine the role vaccination campaigns may play in explaining differences in pandemic trajectories in temperate regions. Our model differentiates between the nature of vaccine- and infection-acquired immunity. In particular, we assume that immunity triggered by infection elicits heterologous cross-protection against viral shedding in addition to long-lasting neutralising antibody, whereas vaccination induces imperfect reduction in susceptibility. We employ an Approximate Bayesian Computation (ABC) framework to calibrate the model using data for pH1N12009 seroprevalence, relative subtype dominance, and annual IARs for Australia and the UK. Heterologous cross-protection substantially suppressed the pandemic IAR over the posterior, with the strength of protection against onward transmission inversely correlated with the initial reproduction number. We show that IAV pandemic timing relative to the usual seasonal influenza cycle influenced the size of the initial waves of pH1N12009 in temperate regions and the impact of vaccination campaigns.

Changing social contact patterns among US workers during the COVID-19 pandemic: April 2020 to December 2021

Non-pharmaceutical interventions minimize social contacts, hence the spread of respiratory pathogens such as influenza and SARS-CoV-2. Globally, there is a paucity of social contact data from the workforce. In this study, we quantified two-day contact patterns among USA employees. Contacts were defined as face-to-face conversations, involving physical touch or proximity to another individual and were collected using electronic self-kept diaries. Data were collected over 4 rounds from 2020 to 2021 during the COVID-19 pandemic. Mean (standard deviation) contacts reported by 1456 participants were 2.5 (2.5), 8.2 (7.1), 9.2 (7.1) and 10.1 (9.5) across round 1 (April-June 2020), 2 (November 2020-January 2021), 3 (June-August 2021), and 4 (November-December 2021), respectively. Between round 1 and 2, we report a 3-fold increase in the mean number of contacts reported per participant with no major increases from round 2-4. We then modeled SARS-CoV-2 transmission at home, work, and community settings. The model revealed reduced relative transmission in all settings in round 1. Subsequently, transmission increased at home and in the community but remained exceptionally low in work settings. To accurately parameterize models of infection transmission and control, we need empirical social contact data that capture human mixing behavior across time.

A mixed-method approach to determining contact matrices in the Cox's Bazar refugee settlement

Contact matrices are an important ingredient in age-structured epidemic models to inform the simulated spread of the disease between subgroups of the population. These matrices are generally derived using resource-intensive diary-based surveys and few exist in the Global South or tailored to vulnerable populations. In particular, no contact matrices exist for refugee settlements-locations under-served by epidemic models in general. In this paper, we present a novel, mixed-method approach for deriving contact matrices in populations, which combines a lightweight, rapidly deployable survey with an agent-based model of the population informed by census and behavioural data. We use this method to derive the first set of contact matrices for the Cox's Bazar refugee settlement in Bangladesh. To validate our approach, we apply it to the UK population and compare our derived matrices with well-known contact matrices collected using traditional methods. Our findings demonstrate that our mixed-method approach successfully addresses some of the challenges faced by traditional and agent-based approaches to deriving contact matrices. It also shows potential for implementation in resource-constrained environments. This work therefore contributes to a broader aim of developing new methods and mechanisms of data collection for modelling disease spread in refugee and internally displaced person (IDP) settlements and better serving these vulnerable communities.

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.

Predictions from standard epidemiological models of consequences of segregating and isolating vulnerable people into care facilities

OBJECTIVES

Since the declaration of the COVID-19 pandemic, many governments have imposed policies to reduce contacts between people who are presumed to be particularly vulnerable to dying from respiratory illnesses and the rest of the population. These policies typically address vulnerable individuals concentrated in centralized care facilities and entail limiting social contacts with visitors, staff members, and other care home residents. We use a standard epidemiological model to investigate the impact of such circumstances on the predicted infectious disease attack rates, for interacting robust and vulnerable populations.

METHODS

We implement a general susceptible-infectious-recovered (SIR) compartmental model with two populations: robust and vulnerable. The key model parameters are the per-individual frequencies of within-group (robust-robust and vulnerable-vulnerable) and between-group (robust-vulnerable and vulnerable-robust) infectious-susceptible contacts and the recovery times of individuals in the two groups, which can be significantly longer for vulnerable people.

RESULTS

Across a large range of possible model parameters including degrees of segregation versus intermingling of vulnerable and robust individuals, we find that concentrating the most vulnerable into centralized care facilities virtually always increases the infectious disease attack rate in the vulnerable group, without significant benefit to the robust group.

CONCLUSIONS

Isolated care homes of vulnerable residents are predicted to be the worst possible mixing circumstances for reducing harm in epidemic or pandemic conditions.

Data-driven models for the risk of infection and hospitalization during a pandemic: Case study on COVID-19 in Nepal

The newly emerging pandemic disease often poses unexpected troubles and hazards to the global health system, particularly in low and middle-income countries like Nepal. In this study, we developed mathematical models to estimate the risk of infection and the risk of hospitalization during a pandemic which are critical for allocating resources and planning health policies. We used our models in Nepal's unique data set to explore national and provincial-level risks of infection and risk of hospitalization during the Delta and Omicron surges. Furthermore, we used our model to identify the effectiveness of non-pharmaceutical interventions (NPIs) to mitigate COVID-19 in various groups of people in Nepal. Our analysis shows no significant difference in reproduction numbers in provinces between the Delta and Omicron surge periods, but noticeable inter-provincial disparities in the risk of infection (for example, during Delta (Omicron) surges, the risk of infection of Bagmati province is: ∼ 98.94 (89.62); Madhesh province: ∼ 12.16 (5.1); Karnali province ∼31.16 (3) per hundred thousands). Our estimates show a significantly low level of hospitalization risk during the Omicron surge compared to the Delta surge (hospitalization risk is: ∼10% in Delta and ∼2.5% in Omicron). We also found significant inter-provincial disparities in the hospitalization rate (for example, ∼ 6% in Madhesh province and ∼ 21% in Sudur Paschim) during the Delta surge. Moreover, our results show that closing only schools, colleges, and workplaces reduces the risk of infection by one-third, while a complete lockdown reduces the infections by two-thirds. Our study provides a framework for the computation of the risk of infection and the risk of hospitalization and offers helpful information for controlling the pandemic.

Rational social distancing policy during epidemics with limited healthcare capacity

Epidemics of infectious diseases posing a serious risk to human health have occurred throughout history. During recent epidemics there has been much debate about policy, including how and when to impose restrictions on behaviour. Policymakers must balance a complex spectrum of objectives, suggesting a need for quantitative tools. Whether health services might be 'overwhelmed' has emerged as a key consideration. Here we show how costly interventions, such as taxes or subsidies on behaviour, can be used to exactly align individuals' decision making with government preferences even when these are not aligned. In order to achieve this, we develop a nested optimisation algorithm of both the government intervention strategy and the resulting equilibrium behaviour of individuals. We focus on a situation in which the capacity of the healthcare system to treat patients is limited and identify conditions under which the disease dynamics respect the capacity limit. We find an extremely sharp drop in peak infections at a critical maximum infection cost in the government's objective function. This is in marked contrast to the gradual reduction of infections if individuals make decisions without government intervention. We find optimal interventions vary less strongly in time when interventions are costly to the government and that the critical cost of the policy switch depends on how costly interventions are.

Cost-effectiveness analysis of COVID-19 variants effects in an age-structured model

This study analyzes the impact of COVID-19 variants on cost-effectiveness across age groups, considering vaccination efforts and nonpharmaceutical interventions in Republic of Korea. We aim to assess the costs needed to reduce COVID-19 cases and deaths using age-structured model. The proposed age-structured model analyzes COVID-19 transmission dynamics, evaluates vaccination effectiveness, and assesses the impact of the Delta and Omicron variants. The model is fitted using data from the Republic of Korea between February 2021 and November 2022. The cost-effectiveness of interventions, medical costs, and the cost of death for different age groups are evaluated through analysis. The impact of different variants on cases and deaths is also analyzed, with the Omicron variant increasing transmission rates and decreasing case-fatality rates compared to the Delta variant. The cost of interventions and deaths is higher for older age groups during both outbreaks, with the Omicron outbreak resulting in a higher overall cost due to increased medical costs and interventions. This analysis shows that the daily cost per person for both the Delta and Omicron variants falls within a similar range of approximately $10-$35. This highlights the importance of conducting cost-effect analyses when evaluating the impact of COVID-19 variants.

A model-based strategy for the COVID-19 vaccine roll-out in the Philippines

COVID-19 has affected millions of people worldwide, causing illness and death, and disrupting daily life while imposing a significant social and economic burden. Vaccination is an important control measure that significantly reduces mortality if properly and efficiently distributed. In this work, an age-structured model of COVID-19 transmission, incorporating an unreported infectious compartment, is developed. Three age groups are considered: young (0-19 years), adult (20-64 years), and elderly (65+ years). The transmission rate and reporting rate are determined for each group by utilizing the number of COVID-19 cases in the National Capital Region in the Philippines. Optimal control theory is employed to identify the best vaccine allocation to different age groups. Further, three different vaccination periods are considered to reflect phases of vaccination priority groups: the first, second, and third account for the inoculation of the elderly, adult and elderly, and all three age groups, respectively. This study could guide in making informed decisions in mitigating a population-structured disease transmission under limited resources.

Circulation of Respiratory Viruses in the City: Towards an Agent-Based Ecosystem model

Mathematical models play an important role in management of outbreaks of acute respiratory infections (ARI). While such models are generally used to study the spread of a solitary virus, in reality multiple viruses co-circulate in the population. These viruses have been studied in detail, including the course of infection and immune defense mechanisms. We developed an agent-based model, called ABM-ARI, assimilating heterogeneous data and theoretical knowledge into a biologically motivated system, that allows to reproduce the seasonal patterns of ARI incidence and simulate interventions. ABM-ARI uses city-specific data to create a synthetic population and to construct realistic contact networks in different activity settings. Characteristics of infection, immune protection and non-specific resistance were varied between individuals to account for the population heterogeneity. For the calibration, we minimised the normalised mean absolute error between simulated and observed epidemic curves. ABM-ARI was built based on the quantitative assessment of features of predominant respiratory viruses and epidemiological characteristics of the population. It provides a good fit to the observed epidemic curves for different age groups and viruses. We also simulated one-week school closures when student absences were at or above 10%, 20% or 30% and found that only 10% and 20% thresholds resulted in a reduction of the incidence. ABM-ARI has a great potential in tackling the challenge of emerging infections by simulating and evaluating the effectiveness of various interventions.

Effective vaccination strategies to control COVID-19 in Korea: a modeling study

OBJECTIVES

In Korea, as immunity levels of the coronavirus disease 2019 (COVID-19) in the population acquired through previous infections and vaccinations have decreased, booster vaccinations have emerged as a necessary measure to control new outbreaks. The objective of this study was to identify the most suitable vaccination strategy for controlling the surge in COVID-19 cases.

METHODS

A mathematical model was developed to concurrently evaluate the immunity levels induced by vaccines and infections. This model was then employed to investigate the potential for future resurgence and the possibility of control through the use of vaccines and antivirals.

RESULTS

As of May 11, 2023, if the current epidemic trend persists without further vaccination efforts, a peak in resurgence is anticipated to occur around mid-October of the same year. Under the most favorable circumstances, the peak number of severely hospitalized patients could be reduced by 43% (n=480) compared to the scenario without vaccine intervention (n=849). Depending on outbreak trends and vaccination strategies, the best timing for vaccination in terms of minimizing this peak varies from May 2023 to August 2023.

CONCLUSIONS

Our findings suggest that if the epidemic persist, the best timing for administering vaccinations would need to be earlier than currently outlined in the Korean plan. It is imperative to continue monitoring outbreak trends, as this is key to determining the best vaccination timing in order to manage potential future surges.