Improving modelling for epidemic responses: reflections from members of the UK infectious disease modelling community on their experiences during the COVID-19 pandemic.

Background: The COVID-19 pandemic both relied and placed significant burdens on the experts involved from research and public health sectors. The sustained high pressure of a pandemic on responders, such as healthcare workers, can lead to lasting psychological impacts including acute stress disorder, post-traumatic stress disorder, burnout, and moral injury, which can impact individual wellbeing and productivity. Methods: As members of the infectious disease modelling community, we convened a reflective workshop to understand the professional and personal impacts of response work on our community and to propose recommendations for future epidemic responses. The attendees represented a range of career stages, institutions, and disciplines. This piece was collectively produced by those present at the session based on our collective experiences. Results: Key issues we identified at the workshop were lack of institutional support, insecure contracts, unequal credit and recognition, and mental health impacts. Our recommendations include rewarding impactful work, fostering academia-public health collaboration, decreasing dependence on key individuals by developing teams, increasing transparency in decision-making, and implementing sustainable work practices. Conclusions: Despite limitations in representation, this workshop provided valuable insights into the UK COVID-19 modelling experience and guidance for future public health crises. Recognising and addressing the issues highlighted is crucial, in our view, for ensuring the effectiveness of epidemic response work in the future.

Aetiology of vaginal discharge, urethral discharge, and genital ulcer in sub-Saharan Africa: A systematic review and meta-regression.

BACKGROUND: Syndromic management is widely used to treat symptomatic sexually transmitted infections in settings without aetiologic diagnostics. However, underlying aetiologies and consequent treatment suitability are uncertain without regular assessment. This systematic review estimated the distribution, trends, and determinants of aetiologies for vaginal discharge, urethral discharge, and genital ulcer in sub-Saharan Africa (SSA). METHODS AND FINDINGS: We searched Embase, MEDLINE, Global Health, Web of Science, and grey literature from inception until December 20, 2023, for observational studies reporting aetiologic diagnoses among symptomatic populations in SSA. We adjusted observations for diagnostic test performance, used generalised linear mixed-effects meta-regressions to generate estimates, and critically appraised studies using an adapted Joanna Briggs Institute checklist. Of 4,418 identified records, 206 reports were included from 190 studies in 32 countries conducted between 1969 and 2022. In 2015, estimated primary aetiologies for vaginal discharge were candidiasis (69.4% [95% confidence interval (CI): 44.3% to 86.6%], n = 50), bacterial vaginosis (50.0% [95% CI: 32.3% to 67.8%], n = 39), chlamydia (16.2% [95% CI: 8.6% to 28.5%], n = 50), and trichomoniasis (12.9% [95% CI: 7.7% to 20.7%], n = 80); for urethral discharge were gonorrhoea (77.1% [95% CI: 68.1% to 84.1%], n = 68) and chlamydia (21.9% [95% CI: 15.4% to 30.3%], n = 48); and for genital ulcer were herpes simplex virus type 2 (HSV-2) (48.3% [95% CI: 32.9% to 64.1%], n = 47) and syphilis (9.3% [95% CI: 6.4% to 13.4%], n = 117). Temporal variation was substantial, particularly for genital ulcer where HSV-2 replaced chancroid as the primary cause. Aetiologic distributions for each symptom were largely the same across regions and population strata, despite HIV status and age being significantly associated with several infection diagnoses. Limitations of the review include the absence of studies in 16 of 48 SSA countries, substantial heterogeneity in study observations, and impeded assessment of this variability due to incomplete or inconsistent reporting across studies. CONCLUSIONS: In our study, syndrome aetiologies in SSA aligned with World Health Organization guidelines without strong evidence of geographic or demographic variation, supporting broad guideline applicability. Temporal changes underscore the importance of regular aetiologic re-assessment for effective syndromic management. PROSPERO NUMBER: CRD42022348045.

The impact of health inequity on spatial variation of COVID-19 transmission in England.

