The effectiveness of early start of Grade III response to dengue in Guangzhou, China: A population-based interrupted time-series study.

In 2019, dengue incidences increased dramatically in many countries. However, the prospective growth in dengue incidence did not occur in Guangzhou, China. We examined the effectiveness of early start of Grade III response to dengue in Guangzhou. We extracted the data on daily number of dengue cases during 2017-2019 in Guangzhou and weekly data for Foshan and Zhongshan from the China National Notifiable Disease Reporting System, while the data on weekly number of positive ovitraps for adult and larval Aedes albopictus were obtained from Guangzhou Center for Disease Control and Prevention. We estimated the number of dengue cases prevented by bringing forward the starting time of Grade III response from September in 2017-2018 to August in 2019 in Guangzhou using a quasi-Poisson regression model and applied the Baron and Kenny's approach to explore whether mosquito vector density was a mediator of the protective benefit. In Guangzhou, early start of Grade III response was associated with a decline in dengue incidence (relative risk [RR]: 0.54, 95% confidence interval [CI]: 0.43-0.70), with 987 (95% CI: 521-1,593) cases averted in 2019. The rate of positive ovitraps also significantly declined (RR: 0.64, 95% CI: 0.53-0.77). Moreover, both mosquito vector density and early start of Grade III response was significantly associated with dengue incidence after adjustment for each other. By comparing with the incidence in Foshan and Zhongshan where the Grade III response has not been taken, benefits from the response starting in August were confirmed but not if starting from September. Early start of Grade III response has effectively mitigated the dengue burden in Guangzhou, China, which might be partially through reducing the mosquito vector density. Our findings have important public health implications for development and implementation of dengue control interventions for Guangzhou and other locations with dengue epidemics.

Introduction

Dengue is the most rapidly spreading mosquito-borne disease throughout the world and it is currently endemic in 128 countries [1]. Globally, around half of the population are at risk of dengue virus infection, with approximately 400 million people infected by dengue viruses annually [2]. Climate change, population growth, increasing human mobility, and unplanned urbanization are expected to facilitate the spread of dengue viruses and thereby aggravate the dengue disease burden [3-6]. In 2019, World health Organization (WHO) listed dengue virus as the world’s top 10 public health threats [1].

Vaccination is an intervention for the prevention of dengue infection. However, dengue vaccination programs are only recommended for areas with a high dengue burden, considering the vaccine safety issues for seronegative individuals and low vaccine efficacy [7-9]. Even as effective vaccines become available, vector control is crucial in dengue prevention and mitigation [10, 11]. Previous studies have reported the effectiveness of interventions which included vector control measures in locations such as two Chinese cities of Guangzhou [12, 13] and Ningbo [14], Sri Lanka [15], and Singapore [16]. In addition, Liu et al. illustrated that the implementation of early rigorous interventions could reduce dengue epidemic through simulating dengue incidences under different scenarios [17]. However, there is a lack of data on comparing the effectiveness of interventions commenced at different times for the identification of appropriate start time of the intervention.

Dengue disease burden distributes disproportionately globally, with Asia representing 75% of the burden [18]. Guangzhou, located in the subtropical climatic zone (latitude: 23°07′N; longitude 113°15′E) and having a permanent population of 14.9 million, is most highly influenced by dengue in China. During 2006–2014, approximately 70% of the dengue cases in China were reported in Guangzhou [19]. During 2017–2018, the peak dengue season was from August to November. Guangzhou Municipal Health and Family Planning Commission (GZ MHFPC) divided the events related to dengue into six grades and formulated a series of responses for events of different grades. Considering the dengue epidemic, funding and manpower required for the initiation, the local government started the Grade III response which targeted the Grade III events related to dengue in September in 2017 and 2018 to curb the proliferation of dengue transmission in the population. However, the initiation of Grade III response mainly triggered by the event may be not timely for adequately controlling the consequent epidemics.

In 2019, large dengue outbreaks occurred in many locations and especially, places such as Bangladesh, Pakistan, and WHO Region of the Americas witnessed their worst dengue outbreaks in recent years [1, 20–22]. Due to the escalation of outbreaks, WHO scaled up the global response to dengue [1]. In Guangzhou, a higher mosquito vector density in early months was a potential signal of a large dengue epidemic in the year [23]. In response to the high risk of increasing severity of dengue epidemic, the start time of Grade III response was advanced in August. This provided us a precious opportunity to examine the impact of early implementation of control program on the additional reduction in dengue disease burden. The present study aims to 1) gauge the effect of early start of Grade III response on dengue incidence and 2) explore the role of vector control in dengue prevention and mitigation, when the beginning time of Grade III response was advanced one month ahead in 2019 in Guangzhou, China.

