Next Article in Journal
Coagulopathy of Dengue and COVID-19: Clinical Considerations
Previous Article in Journal
Serological Investigation and Epidemiological Analysis of Bovine Leptospirosis in Egypt
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Epidemiological Characteristics and the Dynamic Transmission Model of Dengue Fever in Zhanjiang City, Guangdong Province in 2018

1
Guangdong Provincial Center for Disease Control and Prevention, Guangzhou 511430, China
2
State Key Laboratory of Molecular Vaccinology and Molecular Diagnostics, School of Public Health, Xiamen University, Xiamen 361102, China
3
School of Public Health, Southern Medical University, Guangzhou 510515, China
4
Zhanjiang Municipal Center for Disease Control and Prevention, Zhanjiang 524037, China
5
Division of Environmental Health Sciences, School of Public Health, University of California, Berkeley, CA 94704, USA
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Trop. Med. Infect. Dis. 2022, 7(9), 209; https://0-doi-org.brum.beds.ac.uk/10.3390/tropicalmed7090209
Submission received: 22 July 2022 / Revised: 14 August 2022 / Accepted: 15 August 2022 / Published: 25 August 2022

Abstract

:
Background: With the progress of urbanization, the mobility of people has gradually increased, which has led to the further spread of dengue fever. This study evaluated the transmissibility of dengue fever within districts and between different districts in Zhanjiang City to provide corresponding advice for cross-regional prevention and control. Methods: A mathematical model of transmission dynamics was developed to explore the transmissibility of the disease and to compare that between different regions. Results: A total of 467 DF cases (6.38 per 100,000 people) were reported in Zhanjiang City in 2018. In the model, without any intervention, the number of simulated cases in this epidemic reached about 950. The dengue fever transmissions between districts varied within and between regions. When the spread of dengue fever from Chikan Districts to other districts was cut off, the number of cases in other districts dropped significantly or even to zero. When the density of mosquitoes in Xiashan District was controlled, the dengue fever epidemic in Xiashan District was found to be significantly alleviated. Conclusions: When there is a dengue outbreak, timely measures can effectively control it from developing into an epidemic. Different prevention and control measures in different districts could efficiently reduce the risk of disease transmission.

1. Introduction

Dengue fever is an important mosquito-borne (mainly by Aedes mosquitoes [1,2,3]) viral disease that is caused by four different serotypes of the dengue virus (DENV 1-4). According to the World Health Organization (WHO), the incidence rate of DF around the world has grown sharply and dengue fever was also listed as one of the potentially threatening diseases in 2019 by the WHO [4]. Dengue could be widespread in tropical and subtropical regions around the world when given favorable temperature, precipitation and humidity conditions [5]. Studies showed that temperature greatly influences Aedes albopictus mosquitoes, which are the vector of dengue fever. Having a minimum temperature above 18 °C is one of the conditions suitable for the Aedes mosquitoes to breed and thus develop a DF epidemic [6].
Guangdong province is located in the south of China and has a subtropical climate, while Zhanjiang City is one of the prefecture-level cities of Guangdong Province. Since the first outbreak in Guangdong Province in 1978, dengue cases have been reported for years [7], and Guangdong has become the region with the highest incidence rate of DF in China since then [8]. Aedes albopictus mosquitoes, which are densely distributed throughout Guangdong province, are considered the main vector [9]. The dengue serotypes detected in Guangzhou were mainly DENV 1 and DENV 2 [10]. Between 2008 and 2018, more than 85% of indigenous DF cases in China occurred in Guangdong Province [11]. In June 2014, an unprecedented DF outbreak occurred in Guangdong Province. A total of 14,594 clinically diagnosed and 11,387 laboratory-confirmed cases were reported, more than any of the DF epidemics in China in the past 20 years [12,13].
Many mathematical models are available for studying dengue fever. Lan Zhou [14] designed a dengue model with vertical transmission and the results showed that if we could control mosquito population growth, dengue could be eliminated. A dengue transmission model that considered time-dependent entomological parameters and assessed seasonal variations in temperature and precipitation during dengue epidemics was created [15]. However, this research focused on the media and ignored the human impact. There was an SIR model that was created to describe the vector-to-host and host-to-vector transmission dynamics of dengue, as well as the optimal time for pesticide use [16]. However, the impact of host mobility on transmission was not considered. G Bhuju et al. used an SEIR-SEI model to explore the effect of dengue viral load on the transmission dynamics of dengue disease through fuzzy behavior [17]. One study explored the role of asymptomatic infection in the epidemiology of the dengue virus through the susceptibility, asymptomatic, infection and recovery (SAIR) model [18]. In Thailand, a national study that revealed the spatial and seasonal patterns of dengue fever at the regional level was conducted and tested whether the peak of dengue fever in Bangkok was ahead of the peak of dengue fever in other provinces in Thailand [19]. However, it has not yet been considered whether there is a mutual transmission of dengue fever in space. This study used the data on dengue fever in Zhanjiang City in 2018 to investigate the mutual transmission of dengue fever between the host and the vector. At the same time, the mobility of hosts between different districts, which causes disease transmission between regions, was taken into consideration.
In order to provide corresponding suggestions for cross-regional prevention and control in Zhanjiang City and to offer methods for studying the trans-regional transmission of DF in other regions, we developed an ordinary differential equation (ODE) model to simulate the transmission dynamics of the dengue virus in this study. This was the first study that used the SEI-SEIAR model to explore the dynamics of dengue fever spreading between districts.

2. Methods

2.1. Materials

In 2018, the first dengue outbreak occurred in Chikan District, followed by outbreaks in Xiashan District, Mazhang District and Kaifa District, all of which are geographically adjacent to Chikan District. Because of the onset time lag of patients between Chikan and the other three districts, we suspected that the disease was transmitted from one region to another. Therefore, we examined the transmissibility between Chikan District and the other three regions. The second indigenous case in Chikan District occurred on 27 August. Since then, the number of cases has gradually increased and peaked on 27 September. At this time, Zhanjiang CDC implemented intervention measures, but the outbreak in Chikan District continued until 6 October. We identified the index case as the first indigenous case in Chikan District. Since there were only two indigenous cases in Potou District, we omitted this district in our analyses. The last indigenous case in the study districts was reported on 28 October in Xiashan District.
We collected data on the daily cases from 17 August to 28 October for the simulation. In addition, we developed a susceptible–exposed–infectious/asymptomatic–removed (SEIAR) model and a susceptible–exposed–infectious (SEI) model for humans and mosquitoes, respectively, to explore the transmissibility of the dengue virus between Chikan district and other districts.

