Best practices for estimating and reporting epidemiological delay distributions of infectious diseases.

Epidemiological delays are key quantities that inform public health policy and clinical practice. They are used as inputs for mathematical and statistical models, which in turn can guide control strategies. In recent work, we found that censoring, right truncation, and dynamical bias were rarely addressed correctly when estimating delays and that these biases were large enough to have knock-on impacts across a large number of use cases. Here, we formulate a checklist of best practices for estimating and reporting epidemiological delays. We also provide a flowchart to guide practitioners based on their data. Our examples are focused on the incubation period and serial interval due to their importance in outbreak response and modeling, but our recommendations are applicable to other delays. The recommendations, which are based on the literature and our experience estimating epidemiological delay distributions during outbreak responses, can help improve the robustness and utility of reported estimates and provide guidance for the evaluation of estimates for downstream use in transmission models or other analyses.

How effective is the BNT162b2 mRNA vaccine against SARS-CoV-2 transmission and infection? A national programme analysis in Monaco, July 2021 to September 2022.

BACKGROUND: We quantified SARS-CoV-2 dynamics in different community settings and the direct and indirect effect of the BNT162b2 mRNA vaccine in Monaco for different variants of concern (VOC). METHODS: Between July 2021 and September 2022, we prospectively investigated 20,443 contacts from 6320 index cases using data from the Monaco COVID-19 Public Health Programme. We calculated secondary attack rates (SARs) in households (n = 13,877), schools (n = 2508) and occupational (n = 6499) settings. We used binomial regression with a complementary log-log link function to measure adjusted hazard ratios (aHR) and vaccine effectiveness (aVE) for index cases to infect contacts and contacts to be infected in households. RESULTS: In households, the SAR was 55% (95% CI 54-57) and 50% (48-51) among unvaccinated and vaccinated contacts, respectively. The SAR was 32% (28-36) and 12% (10-13) in workplaces, and 7% (6-9) and 6% (3-10) in schools, among unvaccinated and vaccinated contacts respectively. In household, the aHR was lower in contacts than in index cases (aHR 0.68 [0.55-0.83] and 0.93 [0.74-1.1] for delta; aHR 0.73 [0.66-0.81] and 0.89 [0.80-0.99] for omicron BA.1&2, respectively). Vaccination had no significant effect on either direct or indirect aVE for omicron BA.4&5. The direct aVE in contacts was 32% (17, 45) and 27% (19, 34), and for index cases the indirect aVE was 7% (- 17, 26) and 11% (1, 20) for delta and omicron BA.1&2, respectively. The greatest aVE was in contacts with a previous SARS-CoV-2 infection and a single vaccine dose during the omicron BA.1&2 period (45% [27, 59]), while the lowest were found in contacts with either three vaccine doses (aVE - 24% [- 63, 6]) or one single dose and a previous SARS-CoV-2 infection (aVE - 36% [- 198, 38]) during the omicron BA.4&5 period. CONCLUSIONS: Protection conferred by the BNT162b2 mRNA vaccine against transmission and infection was low for delta and omicron BA.1&2, regardless of the number of vaccine doses and previous SARS-CoV-2 infection. There was no significant vaccine effect for omicron BA.4&5. Health authorities carrying out vaccination campaigns should bear in mind that the current generation of COVID-19 vaccines may not represent an effective tool in protecting individuals from either transmitting or acquiring SARS-CoV-2 infection.

Nowcasting the 2022 mpox outbreak in England.

In May 2022, a cluster of mpox cases were detected in the UK that could not be traced to recent travel history from an endemic region. Over the coming months, the outbreak grew, with over 3000 total cases reported in the UK, and similar outbreaks occurring worldwide. These outbreaks appeared linked to sexual contact networks between gay, bisexual and other men who have sex with men. Following the COVID-19 pandemic, local health systems were strained, and therefore effective surveillance for mpox was essential for managing public health policy. However, the mpox outbreak in the UK was characterised by substantial delays in the reporting of the symptom onset date and specimen collection date for confirmed positive cases. These delays led to substantial backfilling in the epidemic curve, making it challenging to interpret the epidemic trajectory in real-time. Many nowcasting models exist to tackle this challenge in epidemiological data, but these lacked sufficient flexibility. We have developed a nowcasting model using generalised additive models that makes novel use of individual-level patient data to correct the mpox epidemic curve in England. The aim of this model is to correct for backfilling in the epidemic curve and provide real-time characteristics of the state of the epidemic, including the real-time growth rate. This model benefited from close collaboration with individuals involved in collecting and processing the data, enabling temporal changes in the reporting structure to be built into the model, which improved the robustness of the nowcasts generated. The resulting model accurately captured the true shape of the epidemic curve in real time.