Methods

Data source

We extracted the data on daily number of dengue cases confirmed by serological and/or virologic methods in Guangzhou during 2017–2019 and weekly data for Foshan and Zhongshan from the China National Notifiable Disease Reporting System. The dengue cases were classified into local cases and imported cases. The mosq-ovitrap is the most widely employed device for the surveillance of Aedes albopictus [24, 25] and it is currently used for routine surveillance of container-breeding Aedes in 172 streets in all of the 11 districts of Guangzhou. Data on the number of positive ovitraps for adult and larval Aedes albopictus and total number of retrieved traps were obtained from Guangzhou Center for Disease Control and Prevention (GZ CDC). Data collection frequency was determined according to dengue epidemic, from once a month between January and March, to once a week between July and November. Mosquito ovitrap index (MOI) was calculated as the number of positive ovitraps for adult and larval Aedes albopictus per 100 traps which were retrieved [26]. Daily meteorological data on ambient mean temperature, relative humidity, and rainfall were downloaded from the China Meteorological Data Sharing Service System [27].

The Grade III response to dengue

GZ MHFPC formulated a series of responses for the events related to dengue, which were divided into six grades (i.e., VI, V, IV, III, II, I) according to the severity of events, with the Grade I events being the most severe ones (S1 Text) [28]. Generally speaking, the responses is comprised mainly of vector surveillance, vector control, raising public awareness of the health risk related to mosquito vectors, public announcement of the outbreak. The most important counterpart of these responses is vector control, including killing adult mosquito vectors, clearing indoor and outdoor mosquito breeding grounds. The core area is defined as the area within a radius of 200 meters from a dengue case, while the warning zone is defined as the area within a radius of 400 meters but outside a radius of 200 meters of a dengue case. Streets are required to kill adult mosquito vectors indoors and outdoors in core areas within three days. Subsequently, the intervention is implemented every three days for three times and then once a week until the end of the epidemic. The intervention is meanwhile conducted in warning zones with different frequency (i.e. once a week if mosquito vector density is reduced to the safety level; every three days if mosquito vector density still exceeds the safety level. Details on the safety level are provided in S2 Text).

Vector control requires efforts from government, communities and individuals. Interventions formulated for the events of six grades related to dengue are in general consistent, but the degree of government involvement differs. The responses for events of Grade VI-IV rely mainly on the efforts from local health authority with limited resources. However, local government dominates the work of dengue control when Grade III-I responses start. The aims of the Grade III response, which targets the Grade III events related to dengue (S1 Text), are quickly controlling dengue outbreak and preventing rapid spread of dengue in the population. The initiation of Grade III response is influenced by various factors, such as dengue epidemic, funding and manpower required for the initiation, national and global events occurring at the time. After commencing the Grade III response, the intervention for dengue control is intensified, additional works are required to be conducted, more resources are allocated and the cross-departmental collaboration is reinforced, with direct supervision of senior officials.

Data analysis

Since the data on mosquito vector density were collected at the end of month during January and March, the weekly data were available only for the last weeks of these months, while the values of mosquito vector density were missing in other weeks of these months and the missing values were not imputed. The values of rainfall were highly right skewed (more than half of the values were zeroes) and there were a few outlier values (S1 Fig). The inclusion of rainfall as a continuous variable in the model would impair the inference for the association between rainfall and dengue incidence. Thus, we categorized rainfall into three levels (i.e. 0, 0.05–3.85, >3.85mm; 0.05 and 3.85mm were the minimum and median of the non-zero values of rainfall). We compared the effect of Grade III response commencing in August 2019 on local dengue incidence with the counterfactual scenario that the Grade III response started in September as usual by associating an indicator variable of the Grade III response in 2019 and daily incidence of dengue in Guangzhou using a quasi-Poisson regression model [15]. The model evaluating the effect of early start of Grade III response on dengue incidence (Model 1) was as follows: where t is sequential days from January 1, 2017 to December 31, 2019; D and Pop are daily number of local dengue cases and population at the time point t, respectively. We included the logarithm transformation of population into the model and constrained the corresponding regression coefficient to be one. A natural cubic spline function ns(.) with four degrees of freedom (dfs) was used for the calendar day in each year to control for the seasonality of dengue incidence. In addition, a linear function of calendar year was also included. Candidate climatic variables included daily mean temperature, relative humidity, and a categorical variable of rainfall. We selected climatic variables and corresponding maximum lags according to the Quasi Akaike Information Criterion (QAIC). In Model 1, Temp and Hum are metrics constructed by applying the distributed lag non-linear model (DLNM) to temperature and relative humidity, with maximum lags of 22 and 69 days, respectively (S1 Table). We constrained the time lag for the indicator variable of the Grade III response in 2019 (1: between August 9 and December 18, 2019; 0: otherwise) to be at least 10 days, considering the intrinsic incubation period (IIP) and time delay from illness onset to reporting cases, and selected the maximum lag based on QAIC. The cross-basis of an indicator variable of the Grade III response in 2019 and time lag (Int) was subsequently included in the model (S2 Table). The cross-basis can indicate the additional effect of early start of Grade III response, since the Grade III response started in September in 2017 and 2018, but in August in 2019 (S3 Text).