Considerable spatial heterogeneity has been observed in COVID-19 transmission across administrative areas of England throughout the pandemic. This study investigates what drives these differences. We constructed a probabilistic case count model for 306 administrative areas of England across 95 weeks, fit using a Bayesian evidence synthesis framework. We incorporate the impact of acquired immunity, of spatial exportation of cases, and 16 spatially-varying socio-economic, socio-demographic, health, and mobility variables. Model comparison assesses the relative contributions of these respective mechanisms. We find that spatially-varying and time-varying differences in week-to-week transmission were definitively associated with differences in: time spent at home, variant-of-concern proportion, and adult social care funding. However, model comparison demonstrates that the impact of these terms is negligible compared to the role of spatial exportation between administrative areas. While these results confirm the impact of some, but not all, static measures of spatially-varying inequity in England, our work corroborates the finding that observed differences in disease transmission during the pandemic were predominantly driven by underlying epidemiological factors rather than aggregated metrics of demography and health inequity between areas. Further work is required to assess how health inequity more broadly contributes to these epidemiological factors.

Temporal variations in international air travel: implications for modelling the spread of infectious diseases.

BACKGROUND: The international flight network creates multiple routes by which pathogens can quickly spread across the globe. In the early stages of infectious disease outbreaks, analyses using flight passenger data to identify countries at risk of importing the pathogen are common and can help inform disease control efforts. A challenge faced in this modelling is that the latest aviation statistics (referred to as contemporary data) are typically not immediately available. Therefore, flight patterns from a previous year are often used (referred to as historical data). We explored the suitability of historical data for predicting the spatial spread of emerging epidemics. METHODS: We analysed monthly flight passenger data from the International Air Transport Association to assess how baseline air travel patterns were affected in outbreaks of MERS, Zika, and SARS-CoV-2 over the past decade. We then used a stochastic discrete time SEIR metapopulation model to simulate global spread of different pathogens, comparing how epidemic dynamics differed in simulations based on historical and contemporary data. RESULTS: We observed local, short-term disruptions to air travel from South Korea and Brazil for the MERS and Zika outbreaks we studied, whereas global and longer-term flight disruption occurred during the SARS-CoV-2 pandemic.For outbreak events that were accompanied by local, small, and short-term changes in air travel, epidemic models using historical flight data gave similar projections of timing and locations of disease spread as when using contemporary flight data. However, historical data were less reliable to model the spread of an atypical outbreak such as SARS-CoV-2 in which there were durable and extensive levels of global travel disruption. CONCLUSIONS: The use of historical flight data as a proxy in epidemic models is an acceptable practice except in rare, large epidemics that lead to substantial disruptions to international travel.

SIR or MADAM? The impact of privilege on careers in epidemic modelling.

As we emerge from what may be the largest global public health crises of our lives, our community of epidemic modellers is naturally reflecting. What role can modelling play in supporting decision making during epidemics? How could we more effectively interact with policy makers? How should we design future disease surveillance systems? All crucial questions. But who is going to be addressing them in 10 years' time? With high burnout and poor attrition rates in academia, both magnified in our field by our unprecedented efforts during the pandemic, and with low wages coinciding with inflation at its highest for decades, how do we retain talent? This is a multifaceted challenge, that I argue is underpinned by privilege. In this perspective, I introduce the notion of privilege and highlight how various aspects of privilege (namely gender, ethnicity, sexual orientation, language and caring responsibilities) may affect the ability of individuals to access to and progress within academic modelling careers. I propose actions that members of the epidemic modelling research community may take to mitigate these issues and ensure we have a more diverse and equitable workforce going forward.

Corrigendum to "Modeling vaccination coverage during the 2022 central Ohio measles outbreak: a cross-sectional study" [The Lancet Regional Health-Americas 2023; 23: 100533].

[This corrects the article DOI: 10.1016/j.lana.2023.100533.].

Burden and Trends of Symptomatic Sexually Transmitted Infections in Malawi From 2000 to 2021: Comparative Analysis of Survey and Case Report Data.

BACKGROUND: In settings without etiologic testing for sexually transmitted infections (STIs), programs rely on STI symptom data to inform priorities. To evaluate whether self-reported STI symptoms in household surveys consistently represent the STI burden, we compared symptomatic infection rates between survey self-reporting and health facility case reporting in Malawi. METHODS: We analyzed self-reported symptoms and treatment seeking in the past year among sexually active adults from 4 Malawi Demographic and Health Surveys between 2000 and 2015. Bayesian mixed-effects models were used to estimate temporal trends, spatial variation, and sociodemographic determinants. Survey reporting was compared with health facility syndromic diagnoses between 2014 and 2021. RESULTS: In surveys, 11.0% (95% confidence interval, 10.7%-11.4%) of adults reported STI or STI-related symptoms in the last year, of whom 54.2% (52.8%-55.7%) sought treatment. In facilities, the mean annual symptomatic case diagnosis rate was 3.3%. Survey-reported treatment in the last year was 3.8% (95% credible interval, 2.3%-6.1%) for genital ulcer, 3.8% (2.0%-6.7%) for vaginal discharge, and 2.6% (1.2%-4.7%) for urethral discharge. Mean annual diagnosis rates at facilities were 0.5% for genital ulcer, 2.2% for vaginal discharge, and 2.0% for urethral discharge. Both data sources indicated a higher burden of symptoms among women, individuals older than 25 years, and those in Southern Malawi. CONCLUSIONS: Survey and facility case reports indicated similar spatial and demographic patterns of STI symptom burden and care seeking, but implied large differences in the magnitude and relative burden of symptoms, particularly genital ulcer, which could affect program priorities. Targeted etiologic surveillance would improve interpretation of these data to enable more comprehensive STI surveillance.