Identifying employee, workplace and population characteristics associated with COVID-19 outbreaks in the workplace: a population-based study.

OBJECTIVES: To identify risk factors that contribute to outbreaks of COVID-19 in the workplace and quantify their effect on outbreak risk. METHODS: We identified outbreaks of COVID-19 cases in the workplace and investigated the characteristics of the individuals, the workplaces, the areas they work and the mode of commute to work, through data linkages based on Middle Layer Super Output Areas in England between 20 June 2021 and 20 February 2022. We estimated population-level associations between potential risk factors and workplace outbreaks, adjusting for plausible confounders identified using a directed acyclic graph. RESULTS: For most industries, increased physical proximity in the workplace was associated with increased risk of COVID-19 outbreaks, while increased vaccination was associated with reduced risk. Employee demographic risk factors varied across industry, but for the majority of industries, a higher proportion of black/African/Caribbean ethnicities and living in deprived areas, was associated with increased outbreak risk. A higher proportion of employees in the 60-64 age group was associated with reduced outbreak risk. There were significant associations between gender, work commute modes and staff contract type with outbreak risk, but these were highly variable across industries. CONCLUSIONS: This study has used novel national data linkages to identify potential risk factors of workplace COVID-19 outbreaks, including possible protective effects of vaccination and increased physical distance at work. The same methodological approach can be applied to wider occupational and environmental health research.

EpiBeds: Data informed modelling of the COVID-19 hospital burden in England.

The first year of the COVID-19 pandemic put considerable strain on healthcare systems worldwide. In order to predict the effect of the local epidemic on hospital capacity in England, we used a variety of data streams to inform the construction and parameterisation of a hospital progression model, EpiBeds, which was coupled to a model of the generalised epidemic. In this model, individuals progress through different pathways (e.g. may recover, die, or progress to intensive care and recover or die) and data from a partially complete patient-pathway line-list was used to provide initial estimates of the mean duration that individuals spend in the different hospital compartments. We then fitted EpiBeds using complete data on hospital occupancy and hospital deaths, enabling estimation of the proportion of individuals that follow the different clinical pathways, the reproduction number of the generalised epidemic, and to make short-term predictions of hospital bed demand. The construction of EpiBeds makes it straightforward to adapt to different patient pathways and settings beyond England. As part of the UK response to the pandemic, EpiBeds provided weekly forecasts to the NHS for hospital bed occupancy and admissions in England, Wales, Scotland, and Northern Ireland at national and regional scales.

Novel methods for estimating the instantaneous and overall COVID-19 case fatality risk among care home residents in England.

The COVID-19 pandemic has had high mortality rates in the elderly and frail worldwide, particularly in care homes. This is driven by the difficulty of isolating care homes from the wider community, the large population sizes within care facilities (relative to typical households), and the age/frailty of the residents. To quantify the mortality risk posed by disease, the case fatality risk (CFR) is an important tool. This quantifies the proportion of cases that result in death. Throughout the pandemic, CFR amongst care home residents in England has been monitored closely. To estimate CFR, we apply both novel and existing methods to data on deaths in care homes, collected by Public Health England and the Care Quality Commission. We compare these different methods, evaluating their relative strengths and weaknesses. Using these methods, we estimate temporal trends in the instantaneous CFR (at both daily and weekly resolutions) and the overall CFR across the whole of England, and dis-aggregated at regional level. We also investigate how the CFR varies based on age and on the type of care required, dis-aggregating by whether care homes include nursing staff and by age of residents. This work has contributed to the summary of measures used for monitoring the UK epidemic.

Understanding the leading indicators of hospital admissions from COVID-19 across successive waves in the UK.

