Skip to main content
  • Research article
  • Open access
  • Published:

Environmental correlates of crimean-congo haemorrhagic fever incidence in Bulgaria



Crimean-Congo Haemorrhagic Fever (CCHF) is a zoonotic viral disease transmitted by ixodid tick bites, mainly of Hyalomma spp., or through contact with blood/tissues from infected people or animals. CCHF is endemic in the Balkan area, including Bulgaria, where it causes both sporadic cases and community outbreaks.


We described trends of CCHF in Bulgaria between 1997 and 2009 and investigated the associations between CCHF incidence and a selection of environmental factors using a zero-inflated modelling approach.


A total of 159 CCHF cases (38 women and 121 men) were identified between 1997 and 2009. The incidence was 0.13 cases per 100,000 population/year with a fatality rate of 26%. An epidemic peak was detected close to the Turkish border in the summer of 2002. Most cases were reported between April and September. Increasing mean temperature, Normalized Difference Vegetation Index (NDVI), savannah-type land coverage or habitat fragmentation increased significantly the incidence of CCHF in the CCHF-affected areas. Similar to that observed in Turkey, we found that areas with warmer temperatures in the autumn prior to the case-reporting year had an increased probability of reporting zero CCHF cases.


We identified environmental correlates of CCHF incidence in Bulgaria that may support the prospective implementation of public health interventions.

Peer Review reports


Crimean-Congo Haemorrhagic Fever (CCHF) is a zoonotic disease caused by a Nairovirus of the Bunyaviridae family that is usually transmitted by the bite of ixodid ticks, particularly of the Hyalomma spp. [1, 2], or through contact with tissues, blood or other bodily fluids from infected people and animals [3, 4].

Hyalomma ticks feed only once at each developmental stage. To serve as a CCHF vector, they must become infected at one stage, transmit the virus transtadially to the next stage or transovarially to the next generation, and then transmit the virus by bite to another vertebrate host [5]. Environmental changes can influence both the survival and reproduction of Hyalomma ticks, as well as trigger community outbreaks [6].

CCHF is one of the most feared viral haemorrhagic fevers because of its high fatality rate in human cases (up to 30%, depending on the transmission route) [1, 7, 8] and the potential to cause both community and nosocomial outbreaks [3]. The ratio of sub-clinical to clinical manifestations is approximately 5:1; thus, up to 80% of infections may be asymptomatic [5, 9]. The seroprevalence of CCHF in people with a history of tick bite may be as high as 20% in endemic areas [10]. High-risk groups include farmers, livestock breeders, veterinarians, border army units, abattoir and healthcare workers [11]. Because of exposure to tick bites, outdoor recreational activities are also a recognized risk factor for CCHF in endemic areas [11].

CCHF is endemic in the Balkan area, including Bulgaria [1, 12]. Most cases occur between April and September in the southeast and the south-central part of the country, with a limited number of cases involved per outbreak [1]. These cases follow by approximately one month the appearence of adult forms of Hyalomma marginatum, the main vector species of CCHF in Bulgaria [1]. H. marginatum is believed to be ubiquitous in Bulgaria (Christova I, personal communication). However, field data on its geographical distribution are not available.

The first documented outbreak in Bulgaria occurred during the agricultural collectivization in 1953 [1]. After the introduction of a vaccination programme of high-risk groups of the population in 1974 [12], cases decreased from 1,105 in 1953–1974 to 279 in 1975–1996 [8, 13], to 196 in 1997–2008 [10]. Although the overall number of cases decreased over time, a major outbreak was observed in 2008, with a cluster of cases in the southwest part of Bulgaria, an area historically considered at low risk for CCHF outbreaks [13].

Poverty, human migrations, lack of vector control and inadequate management of infected patients among healthcare workers have been associated with an increased risk of acquiring CCHF [5, 6, 11, 14]. Furthermore, neglect of agricultural lands, agricultural reforms resulting in landscape alterations and global warming may play a relevant role in CCHF emergence as well, although so far no empirical proof has been provided to support this latter claim [15, 16].

