Modelling the influence of naturally acquired immunity from subclinical infection on outbreak dynamics and persistence of rabies in domestic dogs

A number of mathematical models have been developed for canine rabies to explore dynamics and inform control strategies. A common assumption of these models is that naturally acquired immunity plays no role in rabies dynamics. However, empirical studies have detected rabies-specific antibodies in healthy, unvaccinated domestic dogs, potentially due to immunizing, non-lethal exposure. We developed a stochastic model for canine rabies, parameterised for Laikipia County, Kenya, to explore the implications of different scenarios for naturally acquired immunity to rabies in domestic dogs. Simulating these scenarios using a non-spatial model indicated that low levels of immunity can act to limit rabies incidence and prevent depletion of the domestic dog population, increasing the probability of disease persistence. However, incorporating spatial structure and human response to high rabies incidence allowed the virus to persist in the absence of immunity. While low levels of immunity therefore had limited influence under a more realistic approximation of rabies dynamics, high rates of exposure leading to immunizing non-lethal exposure were required to produce population-level seroprevalences comparable with those reported in empirical studies. False positives and/or spatial variation may contribute to high empirical seroprevalences. However, if high seroprevalences are related to high exposure rates, these findings support the need for high vaccination coverage to effectively control this disease.

A Comparative Analysis of Statistical Methods to Estimate the Reproduction Number in Emerging Epidemics, With Implications for the Current Coronavirus Disease 2019 (COVID-19) Pandemic

BACKGROUND

As the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic continues its rapid global spread, quantification of local transmission patterns has been, and will continue to be, critical for guiding the pandemic response. Understanding the accuracy and limitations of statistical methods to estimate the basic reproduction number, R0, in the context of emerging epidemics is therefore vital to ensure appropriate interpretation of results and the subsequent implications for control efforts.

METHODS

Using simulated epidemic data, we assess the performance of 7 commonly used statistical methods to estimate R0 as they would be applied in a real-time outbreak analysis scenario: fitting to an increasing number of data points over time and with varying levels of random noise in the data. Method comparison was also conducted on empirical outbreak data, using Zika surveillance data from the 2015-2016 epidemic in Latin America and the Caribbean.

RESULTS

We find that most methods considered here frequently overestimate R0 in the early stages of epidemic growth on simulated data, the magnitude of which decreases when fitted to an increasing number of time points. This trend of decreasing bias over time can easily lead to incorrect conclusions about the course of the epidemic or the need for control efforts.

CONCLUSIONS

We show that true changes in pathogen transmissibility can be difficult to disentangle from changes in methodological accuracy and precision in the early stages of epidemic growth, particularly for data with significant over-dispersion. As localized epidemics of SARS-CoV-2 take hold around the globe, awareness of this trend will be important for appropriately cautious interpretation of results and subsequent guidance for control efforts.

Descriptive analysis of surveillance data for Zika virus disease and Zika virus-associated neurological complications in Colombia, 2015-2017

Zika virus (ZIKV) is a mosquito-borne pathogen that recently caused a major epidemic in the Americas. Although the majority of ZIKV infections are asymptomatic, the virus has been associated with birth defects in fetuses and newborns of infected mothers as well as neurological complications in adults. We performed a descriptive analysis on approximately 106,000 suspected and laboratory-confirmed cases of Zika virus disease (ZVD) that were reported during the 2015-2017 epidemic in Colombia. We also analyzed a dataset containing patients with neurological complications and recent febrile illness compatible with ZVD. Females had higher cumulative incidence of ZVD than males. Compared to the general population, cases were more likely to be reported in young adults (20 to 39 years of age). We estimated the cumulative incidence of ZVD in pregnant females at 3,120 reported cases per 100,000 population (95% CI: 3,077-3,164), which was considerably higher than the incidence in both males and non-pregnant females. ZVD cases were reported in all 32 departments. Four-hundred and eighteen patients suffered from ZIKV-associated neurological complications, of which 85% were diagnosed with Guillain-Barré syndrome. The median age of ZIKV cases with neurological complications was 12 years older than that of ZVD cases. ZIKV-associated neurological complications increased with age, and the highest incidence was reported among individuals aged 75 and older. Even though neurological complications and deaths due to ZIKV were rare in this epidemic, better risk communication is needed for people living in or traveling to ZIKV-affected areas.