Following the end of universal testing in the UK, hospital admissions are a key measure of COVID-19 pandemic pressure. Understanding leading indicators of admissions at the National Health Service (NHS) Trust, regional and national geographies help health services plan for ongoing pressures. We explored the spatio-temporal relationships of leading indicators of hospitalisations across SARS-CoV-2 waves in England. This analysis includes an evaluation of internet search volumes from Google Trends, NHS triage calls and online queries, the NHS COVID-19 app, lateral flow devices (LFDs), and the ZOE app. Data sources were analysed for their feasibility as leading indicators using Granger causality, cross-correlation, and dynamic time warping at fine spatial scales. Google Trends and NHS triages consistently temporally led admissions in most locations, with lead times ranging from 5 to 20 days, whereas an inconsistent relationship was found for the ZOE app, NHS COVID-19 app, and LFD testing, which diminished with spatial resolution, showing cross-correlation of leads between -7 and 7 days. The results indicate that novel surveillance sources can be used effectively to understand the expected healthcare burden within hospital administrative areas though the temporal and spatial heterogeneity of these relationships is a key determinant of their operational public health utility.

Replacement dynamics and the pathogenesis of the Alpha, Delta and Omicron variants of SARS-CoV-2.

New SARS-CoV-2 variants causing COVID-19 are a major risk to public health worldwide due to the potential for phenotypic change and increases in pathogenicity, transmissibility and/or vaccine escape. Recognising signatures of new variants in terms of replacing growth and severity are key to informing the public health response. To assess this, we aimed to investigate key time periods in the course of infection, hospitalisation and death, by variant. We linked datasets on contact tracing (Contact Tracing Advisory Service), testing (the Second-Generation Surveillance System) and hospitalisation (the Admitted Patient Care dataset) for the entire length of contact tracing in the England - from March 2020 to March 2022. We modelled, for England, time delay distributions using a Bayesian doubly interval censored modelling approach for the SARS-CoV-2 variants Alpha, Delta, Delta Plus (AY.4.2), Omicron BA.1 and Omicron BA.2. This was conducted for the incubation period, the time from infection to hospitalisation and hospitalisation to death. We further modelled the growth of novel variant replacement using a generalised additive model with a negative binomial error structure and the relationship between incubation period length and the risk of a fatality using a Bernoulli generalised linear model with a logit link. The mean incubation periods for each variant were: Alpha 4.19 (95% credible interval (CrI) 4.13-4.26) days; Delta 3.87 (95% CrI 3.82-3.93) days; Delta Plus 3.92 (95% CrI 3.87-3.98) days; Omicron BA.1 3.67 (95% CrI 3.61-3.72) days and Omicron BA.2 3.48 (95% CrI 3.43-3.53) days. The mean time from infection to hospitalisation was for Alpha 11.31 (95% CrI 11.20-11.41) days, Delta 10.36 (95% CrI 10.26-10.45) days and Omicron BA.1 11.54 (95% CrI 11.38-11.70) days. The mean time from hospitalisation to death was, for Alpha 14.31 (95% CrI 14.00-14.62) days; Delta 12.81 (95% CrI 12.62-13.00) days and Omicron BA.2 16.02 (95% CrI 15.46-16.60) days. The 95th percentile of the incubation periods were: Alpha 11.19 (95% CrI 10.92-11.48) days; Delta 9.97 (95% CrI 9.73-10.21) days; Delta Plus 9.99 (95% CrI 9.78-10.24) days; Omicron BA.1 9.45 (95% CrI 9.23-9.67) days and Omicron BA.2 8.83 (95% CrI 8.62-9.05) days. Shorter incubation periods were associated with greater fatality risk when adjusted for age, sex, variant, vaccination status, vaccination manufacturer and time since last dose with an odds ratio of 0.83 (95% confidence interval 0.82-0.83) (P value < 0.05). Variants of SARS-CoV-2 that have replaced previously dominant variants have had shorter incubation periods. Conversely co-existing variants have had very similar and non-distinct incubation period distributions. Shorter incubation periods reflect generation time advantage, with a reduction in the time to the peak infectious period, and may be a significant factor in novel variant replacing growth. Shorter times for admission to hospital and death were associated with variant severity - the most severe variant, Delta, led to significantly earlier hospitalisation, and death. These measures are likely important for future risk assessment of new variants, and their potential impact on population health.