In animals, including livestock and hares, CCHF is asymptomatic [12]. Domestic ruminants (e.g. cattle, sheep and goats) are considered to be short-lived reservoirs of the virus [1], and hares are important amplifying hosts, as they are commonly infested with immature stages of ticks [1, 8]. The enzootic cycle of CCHF therefore depends on a complex network of ticks and host populations, and it is likely to be much more widespread than that revealed by the occurrence of clinical cases [17]. Generally speaking, the risk of CCHF infection in humans increases suddenly when the potential for enzootic transmission and/or the degree of exposure increase as well; thus, the investigation of those environmental factors that may influence the cycle of CCHF is relevant for outbreak preparedness and response [12].

In this study, we aimed at identifying those environmental determinants that may support the development of risk-based surveillance activities for CCHF in Bulgaria. Specifically, we investigated the association between the incidence of CCHF in humans and a selection of landscape and climatic factors that may be relevant to explain the occurrence of CCHF cases in Bulgaria.


Human data

A retrospective study was carried out to identify cases of CCHF that occurred in the resident population of Bulgaria between 1997 and 2009. Cases were collected from surveillance registry data and by contacting local authorities. Cases were laboratory confirmed by viral detection using virus cultivation AND/OR serology by ELISA and complement fixation assay.

Information on age, gender, month of disease onset, and district of residence of the cases was collected from clinical records. Resident population figures for each district of Bulgaria were obtained from the census data of the year 2006. Cases and respective underlying populations were subsequently allocated to each municipality of the districts they belonged to, by month and year of disease onset.

Environmental data

Environmental variables used in the study had a potential impact on Hyalomma tick biology and therefore on the possible occurrence of human CCHF cases. Environmental data for the study area, i.e. Bulgaria with its 253 municipalities, were collected from several sources for the period 2001–2009.

Data on the monthly mean land surface ambient temperature (°C) and medium resolution NDVI, a measure of photosynthetic activity and large-scale vegetation distribution [18], were obtained from the NEO (NASA Earth Observations) datasets ( These are layers derived from radiances measured by the NOAA-AVHRR orbiting satellites that provide well-calibrated datasets at a resolution of 0.1° already corrected for atmospheric and orbital disturbances. Both temperature and NDVI values are calculated by NEO as the monthly mean of daily maximum values. Municipality-averaged monthly values of temperature and NDVI were calculated for each month and year of the study period.

Landscape information was obtained from the 2004–2006 ESA Globcover Land Cover layer derived from an automatic and regionally-tuned classification of time series of MERIS-FR composites. This land cover classification is broken into 22 land cover classes which are defined according to the UN Land Cover Classification System (LCCS) [19]. For each municipality, the percent coverage of each land cover class was calculated. A landscape fragmentation index was also defined as the total number of different land cover classes (i.e. patch categories) that fell in the same municipality, normalized by the corresponding municipality area.

Data on cattle and sheep density (number of animals per Km2) were obtained from the FAO’s gridded livestock density datasets ( Animal densities are calculated by the FAO from official census and survey data accounting for the amount of land suitable for livestock production. Based on the statistical relationship with some environmental variables in similar agro-ecological zones, livestock densities are predicted to match national census totals for the year 2005 according to the FAOSTAT database and are provided at a resolution of 0.05° [20]. Cattle and sheep mean densities were calculated for each municipality.

All spatial data manipulations were performed using ArcMap in ArcGIS (ESRI Inc. Redlands, CA, USA).

Statistical analysis

As the space-time distribution of CCHF cases was positively skewed with a high proportion of zeros, zero-inflated modelling was used to deal with these count data. Zero-inflated models account for “excess zeros” (i.e. zero counts greater than expected based on the non-zero counts) by assuming that the zeros arise from two different processes [21, 22]. The first process occurs with probability pi and produces only zeros (structural zeros), whereas the second process occurs with probability (1 − pi) Pr(Ki = 0), where Pr(Ki = 0) is a Poisson or Negative Binomial (NB) probability function with the zero events occurring by chance (sampling zeros). The overall probability of zero counts is a combination of the probabilities of the two processes. Positive counts (i.e. greater than zero) are assumed to occur with probability (1 − pi) Pr(Ki = yi), where Pr(Ki = yi) is a Poisson or NB probability function for the positive counts [21, 22]. In our study, the structural zero process accounted for those districts in which there were no CCHF cases at all during the entire study period (apparently “CCHF-free” areas). The sampling zero process accounted for those districts were the counts of CCHF cases were either zero or values greater than zero depending on the month and year in question during the entire study period.

As the overdispersion parameter α was significantly different from zero and the Vuong likelihood ratio test revealed that the zero-inflated model was preferred against the standard negative binomial model, a Zero-Inflated Negative Binomial (ZINB) model was chosen, as it accounted both for the excess zeros and for the extra heterogeneity in the outcome [21, 22].

The probability of the structural zero process and the number of event counts in the sampling zero process may depend on different covariates [21, 22]. This allowed us to explain the occurrence of CCHF cases differently: the covariates in the zero-inflated term of the ZINB model modelled the left tail of the case distribution and explained the zero counts in the “CCHF-free” areas, whereas those in the NB part of the ZINB model modelled the right tail and explained the positive counts in the “CCHF-affected” areas. A logit link was used for the excess zeros and a log link for the mean of the NB distribution.

In the NB part of the ZINB model (sampling zero process), the following covariates were included to explain the positive counts of CCHF cases: 1) monthly mean temperature; 2) monthly mean NVDI; 3) habitat fragmentation index (low, medium, and high, according to tertiles), 4) percent coverage of grasslands and areas with predominant scrub and herbaceous vegetation; 5) livestock density; 6) year of the study period (treated as linear splines [23, 24]); and 7) human density (individuals per Km2).

