Skip to main content

Real-time monitoring the transmission potential of COVID-19 in Singapore, March 2020

Abstract

Background

As of March 31, 2020, the ongoing COVID-19 epidemic that started in China in December 2019 is now generating local transmission around the world. The geographic heterogeneity and associated intervention strategies highlight the need to monitor in real time the transmission potential of COVID-19. Singapore provides a unique case example for monitoring transmission, as there have been multiple disease clusters, yet transmission remains relatively continued.

Methods

Here we estimate the effective reproduction number, Rt, of COVID-19 in Singapore from the publicly available daily case series of imported and autochthonous cases by date of symptoms onset, after adjusting the local cases for reporting delays as of March 17, 2020. We also derive the reproduction number from the distribution of cluster sizes using a branching process analysis that accounts for truncation of case counts.

Results

The local incidence curve displays sub-exponential growth dynamics, with the reproduction number following a declining trend and reaching an estimate at 0.7 (95% CI 0.3, 1.0) during the first transmission wave by February 14, 2020, while the overall R based on the cluster size distribution as of March 17, 2020, was estimated at 0.6 (95% CI 0.4, 1.02). The overall mean reporting delay was estimated at 6.4 days (95% CI 5.8, 6.9), but it was shorter among imported cases compared to local cases (mean 4.3 vs. 7.6 days, Wilcoxon test, p < 0.001).

Conclusion

The trajectory of the reproduction number in Singapore underscores the significant effects of successful containment efforts in Singapore, but it also suggests the need to sustain social distancing and active case finding efforts to stomp out all active chains of transmission.

Peer Review reports

Background

The ongoing Coronavirus Disease 2019 (COVID-19) pandemic started with a cluster of pneumonia cases of unknown etiology in Wuhan, China, back in December 2019 [1, 2]. The initial cases have been linked to a wet market in the city of Wuhan, pointing to an animal source of the epidemic [3]. Subsequently, rapid human-to-human transmission of the disease was confirmed in January 2020, and the etiological agent was identified as severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) due to its genetic similarity to the SARS-CoV (severe acute respiratory syndrome coronavirus) discovered in 2003 [4, 5]. The total global case tally has reached 750,890 infections including 36,405 deaths and involving 199 countries as of March 31, 2020 [6]. As the virus continues to spread in the human population, obtaining an accurate “real-time” picture of the epidemic’s trajectory is complicated by several factors including reporting delays and changes in the case definition [7, 8]. Although the COVID-19 case incidence in China has substantially declined, active transmission is now occurring in multiple countries around the world [2]. Epidemiological data from these countries can help to monitor transmission potential of SARS-CoV-2 in near real time.

Outside of China, Singapore, where the first symptomatic imported case (66 years old Chinese male) was reported on January 23, 2020, has been able to maintain relatively low COVID-19 incidence levels through active case finding and strict social distancing measures. Up until March 31, 2020, Singapore has reported 926 laboratory confirmed cases, including 24 reported case importations from Wuhan China and 501 non-Wuhan-related case importations [9]. Imported cases include six individuals who were evacuated from China between January 30 and February 9, 2020, and multiple citizens and long-term Singapore pass holders returning from Asia, Europe, and North America in late March 2020 [9,10,11]. Moreover, Singapore has reported 3 deaths as of March 31, 2020 [9]. On February 4, 2020, the Ministry of Health of Singapore reported its first local cluster of COVID-19, which was linked to the Yong Thai Hang shop [12]. A total of 18 clusters with 2 or more COVID-19 cases have been reported thus far. Table 1 summarizes the characteristics of the 6 largest clusters in Singapore.

Table 1 Characteristics of the largest COVID-19 outbreak in Singapore as of March 17, 2020

Although large-scale community transmission has not been reported in Singapore, the novel coronavirus can rapidly spread in confined and crowded places, as illustrated by large clusters of COVID-19 cases linked to the Grace Assembly of God Church, the Life Church and Missions Singapore, Wizlearn Technologies, and the SAFRA Jurong cluster [13, 14]. In China, substantial hospital-based transmission of SARS-CoV-2 has been reported, with approximately 3400 cases involving healthcare workers [15]. This pattern aligns well with past outbreaks of severe acute respiratory syndrome (SARS) and Middle Eastern respiratory syndrome (MERS) [16], including substantial nosocomial transmission during the 2003 SARS outbreak in Singapore [17]. To minimize the risk of hospital-based transmission of SARS-CoV-2, the Ministry of Health of Singapore has restricted the movement of patients and staff across hospitals [18]. Also, because multiple unlinked COVID-19 cases have been reported in the community [19] and the recognition that a substantial proportion of asymptomatic cases may be spreading the virus [20,21,22], strict social distancing measures have been put in place including advising the public against large social gatherings in order to mitigate the risk of community transmission [23]. These social distancing measures reduce the risk of onward transmission not only within Singapore, but also beyond the borders of this highly connected nation [24]. A recent influx of imported cases from Asia, Europe, and North America into Singapore has triggered travel bans and restrictions for travelers and citizens [9].

The reproduction number is a key threshold quantity to assess the transmission potential of an emerging disease such as COVID-19 [25, 26]. It quantifies the average number of secondary cases generated per case. If the reproduction number is below 1.0, infections occur in isolated clusters as self-limited chains of transmission, and persistence of the disease would require continued undetected importations. On the other hand, reproduction numbers above 1.0 indicate sustained community transmission [16, 26]. Using epidemiological data and mathematical modeling tools, we are monitoring the effective reproduction number, Rt, of SARS-CoV-2 transmission in Singapore in real-time, and here we report the evolution of Rt by March 17, 2020. Specifically, we characterize the growth profile and the effective reproduction number during the first transmission wave from the daily case series of imported and autochthonous cases by date of symptoms onset after adjusting for reporting delays, and we also derive an estimate of the overall reproduction number based on the characteristics of the clusters of COVID-19 in Singapore.