Effectiveness of the BNT162b2 (Pfizer-BioNTech) and the ChAdOx1 nCoV-19 (Oxford-AstraZeneca) vaccines for reducing susceptibility to infection with the Delta variant (B.1.617.2) of SARS-CoV-2.

BACKGROUND: From January to May 2021 the alpha variant (B.1.1.7) of SARS-CoV-2 was the most commonly detected variant in the UK. Following this, the Delta variant (B.1.617.2) then became the predominant variant. The UK COVID-19 vaccination programme started on 8th December 2020. Prior to the Delta variant, most vaccine effectiveness studies focused on the alpha variant. We therefore aimed to estimate the effectiveness of the BNT162b2 (Pfizer-BioNTech) and the ChAdOx1 nCoV-19 (Oxford-AstraZeneca) vaccines in preventing symptomatic and asymptomatic infection with respect to the Delta variant in a UK setting. METHODS: We used anonymised public health record data linked to infection data (PCR) using the Combined Intelligence for Population Health Action resource. We then constructed an SIR epidemic model to explain SARS-CoV-2 infection data across the Cheshire and Merseyside region of the UK. Vaccines were assumed to be effective after 21 days for 1 dose and 14 days for 2 doses. RESULTS: We determined that the effectiveness of the Oxford-AstraZeneca vaccine in reducing susceptibility to infection is 39% (95% credible interval [34, 43]) and 64% (95% credible interval [61, 67]) for a single dose and a double dose respectively. For the Pfizer-BioNTech vaccine, the effectiveness is 20% (95% credible interval [10, 28]) and 84% (95% credible interval [82, 86]) for a single-dose and a double dose respectively. CONCLUSION: Vaccine effectiveness for reducing susceptibility to SARS-CoV-2 infection shows noticeable improvement after receiving two doses of either vaccine. Findings also suggest that a full course of the Pfizer-BioNTech provides the optimal protection against infection with the Delta variant. This reinforces the need to complete the full course programme to maximise individual protection and reduce transmission.

Evolutionary graph theory derived from eco-evolutionary dynamics.

A biologically motivated individual-based framework for evolution in network-structured populations is developed that can accommodate eco-evolutionary dynamics. This framework is used to construct a network birth and death model. The evolutionary graph theory model, which considers evolutionary dynamics only, is derived as a special case, highlighting additional assumptions that diverge from real biological processes. This is achieved by introducing a negative ecological feedback loop that suppresses ecological dynamics by forcing births and deaths to be coupled. We also investigate how fitness, a measure of reproductive success used in evolutionary graph theory, is related to the life-history of individuals in terms of their birth and death rates. In simple networks, these ecologically motivated dynamics are used to provide new insight into the spread of adaptive mutations, both with and without clonal interference. For example, the star network, which is known to be an amplifier of selection in evolutionary graph theory, can inhibit the spread of adaptive mutations when individuals can die naturally.

Approximating Quasi-Stationary Behaviour in Network-Based SIS Dynamics.

Deterministic approximations to stochastic Susceptible-Infectious-Susceptible models typically predict a stable endemic steady-state when above threshold. This can be hard to relate to the underlying stochastic dynamics, which has no endemic steady-state but can exhibit approximately stable behaviour. Here, we relate the approximate models to the stochastic dynamics via the definition of the quasi-stationary distribution (QSD), which captures this approximately stable behaviour. We develop a system of ordinary differential equations that approximate the number of infected individuals in the QSD for arbitrary contact networks and parameter values. When the epidemic level is high, these QSD approximations coincide with the existing approximation methods. However, as we approach the epidemic threshold, the models deviate, with these models following the QSD and the existing methods approaching the all susceptible state. Through consistently approximating the QSD, the proposed methods provide a more robust link to the stochastic models.

Evolutionary bet-hedging in structured populations.

