Search

The 2014 chikungunya outbreak in the Dominican Republic resulted in intense local transmission, with high postoutbreak seroprevalence. The resulting population immunity will likely minimize risk for another large outbreak through 2035, but changes in population behavior or environmental conditions or emergence of different virus strains could lead to increased transmission.

During 2020–2022, players and staff in the English Premier League in the United Kingdom were tested regularly for SARS-CoV-2 with the aim of creating a biosecure bubble for each team. We found that prevalence and reinfection estimates were consistent with those from other studies and with community infection trends.

Summary: Background: COVID-19 dynamics are driven by a complex interplay of factors including population behaviour, new variants, vaccination and immunity from prior infections. We quantify drivers of SARS-CoV-2 transmission in the Dominican Republic, an upper-middle income country of 10.8 million people. We then assess the impact of the vaccination campaign implemented in February 2021, primarily using CoronaVac, in saving lives and averting hospitalisations. Methods: We fit an age-structured, multi-variant transmission dynamic model to reported deaths, hospital bed occupancy, and seroprevalence data until December 2021, and simulate epidemic trajectories under different counterfactual scenarios. Findings: We estimate that vaccination averted 7210 hospital admissions (95% credible interval, CrI: 6830–7600), 2180 intensive care unit admissions (95% CrI: 2080–2280) and 766 deaths (95% CrI: 694–859) in the first 6 months of the campaign. If no vaccination had occurred, we estimate that an additional decrease of 10–20% in population mobility would have been required to maintain equivalent death and hospitalisation outcomes. We also found that early vaccination with CoronaVac was preferable to delayed vaccination using a product with higher efficacy. Interpretation: SARS-CoV-2 transmission dynamics in the Dominican Republic were driven by a substantial accumulation of immunity during the first two years of the pandemic but, despite this, vaccination was essential in enabling a return to pre-pandemic mobility levels without considerable additional morbidity and mortality. Funding: Medical Research Council, Wellcome Trust, Royal Society, US CDC and Australian National Health and Medical Research Council.

Abstract Background Recurring COVID-19 waves highlight the need for tools able to quantify transmission risk, and identify geographical areas at risk of outbreaks. Local outbreak risk depends on complex immunity patterns resulting from previous infections, vaccination, waning and immune escape, alongside other factors (population density, social contact patterns). Immunity patterns are spatially and demographically heterogeneous, and are challenging to capture in country-level forecast models. Methods We used a spatiotemporal regression model to forecast subnational case and death counts and applied it to three EU countries as test cases: France, Czechia, and Italy. Cases in local regions arise from importations or local transmission. Our model produces age-stratified forecasts given age-stratified data, and links reported case counts to routinely collected covariates (e.g. test number, vaccine coverage). We assessed the predictive performance of our model up to four weeks ahead using proper scoring rules and compared it to the European COVID-19 Forecast Hub ensemble model. Using simulations, we evaluated the impact of variations in transmission on the forecasts. We developed an open-source RShiny App to visualise the forecasts and scenarios. Results At a national level, the median relative difference between our median weekly case forecasts and the data up to four weeks ahead was 25% (IQR: 12–50%) over the prediction period. The accuracy decreased as the forecast horizon increased (on average 24% increase in the median ranked probability score per added week), while the accuracy of death forecasts was more stable. Beyond two weeks, the model generated a narrow range of likely transmission dynamics. The median national case forecasts showed similar accuracy to forecasts from the European COVID-19 Forecast Hub ensemble model, but the prediction interval was narrower in our model. Generating forecasts under alternative transmission scenarios was therefore key to capturing the range of possible short-term transmission dynamics. Discussion Our model captures changes in local COVID-19 outbreak dynamics, and enables quantification of short-term transmission risk at a subnational level. The outputs of the model improve our ability to identify areas where outbreaks are most likely, and are available to a wide range of public health professionals through the Shiny App we developed.

To support the ongoing management of viral respiratory diseases while transitioning out of the acute phase of the COVID-19 pandemic, many countries are moving toward an integrated model of surveillance for SARS-CoV-2, influenza virus, and other respiratory pathogens. Although many surveillance approaches catalyzed by the COVID-19 pandemic provide novel epidemiologic insight, continuing them as implemented during the pandemic is unlikely to be feasible for nonemergency surveillance, and many have already been scaled back. Furthermore, given anticipated cocirculation of SARS-CoV-2 and influenza virus, surveillance activities in place before the pandemic require review and adjustment to ensure their ongoing value for public health. In this report, we highlight key challenges for the development of integrated models of surveillance. We discuss the relative strengths and limitations of different surveillance practices and studies as well as their contribution to epidemiologic assessment, forecasting, and public health decision-making.

Abstract Estimating the impact of vaccination and non-pharmaceutical interventions on COVID-19 incidence is complicated by several factors, including successive emergence of SARS-CoV-2 variants of concern and changing population immunity from vaccination and infection. We develop an age-structured multi-strain COVID-19 transmission model and inference framework to estimate vaccination and non-pharmaceutical intervention impact accounting for these factors. We apply this framework to COVID-19 waves in French Polynesia and estimate that the vaccination programme averted 34.8% (95% credible interval: 34.5–35.2%) of 223,000 symptomatic cases, 49.6% (48.7–50.5%) of 5830 hospitalisations and 64.2% (63.1–65.3%) of 1540 hospital deaths that would have occurred in a scenario without vaccination up to May 2022. We estimate the booster campaign contributed 4.5%, 1.9%, and 0.4% to overall reductions in cases, hospitalisations, and deaths. Our results suggest that removing lockdowns during the first two waves would have had non-linear effects on incidence by altering accumulation of population immunity. Our estimates of vaccination and booster impact differ from those for other countries due to differences in age structure, previous exposure levels and timing of variant introduction relative to vaccination, emphasising the importance of detailed analysis that accounts for these factors.

