key: cord-0763973-0fmgaa0h
authors: Wilasang, Chaiwat; Jitsuk, Natcha C.; Sararat, Chayanin; Modchang, Charin
title: Reconstruction of the transmission dynamics of the first COVID-19 epidemic wave in Thailand
date: 2022-02-07
journal: Sci Rep
DOI: 10.1038/s41598-022-06008-x
sha: cfbd8a6d331fc5e562fcece6a604aa37ead895ed
doc_id: 763973
cord_uid: 0fmgaa0h

Thailand was the first country reporting the first Coronavirus disease 2019 (COVID-19) infected individual outside mainland China. Here we delineated the course of the COVID-19 outbreak together with the timeline of the control measures and public health policies employed by the Thai government during the first wave of the COVID-19 outbreak in Thailand. Based on the comprehensive epidemiological data, we reconstructed the dynamics of COVID-19 transmission in Thailand using a stochastic modeling approach. Our stochastic model incorporated the effects of individual heterogeneity in infectiousness on disease transmission, which allows us to capture relevant features of superspreading events. We found that our model could accurately capture the transmission dynamics of the first COVID-19 epidemic wave in Thailand. The model predicted that at the end of the first wave, the number of cumulative confirmed cases was 3091 (95%CI: 2782–3400). We also estimated the time-varying reproduction number (R(t)) during the first epidemic wave. We found that after implementing the nationwide interventions, the R(t) in Thailand decreased from the peak value of 5.67 to a value below one in less than one month, indicating that the control measures employed by the Thai government during the first COVID-19 epidemic wave were effective. Finally, the effects of transmission heterogeneity and control measures on the likelihood of outbreak extinction were also investigated.

The first wave of COVID-19 in Thailand. On 13 January 2020, Thailand reported the first confirmed case of COVID-19 outside China 4, 5, 22 . The first detected case was a traveler from Wuhan, China, who had symptom onset on 3 January 2020. Not long after that, the first locally transmitted COVID-19 case in Thailand was reported on 22 January 2020. The Thai government then enforced several control measures to prevent and mitigate outbreaks 7-10 . Figure 1 shows the numbers of daily and cumulative confirmed locally transmitted COVID-19 cases along with some key events in Thailand during the first epidemic wave. At the end of the first epidemic wave (13 May 2020), 3,017 cumulative cases and 56 deaths (1.9%) had been recorded. Among 3017 infected individuals, 2492 (82.6%) are adults (aged 20-59 years), and 332 (11.0%) are elderly (aged over 60 years old) 8, 23 .

In Thailand, the COVID-19 transmission during the first epidemic wave was mostly related to large cluster outbreaks in crowded areas 6, 7, 9 . As shown in Fig. 1 , the number of daily new cases increased rapidly after two large cluster outbreaks at the Lumphini boxing stadium (6 March 2020) and at a nightclub in Thong Lo (9 March 2020) 6 . The Thai government then employed several control measures to mitigate the spread of COVID-19. Physical distancing and work-from-home measures were implemented. Nightclubs, restaurants, stores, and schools were ordered to close across the country. Also, the Songkran Festival, i.e., the Thai New Year festival, was canceled to avoid mass gatherings 24, 25 . In addition, the government implemented a nationwide curfew (between 10 p.m. to 4 a.m) on 2 April 2020, and all international flights to Thailand were banned starting from 6 April 2020 10, 23 . The daily confirmed cases were subsequently decreased, leading to the lockdown easing on 1 May 2020 8, 10 . The number of locally transmitted COVID-19 cases during the first epidemic wave in Thailand had reached zero for the first time on 13 May 2020 and had since remained zero for at least 200 consecutive days 8, 10, 11 .

To evaluate the effectiveness of the implemented control measures, we then estimated the time-varying reproduction number, R t , which represents the average number of secondary cases caused by a primary infected individual. A value of R t greater than the threshold value of 1 indicates that the epidemic size is growing at time t, whereas R t < 1 indicates that the epidemic size is shrinking at time t. The goal of the control measures is, therefore, to reduce R t < 1. As shown in Fig. 2 , the R t in the early phase of transmission was around 2.72 (95% CI 1. 56-4.35) . However, after the large outbreaks at the Lumphini boxing stadium and the nightclub in Thong Lo, the R t increased rapidly and peaked at 5.67 (95% CI 4.62-6.87) on 18 March 2020, which was just before nightclubs, restaurants, and stores were ordered to close. After Thailand had implemented several control measures, the R t decreased to a value < 1 for the first time on 6 April 2020 and then had since remained below one until the end of the first epidemic wave.