In the logit part of the ZINB model (structural zero process), the mean temperature in October-November of the previous year and the year of the study period (treated as linear splines) were included to explain the excess of zero cases, as it has been shown that late autumn temperatures before the sudden drop in winter temperatures may differ significantly between areas with and without reported human cases of CCHF [18].

The underlying human population was included as exposure variable in the ZINB model. The model accounted for within-municipality correlation resulting from repeated measurements made over time on the same municipalities using robust variance estimators. Statistical analysis was performed using STATA 10 (StataCorp LP, Lakeway Drive, TX, USA), and the significance level was set to 0.05.


Between 1997 and 2009, 159 cases (38 women and 121 men) were identified (Table 1), including 42 fatal cases (case fatality rate: 26%). The average age (years) of cases was 52.80 ± 17.32 standard deviation (median: 54.50). Cases were mainly reported from April to September

Table 1 Number of cases, mean annual incidence per 100,000 and respective 95% confidence intervals (CI) for CCHF in Bulgaria between 1997 and 2009 by district

Cumulative incidence of CCHF was 2.07 per 100,000 population. The incidence per 100,000/year (1997–2009) was 0.13 (95% confidence interval [CI]: 0.11-0.15). There were large variations in CCHF incidence among the districts (Table 1). Values ranged from 0.02 (95% CI: 0.00-0.18) to 0.59 (95% CI: 0.47-0.74) per 100,000/year (1997–2009). The incidence did not vary substantially between 1997 and 2009: the mean incidence per 100,000/year was 0.13 (95% CI:0.10-0.17) during 1997–2000; 0.20 (95% CI: 0.16-0.25) during 2001–2004; and 0.08 (95% CI: 0.06-0.11) during 2005–2009.

The temporal pattern of CCHF cases from 1997 to 2009 showed recurrent peaks in the summer, with a main peak in the summer of 2002 in the southeast and northeast districts of the country, particularly the districts of Stara Zagora, Burgas, Varna, and Haskovo, where 30 cases were reported within the same month of June (Figure 1). The first cases occurred in the districts of Stara Zagora and Burgas, close to the borders with Turkey. In 2008, another cluster of 12 cases was reported in the southwest part of the country, an area previously considered at low risk for CCHF.

Figure 1
figure 1

Distribution of CCHF cases by month and year of occurrence.

