Figures
Abstract
Background
Visceral Leishmaniasis (VL) is a life threatening neglected infectious disease in the Indian subcontinent, transmitted by the bite of female sand flies. Estimation of the infectivity in the vector population, collected in different seasons, may be useful to better understanding the transmission dynamics of VL as well as to plan vector control measures.
Methodology
We collected sand flies from highly endemic regions of Bihar state, India for one year over three seasons. The species of the sand flies were confirmed by species-specific PCR-RFLP. Leishmania donovani infection was investigated in 1397 female Phlebotomus argentipes using PCR, targeting the Leishmania specific minicircle of the kDNA region. Further, the parasitic load in the infected sand flies was measured using quantitative PCR.
Conclusion
Though sand flies were most abundant in the rainy season, the highest rate of infection was detected in the winter season with 2.84% sand flies infected followed by the summer and rainy seasons respectively. This study can help in vector elimination programmes and to reduce disease transmission.
Citation: Tiwary P, Kumar D, Mishra M, Singh RP, Rai M, Sundar S (2013) Seasonal Variation in the Prevalence of Sand Flies Infected with Leishmania donovani. PLoS ONE 8(4): e61370. https://doi.org/10.1371/journal.pone.0061370
Editor: David Joseph Diemert, The George Washington University Medical Center, United States of America
Received: August 3, 2012; Accepted: March 11, 2013; Published: April 9, 2013
Copyright: © 2013 Tiwary et al. 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 author and source are credited.
Funding: The study was funded by National Institute of Allergy and Infectious Disease (NIAID), DMID funding mechanism: Tropical Medicine Research Center Grant number: P50AI074321. Puja Tiwary and Dinesh Kumar thank CSIR for providing fellowships. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
In vector borne diseases, the seasonality of the vector species directly affects the transmission of the infectious disease. Thus, the study of insect behaviour in natural settings is crucial to understanding the current status of any disease. Leishmaniases are some of the most severe forms of vector borne infectious diseases, affecting 98 countries worldwide. In India, visceral leishmaniasis (VL) is a major cause for morbidity and mortality. It is prevalent in the four states of Bihar, Jharkhand, West Bengal and Uttar Pradesh, with about 84% of all VL cases being reported from Bihar [1].
Sand flies, the vectors for leishmaniasis, belong to the genus Phlebotomus in the old world and Lutzomyia in the new world. In India, VL is caused by infection with the protozoan parasite Leishmania donovani with female Phlebotomus argentipes considered as the vector [2]. The female sand fly needs blood for its eggs to develop, and becomes infected with the Leishmania parasite while biting an infected individual or animal. Over a period of 6 to 9 days, the parasite develops inside the sand fly gut [3]. When an infected female sand fly feeds on an uninfected person or animal, it simultaneously inoculates the parasite, thus completing the transmission cycle. The incubation period from infection to clinical VL is between 2 to 6 months [4].
Data on the population dynamics of P. argentipes are scanty from the Indian subcontinent. Changes in climate affect the density and behavioural bionomics of a vector species [5], [6], [7], [8] and thus disease transmission. In India, three seasons occur throughout the year, each with its distinct temperature and humidity. Furthermore, seasonal fluctuations in the pattern of VL transmission have been reported [9] that indicate seasonal variations in rates of L. donovani infection in the female P. argentipes vector that may in turn have direct impact on VL transmission throughout the year. In this study, seasonal variations in the P. argentipes population, the infection rate and the parasitic load were investigated employing molecular methods.
Materials and Methods
Study area
Sand flies were collected from VL endemic regions of Muzaffarpur district, Bihar, located at 26.07°N 85.45°E that lies in a highly active seismic zone of India. For the sampling, villages were screened on the basis of reported VL cases. Relevant information concerning VL cases was obtained from the Kala Azar Medical Research Centre, Muzaffarpur. Using these data, 50 villages were randomly chosen and 10 houses were selected from each village amounting to a total of 500 houses (Figure 1).
Shaded areas represent villages selected for sand fly collection. Distribution of infected sand flies in the summer (hollow arrow), rainy (solid arrow) and winter (star) seasons, respectively.
Approval for the study was obtained from the Ethical Committee of the Institute of Medical Sciences, Banaras Hindu University. Written informed consent was obtained from the head of each household.
Collection of sand flies
To know about the actual status of sand fly infection, collections were conducted during the three different seasons over the course of one year. A total of 500 selected houses were monitored three times during 2010, in the summer (March-April), rainy (August-September) and winter (November-December) seasons. CDC light traps with standard collection cup assemblies (Model 1012 miniature Incandescent Light Trap, John W. Hock, USA) were used for sand fly collection. One light trap was installed inside each house in the evening and left overnight as sand flies are nocturnally active [5], [10]. A total of 10 light traps were installed in 10 houses of each village in one day. In each season, collections were done for 50 days to cover all the villages. Digital thermometers and hygrometers were also installed inside each house to record maximum and minimum temperatures as well as humidity, respectively. On the following day, early in the morning, collected insects were stored in petri dishes containing chloroform coated cotton balls. Sand flies that were resting on the walls of houses were collected using mouth aspirators.
Morphological identification
Morphological identification was done in two steps. In the first step, sand flies were separated from other insects using an optical microscope while in the second step, sand flies were identified using a stereomicroscope (Carl Zeiss, Singapore) on the basis of morphological characteristics such as the thorax and hairs of abdominal tergites [11]. Collected sand flies belonged to three species: P. argentipes, P. papatasi and Sergentomyia babu. Using morphological characteristics of the anal region, these were further divided into male and female. All the sand flies were stored individually in 70% ethanol at −20°C. Only female P. argentipes were selected for this study as they are the only known vector for L. donovani infection.
Sand fly DNA extraction
Each sand fly was individually processed for DNA isolation using a High Pure PCR Template Preparation Kit (Roche, Mannheim, Germany) after crushing the whole sand fly with a sterile polystyrene probe. DNA was quantitated using a NanoDrop 2000 (Thermo Scientific Inc., Wilmington, DE).
Parasite DNA extraction
Leishmania parasite reference strain LEM138 (MHOM/IN/00/DEVI) was cultured in M199 at 26°C and DNA was extracted using DNeasy Blood & Tissue Kit (Qiagen, USA).
Identification of sand fly by PCR-RFLP
For molecular confirmation of morphologically identified species using PCR, primers were designed targeting the 18s rRNA coding gene: 5′-TAGTGAAACCGCAAAAGGCTCAG-3′ (forward) and 5′-CTCGGATGTGAGTCCTGTATTGT-3′ (reverse) [12]. The PCR reaction was carried out in a volume of 25 µl using the pair of primers (10 pmol each), 1 U normal Taq DNA polymerase (New England Biolabs, United Kingdom) supplemented with MgCl2, 10× buffer, 1× BSA, 0.2 mM of dNTPs. For a template, 50 ng of extracted DNA was used while nuclease-free water (QIAGEN, Hilden, Germany) was used as a negative control. PCR conditions were: initial denaturation at 95°C for 5 min, followed by 35 cycles of denaturation at 95°C for 30 sec, annealing at 58°C for 40 sec, extension at 72°C for 30 sec, and a final extension at 72°C for 10 min. Further, restriction digestion was done to discriminate P. argentipes from other morphologically similar sand flies. The amplified PCR product was subjected to restriction digestion using Hinf I restriction enzyme (recognition site GATTC; time saver, New England Bio labs, UK) at 37°C for 15 min. The entire reaction product was electrophoresed on a 1.5% (w/v) ethidium bromide (Merck, Darmastaft, Germany) stained agarose gel (Sigma-Aldrich, USA) (Figure 2).
18S rRNA gene specific PCR products were digested using the Hinf I restriction enzyme. Lane M: 100 bp marker, Lane 1: P. argentipes, Lane 2: P. papatasi, Lane 3: S. babu.
Leishmania specific PCR
Individual P. argentipes confirmed by both morphological and molecular methods were used to study the natural infection rates with Leishmania using primers specific for the minicircle of the kDNA region of the parasite [13]. The primers used were 5′- CCTATTTTACACCAACCCCCAGT-3′ (forward) and 5′- GGGTAGGGGCGTTCTGCGAAA-3′ (reverse). The DNA of L. donovani reference strain, LEM138 (MHOM/IN/00/DEVI), was used as a positive control whereas DNA extracted from male sand flies was used as a negative control. The reaction was carried out in a volume of 25 µL using the primers (10 pM each), 2 U of Taq DNA polymerase supplemented with 10× buffer, MgCl2, 1× BSA, 25 µM dNTPs and 100 ng of template. The PCR reaction conditions included: initial denaturation at 95°C for 5 min, followed by 40 cycles of denaturation at 95°C for 30 sec, annealing at 58°C for 30 sec, extension at 72°C for 30 sec, and a final extension at 72°C for 10 min. Amplified products were resolved on 1.5% ethidium bromide stained agarose gel (Figure 3).
Lane 1: positive control, Lanes 2–13: female sand flies showing amplification in lanes 3, 8, 9 and 10, Lane 14: male sand fly.
Sequencing
Both P. argentipes and L. donovani specific PCR products were directly used for sequencing. The P. argentipes specific PCR product was observed as a single band on the agarose gel, thus leftover primers and dNTPs from PCR products were removed using QIAquick PCR Purification Kit (Qiagen). For L. donovani specific PCR reactions, non-specific products were visible, thus the desired product bands were cut and extracted using QIAquick Gel Extraction Kit (Qiagen). Further, the eluted PCR products were sequenced by BigDye Terminator 3.1 cycle sequencing kit (Applied Biosystems, Foster City, CA, USA) and run through an ABI 3133 sequencer.
Quantitative real-time PCR
Leishmania specific minicircle kDNA primers as described above were used to quantitate parasite load in the PCR-positive sand flies [13]. Determination of parasite load was done using a standard curve made with 10-fold serial dilutions of Leishmania DNA spiked with a fixed amount of sand fly DNA. Parasites were counted using a Neubauer chamber and DNA was eluted with final concentration of 2.5×105 parasites/µl which was diluted 8 times and a standard curve was plotted. A standard curve was obtained by plotting the Ct against the known target quantity on a common log scale (Figure 4). Parasitic load in the infected sand flies was calculated by comparing Ct values with the standard plot (Table 1). Along with the PCR positive samples, PCR negative samples, male sand fly samples and no template master mix were used as negative controls on each plate.
Ct represents the fractional cycle number reflecting a positive PCR result differentiated from background.
For the P. argentipes specific SYBR Green assay, the target was the 18S rRNA encoding gene of the sand fly (GenBank accession number AJ244360.1). For amplification, primers were designed as 5′-TGTGTGCGTTCACTGTCAAAGGTG-3′ (forward) and 5′- AAGTAACTGTACCGGCCCACA-3′ (reverse) to amplify a 188-bp fragment to check for inter-sample variation in sand fly DNA amplification.
The parasitic load of all Leishmania specific PCR-positive sand flies as well as 150 PCR-negative sand flies (50 from each season randomly selected) was calculated. Q-PCR reactions were set up with 10 µl reaction mixture containing 5 µl of 2× SYBR Green reagent (Applied Biosystems), 5 pM each of forward and reverse primers, and 100 ng of each DNA sample. Reactions were analyzed in MicroAmp optical 96 well reaction plates (Applied Biosystems) using a 7500 Real-Time PCR System (Applied Biosystems). The thermal cycling conditions were similar for both L. donovani and P. argentipes specific Q-PCR: initial denaturation stage at 95°C for 8 min, amplification stage (40 cycles at 95°C for 15 sec and 60°C for 1 min each), melting stage (95°C for 15 sec, 60°C for 1 min and 95°C for 30 sec) followed by cooling at 60°C for 15 sec. Each sample was tested in duplicate.
Data analysis and statistical methods
Temperature and humidity were measured in each house involved in the study. A graph was plotted indicating average temperature and relative humidity corresponding to all the 50 villages for each season (Figure 5). There was marked variation in average temperature and relative humidity between seasons.
Average temperature (Celsius) and relative humidity (%) shown for the summer, rainy and winter seasons.
For Q-PCR, generation of amplification plots, mean values, standard curve analysis and calculation of the melting temperature of each amplicon were done directly using the 7500 software v2.0.1 (Applied Biosystems). The fluorescence threshold limit of the ABI 7500 Real-Time PCR System was set at 0.02 as recommended by the manufacturer.
Results
Collection of female P. argentipes
The total collection of sand flies varied according to season: the number of female P. argentipes were 384 (blood fed 58, unfed 172, gravid 154), 591 (blood fed 59, unfed 378, gravid 154) and 422 (blood fed 54, unfed 208, gravid 160) in th esummer, rainy and winter seasons, respectively. Village-wise variation in the number of collected sand flies is shown in Figure 6. All of the morphologically identified P. argentipes were validated using species-specific PCR-RFLP. Sixty randomly chosen P. argentipes specific PCR products (20 from each season) were sequenced using primers specific to the sand fly 18S rRNA coding gene. Sequence analysis confirmed that all of the samples were similar to known P. argentipes 18s rRNA sequences (Genbank Accession number JX910368).
Total collected sand flies (solid hologram) and infected sand flies (star shaped) from each village, shown for the (A) summer, (B) rainy and (C) winter seasons.
L. donovani infection in sand flies
After amplification of the Leishmania-specific minicircle kDNA region, a 116 bp product was observed in the infected female P. argentipes samples. The amplified PCR products were eluted from the agarose gel and sequenced using the same primers. Sequence analysis confirmed that all of the amplified bands had 100% homology with L. donovani minicircle kDNA.
Using L. donovani specific PCR, seasonalal variation in the infection level of sand flies was observed. The positivity rate for the entire year was 1.50% (21/1397). However, the infection rates were 1.04% (4/384) in summer, 0.85% (5/591) in the rainy season, and 2.84% (12/422) in winter. Female sand flies belonged to three categories – blood fed, unfed and gravid – that showed variation in rates of infection. Amongst the infected sand flies, the highest rate of infection was observed in the blood fed group at 2.34% (4/171), followed by the unfed group at 1.58% (12/758), and the gravid group at 1.07% (5/468). The village-wise distribution of infected sand flies shown in Figure 1 and Figure 6.
Parasitic loads of infected sand flies
Parasitic load was investigated by real-time PCR using the standard curve absolute quantification method to estimate the number of Leishmania parasites in naturally-infected sand flies. Parasitic loads of infected sand flies are shown in Table 1. Parasitic loads were highest in gravid sand flies, followed by unfed and blood fed ones.
Discussion
Although VL is considered a significant public health problem in the Indian subcontinent, very little research has been done on the behaviour of sand fly vectors. Correct identification of the vector population and the rate of Leishmania infection in endemic areas is fundamental for understanding the epidemiology of the disease that may help inform vector control programmes.
PCR based methods were used for the study as we needed a highly robust method to include a large sample size. For sand fly identification, all the collected insects were first morphologically identified and differentiated in various groups. For the present study, only female P. argentipes were selected as sexual differentiation is not possible by molecular methods. Further PCR and restriction digestion were done to confirm the sand fly species. Using molecular methods, various sand fly species that co-exist in same area can be easily differentiated. In endemic areas of Muzaffarpur, there are three prevalent species of sand fly: P. argentipes, P. papatasi and S. babu. With restriction digestion, each species displayed an easily identifiable specific band pattern on agarose gel.
Morphological identification of Leishmania parasites by microscopy in sand flies collected in the field may take a long time, as parasite load in naturally infected sand flies could be very low [14]. For highly sensitive Leishmania parasite detection, the minicircle kDNA region was selected as it is present as approximately 10,000 copies per parasite [15]. Two sets of primers were used for the identification of P. argentipes and L. donovani, respectively. Detection of parasite infections in the insects using morphological methods can be time consuming, as consistent observation under a microscope is needed for large sample sizes and samples cannot be used further as there may be a risk of contamination after dissection. PCR based methods targeting various conserved genes such as minicircle kDNA [16], [17], ITS [18], [19], Cytochrome b gene [20] and 18S rRNA [21] encoding regions have been used for detection of the Leishmania parasite in sand fly vectors. As previously reported, since the parasitic load in naturally-infected sand flies is less, a highly sensitive primer was needed. PCR is preferred over microscopy as a large sample size can be processed simultaneously and it is more sensitive for detecting infection. Also, PCR is time-saving and has been used for similar studies in other countries, allowing for comparability of results [20], [22].
In studies of the prevalence of Leishmania infection in sand flies, collection and processing is usually done by two methods, either by pooling sand flies or by testing them individually. The pooling method is easier, time-saving and efficient [21], [22], [23], [24]. However, the exact number of infected sand flies in a positive pool is impossible to calculate and information on individual sand flies, such as parasitic load or blood meal host, can be lost. Such drawbacks can be overcome by assaying individual sand flies. Previous studies have shown prevalences of Leishmania infection in natural sand fly populations ranging from 0.7 to 2.0% [19], [25], [26], [27]. Considering the limitations inherent in pooling samples, as well as the low expected infection rates, each sand fly was processed individually to determine the exact infection rate in our study.
Based on the seasonal collection of sand flies throughout the year, we found that sand flies are highly prevalent in rainy season. The aim of the present study was to describe the seasonal variation in L. donovani infection in the vector population. The overall percentage of L. donovani infection in female P. argentipes was 1.5% although this varied by season. The reason that the lowest rate of infection was observed in the rainy season might be due to the flushing effect of rainfall on immature, emerging sand flies. With the drainage of water during winters, immature sand flies start developing on the ground and thus during winter, more numbers of mature and potentially infective sand flies exist [28]. The highest parasitic load in summer might be due to the occurrence of the highest number of gravid sand flies as they need to ingest a greater quantity of blood for their eggs to develop. At the same time, there are more active VL patients during summer, so the flies may ingest a higher number of parasites while taking a blood meal from an infected human host.
This observation may be beneficial for understanding VL dynamics as behavioural changes in vectors can have direct impact on disease transmission. Vector control is still the most powerful tool for preventing vector borne diseases. Thus, knowledge of the infection status of the vector population has proved crucial for disease control programmes. With the present study, maximum infection was observed in sand flies collected in the months of November-December with 2.84% infection; consequently, the infectivity of sand flies should be maximum during these months. The time interval between the bite of an infected sand fly and development of disease is between 2 to 6 months. This finding is supported by another study in which monthly variations in the incidence of VL were observed, with the highest incidence occurring in March-April [9]. The present study may be beneficial for vector control programmes in targeting the most effective period for insecticidal spraying.
Author Contributions
Conceived and designed the experiments: SS MR. Performed the experiments: PT RPS MM DK. Analyzed the data: SS MR PT RPS MM DK. Contributed reagents/materials/analysis tools: SS MR PT. Wrote the paper: SS MR PT RPS MM DK.
References
- 1. WHO (2010) Control of the leishmaniases. Report of a WHO Expert Committee. World Health Organ Tech Rep Ser 949: 1–200.
- 2. Swaminath CS, Shortt HE, Anderson LA (2006) Transmission of Indian kala-azar to man by the bites of Phlebotomus argentipes, ann and brun. 1942. Indian J Med Res 123: 473–477.
- 3. Kamhawi S (2006) Phlebotomine sand flies and Leishmania parasites: friends or foes? Trends Parasitol 22: 439–445.
- 4. Chappuis F, Sundar S, Hailu A, Ghalib H, Rijal S, et al. (2007) Visceral leishmaniasis: what are the needs for diagnosis, treatment and control? Nat Rev Microbiol 5: 873–882.
- 5. Dinesh DS, Ranjan A, Palit A, Kishore K, Kar SK (2001) Seasonal and nocturnal landing/biting behaviour of Phlebotomus argentipes (Diptera: Psychodidae). Ann Trop Med Parasitol 95: 197–202.
- 6. Gálvez R, Descalzo MA, Miró G, Jiménez MI, Martín O, et al. (2010) Seasonal trends and spatial relations between environmental/meteorological factors and leishmaniosis sand fly vector abundances in Central Spain. Acta Trop 115: 95–102.
- 7. Boussaa S, Guernaoui S, Pesson B, Boumezzough A (2005) Seasonal fluctuations of phlebotomine sand fly populations (Diptera: Psychodidae) in the urban area of Marrakech, Morocco. Acta Trop 95: 86–91.
- 8. Souza NA, Andrade-Coelho CA, Vilela ML, Peixoto AA, Rangel EF (2002) Seasonality of Lutzomyia intermedia and Lutzomyia whitmani (Diptera: Psychodidae: Phlebotominae), occurring sympatrically in area of cutaneous leishmaniasis in the State of Rio de Janeiro, Brazil. Mem Inst Oswaldo Cruz 97: 759–765.
- 9. Malaviya P, Picado A, Singh SP, Hasker E, Singh RP, et al. (2011) Visceral leishmaniasis in Muzaffarpur district, Bihar, India from 1990 to 2008. PLoS One 6: e14751.
- 10. Morrison AC, Ferro C, Pardo R, Torres M, Wilson ML, et al. (1995) Nocturnal activity patterns of Lutzomyia longipalpis (Diptera: Psychodidae) at an endemic focus of visceral leishmaniasis in Colombia. J Med Entomol 32: 605–617.
- 11. Lewis DJ (1982) A taxonomic review of the genus Phlebotomus. Bull Brit Mus Nat Hist Entomol 45: 121–209.
- 12. Tiwary P, Kumar D, Rai M, Sundar S (2012) PCR-RFLP based method for molecular differentiation of sand fly species Phlebotomus argentipes, Phlebotomus papatasi, and Sergentomyia babu found in India. J Med Entomol 49: 1515–1518.
- 13. Nicolas L, Prina E, Lang T, Milon G (2002) Real-time PCR for detection and quantitation of leishmania in mouse tissues. J Clin Microbiol 40: 1666–1669.
- 14. Killick Kendrick R (1990) Phlebotomine vectors of the leishmaniases: a review. Med Vet Entomol 4: 1–24.
- 15. Brewster S, Aslett M, Barker DC (1998) Kinetoplast DNA minicircle database. Parasitol Today 14: 437–438.
- 16. Aransay AM, Scoulica E, Tselentis Y (2000) Detection and identification of Leishmania DNA within naturally infected sand flies by seminested PCR on minicircle kinetoplastic DNA. Appl Environ Microbiol 66: 1933–1938.
- 17. Barker DC, Gibson LJ, Kennedy WP, Nasser AA, Williams RH (1986) The potential of using recombinant DNA species-specific probes for the identification of tropical Leishmania. Parasitology 92 Suppl: S139–174
- 18. Azpurua J, De La Cruz D, Valderama A, Windsor D (2010) Lutzomyia sand fly diversity and rates of infection by Wolbachia and an exotic Leishmania species on Barro Colorado Island, Panama. PLoS Negl Trop Dis 4: e627.
- 19. Svobodova M, Alten B, Zidkova L, Dvorak V, Hlavackova J, et al. (2009) Cutaneous leishmaniasis caused by Leishmania infantum transmitted by Phlebotomus tobbi. Int J Parasitol 39: 251–256.
- 20. Kato H, Uezato H, Katakura K, Calvopina M, Marco JD, et al. (2005) Detection and identification of Leishmania species within naturally infected sand flies in the andean areas of ecuador by a polymerase chain reaction. Am J Trop Med Hyg 72: 87–93.
- 21. Bhattarai NR, Das ML, Rijal S, van der Auwera G, Picado A, et al. (2009) Natural infection of Phlebotomus argentipes with Leishmania and other trypanosomatids in a visceral leishmaniasis endemic region of Nepal. Trans R Soc Trop Med Hyg 103: 1087–1092.
- 22. Anderson JM, Samake S, Jaramillo-Gutierrez G, Sissoko I, Coulibaly CA, et al. Seasonality and prevalence of Leishmania major infection in Phlebotomus duboscqi Neveu-Lemaire from two neighboring villages in central Mali. PLoS Negl Trop Dis 5: e1139.
- 23. Dinesh DS, Kar SK, Kishore K, Palit A, Verma N, et al. (2000) Screening sandflies for natural infection with Leishmania donovani, using a non-radioactive probe based on the total DNA of the parasite. Ann Trop Med Parasitol 94: 447–451.
- 24. Tiwary P, Kumar D, Singh RP, Rai M, Sundar S (2012) Prevalence of sand flies and Leishmania donovani infection in a natural population of female Phlebotomus argentipes in Bihar State, India. Vector Borne Zoonotic Dis 12: 467–472.
- 25. Paiva BR, Oliveira AG, Dorval ME, Galati EA, Malafronte RS (2010) Species-specific identification of Leishmania in naturally infected sand flies captured in Mato Grosso do Sul State, Brazil. Acta Trop 115: 126–130.
- 26. Rocha LS, Falqueto A, dos Santos CB, Ferreira AL, da Graca GC, et al. (2010) Survey of natural infection by Leishmania in sand fly species collected in southeastern Brazil. Trans R Soc Trop Med Hyg 104: 461–466.
- 27. Sanchez-Garcia L, Berzunza-Cruz M, Becker-Fauser I, Rebollar-Tellez EA (2010) Sand flies naturally infected by Leishmania (L.) mexicana in the peri-urban area of Chetumal city, Quintana Roo, Mexico. Trans R Soc Trop Med Hyg 104: 406–411.
- 28.
Mukhopadhyay AK, Rahaman SJ, Chakravarty AK (1990) Effect of flood on Immature stages of Sandflies in a Flood Prone Kala-azar Endemic Village of North Bihar, India. WHO/VBC/90.986.Geneva: World Health Organization.