We hypothesized that the Grade III response starting one month ahead in August could influence local dengue incidence through affecting mosquito vector density. To explore the mediator effect of mosquito vector density, quasi-Poisson regression models were further applied to evaluate the effect of early start of Grade III response on the rate of positive ovitraps for adult and larval Aedes albopictus (Model 2) and to associate the indicator variable of the Grade III response in 2019, logarithm transformation of MOI, and weekly number of local dengue cases (Model 3), according to the prerequisites for establishing mediation proposed by Baron and Kenny [29, 30] (S2 Fig). Model 2 and Model 3 were as follows: where tw is sequential weeks from 1 to 156; Pos and Trap are the number of positive ovitraps for adult and larval Aedes albopictus and total number of traps at the time point tw, respectively; Week represents the calendar week in each year; D and Pop are the number of local dengue cases and population at the time point tw, respectively. According to Model 1, moving averages of 0–3 and 0–10 weeks were applied to temperature and relative humidity for Model 3, respectively. The time lags between climatic variables and the positive rate in Model 2 and the lag between the logarithm transformation of MOI and dengue incidence in Model 3 were selected jointly, assuming that climatic factors could influence dengue incidence through mosquito vector density. The current-week temperature and a moving average of relative humidity at a lag of 0–7 weeks were included in Model 2, meanwhile the logarithm transformation of MOI at a lag of 1–3 weeks were included in Model 3. In addition, to describe the effects of early start of Grade III response on mosquito vector density and weekly dengue incidence, the indicator variable of the Grade III response in 2019 of current week and at a lag of three weeks were included in Model 2 and Model 3, respectively according to QAIC (S4 Text and S3 Fig).

Subsequently, we used Foshan and Zhongshan, where Grade III response did not initiate during 2017–2019, as negative controls to examine the effects of Grade III response respectively starting in August and September on dengue incidence. The model (Model 4) was as follow: where Intm is a categorical variable indicating the Grade III response starting in August (August 9–December 18, 2019) and starting in September (September 13–December 19, 2017 and September 14–December 17, 2018) in Guangzhou and no Grade III response in Zhongshan and Foshan during the whole study period; GDP and Green represent gross domestic product and areas of green spaces per capita (in parks); α is a random-effect intercept for three cities.

The reduction in daily number of dengue cases due to early start of Grade III response was estimated by subtracting the predicted number of dengue cases given the Grade III response commenced in August 2019 from that under the counterfactual scenario that the Grade III response started in September 2019. Subsequently, the total number of dengue cases averted was estimated by summing up the daily number of avoided cases during the rest period of 2019 after the Grade III response. To investigate the appropriateness of the models used in the main analysis, histograms of residuals, scatter plots of residuals versus predicted values of response variables, and plots of partial autocorrelation function (PACF) of residuals of the fitted models are presented.

We conducted sensitivity analysis to assess the cumulative relative risk (RR) of dengue over 0–21 days associated with early start of Grade III response, checking whether the effect of early start of Grade III response began within 10 days; to check the robustness of the effect estimates of early start of Grade III response and MOI, by changing the parameters and adjusting for the number of imported cases in the models (S5 Text). We conducted all analyses using R software version 3.6.2 (R Foundation for Statistical Computing).

Results

During 2017–2019, there were 3,871 dengue cases reported in Guangzhou, among which 3,454 (89.2%) were local cases (2017: 873; 2018: 1,199; 2019: 1,382), while 417 (10.8%) were imported cases. The first local case in each year under study was consistently identified in June, after which the number of dengue cases increased until October and then declined (Fig 1 and S4 Fig). MOI, temperature, relative humidity, and rainfall during the months between January and June were on average higher in 2019 than in 2017 and 2018 (S4 Fig).

