(C) PLOS One
This story was originally published by PLOS One and is unaltered.
. . . . . . . . . .



Socio-economic factors and its influence on the association between temperature and dengue incidence in 61 Provinces of the Philippines, 2010–2019 [1]

['Xerxes Seposo', 'Department Of Hygiene', 'Hokkaido University', 'Sapporo', 'Hokkaido Japan', 'School Of Tropical Medicine', 'Global Health', 'Nagasaki University', 'Nagasaki', 'Ateneo School Of Medicine']

Date: 2023-11

This study examined the effect of temperature on dengue in the 61 Provinces of the Philippines. We found that as temperature increases, so does the risk of dengue. However, the increase in the risk of dengue plateaus at 27.5 degrees Celsius, and beyond that temperature the risk decreases. The association between temperature and dengue incidence at high temperatures was stronger among provinces with high population density, less people in urban areas with larger household size, high poverty incidence, higher health spending per capita, and in lower latitudes compared to their counterparts. In particular, results showed that people in the rural areas had higher dengue risk than those in urban areas, which is contrary to expectations that dengue is an urban-dwelling related disease. On another hand, we found that high health spending per capita does not equate to better health outcomes. Efficiency of use of resources may come into play in terms of improving the health outcomes and not just the quantity or amount of the difference in the dengue risk by socio-economic factors highlights the need for a tailor-fit dengue response.

Data Availability: The data set for dengue incidence was acquired with permission from the Philippines Department of Health Epidemiological Bureau, while that of the provincial mean weekly temperatures were acquired from the Philippine Atmospheric, Geophysical and Astronomical Services Administration. The data sets may be accessed upon formal request by contacting the data source owners directly at [email protected] +63 (02) 8651 7800 and (+632)8927-1541 or [email protected] , respectively. R codes to perform the analyses are available upon formal request from the Department of Hygiene, Graduate School of Medicine, Hokkaido University at [email protected] .

Several temperature-dengue studies have documented evident heterogeneity in the association, which may potentially be due to socio-demographic and economic factors innate to the specific locations [ 16 – 19 ]. Certain vulnerable sub-populations have been observed to have higher infectious disease risk incidence such as those with higher poverty incidence, urban to population ratios, and smaller government health expenditures, as they have lower adaptive capacity and are more susceptible to negative health outcomes when exposed to environmental hazards [ 20 – 22 ]. While noted, examination of the contributors to heterogeneity of the linear association between temperature and dengue have been relatively scarce, with little to no evidence on how these socio-demographic and economic factors impact the non-linear temperature-dengue association. Against this backdrop, this study sought to determine the association between temperature and dengue in multiple locations in the Philippines, focusing on the role of socio-demographic and economic factors in the temperature-dengue association.

There have been similar climate change-dengue studies conducted in Southeast Asia, observing an increase of dengue risk within a specified range of high temperatures. In Thailand, vector efficiency increased in temperatures from 20–30 degrees Celsius, affecting annual cyclic patterns of dengue hemorrhagic fever epidemics in Bangkok [ 10 ]. Similar results were noted with the studies in Singapore [ 11 – 13 ]. A study in China [ 14 ] observed an inverse-U shape association between temperature and dengue. One study in the Philippines observed significant interaction effects between seasons, temperature, and altitude on the hatch rates (HRs) of A. aegypti larvae and reproductive outputs (ROs) of mosquitoes [ 15 ]. HRs in both seasons were highest at 25 degrees Celsius and lowest at 38 degrees Celsius. ROs were highest at 25 degrees Celsius in the wet season and at 18 degrees Celsius in the dry season. Furthermore, the study noted that the incidence of the dengue-causing mosquito phenotype was found to be highest at 18 degrees Celsius and present even at 38 degrees Celsius in both seasons, suggesting larval adaptation to global warming [ 15 ].