Resurgence of SARS-CoV-2: Detection by community viral surveillance

Surveillance of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic has mainly relied on case reporting, which is biased by health service performance, test availability, and test-seeking behaviors. We report a community-wide national representative surveillance program in England based on self-administered swab results from ~594,000 individuals tested for SARS-CoV-2, regardless of symptoms, between May and the beginning of September 2020. The epidemic declined between May and July 2020 but then increased gradually from mid-August, accelerating into early September 2020 at the start of the second wave. When compared with cases detected through routine surveillance, we report here a longer period of decline and a younger age distribution. Representative community sampling for SARS-CoV-2 can substantially improve situational awareness and feed into the public health response even at low prevalence.

The J-IDEA Pandemic Planner: A Framework for Implementing Hospital Provision Interventions During the COVID-19 Pandemic

BACKGROUND

Planning for extreme surges in demand for hospital care of patients requiring urgent life-saving treatment for coronavirus disease 2019 (COVID-19), while retaining capacity for other emergency conditions, is one of the most challenging tasks faced by health care providers and policymakers during the pandemic. Health systems must be well-prepared to cope with large and sudden changes in demand by implementing interventions to ensure adequate access to care. We developed the first planning tool for the COVID-19 pandemic to account for how hospital provision interventions (such as cancelling elective surgery, setting up field hospitals, or hiring retired staff) will affect the capacity of hospitals to provide life-saving care.

METHODS

We conducted a review of interventions implemented or considered in 12 European countries in March to April 2020, an evaluation of their impact on capacity, and a review of key parameters in the care of COVID-19 patients. This information was used to develop a planner capable of estimating the impact of specific interventions on doctors, nurses, beds, and respiratory support equipment. We applied this to a scenario-based case study of 1 intervention, the set-up of field hospitals in England, under varying levels of COVID-19 patients.

RESULTS

The Abdul Latif Jameel Institute for Disease and Emergency Analytics pandemic planner is a hospital planning tool that allows hospital administrators, policymakers, and other decision-makers to calculate the amount of capacity in terms of beds, staff, and crucial medical equipment obtained by implementing the interventions. Flexible assumptions on baseline capacity, the number of hospitalizations, staff-to-beds ratios, and staff absences due to COVID-19 make the planner adaptable to multiple settings. The results of the case study show that while field hospitals alleviate the burden on the number of beds available, this intervention is futile unless the deficit of critical care nurses is addressed first.

DISCUSSION

The tool supports decision-makers in delivering a fast and effective response to the pandemic. The unique contribution of the planner is that it allows users to compare the impact of interventions that change some or all inputs.

Prevalence of antibody positivity to SARS-CoV-2 following the first peak of infection in England: Serial cross-sectional studies of 365,000 adults

BACKGROUND

The time-concentrated nature of the first wave of the COVID-19 epidemic in England in March and April 2020 provides a natural experiment to measure changes in antibody positivity at the population level before onset of the second wave and initiation of the vaccination programme.

METHODS

Three cross-sectional national surveys with non-overlapping random samples of the population in England undertaken between late June and September 2020 (REACT-2 study). 365,104 adults completed questionnaires and self-administered lateral flow immunoassay (LFIA) tests for IgG against SARS-CoV-2.

FINDINGS

Overall, 17,576 people had detectable antibodies, a prevalence of 4.9% (95% confidence intervals 4.9, 5.0) when adjusted for test characteristics and weighted to the adult population of England. The prevalence declined from 6.0% (5.8, 6.1), to 4.8% (4.7, 5.0) and 4.4% (4.3, 4.5), over the three rounds of the study a difference of -26.5% (-29.0, -23.8). The highest prevalence and smallest overall decline in positivity was in the youngest age group (18-24 years) at -14.9% (-21.6, -8.1), and lowest prevalence and largest decline in the oldest group (>74 years) at -39.0% (-50.8, -27.2). The decline from June to September 2020 was largest in those who did not report a history of COVID-19 at -64.0% (-75.6, -52.3), compared to -22.3% (-27.0, -17.7) in those with SARS-CoV-2 infection confirmed on PCR.

INTERPRETATION