Methods

Data

We obtained the daily series of 247 confirmed COVID-19 cases in Singapore between January 23 and March 17, 2020, from public records of the Ministry of Health, Singapore, as of March 17, 2020 [27]. Individual-level case details including the dates of symptom onset, the date of reporting, and whether the case is autochthonous (local transmission) or imported are publicly available. Clusters consisting of two or more cases according to the infection source were also assembled from case descriptions obtained from field investigations conducted by the Ministry of Health, Singapore [27]. Single imported cases are analyzed as clusters of size 1 whereas unlinked cases were excluded from the cluster analysis.

Transmission clusters

As of March 17, 2020, 18 different clusters of COVID-19 cases with 2–48 cases per cluster have been reported in Singapore. A schematic diagram and characteristics of the COVID-19 clusters in Singapore are given in Fig. 1 and Table 1. The geographic location of the six clusters accounting for 45.3% of the total cases is shown in Fig. 2 whereas the corresponding distribution of cluster sizes is shown in Fig. 3.

Fig. 1
figure1

Cluster network of the cases in Singapore for the COVID-19 global pandemic as of March 17, 2020. The pink circles represent the cases linked to Wuhan, the green circles represent the non-Wuhan-related case importations, and the blue circles represent cases with no travel history to China. The larger dotted circles represent the COVID-19 disease clusters. Each blue arrow represents the direction in which the disease was transmitted. Pink arrows represent immediate family. Dates below the circles are the dates of case reporting

Fig. 2
figure2

Map depicting the spatial distribution of the 6 largest COVID-19 clusters in Singapore; Grand Hyatt cluster, Yong Thai Hang cluster, Seletar Aerospace cluster, Wizlearn Technologies cluster, SAFRA Jurong cluster, and The Grace Assembly of God Church and Life Church and Missions cluster as of March 17, 2020

Fig. 3
figure3

Distribution of COVID-19 cluster sizes in Singapore as of March 17, 2020

Yong Thai Hang cluster

This cluster with 9 cases was the first to be reported in Singapore. It has nine traceable links, including eight Chinese and one Indonesian national associated with the visit of Chinese tourists to the Yong Thai Hang health products store, a shop that primarily serves the Chinese population, on January 23, 2020. Four shop employees and the tour guide were first identified as a cluster on February 4, 2020 [12, 28, 29]. The tour guide subsequently infected her husband, a newborn, and the domestic helper [29]. No further cases have been added to this cluster as of February 8, 2020.

Grand Hyatt hotel

This cluster with 3 local cases was the second cluster to receive international attention, as it originated from a business meeting held at the Grand Hyatt hotel attended by Singaporean locals and the Chinese visitors from Hubei [30]. Four international cases associated with this cluster had left Singapore before the onset of symptoms. All Singaporean residents associated with this cluster have recovered as of February 19, 2020 [30]. No additional cases have been added to this cluster as of February 8, 2020.

Seletar Aerospace Heights cluster

This cluster with 5 Bangladeshi work pass holders was identified on February 9, 2020. No further cases have been added to this cluster as of February 15, 2020.

The Life Church and Missions and The Grace Assembly of God cluster

This Singaporean cluster is composed of 33 cases, including two imported cases and 31 local cases. The cluster started during The Life Church and Missions service event in Paya Lebar on January 19, 2020. This event was apparently seeded by two visitors from Wuhan China who infected a couple with SARS-CoV-2 at the church. The infected couple likely passed the infection to another case during a Lunar New Year’s celebration on January 25, 2020. This case had subsequently infected Grace Assembly of God church staff at the Tanglin branch, generating secondary cases by the time he was reported on February 14, 2020. Two branches of the Grace Assembly of God church at Tanglin and Bukit Batok have been included in this cluster [28, 31]. This church serves an average of 4800 people in attendance over the weekend. While the church has momentarily closed, field investigations have not led to conclusive evidence regarding superspreading transmission. No further cases have been added to this cluster as of March 9, 2020.

SAFRA Jurong cluster

The largest cluster composed of 48 local cases is linked to a private dinner function at SAFRA Jurong restaurant on February 15, 2020. The restaurant was closed for cleaning from February 16 to February 19, 2020, following the dinner function. The latest case was added to this cluster on March 16, 2020.

Wizlearn Technologies cluster

This cluster which comprises of 14 cases was identified on February 26, 2020. Wizlearn Technologies is an e-learning solutions company. The latest case was added to this cluster on March 3, 2020.

Church of Singapore cluster

The first case of this cluster was identified on March 14, 2020, originating as a secondary case from a case in the SAFRA Jurong cluster. This cluster is composed of 3 local cases. No further cases have been added to this cluster since March 16, 2020.

Boulder gym cluster

The first case of this cluster was identified on March 8, 2020, also linked to the SAFRA Jurong cluster. This cluster is composed of 3 local cases. No further cases have been added to this cluster since March 10, 2020.

Cluster A

The first case of this cluster was identified on February 14, 2020. The cluster comprises of 3 local cases. No further cases have been added in this cluster as of February 18, 2020.

Cluster B

The first case of this cluster was identified on February 19, 2020. This cluster is composed of two local cases. No further cases have been added in this cluster since February 21, 2020.

Cluster C

This first case of this cluster was identified on March 3, 2020. This cluster is composed of 2 local cases. No further cases have been added to this cluster since March 6, 2020.

Cluster D

The first case of this cluster was identified on March 7, 2020. The two cases (one imported and one local) in this cluster are related to each other. No further cases have been added in this cluster since March 8, 2020.

Cluster E

The two cases (an imported and a local case) of this cluster were identified on March 11, 2020. No further cases have been added in this cluster since March 11, 2020.

Cluster F