Marburg virus disease outbreaks, mathematical models, and disease parameters: a systematic review.

The 2023 Marburg virus disease outbreaks in Equatorial Guinea and Tanzania highlighted the importance of better understanding this lethal pathogen. We did a systematic review (PROSPERO CRD42023393345) of peer-reviewed articles reporting historical outbreaks, modelling studies, and epidemiological parameters focused on Marburg virus disease. We searched PubMed and Web of Science from database inception to March 31, 2023. Two reviewers evaluated all titles and abstracts with consensus-based decision making. To ensure agreement, 13 (31%) of 42 studies were double-extracted and a custom-designed quality assessment questionnaire was used for risk of bias assessment. We present detailed information on 478 reported cases and 385 deaths from Marburg virus disease. Analysis of historical outbreaks and seroprevalence estimates suggests the possibility of undetected Marburg virus disease outbreaks, asymptomatic transmission, or cross-reactivity with other pathogens, or a combination of these. Only one study presented a mathematical model of Marburg virus transmission. We estimate an unadjusted, pooled total random effect case fatality ratio of 61·9% (95% CI 38·8-80·6; I2=93%). We identify epidemiological parameters relating to transmission and natural history, for which there are few estimates. This systematic review and the accompanying database provide a comprehensive overview of Marburg virus disease epidemiology and identify key knowledge gaps, contributing crucial information for mathematical models to support future Marburg virus disease epidemic responses.

Retrospective evaluation of real-time estimates of global COVID-19 transmission trends and mortality forecasts.

Since 8th March 2020 up to the time of writing, we have been producing near real-time weekly estimates of SARS-CoV-2 transmissibility and forecasts of deaths due to COVID-19 for all countries with evidence of sustained transmission, shared online. We also developed a novel heuristic to combine weekly estimates of transmissibility to produce forecasts over a 4-week horizon. Here we present a retrospective evaluation of the forecasts produced between 8th March to 29th November 2020 for 81 countries. We evaluated the robustness of the forecasts produced in real-time using relative error, coverage probability, and comparisons with null models. During the 39-week period covered by this study, both the short- and medium-term forecasts captured well the epidemic trajectory across different waves of COVID-19 infections with small relative errors over the forecast horizon. The model was well calibrated with 56.3% and 45.6% of the observations lying in the 50% Credible Interval in 1-week and 4-week ahead forecasts respectively. The retrospective evaluation of our models shows that simple transmission models calibrated using routine disease surveillance data can reliably capture the epidemic trajectory in multiple countries. The medium-term forecasts can be used in conjunction with the short-term forecasts of COVID-19 mortality as a useful planning tool as countries continue to relax public health measures.

Bayesian reconstruction of SARS-CoV-2 transmissions highlights substantial proportion of negative serial intervals.

BACKGROUND: The serial interval is a key epidemiological measure that quantifies the time between the onset of symptoms in an infector-infectee pair. It indicates how quickly new generations of cases appear, thus informing on the speed of an epidemic. Estimating the serial interval requires to identify pairs of infectors and infectees. Yet, most studies fail to assess the direction of transmission between cases and assume that the order of infections - and thus transmissions - strictly follows the order of symptom onsets, thereby imposing serial intervals to be positive. Because of the long and highly variable incubation period of SARS-CoV-2, this may not always be true (i.e an infectee may show symptoms before their infector) and negative serial intervals may occur. This study aims to estimate the serial interval of different SARS-CoV-2 variants whilst accounting for negative serial intervals. METHODS: This analysis included 5 842 symptomatic individuals with confirmed SARS-CoV-2 infection amongst 2 579 households from September 2020 to August 2022 across England & Wales. We used a Bayesian framework to infer who infected whom by exploring all transmission trees compatible with the observed dates of symptoms, based on a wide range of incubation period and generation time distributions compatible with estimates reported in the literature. Serial intervals were derived from the reconstructed transmission pairs, stratified by variants. RESULTS: We estimated that 22% (95% credible interval (CrI) 8-32%) of serial interval values are negative across all VOC. The mean serial interval was shortest for Omicron BA5 (2.02 days, 1.26-2.84) and longest for Alpha (3.37 days, 2.52-4.04). CONCLUSIONS: This study highlights the large proportion of negative serial intervals across SARS-CoV-2 variants. Because the serial interval is widely used to estimate transmissibility and forecast cases, these results may have critical implications for epidemic control.

