Next Article in Journal
Lupeol Acetate and α-Amyrin Terpenes Activity against Trypanosoma cruzi: Insights into Toxicity and Potential Mechanisms of Action
Next Article in Special Issue
Integrated Strategies for Aedes aegypti Control Applied to Individual Houses: An Approach to Mitigate Vectorial Arbovirus Transmission
Previous Article in Journal
It Was Not the Perfect Storm: The Social History of the HIV-2 Virus in Guinea-Bissau
Previous Article in Special Issue
Dengue as a Disease Threatening Global Health: A Narrative Review Focusing on Latin America and Brazil
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatial Analysis of Dengue Clusters at Department, Municipality and Local Scales in the Southwest of Colombia, 2014–2019

by
Catalina Marceló-Díaz
1,*,
María Camila Lesmes
1,2,
Erika Santamaría
1,
José Alejandro Salamanca
2,
Patricia Fuya
1,
Horacio Cadena
3,
Paola Muñoz-Laiton
1 and
Carlos Andrés Morales
4
1
Grupo de Entomología, Instituto Nacional de Salud, Bogotá 111321, Colombia
2
Facultad de Ciencias Ambientales y de la Sostenibilidad, Programa de Ingeniería Geográfica y Ambiental, Universidad de Ciencias Aplicadas y Ambientales, UDCA, Bogotá 111166, Colombia
3
Programa de Estudio y Control de Enfermedades Tropicales, PECET, Universidad de Antioquia, Medellín 050010, Colombia
4
Secretaría de Salud Departamental del Cauca, Popayán 190001, Colombia
*
Author to whom correspondence should be addressed.
Trop. Med. Infect. Dis. 2023, 8(5), 262; https://doi.org/10.3390/tropicalmed8050262
Submission received: 7 February 2023 / Revised: 7 April 2023 / Accepted: 26 April 2023 / Published: 2 May 2023
(This article belongs to the Special Issue Emerging Topics in Arbovirus Vectors)

Abstract

:
Dengue is an arbovirus transmitted by mosquitoes of the genus Aedes and is one of the 15 main public health problems in the world, including Colombia. Where limited financial resources create a problem for management, there is a need for the department to prioritize target areas for public health implementation. This study focuses on a spatio-temporal analysis to determine the targeted area to manage the public health problems related to dengue cases. To this end, three phases at three different scales were carried out. First, for the departmental scale, four risk clusters were identified in Cauca (RR ≥ 1.49) using the Poisson model, and three clusters were identified through Getis-Ord Gi* hotspots analysis; among them, Patía municipality presented significantly high incidence rates in the time window (2014–2018). Second, on the municipality scale, altitude and minimum temperature were observed to be more relevant than precipitation; considering posterior means, no spatial autocorrelation for the Markov Chain Monte Carlo was found (Moran test ˂ 1.0), and convergence was reached for b1–b105 with 20,000 iterations. Finally, on the local scale, a clustered pattern was observed for dengue cases distribution (nearest neighbour index, NNI = 0.202819) and the accumulated number of pupae (G = 0.70007). Two neighbourhoods showed higher concentrations of both epidemiological and entomological hotspots. In conclusion, the municipality of Patía is in an operational scenario of a high transmission of dengue.

1. Introduction

Dengue is the most rapidly expanding arboviruses in the world, and its incidence has increased in recent decades in tropical and subtropical countries, becoming a global public health problem [1,2]. In the Americas there was recorded the highest number of dengue cases in history in 2019 with more than 2.7 million cases, including 22,127 serious cases and 1206 deaths [3].
In Colombia, approximately 58% of the population is at risk of having the disease, since it is endemic to approximately 909 municipalities, and 108 municipalities are centres of dengue transmission [1]. In 2021, during the epidemiological week 52, about 53,334 cases were recorded, of which 50% corresponded to dengue with warning signs, and 1.8% corresponded to severe dengue. Although the behaviour of the disease at the national level was within the expected range, the resulting mortality showed a significant increase. Some departments, including Cauca, are in a situation of alert [4], and Cauca itself reported 336 cases over the course of 2022. Between 2020 and 2021, this department was also under alert since its epidemiological behaviour exceeded the expected threshold [4,5].
The main vector of dengue in Colombia is the Ae. aegypti mosquito, an anthropophilic species that is found mainly in urban areas. In general, this vector infests natural or artificial containers in or around homes [6]. Its distribution has been reported in all departments of Colombia up to 2302 m.a.s.l. [7]. Under optimal conditions of food availability and adequate oviposition sites, the average dispersal of a female Aedes spp. mosquito is estimated to be between 50 and 100 m, which limits its visits to two or three dwellings during its adult life [8], although fed females have been recorded to disperse as far as 800 m in 6 days [9].
The dynamics that allow the transmission of dengue respond to a series of general and specific factors known as macrodeterminants and microdeterminants of dengue. Macrodeterminants can be (1) environmental, such as altitude, latitude and climatic conditions that facilitate larval development, vector survival and virus replication in the vector [10] and (2) demographic and socioeconomic conditions, such as intraurban mobility, inadequate urbanization patterns and lack of access to water supply and sewerage services, as well as cultural conditions, such as a low perception of risk. Microdeterminants include host characteristics, the circulation of the various dengue serotypes and entomological factors. These considerations highlight the presence of various conditions that favour the endemic–epidemic transmission of dengue in the different municipalities of the department of Cauca [10].
These variables provide valuable information, which is collected by health institutions and thus forms the early warning components of systematic public health surveillance systems, when this surveillance is based on indicators. Although dengue surveillance is traditionally syndromic or sentinel, there are also alternative sources of information which do not necessarily belong to the health field and could be taken into account to carry out surveillance based on indicators such as meteorological, entomological and environmental data.
Given the multiplicity of factors that affect the transmission of dengue, geographic information systems (GISs) have become a very useful tool for identifying the distribution patterns of events and studying their association with the macrodeterminants mentioned above.
Thus, the use of global and local spatial-temporal analysis tools such as the nearest neighbour index (NNI) [11]; the Getis Ord-Gi* statistic, used to detect hot and cold dengue spots with geographically homogeneous high or low values [12]; and kernel density, applied in several studies on dengue [13,14] among other techniques, allow the generation of maps that highlight the geographic areas and population groups at risk [15,16]. Those generated maps can serve as a basis for responsible institutions to design strategies for preventing and reducing the risk of dengue epidemics.
In order to assist in determining the problematic dengue areas to ensure an effective management of the dengue cases, a spatio-temporal analysis was conducted. This analysis contributes to the process of spatial stratification, allowing the prioritization of control actions, to generate base information for developing contingency plans for vector-borne diseases, such as dengue fever.
The primary objective of this study was to analyse the spatio-temporal distribution of dengue in the department of Cauca to identify high-dengue-risk municipalities through spatial scan statistics and cluster analysis.
The department of Cauca was selected due to its epidemiological behaviour in the last ten years, with the most important epidemic outbreak recorded in this period occurring between the end of 2009 and 2010. The endemic–epidemic characteristics of Cauca provide a suitable study area to carry out a series of spatial analyses since it presents temporal, seasonal, and cyclical behaviours in at-risk populations, enabling a better characterization without presenting biases due to epidemics or hyperendemic trends.
A secondary aim was to contribute to the spatial stratification of dengue in the identified high-risk municipalities, which includes determining the effect of environmental risk factors on the rate of dengue cases and the identification of important transmission foci [8]. With this study, the finding determines the level of risks and helps to determine actions needed to control the risks.
A third objective was to prioritize vector-control activities through the identification and analysis of persistent dengue and entomological clusters by neighbourhood, to reveal correlations and to generate base information for the elaboration of contingency plans for vector-borne diseases, such as dengue, in problem areas which will be targeted for intensive intervention.

2. Materials and Methods

The study comprised three phases. The first phase targeted the locations of risk clusters in the department of Cauca through a spatio-temporal analysis. The department of Cauca was selected because it has been on epidemiological warning for more than three consecutive years since 2019, given the higher-than-expected increase in dengue cases, as reported in the epidemiological reports of the National Institute of Health (INS, acronym in Spanish) [4,5]. In the second phase, environmental variables related to the disease were identified through the development of a Poisson regression model for the municipality of Patía, which was selected because it shows a relatively moderate dengue risk. In the third phase of the study, areas with higher-than-expected incidence of the disease were established and clusters were identified at the neighbourhood scale (an urban spatial unit), using spatial analysis processes in GIS software. These urban spatial units should be targeted for intensive intervention and surveillance. Figure 1 shows the development of each of the phases.

2.1. Study Area

The department of Cauca is located in southwestern Colombia between the coordinates 00°58′ and 3°19′ north latitude and 75°47′ and 77°57′ west longitude. It is divided into 42 municipalities covering an area of 29,308 km2 and has a population of 1,491,937 inhabitants [17]. Its average altitude is 1693 m.a.s.l.; however, it contains all thermal floors, which is reflected in its diversity of ecosystems. It contains mangrove forests; moors; wetlands; warm forests; sub-Andean, Andean and high-Andean forests; inter-Andean valleys; and dry and subxerophytic forests [18].
The climate in the department of Cauca has a great variability depending on its geographic and geomorphologic characteristics. On the Pacific coast of the department, there are warm wet and semi-humid pluvial climates, with rainfall that can range from 3000 mm to 7000 mm per year and average temperatures of 24 °C. In the central subregion, the climate is predominantly temperate wet, with temperatures that oscillate between 18 °C and 24 °C. The Patía Valley stands out as the driest area in the department as it has a warm dry climate and receives less than 1500 mm of rainfall per year. In contrast, the eastern edge of Cauca has a cold dry and very cold dry climate, while the Amazonian foothills has a warm wet climate.
Cauca has seven subregions: the central subregion, where the capital Popayán is located; the northern subregion, divided into 13 municipalities, including Miranda, Puerto Tejada and Villa Rica; the eastern subregion; the Pacific subregion; the southern subregion, where the municipality of Patía (see Section 2.3.2) and a large part of the Patía River valley are located; the subregion of the Massif; and the Amazonian Piedmont subregion, which includes the municipality of Piamonte [18].
For each phase of the study, we considered the results obtained in the previous phase. For Phase 1, Cauca was chosen as we explained before. For Phase 2, Patía municipality was selected considering the findings of spatio-temporal analysis and for Phase 3 we chose its municipal seat, because of its important as an administrative centre, seat of government, or capital city. Additionally, this municipal seat has been positive for pupae sampling throughout several years in several neighbourhoods.

2.2. Data Description

Epidemiological data regarding dengue and severe dengue cases were used in each of the phases of the study, including laboratory-confirmed cases and cases confirmed by epidemiological linkage. Data also included sociodemographic characteristics at the patient level such as age, sex, occupation, and ethnic data. These data are routinely reported by the Unidades Primarias Generadoras de Datos (UPGD in spanish) following the guidelines established by the WHO in 2009 [19].
In Colombia, a laboratory-confirmed case is a probable case of dengue, severe dengue, or dengue mortality confirmed by any of the laboratory criteria for dengue diagnosis. A positive PCR or viral isolation in patients with less than 5 days of onset of fever or Dengue IgM Dengue ELISA in patients with 5 or more days of onset of fever, while a confirmed case by epidemiological linkage involves confirming probable cases of dengue from laboratory-confirmed cases using the association of person, time and space [19]. This information is recorded in a database that contains epidemiological records by epidemiological week, so for the present study it was necessary to group them by year.
The procedures and specifications applied to the data for the development of the phases of the study are described below:

2.2.1. Dengue Incidence Data for the Department Scale

For the spatio-temporal analysis, dengue event notification data from the department of Cauca recorded in the period 2012 to 2018, acquired through the National System of Public Health Surveillance (SIVIGILA, for its Spanish abbreviation) were used. Data available at https://www.ins.gov.co/Direcciones/Vigilancia/Paginas/SIVIGILA.aspx (accessed on 10 January 2020).
The year 2019 was excluded from the analysis because it was an epidemic period; therefore, to avoid overestimating the cluster, that year was not considered in the spatio-temporal analysis. For hot spots detection, the incidence rates of dengue per 10,000 inhabitants from 2014 to 2018 were calculated from the cases reported to the SIVIGILA, and the population statistics of the projection system reported by the National Administrative Department of Statistics in 2020; in both process data were georeferenced. In 2014, a total of 378 cases of dengue fever were registered in the department of Cauca, with an increase of 40% in 2015 with 631 cases. About 803 cases were confirmed in 2016 and only 96 cases in 2017. Finally, in 2018, there were 177 cases registered.

2.2.2. Dengue Data for the Local Scale

For the municipality of Patía the databases of dengue and severe dengue cases from 2015 to 2019 reported by the Epidemiology Surveillance Group of the Department of Health of Cauca were reviewed (data available by request). Of the 246 total cases, 66 were discarded for reasons such as duplication of information, addresses that could not be found or different origins. As a result, 180 cases of dengue were used, including six cases of severe dengue, of whom 167 were georeferenced and 13 were geocode. To explore the data, a dengue endemic channel for the municipality of Patía was developed. An endemic channel represents the number of cases within the expected normal seasonal range; anything above this moving threshold would be considered representative of an unprecedented number of cases, i.e., an outbreak (Figure 2). The zones of success, security, and alarm determine the habitual behavior of the dengue cases.
Events were differentiated when they had been confirmed by a laboratory or clinic, confirmed by an epidemiological linkage, discarded, updated, or adjusted due to typing error. Subsequently, only the variables that were necessary to carry out the geocoding or georeferencing process were selected as appropriate. Whether the case had moved from his or her municipality of residence was also noted, and the location of the address and/or the neighbourhood was taken into account.
The geocoding process, which consists of assigning geographic coordinates to the addresses, was carried out in cases where the detailed description of the location was available, through Google Maps Geocoding API [20]. Georeferencing, which also involves the positioning of the element in a coordinate system, was applied for general descriptions referring to polygon-type elements.
Georeferencing was carried out using the point-radius method [21]. The assignment of coordinates was performed using the centroid of the most specific political divisions, which could correspond to a village, section, or neighbourhood. This centroid was calculated using the Calculate Geometry tool of ArcGIS® 10.8 for regular polygons and the Point on Surface tool of QGIS 3.10.5 for irregular polygons, using the official layers of the rural and urban sections provided by DANE [22]. For the polygons of neighbourhoods in the municipality, the layer was downloaded from the collaborative project Open Street Maps (OSM) due to its greater precision (https://www.openstreetmap.org/#map=5/4.632/-74.299 accessed on 25 June 2020). Data with ambiguous locations and data with neighbourhood or village locations that did not correspond to the information of the municipality or cases that reported a municipality of residence other than Patía were excluded.

2.3. Study Design

A retrospective study based on the epidemiological data collected for the observation window from 2012 to 2019; after performing its georeferencing, the study was developed at three scales: departmental, municipal, and local.

2.3.1. Phase 1: Spatio-Temporal Analysis of Dengue Distribution for Cluster Detection, Targeting High-Dengue-Risk Clusters of Municipalities at the Departmental Scale (Objective 1)

  • Space-time clustering: A space-time scan statistic for the detection of risk clusters following the Poisson model [23] was used for the years 2012 to 2018 when a spatial pattern of dengue cases existed. Targeting was carried out using the reported epidemiological information and included the following steps: (1) grouping of the epidemiological cases of the event reported to SIVIGILA according to the municipality of residence of the case or where the case originated, (2) information regarding the total number of cases of dengue and severe dengue taking into account the time of year and (3) location of the centroids in terms of latitude and longitude for each municipality in the department.
  • The information was processed using Kulldorff’s method [23]. This method uses a model in which the number of events in a geographic area is distributed according to the Poisson model; under the null hypothesis, and without covariates, the expected number of cases in each area is proportional to its population size, or to the person-years in that area. Because the data are aggregated in census districts, the measurement was concentrated in terms of the central coordinates of those districts and was expanded along a third dimension that reflects the size of the population as it changes over time.
  • The retrospective space-time analysis scanning for clusters with high rates using the discrete Poisson model was defined by a cylindrical window with a circular geographic base and with a height corresponding to time. The base was defined exactly as for the purely spatial scan statistic, while the height reflects the time frame of potential clusters. The cylindrical window was then moved in space and time, so that for each possible geographical location and size, it also visited each possible time frame. In effect, we obtained an infinite number of overlapping cylinders of different size and shape, jointly covering the entire study region, where each cylinder reflected a possible cluster [23].
  • Hotspot detection: Once the overall landscape of the Cauca department was defined by clustering risk, a more detailed analysis was carried out to distinguish the year-to-year epidemiological behaviour of each municipality using the High-Low Clustering technique (Getis-Ord General and Getis Ord-Gi* statistic). This technique identifies the areas where the highest incidence rates of dengue are geographically homogeneous and concentrated under statistical significance parameters. These areas, known as hot and cold spots, as they specifically consider the extreme high or low values, have been previously used by authors such as Khormi and Kumar [24] and Mutheneni et al. [12] to study the spatial patterns of dengue and their potential application in disease management programmes.
  • The clusters that presented a pattern of distribution of the incidence rates of dengue per 10,000 inhabitants—that is, data that were not randomly distributed in the department of Cauca for the period between 2014 and 2018—were identified. The conceptualization of inverse distance was determined using the Euclidean distance without defining a distance threshold. Additionally, the statistical significance was based on the False Discovery Rate (FDR) correction.
  • These spatial geostatistical techniques allowed us to characterize the geographic distribution of the incidence of dengue [25] and to thereby identify statistically significant hot spots and cold spots using the Getis-Ord Gi* statistic. This statistic allowed us to examine the local level of spatial clustering to identify and visualize the municipalities whose dengue rate values were extreme and geographically homogeneous.

2.3.2. Phase 2: Spatial Stratification of Dengue Risk Villages by Poisson Regression, to Determine the Effect of Environmental Factors That Influence the Spatial Variation of Dengue at the Municipal Scale (Objective 2)

  • Study site: Patía municipality was chosen for spatial and environmental factors analysis within the department after carrying out Phase 1; the municipality had one of the highest incidences of dengue fever in the time window (2012–2018). In relation to entomological indicators, it has been positive for pupae index throughout the years in several neighbourhoods of its municipal seat.
  • The municipality is located in the department of Cauca and has an area of 723 km2. As of 2021, it had a population of 37,793 inhabitants, of whom 13,598 lived in the urban area. The municipal seat, called Bordo-Patía, is located at the coordinates 02°06′56″ N and 76°59′21″ W and at an altitude of 910 m.a.s.l. The average temperature of the municipality is between 25 and 27 °C (maximum temperature 33–38 °C and minimum temperature 15–19 °C), and the average annual precipitation is 2171 mm. The municipality includes zones of premontane rainforest, premontane dry forest, and tropical dry forest [26].
  • The climate in Patía is conditioned by its geoforms. A warm climate is representative of the depression or valley of Patía, which covers most of its extent. A cold climate is associated with the western mountain range that passes through the north of the municipality. Precipitation in the municipality has a bimodal distribution and is divided into two wet periods (March–May and October–December, the latter of which is more intense) separated by two dry periods (January–February and June–September) [27].
  • Environmental factors: The environmental variables that were considered in the Poisson regression were altitude, minimum and maximum temperature, and precipitation. To acquire the data associated with these variables, the NASA Earth Data platform and the Climate Engine climate database were used.
  • Altitude data were obtained from the NASA website, https://search.asf.alaska.edu/#/ (accessed on 8 October 2020), which has a variety of high-resolution cartographic resources. To download the digital elevation model (DEM) from which the altitude data were derived, a user account was created to access the platform. After the area of interest was located, the dataset “ALOS PALSAR” was selected, which contains the global DEM at a resolution of 12.5 m obtained through the synthetic aperture radar of the Advanced Land Observation Satellite (ALOS).
  • The climate data were acquired through the website https://app.climateengine.org/climateEngine (accessed on 24 March 2021); after logging into the platform, a series of parameters were defined, such as the dataset, the meteorological variables of precipitation, maximum and minimum temperature, and the time interval (2015–2019). The dataset selected for the study was Terraclimate, which has monthly information on all required variables with a spatial resolution of 4 km. It is the most detailed among the various datasets available in the Climate Engine platform, which covers pixel values between 5 × 5 km and 55 × 55 km. Additionally, this dataset has an array of data from WorldClim, CRUTs 4.0 and JRA 55, which were structured and validated using interpolation and reanalysis techniques [28].
  • Poisson regression analysis: Prior to modelling, univariate Poisson regression models were developed to identify the environmental factors that influence the increase of and spatial variation in the dengue burden in the municipality; a Poisson regression model was run using ArcGis Pro© software, taking into consideration the statistically significant variables between altitude, temperature (maximum and minimum) and precipitation.
  • To identify how the environmental factors influence the increase of and spatial variation in the dengue burden in a village within the municipality, a Poisson regression model was run using WinBUGS [29] and GeoBUGS [30] software packages, both useful for making inferences under a Bayesian framework using the Gibbs Sampling method.
  • WinBUGS approaches Bayesian estimation problems by multiplying a priori distribution by the likelihood and then simulating samples from the a posteriori distribution using the Gibbs algorithm. Predictive maps for the risk of infection were obtained taking into consideration the statistically significant variables such as altitude and minimum temperature.
  • Subsequently, a Hierarchical Bayesian Model (HBM) was built in two stages; the HBM uses multiple levels of analysis in an iterative way [31]. As described by Best [32], in the HBA the unexplained extra-variance found in spatial statistics is identified as either spatially correlated effects or heterogeneity effects.
  • A purely spatial modelling in two stages was followed: at the first stage, a likelihood model for the observed and expected dengue disease counts was specified based on the environmental variables. At the second stage, a prior model over the space of possible relative risks (RR) was specified. The data included two covariates measuring the minimum temperature and elevation within the municipality, and a list of adjacent villages, using the intrinsic conditional autoregressive (CAR) prior proposed by Besag, York and Mollie [33]. This model considers the spatial correlation between neighbouring areas. The general model could be written as:
Oi ~ Poisson(μi),
log (µi) = log Ei + α0 + α1xi/10 + bi,
  • Where α0 is the log relative risk for dengue in the study region, xi is one of the covariates with associated regression coefficient α1, and bi is an area that represents the residual or unexplained relative risk of disease. To allow for spatial dependence between the random effects bi in nearby areas, the car.normal distribution was used. The set of posterior means of relative risks was then used to create a map to visualize high- or low-risk segments.
  • The Bayesian interpretation of probability allows (proper prior) probabilities to be assigned subjectively to random events, in accordance with the natural history of the disease. The Markov Chain Monte Carlo (MCMC) implementation ran a sampling chain for 20,000 iterations, and the first 1000 iterations were discarded as pre-convergence “burn-in” [34]. Spatial autocorrelation was evaluated using Moran’s Index in the residuals. After running the HBM, we use the package CODA and R2WinBUGS to analyse the outputs [29,30,31,32].

2.3.3. Phase 3: Spatial Analysis for the Identification of Epidemiological and Entomological Clusters at the Local Scale, to Prioritize Vector-Control Activities (Objective 3)

  • Spatial analysis with Hot Spot analysis: Prior to hot spot analysis, the Getis-Ord General G statistic was used to identify significant risk clusters. The z-score and p-value are measures of statistical significance that lead to acceptance or rejection of the null hypothesis. For this technique, the null hypothesis states that the values associated with features are randomly distributed. For this analysis, a radius of 100 m was established as threshold distance within the neighbourhood area. This radius was selected based on the dispersal area of the vector [9,35,36] and the results of other similar studies [37,38,39,40]. The potential clusters with the highest burden of dengue per unit area were identified, based on the incidence rate in the neighbourhoods that permanently contribute to the cumulative burden of dengue cases in the municipal seat of Patía.
  • Getis-Ord Gi* test for entomological data: The existence of spatial distribution patterns of the pupal stage of Ae. aegypti vector was evaluated through the study of entomological variables at the neighbourhood scale, such as the total number of pupae, the index of pupae per person and the Breteau index (BI) for spatial analysis in GIS software. Those indicators were selected because the sub-department of health of Cauca (Secretaría de Salud) collected systematically and continuously in conducting this type of entomological surveillance.
  • First, an analysis of the spatial distribution of the vector was performed, with a measure of the degree of clustering (high or low) using the general G statistic of Getis-Ord to indicate whether the pattern was uniformly or randomly clustered based on the fixed distance band parameter. Subsequently, neighbourhoods with groupings or high and homogeneous concentrations of clusters were identified with the Getis-Ord Gi* statistic using 100 m as a distance threshold. The statistical significance was based on the False Discovery Rate (FDR) correction. The information on the variables used was obtained through the Basic Entomology Unit of the Department of Health of Cauca (data available by request).
  • For both spatial analyses we consider “local” in the sense of considering neighbourhood areas within the municipality.

2.4. Ethics Statement

Access to patient data, including the home addresses of the patients, was approved by the Ethics Committee of the National Institute of Health, Colombia (CEMIN 13-2019).

3. Results

3.1. Geographic Distribution of Dengue Incidence

Between 2012 and 2018, the Cauca department showed about 3023 dengue cases throughout the territory. A time-series plot of dengue virus cases (DENV,) in Figure 3, shows dengue cases on the vertical axis, and the time period for the analysis on the horizontal axis. The series shown below is non-stationary and non-linear, with trends that could be associated with a seasonal component.
Figure 4 shows the variation of the incidence rate of dengue for the 42 municipalities of the department of Cauca between 2014 and 2018. Global cluster analysis was performed to represent a variable range between 0.10 and 83.23 cases per 10,000 inhabitants.
In general, the department of Cauca had high incidence rates of dengue in 2014, 2015 and 2018. The epidemiological behaviour of the department was characterized by having 376 cases in 2014; a 64% increase in 2015, with a cumulative incidence rate of 47.04 per 100,000 inhabitants; a gradual decrease in the incidence rate in 2016 and 2017 because of the effect of interventions in the Zika outbreak; and an unusual increase in cases in 2018, with a total of 177 accumulated cases (Figure 2 and Figure 3).
The municipalities in the northeast of the department (Miranda, Puerto Tejada, Villa Rica, Padilla, Corinto, Caloto and Santander de Quilichao) had high incidence rates every year, while the municipalities towards the Pacific coast (Guapi, Timbiquí and Algeria, considering that the latter had very low entomological indices of pupae) had variable and non-persistent incidence rates. Towards the southeastern part of the department, there was an unusual increase in the number of cases in Piamonte (IR = 75.3; POP = 9335), while in the southern region of the department, the municipality of Sucre had the highest incidence rate (IR = 83.23; POP = 9748). The municipality of Patía had continuously moderate to high incidence rates during the observation window (IR = 3.8–28.4; POP = 37,793).

3.1.1. Identification of Clusters of High Risk for Dengue in Cauca

In the retrospective spatio-temporal analysis, the search for clusters with high rates using the discrete Poisson model allowed the identification of four statistically significant geographic clusters with a total of 2697 cases and an annual rate of 29.3 cases per 100,000 inhabitants. This model is useful when count data are available, where there is a background population from which the cases arise, and under the null hypothesis that the cases of dengue are independent of each other.
The first cluster was obtained for 2018, with the highest relative risk of 26.25, while the second cluster, which included the municipality of Patía, presented a relative risk of 8.17 (p ˂ 0.05) specifically for the period between January 2014 and December 2016. In this time frame, the observed cases (n = 293) were higher than the expected cases (n = 39.65) for these municipalities. The third cluster had a relative risk of 4.83 between 2013 and 2015 and was statistically significant, while the fourth cluster presented a value of 1.49 for 2016 (p > 0.1). Table 1 shows the results for each of the clusters from the analysis of the number of registered cases relative to the number of expected cases. A graphical representation of the four clusters is shown in Figure 5.

3.1.2. Getis-Ord Gi* Dengue Hotspot Detection

When analysing the spatial distribution of the incidence rate of dengue for the years 2014 to 2018 in the department, a hot spot analysis identified significantly high rates (≥95% confidence) in the municipalities of Miranda and Patía (year 2014), and the municipalities of Patía and Puerto Tejada (year 2015). No municipalities were identified as “cold” areas where the incidence rate was significantly low.
In the global analysis for 2014 and 2015, the presence of clusters was evidenced for the variable incidence rate (p < 0.1; z-score ≥ 1.65). The spatial pattern of the disease was distributed in an aggregated way that formed clusters, as the spatial analysis looks for patterns and meanings. In Figure 6, the dark and light red colours indicate hotspots of dengue cases (z-score Getis-Ord > 2.58 statistically significant). The blue and light blue colours represent cold spot areas (z-score Getis-Ord < 2.58 statistically significant).
The presence of the municipality of Patía in the cluster analysis persisted until 2017, when it presented significantly high incidence rates, along with the municipality of Piamonte. However, the risk clusters identified in that year were not statistically significant (p > 0.1; z-score: 0.023974). For 2018, the municipalities of Piamonte and Santa Rosa showed high incidence rates, although the presence of Ae. aegypti was not reported in the municipal seat of Santa Rosa. This difference was not statistically significant (p > 0.1; z-score ˂ 1.65).
Once the geographical variation of the event was established, the areas within the municipality with a higher than expected incidence of the disease at the section and neighbourhood scales were identified. This analysis confirmed geographic variation with respect to the incidence rate variable, which was evidenced by identifying the neighbourhoods with higher or lower proportions of events.

3.2. Spatial Variation of the Probability of Dengue Incidence as a Function of Environmental Variables in the Municipality of Patía

In the epidemiological targeting for Patía, the municipality reached significant incidence rates of dengue in 2014, 2015, 2016 and 2018. During the latter year, the incidence rate was 6.5 per 10,000 inhabitants, with a total of 24 accumulated cases (Figure 2 and Figure 3). In 2015, 41% (n = 72) of the total number of cases in the period were analysed; meanwhile, 2016 gathered 32% (n = 56), followed by 2018 and 2019 with 11% (n = 19) and 10% (n = 17) respectively; the year with the lowest number of cases was 2017, with 7% (n = 12). Figure 7 shows the location of dengue cases registered in the period 2015–2019.
These incidence fluctuations may reveal a relationship between dengue disease and the explanatory environmental variables, increasing the probability of finding dengue and its vectors in the villages and population centres analysed, which will allow the carrying out of a characterization process, which is secondary to the targeting process and is useful for identifying and describing the main epidemiological and environmental variables that shape the dynamics of dengue transmission in one of the prioritized clusters of the disease at the departmental scale.
Using a count (Poisson) model type for our dependent discrete variable allowed us to identify the elevation (β = 0.0197) and the minimum temperature (β = 45,185) data as statistically significant (p ≤ 0.01). This results indicate that, in the municipality of Patía, dengue cases are positively influenced by altitude and minimum temperature. Both variables had a significant z-statistic (elev. = 11.2427; min.temp. = 10.4163), meaning that they could contribute significantly to the model.
The identification of risky villages or towns within the municipality requires spatial analysis using GIS and Hierarchical Bayesian Modelling (HBM). It is common knowledge among vector-borne disease researchers that elevation is a risky factor for dengue, but the relationship between the number of cases and minimum temperature may still exist within each particular region. Figure 8 shows the changes in the relative risk (mean) of dengue infection through a HBM as a function of environmental variables: altitude and minimum temperature. This posterior relative risk map clearly shows the characteristic Bayesian smoothing of the crude relative risks.
Additionally, the Caldas-Lang Climate Classification shown in Figure 8 is a classification method widely used in Colombia to characterize the climate, since it allows observation of the general behaviour of temperature as a function of the altitude and humidity of a certain region. This figure shows the risk levels divided into seven levels from lower to higher; a high probability of dengue disease is expected in 18 of 104 villages, which are located towards the northwest of the municipality, and in some isolated villages such as Piedra-Sentada, el Estrecho and El Bordo.
The results of the characterization process show that although the probability of dengue infection for the years 2015 and 2019 as a function of the minimum temperature and altitude variables was higher in the northwest part of the department, some population centres, including El Bordo-Patía, were influenced by other variables, such as biological and sociodemographic variables. Table 2 shows the posterior means and posterior standard deviations of the Poisson regression coefficients related to dengue cases. The model run on WinBUGS software with 20,000 iterations starting at 1001.
From the relative risk estimations displayed in Table 2, analysis reveals that RR estimations at the northwest of the municipality are greater than those near the Patía River valley. The data displayed provide a picture about the application of the spatial model, such as BYM model, that is even more realistic than that the frequentist approach.
The spatial autocorrelation was evaluated with the regression residuals to assess if the Poisson model was correctly specified. Figure 9 shows the cross-correlations between variables in Markov Chain Monte Carlo output showing no spatial autocorrelation for the Markov chain (Moran Test ˂ 1.0). Trace and kernel density for priors (Figure 10) are both available as separate plots, but they are available together via the plot method using the CODA package through R. The plot shows convergence in the case of parameter b1 as well as the other parameters of b. However, the lack of convergence in the case of α2, specifically between the range of 10,000 and 15,000 iterations, may be a wrong sign about convergence if one takes into account that the density distribution is rather irregular. Data used for cross-correlations and trace and kernel density for nodes are contained within the supplementary material section.
The BYM-Poisson model is a good hierarchical distribution to model our data because of the nature of the outcome vatiable. Hence, one of the advantages of this model allows the identification of spatially aggregated count data when incidence is not so high. There are several methodological alternatives for generating estimates with Bayesian hierarchical models and the accuracy of the risk map could be increased for more detailed comparisons by including more explanatory variables. The BYM-Poisson model indicates that, in the municipality of Patía, dengue cases are positively and spatially influenced by altitude and minimum temperature. The results of the estimation can be used as the reference to anticipate the spread of dengue in a municipality, i.e., Patía. The strategic decision and action must be implemented in critical months and must be carried out by the ETV Programme and the Inspection, Surveillance and Health Control Process at Cauca to achieve a more optimal prevention.

3.3. Determination of Disease and Entomological Cluster in the Municipal Seat of Patía

Table 3 shows the result of the total data with the cases by origin, discarded cases, and those that were integrated into the geocoding or georeferencing process. Of the 246 cases used for assigning latitude and longitude, 66 cases were discarded for reasons such as duplication of information between neighbourhood and village, addresses that could not be found, or different origins. Of the total registered cases, 13% were displaced, and six of these cases were reported in Popayan as their municipality of residence and notification; this municipality is the capital city of Cauca (economic and administrative centre of the department).
In total, 180 cases of dengue, including six cases of severe dengue in the municipality of Patía that occurred between 2015 and 2019, were georeferenced (n = 167) and geocoded (n = 13). Of the 180 cases of dengue, 59.9% were men and 40.1% were women. The age range of those affected by the disease ranged from 1 to 89 years, with the following breakdown: 1–5 years, 6.2%; 6–11 years, 16.4%; 12–18 years, 15.8%; 19–26 years, 19.2%; 27–59 years, 35%; and 60 years and older, 7.3%.
After carrying out the global analysis, the −1.65 ≥ z-score ≤ 1.65 pattern did not appear to be significantly different than random at the local scale. Although the calculation of the Getis-Ord General G statistic (=0.001909) for the study period did not allow us to establish a clustering pattern (p < 0.1780; z-score: 1.3468), nine neighbourhoods were found to have more than 15 cases per 1000 habitants and four neighbourhoods more than 23 cases per 1000 habitants (Figure 11).

Spatial Pattern Analysis of Entomological Variables

Regarding the immature stage of the species, for the accumulated pupae for the years 2017 to 2019, the observed general G of Getis-Ord was greater than zero (G = 0.700066), showing that the variable has an aggregate/cluster distribution pattern. The null hypothesis test of complete spatial randomness (CSR) was rejected (Z coefficient = 1.808633, p < 0.1). This means that the spatial distribution of the high vector distribution values in the study area was more spatially clustered than would be expected if the underlying spatial distribution processes were random.
Figure 12 shows the neighbourhoods with statistically significant values. The analysis of hot spots by applying the Getis-Ord Gi* statistic to the 26 neighbourhoods showed that the Libertador neighbourhood has the highest spatial clustering pattern for the distribution of pupae applying the FDR correction (Gi p-value = 0.001; nNeighbors = 22). Additionally, when using the critical p-values and z-scores, the Modelo and Olaya Herrera neighbourhoods had a statistically significant pattern, with 90% confidence (p ≤ 0.1).
Figure 13 shows spatial autocorrelation for a series of distances and creates a line graph of those distances and their corresponding z-scores. Z-scores reflect the intensity of spatial clustering, and statistically significant peak z-scores indicate distances where spatial processes promoting clustering are most pronounced. For the accumulated pupae variable, we did not find peak distances statistically significant vis-à-vis the radius parameter established.
However, one peak at 536.74 m showed the lower p-value, p ≤ 0.1 (z-score = 1.3622).
For the Breteau Index (pupae), the cluster analysis by year did not show statistically significant patterns. However, the neighbourhoods with the highest BIs were Galán (BI = 18.18) in 2017, Prados del Norte (BI = 11.21) in 2018 and El Lago (BI = 12.5) in 2019. For the cumulative years (2017–2019), the highest BIs (BI > 8.33) were reported for the neighbourhoods Calle Nueva, Modelo, Libertador, Las Ferias and La Floresta. Table 4 shows a comparison of the cumulative values of pupae, the BIs and the cases reported during that period.

4. Discussion

This study modelled dengue data in the department of Cauca (Colombia) using spatio-temporal models, correlating the number of cases to environmental variables and estimating the relative risk at the municipality scale. Global and local measures of autocorrelation performed better than spatio-temporal measures when analysing the incidence rate of dengue. Dengue infection risk maps indicated the municipalities were above and below the reference risk, which allows the identification of how elevation and temperature can be used to estimate the probability of disease incidence.

4.1. Epidemiological Behaviour of High-Risk Clusters for Dengue in Cauca

Spatio-temporal clusters using the discrete Poisson model and the Getis Ord-Gi* spatial analysis technique have some similarities and differences. The identification of four clusters through retrospective spatio-temporal analysis limited the search for clusters with high rates to the municipalities of (1) Piamonte; (2) Patía and Sucre; (3) Miranda, Puerto Tejada, Villa Rica, Padilla, Corinto, Caloto, Santander de Quilichao and Guachené; and (4) López de Micay, Timbiquí, Suárez, El Tambo and Argelia. While this analysis allowed the relative risk of the clusters to be determined, the analysis of hot spots was not only consistent with the incidence rate patterns and identified the years when the rates were significantly high (≥90% confidence).
In 2018, Mutheneni et al. conducted a spatial analysis with Getis-Ord Gi*, as this was the method that showed the best autocorrelation [41]. The G Index of Getis and Ord [42] help to identify the degree to which the units of analysis with high values (hot spots) or low values (cold spots) are grouped; that is, it prioritized the formation of clusters.
The spatial distribution pattern of dengue cases was significantly clustered and identified dengue hot spots in the department of Cauca. A consistent dengue hot spot during the study period was identified in the municipality of Patía (≥95% confidence). The other significant hot spots (Miranda, Puerto Tejada and Villa Rica) were mainly located in the northeastern regions of the department near the Valle del Cauca, as Cali is the third most important city in Colombia and is a place with hyperendemic transmission of dengue [39]. The size of the region is related to the geographical characteristics of the territory, since the hot spots lie in highly transited areas along the Pan-American Highway, a commercial and human corridor that connects the departments of southwestern Colombia.
Regarding the space-time scan statistic and Getis-Ord technique, most proposed tests for spatial clustering are tests for global clustering. Some of those methods test for clustering throughout the study area without the ability to identify the location of specific clusters. As such, these tests and the spatial scan statistic complement each other, since they are useful for different purposes.
Additionally, it was observed that in the population of the municipality of Patía, the number of cases of dengue infection was slightly higher among men than among women. These results are similar to those of other studies [41,43,44,45] and contrast with those found in Molineros et al. and Restrepo et al. [46,47]. Additionally, the highest number of cases was found among those aged 27 to 59 years, followed by those aged 19 to 26 and 6 to 11 years. This could be because there was a higher number of cases reported in these age groups in this population. Additionally, the greater number of cases in these age groups may be due to the circulating dengue virus serotype and its mobility patterns.

4.2. Importance of the Environmental Variables Elevation and Minimum Temperature in the Prediction of Dengue

Elevation was the variable with the strongest relationship with the incidence of dengue (μ α1 = −0.02165). Mena et al. and Gyawali et al. (2021) showed that elevation is negatively associated with the incidence of dengue in Costa Rica and Nepal, respectively [48,49]. In Colombia, similar results have been reported. Vásquez found that altitude is a major environmental variable in the incidence of dengue, as increased altitude decreases the risk of dengue by between 60 and 89% within the altitudinal range in Cundinamarca, Colombia [50].
In Colombia, the vector Ae. aegypti has been reported at up to 2302 m.a.s.l. However, its infection with dengue virus has been indicated at a maximum height of 1984 m.a.s.l., in Bello, Antioquia [7]. In the municipality of Patía, the altitude ranges from 536 to 3264 m.a.s.l., and the areas with the highest probability of infection are located at El Bordo and El Estrecho, which is consistent with lower elevation values. Altitude has been indicated as a macroenvironmental factor that limits and influences the development of the vector and the virus [51]. Additionally, altitude is a modifying agent of the microclimate [50,52].
In this study, the minimum temperature was significantly associated with the incidence of dengue in the time window (μ α2 = −3.194), and 2016 and 2018 had minimum yearly temperatures of 17.68 and 17.62 °C, respectively. Other studies have associated the minimum temperature with a higher incidence of dengue. Tuladhar et al. found a greater correlation of dengue incidence with the minimum temperature than with the maximum temperature in Chitwan, Nepal [53]. Similar results were obtained in southern Taiwan [54]; in Mexico, where there is a rapid increase in risk when the average minimum temperatures rise above 18 °C [55]; and in Cali, Colombia, where a strong association has been established between minimum temperature and dengue outbreaks. These studies suggest that periods in which extreme daily temperatures are limited to the range of 18–32 °C promote the growth of the vector population, the amplification of the virus, increased vector feeding, and increased contact with the human host [56].
Carrington et al. performed a vector competence experiment that showed that temperature fluctuations between 18 °C and 20 °C promote more rapid dissemination of the dengue virus than constant temperatures of 20 °C, as only 18.9 days are required for 50% of the mosquitoes exposed to the virus to develop a disseminated infection [57]. These results indicate a greater potential for dengue virus transmission at low temperatures with natural fluctuations and an accelerated rate of dissemination. Although constant temperatures of 26 °C showed similar dissemination results, mortality at this temperature was higher than that at the fluctuating temperature. An increase to a temperature of 30 °C becomes harmful for the mosquito because such temperatures can increase mortality and affect the reproductive function of adults [58].
According to the climatic classification established by the Caldas–Lang climate zoning guidelines, the 50.9% of the municipality of Patía is located in the warm semi-humid climate, 17.4% in the temperate semi-humid zone, 12.7% in the warm semi-arid zone, 12.3% in the temperate wet zone, and the remaining percentage is distributed between cold wet, cold superhumid, and highly cold superhumid climates typical of altitudes above 2000 m.a.s.l. Based on these characteristics, it can be suggested that warm semi-arid and warm semi-humid climates, which are in altitude bands below 1000 m.a.s.l, are directly related to a high probability of dengue transmission. Additionally, the 125.6 km2 extension, in which the temperature ranges from a temperate semi-humid climate between 17.5 °C and 24 °C through the altitudinal strip 1000 to 2000 m.a.s.l., could promote dissemination of the virus more effectively. This confirms the altitudinal variation in temperature and the trends of dengue disease transmission in warm areas, which is limited even in temperate areas, as evidenced in the present study. Similarly, the probability of infection in 2016 and 2018 was related to the minimum temperature (17.68 and 17.62 °C), which was not necessarily the highest among the years studied, as would be expected in dengue transmission scenarios. However, it was consistent with the optimal temperature for vector competence.
Finally, although the relationship between precipitation and dengue was not significant, the minimum temperature, maximum temperature, and precipitation values recorded were consistent with the phases of El Niño and La Niña. The 2015 and 2019 years had the highest temperature values, consistent with the warm phase (28.51 and 29.23 °C, respectively), while 2017 had the highest average precipitation value (2160 mm), due to the cold phase, among other factors.
A study conducted in Mérida, Mexico, found that rainfall variation was the main variable that explained seasonal changes in the abundance of Ae. aegypti and in cases of dengue in this locality [59]. Similarly, precipitation is an environmental variable that has been related to an increase in dengue cases in different countries such as Puerto Rico, Thailand, and Venezuela [60,61]. Precipitation can influence the transmission of dengue through its impact on the vector population, since this variable generates a greater abundance of breeding sites and stimulates the hatching of eggs [60].
Regarding the El Niño Southern Oscillation (ENSO) climate phenomenon, which consists of a warm phase known as El Niño and a cold phase known as La Niña, data from the National Oceanic and Atmospheric Administration (NOAA) show that from 2015 to 2019, the year with the most intense and prolonged warm phase (El Niño) was 2015, which began in October 2014 and ended in April 2016 and was followed by another less intense warm phase from September 2018 to June 2019. Furthermore, the cold phases (La Niña) during this period were shorter and less intense. The most significant phase occurred between October 2017 and April 2018, while a less significant phase occurred between August and December 2016. The highest precipitation within the analysed time window was recorded in 2017.
The implementation of the spatial BYM model in estimating the RR resulted in significant differences of dengue transmission in Patía (Cauca). This analysis also offers the advantage of using count of disease cases when incidence is not so high. With this type of model, the residual relative risks in nearby areas are expected to be more similar than in faraway areas. Further research can be conducted by choosing the prior distribution with other parameter values such as sociodemographic or entomological values.

4.3. Relationship between the Dengue Disease in the Municipal Seat of Patía and the Spatial Distribution Pattern of the Vector

The distribution pattern disease does not appear to be significantly different than random. However, the neighbourhoods with the highest incidence rate corresponded to neighbourhoods with the highest number of accumulated dengue cases during the study period: Olaya Herrera, Libertador, Postobón, and Popular. The same was not true for the neighbourhoods Hueco Lindo and Calle Nueva, which, despite not having a significant number of cases, had a high density of cases due to the superposition of the values of the pixels (10 m × 10 m each), which is more noticeable in small and contiguous spaces such as those encompassed by these neighbourhoods.
It is important to mention that the neighbourhoods that composed the risk clusters identified using kernel density (see Supplementary Materials, Figure S2) had a greater number of inhabitants than the neighbourhoods that did not represent any risk, which can affect the obtained results according to the provisions of Khormi and Kumar [62]. The Olaya Herrera neighbourhood features open spaces such as cemeteries and market squares, which could favour clustering and the presence of breeding sites, as could the neighbourhood’s proximity to rural areas. Other explanatory factors, such as access to basic sanitation services, water use and the disposal of water containers and the mobility of the population, should be considered in future studies.
Furthermore, the analysis of hot spots identified neighbourhoods with extreme and geographically homogeneous values in terms of the total number of pupae: Libertador, Olaya Herrera, and Modelo. These neighbourhoods had higher numbers of pupae than the other neighbourhoods; similarly, contiguous neighbourhoods had significant values, especially in the case of Libertador. One of the most important factors in this analysis was the critical distance of 100 m, which was defined according to the dispersal of the vector. Although this analysis focuses on the pupal stage, the threshold is consistent with other similar studies. Garelli et al. [63] found a clustering of pupae was found in a maximum radius of 150 m. The concentration of the high values found is of great relevance, since the spatial grouping of pupae is usually weak, as described by Khormi et al. and Garelli et al. [62,63]. In addition, this result may be associated with the presence of adult mosquitoes, as indicated by LaCon et al. [64].
Although we did not find a peak by incremental spatial autocorrelation, this most often happens in cases where data have been aggregated and the scale of the processes affecting the input field variable are smaller than the aggregation scheme.
When comparing the results of the analyses of the epidemiological and entomological data mentioned above, it is observed that the neighbourhoods of Libertador and Olaya Herrera simultaneously present the highest density of dengue cases and the highest pupae concentration. However, the correlation between the total number of cases and the total number of pupae registered in the time window of 2017–2019 was weak and non-linear, as demonstrated by the Spearman coefficient (0.1893).

4.4. Study Limitations

The findings of this study have to be seen in light of some limitations. First, in Phase 1, some of the cases were georeferenced using the neighbourhood centroid; we recommend to health authorities the use of tools that allow georeferencing the actual residence address of the case. Second, in Phase 2, the effect estimates in the model are based on the basis of a purely spatial model; even though some of the most important environmental variables were included, they could be improved further with entomological data. Third, in Phase 3, due to the fact that the surveillance of immature stages in the department of Cauca started in 2017, we did not use the same timeframe for epidemiological and entomological clusters.
An additional risk worth mentioning is that of security and public order in the department, as this may have an impact on the collection of primary information in the field. In addition, displacement in conflict zones can be a limitation for primary data collection, so it is recommended to have the support of stakeholders at the local level.

5. Conclusions

Previously, the dengue transmission risk stratification process consisted of the formation of socioecological and epidemiological strata, defined according to the distribution and frequency of risk factors and the endemic and hyperendemic transmission patterns at the local scale [8]. Currently, dengue risk stratification involves five operational scenarios that aim to explain and/or predict how environmental, sociodemographic, and entomological–epidemiological variables will increase or reduce the incidence of the disease [65].
This study is the first to report high-risk clusters of dengue for the department of Cauca based on spatial analysis techniques and considering the trends of dengue disease over a period of seven years. It also evaluated the influence of other variables, such as elevation, minimum temperature, maximum temperature, and precipitation over dengue behaviour in Patía, Cauca. Between the years 2015 and 2019, the main variables that were related to the presence of dengue were altitude and minimum temperature—findings that are consistent with the geographical and climatic conditions of the municipality. Although the effects of elevation on the incidence of dengue are widely known, the minimum temperature is emerging as an important variable, as it is positively related to temperature fluctuations and infection in the vector.
At the local scale, the study established an association between dengue conglomerates and the accumulation of Ae. aegypti pupae in the municipal seat of Patía. The neighbourhoods with the highest kernel density (Olaya Herrera, Libertador, Postobón and Popular) correspond to the neighbourhoods with the highest number of cumulative dengue cases during the study period. Similarly, the Libertador neighbourhood showed extreme and geographically homogeneous values in terms of the cumulative total pupae.
To determine the most cost-effective intervention strategy, it is necessary to construct operational scenarios. In this study, as a first step, the localities were stratified, taking the municipality of Patía as an example. In addition, information was obtained for environmental variables and entomological–epidemiological antecedents. Considering that scenarios should be based on the risks of transmission, the municipality of Patía is characterized as an area with high transmission of dengue virus within the department of Cauca. This is because this municipality has: (1) urban areas with a higher-than-average incidence in at least three of the last five years, according to the data generated by the epidemiological surveillance system; (2) there are established populations of the vector, (3) persistent transmission and (4) various outbreaks during the year, with seasonal behaviour, and (5) cases of severe dengue have been reported.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/tropicalmed8050262/s1, Figure S1: Probability of dengue infection based on environmental variables altitude, minimum temperature and precipitation, years 2016–2018 (Patía, Cauca); Figure S2: Dengue fever case density and risk clusters in Bordo-Patía, Cauca (2015–2019). Table S1. Poisson regression based on environmental factors data (original dataset). Code S1. Computer code with the model and data used in GeoBUGS format. Code S2. Computer code for use the package CODA and R2WinBUGS to analyse the spatial autocorrelation and kernel density outputs.

Author Contributions

Conceptualization, C.M.-D. and E.S.; methodology, C.M.-D. and J.A.S.; software, C.M.-D. and M.C.L.; validation, E.S. and J.A.S.; formal analysis, C.M.-D.; E.S. and J.A.S.; investigation, E.S., C.A.M., M.C.L., P.F. and C.M.-D.; resources, E.S, P.F. and C.A.M.; data curation, C.M.-D., M.C.L. and P.M.-L.; writing—original draft preparation, C.M.-D., M.C.L. and P.M.-L.; writing—review and editing, E.S., H.C., J.A.S., P.F. and C.A.M.; visualization, M.C.L. and C.M.-D.; supervision, E.S.; project administration, E.S. and C.M.-D.; funding acquisition, E.S. and C.M.-D. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Institute of Health, Bogotá, Colombia (research project CEMIN 13-2019), the Secretaria de Salud Departamental del Cauca, the Universidad de Ciencias Aplicadas y Ambientales-UDCA and Minciencias (research project 210484467217).

Institutional Review Board Statement

The study was conducted in accordance with the local health authorities and with the presence of local health agents, and had the scientific licenses number provided by INS (CEMIN 13-2019).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Most of the relevant data are within the manuscript and its Supporting Information files. Additional data presented in this study are available on request from the corresponding author. Georeferenced data presented in this study are openly available in [Kaggle] at [Stratification of Dengue in Cauca—Colombia].

Acknowledgments

We wish to express our gratitude to the officials of the Department of Health of Cauca, especially to Hernando Gil Gómez, leader of the Inspection, Surveillance and Health Control Process, for his support in the realization of the research project; to the leader of the Public Health Surveillance Process, Duban Quintero, for providing us with the databases that were used for georeferencing the cases; to the engineer Anderson Hair Piamba in the ETV Programme for his support; and to the ETV technicians in the municipality of Patía for the collection of the immature stages of the Ae. aegypti vector. Additionally, we thank the reviewers for their constructive comments, which helped us to improve the manuscript, and to the epidemiologist Paola Elena León for her advice in spatial analysis.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Padilla, J.C.; Lizarazo, F.E.; Murillo, O.L.; Mendigaña, F.A.; Pachón, E.; Vera, M.J. Transmission Scenarios of Major Vector-Borne Diseases in Colombia, 1990–2016. Biomédica 2017, 37, 27–40. [Google Scholar] [CrossRef]
  2. Ye, J.; Moreno, M. Comparing Different Spatio-Temporal Modeling Methods in Dengue Fever Data Analysis in Colombia during 2012–2015. Spat. Spat. Epidemiol. 2020, 34, 100360. [Google Scholar] [CrossRef] [PubMed]
  3. Organización Panamericana de la Salud/Organización Mundial de la Salud. Actualización Epidemiológica Dengue 11 de Noviembre de 2019. Available online: https://www3.paho.org/hq/index.php?option=com_docman&view=download&category_slug=dengue-2158&alias=50965-11-de-noviembre-de-2019-dengue-actualizacion-epidemiologica-1&Itemid=270&lang=es (accessed on 11 November 2019).
  4. Sistema Nacional de Vigilancia en Salud Pública-SIVIGILA, Instituto Nacional de Salud Boletín Epidemiológico Semanal, Semana Epidemiológica 32: 8 Al 14 de Agosto de 2021. Colombia; 2021. Available online: https://www.ins.gov.co/buscador-eventos/BoletinEpidemiologico/2021_Boletin_epidemiologico_semana_32.pdf (accessed on 20 August 2021).
  5. Sistema Nacional de Vigilancia en Salud Pública-SIVIGILA, Instituto Nacional de Salud. Boletín Epidemiológico Semanal, Semana Epidemiológica 52: 20 Al 26 de Diciembre de 2020. 2020. Available online: https://www.ins.gov.co/buscador-eventos/BoletinEpidemiologico/2020_Boletin_epidemiologico_semana_52.pdf (accessed on 30 December 2021).
  6. Suaza, J.; Barajas, J.; Galeano, E.; Uribe, S. Criaderos Naturales Para Aedes aegypti (Linnaeus, 1762) y Aedes albopictus (Sus, 1895) Vectores de Arbovirus En La Ciudad de Medellín (Antioquia, Colombia). Boletín del Mus. Entomológico Fr. Luis Gall. 2013, 5, 18–24. [Google Scholar]
  7. Ruiz, F.; González, A.; Vélez, A.; Gómez, G.; Zuleta, L.; Uribe, S.; Vélez, I. Presencia de Aedes (Stegomyia) aegypti (Linnaeus, 1762) y Su Infección Natural Con El Virus Del Dengue En Alturas No Registradas Para Colombia. Biomédica 2016, 36, 303–308. [Google Scholar] [CrossRef]
  8. Ministerio de la Protección Social & Instituto Nacional de Salud & Organización Panamericana de la Salud. Gestión Para La Vigilancia Entomológica y Control de La Transmisión Del Dengue. Colombia; 2011. Available online: https://www.minsalud.gov.co/sites/rid/Lists/BibliotecaDigital/RIDE/DE/gestion-vigilancia-entomologica-dengue.pdf (accessed on 10 August 2020).
  9. Honório, N.; Costa, W.; Leite, P.; Gonçalves, J.; Lounibos, L.; Lourenço-de-Oliveira, R. Dispersal of Aedes aegypti and Aedes albopictus (Diptera: Culicidae) in an Urban Endemic Dengue Area in the State of Rio de Janeiro, Brazil. Mem Inst Oswaldo Cruz 2003, 98, 191–198. [Google Scholar] [CrossRef] [PubMed]
  10. Padilla, J.C.; Rojas, D.P.; Roberto, S.-G. Dengue En Colombia: Epidemiología de la Reemergencia a la Hiperendemia. Rev. Salud. Bosque. 2015, 5, 81–83. [Google Scholar]
  11. Balaji, D.; Saravanabavan, V. Geo Spatial Variation of Dengue Risk Zone in Madurai City Using Autocorrelation Techniques. GeoJournal 2021, 86, 1481–1501. [Google Scholar] [CrossRef]
  12. Mutheneni, S.R.; Mopuri, R.; Naish, S.; Gunti, D.; Upadhyayula, S.M. Spatial Distribution and Cluster Analysis of Dengue Using Self Organizing Maps in Andhra Pradesh, India, 2011–2013. Parasite Epidemiol. Control 2018, 3, 52–61. [Google Scholar] [CrossRef]
  13. Mala, S.; Kumar, M. Geographic Information System Based Spatio-Temporal Dengue Fever Cluster Analysis and Mapping. Egypt. J. Remote Sens. Sp. Sci. 2019, 22, 297–304. [Google Scholar] [CrossRef]
  14. Oviedo, M.; Lage, R.; Romero, R.; Fonseca, C.; Amaral, J. Spatial and Statistical Methodologies to Determine the Distribution of Dengue in Brazilian Municipalities and Relate Incidence with the Health Vulnerability Index. Spat. Spat. Epidemiol. 2014, 11, 143–151. [Google Scholar] [CrossRef]
  15. Fonseca, F.; Carvalho, M. Spatial and Temporal Analysis of Epidemiological Data. 1996. In: GIS for Health and the Environment. Available online: http://www.idrc.ca/books/focus/766/nobre.html (accessed on 11 September 2020).
  16. Kitron, U. Landscape Ecology and Epidemiology of Vector-Borne Diseases: Tools for Spatial Analysis. J. Med. Entomol 1998, 35, 435–445. [Google Scholar] [CrossRef] [PubMed]
  17. Departamento Nacional de Planeación Ficha Departamento de Cauca. Terridata. 2020. Available online: https://terridata.dnp.gov.co/ (accessed on 11 September 2020).
  18. Gobernación del Cauca Plan de Desarrollo Departamental 2020–2023. Colombia. 2020. Available online: https://obsgestioneducativa.com/download/plan-de-desarrollo-departamental-cauca-2020-2023/ (accessed on 22 November 2019).
  19. Zambrano, P. Protocolo de Vigilancia En Salud Pública. Dengue Código: 210–220–580. Vigil. Análisis Riesgo en Salud Pública Instituto Nacional de Salud; Instituto Nacional de Salud: Bogotá, Colombia, 2022; Volume 4. Available online: https://www.ins.gov.co/buscador-eventos/Lineamientos/Pro_Dengue.pdf (accessed on 15 June 2020).
  20. Rodrígues, T.; da Silveira, I.; Leite, W. Improving Geocoding Matching Rates of Structured Addresses in Rio de Janeiro, Brazil. Cad. Saúde Pública 2021, 37, e00039321. [Google Scholar] [CrossRef]
  21. Wieczorek, J.; Guo, Q.; Hujmans, R. The Point-Radius Method for Georeferencing Locality Descriptions and Calculating Associated Uncertainty. Int. J. Geogr. Inf. Sci. 2004, 18, 745–767. [Google Scholar] [CrossRef]
  22. Departamento Administrativo Nacional de Estadística Marco Geoestadístico Nacional. Geoportal DANE. Available online: https://geoportal.dane.gov.co/servicios/descarga-y-metadatos/descarga-mgn-marco-geoestadistico-nacional/ (accessed on 30 May 2020).
  23. Kulldorff, M. A Spatial Scan Statistic. Commun. Stat. 1997, 26, 1481–1496. [Google Scholar] [CrossRef]
  24. Khormi, H.M.; Kumar, L. Identifying and Visualizing Spatial Patterns and Hot Spots of Clinically-Confirmed Dengue Fever Cases and Female Aedes aegypti Mosquitoes in Jeddah, Saudi Arabia. Dengue Bull. 2011, 35, 15–34. [Google Scholar]
  25. Castillo, K.C.; Körbl, B.; Stewart, A.; Gonzalez, J.F.; Ponce, F. Application of Spatial Analysis to the Examination of Dengue Fever in Guayaquil, Ecuador. Procedia Environ. Sci. 2011, 7, 188–193. [Google Scholar] [CrossRef]
  26. Asociación Supradepartamental de Municipios de la Región del Alto Patía-Asopatía; Corporación Colombiana de Investigación Agropecuaria-Corpoica. Plan Básico de Ordenamiento Territorial Municipio de Patía (Cauca); Instituto Nacional de Salud: Bogotá, Colombia, 2001. Available online: http://www.asopatia.gov.co/ (accessed on 20 July 2020).
  27. Vergara, H.; Torres, P. Aspectos Generales Del Valle de Patía. Rev. Noved. Colomb. 2017, 12, 11–24. [Google Scholar]
  28. Abatzoglou, J.T.; Dobrowski, S.Z.; Parks, S.A.; Hegewisch, K.C. TerraClimate, a High-Resolution Global Dataset of Monthly Climate and Climatic Water Balance from 1958–2015. Sci. Data 2018, 5, 170191. [Google Scholar] [CrossRef]
  29. Spiegelhalter, D.J.; Thomas, A.; Best, N.G. WinBUGS Examples Volume 2. Imp. Coll. Med. Res. Counc. UK 2004, 6, 1–42. [Google Scholar]
  30. Thomas, A.; Best, N.; Lunn, D.; Spiegelhalter, D.; Way, R. GeoBUGS User Manual. GeoBugs 2004, 1, 45. [Google Scholar]
  31. Lawson, A. Using R for Bayesian Spatial and Spatio-Temporal Health Modeling; Kindle, Ed.; Chapman & Hall/CRC The R Series; CRC Press: Boca Raton, FL, USA, 2021. [Google Scholar]
  32. Best, N.; Arnold, R.; Thomas, A.; Waller, L.; Collon, E. Bayesian Models for Spatially Correlated Disease and Exposure Data. Bayesian Stat. 1999, 6, 131–156. [Google Scholar]
  33. Besag, J.; York, J.; Mollie, A. Bayesian Image Restoration, with Two Applications in Spatial Statistics (with Discussion). Ann. Inst. Stat. Math. 1991, 43, 1–59. [Google Scholar] [CrossRef]
  34. Stevenson, M.A.; Morris, R.S.; Lawson, A.B.; Wilesmith, J.W.; Ryan, J.B.M.; Jackson, R. Area-Level Risks for BSE in British Cattle before and after the July 1988 Meat and Bone Meal Feed Ban. Prev. Vet. Med. 2005, 69, 129–144. [Google Scholar] [CrossRef]
  35. Nelson, M. Aedes Aegypti: Biología y Ecología; Pan American Health Organization: Washington, DC, USA, 1986; pp. 1–62. [Google Scholar]
  36. Lourenço de Oliveira, R.; Castro, M.; Braks, M.; Lounibos, L. The Invasion of Urban Forest by Dengue Vectors in Rio de Janeiro. J. Vector Ecol. 2004, 29, 94–100. [Google Scholar]
  37. Rotela, C. Desarrollo de Modelos e Indicadores Remotos de Riesgo Epidemiológico de Dengue En Argentina; Universidad Nacional de Córdoba: Córdoba, Argentina, 2012; Available online: https://rdu.unc.edu.ar/handle/11086/11609 (accessed on 30 June 2020).
  38. Rojas, R.; Nazarena, M. Usos de Herramientas Geoespaciales En La Detección de Áreas Con Riesgo Epidemiológico a Partir de Variables Biofísicas y Casos de Dengue En Jujuy—Argentina; Universidad Nacional de Córdoba: Córdoba, Argentina, 2016; Available online: https://rdu.unc.edu.ar/handle/11086/2841 (accessed on 30 June 2020).
  39. Cuartas, D.E.; Martínez, G.; Caicedo, D.M.; Garcés, J.; Ariza-Araújo, Y.; Peña, M.; Méndez, F. Distribución Espacial de Criaderos Positivos y Potenciales de Aedes aegypti. Biomédica 2017, 37, 59–66. [Google Scholar] [CrossRef]
  40. Souza, R.; Carvalho, M. Análise Da Distribuição Espacial de Larvas de Aedes aegypti Na Ilha Do Governador, Rio de Janeiro, Brasil. Cad. Saude Publica 2000, 16, 31–42. [Google Scholar] [CrossRef] [PubMed]
  41. Mutheneni, S.R.; Morse, A.P.; Caminade, C.; Upadhyayula, S.M. Dengue burden in India: Recent trends and importance of climatic parameters. Emerg. Microbes Infect. 2019, 6, 1–10. [Google Scholar] [CrossRef]
  42. Getis, A.; Ord, K. The Analysis of Spatial Association by Use of Distance Statistics. Geogr. Anal. 1992, 24, 189–206. [Google Scholar] [CrossRef]
  43. Dzul-Manzanilla, F.; Correa-Morales, F.; Che-Mendoza, A.; Palacio-Vargas, J.; Sánchez-Tejeda, G.; González-Roldan, J.F.; López-Gatell, H.; Flores-Suárez, A.E.; Gómez-Dantes, H.; Coelho, G.E.; et al. Identifying Urban Hotspots of Dengue, Chikungunya, and Zika Transmission in Mexico to Support Risk Stratification Efforts: A Spatial Analysis. Lancet Planet. Health 2021, 5, e277–e285. [Google Scholar] [CrossRef]
  44. Chaparro, P.; León-Quevedo, W.; Castañeda, C. Comportamiento de La Mortalidad Por Dengue En Colombia Entre 1985 y 2012. Biomédica 2016, 36, 125–134. [Google Scholar] [CrossRef]
  45. Hernández, M.; Arboleda, D.; Arce, S.; Benavides, A.; Tejada, P.; Ramírez, S.; Cubides, A. Metodología Para La Elaboración de Canales Endémicos y Tendencia de La Notificación Del Dengue, Valle Del Cauca, Colombia, 2009–2013. Biomédica 2016, 36, 98–107. [Google Scholar] [CrossRef]
  46. Molineros, L.F.; Pinzón, E.M.; Rengifo, N.E.; Daza, C.F.; Hernández-Carrillo, M.; Ortiz, M.E.; Lesmes, M.C. Seroprevalencia de Dengue En Municipios Con Transmisión Hiperendémica y Mesoendémica, Valle Del Cauca, Colombia. Rev. Cuba. Salud Pública 2020, 46, 1–20. [Google Scholar]
  47. Restrepo, B.; Arboleda, M.; Lopera, T. Estudio Seroepidemiológico de Dengue En La Región Del Urabá Antioqueño. Infectio 2004, 8, 255–262. [Google Scholar]
  48. Mena, N.; Troyo, A.; Bonilla-Carrión, R.; Calderón-Arguedas, O. Factores Asociados Con La Incidencia de Dengue En Costa Rica. Rev. Panam. Salud Pública 2011, 29, 234–242. [Google Scholar] [CrossRef]
  49. Gyawali, N.; Johnson, B.; Dixit, S.; Devine, G. Patterns of Dengue in Nepal from 2010–2019 in Relation to Elevation and Climate. Trans. R. Soc. Trop. Med. Hyg. 2021, 115, 741–749. [Google Scholar] [CrossRef]
  50. Vásquez, A. Factores Geográficos, Ecológicos y Sociodemográficos En La Ocurrencia de Dengue En Cundinamarca; Universidad Nacional de Colombia: Bogotá, Colombia, 2019. [Google Scholar]
  51. Organización Panamericana de la Salud/Organización Mundial de la Salud. Dengue y Dengue Hemorrágico En Las Américas: Guías Para Su Prevención y Control. 1995. Available online: https://iris.paho.org/handle/10665.2/36861 (accessed on 30 June 2021).
  52. Eslava, J. Climatología y Diversidad Climática En Colombia. Rev. Acad. Colomb. Cienc. 1993, 18, 508–538. [Google Scholar]
  53. Tuladhar, R.; Singh, A.; Varma, A.; Choudhary, D.K. Climatic Factors Influencing Dengue Incidence in an Epidemic Area of Nepal. BMC Res. 2019, 12, 131. [Google Scholar] [CrossRef]
  54. Chien, L.; Yu, H. Impact of Meteorological Factors on the Spatiotemporal Patterns of Dengue Fever Incidence. Environ. Int. 2014, 73, 46–56. [Google Scholar] [CrossRef] [PubMed]
  55. Colón, F.; Fezzi, C.; Lake, I.; Hunter, P. Los Efectos Del Clima y El Cambio Climático Sobre El Dengue. PLoS Negl. Trop. Dis. 2013, 7, e2503. [Google Scholar] [CrossRef]
  56. Eastin, M.; Delmelle, E.; Casas, I.; Wexler, J.; Self, C. Intra and Interseasonal Autoregressive Prediction of Dengue Outbreaks Using Local Weather and Regional Climate for a Tropical Environment in Colombia. Am. J. Trop. Med. Hyg. 2014, 91, 598–610. [Google Scholar] [CrossRef] [PubMed]
  57. Carrington, L.; Armijos, M.; Lambrechts, L.; Scott, T. Fluctuations at a Low Mean Temperature Accelerate Dengue Virus Transmission by Aedes aegypti. PLoS Negl. Trop. Dis. 2013, 7, e2190. [Google Scholar] [CrossRef]
  58. Carrington, L.; Armijos, M.; Lambrechts, L.; Barker, C.; Scott, T. Effects of Fluctuating Daily Temperatures at Critical Thermal Extremes on Aedes aegypti Life-History Traits. PLoS ONE 2013, 8, e58824. [Google Scholar] [CrossRef]
  59. García, J.; Loroño, M.; Farfán, J.; Flores, L.; Del Pilar Rosado, E.; Rivero, N.; Najera-Vazquez, R.; Gomez-Carro, S.; Lira-Zumbardo, V.; Gonzalez-Martinez, P.; et al. Dengue Virus Infected Aedes aegypti in the Home Environment. Am. J. Trop. Med. Hyg. 2008, 79, 940–950. [Google Scholar] [CrossRef]
  60. Johansson, M.; Cummings, D.; Glass, G. Multi-Year Climate Variability and Dengue-El Niño Southern Oscillation, Weather, and Dengue Transmission in Puerto Rico, Mexico, and Thailand: A Longitudinal Data Analysis. PLoS Med. 2009, 6, e1000168. [Google Scholar] [CrossRef] [PubMed]
  61. Rubio, Y.; Pérez, L.; Infante, M.; Comach, G.; Urdaneta, L. Influencia de Las Variables Climáticas En La Casuística de Dengue y La Abundancia de Aedes aegypti (Díptera: Culicidae) En Maracay, Venezuela. Bol. Malariol. Salud Ambient. 2011, 51, 145–158. [Google Scholar]
  62. Khormi, H.M.; Kumar, L. Modeling Dengue Fever Risk Based on Socioeconomic Parameters, Nationality and Age Groups: GIS and Remote Sensing Based Case Study. Sci. Total Environ. 2011, 409, 4713–4719. [Google Scholar] [CrossRef]
  63. Garelli, F.; Espinosa, M.; Gürtel, R. Spatial Analysis of Aedes aegypti Immatures in Northern Argentina: Clusters and Temporal Instability. Acta Trop. 2013, 128, 461–467. [Google Scholar] [CrossRef]
  64. LaCon, G.; Morrison, A.C.; Astete, H.; Stoddard, S.T.; Paz-Soldan, V.A.; Elder, J.P.; Halsey, E.S.; Scott, T.W.; Kitron, U.; Vazquez-Prokopec, G.M. Shifting Patterns of Aedes aegypti Fine Scale Spatial Clustering in Iquitos, Peru. PLoS Negl. Trop. Dis. 2014, 8, e3038. [Google Scholar] [CrossRef]
  65. Organización Panamericana de la Salud/Organización Mundial de la Salud. Documento Técnico Para La Implementación de Intervenciones Basado En Escenarios Operativos Genéricos Para El Control Del Aedes aegypti; Pan American Health Organization: Washington, DC, USA, 2019; Available online: https://iris.paho.org/handle/10665.2/51654 (accessed on 12 February 2020).
Figure 1. Methodological phases of the study design.
Figure 1. Methodological phases of the study design.
Tropicalmed 08 00262 g001
Figure 2. Dengue endemic channel for the municipality of Patía. Yellow zone represents the number of cases within the expected normal seasonal range. Anything above this moving threshold would be considered representative of an unprecedented number of cases in the municipality (Red zone).
Figure 2. Dengue endemic channel for the municipality of Patía. Yellow zone represents the number of cases within the expected normal seasonal range. Anything above this moving threshold would be considered representative of an unprecedented number of cases in the municipality (Red zone).
Tropicalmed 08 00262 g002
Figure 3. Time series of dengue cases from the department of Cauca from 2012–2018.
Figure 3. Time series of dengue cases from the department of Cauca from 2012–2018.
Tropicalmed 08 00262 g003
Figure 4. Dengue incidence rate × 10,000 population for the department of Cauca (2014–2018).
Figure 4. Dengue incidence rate × 10,000 population for the department of Cauca (2014–2018).
Tropicalmed 08 00262 g004
Figure 5. High-risk spatiotemporal clusters of dengue in the department of Cauca (2012–2018).
Figure 5. High-risk spatiotemporal clusters of dengue in the department of Cauca (2012–2018).
Tropicalmed 08 00262 g005
Figure 6. Hot spots for the department of Cauca considering the incidence rate of dengue (2014–2018).
Figure 6. Hot spots for the department of Cauca considering the incidence rate of dengue (2014–2018).
Tropicalmed 08 00262 g006
Figure 7. Observed number of dengue cases in Patía municipality for 2015–2019 period.
Figure 7. Observed number of dengue cases in Patía municipality for 2015–2019 period.
Tropicalmed 08 00262 g007
Figure 8. Poisson spatial BYM regression model for dengue disease based on environmental variables altitude and minimum temperature, years 2015–2019 (Patía, Cauca).
Figure 8. Poisson spatial BYM regression model for dengue disease based on environmental variables altitude and minimum temperature, years 2015–2019 (Patía, Cauca).
Tropicalmed 08 00262 g008
Figure 9. Cross-correlations for Markov chain in Poisson spatial regression model for dengue disease based on environmental variables of altitude and minimum temperature, years 2015–2019 (Patía, Cauca).
Figure 9. Cross-correlations for Markov chain in Poisson spatial regression model for dengue disease based on environmental variables of altitude and minimum temperature, years 2015–2019 (Patía, Cauca).
Tropicalmed 08 00262 g009
Figure 10. Trace and density for nodes in Markov chain for Poisson spatial regression model (2015–2019).
Figure 10. Trace and density for nodes in Markov chain for Poisson spatial regression model (2015–2019).
Tropicalmed 08 00262 g010
Figure 11. Dengue incidence rate per 1000 habitants in Bordo-Patía, Cauca (2015–2019).
Figure 11. Dengue incidence rate per 1000 habitants in Bordo-Patía, Cauca (2015–2019).
Tropicalmed 08 00262 g011
Figure 12. Accumulated pupae of A. aegypti (2017–2019) in the neighbourhoods of Bordo–Patía (left); hot spots for pupae calculated using the Getis-Ord Gi* statistic (right).
Figure 12. Accumulated pupae of A. aegypti (2017–2019) in the neighbourhoods of Bordo–Patía (left); hot spots for pupae calculated using the Getis-Ord Gi* statistic (right).
Tropicalmed 08 00262 g012
Figure 13. Incremental spatial autocorrelation for the accumulated pupae in Bordo–Patía neighbourhoods (2017–2019).
Figure 13. Incremental spatial autocorrelation for the accumulated pupae in Bordo–Patía neighbourhoods (2017–2019).
Tropicalmed 08 00262 g013
Table 1. Spatio-temporal analysis for the identification of high-risk clusters for dengue in the department of Cauca (2012–2018).
Table 1. Spatio-temporal analysis for the identification of high-risk clusters for dengue in the department of Cauca (2012–2018).
Dengue Clusters in the Department of Cauca
Cluster I *Cluster II *Cluster III *Cluster IV
MunicipalitiesPiamontePatía and SucreCaloto, Miranda, Corinto, Padilla, Puerto Tejada,
Guachené, Santander de
Quilichao and Villa Rica
López de Micay, Tim-
biquí, Suárez, El Tambo and Argelia
Time frame20182014–20162013–20152016
Population734345,103274,554134,841
Registered cases5629385759
Expected cases2.1839.65237.2239.77
Relative risk26.258.174.831.49
Likelihood ratio128.58345.14566.614.11
*p < 0.001.
Table 2. Summary of the results of relative risk estimations using BYM model for dengue.
Table 2. Summary of the results of relative risk estimations using BYM model for dengue.
NodeMeanStandard
Deviation
Monte Carlo
Standard Error
2.5%Median97.5%
alpha07.6953.470.29021.7957.50514.39
alpha1−0.021650.0084916.66 × 10−4−0.03802−0.02159−0.006483
alpha2−3.1941.4830.1237−5.948−3.101−0.6001
b [1]−0.54350.60660.01957−1.777−0.53340.6039
b [2]1.340.40150.018270.55721.3372.132
b [3]−0.53950.20340.012392.6723.0573.464
b [4]−0.31360.58440.01519−1.499−0.30070.7992
b [5]−0.10960.51870.01175−1.165−0.099050.8628
b [6]−0.21940.53980.01347−1.334−0.20080.7841
b [7]−0.5720.58910.01551−1.775−0.55460.5265
b [8]0.60910.6280.01041−0.71120.64081.763
b [9]0.017790.59550.008791−1.2020.035561.121
b [10]0.20880.52010.01189−0.84480.22131.185
b [11]0.33560.7690.02283−1.2130.35751.799
b [12]0.22520.58910.005614−0.97380.24111.339
b [13]0.67860.73780.02756−0.80570.6862.082
b [14]−0.16650.50870.01047−1.227−0.14890.7755
b [15]−0.32460.58740.01131−1.532−0.30950.7713
b [16]0.52090.87190.03325−1.2560.54062.169
b [17]−0.26560.81660.01115−2.038−0.21081.168
b [18]0.27780.54130.004316−0.82490.29231.291
b [19]0.55240.92010.03325−1.3010.57252.294
b [20]0.32770.50340.004141−0.68650.33741.296
b [21]0.59870.5830.00776−0.59360.62081.676
b [22]0.69850.79480.02767−0.91440.71182.201
b [23]0.70030.77460.03044−0.83980.70472.194
b [24]0.64070.84430.03292−1.0520.65592.274
b [25]0.66980.93540.03484−1.1990.68452.46
b [26]0.60660.93970.03487−1.2680.6332.4
b [27]0.49070.80110.03154−1.1270.50832.01
b [28]0.26740.73990.02962−1.2170.28981.654
b [29]0.46550.76950.02498−1.1130.49541.902
b [30]1.0990.57080.01492−0.06591.1192.17
b [31]0.55010.66910.01684−0.8130.56871.815
b [32]0.40660.7780.01396−1.2270.43861.83
b [33]−0.09150.75890.7589−1.718−0.034751.237
b [34]0.47840.59760.009939−0.76410.49661.595
b [35]0.8490.70010.02512−0.56240.86762.169
b [36]0.95510.52680.005398−0.13240.97571.924
b [37]0.6060.62360.02026−0.64530.61391.809
b [38]0.81780.53070.008372−0.26790.82931.818
b [39]0.44210.63460.02222−0.85580.45621.663
b [40]0.27330.65470.02012−1.0510.29351.513
b [41]0.73220.58010.01306−0.46730.75221.81
b [42]0.32730.56430.006039−0.82950.34571.381
b [43]−0.3410.7060.01081−1.842−0.30080.9222
b [44]0.050250.66270.006465−1.3340.073331.279
b [45]−0.12340.61030.01085−1.353−0.10651.033
b [46]−0.09730.59370.008167−1.32−0.079681.005
b [47]−0.0097870.51730.005345−1.0680.0053710.9668
b [48]0.18150.51570.00985−0.8610.18921.171
b [49]0.16970.60750.005576−1.0710.18971.316
b [50]0.31640.56310.006296−0.84120.3361.368
b [51]0.098940.7090.01129−1.3460.12111.418
b [52]0.15120.63550.01506−1.1370.16281.357
b [53]0.079440.56450.004516−1.0740.097021.143
b [54]0.3770.65130.02601−0.90650.38211.643
b [55]−0.67580.53080.0199−1.748−0.6610.3264
b [56]0.19520.57870.01601−1.0030.21111.279
b [57]0.21950.61240.01644−1.0240.2361.374
b [58]0.42430.61110.02113−0.78720.43051.605
b [59]0.65990.8340.02483−1.0750.69532.181
b [60]0.14360.60620.01453−1.0760.15991.293
b [61]0.19390.66350.01304−1.1850.21811.416
b [62]−0.019980.71520.02696−1.478−0.011881.346
b [63]0.28540.61760.009941−0.98630.30861.432
b [64]0.36020.69330.0105−1.060.39341.629
b [65]0.84060.75160.01541−0.78250.90242.152
b [66]−0.35310.52540.01345−1.404−0.34170.6417
b [67]−0.30.64030.01245−1.612−0.28080.8996
b [68]−0.21660.56350.0104−1.381−0.20010.8314
b [69]−0.47950.61480.01522−1.742−0.45330.6633
b [70]−0.73090.65640.01606−2.078−0.71530.5015
b [71]−0.38470.60530.01633−1.613−0.36690.7566
b [72]−0.26050.60730.0147−1.523−0.23390.8679
b [73]−0.082150.63190.009134−1.407−0.047731.057
b [74]1.6960.43050.0058310.8021.7142.492
b [75]−0.07930.56850.009272−1.267−0.054480.9694
b [76]−0.29720.67270.01103−1.692−0.26650.9358
b [77]−0.042710.59720.009422−1.281−0.017511.048
b [78]−0.24960.55430.00934−1.413−0.22060.7576
b [79]−0.12820.55720.01277−1.305−0.10670.8917
b [80]−0.40820.51740.01142−1.498−0.38590.5412
b [81]−0.73110.61990.01116−2.031−0.69840.3958
b [82]−0.74440.54250.01338−1.871−0.71510.2382
b [83]−0.28320.93970.01321−2.357−0.19851.312
b [84]−1.2560.68060.01606−2.659−1.228−0.009965
b [85]−0.88980.6550.01436−2.256−0.86080.3195
b [86]−1.2940.55510.01508−2.444−1.27−0.27
b [87]−1.3480.74380.01577−2.915−1.310.007691
b [88]−0.41580.67080.0184−1.799−0.39650.8323
b [89]−0.053540.61860.01769−1.314−0.041621.119
b [90]−1.0610.57910.01582−2.244−1.0460.03023
b [91]−1.2650.70970.01687−2.74−1.230.02105
b [92]−0.81220.59550.01683−2.022−0.79830.3181
b [93]−0.86150.68130.01718−2.263−0.83140.4095
b [94]−0.83090.58830.01738−2.03−0.81360.2922
b [95]−0.77910.61030.01846−2.038−0.760.3611
b [96]−1.0850.67010.01851−2.488−1.0570.1453
b [97]−0.26660.62030.01904−1.56−0.250.9118
b [98]−1.0020.71060.01762−2.467−0.9760.3117
b [99]−0.71580.5730.01842−1.878−0.70470.3791
b [100]−0.93110.67180.01815−2.343−0.89140.2785
b [101]−0.54150.83620.01919−2.312−0.49270.9519
b [102]−0.30210.60520.01899−1.546−0.28370.8592
b [103]−0.54780.66630.0198−1.923−0.51850.6852
b [104]−0.82850.73590.02025−2.37−0.79770.5246
b [105]3.060.20340.012392.6723.0573.464
tau0.37710.082710.0012860.24070.36880.5622
Table 3. Dengue cases registered in the municipality of Patía by source. Dengue cases from the a national and b department level.
Table 3. Dengue cases registered in the municipality of Patía by source. Dengue cases from the a national and b department level.
Dengue Cases in the Municipality of Patia (2015–2019)
SIVIGILA Notified aCAUCA
Notified b
ProvenienceDismissFinal Data
Year
2015101971091297
2016658691388
2017141423914
2018242431625
2019172228622
Total22124328236246
Table 4. Dengue cases and Pupae Index in Bordo-Patía neighbourhoods, disaggregated by year (2017–2019).
Table 4. Dengue cases and Pupae Index in Bordo-Patía neighbourhoods, disaggregated by year (2017–2019).
NeighbourhoodsPopulationYear 2017Year 2018Year 2019
No. CasesTotal PupaeBreteau
Index
Pupae/PersonNo. CasesTotal PupaeBreteau
Index
Pupae/PersonNo. CasesTotal PupaeBreteau
Index
Pupae/Person
Altillo1850347.220.130182.780.060546.180.22
Aruba y Curazao178000.000.00000.000.00200.000.00
Balboita365092.380.02000.000.001172.080.04
Calle Nueva2910673.570.120351.920.100200.690.02
Centro263021.590.01000.000.000354.170.15
El Campín15501755.930.310182.500.0601728.040.26
El Lago523094.170.040127.140.09110212.500.28
El Peñol4160425.360.100335.410.120604.250.12
Estudiantes3680802.780.090324.170.0601355.090.14
Fundadores4190815.290.220174.780.060884.170.15
Galán19904918.180.54117.55.560.11000.000.00
Hospital8000.000.00000.000.00000.000.00
Hueco Lindo3591102.54.290.190452.380.140965.000.13
Jardín3030631.850.09000.000.000292.780.03
La Floresta28104315.380.410127.140.10000.000.00
La Unión3481986.210.211337.140.100554.760.09
Las Ferias86016314.290.490152.780.0701089.440.21
Libertador583141015.930.4701458.330.3301384.700.15
Limonar163000.000.00000.000.000142.500.05
Modelo36901248.250.220699.520.22028211.550.42
Olaya Herrera14610983.270.093124.57.920.1702247.320.18
Pablo VI–Parte I187000.000.00000.000.000615.950.12
Pablo VI–Parte II513000.000.00000.000.00000.000.00
Popular29421107.250.21000.000.001995.170.14
Postobon71701146.230.162151.430.030603.660.07
Prados del Norte3271332.380.1204211.210.201775.140.14
Puerto Nuevo1980605.810.23000.000.000211.960.06
Versalles1300485.560.170183.330.081334.580.11
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Marceló-Díaz, C.; Lesmes, M.C.; Santamaría, E.; Salamanca, J.A.; Fuya, P.; Cadena, H.; Muñoz-Laiton, P.; Morales, C.A. Spatial Analysis of Dengue Clusters at Department, Municipality and Local Scales in the Southwest of Colombia, 2014–2019. Trop. Med. Infect. Dis. 2023, 8, 262. https://doi.org/10.3390/tropicalmed8050262

AMA Style

Marceló-Díaz C, Lesmes MC, Santamaría E, Salamanca JA, Fuya P, Cadena H, Muñoz-Laiton P, Morales CA. Spatial Analysis of Dengue Clusters at Department, Municipality and Local Scales in the Southwest of Colombia, 2014–2019. Tropical Medicine and Infectious Disease. 2023; 8(5):262. https://doi.org/10.3390/tropicalmed8050262

Chicago/Turabian Style

Marceló-Díaz, Catalina, María Camila Lesmes, Erika Santamaría, José Alejandro Salamanca, Patricia Fuya, Horacio Cadena, Paola Muñoz-Laiton, and Carlos Andrés Morales. 2023. "Spatial Analysis of Dengue Clusters at Department, Municipality and Local Scales in the Southwest of Colombia, 2014–2019" Tropical Medicine and Infectious Disease 8, no. 5: 262. https://doi.org/10.3390/tropicalmed8050262

APA Style

Marceló-Díaz, C., Lesmes, M. C., Santamaría, E., Salamanca, J. A., Fuya, P., Cadena, H., Muñoz-Laiton, P., & Morales, C. A. (2023). Spatial Analysis of Dengue Clusters at Department, Municipality and Local Scales in the Southwest of Colombia, 2014–2019. Tropical Medicine and Infectious Disease, 8(5), 262. https://doi.org/10.3390/tropicalmed8050262

Article Metrics

Back to TopTop