Results from the ZINB model are reported in Table 2. According to the model characteristics, the results are presented for its two constitutive parts separately: the first part (Table 2-A), based on the sampling zeros process, shows the pattern of associations between the explanatory variables and the incidence of CCHF in the “CCHF-affected” areas, while the second part (Table 2-B), based on the structural zero process, presents the association of the selected explanatory variable with the consistent presence of zero CCHF cases in the “CCHF-free” areas.

Table 2 Incidence rate ratios (IRR), Odds ratios (OR) and 95% confidence intervals (CI) from the zero inflated- negative binomial regression model

In those districts where CCHF cases were observed, a unit increase in the mean temperature and NDVI values was significantly associated with an increase in the incidence rate by 5.5% (Incidence Rate Ratio [IRR]: 1.055), and by 1.8% (IRR: 1.018), respectively. Compared to areas with a lowly fragmented land structure, those with a high and medium level of fragmentation were significantly associated with an increase in CCHF incidence (IRR: 1.558; and IRR: 1.402, respectively). In the districts with a high and medium coverage of grassland, scrub and herbaceous vegetation, the incidence was 4.261 and 3.994 times the incidence of the reference group (i.e. the districts where this type of vegetation was less represented). There was no evidence of association between CCHF incidence and human density, cattle density, and sheep density.

The absence of CCHF cases in the “CCHF-free” areas was significantly associated with the mean temperature from October to November (Odds Ratio [OR]: 1.129), meaning that the higher the temperatures in October-November the higher the probability of reporting zero cases in the following year.


Bulgaria has been challenged by CCHF since the Second World War, becoming only the European country in which the disease is so far officially considered to be endemic, with recurrent epidemics particularly in the southeast part of the country. A similar epidemiological situation is increasingly being documented also in neighbouring countries, such as Albania, Kosovo, Russia, Turkey, Ukraine and Greece [25].

In this study, we described recent patterns in CCHF incidence in Bulgaria and investigated the extent to which the considered environmental factors were associated with the observed pattern of CCHF incidence in Bulgaria between 2001 and 2009, using a zero-inflated modelling approach.

Explaining zero counts of CCHF cases

There were many areas with zero cases of CCHF during the entire study period. This occurs frequently with vector-borne diseases, as transmission to humans may depend on several biotic and abiotic factors as well as their interactions that may vary substantially over space and time. Moreover, it has been shown that CCHF does not occur in all the areas characterized by a suitable habitat for H. marginatum, but only in a narrow subset of them, as observed during the Turkish epidemic in 2003–2008 [18].

We observed that increasing mean temperature in late autumn (i.e. October-November) increased the likelihood of reporting zero counts of CCHF cases in the following year. A similar finding has previously been observed in Turkey, where an analysis of the environmental factors that differed between areas with and without CCHF cases showed that the areas with CCHF cases had colder mean temperatures in the late autumn and winter (November-February) compared to those without CCHF cases [18]. As immatures moult in the late autumn at that latitude, high temperatures after moulting would trigger their questing activity for hosts, whereas low temperatures would send them into quiescence. It has therefore been hypothesized that in areas with colder temperatures after moulting, more abundant overwintering tick populations would be synchronously activated in the spring of the following year, thereby increasing the potential for viral circulation [18].

Explaining non-zero counts of CCHF cases

We found that the incidence of CCHF increased with increasing mean temperature in areas with CCHF cases. Temperature has previously been identified as the main factor driving the seasonal pattern and abundance of H. marginatum [26, 27], as high temperatures, especially in the spring and summer, may accelerate H. marginatum cycle by switching on its interstadial development [28] and host questing activity [3, 26, 28]. Moreover, humans may adjust their recreational and subsistence outdoor activities according to the weather, increasing the chance of encountering and being bitten by these ticks [14, 29].

The incidence of CCHF was higher in those areas characterized by a high proportion of grasslands, scrub and herbaceous vegetation (savannah-type environment), where in fact Hyalomma ticks usually find their ecological optima [2, 30].

Livestock density was not significantly associated with CCHF incidence. In South Africa, the prevalence of antibody against CCHF virus has been found to be associated with the handling of sheep [1], suggesting that people with occupational exposure to livestock are particularly at-risk for CCHF. However, veterinarians and farm workers are usually vaccinated against CCHF in Bulgaria, and this may be the reason as to why this association was insignificant in our study.