Estimating the epidemic reproduction number from temporally aggregated incidence data: A statistical modelling approach and software tool.

The time-varying reproduction number (Rt) is an important measure of epidemic transmissibility that directly informs policy decisions and the optimisation of control measures. EpiEstim is a widely used opensource software tool that uses case incidence and the serial interval (SI, time between symptoms in a case and their infector) to estimate Rt in real-time. The incidence and the SI distribution must be provided at the same temporal resolution, which can limit the applicability of EpiEstim and other similar methods, e.g. for contexts where the time window of incidence reporting is longer than the mean SI. In the EpiEstim R package, we implement an expectation-maximisation algorithm to reconstruct daily incidence from temporally aggregated data, from which Rt can then be estimated. We assess the validity of our method using an extensive simulation study and apply it to COVID-19 and influenza data. For all datasets, the influence of intra-weekly variability in reported data was mitigated by using aggregated weekly data. Rt estimated on weekly sliding windows using incidence reconstructed from weekly data was strongly correlated with estimates from the original daily data. The simulation study revealed that Rt was well estimated in all scenarios and regardless of the temporal aggregation of the data. In the presence of weekend effects, Rt estimates from reconstructed data were more successful at recovering the true value of Rt than those obtained from reported daily data. These results show that this novel method allows Rt to be successfully recovered from aggregated data using a simple approach with very few data requirements. Additionally, by removing administrative noise when daily incidence data are reconstructed, the accuracy of Rt estimates can be improved.

Scoring epidemiological forecasts on transformed scales.

Forecast evaluation is essential for the development of predictive epidemic models and can inform their use for public health decision-making. Common scores to evaluate epidemiological forecasts are the Continuous Ranked Probability Score (CRPS) and the Weighted Interval Score (WIS), which can be seen as measures of the absolute distance between the forecast distribution and the observation. However, applying these scores directly to predicted and observed incidence counts may not be the most appropriate due to the exponential nature of epidemic processes and the varying magnitudes of observed values across space and time. In this paper, we argue that transforming counts before applying scores such as the CRPS or WIS can effectively mitigate these difficulties and yield epidemiologically meaningful and easily interpretable results. Using the CRPS on log-transformed values as an example, we list three attractive properties: Firstly, it can be interpreted as a probabilistic version of a relative error. Secondly, it reflects how well models predicted the time-varying epidemic growth rate. And lastly, using arguments on variance-stabilizing transformations, it can be shown that under the assumption of a quadratic mean-variance relationship, the logarithmic transformation leads to expected CRPS values which are independent of the order of magnitude of the predicted quantity. Applying a transformation of log(x + 1) to data and forecasts from the European COVID-19 Forecast Hub, we find that it changes model rankings regardless of stratification by forecast date, location or target types. Situations in which models missed the beginning of upward swings are more strongly emphasised while failing to predict a downturn following a peak is less severely penalised when scoring transformed forecasts as opposed to untransformed ones. We conclude that appropriate transformations, of which the natural logarithm is only one particularly attractive option, should be considered when assessing the performance of different models in the context of infectious disease incidence.

Extending EpiEstim to estimate the transmission advantage of pathogen variants in real-time: SARS-CoV-2 as a case-study.