A large proportion of the population remained susceptible to SARS-CoV-2 infection in England based on naturally acquired immunity from the first wave. Widespread vaccination is needed to confer immunity and control the epidemic at population level.

FUNDING

This work was funded by the Department of Health and Social Care in England.

Interventions targeting non-symptomatic cases can be important to prevent local outbreaks: SARS-CoV-2 as a case study

During infectious disease epidemics, an important question is whether cases travelling to new locations will trigger local outbreaks. The risk of this occurring depends on the transmissibility of the pathogen, the susceptibility of the host population and, crucially, the effectiveness of surveillance in detecting cases and preventing onward spread. For many pathogens, transmission from pre-symptomatic and/or asymptomatic (together referred to as non-symptomatic) infectious hosts can occur, making effective surveillance challenging. Here, by using SARS-CoV-2 as a case study, we show how the risk of local outbreaks can be assessed when non-symptomatic transmission can occur. We construct a branching process model that includes non-symptomatic transmission and explore the effects of interventions targeting non-symptomatic or symptomatic hosts when surveillance resources are limited. We consider whether the greatest reductions in local outbreak risks are achieved by increasing surveillance and control targeting non-symptomatic or symptomatic cases, or a combination of both. We find that seeking to increase surveillance of symptomatic hosts alone is typically not the optimal strategy for reducing outbreak risks. Adopting a strategy that combines an enhancement of surveillance of symptomatic cases with efforts to find and isolate non-symptomatic infected hosts leads to the largest reduction in the probability that imported cases will initiate a local outbreak.

Genetic evidence for the association between COVID-19 epidemic severity and timing of non-pharmaceutical interventions

Unprecedented public health interventions including travel restrictions and national lockdowns have been implemented to stem the COVID-19 epidemic, but the effectiveness of non-pharmaceutical interventions is still debated. We carried out a phylogenetic analysis of more than 29,000 publicly available whole genome SARS-CoV-2 sequences from 57 locations to estimate the time that the epidemic originated in different places. These estimates were examined in relation to the dates of the most stringent interventions in each location as well as to the number of cumulative COVID-19 deaths and phylodynamic estimates of epidemic size. Here we report that the time elapsed between epidemic origin and maximum intervention is associated with different measures of epidemic severity and explains 11% of the variance in reported deaths one month after the most stringent intervention. Locations where strong non-pharmaceutical interventions were implemented earlier experienced much less severe COVID-19 morbidity and mortality during the period of study.

A Quantitative Framework for Defining the End of an Infectious Disease Outbreak: Application to Ebola Virus Disease

The end-of-outbreak declaration is an important step in controlling infectious disease outbreaks. Objective estimation of the confidence level that an outbreak is over is important to reduce the risk of postdeclaration flare-ups. We developed a simulation-based model with which to quantify that confidence and tested it on simulated Ebola virus disease data. We found that these confidence estimates were most sensitive to the instantaneous reproduction number, the reporting rate, and the time between the symptom onset and death or recovery of the last detected case. For Ebola virus disease, our results suggested that the current World Health Organization criterion of 42 days since the recovery or death of the last detected case is too short and too sensitive to underreporting. Therefore, we suggest a shift to a preliminary end-of-outbreak declaration after 63 days from the symptom onset day of the last detected case. This preliminary declaration should still be followed by 90 days of enhanced surveillance to capture potential flare-ups of cases, after which the official end of the outbreak can be declared. This sequence corresponds to more than 95% confidence that an outbreak is over in most of the scenarios examined. Our framework is generic and therefore could be adapted to estimate end-of-outbreak confidence for other infectious diseases.

Better educational signage could reduce disturbance of resting dolphins

Spinner dolphins on Hawai‘i Island’s west coast (Stenella longirostris longirostris) rest by day in protected bays that are increasingly popular for recreation. Because more frequent interactions of people with these dolphins is likely to reduce rest for dolphins and to explain recent decline in dolphin abundance, the National Oceanic and Atmospheric Administration (NOAA) proposed stricter rules regarding interactions with spinner dolphins near the main Hawaiian Islands and plans to increase enforcement. Simultaneous investment in public education about both interaction rules and their biological rationale has been and is likely to be relatively low. To test the hypothesis that more educational signage will reduce human-generated disturbance of dolphins, a paper questionnaire was distributed to 351 land-based, mostly unguided visitors at three dolphin resting bays on Hawai‘i Island’s west coast. Responses indicated that visitors wanted to see dolphins, were ignorant of interaction rules, were likely to read signs explaining rules and their biological rationales, and were likely to follow known rules. Therefore, investment in effective educational signage at dolphin resting bays is recommended as one way to support conservation of spinner dolphins on Hawai‘i Island’s west coast and similar sites in the Hawaiian archipelago.