2.2. Study Site

In 2018, Zhanjiang consisted of 5 municipal districts (Chikan, Xiashan, Potou, Mazhang and Kaifa Districts), 3 county-level cities (Wuchuan, Leizhou and Lianjiang Cities) and 2 counties (Xuwen and Suixi Counties). The permanent resident population of Zhanjiang City was 7.33 million. The population of Chikan District, which is the center of Zhanjiang City, was about 298,000 people and the land area was 79 square kilometers. Chikan District (110°20′ to 110°21′ E, 21°14′ to 21°19′ N) is located in the tropical northern margin south of the Tropic of Cancer, with an annual mean temperature of 22.8 to 23.5 °C, annual cumulative precipitation of 1596 mm and annual mean sunshine duration time of 1927.7 h. Chikan District is surrounded by Mazhang, Xiashan, Potou and Kaifa Districts. In 2018, the population of Xiashan District was approximately 436,000 people, Mazhang had approximately 275,000 people and Kaifa had approximately 334,000 people. Aedes mosquitoes are the only dengue vector species in Guangdong Province.

2.3. Data Collection

Data from DF cases and daily vector surveillance were obtained from Zhanjiang Center for Disease Control and Prevention (CDC). General information for each case included gender, age, address and occupation. For a subset of the cases, information about the type of diagnosis (confirmed and clinical), virus serotype, and onset and diagnosis time were also collected.

2.4. Case Definitions

DF cases were diagnosed according to the Dengue Diagnosis Criteria (WS 216—2008, WS 216—2018) issued by the Chinese Ministry of Health as follows:
Clinically diagnosed cases: suspected cases of leucopenia or thrombocytopenia, or suspected cases with a positive serum test for IgG or IgM.
Laboratory confirmed cases: suspected cases of dengue virus (DENV) RNA in serum detected using real-time fluorescent quantitative PCR; or the virus was isolated from the blood, tissue or cerebrospinal fluid of an acutely infected person; or the IgG titer in the recovery stage was 4 times higher than that in the acute stage.
Imported cases were defined as dengue cases who had traveled to dengue-endemic districts within 2 weeks before symptom onset, while local cases were defined as those who did not leave Zhanjiang City within 2 weeks before the onset of dengue [20].
According to the Dengue Diagnosis Criteria (WS 216—2008, WS 216—2018) issued by the Chinese Ministry of Health, a DF outbreak is defined as the occurrence of 3 or more indigenous laboratory-diagnosed dengue cases during a maximum incubation period (14 days) in a relatively concentrated population location (such as a community, neighborhood committee, village, school or another collective unit)

2.5. Vector Surveillance

During the epidemic season, when imported or local cases of dengue fever infection are found, as an important part of the investigation and treatment of the epidemic, emergency surveillance is activated.

2.5.1. Surveillance Area

Core area: With reference to the activity of Aedes mosquito, multiple core areas were divided centered on the infected person’s residence, households adjacent to it, the infected person’s workplace, and other places where the case has been.Alert zone: A radius of 200 m outside the core zone. In rural areas, natural villages and tunnels around the core area, and if necessary, administrative villages or even townships, are set as alert zones. If necessary, administrative villages or even towns or villages should be considered alert zones. In urban areas, the alert area is usually the core area surrounded by some streets, lanes, neighborhood committees or streets.
Surveillance zone: According to the epidemic size and epidemic season of different dengue risk areas, the surveillance zone will be set up outside the vigilance zone. The surveillance area is defined around the alert area.

2.5.2. Monitoring Methods

All dengue mosquito vector emergency monitoring points were monitored using the Breteau index method.

2.5.3. The Frequency of Monitoring

Breteau index method: After 1–2 days of the dengue epidemic, a full coverage survey and emergency mosquito vector control were conducted in the core area, followed by repeated control and survey every 2–3 days until the BI was less than 5. The alert area was surveyed once a week, while the surveillance area was surveyed once every 2 weeks.

2.5.4. Data Analysis and Feedback

Risk assessment: A Breteau index (BI) of less than 5 is the threshold for controlling dengue transmission. A BI greater than 5 has the risk of transmission, greater than 10 has the risk of an outbreak, and greater than 20 has the risk of a regional epidemic, requiring continued cleanup of breeding sites, which should be cleared of adult mosquitoes.
Breteau index (BI) is an index that describes the density of Aedes mosquitoes and is calculated as follows:
BI = Positive   Containers   Properties   Searched   × 100
Zhanjiang CDC staff followed the technical guidelines for dengue fever control issued by the Chinese Center for Disease Control and Prevention. Four streets/villages were selected in different geographical locations for each surveillance site; at least 100 residential households were surveyed, small water storage containers in indoor and outdoor plots were examined, and then mosquito larvae in positive containers were collected. To avoid the effect of continuous monitoring on mosquito density, two adjacent monitoring sessions were conducted in different households.

2.6. Transmission Model

