Abstract
Since the detection of first case of COVID-19 in Chile on March 3rd, 2020, a total of 301019 cases including 6434 deaths have been reported in Chile as of July 7th, 2020. In this manuscript we estimate the reproduction number during the early transmission phase in Chile and study the effectiveness of control interventions by conducting short-term forecasts based on the early transmission dynamics of COVID-19. We also estimate the reproduction number and conduct short term forecasts for the most recent developments in the epidemic trajectory of COVID-19 in Chile (May 9th-July 7th, 2020) to study the effectiveness of re-imposition of lockdowns in the country. The incidence curve in Chile displays early sub-exponential growth dynamics with the scaling of growth parameter, p, estimated at 0.8 (95% CI: 0.7, 0.8) and the reproduction number, estimated at 1.8 (95% CI: 1.6, 1.9). Our analysis emphasizes that the control measures at the start of the epidemic significantly slowed down the spread of the virus. However, easing of the COVID-19 restrictions and spread of virus to the low income neighborhoods in May led to a new wave of infections, followed by the re-imposition of lockdowns in Santiago and other municipalities. These measures have decelerated the virus spread with R estimated at ~0.87(95% CI: 0.84, 0.89) as of July 7th, 2020. Our current findings point that the sustained transmission of SARS-CoV-2 in Chile is being brought under control. The COVID-19 epidemic followed an early sub-exponential growth trend (p ~0.8) that transformed into a linear growth trend (p ~0.5) as of July 7th, 2020. While the broad scale social distancing interventions have slowed the virus spread, the number of new COVID-19 cases continue to accrue, underscoring the need for persistent social distancing and active case detection and isolation efforts to bring epidemic under control.
Author summary In context of the ongoing COVID-19 pandemic, Chile is one of the hardest hit countries in Latin America, struggling to contain the spread of the virus. In this manuscript we employ renewal equation to estimate the reproduction number for the early ascending phase of the COVID-19 epidemic and the most recent time period to guide the magnitude and intensity of the interventions required to combat the COVID-19 epidemic. We also generate short terms forecasts based on the epidemic trajectory using phenomenological models and assess counterfactual scenarios to understand any additional resources required to contain the spread of virus. Our results indicate early sustained transmission of SARS-CoV-2. However, the initial control measures at the start of the epidemic significantly slowed down the spread of the virus whose effect is visible two weeks after the implementation of interventions. Easing of the COVID-19 restrictions in May led to a new wave of infections, followed by the re-imposition of lockdowns in Santiago and other municipalities. While the broad scale social distancing interventions have slowed the most recent spread of the virus spread, the number of new COVID-19 cases continue to accrue, underscoring the need for persistent social distancing efforts to bring epidemic under control.
Introduction
The coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was declared a global pandemic by the World Health Organization (WHO) on March 11th, 2020 [1]. This highly contagious unprecedented virus has impacted government and public institutions, strained the health care systems, restricted people in their homes and caused country-wide lock downs resulting in a global economic crisis. Moreover, as of July 7th, 2020, more than 11.5 million COVID-19 cases in 213 countries and territories have been reported, including more than 530,000 deaths [2]. The social, economic and psychological impact of this pandemic on much of the world’s population is profound [3].
Soon after its initial rapid spread in China, the first case of novel coronavirus beyond China was reported in Thailand on January 13th, 2020 [4]. The first case in the USA was not identified until January 20th, 2020 followed by the detection of the first cases in the European territory on January 24th, 2020 [5, 6]. The COVID-19 pandemic has since spread to every continent except the Antarctica. While some countries like New Zealand, and the Republic of Korea have steadily suppressed the COVID-19 spread, reporting < 70 cases per day as of July 7th, 2020, other countries like Brazil, India and the USA still struggle to contain the increasing number of cases [7]. Subsequently, considerable COVID-19 outbreaks have occurred in Latin America since late February 2020, with 1603055 cases reported in Brazil, followed by 302718 cases in Peru and 301019 cases in Chile as of July 7th, 2020 [7, 8]. The first case of SARS-CoV-2 in Chile was identified on March 3rd, 2020. While the initial cases were imported from Southeast Asia and Europe, the COVID-19 case counts have expanded in this country, placing Chile in phase 4 of the pandemic on March 25th, 2020 [9, 10].
Chile is the fifth country in Latin America after Brazil, Mexico, Ecuador and Argentina to report COVID-19 cases. The first six imported cases were reported in the city of Talca and the capital of Chile, Santiago [10]. However, since the early phase of the outbreak, Chile has employed an agile public health response by announcing a ban on large public health gatherings of more than 500 people on March 13th, 2020 when the nationwide cumulative case count reached 44 reported cases [11]. Moreover, the Chilean government announced the closure of all day cares, schools and universities on March 16th, 2020. This was followed by the announcement to close country borders on March 18th, 2020, and the declaration of national emergency on the same date, which was accompanied with several concrete interventions to further contain the outbreak in the region [12]. In particular, these included a night time curfew in Chile starting on March 22nd, 2020, and small-area lockdowns (i.e., intermittent lockdowns at the municipality level depending on total cases and case growth) starting on March 28th, 2020 in two municipalities in Southern Chile and seven municipalities in Santiago [13]. These initial containment strategies kept the COVID-19 case counts lower than regional peers; Brazil, Peru and Ecuador until the end of April, 2020. However, the government started to ease the COVID-19 restrictions in early May by reopening the economy under the “Safe Return” plan, including televised opening of some businesses and stores, as new infections had reduced between 350-500 per day by the end of April, implying an only apparent flattening of the COVID-19 curve [14–16]. Moreover, the approach of imposing and lifting the small-area lockdowns by area proved unsuccessful, resulting in a new wave of infections; with the virus spreading from more affluent areas of Chile to poorer, more crowded communities, forcing the government to reimpose lockdown measures in Santiago in mid-May [17]. The case counts have since continued to increase averaging ~4943 cases per day in June 2020. Thus far, Chile has accumulated a total of 301019 reported cases including 6434 deaths as of July 7th, 2020. The majority (~77%) of COVID-19 cases are concentrated in Chile’s capital, Santiago, with 231171 reported cases, followed by Valparaiso (13433) in coastal central Chile, and Antofagasta in Northern Chile (10594) [18]. However, the crude case fatality rate in Chile (2.1%) remains below the global average (4.6%) [7, 19].
In this study we estimate the transmission potential of COVID-19, including the effective reproduction number, R, during the early transmission phase of the COVID-19 epidemic in Chile that can guide the magnitude and intensity of control interventions to combat the epidemic [20, 21]. We can also study the effectiveness of control interventions in Chile (see Table 1) on the transmission rate by conducting short-term forecasts calibrated using the early trajectory of the epidemic [22]. We also generate the reproduction number and timely short term forecasts from the COVID-19 reported cases for the most recent time period, as of July 7th, 2020, to anticipate any additional resources required to combat the growing epidemic in Chile.
Methods
COVID-19 incidence and testing data
We obtained updates on the daily series of COVID-19 cases as of July 7th, 2020 from the publicly available data from the GitHub repository created by the government of Chile [8]. Incidence case data by the date of reporting, confirmed by PCR (polymerase chain reaction) tests from March 3rd–July 7th, 2020 were analyzed. The daily testing and positivity rates available from April 9th–July 7th, 2020 were also analyzed.
Generalized growth model (GGM)
We generate short term forecasts using the generalized growth model (GGM) that characterizes the early ascending phase of the epidemic by estimating two parameters: (1) the intrinsic growth rate, r; and (2) a dimensionless “deceleration of growth” parameter, p. This model allows to capture a range of epidemic growth profiles by modulating parameter p. The GGM model is given by the following differential equation:
In this equation C′(t) describes the incidence curve over time t, C(t) describes the cumulative number of cases at time t and p∈[0,1] is a “deceleration of growth” parameter. This equation becomes constant incidence over time if p=0 and an exponential growth model for cumulative cases if p=1. Whereas if p is in the range 0< p <1, then the model indicates sub-exponential growth dynamics [23, 24].
Generalized logistic growth model (GLM)
The generalized logistic growth model (GLM) is an extension of the simple logistic growth model that allows capture of a range of epidemic growth profiles including sub-exponential (polynomial) and exponential growth dynamics. GLM characterizes epidemic growth by estimating (i) the intrinsic growth rate, r (ii) a dimensionless “deceleration of growth” parameter, p and (iii) the final epidemic size, k0. The deceleration parameter modulates the epidemic growth patterns including the sub-exponential growth (0< p <1), constant incidence (p =0) and exponential growth dynamics (p =1). The GLM model is given by the following differential equation:
Where describes the incidence over time t and the cumulative number of cases at time t is given by C(t) [22].
Calibration of the GGM and GLM model
We calibrate the GGM and the GLM model to the daily incidence curve by dates of reporting in Chile using time series data that is available from March 3rd–March 30th, 2020 and from May 9th – July 7th, 2020, respectively.
Model parameters are estimated by a non-linear least square fitting of model solution to the incidence data by the date of reporting. This is achieved by searching for the set of model parameters that minimizes the sum of squared differences between the observed data yti = yt1, yt2, … ytn and the corresponding mean incidence curve given by : where corresponds to the set of parameters of the GGM model and to the set of parameters of the GLM model. In both cases, the objective function for the best fit solution of is given by : where ti is the time stamp at which the time series data are observed and n is the total number of data points available for inference. The initial condition is fixed to the first observation in the data set. This way, gives the best fit to the time series data Next, we utilize a parametric bootstrapping approach assuming negative binomial error structure for the GGM and GLM model to derive uncertainty in the parameters obtained by non-linear least square fit of the data as previously described [23, 25]. The variance is assumed to be 3 times the mean for GGM and 96 times of the mean for the GLM. The model confidence intervals of parameters and the 95% prediction intervals of model fit are also obtained using the parametric bootstrap approach [23].
Reproduction number, R, from case incidence using GGM
The reproduction number, R, is defined as the average number of secondary cases generated by a primary case at time t during the outbreak. This is a crucial quantity to identify the intensity of interventions required to contain an epidemic [26–28]. Estimates of effective R indicate if the disease transmission continues (R>1) or if the active disease transmission ceases (R<1). Therefore, in order to contain an outbreak, we need to maintain R<1. We estimate the reproduction number by calibrating the GGM to the early growth phase of the epidemic (27 days) [24]. We model the generation interval of SARS-CoV-2 assuming gamma distribution with a mean of 5.2 days and a standard deviation of 1.72 days [29]. We estimate the growth rate parameter, r, and the deceleration of growth parameter, p, as described above. Then we simulate the progression of local incidence cases Ii at calendar time ti, from the calibrated GGM model and apply the discretized probability distribution of the generation interval denoted by ρi to the renewal equation to estimate the reproduction number [20, 21, 30]:
The numerator represents the total new cases Ii, and the denominator represents the total number of cases that contribute to generating the new cases Ii at time ti. This way, Rt, represents the average number of secondary cases generated by a single case at time t. Next, we derive the uncertainty bounds around the curve of Rt directly from the uncertainty associated with the parameter estimates (r, p). We estimate Rt for 300 simulated curves assuming a negative binomial error structure where variance is assumed to be 3 times the mean [23].
Reproduction number, R, from case incidence using GLM
We also estimate the most recent reproduction number by calibrating the GLM from May 9th – July 7th, 2020 [24]. Next we model the generation interval [29], estimate the model parameters (r, p, k0) from GLM model and the reproduction number from the renewal equation as described above [20, 21, 30]. The uncertainty bounds around the curve of Rt are derived directly from the uncertainty associated with the parameter estimates (r, p, k0). We estimate Rt for 300 simulated curves assuming a negative binomial error structure [23] where the variance is assumed to be 96 times of the mean calculated by averaging mean to variance ratio calculated from the data (by binning data points and calculating directly from the data itself).
3. Results
Case incidence data
A total of 301019 cases were reported by the Ministry of Health Chile as of July 7th, 2020 [8]. On average ~443 (SD: 133.6) new cases per day were reported in April 2020, ~2697 (SD:1342) new cases per day were reported in May 2020, ~4943 (SD:972.2) new cases per day were reported in June 2020 and ~3163 (SD:570) new cases per day have been reported from July 1st-7th, 2020. Figure 1 shows the daily incidence data of all confirmed cases in Chile as of July 7th, 2020.
Initial growth dynamics and estimate of the reproduction number using GGM
We estimate the reproduction number for the first 27 epidemic days incorporating the effects of the social distancing interventions as explained in Table 1. The incidence curve displays sub-exponential growth dynamics with the scaling of growth parameter, p, estimated at 0.77 (95% CI: 0.73, 0.81) and the intrinsic growth rate, r, estimated at 0.81 (95% CI: 0.67, 1.0). The reproduction number for the early transmission phase was estimated at 1.8 (95% CI: 1.6, 1.9) (Figure 2).
Recent growth dynamics and estimate of reproduction number using GLM
We also estimate the reproduction number from May 9th-July 7th, 2020 incorporating the effects of the reimplementation of lockdowns in Santiago, Antofagasta and Valparaíso. The incidence curve displays nearly linear growth trend with the scaling of growth parameter, p, estimated at 0.51 (95% CI: 0.47, 0.56), the intrinsic growth rate, r, estimated at 22 (95% CI: 13, 31) and the final epidemic size, k0, estimated at 3.4 e+05 (95% CI: 3.1 e+05, 3.7 e+05). The reproduction number was estimated at 0.87 (95% CI: 0.84, 0.89) as of July 7th, 2020 (Figure 3).
Assessing the impact of social distancing interventions
To assess the impact of social distancing interventions in Chile given in Table 1, we generated a 20-day ahead forecast for Chile based on the daily incidence curve until March 30th, 2020. The 28-day calibration period of the GGM model yields an estimated growth rate, r, at 0.8 (95% CI: 0.6, 1.0) and a scaling of growth rate parameter, p, at 0.8 (95% CI: 0.7, 0.8). The 20-day ahead forecast suggested that the early social distancing measures significantly slowed down the early spread of the virus in Chile, whose effect is noticeable about two weeks after the implementation of an intervention as shown in Figure 4. A case resurgence was observed in Chile in mid-May 2020. As a consequence of this case resurgence a total lockdown was imposed in Santiago (representing ~77% of total COVID-19 cases) on May 15th, 2020, which is still in place [12]. We generated a 20-day ahead forecast based on daily incidence curve from May 9th-July 7th, 2020. The 60 day calibration of the GLM model yields an estimated scaling of the growth parameter, p, at 0.52 (95% CI: 0.47, 0.57) representing an almost linear growth pattern. The 20-day ahead average forecast utilizing the GLM model shows that Chile will accumulate ~46677 cases (95% PI: 28746-70071) within next 20-days (Figure 5).
COVID-19 Testing rates and positivity rate
Daily testing and positivity rates for the time period April 9th–July 7th, 2020 by the date of reporting are shown in Figure 6. The total number of tests performed for this time period were 1149770, amongst which 263995 had positive results. The average number of tests performed daily was estimated at ~5460 for April 2020 and ~12959 for May 2020, a 137% increase. The testing rate in Chile further increased in June 2020, testing on average ~17578 individuals per day. The positivity rate (percentage of positive tests among the total number of tests) has fluctuated from an average of ~9.07% (SD: 2.3) in April 2020 to an average of ~28.3% (SD: 4.9) in June 2020.
4. Discussion
The estimates of the early transmission potential in Chile for the first 27 days of the epidemic indicates sustained local transmission in the country with the estimate of reproduction number R at ~1.8 (95% CI: 1.6, 1.9) which is in accordance with the estimates of the reproduction numbers retrieved from studies conducted in China, Brazil, Korea, Peru, South Africa and Iran that lie in the range of 1.5-7.1 [31–40]. In contrast, even lower estimates of R (R <1) that have been reported in Singapore and Australia can be correlated with the implementation of early strict social distancing interventions in these countries [41, 42].
The initial scaling of the growth parameter in Chile indicates a sub-exponential growth pattern (p~0.8), consistent with sub-exponential growth patterns of COVID-19 that have been observed in Singapore (p~0 7), Korea (p~0 76) and other Chinese provinces excluding Hubei (p~0 67) [36, 41, 43]. In contrast, some recent studies have reported nearly exponential growth pattern of the COVID-19 in Iran and Peru and an exponential growth pattern in China [33, 40, 43].
Although the initial transmission stage of COVID-19 in Chile has been attributed to multiple case importations, Chile was quick to implement control measures against the COVID-19 epidemic including border closures on March 18th, 2020 to prevent further case importations. The 20-day ahead forecast of our GGM model calibrated to 28 days suggest that the social distancing measures including closure of schools, universities and day cares have helped slow down the early virus spread in the country by reducing population mobility [44]. The commixture of interventions including small-area lockdowns, night time curfew and ban on social gatherings in Chile can be attributed to preventing the disease trend from growing exponentially in the early growth phase. However, the significant increase in case incidence observed in mid-May can be attributed to the phased reopening of the country and easing of distancing restrictions in early May, in the context of the “Safe Return” plan [12]. As the virus reached the lower income neighborhoods in Chile, the pandemic quickly exploded [45]. While, the COVID-19 case incidence exhibited a relative stabilization in case trajectory for April 2020 (with an average of ~443 cases per day), highlighting the positive effects of early quarantine and lockdowns, measured reopening of the country’s economy in early May 2020 resulted in the surge of cases that continue to accrue till date. The total lockdown comprised of stay at home orders imposed in Santiago (which accounts for about 77% of cases in the country) on May 15th is now showing an effect in slowing the transmission of the virus. Similar lockdowns were imposed in Antofagasta on May 5th and in Valparaíso on June 12th, though these regions together represent only ~8% of cases in Chile [12]. The scaling of growth parameter, p, has been estimated at ~0.51 (95% CI: 0.47, 0.56) after the re-imposition of lockdowns and social distancing measures in May, consistent with a linear incidence growth trend. Moreover, the reproduction number is estimated at ~0.87 (95% CI: 0.84, 0.89), indicating a decline in transmission of the virus consistent with the stay at home orders. The 20-day ahead forecast calibrating data to the GLM model (from May 9th-July 7th, 2020) is reasonable indicating a declining trend in case incidence and implies that Chile could accumulate between 28746-70071 cases from July 8th-July 27th, 2020. Therefore, implementation of prolonged quarantine and continued reduction in population mobility is required to contain the spread of the virus [46].
Though the number of reported cases in Latin America remains low compared to the USA, official data for many Latin American countries are incomplete. However, Chile has tested a higher percentage of its residents than any other Latin American nation lending confidence to its reliability [47]. For instance, the average number of COVID-19 tests performed in Chile per day per thousand people is 0.75 compared to the neighboring South American country, Peru (~0.09 tests per thousand people) [48]. The average positivity rate of COVID-19 in Chile is estimated at ~28.3% for June, 2020 and ~20% for the whole span of the epidemic. This high positivity rate indicates that the government has failed to cast a wide enough net to test the masses, and there are probably many more active cases than those detected by epidemiological surveillance, underestimating the epidemic growth curve [49–51]. This positivity rate for Chile is also consistent with the high positivity rates derived from Mexico and Nigeria (20-50%) [48]. In comparison countries like Denmark, Germany, Canada, New Zealand, Australia and South Korea have reached very low positivity rates (0-3%) with South Korea’s large testing capacity combined with a strategy that tracks infected people via cell phones [47, 48, 52]. Other countries like India, US and Singapore exhibit moderately high positivity rates (4-8%) for COVID-19, indicating an overall limited testing in these countries [48, 53]. Moreover, studies suggest there is asymptomatic transmission of SARS-CoV-2, which means we could have underestimated our estimates based on the daily incidence’s growth trend from symptomatic cases [54–56]. On the other hand, preliminary results of a study have shown the relative transmission of asymptomatic cases in Santiago to be almost ~3% [57]. While our study highlights the effectiveness of broad scale social distancing and control interventions in Chile, it also underscores the need for persistent isolation and social distancing measures to stomp all active disease transmission chains in Chile.
This study has some limitations. First our study analyzes cases by the dates of reporting while it is ideal to analyze the cases by the dates of onset or after adjusting for reporting delays. On the other hand, a substantial fraction of the COVID-19 infections exhibit very mild or no symptoms at all, which may not be reflected by data [58]. Second, the data is not stratified by local vs. imported cases, therefore, we assume that all cases contribute equally to the transmission dynamics of COVID-19. Finally, the extent of selective underreporting, and its impact on these results, is difficult to assess.
5. Conclusions
In this study we estimate the transmission potential of SARS-CoV-2 in Chile. Our current findings point to sustained transmission of SARS-CoV-2 in the early phase of the outbreak, with our estimate of the reproduction number at R~1.8. The COVID-19 epidemic in Chile followed an early sub-exponential growth trend (p ~0.8) which has transformed into an almost linear growth trend (p ~0.5) with the most recent estimate of reproduction number (R~0.8) as of July 7th, 2020 indicating a decline in the virus transmission in Chile. The implementation of lockdowns and apt social distancing interventions have indeed slowed the spread of the virus. However, the number of new COVID-19 cases continue to accumulate, underscoring the need for persistent social distancing and active case finding efforts to bring the epidemic under control.
Data Availability
All data are publicly available.
Author Contributions
Conceptualization, G.C. and A.T.; methodology, G.C, A.T.; validation, G.C.; formal analysis, A.T., G.C.; investigation, A.T.; resources, G.C.,A.T.; data curation, A.T.; writing— original draft preparation, A.T., G.C.; writing, review and editing, A.T., G.C., E.U, C.C, K.V., R.R, R.L, G.C.; visualization, A.T., G.C.; supervision, G.C.; project administration, G.C.; funding acquisition, G.C. All authors have read and agreed to the published version of the manuscript.
Funding
G.C. is partially supported from NSF grants 1610429 and 1633381 and R01 GM 130900.
Conflicts of Interest
The authors declare no conflict of interest.