Emerg Infect DisEIDEmerging Infectious Diseases1080-60401080-6059Centers for Disease Control and Prevention19523295272733207-166610.3201/eid1506.071666DispatchIncreasing Incidence of Zoonotic Visceral Leishmaniasis on Crete, GreeceZoonotic Visceral Leishmaniasis on CreteAntoniouMariaMessaritakisIppokratisChristodoulouVasilikiAscoksilakiIoannaKanavakisNikosSuttonAndrew J.CarsonConnorCourtenayOrinUniversity of Crete, Heraklion, Greece (M. Antoniou, I. Messaritakis, V. Christodoulou)Veterinary Laboratory of Heraklion, Heraklion (I. Ascoksilai)Veterinary Services, Sitia, Greece (N. Kanavakis)University of Warwick, Coventry, UK (A.J. Sutton, C. Carson, O. Courtenay)Address for correspondence: Orin Courtenay, Populations and Disease Research Group, Department of Biological Sciences, University of Warwick, Coventry CV4 7AL, UK; email: orin.courtenay@warwick.ac.uk62009156932934

To determine whether the incidence of canine leishmaniasis has increased on Crete, Greece, we fitted infection models to serodiagnostic records of 8,848 dog samples for 1990–2006. Models predicted that seroprevalence has increased 2.4% (95% confidence interval 1.61%–3.51%) per year and that incidence has increased 2.2- to 3.8-fold over this 17-year period.

Keywords: Leishmania infantumzoonotic visceral leishmaniasiscanine leishmaniasisincidencemodelszoonosesparasitesCreteGreecedispatch

Zoonotic visceral leishmaniasis (ZVL) caused by Leishmania infantum is a disease of humans and domestic dogs (the reservoir) transmitted by phlebotomine sandflies. According to the World Health Organization (1), ZVL was first recorded on Crete in 1907, after which it featured prominently in medical literature as a serious public health problem. In Chania, Crete, the annual incidence in the 1930s was 50 cases/30,000 population (2), and 21% of 1,115 dogs were positive for ZVL by the formol-gel serologic test; 30% of those were symptomatic (1).

After World War II, use of DDT against malaria vectors and focal destruction of Leishmania spp.–infected dogs is thought to have reduced ZVL on Crete (3). Sandflies were not found in villages systematically sprayed during 1946–1949 compared with unsprayed villages (4). During 1951–1975, only 33 alleged human ZVL cases were recorded on Crete (5), and reports of canine ZVL were scanty. In 1983 a serosurvey of 72 stray dogs in Chania identified only 1 infected dog (1).

By the late 1980s, Phlebotomus neglectus, the putative vector of L. infantum on Crete (6), was abundant in stone walls inside and outside villages around Heraklion (7). During 1999–2004, P. neglectus was found in abundance in human dwellings and rural locations (hollows in olive trees and near rodent burrows) (8). Since 1991, 38 persons who came to hospitals in Heraklion were confirmed as having cases of ZVL. Today, canine infection is confirmed throughout the island (M. Antoniou, unpub. data), and seroprevalences (30%–40%) are some of the highest reported in Europe. These accounts anecdotally suggest that the incidence of ZVL has increased on Crete during the past decade, which, if so, would be relevant to the public and veterinary health sectors.

The Study

We statistically examined diagnostic records of 8,848 dogs sampled in the eastern Crete districts of Lasithi and Heraklion during 1990–2006. Data were supplied by the Faculty of Medicine, University of Crete (n = 1,205 dogs from 97 villages tested by using an indirect immunofluorescent antibody test [IFAT], cutoff titer 160, during 1999–2006, accompanied by demographic and geographic records); and by the Ministry of Agriculture Serology Laboratory, Heraklion (n = 7,643 dogs tested by using an IFAT, cutoff titer 200, during 1999–2005, but without accompanying records). Samples were collected by veterinarians in private practice from any dog initially brought to their clinic for any reason, or by government veterinarians for any dog encountered, irrespective of clinical condition, during prearranged visits to villages. Numbers of villages and dogs sampled in any year depended on available resources at the time.

To reduce potential sampling bias, we first tested serologic data from 6 villages (located 9–45 km apart) consistently sampled annually during 1999–2006 (n = 744 dogs). The age-adjusted annual prevalence increased significantly with sample year (odds ratio [OR] 1.18, 95% confidence interval [CI] 1.089–1.271, p<0.001) and showed no significant differences in slope or intercept (by using robust standard errors) compared with data from 91 less consistently sampled villages with demographic records (n = 461 dogs; slope × village group interaction OR 0.91, 95% CI 0.800–1.023, p not significant; intercepts OR 5.35, 95% CI 0.902–3.688, p not significant).