Fig 1

Time-series of daily number of reported dengue cases, weekly mosquito ovitrap index (MOI), daily temperature, relative humidity, and rainfall in Guangzhou, 2017–2019.

Grey regions indicate the time periods of Grade III response.

It was observed that dengue incidence in general increased with temperature and relative humidity, although the associations were not linear (S5 Fig). We estimated that the RR of dengue associated with early start of Grade III response in 2019 was 0.54 (95% confidence interval [CI]: 0.43–0.70) over 10–21 days (Fig 2). Accordingly, bringing forward the start time of Grade III response in 2019 prevented 987 (95% CI: 521–1,593) dengue cases (Fig 3). Early start of Grade III response was also associated with a reduction in the positive rate of ovitraps, with a RR of 0.64 (95% CI: 0.53–0.77). Mosquito vector density might be a mediator between early start of Grade III response and dengue incidence, given (1) the statistically significant effects of early start of Grade III response on dengue incidence and the positive rate of ovitraps, indicated in Model 1 and Model 2, respectively; and (2) the association between the logarithm transformation of MOI and dengue incidence after controlling for the indicator variable of the Grade III response revealed in Model 3 (RR: 1.36, 95% CI: 1.01–1.84) (Table 1). The effect of early start of Grade III response on dengue incidence reduced slightly but remained statistically significant after controlling for mosquito vector density (RR: 0.68, 95% CI: 0.52–0.90). The findings implied that early start of Grade III response may reduce mosquito vector density which might indirectly curb the transmission of dengue viruses. In addition, early start of Grade III response could reduce the number of dengue cases in other ways. Results of Model 4 indicated that the Grade III responses starting in August were associated with a decline in dengue incidence, with RR of 0.59 (95% CI: 0.45–0.78) but the benefit was not statistically significant if the response started in September (RR: 0.91, 95% CI: 0.74–1.12) when comparing with Foshan and Zhongshan where Grade III response was not taken during the whole study period.

Fig 2

Cumulative relative risk of dengue due to bringing forward the starting time of Grade III response from September to August in 2019 in Guangzhou over 10–21 days.

The grey region indicates the corresponding 95% confidence interval.

Fig 3

Observed and predicted numbers of dengue cases with the Grade III response starting in September and August 2019 in Guangzhou.

The grey region indicates the time period of Grade III response in 2019. Blue and red regions represent the subtraction of the predicted number of dengue cases given the Grade III response commenced in August 2019 from the predicted number of dengue cases under the counterfactual scenario that the Grade III response started in September 2019.

Table 1

Associations between the logarithm transformation of MOI, early start of Grade III response, and dengue incidence in Model 3.

Variable	RR	(95% CI)
log (MOI)	1.36	(1.01–1.84)
Early start of Grade III response	0.68	(0.52–0.90)

Abbreviations: RR, relative risk; 95% CI, 95% confidence interval; log(MOI), logarithm transformation of mosquito ovitrap index.

Results of model diagnostics suggested that residuals of models were nearly normally distributed and that there was no significant autocorrelation in residuals (S6 Fig). Sensitivity analysis indicated that the results were robust to the specification of parameters and whether the number of imported cases was adjusted or not (S5 Text, S7 Fig, S3, S4 and S5 Tables).

Discussion

In this study, we assessed the effectiveness of early start of Grade III response in 2019 on dengue control in Guangzhou. In the early months of 2019 between January and June, temperature, relative humidity, MOI, and the number of reported dengue cases in Guangzhou were higher than that in 2017 and 2018 (S4 Fig), signaling that the dengue outbreak was likely to be more severe in 2019 than that in previous two years given the positive associations between these factors and dengue incidence reported in previous studies [23, 31]. However, the prospective increase in dengue incidence did not occur in Guangzhou, which could be attributed to bringing forward the start time of Grade III response from September to August. We estimated that early start of Grade III response in 2019 prevented 987 dengue cases in Guangzhou, which were equivalent to 71.4% of the total number of reported local dengue cases in the year. It was estimated that approximately five million RMB could be additionally consumed due to early start of Grade III response in August, according to the investigation by GZ CDC. However, around 95,600 RMB could be saved annually if one dengue case was prevented in Guangzhou [32]. Based on our estimates of the averted dengue cases in 2019, around 94,357,200 RMB would be saved due to early start of Grade III response. That is, early start of Grade III response could save around 90 million RMB (i.e. 12.89 million USD) in 2019. In addition, results for comparing with no Grade III response suggested that the effect of Grade III response starting in August was much more significant than that commencing in September, which also highlighted the importance of early dengue control. Future studies on comprehensive cost-effectiveness analysis of the Grade III response initiating in August or even earlier before the peak of mosquito vector density occurs can inform suitable public health interventions for dengue prevention and mitigation. The experience in Guangzhou, especially the effective integrated intervention strategy with the cross-departmental collaboration and supervision from senior officials could also be beneficial for other locations with dengue epidemics.