As ecosystems evolve, species can become extinct due to fluctuations in the environment. This leads to the evolutionary adaption known as bet-hedging, where species hedge against these fluctuations to reduce their likelihood of extinction. Environmental variation can be either within or between generations. Previous work has shown that selection for bet-hedging against within-generational variation should not occur in large populations. However, this work has been limited by assumptions of well-mixed populations, whereas real populations usually have some degree of structure. Using the framework of evolutionary graph theory, we show that through adding competition structure to the population, within-generational variation can have a significant impact on the evolutionary process for any population size. This complements research using subdivided populations, which suggests that within-generational variation is important when local population sizes are small. Together, these conclusions provide evidence to support observations by some ecologists that are contrary to the widely held view that only between-generational environmental variation has an impact on natural selection. This provides theoretical justification for further empirical study into this largely unexplored area.

Methods for approximating stochastic evolutionary dynamics on graphs.

Population structure can have a significant effect on evolution. For some systems with sufficient symmetry, analytic results can be derived within the mathematical framework of evolutionary graph theory which relate to the outcome of the evolutionary process. However, for more complicated heterogeneous structures, computationally intensive methods are required such as individual-based stochastic simulations. By adapting methods from statistical physics, including moment closure techniques, we first show how to derive existing homogenised pair approximation models and the exact neutral drift model. We then develop node-level approximations to stochastic evolutionary processes on arbitrarily complex structured populations represented by finite graphs, which can capture the different dynamics for individual nodes in the population. Using these approximations, we evaluate the fixation probability of invading mutants for given initial conditions, where the dynamics follow standard evolutionary processes such as the invasion process. Comparisons with the output of stochastic simulations reveal the effectiveness of our approximations in describing the stochastic processes and in predicting the probability of fixation of mutants on a wide range of graphs. Construction of these models facilitates a systematic analysis and is valuable for a greater understanding of the influence of population structure on evolutionary processes.

Deterministic epidemic models overestimate the basic reproduction number of observed outbreaks.

The basic reproduction number, R0, is a well-known quantifier of epidemic spread. However, a class of existing methods for estimating R0 from incidence data early in the epidemic can lead to an over-estimation of this quantity. In particular, when fitting deterministic models to estimate the rate of spread, we do not account for the stochastic nature of epidemics and that, given the same system, some outbreaks may lead to epidemics and some may not. Typically, an observed epidemic that we wish to control is a major outbreak. This amounts to implicit selection for major outbreaks which leads to the over-estimation problem. We formally characterised the split between major and minor outbreaks by using Otsu's method which provides us with a working definition. We show that by conditioning a 'deterministic' model on major outbreaks, we can more reliably estimate the basic reproduction number from an observed epidemic trajectory.

Modelling multiplex testing for outbreak control.

During the SARS-CoV-2 pandemic, polymerase chain reaction (PCR) and lateral flow device (LFD) tests were frequently deployed to detect the presence of SARS-CoV-2. Many of these tests were singleplex, and only tested for the presence of a single pathogen. Multiplex tests can test for the presence of several pathogens using only a single swab, which can allow for: surveillance of more pathogens, targeting of antiviral interventions, a reduced burden of testing, and lower costs. Test sensitivity, however, particularly in LFD tests, is highly conditional on the viral concentration dynamics of individuals. To inform the use of multiplex testing in outbreak detection it is therefore necessary to investigate the interactions between outbreak detection strategies and the differing viral concentration trajectories of key pathogens. Viral concentration trajectories are estimated for SARS-CoV-2 and Influenza A/B. Testing strategies for the first five symptomatic cases in an outbreak are then simulated and used to evaluate key performance indicators. Strategies that use a combination of multiplex LFD and PCR tests achieve; high levels of detection, detect outbreaks rapidly, and have the lowest burden of testing across multiple pathogens. Influenza B was estimated to have lower rates of detection due to its modelled viral concentration dynamics.

The real-time infection hospitalisation and fatality risk across the COVID-19 pandemic in England.

