Superspreading in Early Transmissions of COVID-19 in Indonesia

Introduction

The first two confirmed cases of COVID-19 in Indonesia was announced in March 2, 2020 by the president himself. Jakarta and its buffer zones including Depok have become the epicenter of the outbreak. In its early transmission, about 50% of the cases were from the cities¹. The virus quickly spread to all 34 provinces, with more than 93,000 confirmed cases and 4,500 deaths as of July 23, 2020. As a country with the highest death toll in South East Asia and one with the lowest global testing rate, the first government effort to control the disease spread was by establishing a COVID-19 response acceleration task force led by the head of the National Disaster Management Agency. The government further introduced Large-Scale Social Restrictions, which included measures such as closing public places, restricting public transport, and limiting travel to and from restricted regions. On May 19, the first milestone of 10,000 PCR tests per day was reached. In June 2020, WHO announced that only Jakarta meets the minimum requirement of 1 test per 1,000 population per week for a reliable positivity rate calculation².

As the Indonesian government still struggles to mitigate the COVID-19 epidemic in the country, it becomes critical to understand its infection spreads. Here, we present a first quantitative analysis of early transmission dynamics of the disease and evidence of superspreading events in Indonesia based on contact tracing of secondary cases. Infection clusters are difficult to predict and thus difficult to prevent. Nevertheless, understanding their environmental and behavioral drivers has been recently recognized to be extremely important to inform strategies for prevention and control³ as the reproduction number only is no longer enough⁴. In our analysis of Jakarta, we combine data from Depok that is a main buffer city of the capital. Outside the epicentre, we present data from Batam that has received thousand test kits and protective equipment from Singapore as a comparison. There were 218 cases from 4,092 PCR tests, as of June 26⁵.

Results

Based on 1288 confirmed cases in March and April 2020 in the regions Jakarta-Depok¹ (population 12.6 millions) and Batam⁵ (population 1.4 millions), we find a high degree of individual-level variation in the transmission. Early in the outbreak, superspreaders transmit the majority of coronavirus cases, as can be seen from Figure 1. Defining a superspreader as an individual with at least eight transmissions of the disease, there were 44 clusters in Jakarta-Depok region (2 - 31 March, 2020) and 3 clusters in Batam city (19 March to 7 April, 2020), as can be seen from Table 1.

• Download figure
• Open in new tab

Figure 1.

Networks of secondary cases based on local transmissions in Jakarta-Depok (left) and Batam (right) regions.

We estimated the reproduction number ℛ₀ to be at 6.79 and 2.47, while the overdispersion parameter 𝒦 of a negative-binomial distribution is obtained at 0.06 and 0.2 for Jakarta-Depok and Batam, respectively (see Figure 2), suggesting that a small number of infected individuals are responsible for large amounts of the disease transmission. Indeed, between 10-15% of all infections were responsible for 80% of onward transmission events.

• Download figure
• Open in new tab

Figure 2.

Estimates of the overdispersion parameter 𝒦 and the basic reproduction number ℛ₀. Lines in the top figures show maximum-likelihood fits for the negative binomial distribution.

The fact that Indonesia has a moderate number of positive cases, compared to other countries with roughly the same population size such as Russia, USA, and Brazil, can be caused by three reasons: a low number of tests, the success of Large- Scale Social Restrictions, and overdispersion of the COVID-19 transmission. Our result clearly indicates that the transmission is overdispersed, even though it does not exclude the other possibilities. Therefore, close-interaction activities such as traditional markets, religious gathering, and wedding parties need to be adapted if not restricted as they can become transmission hot spots. Effective response factors to reduce the transmission include aggressive implementation of non-pharmaceutical interventions such as rapid identification and isolation of cases. Furthermore, timeliness is critical to prevent or limit their extent since delay of diagnosis is the most common cause of superspreading clusters⁶.

Methods

We use a sequential Bayesian method based on the discrete SIR model to estimate the basic reproduction number ℛ₀⁷. Time series of the active cases at day n is assumed to follow the model I_n+1 = (γ (ℛ_t − 1) + 1) I_n, where ℛ_t is the instantaneous reproduction number and 1/γ is the infectious time taken to be 9 days herein⁸. Assuming a Poisson distribution for the arrival of active cases and using Bayes’ theorem, the posterior probability of ℛ_t given a confirmed case I_n+1 is computed through

The formula comes from the assumption that the prior distribution of ℛ_t at day n is taken from the posterior of day (n − 1) with a uniform prior at day 1. Estimated ℛ₀ is then the maximum of ℛ_t, see Figure 2.

The parameter 𝒦 is then defined as the overdispersion parameter in the negative binomial distribution⁹, which can be found by fitting the distribution to the secondary cases data. The number of secondary cases by one person x can be expressed in terms of the overdispersion parameter 𝒦 and the probability p,

We use a MATLAB^® function nbinfit to estimate the overdispersion parameter 𝒦. The function uses derivative-free method to find the maximum of the likelihood function.

Data and Code Availability

Secondary cases data (both in Excel and Power BI) and a MATLAB^® code to produce Figure 1 and Figure 2 are available at: https://github.com/agusisma/covidindonesia.

Author contributions statement

A.H., H.S., and M.F.K. analysed the data, calculated ℛ₀ and 𝒦, and discussed the results, N.N., B.L., D.T., and W.W. collected the data. All authors reviewed the manuscript.