COVID-19 in Japan: insights from the first three months of the epidemic

BackgroundUnderstanding the characteristics and natural history of novel pathogens is crucial to inform successful control measures. Japan was one of the first affected countries in the COVID-19 pandemic reporting their first case on 14 January 2020. Interventions including airport screening, contact tracing, and cluster investigations were quickly implemented. Here we present insights from the first 3 months of the epidemic in Japan based on detailed case data. MethodsWe conducted descriptive analyses based on information systematically extracted from individual case reports from 13 January to 31 March 2020 including patient demographics, date of report and symptom onset, symptom progression, travel history, and contact type. We analysed symptom progression and estimated the time-varying reproduction number, Rt, correcting for epidemic growth using an established Bayesian framework. Key delays and the age-specific probability of transmission were estimated using data on exposures and transmission pairs. ResultsThe corrected fitted mean onset-to-reporting delay after the peak was 4 days (standard deviation: {+/-}2 days). Early transmission was driven primarily by returning travellers with Rt peaking at 2.4 (95%CrI:1.6, 3.3) nationally. In the final week of the trusted period, Rt accounting for importations diverged from overall Rt at 1.1 (95% CrI: 1.0, 1.2) compared to 1.5 (95% CrI: 1.3, 1.6) respectively. Household (39.0%) and workplace (11.6%) exposures were the most frequently reported potential source of infection. The estimated probability of transmission was assortative by age. Across all age groups, cases most frequently onset with cough, fever, and fatigue. There were no reported cases of patients <20 years old developing pneumonia or severe respiratory symptoms. ConclusionsInformation collected in the early phases of an outbreak are important in characterising any novel pathogen. Timely recognition of key symptoms and high-risk settings for transmission can help to inform response strategies. The data analysed here were the result of robust and timely investigations and demonstrate the improvements to epidemic control as a result of such surveillance.

A generic method and software to estimate the transmission advantage of pathogen variants in real-time : SARS-CoV-2 as a case-study

Recent months have demonstrated that emerging variants may set back the global COVID-19 response. The ability to rapidly assess the threat of new variants in real-time is critical for timely optimisation of control strategies. We extend the EpiEstim R package, designed to estimate the time-varying reproduction number (Rt), to estimate in real-time the effective transmission advantage of a new variant compared to a reference variant. Our method can combine information across multiple locations and over time and was validated using an extensive simulation study, designed to mimic a variety of real-time epidemic contexts. We estimate that the SARS-CoV-2 Alpha variant is 1.46 (95% Credible Interval 1.44-1.47) and 1.29, (95% CrI 1.29-1.30) times more transmissible than the wild type, using data from England and France respectively. We further estimate that Beta and Gamma combined are 1.25 (95% CrI 1.24-1.27) times more transmissible than the wildtype (France data). All results are in line with previous estimates from literature, but could have been obtained earlier and more easily with our off-the-shelf open-source tool. Our tool can be used as an important first step towards quantifying the threat of new variants in real-time. Given the popularity of EpiEstim, this extension will likely be used widely to monitor the co-circulation and/or emergence of multiple variants of infectious pathogens. Significance StatementEarly assessment of the transmissibility of new variants of an infectious pathogen is critical for anticipating their impact and designing appropriate interventions. However, this often requires complex and bespoke analyses relying on multiple data streams, including genomic data. Here we present a novel method and software to rapidly quantify the transmission advantage of new variants. Our method is fast and requires only routinely collected disease surveillance data, making it easy to use in real-time. The ongoing high level of SARS-CoV-2 circulation in a number of countries makes the emergence of new variants highly likely. Our work offers a powerful tool to help public health bodies monitor such emerging variants and rapidly detect those with increased transmissibility.

National and subnational short-term forecasting of COVID-19 in Germany and Poland, early 2021

