Spatial cluster analysis of human cases of Crimean Congo hemorrhagic fever reported in Pakistan
© Abbas et al.; licensee BioMed Central. 2015
Received: 30 July 2014
Accepted: 15 January 2015
Published: 2 March 2015
Crimean Congo hemorrhagic fever (CCHF) is a tick-borne viral zoonotic disease that has been reported in almost all geographic regions in Pakistan. The aim of this study was to identify spatial clusters of human cases of CCHF reported in country.
Kulldorff’s spatial scan statisitc, Anselin’s Local Moran’s I and Getis Ord Gi* tests were applied on data (i.e. number of laboratory confirmed cases reported from each district during year 2013).
The analyses revealed a large multi-district cluster of high CCHF incidence in the uplands of Balochistan province near it border with Afghanistan. The cluster comprised the following districts: Qilla Abdullah; Qilla Saifullah; Loralai, Quetta, Sibi, Chagai, and Mastung. Another cluster was detected in Punjab and included Rawalpindi district and a part of Islamabad.
We provide empirical evidence of spatial clustering of human CCHF cases in the country. The districts in the clusters should be given priority in surveillance, control programs, and further research.
KeywordsCrimean Congo hemorrhagic fever Pakistan Spatial autocorrelation Clusters Adjusted cumulative incidence
Please see Additional file 1 for translations of the abstract into the six official working languages of the United Nations.
Crimean-Congo hemorrhagic fever (CCHF) is an arboviral zoonotic infection with potential for human-to-human transmission. The disease is widely distributed across Africa, Asia, and Europe. It is generally asymptomatic in infected animals (usually domestic livestock), but highly fatal in humans with a 10 to 50% case fatality . The infection is transmitted to humans by infected ticks, direct contact with meat or blood of infected animals, or direct contact with the blood or secretions of an infected person . The causative virus of CCHF is a member of the genus Nairovirus in the family Bunyaviridae. It is an enveloped virus containing negative-sense; single-stranded RNA composed of three segments. The virus has been classified into seven genotypes (Asia-1, Asia-2, Africa-1, Africa-2, Africa-3, Euro-1, and Euro-2) that coincide well with their geographical regions .
The first confirmed case of CCHF in Pakistan was reported in 1976 followed by outbreaks from almost all areas of the country mainly Balochistan province where the disease still remains entrenched. Genetic analysis of the viruses has confirmed the presence of Asia-1 and Asia-2 genotypes in this province. In Pakistan, CCHF exhibits two annual peaks from March to October . Since 2000, there is an increasing trend in the number of cases, for instance during 2012–2013, a total of 162 cases including 38 deaths reported in different areas . The disease incidence continues to expand and cover previously uninfected areas across all four provinces. Locating areas with relatively high incidence is essential in order to target surveillance and prevention programmes. The purpose of this study was to explore data about human CCHF cases for presence of spatial clusters of the disease in Pakistan.
Information about the occurrence of the disease was extracted from the Weekly Epidemiological Bulletin(s) jointly published by the National Institute of Health (NIH), Islamabad and the country office of the World Health Organization (WHO), Pakistan. The dataset included total count of laboratory confirmed human cases of CCHF reported from different country areas during January to December, 2013. SaTScan™ v9.3.1 and GeoDa™ v.1.6.6 software were used to execute the analyses. Population of each district was extrapolated from the linear growth rate between the censuses of 1981 and 1998. Incidence per million of population was calculated for each district, with the total number of cases as the event variable and the estimated population for 2013 as the base variable. The spatial empirical Bayes (EB) smoothing method was applied to adjust incidence proportion estimates using first order queen contiguity weights. The spatial EB method adjusts rates toward the average of the neighbors rather than the overall mean of the study area. There were no reliable population data for northern areas and disputed Kashmir; therefore these regions were excluded from the analyses.
Global Moran’s I statistic was used to explore data for evidence of spatial clustering at the country level. The results of the test vary between + 1.0 and -1.0. The Moran’s I > 0, = 0, and < 0 indicate the positive spatial autocorrelation, random distribution, and negative spatial autocorrelation, respectively. In order to identify the location of the clusters, spatial scan statistic and local indicators of spatial autocorrelation, local Moran’s I and Getis ord Gi* were also used. The two techniques have been used simultaneously to complement the findings. The number of the permutation test was set to 999 and the significance level was set to 0.05. The formulas of the techniques are given in the Appendix.
where C is the total number of cases, c is the observed number of cases within the window, E [c] is the expected number of cases within the window under the null hypothesis (assuming constant risk in the study domain), and I() is the indicator function which is equal to 1 when the window is more than expected . A retrospective purely spatial scan statistic was applied to detect high -risk clusters of CCHF with a discrete Poisson model. The maximum spatial cluster size of 50% of the population at risk was applied to avoid pre-selection bias.
Studies conducted in Iran and Turkey have also found spatial clustering of CCHF [6, 8]. Khushal et al. studied the residential location of CCHF cases admitted to three public sector hospitals of Quetta –the capital city of Balochistan. The districts with a relatively high prevalence were located in cluster detected in the province. Alam et al.  published a map on the distribution of IgM confirmed CCHF cases reported to National Institute of Health, Islamabad (2003 to 2008). All affected districts were located in clusters identified in our study.
The reasons for detection of disease cluster in Balochistan are unknown. The disease is endemic in this province since 1987. Many cases from Afghanistan are referred to tertiary care hospitals in Quetta city . The province lies on migration routes of the nomads . Moreover, there is cross-border trade of animal’s skins and hides through Balochistan . More than 60% of all human cases of CCHF in Iran have been reported in Sistan and Baluchistan Province. This has been attributed to excessive livestock import and high frequency of Afghani traveling to these areas . In Pakistan, livestock trade and movement during eid ul azha has potential to spread infection in the country. The animals are brought for sale in cities such as Rawalpindi . Relationship between altitude and CCHF occurrence has been investigated. The altitude has been supposed to affect tick distribution and animal husbandry practices . Estrada-Peña et al. modeled the spatial distribution of crimean-congo hemorrhagic fever outbreaks in Turkey and found that areas with higher risk of CCHF were correlated with the zones of high climate suitability for the tick together with a high rate of fragmentation of agricultural land interspersed between forest and shrub-type vegetation. Yagci-Caglayik et al. studied seroprevalence and risk factors of Crimean-Congo hemorrhagic fever in selected seven provinces of Turkey. The most important risk factors for CCHF seropositivity were older age, male gender, illiterate, farmer, animal husbandry, living in rural residence in adobe houses, and a previous tick bite history.
It is important to note that each of the techniques used in this study has its own advantages and limitations. The spatial scan statistic has several advantages over other clustering detection methods because it: i) accounts for problem of multiple testing, ii) adjusts for the heterogeneous population and any number of categorical covariates, and iii) searches for clusters without specifying their size or location thereby overcoming the problem of pre-selection bias. However, SaTScan, often reports heterogeneous clusters i.e. a large cluster may contain a number of low- risk areas. Moreover, it is sensitive to user-defined settings . The shapes of clusters may not always be circular or elliptical. The Getis Ord Gi* identifies clusters of high or low values but cannot capture negative spatial autocorrelation (spatial outliers). Local Moran’s I statistic, on the other hand, has the ability to detect both positive and negative spatial correlations. The Gi* statistic has been preferred in some studies as it matches the usual definition of cluster and is indicated for use with variables that posse natural origin . Geary’s C is another measure of spatial autocorrelation. It emphasizes on differences in values between pairs of observations. sensitive to differences in small neighborhoods and statistically less well behaved compared to local Moran’s I. Therefore Moran’s I is usually preferred to Geary’s C.
Limitations of this study included the use of secondary data and the problems inherent to data quality. Spatial health data might have uncertainty in location (i.e. place of residence and location where the attributes recorded were not same), underreporting and mis-diagnosis. The level of those uncertainties in disease data is unknown and cannot be controlled. The disease patterns are scale dependent and subjected to temporal changes. These facts should be considered in interpretation of the findings. We extrapolated population from inter census growth rate from 1981 to 1998. Population count data are generally stable so it was quite acceptable to extrapolate.
In short, CCHF is an emerging zoonotic disease in Pakistan. Understanding the spatial patterns of the disease is important in order to make rational use of the limited resources. Our study provides empirical evidence for the existence of statistically significant spatial clusters of high CCHF incidence in few districts of Balochistan, Rawalpindi and Islamabad. Those districts may be targeted for research and healthcare interventions on priority such as public awareness programmes, one health initiative etc. Our findings emphasize to initiate sustained and comprehensive surveillance activities focusing timely detection and disease control strategies in these areas. The location of a cluster close to Afghanistan border highlights the importance of the findings for organizations dealing with the disease at the national, regional and global levels.
where x represents CCHF incidence within a given district, W ij is a spatial weight which defines neighboring districts j to I, W i is the sum of weights W ij , is the mean of the country level CCHF incidence, , s is the standard deviation of the x values.
The authors acknowledge NIH, Islamabad for providing the disease data used in this analysis.
- Appannanavar SB, Mishra B: An update on crimean congo hemorrhagic fever. J Global Infect Dis 2011, 3:285–92. 10.4103/0974-777X.83537View ArticleGoogle Scholar
- Sharifi-Mood B, Metanat M, Alavi-Naini R: Prevalence of crimean-congo hemorrhagic fever among high risk human groups. Int J High Risk Behav Addict 2014, 3:e11520.View ArticlePubMed CentralPubMedGoogle Scholar
- Alam MM, Khurshid A, Sharif S, Shaukat S, Suleman RM, Angez M, et al.: Crimean-Congo hemorrhagic fever Asia-2 genotype, Pakistan. Emerg Infect Dis 2013, 19:1017–9. 10.3201/eid1906.120771View ArticlePubMed CentralPubMedGoogle Scholar
- Sheikh AS, Sheikh AA, Sheikh NS, Rafi US, Asif M, Afridi F, et al.: Bi-annual surge of Crimean-Congo haemorrhagic fever (CCHF): a five-year experience. Int J Infect Dis 2005, 9:37–42. 10.1016/j.ijid.2004.02.007View ArticlePubMedGoogle Scholar
- WHO: Weekly epidemiological bulletin - disease early warning system and response in Pakistan. In Book weekly epidemiological bulletin - disease early warning system and response in Pakistan (Editor ed.^eds.). Pakistan: City: National Institute of Health, Islamabad and World Health Organization (WHO); 2014:1–7.Google Scholar
- Estrada-Pena A, Zatansever Z, Gargili A, Aktas M, Uzun R, Ergonul O, et al.: Modeling the spatial distribution of crimean-congo hemorrhagic fever outbreaks in Turkey. Vector Borne Zoonotic Dis 2007, 7:667–78. 10.1089/vbz.2007.0134View ArticlePubMedGoogle Scholar
- Warden CR: Comparison of Poisson and Bernoulli spatial cluster analyses of pediatric injuries in a fire district. Int J Health Geogr 2008, 7:51. 10.1186/1476-072X-7-51View ArticlePubMed CentralPubMedGoogle Scholar
- Atkinson B, Chamberlain J, Jameson LJ, Logue CH, Lewis J, Belobrova EA, et al.: Identification and analysis of Crimean-Congo hemorrhagic fever virus from human sera in Tajikistan. Int J Infect Dis 2013, 17:e1031–7. 10.1016/j.ijid.2013.04.008View ArticlePubMedGoogle Scholar
- Khushal , Athar M, Kashif , Dawood M, Ahmed MD, Saleem MH: Epidemiological trends, clinical profile and geographical distribution of Crimean Congo Hemorrhagic Fever in Quetta, Balochistan. Pakistan J Sci 2014, 66:1–6.Google Scholar
- Ahmad S, Aslam M, Mirza SN: Rangeland degradation and management approaches in Baluchistan, Pakistan. Pak J Bot 2012, 44:127–36.Google Scholar
- Shafique M, Kakar MA: Current livestock marketing and its future prospects for the economic development of Balochistan–Pakistan. Int J Agric Biol 2006, 6:885–95.Google Scholar
- Mostafavi E, Haghdoost A, Khakifirouz S, Chinikar S: Spatial analysis of Crimean Congo hemorrhagic fever in Iran. AmJTrop Med Hyg 2013, 89:1135–41. 10.4269/ajtmh.12-0509View ArticleGoogle Scholar
- Saleem J, Usman M, Nadeem A, Sethi SA, Salman M: Crimean–Congo hemorrhagic fever: a first case from Abbottabad, Pakistan. Int J Infect Dis 2009, 13:e121–3. 10.1016/j.ijid.2008.07.023View ArticlePubMedGoogle Scholar
- Sargianou M, Panos G, Tsatsaris A, Gogos C, Papa A: Crimean-Congo hemorrhagic fever: seroprevalence and risk factors among humans in Achaia, western Greece. Int J Infect Dis 2013, 17:e1160–5. 10.1016/j.ijid.2013.07.015View ArticlePubMedGoogle Scholar
- Yagci-Caglayik D, Korukluoglu G, Uyar Y: Seroprevalence and risk factors of Crimean-Congo hemorrhagic fever in selected seven provinces in Turkey. J Med Virol 2014, 86:306–14. 10.1002/jmv.23699View ArticlePubMedGoogle Scholar
- Osei FB, Duker AA, Stein A: Evaluating spatial and space-time clustering of cholera in Ashanti-Region-Ghana. In InTech Edited by: Gowder S. 2012. doi:10.5772/36316. Available from: http://cdn.intechopen.com/pdfs/34049/InTech-Evaluating_spatial_and_space_time_clustering_of_cholera_in_ashanti_region_ghana.pdf Google Scholar
- Chainey S, Ratcliffe J: GIS and crime mapping. West Sussex, England: John Wiley & Sons; 2005.View ArticleGoogle Scholar
- Shekhar S, Xiong H: Encyclopedia of GIS. New York: Springer; 2008.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.