Vector control is of great importance in preventing and mitigating dengue burden [10, 11]. A previous study conducted in Sri Lanka revealed that mosquito vector density could be a mediator of intensified dengue control measures and dengue incidence [15]. Interestingly, we found that mosquito vector density could also be a mediator of early start of Grade III response and dengue incidence. Killing adult mosquito vectors and eliminating mosquito breeding grounds are counterparts of the integrated intervention measures included in the Grade III response to dengue. Therefore, the Grade III response can contribute to a decline in mosquito vector density. Early start of Grade III response reduced the duration of potential exposure to infected mosquito vectors for local residents and subsequently mitigated the dengue burden. Furthermore, it was observed that the effect of bringing forward the start time of Grade III response one month ahead remained statistically significant, after controlling for the effect of the logarithm transformation of MOI (Table 1). The findings implied that early start of Grade III response may have an impact on dengue incidence not only through reducing mosquito vector density, but also in other ways, such as controlling imported dengue cases, strengthening social preparedness and raising public awareness against dengue fever.

Previous studies have showed that dengue incidence was positively associated with household index, container index, and Breteau index in Vietnam [33] and Guangzhou, China [34]. In the present study, we found a positive association between MOI and dengue incidence. Undoubtedly, a higher mosquito vector density could indicate an elevation in the number of mosquitoes infected by dengue viruses. In this situation, the population could be exposed to more infected mosquitoes and therefore at a higher risk of dengue infection. These findings highlighted the importance of reducing mosquito vector density for dengue control. The result of the association between the index of mosquito vector density and dengue incidence can be used for constructing the predictive model for dengue disease burden [4].

The current study had several limitations. First, the number of reported dengue cases was an underestimate of dengue infections in Guangzhou, since almost all reported cases were hospitalized and many asymptomatic people infected with dengue were not identified or reported [2]. Consequently, the number of prevented cases associated with early start of Grade III response would be underestimated. Second, we analyzed 3-year data, while the modification of the start time of Grade III response was observed only for one year. Cumulative data for different start time of response are required to identify an optimal stating time of Grade III response. In addition, because of data availability, the impact of early start of the Grade III response on mosquito vector density and the association between the logarithm transformation of MOI and dengue incidence were assessed on a weekly basis instead of a daily scale, leading to less data points. Next, the intensity of intervention for dengue control might fluctuate during the period of Grade III response. However, the intervention intensity of Grade III response should be within a certain range and should be significantly stronger than that of Grade VI-IV responses given the high degree of government involvement [28]. Since it is extremely difficult to quantify the intervention intensity, we only used a dummy variable indicating the period of Grade III response in 2009 in the model to examine the average effect of early start of Grade III response on dengue incidence. We explored whether mosquito vector density was a mediator by using Baron and Kenny’s approach instead of structural equation modeling (SEM). Although SEM is a powerful technique to distinguish the direct and indirect effects in causal relationships, unfortunately, the techniques currently available for SEM are not applicable to examine the lagged effects for such non-Gaussian time-series data in this study. Further, we considered population, GDP and green coverage when comparing the dengue incidence among Guangzhou, Foshan and Zhongshan, while other variables such as a greater degree of population mobility may lead to higher dengue incidence in Guangzhou, which may be likely to underestimate the impacts of Grade III response.

In conclusion, bringing forward the starting time of Grade III response one month ahead to August have effectively reduced the dengue burden in Guangzhou, China, which might be partially through controlling the mosquito vector density. Our findings have important public health implications on the dengue control interventions for Guangzhou and other locations with dengue epidemics.

Histogram of rainfall.

(TIF)

Click here for additional data file.

Models used to test whether mosquito vector density was a mediator between early start of Grade III response and dengue incidence.

(TIF)

Click here for additional data file.

Combinations of time lags between climatic variables and mosquito vector density in Model 2 and the lag between logarithm transformation of MOI and dengue incidence in Model 3.

Six combinations of time lags (C1-C6) were shown in the dashed line rectangle. H→M, T→M, and M→D represent the time lags between relative humidity and mosquito vector density, between temperature and mosquito vector density, between mosquito vector density and dengue incidence, respectively.