BackgroundDuring the COVID-19 pandemic there has been a strong interest in forecasts of the short-term development of epidemiological indicators to inform decision makers. In this study we evaluate probabilistic real-time predictions of confirmed cases and deaths from COVID-19 in Germany and Poland for the period from January through April 2021. MethodsWe evaluate probabilistic real-time predictions of confirmed cases and deaths from COVID-19 in Germany and Poland. These were issued by 15 different forecasting models, run by independent research teams. Moreover, we study the performance of combined ensemble forecasts. Evaluation of probabilistic forecasts is based on proper scoring rules, along with interval coverage proportions to assess forecast calibration. The presented work is part of a pre-registered evaluation study and covers the period from January through April 2021. ResultsWe find that many, though not all, models outperform a simple baseline model up to four weeks ahead for the considered targets. Ensemble methods (i.e., combinations of different available forecasts) show very good relative performance. The addressed time period is characterized by rather stable non-pharmaceutical interventions in both countries, making short-term predictions more straightforward than in previous periods. However, major trend changes in reported cases, like the rebound in cases due to the rise of the B.1.1.7 (alpha) variant in March 2021, prove challenging to predict. ConclusionsMulti-model approaches can help to improve the performance of epidemiological forecasts. However, while death numbers can be predicted with some success based on current case and hospitalization data, predictability of case numbers remains low beyond quite short time horizons. Additional data sources including sequencing and mobility data, which were not extensively used in the present study, may help to improve performance. Plain language summaryThe goal of this study is to assess the quality of forecasts of weekly case and death numbers of COVID-19 in Germany and Poland during the period of January through April 2021. We focus on real-time forecasts at time horizons of one and two weeks ahead created by fourteen independent teams. Forecasts are systematically evaluated taking uncertainty ranges of predictions into account. We find that combining different forecasts into ensembles can improve the quality of predictions, but especially case numbers proved very challenging to predict beyond quite short time windows. Additional data sources, in particular genetic sequencing data, may help to improve forecasts in the future.

Global predictions of short- to medium-term COVID-19 transmission trends : a retrospective assessment

BackgroundAs of July 2021, more than 180,000,000 cases of COVID-19 have been reported across the world, with more than 4 million deaths. Mathematical modelling and forecasting efforts have been widely used to inform policy-making and to create situational awareness. Methods and FindingsFrom 8th March to 29th November 2020, we produced weekly estimates of SARS-CoV-2 transmissibility and forecasts of deaths due to COVID-19 for countries with evidence of sustained transmission. The estimates and forecasts were based on an ensemble model comprising of three models that were calibrated using only the reported number of COVID-19 cases and deaths in each country. We also developed a novel heuristic to combine weekly estimates of transmissibility and potential changes in population immunity due to infection to produce forecasts over a 4-week horizon. We evaluated the robustness of the forecasts using relative error, coverage probability, and comparisons with null models. ConclusionsDuring the 39-week period covered by this study, we produced short- and medium-term forecasts for 81 countries. Both the short- and medium-term forecasts captured well the epidemic trajectory across different waves of COVID-19 infections with small relative errors over the forecast horizon. The model was well calibrated with 56.3% and 45.6% of the observations lying in the 50% Credible Interval in 1-week and 4-week ahead forecasts respectively. We could accurately characterise the overall phase of the epidemic up to 4-weeks ahead in 84.9% of country-days. The medium-term forecasts can be used in conjunction with the short-term forecasts of COVID-19 mortality as a useful planning tool as countries continue to relax stringent public health measures that were implemented to contain the pandemic.

The 2020 SARS-CoV-2 epidemic in England: key epidemiological drivers and impact of interventions

We fitted a model of SARS-CoV-2 transmission in care homes and the community to regional surveillance data for England. Among control measures implemented, only national lockdown brought the reproduction number below 1 consistently; introduced one week earlier it could have reduced first wave deaths from 36,700 to 15,700 (95%CrI: 8,900-26,800). Improved clinical care reduced the infection fatality ratio from 1.25% (95%CrI: 1.18%-1.33%) to 0.77% (95%CrI: 0.71%-0.84%). The infection fatality ratio was higher in the elderly residing in care homes (35.9%, 95%CrI: 29.1%-43.4%) than those residing in the community (10.4%, 95%CrI: 9.1%-11.5%). England is still far from herd immunity, with regional cumulative infection incidence to 1st December 2020 between 4.8% (95%CrI: 4.4%-5.1%) and 15.4% (95%CrI: 14.9%-15.9%) of the population. One-sentence summaryWe fit a mathematical model of SARS-CoV-2 transmission to surveillance data from England, to estimate transmissibility, severity, and the impact of interventions

Short-term forecasting of COVID-19 in Germany and Poland during the second wave - a preregistered study

We report insights from ten weeks of collaborative COVID-19 forecasting for Germany and Poland (12 October - 19 December 2020). The study period covers the onset of the second wave in both countries, with tightening non-pharmaceutical interventions (NPIs) and subsequently a decay (Poland) or plateau and renewed increase (Germany) in reported cases. Thirteen independent teams provided probabilistic real-time forecasts of COVID-19 cases and deaths. These were reported for lead times of one to four weeks, with evaluation focused on one- and two-week horizons, which are less affected by changing NPIs. Heterogeneity between forecasts was considerable both in terms of point predictions and forecast spread. Ensemble forecasts showed good relative performance, in particular in terms of coverage, but did not clearly dominate single-model predictions. The study was preregistered and will be followed up in future phases of the pandemic.

Quantifying the dynamics of COVID-19 burden and impact of interventions in Java, Indonesia

BackgroundAs in many countries, quantifying COVID-19 spread in Indonesia remains challenging due to testing limitations. In Java, non-pharmaceutical interventions (NPIs) were implemented throughout 2020. However, as a vaccination campaign launches, cases and deaths are rising across the island. MethodsWe used modelling to explore the extent to which data on burials in Jakarta using strict COVID-19 protocols (C19P) provide additional insight into the transmissibility of the disease, epidemic trajectory, and the impact of NPIs. We assess how implementation of NPIs in early 2021 will shape the epidemic during the period of likely vaccine roll-out. ResultsC19P burial data in Jakarta suggest a death toll approximately 3.3 times higher than reported. Transmission estimates using these data suggest earlier, larger, and more sustained impact of NPIs. Measures to reduce sub-national spread, particularly during Ramadan, substantially mitigated spread to more vulnerable rural areas. Given current trajectory, daily cases and deaths are likely to increase in most regions as the vaccine is rolled-out. Transmission may peak in early 2021 in Jakarta if current levels of control are maintained. However, relaxation of control measures is likely to lead to a subsequent resurgence in the absence of an effective vaccination campaign. ConclusionSyndromic measures of mortality provide a more complete picture of COVID-19 severity upon which to base decision-making. The high potential impact of the vaccine in Java is attributable to reductions in transmission to date and dependent on these being maintained. Increases in control in the relatively short-term will likely yield large, synergistic increases in vaccine impact. Key questionsO_ST_ABSWhat is already known?C_ST_ABSO_LIIn many settings, limited SARS-CoV-2 testing makes it difficult to estimate the true trajectory and associated burden of the virus. C_LIO_LINon-pharmaceutical interventions (NPIs) are key tools to mitigate SARS-CoV-2 transmission. C_LIO_LIVaccines show promise but effectiveness depends upon prioritization strategies, roll-out and uptake. C_LI What are the new findings?O_LIThis study gives evidence of the value of syndrome-based mortality as a metric, which is less dependent upon testing capacity with which to estimate transmission trends and evaluate intervention impact. C_LIO_LINPIs implemented in Java earlier in the pandemic have substantially slowed the course of the epidemic with movement restrictions during Ramadan preventing spread to more vulnerable rural populations. C_LIO_LIPopulation-level immunity remains below proposed herd-immunity thresholds for the virus, though it is likely substantially higher in Jakarta. C_LI What do the new findings imply?O_LIGiven current levels of control, upwards trends in deaths are likely to continue in many provinces while the vaccine is scheduled to be rolled out. A key exception is Jakarta where population-level immunity may increase to a level where the epidemic begins to decline before the vaccine campaign has reached high coverage. C_LIO_LIFurther relaxation of measures would lead to more rapidly progressing epidemics, depleting the eventual incremental effectiveness of the vaccine. Maintaining adherence to control measures in Jakarta may be particularly challenging if the epidemic enters a decline phase but will remain necessary to prevent a subsequent large wave. Elsewhere, higher levels of control with NPIs are likely to yield high synergistic vaccine impact. C_LI

A sub-national analysis of the rate of transmission of COVID-19 in Italy