The first case of this cluster was identified on March 11, 2020. This cluster is composed of 3 local cases. No further cases have been added to this cluster since March 13, 2020.

Cluster G

This cluster is composed of 5 cases, including 4 imported cases. The first case of this cluster was identified on March 10, 2020. No further cases have been added to this cluster since March 13, 2020.

Cluster H

The first case of this cluster was identified on March 14, 2020, a secondary case generated from a case at SAFRA Jurong cluster. This cluster is composed of 3 local cases. No further cases have been added to this cluster since March 16, 2020.

Cluster I

The first case of this cluster was identified on March 14, 2020. This cluster is composed of one local and one imported case. No further cases have been added to this cluster since March 15, 2020.

Cluster J

The first case of this cluster was identified on March 15, 2020. This cluster is composed of one imported and one local case. No further cases have been added to this cluster since March 16, 2020.

Adjusting for reporting delays

As an outbreak progresses in real time, epidemiological curves can be distorted by reporting delays arising from several factors that include (i) delays in case detection during field investigations, (ii) delays in symptom onset after infection, (iii) delays in seeking medical care, (iv) delays in diagnostics, and (v) delays in processing data in surveillance systems [32]. However, it is possible to generate reporting-delay-adjusted incidence curves using standard statistical methods [33]. Briefly, the reporting delay for a case is defined as the time lag in days between the date of onset and date of reporting. Here we adjusted the COVID-19 epidemic curve of local cases by reporting delays using a non-parametric method that employs survival analysis known as the Actuaries method for use with right truncated data, employing reverse time hazards to adjust for reporting delays as described in a previous publication [34,35,36]. The 95% prediction limits are derived according to Lawless and Kalbfleisch [37]. For this analysis, we exclude 7 imported cases and 5 local cases for which dates of symptoms onset are unavailable.

Effective reproduction number from case incidence

We assess the effective reproduction number over the course of the outbreak, Rt, which quantifies the temporal variation in the average number of secondary cases generated per case during the course of an outbreak after considering multiple factors including behavior changes, cultural factors, and the implementation of public health measures [16, 26, 38]. Estimates of Rt > 1 indicate sustained transmission, whereas, Rt < 1 implies that the outbreak is slowing down and the incidence trend is declining. Hence, maintaining Rt < 1 is required to bring an outbreak under control. Using the reporting delay adjusted incidence curve, we estimate the most recent estimate of Rt for COVID-19 in Singapore by characterizing the early transmission phase using a phenomenological growth model as described in previous publications [39,40,41,42]. Specifically, we first characterize daily incidence of local cases for the first transmission wave (January 21–February 14, 2020) using the generalized logistic growth model (GLM) after adjusting for imported cases. This model characterizes the growth profile via three parameters: the growth rate (r), the scaling of the growth parameter (p), and the final epidemic size (K). The GLM can reproduce a range of early growth dynamics, including constant growth (p = 0), sub-exponential or polynomial growth (0 < p < 1), and exponential growth (p = 1) [40, 42]. We denote the local incidence at calendar time ti by Ii, the raw incidence of imported cases at calendar time ti by Ji, and the discretized probability distribution of the generation interval by ρi. The generation interval is assumed to follow a gamma distribution with a mean of 4.41 days and a standard deviation of 3.17 days based on refs. [43, 44]. Then, we can estimate the effective reproduction number by employing the renewal equation given by [45, 46]

$$ {R}_{t_i}=\frac{I_i}{\sum_{j=0}^i\left({I}_{i-j}+\alpha\ {J}_{i-j}\right){\rho}_i} $$

In this equation, the numerator represents the new cases Ii, and the denominator represents the total number of cases that contribute to the new cases Ii at time ti. Parameter 0 ≤α ≤ 1 represents the relative contribution of imported cases to the secondary disease transmission. We perform a sensitivity analyses by setting α = 0.15 and α = 1.0 [47]. Next, in order to derive the uncertainty bounds around the curve of Rt directly from the uncertainty associated with the parameters estimates (r, p, K), we estimate Rt for 300 simulated curves assuming a Poisson error structure [48].

Reproduction number (R) from the analysis of cluster sizes

A second method of inferring the reproduction number applies branching process theory to cluster size data to infer the degree of transmission heterogeneity [49, 50]. Simultaneous inference of heterogeneity and the reproduction number has been shown to improve the reliability of confidence intervals for the reproduction number [51]. In the branching process analysis, the number of transmissions caused by each new infection is modeled as a negative binomial distribution. This is parameterized by the effective reproduction number, R, and the dispersion parameter, k. The reproduction number provides the average number of secondary cases per index case, and the dispersion parameter varies inversely with the heterogeneity of the infectious disease. In this parameterization, a lower dispersion parameter indicates higher transmission heterogeneity.

Branching process theory provides an analytic representation of the size distribution of cluster sizes as a function of R, k, and the number of primary infections in a cluster (as represented in equation of 6 of the supplement of [52]). This permits direct inference of the maximum likelihood estimate and confidence interval for R and k. In this manuscript, we modify the calculation of the likelihood of a cluster size to account for the possibility that truncation of case counts at a specific time point (i.e., March 17, 2020) may result in some infections being unobserved. This is accomplished by denoting x as the sum of the observed number of serial intervals in a cluster. Then the likelihood that an observed cluster of size j containing m imported cases is generated by x infectious intervals is given by:

$$ {l}_{m\to j}^C\left({R}_{\mathrm{eff}},k,x\right)=\frac{m}{j}{l}_{x\to \left(j-m\right)}\left({R}_{\mathrm{eff}},k\right) $$
(1)

where the likelihood of i infections causing j infections is given by:

$$ {l}_{i\to j}\left({R}_{\mathrm{eff}},k\right)=\frac{\Gamma \left(j+ ki\right)}{\Gamma \left(j+1\right)\Gamma (ki)}\ {\left(\frac{k}{R_{\mathrm{eff}}+k}\right)}^{k\mathrm{i}}{\left(\frac{R_{\mathrm{eff}}}{R_{\mathrm{eff}}+k}\right)}^j $$
(2)