The evolution of SARS-CoV-2 has demonstrated that emerging variants can set back the global COVID-19 response. The ability to rapidly assess the threat of new variants is critical for timely optimisation of control strategies. We present a novel method to estimate the effective transmission advantage of a new variant compared to a reference variant combining information across multiple locations and over time. Through an extensive simulation study designed to mimic real-time epidemic contexts, we show that our method performs well across a range of scenarios and provide guidance on its optimal use and interpretation of results. We also provide an open-source software implementation of our method. The computational speed of our tool enables users to rapidly explore spatial and temporal variations in the estimated transmission advantage. We estimate that the SARS-CoV-2 Alpha variant is 1.46 (95% Credible Interval 1.44-1.47) and 1.29 (95% CrI 1.29-1.30) times more transmissible than the wild type, using data from England and France respectively. We further estimate that Delta is 1.77 (95% CrI 1.69-1.85) times more transmissible than Alpha (England data). Our approach can be used as an important first step towards quantifying the threat of emerging or co-circulating variants of infectious pathogens in real-time.

Modeling vaccination coverage during the 2022 central Ohio measles outbreak: a cross-sectional study.

Background: Of the eight large (>50 cases) US postelimination outbreaks, the first and last occurred in Ohio. Ohio's vaccination registry is incomplete. Community-level immunity gaps threaten more than two decades of measles elimination in the US. We developed a statistical model, VaxEstim, to rapidly estimate the early-phase vaccination coverage and immunity gap in the exposed population during the 2022 Central Ohio outbreak. Methods: We used reconstructed daily incidence (from publicly available data) and assumptions about the distribution of the serial interval, or the time between symptom onset in successive measles cases, to estimate the effective reproduction number (i.e., the average number of secondary infections caused by an infected individual in a partially immune population). We estimated early-phase measles vaccination coverage by comparing the effective reproduction number to the basic reproduction number (i.e., the average number of secondary infections caused by an infected individual in a fully susceptible population) while accounting for vaccine effectiveness. Finally, we estimated the early-phase immunity gap as the difference between the estimated critical vaccination threshold and vaccination coverage. Findings: VaxEstim estimated the early-phase vaccination coverage as 53% (95% credible interval, 21%-77%), the critical vaccination threshold as 93%, and the immunity gap as 42% (95% credible interval, 18%-74%). Interpretation: This study estimates a significant immunity gap in the exposed population during the early phase of the 2022 Central Ohio measles outbreak, suggesting a robust public health response is needed to identify the susceptible community and develop community-specific strategies to close the immunity gap. Funding: This work was supported in part by the National Institute of General Medical Sciences, National Institutes of Health; the UK Medical Research Council (MRC); the Foreign, Commonwealth and Development Office; the National Institute for Health Research (NIHR) Health Protection Research Unit in Modelling Methodology; Imperial College London, and the London School of Hygiene & Tropical Medicine, Community Jameel; the EDCTP2 programme, supported by the EU; and the Sergei Brin Foundation.

Epidemiological drivers of transmissibility and severity of SARS-CoV-2 in England.

As the SARS-CoV-2 pandemic progressed, distinct variants emerged and dominated in England. These variants, Wildtype, Alpha, Delta, and Omicron were characterized by variations in transmissibility and severity. We used a robust mathematical model and Bayesian inference framework to analyse epidemiological surveillance data from England. We quantified the impact of non-pharmaceutical interventions (NPIs), therapeutics, and vaccination on virus transmission and severity. Each successive variant had a higher intrinsic transmissibility. Omicron (BA.1) had the highest basic reproduction number at 8.3 (95% credible interval (CrI) 7.7-8.8). Varying levels of NPIs were crucial in controlling virus transmission until population immunity accumulated. Immune escape properties of Omicron decreased effective levels of immunity in the population by a third. Furthermore, in contrast to previous studies, we found Alpha had the highest basic infection fatality ratio (2.9%, 95% CrI 2.7-3.2), followed by Delta (2.2%, 95% CrI 2.0-2.4), Wildtype (1.2%, 95% CrI 1.1-1.2), and Omicron (0.7%, 95% CrI 0.6-0.8). Our findings highlight the importance of continued surveillance. Long-term strategies for monitoring and maintaining effective immunity against SARS-CoV-2 are critical to inform the role of NPIs to effectively manage future variants with potentially higher intrinsic transmissibility and severe outcomes.

Lessons from COVID-19 for rescalable data collection.

Novel data and analyses have had an important role in informing the public health response to the COVID-19 pandemic. Existing surveillance systems were scaled up, and in some instances new systems were developed to meet the challenges posed by the magnitude of the pandemic. We describe the routine and novel data that were used to address urgent public health questions during the pandemic, underscore the challenges in sustainability and equity in data generation, and highlight key lessons learnt for designing scalable data collection systems to support decision making during a public health crisis. As countries emerge from the acute phase of the pandemic, COVID-19 surveillance systems are being scaled down. However, SARS-CoV-2 resurgence remains a threat to global health security; therefore, a minimal cost-effective system needs to remain active that can be rapidly scaled up if necessary. We propose that a retrospective evaluation to identify the cost-benefit profile of the various data streams collected during the pandemic should be on the scientific research agenda.