The COVID-19 pandemic led to 231,841 deaths and 940,243 hospitalisations in England, by the end of March 2023. This paper calculates the real-time infection hospitalisation risk (IHR) and infection fatality risk (IFR) using the Office for National Statistics Coronavirus Infection Survey (ONS CIS) and the Real-time Assessment of Community Transmission Survey between November 2020 to March 2023. The IHR and the IFR in England peaked in January 2021 at 3.39% (95% Credible Intervals (CrI): 2.79, 3.97) and 0.97% (95% CrI: 0.62, 1.36), respectively. After this time, there was a rapid decline in the severity from infection, with the lowest estimated IHR of 0.32% (95% CrI: 0.27, 0.39) in December 2022 and IFR of 0.06% (95% CrI: 0.04, 0.08) in April 2022. We found infection severity to vary more markedly between regions early in the pandemic however, the absolute heterogeneity has since reduced. The risk from infection of SARS-CoV-2 has changed substantially throughout the COVID-19 pandemic with a decline of 86.03% (80.86, 89.35) and 89.67% (80.18, 93.93) in the IHR and IFR, respectively, since early 2021. From April 2022 until March 2023, the end of the ONS CIS study, we found fluctuating patterns in the severity of infection with the resumption of more normative mixing, resurgent epidemic waves, patterns of waning immunity, and emerging variants that have shown signs of convergent evolution.

Understanding the infection severity and epidemiological characteristics of mpox in the UK.

In May 2022, individuals infected with the monkeypox virus were detected in the UK without clear travel links to endemic areas. Understanding the clinical characteristics and infection severity of mpox is necessary for effective public health policy. The study period of this paper, from the 1st June 2022 to 30th September 2022, included 3,375 individuals that tested positive for the monkeypox virus. The posterior mean times from infection to hospital admission and length of hospital stay were 14.89 days (95% Credible Intervals (CrI): 13.60, 16.32) and 7.07 days (95% CrI: 6.07, 8.23), respectively. We estimated the modelled Infection Hospitalisation Risk to be 4.13% (95% CrI: 3.04, 5.02), compared to the overall sample Case Hospitalisation Risk (CHR) of 5.10% (95% CrI: 4.38, 5.86). The overall sample CHR was estimated to be 17.86% (95% CrI: 6.06, 33.11) for females and 4.99% (95% CrI: 4.27, 5.75) for males. A notable difference was observed between the CHRs that were estimated for each sex, which may be indicative of increased infection severity in females or a considerably lower infection ascertainment rate. It was estimated that 74.65% (95% CrI: 55.78, 86.85) of infections with the monkeypox virus in the UK were captured over the outbreak.

Forecasting influenza hospital admissions within English sub-regions using hierarchical generalised additive models.

BACKGROUND: Seasonal influenza places a substantial burden annually on healthcare services. Policies during the COVID-19 pandemic limited the transmission of seasonal influenza, making the timing and magnitude of a potential resurgence difficult to ascertain and its impact important to forecast. METHODS: We have developed a hierarchical generalised additive model (GAM) for the short-term forecasting of hospital admissions with a positive test for the influenza virus sub-regionally across England. The model incorporates a multi-level structure of spatio-temporal splines, weekly cycles in admissions, and spatial correlation. Using multiple performance metrics including interval score, coverage, bias, and median absolute error, the predictive performance is evaluated for the 2022-2023 seasonal wave. Performance is measured against autoregressive integrated moving average (ARIMA) and Prophet time series models. RESULTS: Across the epidemic phases the hierarchical GAM shows improved performance, at all geographic scales relative to the ARIMA and Prophet models. Temporally, the hierarchical GAM has overall an improved performance at 7 and 14 day time horizons. The performance of the GAM is most sensitive to the flexibility of the smoothing function that measures the national epidemic trend. CONCLUSIONS: This study introduces an approach to short-term forecasting of hospital admissions for the influenza virus using hierarchical, spatial, and temporal components. The methodology was designed for the real time forecasting of epidemics. This modelling framework was used across the 2022-2023 winter for healthcare operational planning by the UK Health Security Agency and the National Health Service in England.

Seasonal influenza causes a burden for hospitals and therefore it is useful to be able to accurately predict how many patients might be admitted with the disease. We attempted to predict influenza admissions up to 14 days in the future by creating a computational model that incorporates how the disease is reported and how it spreads. We evaluated our optimised model on data acquired during the winter of 2022-2023 data in England and compared it with previously developed models. Our model was better at modelling how influenza spreads and predicting future hospital admissions than the models we compared it to. Improving how influenza admissions are forecast can enable hospitals to prepare better for increased admissions, enabling improved treatment and reduced death for all patients in hospital over winter.

The rapid replacement of the SARS-CoV-2 Delta variant by Omicron (B.1.1.529) in England.