(TIF)

Click here for additional data file.

Time-series of monthly number of reported dengue cases, average mosquito ovitrap index (MOI), temperature, absolute humidity, and rainfall in Guangzhou, China, during 2017–2019.

(TIF)

Click here for additional data file.

Cumulative relative risk of dengue across temperature and relative humidity.

(A) Cumulative relative risk of dengue across temperature; (B) Cumulative relative risk of dengue across relative humidity. Medians of temperature and relative humidity were treated as the reference levels for the two variables. Grey regions indicate the corresponding 95% confidence intervals.

(TIF)

Click here for additional data file.

Residual analysis of the models used in the main analysis.

(A) Histogram of residuals of Model 1; (B) Histogram of residuals of Model 2; (C) Histogram of residuals of Model 3; (D) Histogram of residuals of Model 4; (E) Scatter plot of residuals versus predicted daily number of reported dengue cases of Model 1; (F) Scatter plot of residuals versus predicted number of positive ovitraps for adult and larval Aedes albopictus of Model 2; (G) Scatter plot of residuals versus predicted weekly number of reported dengue cases of Model 3; (H) Scatter plot of residuals versus predicted weekly number of reported dengue cases of Model 4; (I) Partial autocorrelation function (PACF) of residuals of Model 1; (J) Partial autocorrelation function of residuals of Model 2; (K) Partial autocorrelation function of residuals of Model 3; (L) Partial autocorrelation function of residuals of Model 4.

(TIF)

Click here for additional data file.

Relative risk of dengue due to bringing forward the starting time of Grade III response from September to August in 2019 in Guangzhou over 0–21 days.

Grey regions indicate the corresponding 95% confidence intervals.

(TIF)

Click here for additional data file.

Values of Quasi Akaike Information Criterion (QAIC) corresponding to the models which incorporate natural cubic splines of temperature or relative humidity with different degrees of freedom.

(DOCX)

Click here for additional data file.

Values of Quasi Akaike Information Criterion (QAIC) corresponding to the models which incorporate the cross-basis of an indicator variable of the Grade III response in 2019 and time lag with different maximum lags.

(DOCX)

Click here for additional data file.

Estimates of averted number of dengue cases due to early start of Grade III response in the sensitivity analysis.

(DOCX)

Click here for additional data file.

The effect of early start of Grade III response on the rate of positive ovitraps for adult and larval Aedes albopictus in the sensitivity analysis.

(DOCX)

Click here for additional data file.

The effect of the logarithm transformation of mosquito ovitrap index (MOI) on dengue incidence in the sensitivity analysis.

(DOCX)

Click here for additional data file.

Events related to dengue of six grades and the corresponding responses.

(DOCX)

Click here for additional data file.

The safety level of mosquito vector density.

(DOCX)

Click here for additional data file.

Model 1.

(DOCX)

Click here for additional data file.

Model 2 and Model 3.

(DOCX)

Click here for additional data file.

Sensitivity analysis.

(DOCX)

Click here for additional data file.

4 May 2020

Dear Dr. Ou,

Thank you very much for submitting your manuscript "The effectiveness of early start of Grade III response to dengue in Guangzhou, China: a population-based interrupted time-series study" for consideration at PLOS Neglected Tropical Diseases. As with all papers reviewed by the journal, your manuscript was reviewed by members of the editorial board and by several independent reviewers. In light of the reviews (below this email), we would like to invite the resubmission of a significantly-revised version that takes into account the reviewers' comments.

We cannot make any decision about publication until we have seen the revised manuscript and your response to the reviewers' comments. Your revised manuscript is also likely to be sent to reviewers for further evaluation.

When you are ready to resubmit, please upload the following:

[1] A letter containing a detailed list of your responses to the review comments and a description of the changes you have made in the manuscript. Please note while forming your response, if your article is accepted, you may have the opportunity to make the peer review history publicly available. The record will include editor decision letters (with reviews) and your responses to reviewer comments. If eligible, we will contact you to opt in or out.

[2] Two versions of the revised manuscript: one with either highlights or tracked changes denoting where the text has been changed; the other a clean version (uploaded as the manuscript file).

Important additional instructions are given below your reviewer comments.

Please prepare and submit your revised manuscript within 60 days. If you anticipate any delay, please let us know the expected resubmission date by replying to this email. Please note that revised manuscripts received after the 60-day due date may require evaluation and peer review similar to newly submitted manuscripts.

Thank you again for your submission. We hope that our editorial process has been constructive so far, and we welcome your feedback at any time. Please don't hesitate to contact us if you have any questions or comments.