These combined datasets were then compared with crude prevalence data for 1990–2005 (n = 7,643 dogs) calculated from the ministry records; no difference was detected in prevalence slopes (slope × data source interaction OR 1.06, 95% CI 0.995–1.30, p not significant; Figure). In univariable or multivariable logistic regression that controlled for dog age and clustering on villages (1999–2006, n = 1,205 dogs), no statistical confounding of the probability of a dog being seropositive was attributed to dog use (companion, guard, or hunting dog), sex, crude habitat type, or village altitude (p>0.215; see Technical Appendix Table). The final fit of seroprevalence against time was significant (OR 1.24, 95% CI 1.221–1.264, p<0.001), and the linearized difference in model prevalence (and binomial confidence limits) over the 17-year study equated to a mean prevalence increase of –ln(1 – 0.321)/16 = 2.4% (95% CI 1.61%–3.51%) per year.

To assess the change in infection incidence, our principal aim, we used 3 standard epidemiologic models (911; Technical Appendix) to calculate infection rates accounting for time, dog age, and potential loss of infection. The first method (model 1) used the cross-sectional age-prevalence data (IFAT cutoff titer 160), in which the proportion of seropositive dogs in each age class is fitted by varying the rates of infection and recovery. A second method (model 2) used these same data to describe the infection rate as it varied with both time and age until reaching the best fit. The third method (model 3) estimated the infection rate from longitudinal data of previously unexposed dogs <12 months of age that were followed up during 1 transmission season.

Results from the 3 models were consistent (Tables 1, 2) and showed a relative increase in the mean infection rate estimated to be 2.20–3.78-fold higher during 2005–2006 than during 1999–2000. The models differed in approach and age of dogs considered by necessity of the model, number of estimated parameters, or model reduction. Inclusion of a parameter describing loss of infection (Table 1, model 1) did not significantly lower the infection rate estimates as might be expected compared with a single parameter (Table 1, model 2) or longitudinal (Table 2, model 3) model, both of which identified younger (<2 years of age) dogs to be at substantially greater risk for infection (p<0.0001).

Variation in incidence over time estimated from cross-sectional data for 1,205 dogs with accompanying demographic records, Crete, Greece, 1990–2006*
PeriodModel 1
Model 2
No. dogs
Incidence/mo95% CILoss of infection/mo95% CIIncidence/mo95% CI
1999–20000.0160.0107–0.02060.0450.0260–0.06450.0150.0093–0.0213237
2001–20020.0290.0114–0.04550.0710.0201–0.12110.0230.0166–0.0298219
2003–20040.0300.0216–0.03810.0490.0320–0.06670.0290.0221–0.0367401
2005–20060.0590.0233–0.09460.1060.0383–0.174300.0320.0205–0.0477348

*Only model 1 is designed to estimate loss of infection (serorecovery). CI, confidence interval. A full description of these models is available in the Technical Appendix.

Variation in incidence over time estimated from longitudinal data for 179 dogs with accompanying demographic records, Crete, Greece, 1990–2006*
PeriodModel 3
No positive/ no. testedIncidence/ mo95% CI
1999–20006/560.0140.0061–0.0280
2001–20027/300.0330.0153–0.0612
2003–200412/560.0300.0172–0.0490
2005–200610/370.0390.0210–0.0670

*CI, confidence interval. A full description of this model is available in the Technical Appendix.

Conclusions

The potential contribution of any improvements in diagnostic test sensitivity or vigilance to the increasing incidence of ZVL infection is unclear. The difference in cutoff titers between data sources minimally shifted the absolute prevalence values, but not relative prevalence slopes, with time (Figure). Any loss of infection with age (Table 1, model 1) could result from nonmutually exclusive biologic processes including recovery from infection, death, or reduced past exposure (9,11). The latter possibility is unlikely on Crete because of the higher risk identified in young dogs in all biannual periods. Disproportionate numbers of deaths of seropositive dogs is not suggested by a decline in ZVL clinical signs in older dogs in this study (data not shown) or elsewhere (11,12). Loss of detectable Leishmania-specific antibody is the more likely explanation because the observed rates of serorecovery are not dissimilar to those (e.g., 0.062/month) estimated by cohort studies elsewhere (12).

Annual seroprevalence of zoonotic visceral leishmaniasis in dogs on Crete, Greece, 1990–2006. Shown are logistic fits of age-adjusted prevalences (line and squares) for dogs from 97 villages (indirect immunofluoresecent antibody test [IFAT] cutoff titer 160) and crude seroprevalences (line and triangles) calculated from records of the Veterinary Laboratory of Heraklion, Crete, Hellenic Ministry of Rural Development and Food (www.minagric.gr) (IFAT cutoff titer 200). Binomial standard error bars are shown.

Actual infection rates are likely to be higher than those shown here because IFAT sensitivity is <100%. Similarly, absolute prevalences, particularly low values for 1990–1991, should be treated with caution because the official leishmaniasis control program on Crete (1984–1995) began before this period when infection was presumably sufficient on the island to warrant intervention. The intervention comprised elimination of IFAT-seropositive dogs (cutoff titer 400) but did not include insecticide spraying (DDT spraying ceased in 1950; V. Chatzistefanou, pers. comm.).

The continual increase in canine seroprevalence during the latter part of the intervention (Figure) suggests that the culling policy was unsuccessful in reducing transmission; likely reasons for the low efficacy of dog culling in other leishmaniasis-endemic regions have been described (13). Officially, destruction of seropositive dogs is still required today unless the owner agrees to veterinary treatment of the dog or to keep the dog under sandfly-proof netting. However, there is no current policy on Crete to combat vectors. We conclude that the results of our study are consistent with a postwar reemergence and current increasing incidence of ZVL infection on Crete.

Supplementary MaterialTechnical Appendix

Analysis of Demographic and Environmental Associations with Canine Prevalence of Zoonotic Visceral Leishmaniasis

Suggested citation for this article: Antoniou M, Messaritakis I, Christodoulou V, Ascoksilaki I, Kanavakis N, Sutton A, et al. Increasing incidence of zoonotic visceral leishmaniasis on Crete, Greece. Emerg Infect Dis [serial on the Internet]. 2009 Jun [date cited]. Available from http://www.cdc.gov/EID/content/15/6/932.htm

Acknowledgments

We thank E. Svirinaki for technical assistance; M. Spyridakis, G. Parasyris, C. Parasyris, and K. Milathianakis for contributing dog samples; R.J. Quinnell for comments on an earlier version of the manuscript; and P. Ready and J.P. Dedet for valuable discussions.

Dr Antoniou is an assistant professor of parasitology at the Medical School of the University of Crete. Her research interests are Leishmania and Toxoplasma parasites, the epidemiology of zoonoses, and development of strategies for disease control.

ReferencesChaniotis B Leishmaniasis, sandfly fever and phlebotomine sandflies in Greece: an annotated bibliography. Geneva: World Health Organization; 1994 WHO/LEISH/94.34.Adler S, Theodor O, Witenberg G Investigations on Mediterranean kala azar. XI. A study of leishmaniasis in Canea (Crete).Proc R Soc Lond B Biol Sci 1938;125:491516 10.1098/rspb.1938.0039Hertig M Phlebotomus and residual DDT in Greece and Italy.Am J Trop Med Hyg 1949;29:77380918141397Hadjinicolaou J Present status of Phlebotomus in certain areas of Greece.Bull World Health Organ 1958;19:96779Leger N, Saratsiotis A, Pesson B, Leger P Leishmaniasis in Greece. Results of an entomological study during June 1977.Ann Parasitol Hum Comp 1979;54:1129485036Leger N, Gramiccia M, Gradoni L, Madulo-Leblond G, Pesson B, Ferte H, Isolation and typing of Leishmania infantum from Phlebotomus neglectus on the island of Corfu, Greece.Trans R Soc Trop Med Hyg 1988;82:41920 10.1016/0035-9203(88)90145-93232178Leger N, Pesson B, Madulo-Leblond G, Ferte H, Tselentis I, Antoniou M The phlebotomes of Crete[in French] Biologia Gallo-Hellenica.1993;20:13543Ivovic V, Patakakis M, Tselentis Y, Chaniotis B Faunistic study of sandflies in Greece.Med Vet Entomol 2007;21:1214 10.1111/j.1365-2915.2006.00649.x17373955Williams BG, Dye C Maximum likelihood for parasitologists.Parasitol Today 1994;10:48993 10.1016/0169-4758(94)90163-515275519Sutton AJ, Gay NJ, Edmunds WJ, Hope VD, Gill ON, Hickman M Modelling the force of infection for hepatitis B and hepatitis C in injecting drug users in England and Wales.BMC Infect Dis 2006;6:93 10.1186/1471-2334-6-9316762050Courtenay O, MacDonald DW, Lainson R, Shaw JJ, Dye C Epidemiology of canine leishmaniasis: a comparative serological study of dogs and foxes in Amazon Brazil.Parasitology 1994;109:2739 10.1017/S00311820000783067970884Quinnell RJ, Courtenay O, Garcez L, Dye C The epidemiology of canine leishmaniasis: transmission rates estimated from a cohort study in Amazonian Brazil.Parasitology 1997;115:14356 10.1017/S003118209700120010190170Courtenay O, Quinnell RJ, Garcez LM, Shaw JJ, Dye C Infectiousness in a cohort of Brazilian dogs: why culling fails to control visceral leishmaniasis in areas of high transmission.J Infect Dis 2002;186:131420 10.1086/34431212402201