To reconstruct the transmission dynamics of COVID-19 in Thailand, we used the estimated R t to calculate the time-dependent transmission rate of SARS-CoV-2 in Thailand (see Material and Methods). The comparison between the modeling results and the reported data is shown in Fig. 3 . We found that the trends of both daily new cases and cumulative cases generated from the model agree well with the corresponding observed data (R-square = 0.9996). The mean time delay from onset of infectiousness to isolation, representing the effective infectious period of isolated individuals, was estimated at 2.26 days. The model predicted that at the end of the first wave, the number of cumulative confirmed cases was 3,091 (95%CI 2782-3400), which is close to the number of reported cases. mission dynamics, we employed the constructed model to simulate the SARS-CoV-2 transmission and estimated the probability of outbreak extinction after the introduction of a single infected individual. We found that 64.1% of the model simulations exhibit stochastic extinction (defined as a realization that has no latently infected, asymptomatic, and symptomatic infectious individuals in 30 days after the introduction), with an average time to extinction of 9.6 days. Moreover, of 77.84% of the extinct outbreaks, the index case did not generate any secondary infection. In this work, following Lloyd-Smith et al. 26 , a superspreading event (SSE) is defined as any outbreak event that involves infected individuals who, on average, infect others more than the 99th percentile of the Poisson (R 0 ) distribution. Based on this criterion, we found that 11.1% of the model realizations involve SSEs. Also, as can be seen in Fig. 4 , an outbreak that involves SSEs usually spreads faster than the others.

Evaluating the impact of interventions and heterogeneity in infectiousness on SARS-CoV-2 transmission dynamics. Next, we employed the model to evaluate the impact of interventions and individual variation in infectiousness on the SARS-CoV-2 transmission dynamics and outbreak extinction. Specifically, we explored the intervention scenarios where the mean time delay from onset of infectiousness to isolation (T Q ) is 0.5, 1.0, 2.0, 3.0, and 5.0 days together with the symptomatic case isolation ( q S ) ranging from 10 to 100%. Our modeling results indicated that the outbreak would be more likely to become extinct if more symptomatic cases are isolated or case isolation is performed sooner (Fig. 5A ). For example, more than 80% of the outbreaks will become extinct if all symptomatic cases are isolated within the first day of their infectiousness. Moreover, the isolation of symptomatic cases can also affect the likelihood of SSEs. Fast and efficient case isolation can almost eliminate the likelihood of SSEs (Fig. 5B ). We also explored the intervention scenarios where the asymptomatic case isolation (q A ) is different from the default value of 60%. In this case, we also found similar results in which the outbreak is more likely to become extinct when more asymptomatic cases are isolated ( Fig. S1 in the supplementary information). Figure 6A depicts the role of individual variation in infectiousness on the outbreak extinction. Here, the dispersion parameter k represents the individual heterogeneity in disease transmission. A very high value of k indicates that each primary case generates roundly the same number of secondary cases. In contrast, a low value of k indicates that only a small fraction of infected individuals disproportionately infects a large number of individuals 26 . We found that individual variation on disease transmission favors the outbreak extinction, which is consistent with the previous studies 20, [26] [27] [28] . Moreover, R 0 can also affect the likelihood of outbreak extinction (Fig. 6B ). An outbreak with a lower value of R 0 is more likely to become extinct. However, lowering R 0 could also reduce the impact of individual heterogeneity on SARS-CoV-2 transmission. Finally, to investigate the impact of transmission heterogeneity on the speed of disease transmission, we calculated the first date at which the cumulative number of infected individuals exceeds 100 (T 100 ) (Fig. 6C ). We found that an outbreak with a higher level of transmission heterogeneity spreads slightly faster. Specifically, T 100 decreases from about 28 days for the homogeneous transmission ( k → ∞ ) to approximately 26 days for the transmission with a very high degree of individual heterogeneity (k = 0.01). Also, we found no SSE-related outbreak when k ≥ 10 . In addition, for non-SSE outbreaks with k ≤ 0.05 , the cumulative number of infected individuals could not reach 100 within the simulation time of 30 days. 

In this work, we investigated the first wave of the COVID-19 outbreak in Thailand together with the timeline of the control measures employed by the Thai government. We also estimated the time-varying reproduction number during the first epidemic wave and used stochastic simulations to reconstruct the full transmission dynamics of COVID-19 in Thailand. Recent studies have revealed important transmission features of COVID-19, including the infectiousness of asymptomatic 29, 30 , individual heterogeneity in infectiousness, and SSEs 14, 16, 19, 31 . To accurately capture the COVID-19 transmission dynamics in Thailand, these features were therefore incorporated in our stochastic event-based model. Specifically, our model assumes that both symptomatic and asymptomatic infectious individuals can transmit the disease but with different degrees of infectiousness. In addition, the expected number of secondary cases ( ν ) caused by a primary infectious individual is assumed to follow a gamma distribution with a mean reproduction number R and a dispersion parameter k 26 .

After declaring the state of emergency and implementing the nationwide interventions, the significant decline in both the number of confirmed cases and the positive detection rate in Thailand was observed 32, 33 . To evaluate the effectiveness of the control measures implemented by the Thai government during the first epidemic wave, the time-varying reproduction number, R t , was estimated (Fig. 2) . We found that after implementing the nationwide interventions, the R t in Thailand was decreased from the peak value of 5.67 to a value below one in less than one month. This indicated that the control measures employed by the Thai government during the first COVID-19 epidemic wave were effective. Singapore is another country in Southeast Asia that was successful in controlling the first wave of COVID-19 34 . It was found that the mean time delay from onset of infectiousness to isolation in Singapore was about 2.6 days 35 , which is close to our estimate of 2.26 days in Thailand.

In general, mathematical modeling and computer simulation can be used either as predictive tools or as a means of understanding infectious transmission dynamics 36 . Here, the model was employed to understand the To this end, we incorporated the effects of the implemented interventions during the first COVID-19 epidemic wave by using the estimated R t to estimate the time-varying individual reproduction number ( ν t ) in our model 26 . We found that our model could www.nature.com/scientificreports/ accurately capture the transmission dynamics of the first COVID-19 epidemic wave in Thailand (Fig. 3) . The model predicted that at the end of the first wave, the number of cumulative confirmed cases was 3,091 (95%CI 2782-3400), which is consistent with the reported data. In fact, several modeling studies have been conducted to predict the incidence of COVID-19 and evaluate the effectiveness of control measures 24, 33, [37] [38] [39] . However, most previous models usually assumed a Poisson distribution of the number of secondary infections per infected individual, which could not capture the relevant features of superspreading events (SSEs) 12, 26 . An SSE is usually defined as an event in which an infected individual transmits the disease to a large number of secondary individuals, while most infected individuals infect a few or none 20, 21 . Recent studies have pointed out that only a small fraction of SARS-CoV-2 infected individuals disproportionately infect a large group of people while most infected individuals infect a few or none 20, 21 . It was estimated that only 10-20% of the COVID-19 infected individuals were responsible for about 80% of transmission 12, 16, 26, [40] [41] [42] . To accurately reconstruct the full dynamics of COVID-19 transmission, these features of disease transmission heterogeneity were therefore incorporated in our model. Specially, we assumed that the expected number of secondary cases caused by an infectious individual varies from person to person 26 .

We found that the dynamics of COVID-19 transmission depend on the degree of the individual heterogeneity in disease transmission, represented by the dispersion parameter k. Our modeling results indicated that individual variation on disease transmission favors the outbreak extinction, which is in line with previous studies 20, [26] [27] [28] (Fig. 6) . This is because when individual heterogeneity in disease transmission is high, only a few infected individuals generate a large number of secondary infections, while most infected individuals give rise to a few or none secondary infections, leading to the higher likelihood of the outbreak extinction. However, although individual variation on disease transmission could increase the chance of outbreak extinction, it could also accelerate the transmission dynamics of the non-extinct outbreaks, which is in agreement with other studies 12, 16, 26 .

Our study, as with all modeling studies, has several limitations. The estimation of R t was based on the reported number of confirmed cases. Although there might be potential biases due to delays in case reporting, if the delay from infection to confirmation does not change over time, this will not affect the conclusion of this study. In our study, we also assumed that the testing and reporting efforts are constant over time. However, if the testing effort increases and decreases during a particular time interval, this will increase and decrease the estimated values of R t , respectively 3 . In addition, although recent studies revealed that age-specific contact patterns also play a role in the COVID-19 transmission 18,43,44 , we did not incorporate them into our model; such data for Thailand is not yet available. Finally, our model did not consider migrant movements and human mobility; incorporating these data into the epidemic model might improve the modeling accuracy 45, 46 .

Data sources. The number of laboratory-confirmed COVID-19 cases in Thailand was retrieved from the Center for Systems Science and Engineering (CSSE) at Johns Hopkins University 47 . The information of implemented control measures and large outbreak events was collected from the Centre for COVID-19 Situation Administration (CCSA) of Thailand 8 .

Duration of the first epidemic wave. The first wave of COVID-19 outbreak in Thailand was assumed to start from 22 January 2020, when the first locally transmitted COVID-19 case was reported, to 13 May 2020, which was the first date at which Thailand had no reported on COVID-19 local transmission for at least 200 consecutive days 8,11 .

Estimation of the effective reproduction number. The effective reproduction number, R t , is defined as the average number of secondary cases produced by an infected individual at time t 3 . A value of R t greater than the critical value of 1 indicates that the epidemic size is growing and the infection could be able to spread in the population at time t , whereas a value of R t less than 1 indicates that the epidemic size is shrinking at time t 3 . In this study, we employed a statistical method developed by Cori et al. to estimate R t during the first epidemic wave in Thailand 48, 49 . This R t estimation method only requires the number of daily new cases and the distribution of the corresponding serial interval, which here was assumed to be a discretized Gamma distribution with a mean and the standard deviation of 6.5 days and 4.2 days, respectively 3, 35, 50 . Details of the R t estimation are presented in ref 3 .

A stochastic event-based modeling approach was employed in this study. The schematic of the model specification is shown in Fig. 7 . The model divides the population into six separated compartments; namely, susceptible (S), latently infected (L), symptomatic infectious (I S ), asymptomatic infectious (I A ), quarantined (Q), and recovered (R) compartments. A susceptible individual can get an infection from either a symptomatic or an asymptomatic infectious individual at rates β S and β A , respectively. After being infected, the susceptible individual progresses to the latently infected class. Individuals in this class have already acquired the SARS-CoV-2 infection but are not yet infectious and cannot transmit the disease to other susceptible individuals. Latently infected individuals become infectious at a rate that is inversely proportional to the mean latent period (T L ). A proportion r S of the infected individuals becomes symptomatic infectious individuals, while the remaining 1 − r S becomes asymptomatic infectious individuals. In addition, a proportion q S of symptomatic infectious individuals and a proportion q A of asymptomatic infectious individuals is isolated and quarantined at a rate that is inversely propositional to the mean time delay from onset of infectiousness to isolation (T Q ) while the remaining is recovered at a rate that is inversely proportional the mean infectious period (T I ). We assumed that all isolated individuals could not make a further transmission. 26 . In our work, the individual reproduction number ( ν ) of each infectious individual was assumed to follow a gamma distribution with mean R 0 and dispersion parameter k , ν ∼ gamma(R 0 , k) , where R 0 is the basic reproduction number 12, 26 . The transmission rates due to symptomatic infectious individuals ( β S ) and asymptomatic infectious individuals ( β A ) are given by where ν Q = T Q /T I ν is the expected number of secondary cases generated by a quarantined individual before isolation. Here we assumed that ν Q is linearly proportional to the mean time delay from onset of infectiousness to isolation (T Q ). γ S and γ A are the symptomatic and asymptomatic recovery rates that are inversely proportional to the mean infectious period.

The total number of infection events that occurred during the time interval dt is given by where n I is the total number of infectious individuals at time t, and i = S or A representing symptomatic or asymptomatic infectious individual, respectively. The numbers of all other transition events during the time interval dt were calculated using the standard event-based modeling algorithm and τ-leap method 51 . The model simulations begin with a single infectious individual at time t = 0. For each set of parameters, we simulated 10 batches of 1,000 model realizations. Following Lloyd-Smith et al. 26 , a model realization that involves infected individuals who, on average, infect others more than the 99th percentile of the Poisson (R 0 ) distribution is classified as an SSE-related outbreak. An outbreak extinction was defined as a realization that has no latently infected, asymptomatic, and symptomatic infectious individuals in 30 days after the introduction of the index 

To reconstruct the transmission dynamics of COVID-19 in Thailand, we used the estimated R t during the first epidemic wave in Thailand to estimate the time-varying individual reproduction number ( ν t ). Specifically, ν t was drawn from the gamma distribution with mean R t and dispersion parameter k , ν t ∼ gamma(R t , k) . The model structure employed in this section is the same as the one described in the Model structure section, except that the individual reproduction number here is time-dependent. In addition, based on the surveillance information, the initial number of symptomatic infectious individuals was set at five, while the initial numbers of latently infected, quarantined, recovered, and asymptomatic infectious individuals were set at zero. 

A novel coronavirus from patients with pneumonia in China

Early dynamics of transmission and control of COVID-19: A mathematical modelling study

Reduction in effective reproduction number of COVID-19 is higher in countries employing active case detection with prompt isolation

Timeline: WHO's COVID-19 response

Coronavius disease 2019 (COVID-19) WHO Thailand Situation Report-22

COVID-19-WHO Thailand Situation Reports

Effects of public health interventions on the epidemiological spread during the first wave of the COVID-19 outbreak in Thailand

Coronavius disease 2019 (COVID-19) WHO Thailand Situation Report

Superspreading events in the transmission dynamics of SARS-CoV-2: Opportunities for interventions and control

Identifying and interrupting superspreading events-implications for control of severe acute respiratory syndrome coronavirus 2

Evaluating transmission heterogeneity and super-spreading event of COVID-19 in a metropolis of China

What settings have been linked to SARS-CoV-2 transmission clusters?

Estimating the overdispersion in COVID-19 transmission using outbreak sizes outside China

Pattern of early human-to-human transmission of Wuhan

Characterizing superspreading events and age-specific infectiousness of SARS-CoV-2 transmission in Georgia

Superspreading in early transmissions of COVID-19 in Indonesia

Super-spreaders in infectious diseases

COVID-19 and SARS-CoV-2. Modeling the present, looking at the future

Early transmission dynamics in Wuhan, China, of novel coronavirus-infected pneumonia

Corona Virus Disease

Effect of the Songkran festival on COVID-19 transmission in Thailand. Asian Pac

Superspreading and the effect of individual variation on disease emergence

An event-based model of superspreading in epidemics

The effect of superspreading on epidemic outbreak size distributions

Temporal dynamics in viral shedding and transmissibility of COVID-19

Presumed asymptomatic carrier transmission of COVID-19

Transmission heterogeneities, kinetics, and controllability of SARS-CoV-2

Positive impact of lockdown on COVID-19 outbreak in Thailand

Evaluating the impact of intervention strategies on the first wave and predicting the second wave of COVID-19 in Thailand: A mathematical modeling study

Singapore's pandemic preparedness: an overview of the first wave of COVID-19

Effectiveness of isolation, testing, contact tracing, and physical distancing on reducing transmission of SARS-CoV-2 in different settings: A mathematical modelling study

Dynamics of Zika virus outbreaks: an overview of mathematical modeling approaches

Nowcasting and forecasting the potential domestic and international spread of the 2019-nCoV outbreak originating in Wuhan, China: A modelling study

Modeling and forecasting the early evolution of the Covid-19 pandemic in Brazil

Modeling, state estimation, and optimal control for the US COVID-19 outbreak

Clustering and superspreading potential of SARS-CoV-2 infections in Hong Kong

Why do some COVID-19 patients infect many others, whereas most don't spread the virus at all

High SARS-CoV-2 attack rate following exposure at a choir practice

Evaluating how smartphone contact tracing technology can reduce the spread of infectious diseases: The case of COVID-19

Changes in contact patterns shape the dynamics of the COVID-19 outbreak in China

How the individual human mobility spatio-temporally shapes the disease transmission dynamics

Human mobility: Models and applications

A new framework and software to estimate time-varying reproduction numbers during epidemics

Practical considerations for measuring the effective reproductive number, R t

Epidemiology and Transmission of COVID-19 in Shenzhen China: Analysis of 391 cases and 1,286 of their close contacts

Modeling Infectious Diseases in Humans and Animals

The incubation period of coronavirus disease 2019 (COVID-19) from publicly reported confirmed cases: Estimation and application

Comparison of molecular testing strategies for COVID-19 control: A mathematical modelling study

Reconstruction of the full transmission dynamics of COVID-19 in Wuhan

Estimating the infection and case fatality ratio for coronavirus disease (COVID-19) using age-adjusted data from the outbreak on the Diamond Princess cruise ship

The reproductive number of COVID-19 is higher compared to SARS coronavirus

The authors declare that the data supporting the findings of this study are available within the paper. 

C.M. designed the research. C.W., N.C.J., C.S., and C.M. analyzed the results, edited and draft the manuscript. C.W. performed the simulations.

The authors declare no competing interests.

The online version contains supplementary material available at https:// doi. org/ 10. 1038/ s41598-022-06008-x.Correspondence and requests for materials should be addressed to C.M.

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.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:// creat iveco mmons. org/ licen ses/ by/4. 0/.