where Γ is the gamma function.

To determine the number of observed serial intervals in each cluster, we first estimate the cumulative probability distribution of the serial interval. We assume the serial interval is a gamma distribution, with a mean of 4.7 days and a standard deviation of 2.9 days [43]. This translates to a shape parameter of 2.63 and a scale parameter of 1.79. We then use the difference between the onset data and the end of our study (March 17, 2020) to determine how much of the infectious period was observed. For cases that only have a report date, but no onset date, we assume an onset date that is 6 days earlier than the reporting date. This is based on the average duration between onset date and report date that was observed in the data. When applied to the case series, we are able to assign a total size, the number of imported cases, and the observed number of infectious periods for each cluster in the case series. When no imported cases are known to be in a cluster, we assign the number of imported cases to be one as the cluster must have been initiated by someone (e.g., the index case had contact with a foreign visitor).

When Eq. (1) is applied to the table of cluster size characteristics, the likelihood of the data can be calculated as a function of R and k. Minimizing the likelihood produces the maximum likelihood estimates of R and k. Applying the likelihood ratio test by profiling and R and k produces confidence intervals [53]. Code was run in R version 3.6.1.

Results

Incidence data and reporting delays

The COVID-19 epidemic curve by the date of reporting, stratified for local and imported incidence case counts, is shown in Fig. 4. It shows that the majority of the imported cases are concentrated at the beginning of the outbreak (January 23, 2020 to February 3, 2020) and after March 10, 2020, in Singapore, with an average of ~12 new cases reported per day between March 1, 2020, and March 17, 2020 (Fig. 4). Out of 88 imported cases, only 14 cases have been linked to secondary cases. Meanwhile, a total of 159 autochthonous cases have been reported as of March 17, 2020, including 27 cases that are unlinked to any known transmission chains.

Fig. 4
figure4

Local and imported incidence cases by date of reporting as of March 17, 2020. The solid blue line represents the cumulative case count for the COVID-19 cases in Singapore

The reporting-delay adjusted epidemic curve of local cases by date of symptoms onset roughly displays two small waves of transmission reflecting the occurrence of asynchronous case clusters (Fig. 5). Moreover, the gamma distribution provided a reasonable fit to the distribution of reporting delays for all cases, with a mean reporting delay at 6.4 days (95% CI 5.8, 6.9) (Fig. 6). We also found that imported cases tend to have shorter reporting delays compared to local cases (mean 4.3 vs. 7.6 days, Wilcoxon test, p < 0.001), as imported cases tend to be identified more quickly. The mean of reporting delays for the six large clusters ranged from 4.8 to 13.6 days (Fig. 7).

Fig. 5
figure5

Reporting delay adjusted local incidence for the COVID-19 outbreak in Singapore as of March 17, 2020. Blue bars represent the raw incidence, red solid line represents the adjusted incidence, and red dotted lines represent the 95% lower and upper bound of the adjusted incidence

Fig. 6
figure6

The distribution of reporting delays for all cases as of March 17, 2020. The red line represents the fit of a gamma distribution to the data. The red circle represents the mean of gamma distribution and the horizontal line represents the 95% CI

Fig. 7
figure7

Reporting delay distribution with mean (blue circle) and 95% CI (vertical lines) for each big cluster in Singapore; Grand Hyatt cluster, Yong Thai Hang cluster, Seletar Aerospace cluster, Wizlearn Technologies cluster, SAFRA Jurong cluster, and The Grace Assembly of God Church and Life Church and Missions cluster as of March 17, 2020

Reproduction numbers

For the first small wave of transmission comprising the first 25 epidemic days, the delay-adjusted local incidence curve displays sub-exponential growth dynamics with the scaling of growth parameter p at 0.7 (95% CI 0.4, 1.0), the intrinsic growth rate r estimated at 0.6 (95% CI 0.3, 1.1) and parameter K estimating the wave size estimated at 95 (95% CI 56, 230). Because of the sub-exponential growth dynamics, the effective reproduction number followed a declining trend with the latest estimate at 0.7 (95% CI 0.3, 1.0) when α = 0.15 (Fig. 8). This estimate was not sensitive to changes in parameter α.

Fig. 8
figure8

The effective reproduction number reproduction number with 95% CI estimated by adjusting for the imported cases α = 0.15 during the first transmission wave by February 14, 2020. The effective reproduction number followed a declining trend with the latest estimate at 0.7 (95% CI 0.3, 1.0) by February 14, 2020

Based on the entire distribution of cluster sizes, we jointly estimated the overall reproduction number R and the dispersion parameter k as of March 17, 2020. Fitting the negative binomial distribution to the cluster data in the empirical distributions of the realizations during the early stages of the outbreak in Singapore, the reproduction number is estimated at 0.61 (95% CI 0.39, 1.02) after adjusting for the truncation of the time series leading to the possibility that some infected cases might still cause new infections after March 17, 2020. The dispersion parameter is estimated at 0.11 (95% CI 0.05, 0.25) consistent with SARS-CoV-2 transmission heterogeneity.

Discussion