The emergence of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) B.1.1.529 (Omicron) variant caused international concern due to its rapid spread in Southern Africa. It was unknown whether this variant would replace or coexist with (either transiently or long term) the then-dominant Delta variant on its introduction to England. We developed a set of hierarchical logistic growth models to describe changes in the frequency of S gene target failure (SGTF) PCR tests, a proxy for Omicron. The doubling time of SGTF cases peaked at 1.56 days (95% CI: 1.49 to 1.63) on 5 December, whereas triple-positive cases were halving every 5.82 days (95% CI: 5.11 to 6.67) going into Christmas 2021. We were unable to characterize the replacement of Delta by Omicron with a single rate. The replacement rate decreased by 53.56% (95% CrI: 45.38 to 61.01) between 14 and 15 December, meaning the competitive advantage of Omicron approximately halved. Preceding the changepoint, Omicron was replacing Delta 16.24% (95% CrI: 9.72 to 23.41) faster in those with two or more vaccine doses, indicative of vaccine escape being a substantial component of competitive advantage. Despite the slowdown, Delta was almost entirely replaced in England within a month of the first sequenced domestic case. The synchrony of changepoints across regions at various stages of Omicron epidemics suggests that the growth rate advantage was not attenuated because of biological mechanisms related to strain competition. The step change in replacement could have resulted from behavioral changes, potentially elicited by public health messaging or policies, that differentially affected Omicron.

Transmission dynamics of monkeypox in the United Kingdom: contact tracing study.

OBJECTIVE: To analyse the transmission dynamics of the monkeypox outbreak in the UK, declared a Public Health Emergency of International Concern in July 2022. DESIGN: Contact tracing study, linking data on case-contact pairs and on probable exposure dates. SETTING: Case questionnaires from the UK Health Security Agency (UKHSA), United Kingdom. PARTICIPANTS: 2746 people with polymerase chain reaction confirmed monkeypox virus in the UK between 6 May and 1 August 2022. MAIN OUTCOME MEASURES: The incubation period and serial interval of a monkeypox infection using two bayesian time delay models-one corrected for interval censoring (ICC-interval censoring corrected) and one corrected for interval censoring, right truncation, and epidemic phase bias (ICRTC-interval censoring right truncation corrected). Growth rates of cases by reporting date, when monkeypox virus was confirmed and reported to UKHSA, were estimated using generalised additive models. RESULTS: The mean age of participants was 37.8 years and 95% reported being gay, bisexual, and other men who have sex with men (1160 out of 1213 reporting). The mean incubation period was estimated to be 7.6 days (95% credible interval 6.5 to 9.9) using the ICC model and 7.8 days (6.6 to 9.2) using the ICRTC model. The estimated mean serial interval was 8.0 days (95% credible interval 6.5 to 9.8) using the ICC model and 9.5 days (7.4 to 12.3) using the ICRTC model. Although the mean serial interval was longer than the incubation period for both models, short serial intervals were more common than short incubation periods, with the 25th centile and the median of the serial interval shorter than the incubation period. For the ICC and ICRTC models, the corresponding estimates ranged from 1.8 days (95% credible interval 1.5 to 1.8) to 1.6 days (1.4 to 1.6) shorter at the 25th centile and 1.6 days (1.5 to 1.7) to 0.8 days (0.3 to 1.2) shorter at the median. 10 out of 13 linked patients had documented pre-symptomatic transmission. Doubling times of cases declined from 9.07 days (95% confidence interval 12.63 to 7.08) on the 6 May, when the first case of monkeypox was reported in the UK, to a halving time of 29 days (95% confidence interval 38.02 to 23.44) on 1 August. CONCLUSIONS: Analysis of the instantaneous growth rate of monkeypox incidence indicates that the epidemic peaked in the UK as of 9 July and then started to decline. Short serial intervals were more common than short incubation periods suggesting considerable pre-symptomatic transmission, which was validated through linked patient level records. For patients who could be linked through personally identifiable data, four days was the maximum time that transmission was detected before symptoms manifested. An isolation period of 16 to 23 days would be required to detect 95% of people with a potential infection. The 95th centile of the serial interval was between 23 and 41 days, suggesting long infectious periods.