The distribution of the DF cases in Zhanjiang City is shown in Figure 1. The darkest district was Chikan District.
In the ODE model based on our previous study [9,21], we found that the outbreak was locally transmitted by an imported case; therefore, we also included the index case in the model. We assumed that an imported case transmitted the virus to the outbreak source, where susceptible mosquitoes in Chikan District became infected with the dengue virus by biting the outbreak source and transmitted the virus by biting susceptible people in Chikan District; then, susceptible individuals may become asymptomatic or symptomatic through an incubation period or latent period. Next, these infected people may be bitten by susceptible mosquitoes in other districts, causing the disease to circulate in different districts. At last, the infected individuals would all recover with immunity. We divided individuals in Chikan District into the following five groups: Spc—susceptible, Epc—exposed, Ipc—infectious, Apc—asymptomatic and Rpc—removed. The subscript pc represents people in Chikan District. IIS represents the imported case that initiated the outbreak. Mosquitoes in Chikan District were divided into three groups: Smc—susceptible, Emc—exposed and Imc—infectious. The subscript mc represents mosquitoes in Chikan District. Nmc represents the sum of the mosquitoes in the three mosquito groups. Similarly, individuals in other districts were also divided into Spo, Epo, Ipo, Apo and Rpo, and mosquitoes were divided into Smo, Emo and Imo. The subscript po and mo represent people and mosquitoes in other districts, respectively. The definition model parameters are summarized in Table 1.
The mathematical model was described using the following ordinary differential equations:
d I I s d t = ε I I s d S m c d t = a c c c N m c n c I m c β I s S m c I I s β p c m c S m c A p c + I p c b c S m c d E m c d t = β I s S m c I I s + β p c m c S m c A p c + I p c b c + ω m E m c d I m c d t = a c c c n c I m c + ω m E m c b c I m c N m c = S m c + E m c + I m c d S p c d t = β m c p c S p c I m c β m c p c S p c I m o d E p c d t = β m c p c S p c I m c + β m c p c S p c I m o ω p E p c d I p c d t = 1 q ω p E p c γ I p c d A p c d t = q ω p E p c γ A p c d R p c d t = γ I p c + γ A p c d S m o d t = a o c o N m o n o I m o β p c m o S m o A p c + I p c β p o m o S m o A p o + I p o b o S m o d E m o d t = β p c m o S m o A p c + I p c + β p o m o S m o A p o + I p o b o + ω m E m o d I m o d t = a o c o n o I m o + ω m E m o b o I m o N m o = S m o + E m o + I m o d S p o d t = β m o p o S p o I m o d E p o d t = β m o p o S p o I m o ω p E p o d I p o d t = 1 q ω p E p o γ I p o d A p o d t = q ω p E p o γ A p o d R p o d t = γ I p o + γ A p o
Figure 2 shows the schematic of the transmission model.

2.7. Parameter Estimation

The values of the model parameters are summarized in Table 1. The season is one of the important conditions that affect the density of Aedes albopictus mosquitoes due to the periodicity of DF occurrence; therefore, we introduced the seasonal parameter c according to the previous study and used the trigonometric function to simulate the transmission of infection [22]. The incubation period of the dengue virus in humans usually spans from 4 to 8 days; 6 days is considered appropriate, and thus, we used ωp to represent the incubation relative rate of human infection and defined it as 0.1667 [23]. The extrinsic incubation period (EIP), i.e., the time interval between mosquito infection and when its bites become infectious, is 8–12 days; therefore, we considered 10 days in the simulation, i.e., we used ωm = 0.1000 [24] to represent the relative incubation rate of mosquito infection. The ratio of symptomatic to asymptomatic infection for DF is generally 2.2:1 [25]; therefore, the proportion of asymptomatic infections was q = 0.6875. The infection period is approximately 3–14 days [24]; therefore, the relative removal rate was γ = γ’ = ε = 0.1429. Due to the existing similar studies, we referred to the previously published literature to determine that ac = ao = bc = bo = 0.0714 [21]. The vertical infection rate of individual positive families (as a percentage of the vertical infection rate of offspring) ranges from 1.4% to 17.4% [26], which was simulated as 10.0% in the model; therefore, no = nc = 0.1000. The calculation of parameters cc and co were as follows, where τ and T refer to the simulation delay of the initial time in the whole season and the period of the season cycle, respectively:
c c = c o = c o s 2 π t τ T
On 13 September, the Zhanjiang Municipal Government launched a level III emergency response for this public health emergency, set up an emergency command, and established a joint prevention and control mechanism. On 28 September, Zhanjiang City held a dengue fever epidemic analysis meeting to discuss the situation at each epidemic site in Zhanjiang City and the measures taken, study and judge the development trend of the dengue fever epidemic in Zhanjiang City, discuss the implementation of prevention and control measures and the problems at each risk point, and propose a classification for the next step. The government increased its efforts to prevent and control the epidemic. Due to the interventions on 13 September and 28 September, we divided the dengue fever outbreak into three periods (17 August–12 September, 13–27 September and 28 September–28 October).

2.8. Scenarios

We further simulated two scenarios based on the fitted model. Scenario 1 assumed that the interventions on 13 September and 28 September did not occur and that dengue fever spread at the same rate in the second and third periods as in the first period. Scenario 2 assumed no intervention on 28 September and the rate of dengue fever transmission in the third period was the same as that in the 13 to 28 September period.

2.9. Simulation Method

Coefficients βIS, βpcmc, βmcpc, βpcmo, βmopc βpomo and βmopo were obtained by fitting the transmission models to the collected data with Berkeley Madonna 8.3.18 (developed by Robert Macey and George Oster of the University of California at Berkeley, CA, USA). The goodness of fit was assessed using the least root-mean-square error between the simulated and observed number of new cases per day. The coefficient of determination (R2), was assessed using SPSS 13.0 (IBM Corp, Armonk, NY, USA) to quantify the significance of the fit. In general, the fit was good when R2 is greater than 0.7. Make a map using Google Earth (Google Inc., Mountain View, CA, USA) software.

2.10. Normalization

In each model, we considered 7 transmission relative rates (βIS and so on); therefore, we fit each propagation coefficient separately for the case when the second and third cycles were zero, and then normalized the fitting results.
Constituent   ratio   ( y ) = x o t h e r x o t h e r + x C h i k a n
Normalized = y min ( y ) max ( y ) min ( y )
where Xother is the number of cases in other districts (Mazhang, Xiashan and Potou district), XChikan is the number of cases in Chikan district and y is the constituent ratio.
Based on the normalized results, we could compare the ability of different regions to spread the risk of disease and identify the districts with the highest risk of transmission. We defined that when the relative transmission rate was higher, the risk of disease transmission was higher.

