|Year : 2020 | Volume
| Issue : 1 | Page : 79-86
Exploration of population ecological factors related to the spatial heterogeneity of dengue fever cases diagnosed through a national network of laboratories in India, 2017
Vasna Joshua1, Kanagasabai Kaliaperumal1, Kirubakaran Bajji Krishnamurthy1, Ravi Muthusamy1, Ramachandran Venkatachalam1, Kethara A Gowri1, Vishal C Shete1, Sabarinathan Ramasamy1, Nivedita Gupta2, Manoj Vasant Murhekar3, for VRDL Team4
1 VRDLN Data Management Group, Chennai, Tamil Nadu, India
2 Division of Epidemiology & Communicate Diseases, Indian Council of Medical Research, New Delhi, India
3 ICMR-National Institute of Epidemiology, Chennai, Tamil Nadu, India
|Date of Submission||12-Jun-2018|
|Date of Web Publication||24-Feb-2020|
Dr Vasna Joshua
Indian Council of Medical Research-National Institute of Epidemiology, Second Main Road, Tamil Nadu Housing Board, Ayapakkam, Near Ambattur, Chennai 600 077, Tamil Nadu
Source of Support: None, Conflict of Interest: None
| Abstract|| |
Background & objectives: Dengue virus (DENV) transmission is known to be influenced by the environmental conditions. During 2017, the Viral Research and Diagnostic Laboratories (VRDLs) tested 78,744 suspected dengue fever (DF) patients, of whom, 21,260 were laboratory confirmed. The objectives of the study were to evaluate the hypothesis that spatial heterogeneity existed for DF patients and to identify significant determinants of DENV transmission in various districts across the Indian States during 2017.
Methods: Laboratory confirmed DF cases were analysed from 402 districts spread across the Indian States. The determinants for DF transmission included in the model were population density, proportion of population living in rural areas, proportion o f forest cover area to the total geographical area, proportion of persons not able to read and write and who were aged greater than seven years; the climatic variables considered were minimum, maximum and average temperature, precipitation and cumulative rainfall. The spatial heterogeneity was assessed using spatial regression analysis.
Results: DF cases showed strong spatial dependency, with Moran's I=4.44 (P <0.001). The robust measure for spatial lag (6.55; P=0.01) was found to be the best model fit for the data set. Minimum temperature and cumulative rainfall were significant predictors.
Interpretation & conclusions: A significant increase in the number of dengue cases has occurred when the minimum temperature was 23.0-25.8°C and the cumulative rainfall 118.14-611.64 mm across the Indian districts. Further in-depth investigations incorporating more number of demographic, ecological and socio-economic factors would be needed for robust conclusions.
Keywords: Dengue - ecological factors - humidity - India - spatial regression - temperature - transmission
|How to cite this article:|
Joshua V, Kaliaperumal K, Krishnamurthy KB, Muthusamy R, Venkatachalam R, Gowri KA, Shete VC, Ramasamy S, Gupta N, Murhekar MV, for VRDL Team. Exploration of population ecological factors related to the spatial heterogeneity of dengue fever cases diagnosed through a national network of laboratories in India, 2017. Indian J Med Res 2020;151:79-86
|How to cite this URL:|
Joshua V, Kaliaperumal K, Krishnamurthy KB, Muthusamy R, Venkatachalam R, Gowri KA, Shete VC, Ramasamy S, Gupta N, Murhekar MV, for VRDL Team. Exploration of population ecological factors related to the spatial heterogeneity of dengue fever cases diagnosed through a national network of laboratories in India, 2017. Indian J Med Res [serial online] 2020 [cited 2020 Nov 30];151:79-86. Available from: https://www.ijmr.org.in/text.asp?2020/151/1/79/279129
VRDLs Team (arranged alphabetically by city): Majumdar T, Government Medical College, Agartala; Vegad MM, BJ Medical College, Ahmedabad; Anukumar B, ICMR-National Institute of Virology (Field Unit), Alappuzha; Shukla I, Jawaharlal Nehru Medical College, Aligarh; Dhingra K, Government Medical College, Amristar; Swarnalatha, Government Medical College, Anantpur; Shantala GB, Bengaluru Medical College and Research Institute, Bengaluru; Biswas D, All India Institute of Medical Sciences, Bhopal; Dwibedi B, ICMR-Regional Medical Research Centre, Bhubaneswar; Ratho RK, Postgraduate Institute of Medical Education and Research, Chandigarh; Kaveri K, King Institute of Preventive Medicine and Research, Chennai; Borkakoty BJ, ICMR-Regional Medical Research Centre for NE Region, Dibrugarh; Verma RK, Uttar Pradesh Rural Institute of Medical Sciences, Etawah; Sharma A, Gauhati Medical College, Guwahati; Venkatesha DT, Hassan Institution of Medical Sciences, Hassan; Reddy S, Osmania Medical College, Hyderabad; Singh L, Jawaharlal Nehru Institute of Medical Sciences Hospital, Imphal-East; Barde PV, ICMR-National Institute for Research in Tribal Health, Jabalpur; Majumdar D, Late Sri Baliram Kashyap Memorial Government Medical College, Jagdalpur; Malhotra B, SMS Medical College, Jaipur; Sharma SS, Government Medical College, Jammu; Sinha M, M.P. Shah Government Medical College, Jamnagar; Tiwari YK, Jhalawar Medical College, Jhalawar; Khatri PK, Dr. S.N. Medical College, Jodhpur; Barua P, Jorhat Medical College & Hospital, Jorhat; Jaryal SC, Dr Rajendra Prasad Government Medical College, Kangra; Srinivasan S, Indira Gandhi Medical College and Research Institute, Kathirkamam; Ray R, Institute of Post Graduate Medical Education and Research, Kolkata; Dutta S, ICMR-National Institute of Cholera and Enteric Disease (NICED), Kolkata; Beena PJ, Government Medical College, Kozikhode; Jain A, King George’s Medical University, Lucknow; Sugumari C, Madurai Medical College, Madurai; Arunkumar G, Manipal Centre for Virus Research, Manipal University, Manipal; Amrutha KB, Mysore Medical College and Research Institute, Mysuru; Raut S, Indira Gandhi Government Medical College, Nagpur; Bakshi R, Government Medical College, Patiala; Das P, ICMR-Rajendra Memorial Research Institute of Medical Sciences, Patna; Srivastava RK, Patna Medical College and Hospital, Patna; Vijayachari P, ICMR-Regional Medical Research Centre, Port Blair; Kumar M, Rajendra Institute of Medical Sciences, Ranchi; Gil PS, Pt. B.D. Sharma Post Graduate Institute of Medical Sciences, Rohtak; Phukan AC, North Eastern Indira Gandhi Regional Institute of Health and Medical Sciences, Shillong; Singh D, Indira Gandhi Medical College, Shimla; Kakru DK, Sher-I-Kashmir Institute of Medical Sciences, Srinagar; Mythreyee KM, Government Theni Medical College, Theni; Devi S, Government Medical College, Thiruvananthapuram; Pillai RM, Rajiv Gandhi Centre for Biotechnology, Thiruvananthapuram; Nair PK, Govt. Medical College, Thrissur; Kalawat U, Sri VenKetswara Institute of Medical Science, Tirupati; Kumari LR, Govt. Siddhartha Medical College, Vijayawada; Kamala P, Andhra Medical College, Vishakhapatnam
Dengue disease is an important public health problem in India. The incidence of dengue has increased 30-fold in the last 50 yr with increasing geographic expansion to new places and in the current decade, from urban to rural settings. In the early 20th century, dengue was endemic in a few southern (Maharashtra, Karnataka, Tamil Nadu and Puducherry) and northern States (Delhi, Rajasthan, Haryana, Punjab and Chandigarh). It has spread to many other States in India, including the Union Territories.
The transmission of dengue virus (DENV) is influenced by the environmental conditions; short-term changes in weather, particularly temperature, precipitation and humidity are correlated with dengue incidence. The intensity and pattern of occurrence of dengue are dependent on interrelated human, vector, environment and virus-related factors,,. These factors have been considered as drivers of the spatial patterns of transmission. Social factors such as educational status of the household members and demographics have been found to be negatively correlated with increased dengue burden. It has been reported that DENV evolution first started with Aedes albopictus in the forests of Southeast Asia, and later as Ae. aegyti and its transmission was more in forest cover areas,. Further, weather variables predicted the intensity and timing of outbreaks (OB), which included minimum, maximum and mean temperature; relative humidity; wind velocity and precipitation. Considering the rising incidence of DENV infection in India, this study was conducted to find out environmental factors associated with the DENV transmission. The main objectives of this study were to evaluate the hypothesis that spatial heterogeneity existed in distribution of dengue fever (DF) cases and to identify significant determinants of DF transmission in different districts in India. This study was undertaken at 51 Viral Research and Diagnostic Laboratories (VRDLs) in 26 States across the country established under the Department of Health Research and Indian Council of Medical Research (DHR/ICMR), Government of India, New Delhi, India.
| Material & Methods|| |
The study was conducted after obtaining the ethical clearance from the Institutional Ethics Committee of the ICMR-National Institute of Epidemiology, Chennai, India. During 2017, 51 VRDLs (41 medical college level, 5 State level and 5 regional level) were functional. Besides the 26 States where these VRDLs were located, these laboratories also provided diagnosis to suspected DF patients from five neighbouring States. Thus, 402 districts from 31 Indian States were included in developing the model.
Data source:During 2017, the Viral Research and Diagnostic Laboratory Network (VRDLN) investigated 78,744 suspected DF patients for NS1 antigen and/or IgM antibodies, of whom, 21,260 (27%) were found to be positive. It was also observed that the median time from the onset of symptom to the presentation was five days [Inter quartile range (IQR) of 3-7 days]. Information about time (date of onset of illness and date of seeking care), place (village, sub district/tehsil and district) and person (age and sex) collected from each suspected DF patient was extracted and used for developing the model.
Geocoding: Four hundred and two districts were converted into point features that contained location (latitude and longitude) and spatial geometry using ArcGIS software (Environmental Systems Research Institute, ArcGIS Desktop: Release 10, 2010). Since the data pertained to only 402 districts, using the point shape file and the software GeoDa version 1.8.14 software, the Thiessen polygons were formed by drawing boundaries according to the distribution of the available data points for 2017, with one polygon per district [Figure 1] covering the entire/maximum area of VRDLN data availability.
|Figure 1: Thiessen polygons for suspected dengue fever reported by the Virus Research and Diagnostic Laboratories in 402 districts across Indian States, 2017.|
Click here to view
Covariate selection was based on review of literature for socio-demographic risk factors associated with DF and essentially with the available data. The determinants included in the model were population density (POPDENS) (i.e., number of inhabitants/km); proportion of population who lived in rural areas (RURAL), proportion of forest cover area to the total geographical area (FCA); proportion of persons who were not able to read and write and who were aged greater than seven years (ILLIT). The data for the above determinants were obtained from 'India population 2017 portal'. The four climatic variables considered were minimum temperature (MIN), maximum temperature (MAX), average temperature (AVG) and precipitation (PRECIP). These were downloaded from AccuWeather website for the respective time point and districts. The cumulative rainfall (CRAIN) was obtained from India Meteorological Department website/local government website dashboards. Since dengue positivity was higher during the months of August to November, with 18,265 (79%) patients occurring during this period [Figure 2], average values for the four climatic variables for the months of August to November 2017 were considered for analysis. The summary statistics of the variables are shown in [Table 1].
|Figure 2: Number of patients tested for dengue virus per month and percentage of positive tests.|
Click here to view
Statistical analysis: The software tool GeoDa was used for exploratory spatial data analysis including data manipulation, mapping and spatial regression analysis. The spatial neighbours were obtained using the Thiessen polygons and rook's weight of 1st order. Using the software, ordinary least square (OLS) regression and spatial regression models were constructed considering district-wise proportion of dengue cases (DENG) as dependent variable and POPDEN, RURAL, ILLIT, FCA, MIN, MAX, AVG, CRAIN and PRECIP as independent variables.
As a first step, the OLS model and its fit under classical Best Linear Unbiased Estimator assumptions was examined (the random errors have a mean of zero, constant variance and follows normal distribution). The step-wise technique was used to select a minimal adequate model with the smallest Akaike's information criterion (AIC) value, with P=0.10 and 0.05 as the entry and removal criteria, respectively. The error variance could well be affected by the spatial dependence in the data set. The spatial independency of residuals was evaluated using spatial autocorrelation coefficient, Moran's I, which was expressed as:
Where n was the total number of districts in the study; i and j represented different districts; x i was the residual of i and x̄ was the mean of residuals; w ij was a measure of spatial weights of i and j.
The values of Moran's I would be approximately between +1 (positive autocorrelation) and −1 (negative autocorrelation), and the expected value in the absence of autocorrelation was (−1)/(n−1). Positive spatial autocorrelation meant similar values tended to occur in adjacent areas, while negative autocorrelation implied nearby locations tended to have dissimilar values. If no spatial autocorrelation was found, then the spatial arrangement would be completely at random.
The relationship between proportion of dengue cases (DENG) by district and set of determinants were explored using the spatial regression approach. Two distinct spatial regression models, i.e., spatial lag model and spatial error models were examined. The spatial lag model is appropriate when the value of an event in one region is directly influenced by the values of the same event in its neighbouring regions. Spatial error model is appropriate when the concern being to adjust for the bias occurring from spatial autocorrelation or spatial dependency. A diagnostic test for spatial lag and error models, namely, Lagrange's multiplier (LM) test for both the models was applied. Further a modified LM test (i.e., Robust LM test) for both spatial lag and error model was attempted if both the models were significant. Finally, the best model fit was with the higher values of R values (significant P value), higher Log likelihood and the lower AIC value. Spatial regression diagnostics were examined using Jarque-Bera test (a goodness of fit test to check for normality of errors), Breusch-Pagan test (test for heteroskedasticity which measures the normality of the error terms) and Lagrange Multiplier (a diagnostic test for spatial lag and error models).
The AIC is the measure of relative goodness of fit of a statistical model. In the general case, AIC=2k–2ln (L); where k denotes the number of parameters in the statistical model, and L, the maximized value of the likelihood function for the estimated model. Given a set of models for the data, the preferred model is the one with the minimum AIC value.
| Results|| |
Of the nine determinants, only four determinants viz., minimum temperature, maximum temperature, precipitation and cumulative rainfall were found suitable for the model fit. The OLS model explained 53 per cent of variation in the dataset [Table 2]. It was also observed that the correlation coefficients between the proportion of dengue cases within district (DENG) and with the four determinants were significantly high (P <0.05).
|Table 2: Summary of ordinary least square, spatial lag and spatial error model|
Click here to view
Other diagnostics of the model were also examined. The Jarque-Bera test score (2.763; P=0.251) indicated normality of the error term and the low probability of the Breusch-Pagan test (27.47; P=0.053) and non-significant points showed the non-existence of heteroskedasticity. A significant Moran's I value of 4.44 (P <0.001) indicated a strong positive spatial autocorrelation of the residuals [Table 2]. It was also observed that the OLS model overestimated the coefficients; hence, the spatial lag and error models were considered better for capturing the spatial effect in the model.
The LM test for both the models showed [Table 2] that the test statistics was significant for both LM Spatial lag model (20.86; <0.001) and LM Spatial error model (17.20; <0.001). Both the models explained 68 per cent of the variation in the data set with an increased percentage of 15 per cent variation than the OLS model. Coefficients of the terms capturing spatial effects, viz., P in spatial lag model (P=029; P <0.001) and λ in spatial error model (λ=0.30; P <0.001)were also significant. The positive value of λ and P indicated substantial spatial dependence in dengue cases across the neighbouring districts. The Robust LM of spatial lag model (6.55; P=0.01) was significant, whereas the Robust LM error model (2.89; P=0.09) was insignificant, which meant that in the presence of error-dependent variable, the spatial lag dependency disappeared. As a result, the general model fit improved for Robust LM spatial lag model, as indicated with the higher values of R2(0.68), significant P=0.01, higher Log likelihood value of 28.5 and lower AIC value of −47.0. Hence, the spatial lag model was considered to be the best model fit rather than OLS and spatial error model. The determinants, minimum temperature and cumulative rainfall turned out to be significant predictors for the spatial lag model.
These analyses indicated that DF cases had a strong spatial clustering and were significantly higher in districts with higher minimum temperature and higher cumulative rainfall across the various districts studied.
| Discussion|| |
Our analyses indicated that DF cases diagnosed through the VRDLs were spatially clustered across the districts in India. The significant Moran's I value showed that dengue cases occurred simultaneously in the same district or adjacent districts in India during 2017. This may also be due to the fact that spatial clustering of any disease is inevitable since human population generally live in spatial clusters rather than random distribution in space.
Exploring the 21,260 serologically positive dengue cases in the study area and their location ecological factors, it was observed from the data (data not shown) that significantly increased number of cases occurred with a minimum temperature ranging between 23.0 and 25.8°C (χ2=11.92; P <0.001) and the cumulative rainfall of 118.14-611.64 mm (χ2=10.09; P <0.001) in the districts.
Many studies have examined the relationship between climatic factors and dengue at different geographic locations. High rainfall has been reported to be associated with dengue cases in Mexico, Puerto Rico, Taiwan, Barbados, Indonesia, Thailand, Trinidad and Venezuela. Another study has shown that temperatures in the range of 22 to 31°C has provided a suitable environment for breeding and abundance of Aedes mosquito species and thereby increasing the risk of dengue cases. A study from Lahore showed that minimum temperature had a significant positive effect whereas maximum temperature and wind showed a significant negative effect. In a study conducted in Bhopal, India most dengue cases occurred in the period followed by peak rainfall, when mean minimum temperature had started falling, while mean maximum temperatures were still high.
The present study has supported the fact that high minimum temperature and increased cumulative rainfall have been conducive for the propagation of dengue virus transmission and led to a significantly increased dengue cases across the various districts in India. The maximum temperature in the present study reached to a maximum of 36.5°C. It has been shown that presumably high temperatures that exceed the optimum may lead to a lower transmission probability of dengue virus than expected under a constant temperature model. The precipitation in the study area ranged between 11 and 391 mm in a month. Heavy precipitation events may sometimes wash away the breeding sites of dengue. These might be the reasons for non-significance of the dengue cases in the study area.
Our analysis had certain limitations. First, the analysis was based on the data collected from patients seeking care at the medical college hospitals and reported by the health authorities covered under VRDLN and might not include all the cases of DF spread over the entire country. Places where VRDLs are located are likely to have higher proportion of patients than the neighbouring districts or districts within neighbouring States. Second, relative humidity was found to be strongly associated with mosquito proliferation in several studies,. However, data about relative humidity were available for only 82 districts for 2017. In these 82 districts, relative humidity and rainfall were highly correlated (correlation coefficient=0.71, P <0.05). Therefore rainfall data were included in the model. The data analysis was based on climate change variables for the four high positivity dengue months (August to December 2017) for the 402 districts; hence it might not represent the entire country. The patient's information was geocoded using the address reported by the patient during the time of investigation and there might be a little shift from their exact location of residence. The study facts pertained to the ecological conditions for the year 2017. Further in-depth investigation by including more number of demographic, ecological and socio-economic factors in the analysis and for a longer period would be more precise for robust conclusions.
Acknowledgment: Authors thank the principal investigators, scientists and staff from Resource Centre, Virus Diagnostic Laboratory (VDL), ICMR-National Institute of Virology, Pune. Authors also thank all consultants of Virus Research and Diagnostic Laboratory Network.
Financial support & sponsorship: Authors acknowledge the Indian Council of Medical Research, New Delhi, for financial support.
Conflict of Interest: None.
| References|| |
Chakravarti A, Arora R, Luxemburger C. Fifty years of dengue in India. Trans R Soc Trop Med Hyg
Scott TW, Morrison AC. Aedes aegypti
density and the risk of dengue virus transmission. In: Takken W, Scott TW, editors. Ecological aspects for application of genetically modified mosquitoes
. The Netherlands: FRONTIS, Dordrecht; 2003. p. 187-206.
Kuno G. Review of the factors modulating dengue transmission. Epidemiol Rev
Zellweger RM, Cano J, Mangeas M, Taglioni F, Mercier A, Despinoy M, et al.
Socioeconomic and environmental determinants of dengue transmission in an urban setting: An ecological study in Nouméa, New Caledonia. PLoS Negl Trop Dis
Ebi KL, Nealon J. Dengue in a changing climate. Environ Res
Lippi CA, Stewart-Ibarra AM, Muñoz ÁG, Borbor-Cordova MJ, Mejía R, Rivero K, et al
. The social and spatial ecology of dengue presence and burden during an outbreak in Guayaquil, Ecuador, 2012. Int J Environ Res Public Health
Rudnick A. Studies of the ecology of dengue in Malaysia: A preliminary report. J Med Entomol
Morin CW, Comrie AC, Ernst K. Climate and dengue transmission: Evidence and implications. Environ Health Perspect
Environmental Systems Research Institute, Inc. ArcGIS desktop
: Release 10
. Redlands, CA, USA: ESRI, Inc.; 2010.
Syabri AL, Kho Y. GeoDa: An introduction to spatial data analysis. Geogr Anal
India Meteorological Department, Ministry of Earth Sciences. Government of India; 2015. Available from: https://mausam.imd.gov.in/
, accessed on February 2, 2018.
Chakravarti IM, Laha RG, Roy J. Handbook of methods of applied statistics
, vol. 1. Hoboken, NJ, USA: John Wiley & Sons; 1967.
Moran PAP. The interpretation of statistical maps. J Roy Stat Soc
Lawson AB, Williams FLR. An introductory guide to disease mapping
. Hoboken, NJ, USA: John Wiley & Son; 2001. p. 1-75.
Chowell G, Sanchez F. Climate-based descriptive models of dengue fever: the 2002 epidemic in Colima, Mexico. J Environ Health
: 40-4, 55.
Johansson MA, Dominici F, Glass GE. Local and global effects of climate on dengue transmission in Puerto Rico. PLoS Negl Trop Dis
Chen SC, Liao CM, Chio CP, Chou HH, You SH, Cheng YH. Lagged temperature effect with mosquito transmission potential explains dengue variability in southern Taiwan: insights from a statistical analysis. Sci Total Environ
Depradine C, Lovell E. Climatological variables and the incidence of Dengue fever in Barbados. Int J Environ Health Res
Corwin AL, Larasati RP, Bangs MJ, Wuryadi S, Arjoso S, Sukri N, et al
. Epidemic dengue transmission in southern Sumatra, Indonesia. Trans R Soc Trop Med Hyg
Wiwanitkit V. An observation on correlation between rainfall and the prevalence of clinical cases of dengue in Thailand. J Vector Borne Dis
Chadee DD, Shivnauth B, Rawlins SC, Chen AA. Climate, mosquito indices and the epidemiology of dengue fever in Trinidad (2002-2004). Ann Trop Med Parasitol
Cecilia D. Current status of dengue and chikungunya in India. WHO South East Asia J Public Health
Ahmed SA, Siddiqi JS, Quaiser S, Kamal S. Using PCA poisson and negative binomial model to study the climatic factor and dengue fever outbreak in Lahore. J Basic Appl Sci
Pakhare A, Sabde Y, Joshi A, Jain R, Kokane A, Joshi R. A study of spatial and meteorological determinants of dengue outbreak in Bhopal City in 2014. J Vector Borne Dis
Lambrechts L, Paaijmans KP, Fansiri T, Carrington LB, Kramer LD, Thomas MB, et al
. Impact of daily temperature fluctuations on dengue virus transmission by Aedes aegypti
. Proc Natl Acad Sci U S A
Viana DV, Ignotti E. The ocurrence of dengue and weather changes in Brazil: a systematic review. Rev Bras Epidemiol
Xu HY, Fu X, Lee LK, Ma S, Goh KT, Wong J, et al
. Statistical modeling reveals the effect of absolute humidity on dengue in Singapore. PLoS Negl Trop Dis
[Figure 1], [Figure 2]
[Table 1], [Table 2]