Sincerely,

Elvina Viennet, PhD

Deputy Editor

PLOS Neglected Tropical Diseases

David Harley

Deputy Editor

PLOS Neglected Tropical Diseases

***********************

Reviewer's Responses to Questions

Key Review Criteria Required for Acceptance?

As you describe the new analyses required for acceptance, please consider the following:

Methods

-Are the objectives of the study clearly articulated with a clear testable hypothesis stated?

-Is the study design appropriate to address the stated objectives?

-Is the population clearly described and appropriate for the hypothesis being tested?

-Is the sample size sufficient to ensure adequate power to address the hypothesis being tested?

-Were correct statistical analysis used to support conclusions?

-Are there concerns about ethical or regulatory requirements being met?

Reviewer #1: (No Response)

Reviewer #2: 1. The response for the dengue includes six grades, which grade is chosen based on the severity of event. The question is whether the grade will be upgraded to a higher grade when the cases are increasing during an outbreak. For example, the cases were relatively lower when Grade III was launched in 2019, whether the intervention was strengthened with cases increasing. As the authors mentioned in the paper, the responses are comprised mainly of vector surveillance, vector control, raising public awareness etc., which are not quantified measurements, how to sustain the consistency of the intervention during an outbreak, and the comparability among 2017, 2018 and 2019.

2. The study indicated mosquito density was a mediator of intensified dengue control measures, climatic factors and dengue incidence. Therefore, it may be not appropriate that mosquito density, and climatic factors and dengue control measures were recognized as the independent variables simultaneously to interpret dengue incidence as in model 3, but path analysis or structural equation model can.

3. It’ known that rainfall influences the mosquito density and then influences dengue incidence. I notice that rainfall was not included into the model based on the study’s statistical strategy. Daily data on rainfall were downloaded, why rainfall was arbitrarily categorized into three levels not as the continuous variable was included in the model.

4. Was lag (Intt, 10-21) or lag (Intt, 10-22) included in model 1?

5. The data on mosquito density was collected once a moth from Jan to Mar, but once a week from Jul to Nov. The variable was included into the model using weekly data, how was the monthly data manipulated?

--------------------

Results

-Does the analysis presented match the analysis plan?

-Are the results clearly and completely presented?

-Are the figures (Tables, Images) of sufficient quality for clarity?

Reviewer #1: (No Response)

Reviewer #2: 1. The results of sensitivity analysis should be submitted in the Sup. Materials.

2. The lags of climatic factors were selected based on QAIC. The exposure-response association between climatic factors and dengue incidence should be shown, which is helpful for reader to assess the model, not only using residuals.

--------------------

Conclusions

-Are the conclusions supported by the data presented?

-Are the limitations of analysis clearly described?

-Do the authors discuss how these data can be helpful to advance our understanding of the topic under study?

-Is public health relevance addressed?

Reviewer #1: (No Response)

Reviewer #2: The study used an interrupted time series method to explore the effectiveness of early start of grade III response to dengue. The intervention of response as binary variable, yes or no, was entered into the model. However, i still have concerns on the consistency of the interventions. Responses are separated into six grades, but the core may be the same, what difference may be the intensity. How is the response measured?

--------------------

Editorial and Data Presentation Modifications?

Use this section for editorial suggestions as well as relatively minor modifications of existing data that would enhance clarity. If the only modifications needed are minor and/or editorial, you may wish to recommend “Minor Revision” or “Accept”.

Reviewer #1: (No Response)

Reviewer #2: (No Response)

--------------------

Summary and General Comments

Use this section to provide overall comments, discuss strengths/weaknesses of the study, novelty, significance, general execution and scholarship. You may also include additional comments for the author, including concerns about dual publication, research ethics, or publication ethics. If requesting major revision, please articulate the new experiments that are needed.

Reviewer #1: 1. Up to date, dengue is not endemic in China. So dengue cases consist of imported cases and local cases in China. I recommend the authors divide dengue cases into imported cases and local cases prior to analysis. Grade III response mainly reduce the number of local cases. Of note, the number of local cases not only influenced by Grade III response but also influenced by the number of imported cases.

2. Grade III response can reduce the number of dengue cases. First, Grade III response can reduce mosquito density which can indirect reduce the transmission of dengue viruses. Second, Grade III response may change the activities of humans or improve humans ability of preventing dengue which direct reduce the number of dengue. SEM should be conducted to indentify direct and indirect impacts of Grade III response.