3. Results

3.1. Reported Dengue Fever Cases in Zhanjiang City in 2018

In 2018, a total of 467 cases of DF were reported in Zhanjiang, including 447 indigenous cases and 20 imported cases, of which 15 and 5 were imported from domestic districts and overseas, respectively. In addition to the unrecorded serotypes, there were 222 cases with DENV 1, 7 cases with DENV 2 and 1 case with DENV 3. Overseas imports were from Cambodia (2 cases), Angola (1 case), Bangkok (1 case) and the Philippines (1 case). No deaths were reported. The annual incidence rate was 6.38 per 100,000 population. There were 234 clinically diagnosed and 233 laboratory-confirmed cases. The only DF outbreak appeared from late August to mid-October when cases spread from the construction sites to communities, and the outbreak was quickly contained after measures were taken to reduce the mosquito density.

3.2. Disease Distribution

On 22 April, the first case of imported dengue fever in Zhanjiang occurred in Chikan District, and on 19 July, the first case of local dengue fever occurred in Xiashan District, and on 17 August, the first case of local dengue fever occurred in Chikan District. The number of cases was sporadically distributed in many places. The outbreak ended on 28 October (Figure 3).
After 17 August, the dengue fever outbreak spread to all areas of Zhanjiang, except Xuwen County. There were three main hotspots with 333 cases in Chikan District, 40 cases in Xiashan District and 23 cases in Ma Zhang District. The distribution of cases in other areas was 22 cases in Leizhou City, 15 cases in Wuchuan City, 13 cases in Suixi County, 11 cases in Kaihua District, 5 cases in Potou District and 1 case in Lianjiang City. Xiashan District and Chikan District had the longest presence of dengue fever with 56 days (Table 2).
The 467 cases of dengue in 2018 involved multiple occupations. Workers (110 cases) accounted for the largest proportion (23.6%) of all cases, followed by household workers and the unemployed (80 cases, 17.1%), students (44 cases, 9.4%), the retired (41 cases, 8.8%) and farmers (41 cases, 8.8%) (Figure 4). The age distribution spanned from 6 months to 91 years old. Adults between the ages of 20 and 40 had the highest incidence (287 cases, 61%) (Figure 5).

3.3. Vector Surveillance

Many districts in Zhanjiang City were monitored from 22 September to 1 November. The data showed that the overall level of the BI in Zhanjiang City was over 10 from 22 September to 8 October. The BI dropped quickly from 9 October. When the BI decreased below 5, the number of new cases was almost zero. The last case in this region appeared on 28 October (Figure 6). Daily monitor points of the Breteau index are also shown in Figure 6. According to the Breteau index, we divided the risk levels and present them in the form of a map in Figure 7.

3.4. Curve Fitting

The results of the curve fitting of the reported data showed that the simulated results agreed well with the reported data: RChikan2 = 0.983 (p < 0.001), RKaifa2 = 0.791 (p < 0.001), RMazhang2 = 0.923 (p < 0.001) and RXiashan2 = 0.956 (p < 0.001). The values of βIS, βpcmc, βmcpc, βpcmo, βmopc βpomo and βmopo were generated via curve fitting and are summarized in Table 3.

3.5. Effectiveness of the Interventions

In scenario 1 (Figure 8), the simulation results showed that in Chikan District, the number of existing cases continued to increase, peaking on 7 October and then slowly decreasing. In Kaifa District, there was no increase in cases after 13 September. In Ma Zhang District, cases slowly increased after 13 September and reached a peak on 7 October. In Xiashan District, the cases increased slowly after 13 September and declined slowly after 7 October.
In scenario 2, existing cases in Chikan District peaked on 29 September and then slowly declined. In Kaifa district, cases continued to increase after 28 September and peaked on 14 October. In Ma Zhang, cases began to decline after 2 October. In Xiashan, cases continued to rise after 28 September and declined after peaking on 12 October.
The simulation results showed that the number of new cases in each model was at least twice the actual number without any intervention, yielding an epidemic duration [27] of 113 days. In scenario 2, the number of new cases in each model was at least 1.4 times the observed data, yielding a DO of 95 (96) days (Table 4).

3.6. Comparison of Transmission Relative Rate

The normalization results showed that the NβIS, Nβpcmc, Nβmopc and Nβpomo for the Chikan–Xiashan District transmission had the biggest values and the value of Nβpcmo for the Chikan–Kaifa District transmission was the smallest (Table 5).

4. Discussion