Overall, current estimates of transmission potential in Singapore, based on two different data sources and different methods, suggest that temporary local transmission potential of SARS-CoV-2 has occurred in Singapore while our most recent estimate of the effective reproduction number is below the epidemic threshold of 1.0 whereas the overall reproduction number derived from the distribution of cluster sizes just barely crosses 1.0 (R = 0.61 (95% CI 0.39, 1.02). Temporary sustained transmission in the beginning of the epidemic can be partly attributed to multiple case importations and initiation of local transmission in the region. While large-scale local transmission has not been reported in Singapore, the fact that asymptomatic and subclinical cases are now well documented for COVID-19 [54] suggests that our estimates could be underestimated [55]. On the other hand, it is not clear if asymptomatic or subclinical cases are as infectious as symptomatic cases. Indeed, we have reported that multiple local cases have yet to be traced to existing transmission chains [9]. Additional data collected during the course of the outbreak will help obtain an improved picture of the transmission dynamics. These findings emphasize the need to strengthen public health interventions including active case contact tracing activities in countries with emerging transmission of SARS-CoV-2 [56]. It is worth noting that imported cases have minor contribution to secondary cases in Singapore, with most of the imported cases dating back to the early phase of the epidemic and between March 10 and 17, 2020. However, there are examples such as the Grant Hyatt Singapore cluster and the Yong Thai Hang cluster that were linked to imported sources, and the original sources had left Singapore before these local clusters emerged [30].

Our Rt estimates for Singapore are substantially lower than mean estimates reported for the COVID-19 epidemic in other parts of the world [57,58,59,60,61,62,63,64,65]. This indicates that containment efforts have a significant impact in Singapore (Table 2). However, some differences in the reproduction numbers reported for the epidemic in China may result from different methods, differences in data sources and time periods used to estimate the reproduction number. Similarly, a recent study has shown an average reporting delay of 6.1 days in China [66] which agrees with our mean estimate for cases in Singapore (6.4 days). Moreover, the scaling parameter for growth rate (p) indicates a sub-exponential growth pattern in Singapore, reflecting the effective isolation and control strategies in the region. This is consistent with a sub-exponential growth pattern for Chinese provinces excluding Hubei (p ~ 0.67), as estimated by a recent study [67] whereas an exponential growth pattern was estimated for Hubei (p ~ 1.0) [67].

Table 2 Timeline of COVID-19 epidemic in Singapore as of March 31, 2020 [9]

A previous study on the 2015 MERS outbreak in South Korea reported substantial potential for superspreading transmission despite a subcritical reproduction number [68]. The lower estimate of the dispersion parameter in our study also indicates significant transmission heterogeneity in Singapore. Superspreading events of MERS-CoV and SARS-CoV associated with nosocomial outbreaks are well documented and driven largely by substantial diagnostic delays [16]. Although the average delay from onset of symptoms to diagnosis/reporting for COVID-19 patients in Singapore is at 6.4 days and no superspreading events has been observed yet, the dispersion parameter, k < 1, indicates the probability of observing large clusters and the potential for superspreading such as the SAFRA Jurong cluster [50, 68]. Therefore, public health measures enacted by public health authorities in Singapore that advise the public to avoid mass gatherings and confined places are crucial to prevent disease amplification events. However, the presence of asymptomatic cases in the community represents an ongoing threat although it is not currently known if subclinical cases are less infectious [20, 69]. This highlights the need for rapid testing of suspected cases to quickly isolate those that test positive for the novel coronavirus. To achieve this goal, public health authorities in Singapore are reactivating 900 general practitioner clinics [70]. While new clusters emerge in Singapore, some clusters including the Yong Thai Hang cluster, Seletar Aerospace cluster, Wizlearn Technologies, and the Grand Hyatt cluster have stabilized (no recent additional cases in most clusters). The “Grace Assembly of God and the Life Church and Missions” cluster and the “SAFRA Jurong” cluster continue to be consolidated .

Beyond Singapore, COVID-19 cases are now being reported in 204 countries including identifiable clusters in many parts of the world [2, 71,72,73,74]. Moreover, Singapore has also produced secondary chains of disease transmission beyond its borders [24]. Although Singapore has been detecting and isolating cases with diligence, our findings underscore the need for continued and sustained containment efforts to prevent large-scale community transmission including nosocomial outbreaks. Overall, the current situation in Singapore highlights the need to investigate the imported, unlinked, and asymptomatic cases that could be a potential source of secondary cases and amplified transmission in confined settings. Although Singapore has a world-class health system including a highly efficient contact tracing mechanism in place that has prevented the outbreak from getting out of control [24, 75], continued epidemiological investigations and active case finding efforts are needed to contain the outbreak.

Our study is not exempt from limitations. First, the outbreak is still ongoing and we continue to monitor the transmission potential of COVID-19 in Singapore. Second, onset dates are missing for 12 cases, which were excluded from our analyses. Third, we cannot rule out that additional cases will be added to existing clusters, which may lead to underestimating the reproduction number based on the cluster size distribution. Fourth, some of the cases are associated with generating secondary chains in more than one cluster, which were included in the most relevant cluster.

Conclusion

This is a real-time study to estimate the evolving transmission potential of SARS-CoV-2 in Singapore. Our current findings point to temporary sustained transmission of SARS-CoV-2, with our most recent estimate of the effective reproduction number lying below 1.02. These estimates highlight the significant impact of containment efforts in Singapore while at the same time suggest the need to maintain social distancing and active case finding efforts to stomp out all active or incoming chains of transmission.

Availability of data and materials

The data supporting the conclusion of this article is available from the Ministry of Health, Singapore website [27].

Abbreviations

COVID-19:

Coronavirus Disease 2019

MERS:

Middle Eastern Respiratory Syndrome

SARS:

Severe acute respiratory syndrome

SARS-CoV-2:

Severe acute respiratory syndrome-related coronavirus 2

References

  1. 1.

    World Health Organization. Coronavirus disease 2019 (COVID-19) Situation Report – 27. World Health Organization; 2020. Available from: https://www.who.int/docs/default-source/coronaviruse/situation-reports/20200216-sitrep-27-covid-19.pdf?sfvrsn=78c0eb78_4.

  2. 2.

    Worldometer. COVID-19 coronavirus / cases 2020. Available from: https://www.worldometers.info/coronavirus/coronavirus-cases/.

  3. 3.

    Nishiura H, Jung S, Linton N, Kinoshita R, Yang Y, Hayashi K, et al. The extent of transmission of novel Coronavirus in Wuhan, China, 2020. J Clin Med. 2020;9(2):330.

    Article  Google Scholar 

  4. 4.

    Tian H, Li Y, Liu Y, Kraemer M, Chen B, Cai J, et al. Early evaluation of Wuhan City travel restrictions in response to the 2019 novel coronavirus outbreak. medRxiv. 2020:2020.01.30.20019844.

  5. 5.

    World Health Organization. Novel Coronavirus (2019-nCoV), Situation report - 1. 2020. Available from: https://www.who.int/docs/default-source/coronaviruse/situation-reports/20200121-sitrep-1-2019-ncov.pdf?sfvrsn=20a99c10_4. Accessed 25 Feb 2020.

  6. 6.

    World Health Organization. Situation Reports Coronavirus World Health Organization 2020. Available from: https://www.who.int/emergencies/diseases/novel-coronavirus-2019/situation-reports. Accessed 20 Feb 2020.

  7. 7.

    Olowski LJ. Coronavirus worse than reported – here's how China is catching up to reality. Fox News; 2020.

  8. 8.

    Kottasová I. China's coronavirus numbers highlight the challenges of an evolving epidemic. CNN; 2020.

  9. 9.

    Ministry Of Health. News Highlights Ministry of Health Singapore. Available from: https://www.moh.gov.sg/news-highlights/. Accessed 31 Mar 2020.

  10. 10.

    Ang D. The Big Story: 92 Singaporeans evacuated from Wuhan; Hikikomori in Singapore. The Straits Times; 2020.

  11. 11.

    Tan A. Coronavirus: Flight with 174 Singaporeans evacuated from Wuhan lands at Changi airport on Sunday morning. The Straits Times; 2020.

  12. 12.

    Ministry Of Health. Confirmed cases of local transmission of novel coronavirus infection in Singapore Minstry of Health Singapore. Available from: https://www.moh.gov.sg/news-highlights/details/confirmed-cases-of-local-transmission-of-novel-coronavirus-infection-in-singapore. Accessed 20 Feb 2020.

  13. 13.

    Ministry Of Health. Five more cases discharged, three new cases of covid-19 infection confirmed Ministry of Health. Available from: https://www.moh.gov.sg/news-highlights/details/five-more-cases-discharged-three-new-cases-of-covid-19-infection-confirmed. Accessed 20 Feb 2020.

  14. 14.

    Ministry Of Health. Two more cases discharged; five new cases of COVID-19 infection confirmed. Ministry of Health; 2020. Available from: https://www.moh.gov.sg/news-highlights/details/two-more-cases-discharged-five-new-cases-of-covid-19-infection-confirmed. Accessed 6 Mar 2020.

  15. 15.

    Secon H. Nearly 3,400 Chinese healthcare workers have gotten the coronavirus, and 13 have died. Business Insider; 2020.

  16. 16.

    Chowell G, Abdirizak F, Lee S, Lee J, Jung E, Nishiura H, et al. Transmission characteristics of MERS and SARS in the healthcare setting: a comparative study. BMC Med. 2015;13(1):210.

    Article  Google Scholar 

  17. 17.

    Gopalakrishna G, Choo P, Leo Y, Tay B, Lim Y, Khan A, et al. SARS transmission and hospital containment. Emerg Infect Dis. 2004;10(3):395–400.

    Article  Google Scholar 

  18. 18.

    Lim J. Coronavirus: Doctors, staff and patients to restrict movements to within one hospital. The Straits Times Singapore; 2020.

  19. 19.

    Aravindan A, Geddie J. Singapore lifts virus alert to SARS level, sparking panic buying. Rueters; 2020.

  20. 20.

    Mizumoto K, Kagaya K, Zarebski A, Chowell G. Estimating the asymptomatic proportion of coronavirus disease 2019 (COVID-19) cases on board the Diamond Princess cruise ship, Yokohama, Japan, 2020. Euro Surveill. 2020;25(10):2000180.

    Article  Google Scholar 

  21. 21.

    Nishiura H, Kobayashi T, Miyama T, Suzuki A, Jung SM, Hayashi K, et al. Estimation of the asymptomatic ratio of novel coronavirus infections (COVID-19). Int J Infect Dis. 2020;94:154-5.

  22. 22.

    Hu Z, Song C, Xu C, Jin G, Chen Y, Xu X, et al. Clinical characteristics of 24 asymptomatic infections with COVID-19 screened among close contacts in Nanjing, China. Sci China Life Sci. 2020;63(5):706-11.

  23. 23.

    Sim D. Coronavirus: why did Singapore have more cases than Hong Kong – until now? This Week in Asia; 2020.

  24. 24.

    Vaswani K. Coronavirus: Why Singapore is so vulnerable to coronavirus spread. BBC News; 2020.

  25. 25.

    Diekmann O, Heesterbeek JAP. Mathematical Epidemiology of Infectious Diseases: Model Building, Analysis and Interpretation. 1st ed. Wiley; 2000. ISBN 0-471-49241-8.

  26. 26.

    Anderson RM, May RM. Infectious Diseases of Humans. Oxford: Oxford Univeristy Press; 1991. ISBN 978-0198540403.

  27. 27.

    Ministry Of Health. Data retrieved from updates on COVID-19, Press Release. Ministry of Health Singapore; 2020. Available  from: https://www.moh.gov.sg/covid-19. Accessed 29 Mar 2020.

  28. 28.

    Khalik S. Coronavirus in Singapore: of the 50 cases so far, local cases exceed imported ones. The Straits Times; 2020.

  29. 29.

    CNA. Coronavirus cases in Singapore: trends, clusters and key numbers to watch. CNA; 2020.

  30. 30.

    Wong L. Coronavirus: Full recovery for all three local cases linked to Grand Hyatt event. The Straits Times; 2020.

  31. 31.

    MalayMail. COVID-19: SAF regular among three new victims in Singapore as Grace Assembly of God cluster grows to 18 cases. Malaymail; 2020.

  32. 32.

    Yan P, Chowell G. Quantitative methods for infectious disease outbreak investigations. Springer; 2019. ISBN 978-3-030-21923-9.

  33. 33.

    Yan P. Estimation for the infection curves for the spread of Severe Acute Respiratory Syndrome (SARS) from a back-calculation approach Center for Discrete Mathematics & Theoratical Computer Science Founded as a National Science Foundation Science and Technology Center 2018. Available from: http://archive.dimacs.rutgers.edu/Workshops/Modeling/slides/Yan.pdf. Accessed 15 Feb 2020.

  34. 34.

    Lawless JF. Adjustments for reporting delays and the prediction of occurred but not reported events. Can J Stat. 1994;22(1):15-31.

  35. 35.

    Taylor GC. Claims Reserving In Non Life Insurance: Elsevier; 1985. ISBN 978-0444878465.

  36. 36.

    Tariq A, Roosa K, Mizumoto K, Chowell G. Assessing reporting delays and the effective reproduction number: the Ebola epidemic in DRC, May 2018-January 2019. Epidemics. 2019;26:128–33.

    CAS  Article  Google Scholar 

  37. 37.

    Kalbfleisch JD, Lawless JF. Estimating the incubation time distribution and expected number of cases of transfusion-associated acquired immune deficiency syndrome. Transfusion. 1989;29(8):672–6.

    CAS  Article  Google Scholar 

  38. 38.

    Nishiura H, Chowell G, Heesterbeek H, Wallinga J. The ideal reporting interval for an epidemic to objectively interpret the epidemiological time course. J R Soc Interface. 2010;7(43):297–307.

    Article  Google Scholar 

  39. 39.

    Chowell G, Hincapie-Palacio D, Ospina J, Pell B, Tariq A, Dahal S, et al. Using phenomenological models to characterize transmissibility and forecast patterns and final burden of Zika epidemics. PLoS Curr. 2016;8:ecurrents.outbreaks.f14b2217c902f453d9320a43a35b583.

    Google Scholar 

  40. 40.

    Shanafelt DW, Jones G, Lima M, Perrings C, Chowell G. Forecasting the 2001 foot-and-mouth disease epidemic in the UK. Ecohealth. 2018;15(2):338–47.

    Article  Google Scholar 

  41. 41.

    Pell B, Kuang Y, Viboud C, Chowell G. Using phenomenological models for forecasting the 2015 Ebola challenge. Epidemics. 2018;22:62–70.

    Article  Google Scholar 

  42. 42.

    Chowell G, Luo R, Sun K, Roosa K, Tariq A, Viboud C. Real-time forecasting of epidemic trajectories using computational dynamic ensembles. Epidemics. 2020;30:100379.

    Article  Google Scholar 

  43. 43.

    Nishiura H, Linton NM, Akhmetzhanov AR. Serial interval of novel coronavirus (COVID-19) infections. Int J Infect Dis. 2020;93:284-6.

  44. 44.

    You C, Deng Y, Hu W, Sun J, Lin Q, Zhou F, et al. Estimation of the time-varying reproduction number of COVID-19 outbreak in China. Int J Hyg Environ Health. 2020:113555.

  45. 45.

    Nishiura H, Chowell G. The Effective Reproduction Number as a Prelude to Statistical Estimation of Time-Dependent Epidemic Trends. Springer D, editor 2009. 103-12 p. Mathematical and Statistical Estimation Approaches in Epidemiology. ISBN 978-90-481-2313-1.

  46. 46.

    Paine S, Mercer G, Kelly P, Bandaranayake D, Baker M, Huang Q, et al. Transmissibility of 2009 pandemic influenza A(H1N1) in New Zealand: effective reproduction number and influence of age, ethnicity and importations. Euro surveill. 2010;15(24). https://doi.org/10.2807/ese.15.24.19591-en.

  47. 47.

    Nishiura H, Roberts MG. Estimation of the reproduction number for 2009 pandemic influenza a(H1N1) in the presence of imported cases. Eurosurveillance. 2010;15(29):19622.

    Article  Google Scholar 

  48. 48.

    Chowell G. Fitting dynamic models to epidemic outbreaks with quantified uncertainty: a primer for parameter uncertainty, identifiability, and forecasts. Infect Dis Model. 2017;2(3):379–98.

    PubMed  PubMed Central  Google Scholar 

  49. 49.

    Blumberg S, Lloyd-Smith JO. Inference of R0 and transmission heterogeneity from the size distribution of stuttering chains. PLoS Comput Biol. 2013;9(5):e1002993.

    CAS  Article  Google Scholar 

  50. 50.

    Lloyd-Smith JO, Schreiber SJ, Kopp PE, Getz WM. Superspreading and the effect of individual variation on disease emergence. Nature. 2005;438(7066):355–9.

    CAS  Article  Google Scholar 

  51. 51.

    Blumberg S, Lloyd-Smith JO. Comparing methods for estimating R0 from the size distribution of subcritical transmission chains. Epidemics. 2013;5(3):131–45.

    CAS  Article  Google Scholar 

  52. 52.

    Blumberg S, Funk S, Pulliam JRC. Detecting differential transmissibilities that affect the size of self-limited outbreaks. PLoS Pathog. 2014;10(10):e1004452.

    Article  Google Scholar 

  53. 53.

    Bolker BM. Ecological models and data in R. Press PU, editor 2008.

  54. 54.

    The Novel Coronavirus Pneumonia Emergency Response Epidemiology Team. The Epidemiological Characteristics of an Outbreak of 2019 Novel Coronavirus Diseases (COVID-19) — China, 2020[J]. China CDC Weekly. 2020;2(8):113-22. https://doi.org/10.46234/ccdcw2020.032.

  55. 55.

    Belluz J. Why the coronavirus outbreak might be much bigger than we know. Vox; 2020.

  56. 56.

    Pung R, Chiew CJ, Young BE, Chin S, Chen MIC, Clapham HE, et al. Investigation of three clusters of COVID-19 in Singapore: implications for surveillance and response measures. Lancet. 2020;395(10229):1039–46.

    CAS  Article  Google Scholar 

  57. 57.

    Shen M, Peng Z, Xiao Y, Zhang L. Modelling the epidemic trend of the 2019 novel coronavirus outbreak in China. bioRxiv.2020.01.23.916726.

  58. 58.

    Mizumoto K, Kagaya K, Chowell G. Early epidemiological assessment of the transmission potential and virulence of 2019 Novel Coronavirus in Wuhan City: China, 2019–2020. medRxiv. 2020:2020.02.12.20022434.

  59. 59.

    Zhao S, Lin Q, Ran J, Musa SS, Yang G, Wang W, et al. Preliminary estimation of the basic reproduction number of novel coronavirus (2019-nCoV) in China, from 2019 to 2020: a data-driven analysis in the early phase of the outbreak. Int J Infect Dis. 2020;92:214–7.

    CAS  Article  Google Scholar 

  60. 60.

    Liu T, Hu J, Kang M, Lin L, Zhong H, Xiao J, et al. Transmission dynamics of 2019 novel coronavirus (2019-nCoV). bioRxiv.2020.01.25.919787.

  61. 61.

    Wu JT, Leung K, Leung GM. Nowcasting and forecasting the potential domestic and international spread of the 2019-nCoV outbreak originating in Wuhan, China: a modelling study. The Lancet. 2020;395(10225):689-97.

  62. 62.

    Read JM, Bridgen JR, Cummings DA, Ho A, Jewell CP. Novel coronavirus 2019-nCoV: early estimation of epidemiological parameters and epidemic predictions. medRxiv. 2020:2020.01.23.20018549.

  63. 63.

    Li Q, Guan X, Wu P, Wang X, Zhou L, Tong Y, et al. Early Transmission Dynamics in Wuhan, China, of Novel Coronavirus–Infected Pneumonia. N Engl J Med. 2020:382:1199-207. https://doi.org/10.1056/NEJMoa2001316.

  64. 64.

    Sanche S, Lin YT, Xu C, Romero-Severson E, Hengartner N, Ke R. High contagiousness and rapid spread of severe acute respiratory syndrome coronavirus 2. Emerg Infect Dis. 2020. https://doi.org/10.3201/eid2607.200282.

  65. 65.

    Shim E, Tariq A, Choi W, Lee Y, Chowell G. Transmission potential and severity of COVID-19 in South Korea. Int J Infect Dis. 2020;93:339-44.

  66. 66.

    Kucharski AJ, Russell TW, Diamond C, Liu Y, Edmunds J, Funk S, et al. Early dynamics of transmission and control of COVID-19: a mathematical modelling study. Lancet Infect Dis. 2020;20(5):553-8.

  67. 67.

    Roosa K, Lee Y, Luo R, Kirpich A, Rothenberg R, Hyman JM, et al. Real-time forecasts of the COVID-19 epidemic in China from February 5th to February 24th, 2020. Infect Dis Model. 2020;5:256–63.

    CAS  PubMed  PubMed Central  Google Scholar 

  68. 68.

    Kucharski AJ, Althaus CL. The role of superspreading in Middle East respiratory syndrome coronavirus (MERS-CoV) transmission. Eurosurveillance. 2015;20(25):21167.

    Article  Google Scholar 

  69. 69.

    Sim D, Xinghui K. Singapore: nine new coronavirus cases, no plans to raise outbreak alert to red. This Week in Asia; 2020.

  70. 70.

    Chang N, Yong M. Public health preparedness clinics reactivated to reduce risk of COVID-19 spread. CNA; 2020.

  71. 71.

    Marcus I. Chronology: Germany and Covid 19 (‘Coronavirus’). The Berlin Spectator; 2020.

  72. 72.

    Rueters. Germany Confirms Seventh Coronavirus Case. Rueters; 2020.

  73. 73.

    McCurry J. Japan reports first coronavirus death as 44 more cases confirmed on cruise ship. The Guardian; 2020.

  74. 74.

    Woods A. 44 new coronavirus cases reported on quarantined Diamond princess cruise. Newyork Post; 2020.

  75. 75.

    Ng Y, Li Z, Chua YX, Chaw WL, Zhao Z, Er B, et al. Evaluation of the Effectiveness of Surveillance and Containment Measures for the First 100 Patients with COVID-19 in Singapore — January 2–February 29, 2020. MMWR Early Release; March 20, 2020. Centers for Disease Control and Prevention.

Download references

Acknowledgements

Not applicable.

Funding

G.C. is supported by NSF grants 1610429 and 1633381. G.C. and S.B. are partially supported by R01 GM 130900.

Author information

Affiliations

Authors

Contributions

A.T, S.B., P.Y., and G.C. analyzed the data. A.T., Y.L, P.Y, and S.M. retrieved and managed data; A.T and G.C wrote the first draft of the manuscript. All authors contributed to writing and revising subsequent versions of the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Amna Tariq.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

G.C. is a Board Member for the journal. The other authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Tariq, A., Lee, Y., Roosa, K. et al. Real-time monitoring the transmission potential of COVID-19 in Singapore, March 2020. BMC Med 18, 166 (2020). https://doi.org/10.1186/s12916-020-01615-9

Download citation

Keywords

  • SARS-CoV-2
  • COVID-19
  • Singapore
  • Transmission potential
  • Transmission heterogeneity
  • Reproduction number
  • Cluster distribution
  • Reporting delay