Quantifying the effect of delaying the second COVID-19 vaccine dose in England: a mathematical modelling study.

BACKGROUND: The UK was the first country to start national COVID-19 vaccination programmes, initially administering doses 3 weeks apart. However, early evidence of high vaccine effectiveness after the first dose and the emergence of the SARS-CoV-2 alpha variant prompted the UK to extend the interval between doses to 12 weeks. In this study, we aimed to quantify the effect of delaying the second vaccine dose in England. METHODS: We used a previously described model of SARS-CoV-2 transmission, calibrated to COVID-19 surveillance data from England, including hospital admissions, hospital occupancy, seroprevalence data, and population-level PCR testing data, using a Bayesian evidence-synthesis framework. We modelled and compared the epidemic trajectory in the counterfactual scenario in which vaccine doses were administered 3 weeks apart against the real reported vaccine roll-out schedule of 12 weeks. We estimated and compared the resulting numbers of daily infections, hospital admissions, and deaths. In sensitivity analyses, we investigated scenarios spanning a range of vaccine effectiveness and waning assumptions. FINDINGS: In the period from Dec 8, 2020, to Sept 13, 2021, the number of individuals who received a first vaccine dose was higher under the 12-week strategy than the 3-week strategy. For this period, we estimated that delaying the interval between the first and second COVID-19 vaccine doses from 3 to 12 weeks averted a median (calculated as the median of the posterior sample) of 58 000 COVID-19 hospital admissions (291 000 cumulative hospitalisations [95% credible interval 275 000-319 000] under the 3-week strategy vs 233 000 [229 000-238 000] under the 12-week strategy) and 10 100 deaths (64 800 deaths [60 200-68 900] vs 54 700 [52 800-55 600]). Similarly, we estimated that the 3-week strategy would have resulted in more infections compared with the 12-week strategy. Across all sensitivity analyses the 3-week strategy resulted in a greater number of hospital admissions. In results by age group, the 12-week strategy led to more hospitalisations and deaths in older people in spring 2021, but fewer following the emergence of the delta variant during summer 2021. INTERPRETATION: England's delayed-second-dose vaccination strategy was informed by early real-world data on vaccine effectiveness in the context of limited vaccine supplies in a growing epidemic. Our study shows that rapidly providing partial (single-dose) vaccine-induced protection to a larger proportion of the population was successful in reducing the burden of COVID-19 hospitalisations and deaths overall. FUNDING: UK National Institute for Health Research; UK Medical Research Council; Community Jameel; Wellcome Trust; UK Foreign, Commonwealth and Development Office; Australian National Health and Medical Research Council; and EU.

Gaps in mobility data and implications for modelling epidemic spread: A scoping review and simulation study.

Reliable estimates of human mobility are important for understanding the spatial spread of infectious diseases and the effective targeting of control measures. However, when modelling infectious disease dynamics, data on human mobility at an appropriate temporal or spatial resolution are not always available, leading to the common use of model-derived mobility proxies. In this study we reviewed the different data sources and mobility models that have been used to characterise human movement in Africa. We then conducted a simulation study to better understand the implications of using human mobility proxies when predicting the spatial spread and dynamics of infectious diseases. We found major gaps in the availability of empirical measures of human mobility in Africa, leading to mobility proxies being used in place of data. Empirical data on subnational mobility were only available for 17/54 countries, and in most instances, these data characterised long-term movement patterns, which were unsuitable for modelling the spread of pathogens with short generation times (time between infection of a case and their infector). Results from our simulation study demonstrated that using mobility proxies can have a substantial impact on the predicted epidemic dynamics, with complex and non-intuitive biases. In particular, the predicted times and order of epidemic invasion, and the time of epidemic peak in different locations can be underestimated or overestimated, depending on the types of proxies used and the country of interest. Our work underscores the need for regularly updated empirical measures of population movement within and between countries to aid the prevention and control of infectious disease outbreaks. At the same time, there is a need to establish an evidence base to help understand which types of mobility data are most appropriate for describing the spread of emerging infectious diseases in different settings.