Guangdong Province has been a high prevalence area for DF. We described the dengue fever epidemic in Zhanjiang City, Guangdong Province, in 2018 and developed an ODE model to fit it to the observed data in Zhanjiang County. The R2 results showed that our model fit the observed data well. The normalized results showed that dengue fever spread in the region. Government interventions brought the outbreak under control.
Young adults aged 20 to 40 years accounted for the majority of cases, which was probably because the outbreak started at construction sites where most of the workers were young adults. Occupationally, workers accounted for the largest proportion, followed by domestic workers and the unemployed, suggesting that more emphasis should be placed on educating them about dengue-related health issues. In addition, workers often work around the clock without protective measures or while wearing light clothing during the hot summer months, which may contribute to mosquito-borne diseases [27]. In subtropical cities, the minimum temperature and humidity during a month are positively correlated with the number of mosquitoes [28,29]. This explains why in our study, the only peak of dengue incidence occurred from late August to mid-October. The appropriate climate and high density of Aedes mosquitoes provide the conditions for epidemics to occur. Socioeconomic factors also played an important role in the epidemic [30,31]. The sanitary conditions at construction sites are poor, with domestic waste and construction materials scattered everywhere. Workers are accustomed to living together, which can cause rapid transmission in the event of a case. Residents of Zhanjiang are accustomed to storing water in tanks, which provides conditions for Aedes mosquitoes to breed and may contribute to the spread of the outbreak from the construction site to the community.
BI surveillance was conducted throughout the epidemic period. The high vector density provided a high receptivity, and therefore, a high probability for DF transmission. Since DF is a mosquito-borne disease, killing mosquitoes can lead to a decrease in the disease burden. When control measures, such as mosquito eradication, were taken, the BI was kept below the safety line and the epidemic was controlled, suggesting that the timely reduction of the Aedes mosquitoes density can contribute to effective prevention and control of the epidemic, which is a result that is also supported by our model.
In this study, the SEIAR model was used to fit the epidemic curves of the outbreak in Zhanjiang City. The results of the coefficient of determination showed a good fit of our model to the reported data, which indicated that the model could fit a similar situation, i.e., comparing the transmissibility of DF between different regions. In scenario 1, the number of dengue cases in Zhanjiang City would have more than doubled and the duration of the outbreak would have more than doubled if the government had not implemented any interventions and instead allowed the disease to spread naturally. When the government took control measures on 13 September (scenario 2), there was some relief. However, the number of cases and the duration were also significantly higher than the actual situation. As we can see, the timely intervention by the government effectively mitigated the duration of the epidemic and the number of infected people. This shows that timely and multiple preventive and control measures are necessary in the event of a dengue fever outbreak.
Among the four districts considered, Chikan district had the strongest internal transmission, which indicated that the transmission of dengue fever was strongest within Chikan district. This was because Chikan District is located in the central area of Zhanjiang City, with large population mobility and a huge population. Karl et al. said human activity can promote the spread of dengue [32,33]. A study by Poter et al. [34] in two textile mills showed that some workers were infected at the workplace and some at home, suggesting that controlling human activity is an important condition for controlling dengue. One potential problem could be the weather conditions in different areas of Zhanjiang City. Different local weather also has an impact on the local mosquito activity profile in different parts of the city [35].
Zhanjiang has a permanent resident population of 7.33 million with high mobility. The epidemic first occurred in Chikan District, then in Mazhang, Kaifa and Xiashan Districts. This chronological order suggested that there may have been inter-regional transmission of the dengue virus. The normalized results showed that dengue fever was transmitted within the districts. This suggested that the detection of a case in one district can serve as an early warning signal and that other districts should implement appropriate interventions. When the transmission of the dengue virus from Chikan to other districts was cut off, the number of cases in other districts was significantly reduced or even zero. When measures were taken to control the epidemic diseases in Chikan District in time, such as controlling the population movement between regions, case isolation, adult vector control and larvae control, the epidemic in other districts was alleviated or even stopped from breaking out. In addition, in the analysis of inter-regional transmission, the greatest human-to-mosquito transmissibility was from Chikan District to Mazhang District, where infected people moved from Chikan District to Mazhang District and infected local mosquitoes, suggesting that the risk of DF transmission between these two districts was greater than in other districts. Once an outbreak occurs, the movement of people between these two districts must be more strictly controlled. Furthermore, we found that the greatest transmissibility from mosquitoes to people was from people who were infected with Xiashan District mosquitoes to people in Chikan District, suggesting that stricter prevention and intervention measures must be implemented in these two areas in terms of controlling the density of adults and larvae.

Limitations

Although we gained some insights from this research, we must recognize that the data used in this study were collected from the passive public health surveillance system. Reporting methods, under-reporting or misdiagnosis, and the availability of health facilities and laboratory diagnostics may have all affected the data quality. Moreover, some asymptomatic cases may not have visited the hospital, which was also not included in the surveillance system. Since there have been few dengue fever cases in Zhanjiang in the past few years, we used only one year of data to elucidate the pattern of transmission between districts. If data become available in the future, some more multidimensional and multi-year regional analyses can be done.

5. Conclusions

When a dengue outbreak occurs in one district, it spreads to other districts. Thus, a dengue outbreak in one district can serve as an early warning signal for other districts. The OED model can effectively model the spread of dengue between districts. Adopting the same level of control in all districts is not the best approach, and concentrating different control measures in different districts is the best way to control dengue outbreaks.

Author Contributions

M.Z.: Conceptualization, investigation, resources, writing—original draft; J.-F.H.: methodology, data curation, formal analysis, writing—original draft, software; M.K.: conceptualization, investigation, resources, writing—original draft; X.-C.L.: data curation, methodology, formal analysis, writing—original draft; H.-Y.L.: data curation, visualization; Z.-Y.Z.: methodology, formal analysis; G.-Q.Y.: investigation, resources; S.-N.L.: data curation; J.R.: data curation, visualization; J.-W.X.: data curation; Y.-Z.Z.: visualization; Y.W.: data curation, visualization; M.Y.: visualization; S.-X.T.: writing—review and editing, project administration, supervision, conceptualization, methodology; Q.C.: writing—review and editing, project administration, supervision, conceptualization; T.-M.C.: conceptualization, methodology, writing—review and editing, project administration, supervision. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data supporting the conclusions of this article are included within the article.

Acknowledgments