Reduction in mobility and COVID-19 transmission

In response to the COVID-19 pandemic, countries have sought to control SARS-CoV-2 transmission by restricting population movement through social distancing interventions, thus reducing the number of contacts. Mobility data represent an important proxy measure of social distancing, and here, we characterise the relationship between transmission and mobility for 52 countries around the world. Transmission significantly decreased with the initial reduction in mobility in 73% of the countries analysed, but we found evidence of decoupling of transmission and mobility following the relaxation of strict control measures for 80% of countries. For the majority of countries, mobility explained a substantial proportion of the variation in transmissibility (median adjusted R-squared: 48%, interquartile range - IQR - across countries [27-77%]). Where a change in the relationship occurred, predictive ability decreased after the relaxation; from a median adjusted R-squared of 74% (IQR across countries [49-91%]) pre-relaxation, to a median adjusted R-squared of 30% (IQR across countries [12-48%]) post-relaxation. In countries with a clear relationship between mobility and transmission both before and after strict control measures were relaxed, mobility was associated with lower transmission rates after control measures were relaxed indicating that the beneficial effects of ongoing social distancing behaviours were substantial.

SARS-CoV-2 antibody prevalence in England following the first peak of the pandemic

England has experienced a large outbreak of SARS-CoV-2, disproportionately affecting people from disadvantaged and ethnic minority communities. It is unclear how much of this excess is due to differences in exposure associated with structural inequalities. Here, we report from the REal-time Assessment of Community Transmission-2 (REACT-2) national study of over 100,000 people. After adjusting for test characteristics and re-weighting to the population, overall antibody prevalence is 6.0% (95% CI: 5.8-6.1). An estimated 3.4 million people had developed antibodies to SARS-CoV-2 by mid-July 2020. Prevalence is two- to three-fold higher among health and care workers compared with non-essential workers, and in people of Black or South Asian than white ethnicity, while age- and sex-specific infection fatality ratios are similar across ethnicities. Our results indicate that higher hospitalisation and mortality from COVID-19 in minority ethnic groups may reflect higher rates of infection rather than differential experience of disease or care.

Early Transmission Dynamics, Spread, and Genomic Characterization of SARS-CoV-2 in Panama

We report an epidemiologic analysis of 4,210 cases of infection with severe acute respiratory syndrome coronavirus 2 and genetic analysis of 313 new near-complete virus genomes in Panama during March 9-April 16, 2020. Although containment measures reduced R0 and Rt, they did not interrupt virus spread in the country.

Database of epidemic trends and control measures during the first wave of COVID-19 in mainland China

OBJECTIVES

In this data collation study, we aimed to provide a comprehensive database describing the epidemic trends and responses during the first wave of coronavirus disease 2019 (COVID-19) throughout the main provinces in China.

METHODS

From mid-January to March 2020, we extracted publicly available data regarding the spread and control of COVID-19 from 31 provincial health authorities and major media outlets in mainland China. Based on these data, we conducted descriptive analyses of the epidemic in the six most-affected provinces.

RESULTS

School closures, travel restrictions, community-level lockdown, and contact tracing were introduced concurrently around late January but subsequent epidemic trends differed among provinces. Compared with Hubei, the other five most-affected provinces reported a lower crude case fatality ratio and proportion of critical and severe hospitalised cases. From March 2020, as the local transmission of COVID-19 declined, switching the focus of measures to the testing and quarantine of inbound travellers may have helped to sustain the control of the epidemic.

CONCLUSIONS

Aggregated indicators of case notifications and severity distributions are essential for monitoring an epidemic. A publicly available database containing these indicators and information regarding control measures is a useful resource for further research and policy planning in response to the COVID-19 epidemic.

State-level tracking of COVID-19 in the United States