Reviewer #2: The aim of the study is to evaluate the effectiveness of early Grade III response to dengue in the high risk area of China, Guangzhou, which is interesting and is benefit for the local health authority and government. The study did a lot of analyses, and the main results were shown. I think other results should also be submitted in Sup. Materials to support the conclusion.

--------------------

PLOS authors have the option to publish the peer review history of their article (what does this mean?). If published, this will include your full peer review and any attached files.

If you choose “no”, your identity will remain anonymous but your review may still be made public.

Do you want your identity to be public for this peer review? For information about this choice, including consent withdrawal, please see our Privacy Policy.

Reviewer #1: No

Reviewer #2: No

Figure Files:

While revising your submission, please upload your figure files to the Preflight Analysis and Conversion Engine (PACE) digital diagnostic tool, https://pacev2.apexcovantage.com. PACE helps ensure that figures meet PLOS requirements. To use PACE, you must first register as a user. Then, login and navigate to the UPLOAD tab, where you will find detailed instructions on how to use the tool. If you encounter any issues or have any questions when using PACE, please email us at figures@plos.org.

Data Requirements:

Please note that, as a condition of publication, PLOS' data policy requires that you make available all data used to draw the conclusions outlined in your manuscript. Data must be deposited in an appropriate repository, included within the body of the manuscript, or uploaded as supporting information. This includes all numerical values that were used to generate graphs, histograms etc.. For an example see here: http://www.plosbiology.org/article/info%3Adoi%2F10.1371%2Fjournal.pbio.1001908#s5.

Reproducibility:

To enhance the reproducibility of your results, PLOS recommends that you deposit laboratory protocols in protocols.io, where a protocol can be assigned its own identifier (DOI) such that it can be cited independently in the future. For instructions see https://journals.plos.org/plosntds/s/submission-guidelines#loc-methods

26 Jun 2020

Submitted filename: Responses to Reviewers comments.docx

Click here for additional data file.

30 Jun 2020

Dear Dr. Ou,

We are pleased to inform you that your manuscript 'The effectiveness of early start of Grade III response to dengue in Guangzhou, China: a population-based interrupted time-series study' has been provisionally accepted for publication in PLOS Neglected Tropical Diseases.

Before your manuscript can be formally accepted you will need to complete some formatting changes, which you will receive in a follow up email. A member of our team will be in touch with a set of requests.

Please note that your manuscript will not be scheduled for publication until you have made the required changes, so a swift response is appreciated.

IMPORTANT: The editorial review process is now complete. PLOS will only permit corrections to spelling, formatting or significant scientific errors from this point onwards. Requests for major changes, or any which affect the scientific understanding of your work, will cause delays to the publication date of your manuscript.

Should you, your institution's press office or the journal office choose to press release your paper, you will automatically be opted out of early publication. We ask that you notify us now if you or your institution is planning to press release the article. All press must be co-ordinated with PLOS.

Thank you again for supporting Open Access publishing; we are looking forward to publishing your work in PLOS Neglected Tropical Diseases.

Best regards,

Elvina Viennet, PhD

Deputy Editor

PLOS Neglected Tropical Diseases

David Harley

Deputy Editor

PLOS Neglected Tropical Diseases

***********************************************************

29 Jul 2020

Dear Dr. Ou,

We are delighted to inform you that your manuscript, "The effectiveness of early start of Grade III response to dengue in Guangzhou, China: a population-based interrupted time-series study," has been formally accepted for publication in PLOS Neglected Tropical Diseases.

We have now passed your article onto the PLOS Production Department who will complete the rest of the publication process. All authors will receive a confirmation email upon publication.

The corresponding author will soon be receiving a typeset proof for review, to ensure errors have not been introduced during production. Please review the PDF proof of your manuscript carefully, as this is the last chance to correct any scientific or type-setting errors. Please note that major changes, or those which affect the scientific understanding of the work, will likely cause delays to the publication date of your manuscript. Note: Proofs for Front Matter articles (Editorial, Viewpoint, Symposium, Review, etc...) are generated on a different schedule and may not be made available as quickly.

Soon after your final files are uploaded, the early version of your manuscript will be published online unless you opted out of this process. The date of the early version will be your article's publication date. The final article will be published to the same URL, and all versions of the paper will be accessible to readers.

Thank you again for supporting open-access publishing; we are looking forward to publishing your work in PLOS Neglected Tropical Diseases.

Best regards,

Shaden Kamhawi

co-Editor-in-Chief

PLOS Neglected Tropical Diseases

Paul Brindley

co-Editor-in-Chief

PLOS Neglected Tropical Diseases