Areas with a higher incidence of CCHF were those characterized by a highly fragmented habitat. A similar association has been found in previous studies [18, 31, 32], supporting the hypothesis that a fragmented land structure may increase the risk of acquiring CCHF by favouring viral circulation and amplification though frequent at-risk contacts between ticks, humans, livestock and wildlife [18]. Habitat fragmentation has also been identified as a critical factor for the distribution of Lyme disease [33], and recent modelling studies have suggested that increasing host biodiversity (species richness) may decrease the risk of acquiring Lyme disease through a dilution effect, whereas habitat fragmentation may increase this risk by reducing wildlife biodiversity [34]. However, it is difficult to conceptualize how the dilution effect and the impact of reduced biodiversity actually apply to CCHF, a disease characterized by a variety of different transmission and tick hosts. Possibly, the small fragments would increase appropriate habitat for ticks and their hosts, such as high densities of small mammals due to lower densities of predators and/or competitors, resulting in a critical rise in the rate of at-risk contacts.

We found a significant positive association between NDVI and CCHF incidence. A similar correlation has previously been found in Turkey [35], suggesting that the NDVI may proxy for tick seasonal activity. The NDVI reflects the amount and vigour of healthy green biomass during the seasonal cycle of vegetation as a surrogate of the seasonal fluctuations in weather and has long been recognized as one of the most important environmental variables influencing habitat suitability for ticks [32, 36, 37], as well as in regulating tick survival rates [27, 38]. This is also partly due to the growth and dispersal rates of vertebrate host populations, particularly wild herbivores, which are strongly influenced by the quality and quantity of the available forage vegetation. Consistent with previous findings [39], the largest increase in NDVI occurred in the early growing season, i.e. the spring, but a relatively large NDVI was also found in the summer, which apparently accounted for the significant positive association with CCHF incidence in the “CCHF-affected” areas we found. As the effect of NDVI was still significant even after controlling for the land cover that best characterizes H. marginatum in the ZINB model, it is likely that the observed association with NDVI reflected the natural seasonal fluctuations of environmental exposure to abiotic factors (e.g. intensity and length to sunlight exposure) that can influence tick phenology and fitness, and consequently CCHF occurrence.


Infection rates of ticks and vertebrate hosts as well as frequency of at-risk contacts between humans and ticks are major determinants for CCHF infection in humans [40]. However, this kind of information and quantitative data on the abundance of ticks and vertebrate reservoirs (i.e. non-migratory terrestrial mammals, such as hares, small rodents, hedgehogs and resident bird species) were not available in our study. This may have limited our capacity to fully explain the observed pattern of CCHF in Bulgaria.

Another limitation is the use of administrative boundaries (i.e. districts and municipalities) as territorial units, as the lack of coordinates of cases did not allow for the use of other geographical resolutions. However, the use of such territorial units allowed for intra-unit variability in both biotic and abiotic factors, making it possible to generalize model predictions. This is because infection may have been acquired away from home, but possibly within the same municipality/district of residence location. Moreover, quality of data for probable sites of exposure (point locations) may vary according to differences in information gathering practice and interest of individual physicians, especially regarding asymptomatic or mild cases who have not been recognized as such, thereby precluding reliable use of any point locations as the basis for determining environmental correlates of CCHF incidence in Bulgaria.


This study employed both GIS-derived information and classical statistical techniques to investigate CCHF patterns in Bulgaria and found that increasing mean temperature, NDVI, savannah-type land coverage and habitat fragmentation increased significantly the risk of CCHF in the “CCHF-affected” areas. Similar to that observed in Turkey, the temperature in the late autumn prior to the case-reporting year seemed to influence the probability of occurrence of zero CCHF cases in the “CCHF-free” areas. However, to further investigate the driving environmental factors behind CCHF patterns, data on tick abundance and infection rates, wildlife populations and societal factors are required.