As of 1st June 2020, the US Centres for Disease Control and Prevention reported 104,232 confirmed or probable COVID-19-related deaths in the US. This was more than twice the number of deaths reported in the next most severely impacted country. We jointly model the US epidemic at the state-level, using publicly available death data within a Bayesian hierarchical semi-mechanistic framework. For each state, we estimate the number of individuals that have been infected, the number of individuals that are currently infectious and the time-varying reproduction number (the average number of secondary infections caused by an infected person). We use changes in mobility to capture the impact that non-pharmaceutical interventions and other behaviour changes have on the rate of transmission of SARS-CoV-2. We estimate that R_t was only below one in 23 states on 1st June. We also estimate that 3.7% [3.4%-4.0%] of the total population of the US had been infected, with wide variation between states, and approximately 0.01% of the population was infectious. We demonstrate good 3 week model forecasts of deaths with low error and good coverage of our credible intervals.

State-level tracking of COVID-19 in the United States

As of 1st June 2020, the US Centres for Disease Control and Prevention reported 104,232 confirmed or probable COVID-19-related deaths in the US. This was more than twice the number of deaths reported in the next most severely impacted country. We jointly model the US epidemic at the state-level, using publicly available death data within a Bayesian hierarchical semi-mechanistic framework. For each state, we estimate the number of individuals that have been infected, the number of individuals that are currently infectious and the time-varying reproduction number (the average number of secondary infections caused by an infected person). We use changes in mobility to capture the impact that non-pharmaceutical interventions and other behaviour changes have on the rate of transmission of SARS-CoV-2. We estimate that Rt was only below one in 23 states on 1st June. We also estimate that 3.7% [3.4%-4.0%] of the total population of the US had been infected, with wide variation between states, and approximately 0.01% of the population was infectious. We demonstrate good 3 week model forecasts of deaths with low error and good coverage of our credible intervals.

State-level tracking of COVID-19 in the United States

As of 1st June 2020, the US Centres for Disease Control and Prevention reported 104,232 confirmed or probable COVID-19-related deaths in the US. This was more than twice the number of deaths reported in the next most severely impacted country. We jointly model the US epidemic at the state-level, using publicly available death data within a Bayesian hierarchical semi-mechanistic framework. For each state, we estimate the number of individuals that have been infected, the number of individuals that are currently infectious and the time-varying reproduction number (the average number of secondary infections caused by an infected person). We use changes in mobility to capture the impact that non-pharmaceutical interventions and other behaviour changes have on the rate of transmission of SARS-CoV-2. We estimate that R_t was only below one in 23 states on 1st June. We also estimate that 3.7% [3.4%-4.0%] of the total population of the US had been infected, with wide variation between states, and approximately 0.01% of the population was infectious. We demonstrate good 3 week model forecasts of deaths with low error and good coverage of our credible intervals.

Endemic and Epidemic Human Alphavirus Infections in Eastern Panama: An Analysis of Population-Based Cross-Sectional Surveys

Madariaga virus (MADV) has recently been associated with severe human disease in Panama, where the closely related Venezuelan equine encephalitis virus (VEEV) also circulates. In June 2017, a fatal MADV infection was confirmed in a community of Darien Province. We conducted a cross-sectional outbreak investigation with human and mosquito collections in July 2017, where sera were tested for alphavirus antibodies and viral RNA. In addition, by applying a catalytic, force-of-infection (FOI) statistical model to two serosurveys from Darien Province in 2012 and 2017, we investigated whether endemic or epidemic alphavirus transmission occurred historically. In 2017, MADV and VEEV IgM seroprevalences were 1.6% and 4.4%, respectively; IgG antibody prevalences were MADV: 13.2%, VEEV: 16.8%, Una virus (UNAV): 16.0%, and Mayaro virus: 1.1%. Active viral circulation was not detected. Evidence of MADV and UNAV infection was found near households, raising questions about its vectors and enzootic transmission cycles. Insomnia was associated with MADV and VEEV infections, depression symptoms were associated with MADV, and dizziness with VEEV and UNAV. Force-of-infection analyses suggest endemic alphavirus transmission historically, with recent increased human exposure to MADV and VEEV in Aruza and Mercadeo, respectively. The lack of additional neurological cases suggests that severe MADV and VEEV infections occur only rarely. Our results indicate that over the past five decades, alphavirus infections have occurred at low levels in eastern Panama, but that MADV and VEEV infections have recently increased-potentially during the past decade. Endemic infections and outbreaks of MADV and VEEV appear to differ spatially in some locations of eastern Panama.

An exact method for quantifying the reliability of end-of-epidemic declarations in real time

We derive and validate a novel and analytic method for estimating the probability that an epidemic has been eliminated (i.e. that no future local cases will emerge) in real time. When this probability crosses 0.95 an outbreak can be declared over with 95% confidence. Our method is easy to compute, only requires knowledge of the incidence curve and the serial interval distribution, and evaluates the statistical lifetime of the outbreak of interest. Using this approach, we show how the time-varying under-reporting of infected cases will artificially inflate the inferred probability of elimination, leading to premature (false-positive) end-of-epidemic declarations. Contrastingly, we prove that incorrectly identifying imported cases as local will deceptively decrease this probability, resulting in delayed (false-negative) declarations. Failing to sustain intensive surveillance during the later phases of an epidemic can therefore substantially mislead policymakers on when it is safe to remove travel bans or relax quarantine and social distancing advisories. World Health Organisation guidelines recommend fixed (though disease-specific) waiting times for end-of-epidemic declarations that cannot accommodate these variations. Consequently, there is an unequivocal need for more active and specialised metrics for reliably identifying the conclusion of an epidemic.

Adaptive Estimation for Epidemic Renewal and Phylogenetic Skyline Models

Estimating temporal changes in a target population from phylogenetic or count data is an important problem in ecology and epidemiology. Reliable estimates can provide key insights into the climatic and biological drivers influencing the diversity or structure of that population and evidence hypotheses concerning its future growth or decline. In infectious disease applications, the individuals infected across an epidemic form the target population. The renewal model estimates the effective reproduction number, R, of the epidemic from counts of observed incident cases. The skyline model infers the effective population size, N, underlying a phylogeny of sequences sampled from that epidemic. Practically, R measures ongoing epidemic growth while N informs on historical caseload. While both models solve distinct problems, the reliability of their estimates depends on p-dimensional piecewise-constant functions. If p is misspecified, the model might underfit significant changes or overfit noise and promote a spurious understanding of the epidemic, which might misguide intervention policies or misinform forecasts. Surprisingly, no transparent yet principled approach for optimizing p exists. Usually, p is heuristically set, or obscurely controlled via complex algorithms. We present a computable and interpretable p-selection method based on the minimum description length (MDL) formalism of information theory. Unlike many standard model selection techniques, MDL accounts for the additional statistical complexity induced by how parameters interact. As a result, our method optimizes p so that R and N estimates properly and meaningfully adapt to available data. It also outperforms comparable Akaike and Bayesian information criteria on several classification problems, given minimal knowledge of the parameter space, and exposes statistical similarities among renewal, skyline, and other models in biology. Rigorous and interpretable model selection is necessary if trustworthy and justifiable conclusions are to be drawn from piecewise models. [Coalescent processes; epidemiology; information theory; model selection; phylodynamics; renewal models; skyline plots].

Inference of COVID-19 epidemiological distributions from Brazilian hospital data

Knowing COVID-19 epidemiological distributions, such as the time from patient admission to death, is directly relevant to effective primary and secondary care planning, and moreover, the mathematical modelling of the pandemic generally. We determine epidemiological distributions for patients hospitalized with COVID-19 using a large dataset (N = 21 000 - 157 000) from the Brazilian Sistema de Informação de Vigilância Epidemiológica da Gripe database. A joint Bayesian subnational model with partial pooling is used to simultaneously describe the 26 states and one federal district of Brazil, and shows significant variation in the mean of the symptom-onset-to-death time, with ranges between 11.2 and 17.8 days across the different states, and a mean of 15.2 days for Brazil. We find strong evidence in favour of specific probability density function choices: for example, the gamma distribution gives the best fit for onset-to-death and the generalized lognormal for onset-to-hospital-admission. Our results show that epidemiological distributions have considerable geographical variation, and provide the first estimates of these distributions in a low and middle-income setting. At the subnational level, variation in COVID-19 outcome timings are found to be correlated with poverty, deprivation and segregation levels, and weaker correlation is observed for mean age, wealth and urbanicity.