Abstract Background Understanding the overall effectiveness of non-pharmaceutical interventions to control the COVID-19 pandemic and reduce the burden of disease is crucial for future pandemic planning. However, quantifying the effectiveness of specific control measures and the extent of missed infections, in the absence of early large-scale serological surveys or random community testing, has remained challenging. Methods Combining data on notified local COVID-19 cases with known and unknown sources of infections in Singapore with a branching process model, we reconstructed the incidence of missed infections during the early phase of the wild-type SARS-CoV-2 and Delta variant transmission. We then estimated the relative effectiveness of border control measures, case finding and contact tracing when there was no or low vaccine coverage in the population. We compared the risk of ICU admission and death between the wild-type SARS-CoV-2 and the Delta variant in notified cases and all infections. Results We estimated strict border control measures were associated with 0.2 (95% credible intervals, CrI 0.04–0.8) missed imported infections per notified case between July and December 2020, a decline from around 1 missed imported infection per notified case in the early phases of the pandemic. Contact tracing was estimated to identify 78% (95% CrI 62–93%) of the secondary infections generated by notified cases before the partial lockdown in Apr 2020, but this declined to 63% (95% CrI 56–71%) during the lockdown and rebounded to 78% (95% CrI 58–94%) during reopening in Jul 2020. The contribution of contact tracing towards overall outbreak control also hinges on ability to find cases with unknown sources of infection: 42% (95% CrI 12–84%) of such cases were found prior to the lockdown; 10% (95% CrI 7–15%) during the lockdown; 47% (95% CrI 17–85%) during reopening, due to increased testing capacity and health-seeking behaviour. We estimated around 63% (95% CrI 49–78%) of the wild-type SARS-CoV-2 infections were undetected during 2020 and around 70% (95% CrI 49–91%) for the Delta variant in 2021. Conclusions Combining models with case linkage data enables evaluation of the effectiveness of different components of outbreak control measures, and provides more reliable situational awareness when some cases are missed. Using such approaches for early identification of the weakest link in containment efforts could help policy makers to better redirect limited resources to strengthen outbreak control.

Here, the authors simulate COVID-19 outbreaks on an empirical contact network derived from digital contact data collected on cruise ships. They model impacts of different control measures and find that combinations of measures, particularly vaccination and rapid antigen testing, are important for mitigating outbreaks.

Abstract Background Subnational heterogeneity in immunity to measles can create pockets of susceptibility and result in long-lasting outbreaks despite high levels of national vaccine coverage. The elimination status defined by the World Health Organization aims to identify countries where the virus is no longer circulating and can be verified after 36 months of interrupted transmission. However, since 2018, numerous countries have lost their elimination status soon after reaching it, showing that the indicators defining elimination may not be associated with lower risks of outbreaks. Methods We quantified the impact of local vaccine coverage and recent levels of incidence on the dynamics of measles in each French department between 2009 and 2018, using mathematical models based on the “Endemic-Epidemic” regression framework. After fitting the models using daily case counts, we simulated the effect of variations in the vaccine coverage and recent incidence on future transmission. Results High values of local vaccine coverage were associated with fewer imported cases and lower risks of local transmissions, but regions that had recently reported high levels of incidence were also at a lower risk of local transmission. This may be due to additional immunity accumulated during recent outbreaks. Therefore, the risk of local transmission was not lower in areas fulfilling the elimination criteria. A decrease of 3% in the 3-year average vaccine uptake led to a fivefold increase in the average annual number of cases in simulated outbreaks. Conclusions Local vaccine uptake was a reliable indicator of the intensity of transmission in France, even if it only describes yearly coverage in a given age group, and ignores population movements. Therefore, spatiotemporal variations in vaccine coverage, caused by disruptions in routine immunisation programmes, or lower trust in vaccines, can lead to large increases in both local and cross-regional transmission. The incidence indicator used to define the elimination status was not associated with a lower number of local transmissions in France, and may not illustrate the risks of imminent outbreaks. More detailed models of local immunity levels or subnational seroprevalence studies may yield better estimates of local risk of measles outbreaks.