The authors thank all study participants for providing the data and field investigators for collecting the data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Liu, K.; Hou, X.; Wang, Y.; Sun, J.; Xiao, J.; Li, R.; Lu, L.; Xu, L.; Sang, S.; Hu, J.; et al. The driver of dengue fever incidence in two high-risk areas of China: A comparative study. Sci. Rep. 2019, 9, 19510. [Google Scholar] [CrossRef] [PubMed]
  2. Lai, S.; Huang, Z.; Zhou, H.; Anders, K.L.; Perkins, T.A.; Yin, W.; Li, Y.; Mu, D.; Chen, Q.; Zhang, Z.; et al. The changing epidemiology of dengue in China, 1990–2014: A descriptive analysis of 25 years of nationwide surveillance data. BMC Med. 2015, 13, 100. [Google Scholar] [CrossRef] [PubMed]
  3. Bhatt, S.; Gething, P.W.; Brady, O.J.; Messina, J.P.; Farlow, A.W.; Moyes, C.L.; Drake, J.M.; Brownstein, J.S.; Hoen, A.G.; Sankoh, O.; et al. The global distribution and burden of dengue. Nature 2013, 496, 504–507. [Google Scholar] [CrossRef] [PubMed]
  4. World Health Organization. Dengue and Severe Dengue: World Health Organization. 2020. Available online: https://www.who.int/health-topics/dengue-and-severe-dengue#tab=tab_1 (accessed on 1 December 2020).
  5. Kraemer, M.U.; Sinka, M.E.; Duda, K.A.; Mylne, A.Q.; Shearer, F.M.; Barker, C.M.; Moore, C.G.; Carvalho, R.G.; Coelho, G.E.; Bortel, W.V.; et al. The global distribution of the arbovirus vectors Aedes aegypti and Ae. albopictus. Elife 2015, 4, e08347. [Google Scholar] [CrossRef]
  6. Zhu, B.; Wang, L.; Wang, H.; Cao, Z.; Zha, L.; Li, Z.; Ye, Z.; Zhang, J.; Song, H.; Sun, Y. Prediction model for dengue fever based on interactive effects between multiple meteorological factors in Guangdong, China (2008–2016). PLoS ONE 2019, 14, e0225811. [Google Scholar] [CrossRef]
  7. Yang, L.; Chen, Y.; Yan, H.; Zhang, P.; Xu, X.; Tang, B.; Zhao, P.; Ren, R. A survey of the 2014 dengue fever epidemic in Guangzhou, China. Emerg. Microbes Infect. 2015, 4, e57. [Google Scholar] [CrossRef]
  8. Fan, J.; Lin, H.; Wang, C.; Bai, L.; Yang, S.; Chu, C.; Yang, W.; Liu, Q. Identifying the high-risk areas and associated meteorological factors of dengue transmission in Guangdong Province, China from 2005 to 2011. Epidemiol. Infect. 2014, 142, 634–643. [Google Scholar] [CrossRef]
  9. Liu, X.; Zhang, M.; Cheng, Q.; Zhang, Y.; Ye, G.; Huang, X.; Zhao, Z.; Rui, J.; Hu, Q.; Frutos, R.; et al. Dengue fever transmission between a construction site and its surrounding communities in China. Parasit. Vectors 2021, 14, 22. [Google Scholar] [CrossRef]
  10. Sang, S.; Chen, B.; Wu, H.; Yang, Z.; Di, B.; Wang, L.; Tao, X.; Liu, X.; Liu, Q. Dengue is still an imported disease in China: A case study in Guangzhou. Infect. Genet. Evol. 2015, 32, 178–190. [Google Scholar] [CrossRef]
  11. Zhao, S.; Musa, S.S.; Meng, J.; Qin, J.; He, D. The long-term changing dynamics of dengue infectivity in Guangdong, China, from 2008–2018: A modelling analysis. Trans. R. Soc. Trop. Med. Hyg. 2020, 114, 62–71. [Google Scholar] [CrossRef]
  12. Wu, J.Y.; Lun, Z.R.; James, A.A.; Chen, X.G. Dengue Fever in mainland China. Am. J. Trop. Med. Hyg. 2010, 83, 664–671. [Google Scholar] [CrossRef] [PubMed]
  13. Huang, L.; Luo, X.; Shao, J.; Yan, H.; Qiu, Y.; Ke, P.; Zheng, W.; Xu, B.; Li, W.; Sun, D.; et al. Epidemiology and characteristics of the dengue outbreak in Guangdong, Southern China, in 2014. Eur. J. Clin. Microbiol. Infect. Dis. 2016, 35, 269–277. [Google Scholar] [CrossRef] [PubMed]
  14. Zou, L.; Chen, J.; Feng, X.M.; Ruan, S.G. Analysis of a Dengue Model with Vertical Transmission and Application to the 2014 Dengue Outbreak in Guangdong Province, China. Bull. Math. Biol. 2018, 80, 2633–2651. [Google Scholar] [CrossRef]
  15. Yang, H.M.; Boldrini, J.L.; Fassoni, A.C.; Freitas, L.F.; Gomez, M.C.; de Lima, K.K.; Andrade, V.R.; Freitas, A.R.R. Fitting the Incidence Data from the City of Campinas, Brazil, Based on Dengue Transmission Modellings Considering Time-Dependent Entomological Parameters. PLoS ONE 2016, 11, e0152186. [Google Scholar] [CrossRef] [PubMed]
  16. Chavez, J.P.; Gotz, T.; Siegmund, S.; Wijaya, K.P. An SIR-Dengue transmission model with seasonal effects and impulsive control. Math. Biosci. 2017, 289, 29–39. [Google Scholar] [CrossRef]
  17. Bhuju, G.; Phaijoo, G.R.; Gurung, D.B. Fuzzy Approach Analyzing SEIR-SEI Dengue Dynamics. Biomed. Res. Int. 2020, 2020, 1508613. [Google Scholar] [CrossRef] [PubMed]
  18. Grunnill, M. An exploration of the role of asymptomatic infections in the epidemiology of dengue viruses through susceptible, asymptomatic, infected and recovered (SAIR) models. J. Theor. Biol. 2018, 439, 195–204. [Google Scholar] [CrossRef]
  19. Xu, Z.; Bambrick, H.; Pongsumpun, P.; Ming Tang, I.; Yakob, L.; Devine, G.; Frentiu, F.D.; Williams, G.; Hu, W. Does Bangkok have a central role in the dengue dynamics of Thailand? Parasit. Vectors 2020, 13, 22. [Google Scholar] [CrossRef]
  20. Sun, J.; Zhang, H.; Tan, Q.; Zhou, H.; Guan, D.; Zhang, X.; Duan, J.; Cai, S.; Peng, Z.; He, J.; et al. The epidemiological characteristics and molecular phylogeny of the dengue virus in Guangdong, China, 2015. Sci. Rep. 2018, 8, 9976. [Google Scholar] [CrossRef]
  21. Yi, B.; Chen, Y.; Ma, X.; Rui, J.; Cui, J.A.; Wang, H.; Li, J.; Chan, S.F.; Wang, R.; Ding, K.; et al. Incidence dynamics and investigation of key interventions in a dengue outbreak in Ningbo City, China. PLoS Negl. Trop. Dis. 2019, 13, e0007659. [Google Scholar] [CrossRef] [Green Version]
  22. Yang, H.M.; Macoris, M.L.; Galvani, K.C.; Andrighetti, M.T.; Wanderley, D.M. Assessing the effects of temperature on the population of Aedes aegypti, the vector of dengue. Epidemiol. Infect. 2009, 137, 1188–1202. [Google Scholar] [CrossRef] [PubMed]
  23. Chan, M.; Johansson, M.A. The incubation periods of Dengue viruses. PLoS ONE 2012, 7, e50972. [Google Scholar] [CrossRef] [PubMed]
  24. Andraud, M.; Hens, N.; Marais, C.; Beutels, P. Dynamic epidemiological models for dengue transmission: A systematic review of structural approaches. PLoS ONE 2012, 7, e49085. [Google Scholar] [CrossRef]
  25. Wang, T.; Wang, M.; Shu, B.; Chen, X.Q.; Luo, L.; Wang, J.Y.; Cen, Y.Z.; Anderson, B.D.; Merrill, M.M.; Merrill, H.R.; et al. Evaluation of inapparent dengue infections during an outbreak in Southern China. PLoS Negl. Trop. Dis. 2015, 9, e0003677. [Google Scholar] [CrossRef] [PubMed]
  26. Grunnill, M.; Boots, M. How Important is Vertical Transmission of Dengue Viruses by Mosquitoes (Diptera: Culicidae)? J. Med. Entomol. 2016, 53, 1–19. [Google Scholar] [CrossRef]
  27. Sang, S.; Yin, W.; Bi, P.; Zhang, H.; Wang, C.; Liu, X.; Chen, B.; Yang, W.; Liu, Q. Predicting local dengue transmission in Guangzhou, China, through the influence of imported cases, mosquito density and climate variability. PLoS ONE 2014, 9, e102755. [Google Scholar] [CrossRef]
  28. Liang, L.; Hualiang, L.; Linwei, T.; Weizhong, Y.; Jimin, S.; Qiyong, L. Time series analysis of dengue fever and weather in Guangzhou, China. BMC Public Health 2009, 9, 395. [Google Scholar] [CrossRef]
  29. Zonetti, L.F.C.; Coutinho, M.C.; de Araujo, A.S. Molecular Aspects of the Dengue Virus Infection Process: A Review. Protein Pept. Lett. 2018, 25, 712–719. [Google Scholar] [CrossRef]
  30. Gubler, D.J. Dengue, Urbanization and Globalization: The Unholy Trinity of the 21(st) Century. Trop. Med. Health 2011, 39 (Suppl. S4), 3–11. [Google Scholar] [CrossRef]
  31. Hales, S.; de Wet, N.; Maindonald, J.; Woodward, A. Potential effect of population and climate changes on global distribution of dengue fever: An empirical model. Lancet 2002, 360, 830–834. [Google Scholar] [CrossRef]
  32. Stoddard, S.T.; Forshey, B.M.; Morrison, A.C.; Paz-Soldan, V.A.; Vazquez-Prokopec, G.M.; Astete, H.; Reiner, R.C., Jr.; Vilcarromero, S.; Elder, J.P.; Halsey, E.S.; et al. House-to-house human movement drives dengue virus transmission. Proc. Natl. Acad. Sci. USA 2013, 110, 994–999. [Google Scholar] [CrossRef] [PubMed]
  33. Karl, S.; Halder, N.; Kelso, J.K.; Ritchie, S.A.; Milne, G.J. A spatial simulation model for dengue virus infection in urban areas. BMC Infect. Dis. 2014, 14, 447. [Google Scholar] [CrossRef] [PubMed]
  34. Nevai, A.L.; Soewono, E. A model for the spatial transmission of dengue with daily movement between villages and a city. Math. Med. Biol. 2014, 31, 150–178. [Google Scholar] [CrossRef] [PubMed]
  35. Mincham, G.; Baldock, K.L.; Rozilawati, H.; Williams, C.R. Development of a mechanistic dengue simulation model for Guangzhou. Epidemiol. Infect. 2019, 147, e125. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Cumulative number of cases and cumulative incidence of dengue fever in Zhanjiang City in 2018.
Figure 1. Cumulative number of cases and cumulative incidence of dengue fever in Zhanjiang City in 2018.
Tropicalmed 07 00209 g001
Figure 2. Schematic of the dengue fever transmission model.
Figure 2. Schematic of the dengue fever transmission model.
Tropicalmed 07 00209 g002
Figure 3. The outbreak of Dengue fever in Zhanjiang in 2018.
Figure 3. The outbreak of Dengue fever in Zhanjiang in 2018.
Tropicalmed 07 00209 g003
Figure 4. Profession distribution of dengue fever in Zhanjiang, Guangdong Province, China, in 2018 (HAU*: Housework an unemployed; KC*: Kindergarten children).
Figure 4. Profession distribution of dengue fever in Zhanjiang, Guangdong Province, China, in 2018 (HAU*: Housework an unemployed; KC*: Kindergarten children).
Tropicalmed 07 00209 g004
Figure 5. Gender distribution of dengue fever in Zhanjiang, Guangdong Province, China, in 2018. There were 272 males and 195 females, with a sex ratio of 1.39:1 and there was no significant difference in the sex ratio of cases (p > 0.05).
Figure 5. Gender distribution of dengue fever in Zhanjiang, Guangdong Province, China, in 2018. There were 272 males and 195 females, with a sex ratio of 1.39:1 and there was no significant difference in the sex ratio of cases (p > 0.05).
Tropicalmed 07 00209 g005
Figure 6. The number of new cases and the Breteau index in Zhanjiang in 2018.
Figure 6. The number of new cases and the Breteau index in Zhanjiang in 2018.
Tropicalmed 07 00209 g006
Figure 7. Distribution map of the Breteau index for dengue fever in Zhanjiang City in 2018.
Figure 7. Distribution map of the Breteau index for dengue fever in Zhanjiang City in 2018.
Tropicalmed 07 00209 g007
Figure 8. Curve fitting of dengue fever in Zhanjiang City, 2018.
Figure 8. Curve fitting of dengue fever in Zhanjiang City, 2018.
Tropicalmed 07 00209 g008
Table 1. Parameter definitions and values.
Table 1. Parameter definitions and values.
ParameterDescriptionUnitValueRangeMethod
ac/aoDaily birth rate of mosquitoes in Chikan district or other districts *1day−10.07140.0200–0.2500References
cc/coSeasonality parameter of the mosquitoes’ population in Chikan district or other districts1See Equation (2) *20–1Curve fitting
nc/noProportion of transovarial transmission in Chikan District or other districts10.10.0140–0.1740Reference
bc/boDaily death rate of mosquitoes in Chikan District or other districtsday−10.07140.0200–0.2500Reference
ωm/ωpRelative incubation rate of mosquito infection or human infectionday−10.1000/0.16670.0833–0.1250/Reference
0.1250–0.2500
qProportion of asymptomatic infections10.68750–1Reference
γRelative removal rate of infectious individualsday−10.14290.0714–0.3333Reference
γ’Relative removal rate of asymptomatic individualsday−10.14290.0714–0.3333Reference
τSimulation delay of the initial time in the whole seasonday279≥0Analysis of the reported data
TDuration of the cycleday365≥0Analysis of the reported data
βISRelative transmission rate from imported mosquitoes to mosquitoes in Chikan District1See Table 2≥0Curve fitting
βmcpcRelative transmission rate from mosquitoes in Chikan District to humans in Chikan District1See Table 2≥0Curve fitting
βmopcRelative transmission rate from mosquitoes in other districts to humans in Chikan District1See Table 2≥0Curve fitting
βmopoRelative transmission rate from mosquitoes in other districts to humans in other districts1See Table 2≥0Curve fitting
βpcmcRelative transmission rate from humans in Chikan District to mosquitoes in Chikan District1See Table 2≥0Curve fitting
βpcmoRelative transmission rate from humans in Chikan District to mosquitoes in other districts1See Table 2≥0Curve fitting
βpomoRelative transmission rate from humans in other districts to humans in other districts1See Table 2≥0Curve fitting
*1 Subscript c indicates Chikan District, subscript o indicates other districts (Kaifa, Xiasha and Mazhang District). *2 The parameter c was simulated with the trigonometric function using the parameters τ and T.
Table 2. Regional distribution of dengue fever outbreak in Zhanjiang.
Table 2. Regional distribution of dengue fever outbreak in Zhanjiang.
VariableMaleFemaleTotalDuration (Days)
District
Chikan20013333356
Xiashan22184056
Leizhou10122225
Mzhang13102342
Wuchuan691545
Suixi581334
Kaifa831123
Potou41550
Lianjiang0111
Total268195463-
Table 3. Parameter values.
Table 3. Parameter values.
ParameterChikan–KaifaChikan–MazhangChikan–Xiashan
βIS1.73 × 10−92.23 × 10−91.53 × 10−9
βpcmc3.19 × 10−84.13 × 10−83.10 × 10−8
βmcpc1.86 × 10−51.49 × 10−51.94 × 10−5
βmopc7.88 × 10−101.18 × 10−97.31 × 10−10
βpomo5.16 × 10−166.72 × 10−165.30 × 10−16
βpcmo5.21 × 10−82.55 × 10−82.57 × 10−8
βmopo5.67 × 10−71.66 × 10−61.70 × 10−6
Table 4. Scenario fitting results.
Table 4. Scenario fitting results.
DistrictFitScenario 1Scenario 2
Number of New CasesDurationNumber of New CasesProportion *DurationNumber of New CasesProportion *Duration
Chikan41056974137.56%11458141.71%95
Kaifa1223---43258.33%37
Chikan40656932129.56%11457942.61%95
Xiashan3455---97185.29%54
Chikan40656917125.86%11358644.33%96
Mazhang2741---4359.26%33
*: The increase ratio of the number of new cases in the two scenarios to the fitted data.
Table 5. Normalized results.
Table 5. Normalized results.
Relative Transmission Chikan–Kaifa DistrictChikan–Xiashan DistrictChikan–Mazhang District
NβIS0.4414410.4632940.389146
Nβpcmc0.4491670.4644210.39408
Nβmcpc111
Nβmopc0.4414410.4632940.389146
Nβpomo0.4414410.4632940.389146
Nβpcmo0.0640070.2113370.389146
Nβpcmo000
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, M.; Huang, J.-F.; Kang, M.; Liu, X.-C.; Lin, H.-Y.; Zhao, Z.-Y.; Ye, G.-Q.; Lin, S.-N.; Rui, J.; Xu, J.-W.; et al. Epidemiological Characteristics and the Dynamic Transmission Model of Dengue Fever in Zhanjiang City, Guangdong Province in 2018. Trop. Med. Infect. Dis. 2022, 7, 209. https://0-doi-org.brum.beds.ac.uk/10.3390/tropicalmed7090209

AMA Style

Zhang M, Huang J-F, Kang M, Liu X-C, Lin H-Y, Zhao Z-Y, Ye G-Q, Lin S-N, Rui J, Xu J-W, et al. Epidemiological Characteristics and the Dynamic Transmission Model of Dengue Fever in Zhanjiang City, Guangdong Province in 2018. Tropical Medicine and Infectious Disease. 2022; 7(9):209. https://0-doi-org.brum.beds.ac.uk/10.3390/tropicalmed7090209

Chicago/Turabian Style

Zhang, Meng, Jie-Feng Huang, Min Kang, Xing-Chun Liu, Hong-Yan Lin, Ze-Yu Zhao, Guo-Qiang Ye, Sheng-Nan Lin, Jia Rui, Jing-Wen Xu, and et al. 2022. "Epidemiological Characteristics and the Dynamic Transmission Model of Dengue Fever in Zhanjiang City, Guangdong Province in 2018" Tropical Medicine and Infectious Disease 7, no. 9: 209. https://0-doi-org.brum.beds.ac.uk/10.3390/tropicalmed7090209

Article Metrics

Back to TopTop