Literature identifies ambient temperature as an important meteorological risk factor for dengue transmission [ 5 ]. A systematic review by Li et. al. [ 5 ] revealed that the minimum temperature of 18 degrees Celsius is a critical limiting risk factor for dengue transmission, with high ambient temperatures promoting vector proliferation. Other global meta-analyses pose that dengue epidemics have a critical range of 22–29 degrees Celsius and that extremely higher temperatures reduce vector survival time [ 1 , 6 ]. In the advent of climate change, several studies have shown that current projections of 2 degrees Celsius by the end of the next century will result in an expansion of latitude and altitude range and duration of transmission of dengue [ 5 , 7 , 8 ] or even a shift in geographical regions affected by vector-borne disease outbreaks [ 9 ].

Dengue, transmitted by the mosquito vector Aedes aegypti, remains a major public health concern in tropical areas [ 1 ]. Over the last 50 years, dengue incidence has increased globally by 30-fold and 70% of the global burden of disease comes from Asia [ 2 ]. In 2019, the Philippines declared dengue a national epidemic of concern with 420,453 cases. Similar increases in dengue morbidity have been noted in Vietnam and Malaysia [ 3 ]. Heightened dengue transmission may be partly driven by the increase in international travel, unplanned urbanization, lack of effective vector management, climate change, and poor socio-economic status, among others [ 1 , 4 ].

Effect modification by socioeconomic factors was also examined through an interaction analysis. In brief, the 5 th and 95 th percentile values for the specific effect modifier represented the low and high percentiles. The respective values of the effect modifier were then deducted from the 5 th or the 95 th percentile values, which effectively centers these values around the low and high percentiles. The temperature crossbasis term is then multiplied by the centered effect modifier in either (low or high) percentiles. The resultant crossbasis represents the exposure-lag-response dependency corresponding to the percentile centering point of the effect modifier [ 38 ].

Y t,i is the weekly dengue incidence in time (t) of province (i). The weekly dengue incidence is assumed to follow a negative binomial distribution, which is similar to previous studies [ 35 ]. α is the intercept. logpop is the log of the population, which is an offset of the changing population sizes [ 5 ]. cbtemp is a crossbasis function of the exposure and lag dimensions of temperature, with both dimensions parameterized with 4 degrees of freedom (df) using natural cubic splines, respectively. The maximum lag was set at 18 weeks (or 4.5 months), based on the maximum lag sensitivity analyses ( S2 Fig ) and loosely similar to other studies whose maximum lag was set around 20 weeks [ 13 , 36 ]. After further examination, a seasonal pattern can be noted in Fig 1 , with a peak of the dengue cases around 32 nd /33 rd week. Thus, we adjusted for the seasonality of dengue cases through a 4 df, based on a sensitivity analysis in S3 Fig . Owing to the variations in the annual distribution of cases per year ( S4 Fig ), we added a random effect term of year; (1|Year). A random effect term of Province (1|Province) was also added to emulate a (partial) pooling effect in estimating the nationally representative association. Several studies have noted the potential immunity in the population, which could inflate the susceptible population in the subsequent time scale. To model this phenomenon, Imai et al. [ 37 ] utilized the previous time scale’s case as a potential adjustment. Here, we adjusted for the previous week’s case; Y t−1,i . ε is the error term.

In the absence of consistent background monitoring data for each of the provinces, hourly temperature was obtained from ERA5-land, which has a spatial resolution of 9km on a reduced Gaussian grid. Developed by the Copernicus Climate Change Service (C3S) on behalf of the European Union, ERA5-land is a publicly available enhanced global dataset that has shown accurate temperature readings in Europe and Central Asia [ 23 – 25 ]. Hourly temperatures were then aggregated to daily 24 hour mean temperature, and thereafter was aggregated into a weekly measure. Statistical correlation of temperature obtained from ERA5-land and those locations with existing background monitoring stations is high. For the case of Manila City, Pearson’s correlation coefficient was at 0.88 (p-val < 0.05), whereas for Cebu City it was at 0.86 (p-val < 0.05) ( S1 Fig ). Over the same period, daily temperature was aggregated by epidemiological week (EpiWeek) to match the dengue incidence data. On a similar note, several studies have utilized ERA5-land in examining the effect of temperature on specific health outcomes [ 26 – 28 ].