Italy was the first European country to experience sustained local transmission of COVID-19. As of 1st May 2020, the Italian health authorities reported 28,238 deaths nationally. To control the epidemic, the Italian government implemented a suite of non-pharmaceutical interventions (NPIs), including school and university closures, social distancing and full lockdown involving banning of public gatherings and non essential movement. In this report, we model the effect of NPIs on transmission using data on average mobility. We estimate that the average reproduction number (a measure of transmission intensity) is currently below one for all Italian regions, and significantly so for the majority of the regions. Despite the large number of deaths, the proportion of population that has been infected by SARS-CoV-2 (the attack rate) is far from the herd immunity threshold in all Italian regions, with the highest attack rate observed in Lombardy (13.18% [10.66%-16.70%]). Italy is set to relax the currently implemented NPIs from 4th May 2020. Given the control achieved by NPIs, we consider three scenarios for the next 8 weeks: a scenario in which mobility remains the same as during the lockdown, a scenario in which mobility returns to pre-lockdown levels by 20%, and a scenario in which mobility returns to pre-lockdown levels by 40%. The scenarios explored assume that mobility is scaled evenly across all dimensions, that behaviour stays the same as before NPIs were implemented, that no pharmaceutical interventions are introduced, and it does not include transmission reduction from contact tracing, testing and the isolation of confirmed or suspected cases. New interventions, such as enhanced testing and contact tracing are going to be introduced and will likely contribute to reductions in transmission; therefore our estimates should be viewed as pessimistic projections. We find that, in the absence of additional interventions, even a 20% return to pre-lockdown mobility could lead to a resurgence in the number of deaths far greater than experienced in the current wave in several regions. Future increases in the number of deaths will lag behind the increase in transmission intensity and so a second wave will not be immediately apparent from just monitoring of the daily number of deaths. Our results suggest that SARS-CoV-2 transmission as well as mobility should be closely monitored in the next weeks and months. To compensate for the increase in mobility that will occur due to the relaxation of the currently implemented NPIs, adherence to the recommended social distancing measures alongside enhanced community surveillance including swab testing, contact tracing and the early isolation of infections are of paramount importance to reduce the risk of resurgence in transmission. SUGGESTED CITATIONMichaela A. C. Vollmer, Swapnil Mishra, H Juliette T Unwin, Axel Gandy et al. Using mobility to estimate the transmission intensity of COVID-19 in Italy: a subnational analysis with future scenarios. Imperial College London (2020) doi:https://doi.org/10.25561/78677 This work is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License.

Report 21: Estimating COVID-19 cases and reproduction number in Brazil

1Brazil is currently reporting the second highest number of COVID-19 deaths in the world. Here we characterise the initial dynamics of COVID-19 across the country and assess the impact of non-pharmaceutical interventions (NPIs) that were implemented using a semi-mechanistic Bayesian hierarchical modelling approach. Our results highlight the significant impact these NPIs had across states, reducing an average Rt > 3 to an average of 1.5 by 9-May-2020, but that these interventions failed to reduce Rt < 1, congruent with the worsening epidemic Brazil has experienced since. We identify extensive heterogeneity in the epidemic trajectory across Brazil, with the estimated number of days to reach 0.1% of the state population infected since the first nationally recorded case ranging from 20 days in Sao Paulo compared to 60 days in Goias, underscoring the importance of sub-national analyses in understanding asynchronous state-level epidemics underlying the national spread and burden of COVID-19.

Estimates of the severity of COVID-19 disease

BackgroundA range of case fatality ratio (CFR) estimates for COVID-19 have been produced that differ substantially in magnitude. MethodsWe used individual-case data from mainland China and cases detected outside mainland China to estimate the time between onset of symptoms and outcome (death or discharge from hospital). We next obtained age-stratified estimates of the CFR by relating the aggregate distribution of cases by dates of onset to the observed cumulative deaths in China, assuming a constant attack rate by age and adjusting for the demography of the population, and age- and location-based under-ascertainment. We additionally estimated the CFR from individual line-list data on 1,334 cases identified outside mainland China. We used data on the PCR prevalence in international residents repatriated from China at the end of January 2020 to obtain age-stratified estimates of the infection fatality ratio (IFR). Using data on age-stratified severity in a subset of 3,665 cases from China, we estimated the proportion of infections that will likely require hospitalisation. FindingsWe estimate the mean duration from onset-of-symptoms to death to be 17.8 days (95% credible interval, crI 16.9-19.2 days) and from onset-of-symptoms to hospital discharge to be 22.6 days (95% crI 21.1-24.4 days). We estimate a crude CFR of 3.67% (95% crI 3.56%-3.80%) in cases from mainland China. Adjusting for demography and under-ascertainment of milder cases in Wuhan relative to the rest of China, we obtain a best estimate of the CFR in China of 1.38% (95% crI 1.23%-1.53%) with substantially higher values in older ages. Our estimate of the CFR from international cases stratified by age (under 60 / 60 and above) are consistent with these estimates from China. We obtain an overall IFR estimate for China of 0.66% (0.39%-1.33%), again with an increasing profile with age. InterpretationThese early estimates give an indication of the fatality ratio across the spectrum of COVID-19 disease and demonstrate a strong age-gradient in risk.