Abstract Background Stratifying dengue risk within endemic countries is crucial for allocating limited control interventions. Current methods of monitoring dengue transmission intensity rely on potentially inaccurate incidence estimates. We investigated whether incidence or alternate metrics obtained from standard, or laboratory, surveillance operations represent accurate surrogate indicators of the burden of dengue and can be used to monitor the force of infection (FOI) across urban centres. Methods Among those who reported and resided in 13 cities across the Philippines, we collected epidemiological data from all dengue case reports between 2014 and 2017 (N 80,043) and additional laboratory data from a cross-section of sampled case reports (N 11,906) between 2014 and 2018. At the city level, we estimated the aggregated annual FOI from age-accumulated IgG among the non-dengue reporting population using catalytic modelling. We compared city-aggregated FOI estimates to aggregated incidence and the mean age of clinically and laboratory diagnosed dengue cases using Pearson’s Correlation coefficient and generated predicted FOI estimates using regression modelling. Results We observed spatial heterogeneity in the dengue average annual FOI across sampled cities, ranging from 0.054 [0.036–0.081] to 0.249 [0.223–0.279]. Compared to FOI estimates, the mean age of primary dengue infections had the strongest association (ρ −0.848, p value

Abstract Background Routine asymptomatic testing using RT-PCR of people who interact with vulnerable populations, such as medical staff in hospitals or care workers in care homes, has been employed to help prevent outbreaks among vulnerable populations. Although the peak sensitivity of RT-PCR can be high, the probability of detecting an infection will vary throughout the course of an infection. The effectiveness of routine asymptomatic testing will therefore depend on testing frequency and how PCR detection varies over time. Methods We fitted a Bayesian statistical model to a dataset of twice weekly PCR tests of UK healthcare workers performed by self-administered nasopharyngeal swab, regardless of symptoms. We jointly estimated times of infection and the probability of a positive PCR test over time following infection; we then compared asymptomatic testing strategies by calculating the probability that a symptomatic infection is detected before symptom onset and the probability that an asymptomatic infection is detected within 7 days of infection. Results We estimated that the probability that the PCR test detected infection peaked at 77% (54–88%) 4 days after infection, decreasing to 50% (38–65%) by 10 days after infection. Our results suggest a substantially higher probability of detecting infections 1–3 days after infection than previously published estimates. We estimated that testing every other day would detect 57% (33–76%) of symptomatic cases prior to onset and 94% (75–99%) of asymptomatic cases within 7 days if test results were returned within a day. Conclusions Our results suggest that routine asymptomatic testing can enable detection of a high proportion of infected individuals early in their infection, provided that the testing is frequent and the time from testing to notification of results is sufficiently fast.

Dengue and Zika virus are closely related flaviviruses but can have contrasting transmission dynamics in the same populations. Here, the authors use a model combining serological, surveillance and viral sequence data to explain differences in transmission dynamics in Fiji.

Abstract Predicting the likelihood of wildlife presence at potential wildlife–livestock interfaces is challenging. These interfaces are usually relatively small geographical areas where landscapes show large variation over small distances. Models of wildlife distribution based on coarse data over wide geographical ranges may not be representative of these interfaces. High‐resolution data can help identify fine‐scale predictors of wildlife habitat use at a local scale and provide more accurate predictions of species habitat use. These data may be used to inform knowledge of interface risks, such as disease transmission between wildlife and livestock, or human–wildlife conflict. This study uses fine‐scale habitat use data from wild boar (Sus scrofa) based on activity signs and direct field observations in and around the Forest of Dean in Gloucestershire, England. Spatial logistic regression models fitted using a variant of penalized quasi‐likelihood were used to identify habitat‐based and anthropogenic predictors of wild boar signs. Our models showed that within the Forest of Dean, wild boar signs were more likely to be seen in spring, in forest‐type habitats, closer to the center of the forest and near litter bins. In the area surrounding the Forest of Dean, wild boar signs were more likely to be seen in forest‐type habitats and near recreational parks and less likely to be seen near livestock. This approach shows that wild boar habitat use can be predicted using fine‐scale data over comparatively small areas and in human‐dominated landscapes, while taking account of the spatial correlation from other nearby fine‐scale data‐points. The methods we use could be applied to map habitat use of other wildlife species in similar landscapes, or of movement‐restricted, isolated, or fragmented wildlife populations.

Abstract Background In dengue-endemic countries, targeting limited control interventions to populations at risk of severe disease could enable increased efficiency. Individuals who have had their first (primary) dengue infection are at risk of developing more severe secondary disease, thus could be targeted for disease prevention. Currently, there is no reliable algorithm for determining primary and post-primary (infection with more than one flavivirus) status from a single serum sample. In this study, we developed and validated an immune status algorithm using single acute serum samples from reporting patients and investigated dengue immuno-epidemiological patterns across the Philippines. Methods During 2015/2016, a cross-sectional sample of 10,137 dengue case reports provided serum for molecular (anti-DENV PCR) and serological (anti-DENV IgM/G capture ELISA) assay. Using mixture modelling, we re-assessed IgM/G seroprevalence and estimated functional, disease day-specific, IgG:IgM ratios that categorised the reporting population as negative, historical, primary and post-primary for dengue. We validated our algorithm against WHO gold standard criteria and investigated cross-reactivity with Zika by assaying a random subset for anti-ZIKV IgM and IgG. Lastly, using our algorithm, we explored immuno-epidemiological patterns of dengue across the Philippines. Results Our modelled IgM and IgG seroprevalence thresholds were lower than kit-provided thresholds. Individuals anti-DENV PCR+ or IgM+ were classified as active dengue infections (83.1%, 6998/8425). IgG− and IgG+ active dengue infections on disease days 1 and 2 were categorised as primary and post-primary, respectively, while those on disease days 3 to 5 with IgG:IgM ratios below and above 0.45 were classified as primary and post-primary, respectively. A significant proportion of post-primary dengue infections had elevated anti-ZIKV IgG inferring previous Zika exposure. Our algorithm achieved 90.5% serological agreement with WHO standard practice. Post-primary dengue infections were more likely to be older and present with severe symptoms. Finally, we identified a spatio-temporal cluster of primary dengue case reporting in northern Luzon during 2016. Conclusions Our dengue immune status algorithm can equip surveillance operations with the means to target dengue control efforts. The algorithm accurately identified primary dengue infections who are at risk of future severe disease.

Abstract Background The health impact of COVID-19 may differ in African settings as compared to countries in Europe or China due to demographic, epidemiological, environmental and socio-economic factors. We evaluated strategies to reduce SARS-CoV-2 burden in African countries, so as to support decisions that balance minimising mortality, protecting health services and safeguarding livelihoods. Methods We used a Susceptible-Exposed-Infectious-Recovered mathematical model, stratified by age, to predict the evolution of COVID-19 epidemics in three countries representing a range of age distributions in Africa (from oldest to youngest average age: Mauritius, Nigeria and Niger), under various effectiveness assumptions for combinations of different non-pharmaceutical interventions: self-isolation of symptomatic people, physical distancing and ‘shielding’ (physical isolation) of the high-risk population. We adapted model parameters to better represent uncertainty about what might be expected in African populations, in particular by shifting the distribution of severity risk towards younger ages and increasing the case-fatality ratio. We also present sensitivity analyses for key model parameters subject to uncertainty. Results We predicted median symptomatic attack rates over the first 12 months of 23% (Niger) to 42% (Mauritius), peaking at 2–4 months, if epidemics were unmitigated. Self-isolation while symptomatic had a maximum impact of about 30% on reducing severe cases, while the impact of physical distancing varied widely depending on percent contact reduction and R 0. The effect of shielding high-risk people, e.g. by rehousing them in physical isolation, was sensitive mainly to residual contact with low-risk people, and to a lesser extent to contact among shielded individuals. Mitigation strategies incorporating self-isolation of symptomatic individuals, moderate physical distancing and high uptake of shielding reduced predicted peak bed demand and mortality by around 50%. Lockdowns delayed epidemics by about 3 months. Estimates were sensitive to differences in age-specific social mixing patterns, as published in the literature, and assumptions on transmissibility, infectiousness of asymptomatic cases and risk of severe disease or death by age. Conclusions In African settings, as elsewhere, current evidence suggests large COVID-19 epidemics are expected. However, African countries have fewer means to suppress transmission and manage cases. We found that self-isolation of symptomatic persons and general physical distancing are unlikely to avert very large epidemics, unless distancing takes the form of stringent lockdown measures. However, both interventions help to mitigate the epidemic. Shielding of high-risk individuals can reduce health service demand and, even more markedly, mortality if it features high uptake and low contact of shielded and unshielded people, with no increase in contact among shielded people. Strategies combining self-isolation, moderate physical distancing and shielding could achieve substantial reductions in mortality in African countries. Temporary lockdowns, where socioeconomically acceptable, can help gain crucial time for planning and expanding health service capacity.

Abstract Background Asymptomatic or subclinical SARS-CoV-2 infections are often unreported, which means that confirmed case counts may not accurately reflect underlying epidemic dynamics. Understanding the level of ascertainment (the ratio of confirmed symptomatic cases to the true number of symptomatic individuals) and undetected epidemic progression is crucial to informing COVID-19 response planning, including the introduction and relaxation of control measures. Estimating case ascertainment over time allows for accurate estimates of specific outcomes such as seroprevalence, which is essential for planning control measures. Methods Using reported data on COVID-19 cases and fatalities globally, we estimated the proportion of symptomatic cases (i.e. any person with any of fever ≥ 37.5 °C, cough, shortness of breath, sudden onset of anosmia, ageusia or dysgeusia illness) that were reported in 210 countries and territories, given those countries had experienced more than ten deaths. We used published estimates of the baseline case fatality ratio (CFR), which was adjusted for delays and under-ascertainment, then calculated the ratio of this baseline CFR to an estimated local delay-adjusted CFR to estimate the level of under-ascertainment in a particular location. We then fit a Bayesian Gaussian process model to estimate the temporal pattern of under-ascertainment. Results Based on reported cases and deaths, we estimated that, during March 2020, the median percentage of symptomatic cases detected across the 84 countries which experienced more than ten deaths ranged from 2.4% (Bangladesh) to 100% (Chile). Across the ten countries with the highest number of total confirmed cases as of 6 July 2020, we estimated that the peak number of symptomatic cases ranged from 1.4 times (Chile) to 18 times (France) larger than reported. Comparing our model with national and regional seroprevalence data where available, we find that our estimates are consistent with observed values. Finally, we estimated seroprevalence for each country. As of 7 June, our seroprevalence estimates range from 0% (many countries) to 13% (95% CrI 5.6–24%) (Belgium). Conclusions We found substantial under-ascertainment of symptomatic cases, particularly at the peak of the first wave of the SARS-CoV-2 pandemic, in many countries. Reported case counts will therefore likely underestimate the rate of outbreak growth initially and underestimate the decline in the later stages of an epidemic. Although there was considerable under-reporting in many locations, our estimates were consistent with emerging serological data, suggesting that the proportion of each country’s population infected with SARS-CoV-2 worldwide is generally low.

Various interventions for live poultry markets (LPMs) have emerged to control outbreaks of avian influenza A(H7N9) virus in mainland China since March 2013. We assessed the effectiveness of various LPM interventions in reducing transmission of H7N9 virus across 5 annual waves during 2013–2018, especially in the final wave. With the exception of waves 1 and 4, various LPM interventions reduced daily incidence rates significantly across waves. Four LPM interventions led to a mean reduction of 34%–98% in the daily number of infections in wave 5. Of these, permanent closure provided the most effective reduction in human infection with H7N9 virus, followed by long-period, short-period, and recursive closures in wave 5. The effectiveness of various LPM interventions changed with the type of intervention across epidemics. Permanent LPM closure should be considered to maintain sufficient effectiveness of interventions and prevent the recurrence of H7N9 epidemics.

Background: The duration of immunity against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is still uncertain, but it is of key clinical and epidemiological importance. Seasonal human coronaviruses (HCoV) have been circulating for longer and, therefore, may offer insights into the long-term dynamics of reinfection for such viruses. Methods: Combining historical seroprevalence data from five studies covering the four circulating HCoVs with an age-structured reverse catalytic model, we estimated the likely duration of seropositivity following seroconversion. Results: We estimated that antibody persistence lasted between 0.9 (95% Credible interval: 0.6 - 1.6) and 3.8 (95% CrI: 2.0 - 7.4) years. Furthermore, we found the force of infection in older children and adults (those over 8.5 [95% CrI: 7.5 - 9.9] years) to be higher compared with young children in the majority of studies. Conclusions: These estimates of endemic HCoV dynamics could provide an indication of the future long-term infection and reinfection patterns of SARS-CoV-2.

Background Targeting interventions to areas that have recently experienced cases of disease is one strategy to contain outbreaks of infectious disease. Such case-area targeted interventions (CATI) have become an increasingly popular approach for dengue control but there is little evidence to suggest how precisely targeted or how recent cases need to be, to mount an effective response. The growing interest in the development of prophylactic and therapeutic drugs for dengue has also given new relevance for CATI strategies to interrupt transmission or deliver early treatment. Methods/Principal findings Here we develop a patch-based mathematical model of spatial dengue spread and fit it to spatiotemporal datasets from Singapore. Simulations from this model suggest CATI strategies could be effective, particularly if used in lower density areas. To maximise effectiveness, increasing the size of the radius around an index case should be prioritised even if it results in delays in the intervention being applied. This is partially because large intervention radii ensure individuals receive multiple and regular rounds of drug dosing or vector control, and thus boost overall coverage. Given equivalent efficacy, CATIs using prophylactic drugs are predicted to be more effective than adult mosquito-killing vector control methods and may even offer the possibility of interrupting individual chains of transmission if rapidly deployed. CATI strategies quickly lose their effectiveness if baseline transmission increases or case detection rates fall. Conclusions/Significance These results suggest CATI strategies can play an important role in dengue control but are likely to be most relevant for low transmission areas where high coverage of other non-reactive interventions already exists. Controlled field trials are needed to assess the field efficacy and practical constraints of large operational CATI strategies. Author summary In resource limited settings there is a pressing need for more efficient, more targeted ways of controlling transmission and preventing outbreaks. One option is to use case-area targeted interventions (CATI) that are focussed on areas that have recently reported disease cases. The effectiveness of such CATI strategies is highly dependent on how the disease spreads. Despite CATI strategies being widely used to control the vector-transmitted disease dengue, little evidence underpins its effectiveness. In this analysis we formulate a mathematical model designed to test the effectiveness of CATI strategies for dengue control in Singapore- a best case test scenario for the approach. Simulation from this model suggested CATI are likely to be effective for dengue, but need to have large (250m+) radii around index cases and may not be suitable in higher transmission areas. These results, when combined with limited field evidence of efficacy, suggest that CATI strategies are unlikely to be universally applicable dengue control tools. Only once high coverage with other (non-reactive) interventions has been achieved and comprehensive rapid disease surveillance has been established are CATI strategies likely to become efficient methods of disease control.

Objectives: In Fiji, autochthonous chikungunya virus (CHIKV) infection was first detected in March 2015. In a previous serosurvey conducted during October–November 2015, we reported a prevalence of anti-CHIKV IgG antibodies of 0.9%. In the present study, we investigated the seroprevalence of CHIKV two years after its emergence in Fiji. Methods: Sera from 320 residents of Fiji recruited in June 2017, from the same cohort of individuals that participated in the serosurvey in 2015, were tested for the presence of IgG antibodies against CHIKV using a recombinant antigen-based microsphere immunoassay. Results: Between 2015 and 2017, CHIKV seroprevalence among residents increased from 0.9% (3/333) to 12.8% (41/320). Of the participants with available serum samples collected in both 2015 and 2017 (n = 200), 31 (15.5%) who were seronegative in 2015 had seroconverted to CHIKV in 2017. Conclusions: Our findings suggest that low-level transmission of CHIKV occurred during the two years following the emergence of the virus in Fiji. No CHIKV infection has been reported in Fiji since 2017, but due to the presumed low herd immunity of the population, the risk of CHIKV re-emergence is high. Consequently, chikungunya should be considered in the differential diagnosis of acute febrile diseases in Fiji. Keywords: Chikungunya, Ross River, Arbovirus, Seroprevalence, Fiji, Pacific

Abstract Background Wolbachia-infected mosquitoes reduce dengue virus transmission, and city-wide releases in Yogyakarta city, Indonesia, are showing promising entomological results. Accurate estimates of the burden of dengue, its spatial distribution and the potential impact of Wolbachia are critical in guiding funder and government decisions on its future wider use. Methods Here, we combine multiple modelling methods for burden estimation to predict national case burden disaggregated by severity and map the distribution of burden across the country using three separate data sources. An ensemble of transmission models then predicts the estimated reduction in dengue transmission following a nationwide roll-out of wMel Wolbachia. Results We estimate that 7.8 million (95% uncertainty interval [UI] 1.8–17.7 million) symptomatic dengue cases occurred in Indonesia in 2015 and were associated with 332,865 (UI 94,175–754,203) lost disability-adjusted life years (DALYs). The majority of dengue’s burden was due to non-severe cases that did not seek treatment or were challenging to diagnose in outpatient settings leading to substantial underreporting. Estimated burden was highly concentrated in a small number of large cities with 90% of dengue cases occurring in 15.3% of land area. Implementing a nationwide Wolbachia population replacement programme was estimated to avert 86.2% (UI 36.2–99.9%) of cases over a long-term average. Conclusions These results suggest interventions targeted to the highest burden cities can have a disproportionate impact on dengue burden. Area-wide interventions, such as Wolbachia, that are deployed based on the area covered could protect people more efficiently than individual-based interventions, such as vaccines, in such dense environments.

Zika and chikungunya viruses were first detected in Fiji in 2015. Examining surveillance and phylogenetic and serologic data, we found evidence of low-level transmission of Zika and chikungunya viruses during 2013–2017, in contrast to the major outbreaks caused by closely related virus strains in other Pacific Island countries.

A unique outbreak of Ross River virus (RRV) infection was reported in Fiji in 1979. In 2013, RRV seroprevalence among residents was 46.5% (362/778). Of the residents who were seronegative in 2013 and retested in 2015, 10.9% (21/192) had seroconverted to RRV, suggesting ongoing endemic circulation of RRV in Fiji.

Abstract Background Between August and December 2017, more than 625,000 Rohingya from Myanmar fled into Bangladesh, settling in informal makeshift camps in Cox’s Bazar district and joining 212,000 Rohingya already present. In early November, a diphtheria outbreak hit the camps, with 440 reported cases during the first month. A rise in cases during early December led to a collaboration between teams from Médecins sans Frontières—who were running a provisional diphtheria treatment centre—and the London School of Hygiene and Tropical Medicine with the goal to use transmission dynamic models to forecast the potential scale of the outbreak and the resulting resource needs. Methods We first adjusted for delays between symptom onset and case presentation using the observed distribution of reporting delays from previously reported cases. We then fit a compartmental transmission model to the adjusted incidence stratified by age group and location. Model forecasts with a lead time of 2 weeks were issued on 12, 20, 26 and 30 December and communicated to decision-makers. Results The first forecast estimated that the outbreak would peak on 19 December in Balukhali camp with 303 (95% posterior predictive interval 122–599) cases and would continue to grow in Kutupalong camp, requiring a bed capacity of 316 (95% posterior predictive interval (PPI) 197–499). On 19 December, a total of 54 cases were reported, lower than forecasted. Subsequent forecasts were more accurate: on 20 December, we predicted a total of 912 cases (95% PPI 367–2183) and 136 (95% PPI 55–327) hospitalizations until the end of the year, with 616 cases actually reported during this period. Conclusions Real-time modelling enabled feedback of key information about the potential scale of the epidemic, resource needs and mechanisms of transmission to decision-makers at a time when this information was largely unknown. By 20 December, the model generated reliable forecasts and helped support decision-making on operational aspects of the outbreak response, such as hospital bed and staff needs, and with advocacy for control measures. Although modelling is only one component of the evidence base for decision-making in outbreak situations, suitable analysis and forecasting techniques can be used to gain insights into an ongoing outbreak.

Background: The duration of immunity against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is still uncertain, but it is of key clinical and epidemiological importance. Seasonal human coronaviruses (HCoV) have been circulating for longer and, therefore, may offer insights into the long-term dynamics of reinfection for such viruses. Methods: Combining historical seroprevalence data from five studies covering the four circulating HCoVs with an age-structured reverse catalytic model, we estimated the likely duration of seropositivity following seroconversion. Results: We estimated that antibody persistence lasted between 0.9 (95% Credible interval: 0.6 - 1.6) and 3.8 (95% CrI: 2.0 - 7.4) years. Furthermore, we found the force of infection in older children and adults (those over 8.5 [95% CrI: 7.5 - 9.9] years) to be higher compared with young children in the majority of studies. Conclusions: These estimates of endemic HCoV dynamics could provide an indication of the future long-term infection and reinfection patterns of SARS-CoV-2.

Introduction: Contact tracing has the potential to control outbreaks without the need for stringent physical distancing policies, e.g. civil lockdowns. Unlike forward contact tracing, backward contact tracing identifies the source of newly detected cases. This approach is particularly valuable when there is high individual-level variation in the number of secondary transmissions (overdispersion). Methods: By using a simple branching process model, we explored the potential of combining backward contact tracing with more conventional forward contact tracing for control of COVID-19. We estimated the typical size of clusters that can be reached by backward tracing and simulated the incremental effectiveness of combining backward tracing with conventional forward tracing. Results: Across ranges of parameter values consistent with dynamics of SARS-CoV-2, backward tracing is expected to identify a primary case generating 3-10 times more infections than a randomly chosen case, typically increasing the proportion of subsequent cases averted by a factor of 2-3. The estimated number of cases averted by backward tracing became greater with a higher degree of overdispersion. Conclusion: Backward contact tracing can be an effective tool for outbreak control, especially in the presence of overdispersion as is observed with SARS-CoV-2.

Introduction: Contact tracing has the potential to control outbreaks without the need for stringent physical distancing policies, e.g. civil lockdowns. Unlike forward contact tracing, backward contact tracing identifies the source of newly detected cases. This approach is particularly valuable when there is high individual-level variation in the number of secondary transmissions (overdispersion). Methods: By using a simple branching process model, we explored the potential of combining backward contact tracing with more conventional forward contact tracing for control of COVID-19. We estimated the typical size of clusters that can be reached by backward tracing and simulated the incremental effectiveness of combining backward tracing with conventional forward tracing. Results: Across ranges of parameter values consistent with dynamics of SARS-CoV-2, backward tracing is expected to identify a primary case generating 3-10 times more infections than average, typically increasing the proportion of subsequent cases averted by a factor of 2-3. The estimated number of cases averted by backward tracing became greater with a higher degree of overdispersion. Conclusion: Backward contact tracing can be an effective tool for outbreak control, especially in the presence of overdispersion as was observed with SARS-CoV-2.

Wild animals are the source of many pathogens of livestock and humans. Concerns about the potential transmission of economically important and zoonotic diseases from wildlife have led to increased surveillance at the livestock-wildlife interface. Knowledge of the types, frequency and duration of contacts between livestock and wildlife is necessary to identify risk factors for disease transmission and to design possible mitigation strategies. Observing the behaviour of many wildlife species is challenging due to their cryptic nature and avoidance of humans, meaning there are relatively few studies in this area. Further, a consensus on the definition of what constitutes a ‘contact’ between wildlife and livestock is lacking. A systematic review was conducted to investigate which livestock-wildlife contacts have been studied and why, as well as the methods used to observe each species. Over 30,000 publications were screened, of which 122 fulfilled specific criteria for inclusion in the analysis. The majority of studies examined cattle contacts with badgers or with deer; studies involving wild pig contacts with cattle or with domestic pigs were the next most frequent. There was a range of observational methods including motion-activated cameras and global positioning system collars. As a result of the wide variation and lack of consensus in the definitions of direct and indirect contacts, we developed a unified framework to define livestock-wildlife contacts that is sufficiently flexible to be applied to most wildlife and livestock species for non-vector-borne diseases. We hope this framework will help standardise the collection and reporting of contact data; a valuable step towards being able to compare the efficacy of wildlife-livestock observation methods. In doing so, it may aid the development of better disease transmission models and improve the design and effectiveness of interventions to reduce or prevent disease transmission.

Background: A novel coronavirus disease (COVID-19) outbreak has now spread to a number of countries worldwide. While sustained transmission chains of human-to-human transmission suggest high basic reproduction number R0, variation in the number of secondary transmissions (often characterised by so-called superspreading events) may be large as some countries have observed fewer local transmissions than others. Methods: We quantified individual-level variation in COVID-19 transmission by applying a mathematical model to observed outbreak sizes in affected countries. We extracted the number of imported and local cases in the affected countries from the World Health Organization situation report and applied a branching process model where the number of secondary transmissions was assumed to follow a negative-binomial distribution. Results: Our model suggested a high degree of individual-level variation in the transmission of COVID-19. Within the current consensus range of R0 (2-3), the overdispersion parameter k of a negative-binomial distribution was estimated to be around 0.1 (median estimate 0.1; 95% CrI: 0.05-0.2 for R0 = 2.5), suggesting that 80% of secondary transmissions may have been caused by a small fraction of infectious individuals (~10%). A joint estimation yielded likely ranges for R0 and k (95% CrIs: R0 1.4-12; k 0.04-0.2); however, the upper bound of R0 was not well informed by the model and data, which did not notably differ from that of the prior distribution. Conclusions: Our finding of a highly-overdispersed offspring distribution highlights a potential benefit to focusing intervention efforts on superspreading. As most infected individuals do not contribute to the expansion of an epidemic, the effective reproduction number could be drastically reduced by preventing relatively rare superspreading events.

Background: A novel coronavirus disease (COVID-19) outbreak has now spread to a number of countries worldwide. While sustained transmission chains of human-to-human transmission suggest high basic reproduction number R0, variation in the number of secondary transmissions (often characterised by so-called superspreading events) may be large as some countries have observed fewer local transmissions than others. Methods: We quantified individual-level variation in COVID-19 transmission by applying a mathematical model to observed outbreak sizes in affected countries. We extracted the number of imported and local cases in the affected countries from the World Health Organization situation report and applied a branching process model where the number of secondary transmissions was assumed to follow a negative-binomial distribution. Results: Our model suggested a high degree of individual-level variation in the transmission of COVID-19. Within the current consensus range of R0 (2-3), the overdispersion parameter k of a negative-binomial distribution was estimated to be around 0.1 (median estimate 0.1; 95% CrI: 0.05-0.2 for R0 = 2.5), suggesting that 80% of secondary transmissions may have been caused by a small fraction of infectious individuals (~10%). A joint estimation yielded likely ranges for R0 and k (95% CrIs: R0 1.4-12; k 0.04-0.2); however, the upper bound of R0 was not well informed by the model and data, which did not notably differ from that of the prior distribution. Conclusions: Our finding of a highly-overdispersed offspring distribution highlights a potential benefit to focusing intervention efforts on superspreading. As most infected individuals do not contribute to the expansion of an epidemic, the effective reproduction number could be drastically reduced by preventing relatively rare superspreading events.

We estimate the number of COVID-19 cases from newly reported deaths in a population without previous reports. Our results suggest that by the time a single death occurs, hundreds to thousands of cases are likely to be present in that population. This suggests containment via contact tracing will be challenging at this point, and other response strategies should be considered. Our approach is implemented in a publicly available, user-friendly, online tool.

Real-time forecasts of infectious diseases can help public health planning, especially during outbreaks. If forecasts are generated from mechanistic models, they can be further used to target resources or to compare the impact of possible interventions. However, paremeterising such models is often difficult in real time, when information on behavioural changes, interventions and routes of transmission are not readily available. Here, we present a semi-mechanistic model of infectious disease dynamics that was used in real time during the 2013–2016 West African Ebola epidemic, and show fits to a Ebola Forecasting Challenge conducted in late 2015 with simulated data mimicking the true epidemic. We assess the performance of the model in different situations and identify strengths and shortcomings of our approach. Models such as the one presented here which combine the power of mechanistic models with the flexibility to include uncertainty about the precise outbreak dynamics may be an important tool in combating future outbreaks. Keywords: Forecasting, Real-time modelling, Infectious disease dynamics, Outbreak

Evidence is increasing that Zika virus–related adverse outcomes can occur throughout pregnancy. Mathematical modeling analysis using reported outcome data suggests that surveillance for these outcomes should begin as soon as an outbreak is detected and should continue for 40 weeks after the outbreak ends.

Obesity is a key correlate of severe SARS-CoV-2 outcomes while the role of obesity on risk of SARS-CoV-2 infection, symptom phenotype, and immune response remain poorly defined. We examined data from a prospective SARS-CoV-2 cohort study to address these questions. Serostatus, body mass index, demographics, comorbidities, and prior COVID-19 compatible symptoms were assessed at baseline and serostatus and symptoms monthly thereafter. SARS-CoV-2 immunoassays included an IgG ELISA targeting the spike RBD, multiarray Luminex targeting 20 viral antigens, pseudovirus neutralization, and T cell ELISPOT assays. Our results from a large prospective SARS-CoV-2 cohort study indicate symptom phenotype is strongly influenced by obesity among younger but not older age groups; we did not identify evidence to suggest obese individuals are at higher risk of SARS-CoV-2 infection; and remarkably homogenous immune activity across BMI categories suggests immune protection across these groups may be similar.

Zika virus (ZIKV) exposure across flavivirus-endemic countries, including the Philippines, remains largely unknown despite sporadic case reporting and environmental suitability for transmission. Using laboratory surveillance data from 2016, 997 serum samples were randomly selected from suspected dengue (DENV) case reports across the Philippines and assayed for serological markers of short-term (IgM) and long-term (IgG) ZIKV exposure. Using mixture models, we re-evaluated ZIKV IgM/G seroprevalence thresholds and used catalytic models to quantify the force of infection (attack rate, AR) from age-accumulated ZIKV exposure. While we observed extensive ZIKV/DENV IgG cross-reactivity, not all individuals with active DENV presented with elevated ZIKV IgG, and a proportion of dengue-negative cases (DENV IgG-) were ZIKV IgG-positive (14.3%, 9/63). We identified evidence of long-term, yet not short-term, ZIKV exposure across Philippine regions (ZIKV IgG+: 31.5%, 314/997) which was geographically uncorrelated with DENV exposure. In contrast to the DENV AR (12.7% (95%CI: 9.1–17.4%)), the ZIKV AR was lower (5.7% (95%CI: 3–11%)) across the country. Our results provide evidence of widespread ZIKV exposure across the Philippines and suggest the need for studies to identify ZIKV infection risk factors over time to better prepare for potential future outbreaks.

Using an Ebola virus disease transmission model, we found that addition of ring vaccination at the outset of the West Africa epidemic might not have led to containment of this disease. However, in later stages of the epidemic or in outbreaks with less intense transmission or more effective control, this strategy could help eliminate the disease.

To determine transmission potential of influenza A(H7N9) virus, we used symptom onset data to compare 2 waves of infection in China during 2013–2014. We found evidence of increased transmission potential in the second wave and showed that live bird market closure was significantly less effective in Guangdong than in other regions.

In some parts of western Africa, Ebola treatment centers (ETCs) have reached capacity. Unless capacity is rapidly scaled up, the chance to avoid a generalized Ebola epidemic will soon diminish. The World Health Organization and partners are considering additional Ebola patient care options, including community care centers (CCCs), small, lightly staffed units that could be used to isolate patients outside the home and get them into care sooner than otherwise possible. Using a transmission model, we evaluated the benefits and risks of introducing CCCs into Sierra Leone’s Western Area, where most ETCs are at capacity. We found that use of CCCs could lead to a decline in cases, even if virus transmission occurs between CCC patients and the community. However, to prevent CCC amplification of the epidemic, the risk of Ebola virus–negative persons being exposed to virus within CCCs would have to be offset by a reduction in community transmission resulting from CCC use.

ABSTRACT During 2015 to 2016, Brazil reported more Zika virus (ZIKV) cases than any other country, yet population exposure remains unknown. Serological studies of ZIKV are hampered by cross-reactive immune responses against heterologous viruses. We conducted serosurveys for ZIKV, dengue virus (DENV), and Chikungunya virus (CHIKV) in 633 individuals prospectively sampled during 2015 to 2016, including microcephaly and non-microcephaly pregnancies, HIV-infected patients, tuberculosis patients, and university staff in Salvador in northeastern Brazil using enzyme-linked immunosorbent assays (ELISAs) and plaque reduction neutralization tests. Sera sampled retrospectively during 2013 to 2015 from 277 HIV-infected patients were used to assess the spread of ZIKV over time. Individuals were georeferenced, and sociodemographic indicators were compared between ZIKV-positive and -negative areas and areas with and without microcephaly cases. Epidemiological key parameters were modeled in a Bayesian framework. ZIKV seroprevalence increased rapidly during 2015 to 2016, reaching 63.3% by 2016 (95% confidence interval [CI], 59.4 to 66.8%), comparable to the seroprevalence of DENV (75.7%; CI, 69.4 to 81.1%) and higher than that of CHIKV (7.4%; CI, 5.6 to 9.8%). Of 19 microcephaly pregnancies, 94.7% showed ZIKV IgG antibodies, compared to 69.3% of 257 non-microcephaly pregnancies (P = 0.017). Analyses of sociodemographic data revealed a higher ZIKV burden in low socioeconomic status (SES) areas. High seroprevalence, combined with case data dynamics allowed estimates of the basic reproduction number R0 of 2.1 (CI, 1.8 to 2.5) at the onset of the outbreak and an effective reproductive number Reff of

Reconstructing the history of individual transmission events between cases is key to understanding what factors facilitate the spread of an infectious disease. Since conducting extended contact-tracing investigations can be logistically challenging and costly, statistical inference methods have been developed to reconstruct transmission trees from onset dates and genetic sequences. However, these methods are not as effective if the mutation rate of the virus is very slow, or if sequencing data is sparse. We developed the package o2geosocial to combine variables from routinely collected surveillance data with a simple transmission process model. The model reconstructs transmission trees when full genetic sequences are unavailable, or uninformative. Our model incorporates the reported age-group, onset date, location and genotype of infected cases to infer probabilistic transmission trees. The package also includes functions to summarise and visualise the inferred cluster size distribution. The results generated by o2geosocial can highlight regions where importations repeatedly caused large outbreaks, which may indicate a higher regional susceptibility to infections. It can also be used to generate the individual number of secondary transmissions, and show the features associated with individuals involved in high transmission events. The package is available for download from the Comprehensive R Archive Network (CRAN) and GitHub.

Reconstructing the history of individual transmission events between cases is key to understanding what factors facilitate the spread of an infectious disease. Since conducting extended contact-tracing investigations can be logistically challenging and costly, statistical inference methods have been developed to reconstruct transmission trees from onset dates and genetic sequences. However, these methods are not as effective if the mutation rate of the virus is very slow, or if sequencing data is sparse. We developed the package o2geosocial to combine variables from routinely collected surveillance data with a simple transmission process model. The model reconstructs transmission trees when full genetic sequences are not available, or uninformative. Our model incorporates the reported age-group, onset date, location and genotype of infected cases to infer probabilistic transmission trees. The package also includes functions to summarise and visualise the inferred cluster size distribution. The results generated by o2geosocial can highlight regions where importations repeatedly caused large outbreaks, which may indicate a higher regional susceptibility to infections. It can also be used to generate the individual number of secondary transmissions, and show the features associated with individuals involved in high transmission events. The package is available for download from the Comprehensive R Archive Network (CRAN) and GitHub.