Nationwide, weekly notifiable dengue cases per province (n = 81) from 2010–2019 were obtained from the Philippines’ Department of Health–Epidemiological Bureau. In the Philippines, case surveillance is coursed through the Philippine Integrated Disease Surveillance and Response, which is subsequently managed by the Philippine Department of Health’s Epidemiological Bureau. In this study, provinces (n = 61) with more than 300,000 population were utilized for analysis. Provinces with less than 300,000 population had a higher percentage of 0 observations, affecting statistical power and the subsequent uncertainty in the estimation, and were thus excluded from the current study.

Socio-economic factors had varying impact on the temperature-dengue association (in Fig 5 ). Operationalization of the low and high for each socioeconomic variable is relative to each socioeconomic variable’s 5 th and 95 th percentiles, respectively. Provinces with high population density, less people in urban areas with larger household size, high poverty incidence, higher health spending per capita, and in lower latitudes were noted to exhibit statistically higher dengue risk compared to their counterparts at the upper temperature range (Figs 5 and S7 ). Modifications in the lower temperature range were also observed for all effect modifiers, except for latitude. Main effects of the effect modifiers are shown in S3 Table .

Fig 3 depicts the overall exposure-response association. Dengue-related risk is linearly increasing from the minimum temperature of 15.8 degrees Celsius and peaking at 27.5 degrees Celsius before declining. The inverted U-shaped risk pattern between temperature and dengue remains to be robust even after adjusting for precipitation as shown in S5 Fig . On the other hand, cumulative lag-response association at the 90 th percentile, relative to the 50 th percentile of the prediction interval, shows a reduced risk at the earlier lags but with a gradual increase in latter lags peaking at around 16 weeks (or approximately 4 months), followed by a subsequent reduction in the risks ( Fig 4 ). Further examination of the exposure- and lag-response associations by specific lag period (Weeks 0, 3, 6, 9, 12 and 18) and temperature percentiles (1 st , 5 th , 10 th , 90 th , 95 th , and 99 th ), revealed nearly consistent patterns with the overall exposure and lag associations (in S6 Fig ).

Fig 2A shows the distribution of dengue incidence per 100,000 population from 2010–2019, in all study provinces across the three major island groups of Luzon (north), Visayas (middle), and Mindanao (south). The top three provinces recording the highest dengue incidence were the highly urbanized provinces of Davao del Sur (11.14 per 100,000 population), South Cotabato (8.19 per 100,000 population), and Misamis Oriental (6.61 per 100,000 population). Provinces reporting the lowest dengue incidence were Lanao del Sur (0.29 per 100,000 population), Tawi-tawi (0.26 per 100,000 population, and Sulu (0.10 per 100,000 population), which are high-conflict areas in Mindanao [ 44 , 45 ]. Fig 2B indicates the average reported temperature per province in the same time period. The provinces that reported the highest temperatures between 26–27 degrees Celsius clustered around central and southern Luzon, Visayas, and the northern areas of Mindanao. The highest recorded daily mean temperatures were in Tawi-tawi (27.04 degrees Celsius), Sulu (26.88 degrees Celsius), and Basilan (26.67 degrees Celsius).

Discussion

Nationally representative exposure-response function between temperature and dengue in the Philippines was noted to have an inverse U-shaped, non-linear association, with the relative risk peaking at 27.5 degrees Celsius. In Guangzhou, China, both Wu et al. [14] and Xiang et al. [46], also noted an inverse U-shaped temperature-dengue association. The decline in the risk beyond 27.5 degrees Celsius can be partially explained by the physiological mechanisms affecting vector survival. Extremely hot temperatures beyond maximum temperatures have an inverse effect on fecundity, blood size meal, and vector viability, reducing mosquito survival time [6,47]. Dengue fever risk was seen to increase after extremely high temperatures due to accelerated incubation and increased mosquito population, although dengue lifespan and egg survival rate was decreased [48]. Hawley [49] observed that the optimal thermal conditions for A. albopictus survival is around 20 to 27 degrees Celsius. Another study by Servadio et. al. [9,13], presented similar results of a similar parabolic curve, with maximum development rate in vector-borne diseases at temperatures between 28 and 32°C and none under 18C or over 34°C in South and Southeast Asian countries. This is further supported by other studies [50,51], which observed a decline in vector’s survival rate beyond 28 degrees Celsius.

Our results have shown a 16-week (or approximately 4 months) lag association between temperature and dengue incidence as shown in Fig 4. Previous studies have also utilized longer lag periods to examine the lag association. In Singapore, Hii et al. [13] utilized a 20-week maximum lag to examine the effects of weather variables (including temperature) on dengue incidence. Results showed higher relative risks of dengue cases at a time lag of 3–4 months, which is consistent with our study findings. Similarly, Hii et al. [36] also utilized a 20-week maximum lag in examining the optimal lead time for dengue forecasting (parameterized by weather variables), with an optimal lead time at around 3 months. Koh et al. [52] also found that dengue incidence in 2005 was highly associated with weekly mean temperature at a lag time of 18 weeks. Heng et al. [53] noted that dengue incidence occurred at a lag time of 8–20 weeks following elevated temperature. Similar 3-month lag associations have been observed in Bangkok [54] and Taiwan [14]. The delayed association between temperature and dengue incidence may potentially be related to the incubation periods in host–vector–pathogen transmissions cycle plus reproduction, maturation and survival rate of vector [13].

Similar to previous studies, higher population density was associated with a visually higher temperature-dengueran association than low population dense areas [33,55]. In Indonesia, Nuranisa, et al. [56] noted a positive correlational association between dengue incidence and population density. These highly dense areas provide suitable conditions for vector proliferation and subsequent disease transmission [57]. Noteworthily, urbanization revealed a reduced risk for temperature-related dengue incidence, contrary to prior literature, which associates them as well-documented risk factors [17,19,54,58–61]. Unplanned urbanization promotes close proximity to larval habitats, dengue-prone water storage practices, and rapid population movement [1,59]. Notably, the top three provinces for dengue incidence, Davao del Sur, South Cotabato, and Misamis Oriental, are some of the more urbanized cities and more populated areas in the country [45]. These results suggest that the connection between dengue, human behavior, and temporal relationships is more complex than previously thought. As to whether dengue remains primarily a rural disease as a result of increased rural transmission or due to diminishing incidence rates in urban areas due to concentrated vector control efforts and increased health spending in urban areas, remains to be elucidated. Dengue transmission may indeed be increasing in rural areas due to the interplay of additional modifying variables that were not included in the study due to lack of geographical and behavioral data, such as forest cover or the presence of vegetation, open water containers, and degree of intrahousehold crowding, which are known dengue risk factors [1,17,19,54,58,59,61,62]. One study in Vietnam [34] found that risk of dengue was higher in rural areas due to the lack of piped water supply. Another study in Rio De Janeiro noted that while highly-populated, highly-urbanized areas displayed increased dengue incidence, so did other areas with low human density but with medium to highly-vegetated neighborhoods [62]. Further investigations are required to prove these hypotheses. While our analysis provides hindsight on the role of province-specific characteristics and its influence on the relationship between dengue and temperature, we cannot ascertain the possible mechanisms that are related to the observed associations.

Average household size was found to have higher temperature-related dengue risk than small households. Similar results were seen in New Caledonia [16]. Its known vectors, Aedes aegypti or Aedes albopictus, exhibit limited flight capacity and remain in close proximity to their breeding water containers [60,63]. Domestic workers and housewives are also more vulnerable to dengue compared to those who spend less time at home [17,59]. Larger families may also have larger dwellings, which may have more water-storage containers, such as flowerpots, vases, and lush areas that can serve as breeding sites [60]. As dengue cases have also been observed to cluster within households and other countries [63–65], intrahousehold dengue may be more prevalent than interhousehold dengue transmission. In contrast, a study in Brazil found large household sizes were protective against dengue risk, as smaller families (<7 people) were found to be living in more cramped spaces than their larger counterparts in the studied area [17]. This is in concurrence to the results observed by Zambrano, et. al. [59] who argued that the increased ratio of household members to the number of rooms, otherwise known as “crowding”, was more predictive of dengue incidence than the absolute number of household members. Risks apparent in the lower temperature range among large household size families may potentially also be related to the persistence of disease transmission whose range falls between 21.6–32.9 degrees Celsius [46] coupled by the crowding of the household members.

Poverty incidence was also found to modify the temperature-dengue association in both low and high temperature ranges. Several studies associate poverty with dengue risk due to substandard housing, solid waste mismanagement, and sewerage systems that favor close proximity to larval habitats [17,58,59]. Mulligan et. al. [61] argues that dengue’s epidemiological and transmission characteristics in relation to poverty should be interpretated in the context of how these specific poverty-related factors (i.e. household overcrowding, absence of air conditioning, low education, low income, poor housing quality, self-assessment of economic situation as poor, low family income, higher population density, and inappropriate use of covers on water containers, among others) affect the biology and ecology of the vector Aedes aegypti. In addition, poverty incidence may simply represent a variety of behavioral and exposure differences of the population. Resource-poor households/populations tend to work mostly during extreme conditions, including that during low temperatures and high temperatures in contrast to those resource-rich households/populations. Populations in resource-poor locations would prioritize more food and other basic necessities than personal risk prevention of dengue (e.g. mosquito repelling lotion, long sleeves, etc.). Since dengue transmission may persist in a range of 21.6–32.9 degrees Celsius [46]. If we chart against Fig 4 (poverty incidence panel), we note that the range 21 to 29 degrees Celsius, which is well within the range of dengue transmission, hence the potential high risk at these extremes for the high poverty incidence populations. The low risk at the middle may also be similarly related to low exposure of the population since these are temperatures that are conducive for most less labor-intensive jobs, whereby these populations don’t engage with. These are mostly assumptions that remain to be elucidated in future studies.

Contrary to the expected outcome that provinces with higher health expenditures would exhibit a reduced risk of temperature-related dengue risk due to the potential focus of these resources on overall health enhancement, our findings indicate a higher temperature-related dengue risk in areas with greater health spending as opposed to those with lower health expenditures The increased dengue risk amidst high health spending may seem counterintuitive, however, health spending does not necessarily equate to improved health. Dieleman et. al. [66] conducted a study analyzing factors associated with increases in US health care spending from 1996 to 2013. The study found that despite the substantial increase in health care spending during this period, there was no significant improvement in health outcomes, such as mortality rates or disease burden. The linkage between health care spending and improvement of health outcome is not direct and may be mediated by several factors. In Sub-Saharan Africa (SSA), Makuta and O’Hare [67] examined the relationship between public spending on health, quality of governance, and health outcomes. Their findings indicated that health spending has a statistically significant impact on improving health outcomes, such as under-five mortality and life expectancy. However, the impact of health spending on health outcomes is mediated by the quality of governance. Grigoli and Kapsoli [68] also emphasized the need to consider the efficiency of health spending. Simply increasing public expenditure in the health sector may not significantly affect health outcomes if the efficiency of spending is low. Therefore, it is crucial to ensure that health spending is allocated and utilized efficiently to achieve desired improvements in health outcomes. Similarly, the increased in temperature-related dengue risk may also be due to the higher health-seeking behavior of the population, especially in resource-rich locations. Consistent with a previous literature [69], hydrometeorological changes were found to substantially affect the risk of dengue in regions with mid-to-low latitudes in China. Patz et. al. [8] assert that increased incidence of dengue may first occur in regions bordering endemic zones in latitude or altitude. In summary, from a programmatic, local risk standpoint, our study shows that the association can be modified by area specific factors, thus merits the need for a tailor-fit dengue-control response. The Philippines may benefit from devolved decentralized and localized dengue control and early warning systems in order to account for these variabilities.

[END]
---
[1] Url: https://journals.plos.org/plosntds/article?id=10.1371/journal.pntd.0011700

Published and (C) by PLOS One
Content appears here under this condition or license: Creative Commons - Attribution BY 4.0.

via Magical.Fish Gopher News Feeds:
gopher://magical.fish/1/feeds/news/plosone/

You are viewing proxied material from magical.fish. The copyright of proxied material belongs to its original authors. Any comments or complaints in relation to proxied material should be directed to the original authors of the content concerned. Please see the disclaimer for more details.