There is urgent need to improve both passive and active surveillance activities for CCHF in Bulgaria and to evaluate vaccine coverage in high-risk areas. A holistic approach to risk assessment of CCHF under global climate change is also required to estimate the impact on, and the vulnerability of, human populations in relation to ticks and the prevalence of infection in vertebrate hosts under different environmental scenarios [41].

Despite the wealth of research documenting the widespread distribution, diversity of possible vectors and potential reservoirs for CCHF, our understanding of its epidemiology remains incomplete and at times we can only speculate on the importance of the different variables that can influence CCHF transmission cycle and stimulate outbreaks. In this regard, the factors underlying the spread of CCHF in Bulgaria we found, although they have not yet fully clarified, can make a substantial contribution towards fulfilling the information needs of public health policy-makers when challenging CCHF in varying eco-epidemiological dimensions.

Key points

CCHF is endemic in Bulgaria where it causes either sporadic cases and community outbreaks.

A total of 159 CCHF cases (38 women and 121 men) were identified between 1997 and 2009 (incidence: 0.13 cases per 100,000 population/year; fatality rate: 26%).

Increasing temperature, NDVI, savannah-type land coverage, and habitat fragmentation increased the risk of CCHF in areas affected from the disease.

Areas with lower winter temperatures were associated with an increased probability of reporting zero CCHF.


  1. Hoogstraal H: The epidemiology of tick-borne Crimean-Congo hemorrhagic fever in Asia, Europe, and Africa. J Med Entomol. 1979, 22: 307-417.

    Article  Google Scholar 

  2. Manilla G: Acari Ixodidae, Fauna d’Italia. Parte speciale: Ixodidae. Acari, Ixodida. 1998, 193-195.

    Google Scholar 

  3. Ergonul O: Crimean-Congo haemorrhagic fever. Lancet Infect Dis. 2006, 6: 203-214. 10.1016/S1473-3099(06)70435-2.

    Article  PubMed  Google Scholar 

  4. Whitehouse CA: Crimean-Congo hemorrhagic fever. Antiviral Res. 2004, 64: 145-160.

    Article  CAS  PubMed  Google Scholar 

  5. Gale P, Estrada-Pena A, Martinez M, Ulrich RG, Wilson A, Capelli G, et al: The feasibility of developing a risk assessment for the impact of climate change on the emergence of Crimean-Congo haemorrhagic fever in livestock in Europe: a Review. J Appl Microbiol. 2010, 108: 1859-1870.

    CAS  PubMed  Google Scholar 

  6. Randolph S, Ergonul O: Crimean-Congo hemorrhagic fever: exceptional epidemic of viral hemorrhagic fever in Turkey. Future Medicine. 2008, 3: 303-306.

    Google Scholar 

  7. Papa A, Christova I, Papadimitriou E, Antoniadis A: Crimean-Congo hemorrhagic fever in Bulgaria. Emerg Infect Dis. 2004, 10: 1465-1467. 10.3201/eid1008.040162.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Charrel RN, de’Lamballerie X: [West Nile virus, an emerging arbovirus]. Presse Med. 2004, 33: 1521-1526. 10.1016/S0755-4982(04)98977-4.

    Article  PubMed  Google Scholar 

  9. Goldfarb LG, Chumakov MP, Myskin AA, Kondratenko VF, Reznikova OY: An epidemiological model of Crimean hemorrhagic fever. AmJTrop Med Hyg. 1980, 29: 260-264.

    CAS  Google Scholar 

  10. Leblebicioglu H: Crimean-Congo haemorrhagic fever in Eurasia. Int J Antimicrob Agents. 2010, 36 (Suppl 1): S43-S46.

    Article  CAS  PubMed  Google Scholar 

  11. Vorou R, Pierroutsakos IN, Maltezou HC: Crimean-Congo hemorrhagic fever. Curr Opin Infect Dis. 2007, 20: 495-500. 10.1097/QCO.0b013e3282a56a0a.

    Article  PubMed  Google Scholar 

  12. Maltezou HC, Andonova L, Andraghetti R, Bouloy M, Ergonul O, Jongejan F, et al: Crimean-Congo hemorragic fever in Europe: current situation calls for preparedness. Euro Surveill. 2010, 15: pii=19504

    Google Scholar 

  13. Christova I, Di CA, Papa A, Castilletti C, Andonova L, Kalvatchev N, et al: Crimean-Congo hemorrhagic fever, southwestern Bulgaria. Emerg Infect Dis. 2009, 15: 983-985. 10.3201/eid1506.081567.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Randolph SE: To what extent has climate change contributed to the recent epidemiology of tick-borne diseases?. Vet Parasitol. 2010, 167: 92-94. 10.1016/j.vetpar.2009.09.011.

    Article  PubMed  Google Scholar 

  15. Estrada-Pena A: Forecasting habitat suitability for ticks and prevention of tick-borne diseases. Vet Parasitol. 2001, 98: 111-132. 10.1016/S0304-4017(01)00426-5.

    Article  CAS  PubMed  Google Scholar 

  16. Gubler DJ, Reiter P, Ebi KL, Yap W, Nasci R, Patz JA: Climate variability and change in the United States: potential impacts on vector- and rodent-borne diseases. Environ Health Perspect. 2001, 109 (Suppl 2): 223-233.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Randolph SE, Rogers DJ: The arrival, establishment and spread of exotic diseases: patterns and predictions. Nat Rev Microbiol. 2010, 8: 361-371. 10.1038/nrmicro2336.

    Article  CAS  PubMed  Google Scholar 

  18. Estrada-Pena A, Vatansever Z, Gargili A, Ergonul O: The trend towards habitat fragmentation is the key factor driving the spread of Crimean-Congo haemorrhagic fever. Epidemiol Infect. 2010, 138: 1194-1203. 10.1017/S0950268809991026.

    Article  CAS  PubMed  Google Scholar 

  19. Bontempsm S, Defourny P, Van Bogaert E, Arino O, Kalogirou V, Ramos Perez J: GlobCover Project. 2009: Products description and validation report, European Space Agency and Universitè catholique Louvain, Brussels, Belgium. 2011, Ref Type: Report

    Google Scholar 

  20. Robinson TP, Franceschini G, Wint W: The Food and Agriculture Organization's Gridded Livestock of the World. Vet Ital. 2007, 43: 745-751.

    PubMed  Google Scholar 

  21. Lambert D: Zero-inflated Poisson regression with application to defects in manufacturing. Technometrics. 1992, 34: 1-14. 10.2307/1269547.

    Article  Google Scholar 

  22. Heilbron DC: Zero-altered and other regression models for count data withadded zeros. Biom J. 1994, 36: 531-547. 10.1002/bimj.4710360505.

    Article  Google Scholar 

  23. Ahlberg JH, Nilson EN, Walsh JL: Fundamental properties of generalized splines. Proc Natl Acad Sci U S A. 1964, 52: 1412-1419. 10.1073/pnas.52.6.1412.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Rothman KJ, Greenland S: Regression splines. 1998, Philadelphia: Lippincott Williams andWilkins Publishers, 392-394.

    Google Scholar 

  25. ECDC: Consultation on Crimean-Congo haemorragic fever prevention and control. 2008, Ref Type: Report

    Google Scholar 

  26. Estrada-Pena A, Martinez AM, Munoz Reoyo MJ: A population model to describe the distribution and seasonal dynamics of the tick Hyalomma marginatum in the Mediterranean Basin. Transbound Emerg Dis. 2011, 58: 213-223. 10.1111/j.1865-1682.2010.01198.x.

    Article  CAS  PubMed  Google Scholar 

  27. Estrada-Pena A: Climate, niche, ticks, and models: what they are and how we should interpret them. Parasitol Res. 2008, 103 (Suppl 1): S87-S95.

    Article  PubMed  Google Scholar 

  28. Gray JS, Dautel H, Estrada-Pena A, Kahl O, Lindgren E: Effects of climate change on ticks and tick-borne diseases in europe. Interdiscip Perspect Infect Dis. 2009, 200 (9): 593232-

    Google Scholar 

  29. Randolph SE, Asokliene L, vsic-Zupanc T, Bormane A, Burri C, Gern L, et al: Variable spikes in tick-borne encephalitis incidence in 2006 independent of variable tick abundance but related to weather. Parasit Vectors. 2008, 1: 44-10.1186/1756-3305-1-44.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Emelyanova N: Seasonal changes and host adaptability of ticks of the species Hyalomma marginatum in the Stavropol territory. Zurnal Mikrobiologii, epidemiologii I immunobiologii. 2005, 4: 115-118.

    Google Scholar 

  31. Estrada-Pena A, Venzal JM: Climate niches of tick species in the Mediterranean region: modeling of occurrence data, distributional constraints, and impact of climate change. J Med Entomol. 2007, 44: 1130-1138. 10.1603/0022-2585(2007)44[1130:CNOTSI]2.0.CO;2.

    Article  PubMed  Google Scholar 

  32. 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-678. 10.1089/vbz.2007.0134.

    Article  PubMed  Google Scholar 

  33. Jackson LE, Hilborn ED, Thomas JC: Towards landscape design guidelines for reducing Lyme disease risk. Int J Epidemiol. 2006, 35: 315-322.

    Article  PubMed  Google Scholar 

  34. Schmidt KA, Ostfeld RS: Biodiversity and the dilution effect in disease ecology. Ecology. 2001, 82: 609-619. 10.1890/0012-9658(2001)082[0609:BATDEI]2.0.CO;2.

    Article  Google Scholar 

  35. Estrada-Pena A, Vatansever Z, Gargili A, Buzgan T: An early warning system for Crimean-Congo haemorrhagic fever seasonality in Turkey based on remote sensing technology. Geospat Health. 2007, 2: 127-135.

    Article  PubMed  Google Scholar 

  36. Sumilo D, Bormane A, Asokliene L, Lucenko I, Vasilenko V, Randolph S: Tick-borne encephalitis in the Baltic States: identifying risk factors in space and time. Int J Med Microbiol. 2006, 296 (Suppl 40): 76-79.

    Article  PubMed  Google Scholar 

  37. Estrada-Pena A, Quiez J, Sanchez AC: Species composition, distribution, and ecological preferences of the ticks of grazing sheep in north-central Spain. Med Vet Entomol. 2004, 18: 123-133. 10.1111/j.0269-283X.2004.00486.x.

    Article  CAS  PubMed  Google Scholar 

  38. Rogers DJ: Satellites, space, time and the African trypanosomiases. Adv Parasitol. 2000, 47: 129-171.

    Article  CAS  PubMed  Google Scholar 

  39. Piao S, Fang J, Ji W, Guo Q, Ke J, Tao S: Variation in a satellite-based vegetation index in relation to climate in China. J Veg Sci. 2004, 15: 219-226. 10.1658/1100-9233(2004)015[0219:VIASVI]2.0.CO;2.

    Article  Google Scholar 

  40. Randolph SE: The shifting landscape of tick-borne zoonoses: tick-borne encephalitis and Lyme borreliosis in Europe. Philos Trans R Soc Lond B Biol Sci. 2001, 356: 1045-1056. 10.1098/rstb.2001.0893.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Rezza G: Re-emergence of Chikungunya and other scourges: the role of globalization and climate change. Ann Ist Super Sanità. 2008, 44: 315-318.

    PubMed  Google Scholar 

Pre-publication history

Download references


This study was funded by the Italian National Public Health Institute. Retrospective study of anonymous data. No ethics approval required.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Fenicia M Vescio.

Additional information

Competing interests

The authors’ declare that they have no competing interests.

Authors’ contribution

FV, LB, LMG and GR conceived and designed the study, performed/interpreted the analysis, and prepared the manuscript for publication. IC and ET carried out the laboratory work, collected the data and contributed to the preparation of the manuscript. CK and LA supported the interpretation of results and contributed to the preparation of the manuscript. All authors have read and approved the final version of the manuscript.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Rights and permissions

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 (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Vescio, F.M., Busani, L., Mughini-Gras, L. et al. Environmental correlates of crimean-congo haemorrhagic fever incidence in Bulgaria. BMC Public Health 12, 1116 (2012).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: