Phlebotomus (Paraphlebotomus) sergenti is at least one of the confirmed vectors for the transmission of cutaneous leishmaniasis caused by Leishmania tropica and distributed widely in Morocco. This form of leishmaniasis is considered largely as anthroponotic, although dogs were found infected with Leishmania tropica, suggestive of zoonosis in some rural areas. This survey aimed at (i) studying the presence of Leishmania in field caught Phlebotomus sergenti, (ii) investigating genetic diversity within Leishmania tropica and (iii) identifying the host-blood feeding preferences of Phlebotomus sergenti. A total of 4,407 sand flies were collected in three rural areas of Azilal province, using CDC miniature light traps. Samples collected were found to consist of 13 species: Phlebotomus spp. and 3 Sergentomyia spp. The most abundant species was Phlebotomus sergenti, accounting for 45. 75 % of the total. 965 female Phlebotomus sergenti were screened for the presence of Leishmania by ITS1-PCR-RFLP, giving a positive rate of 5. 7% (55/965), all being identified as Leishmania tropica. Nucleotide heterogeneity of PCR-amplified ITS1-5. 8S rRNA gene-ITS2 was noted. Analyses of 31 sequences obtained segregated them into 16 haplotypes, of which 7 contain superimposed peaks at certain nucleotide positions, suggestive of heterozygosity. Phlebotomus sergenti collected were found to feed on a large variety of vertebrate hosts, as determined by Cytochrome b sequencing of the DNA from the blood meals of 64 engorged females. Our findings supported the notion that Phlebotomus sergenti is the primary vector of Leishmania tropica in this focus, and that the latter is genetically very heterogeneous. Furthermore, our results might be suggestive of a certain level of heterozygosity in Leishmania tropica population. This finding, as well as the feeding of the vectors on different animals are of interest for further investigation.
Leishmaniases are complex diseases of worldwide distribution caused by >20 Leishmania species, which are parasitic protozoa transmitted by the bite of infected female sand flies. The disease affects 98 Mediterranean and other endemic countries putting a population of 350 million people at risk of infection and causing 1. 3 million new cases and 20,000 to 30,000 deaths annually [1,2]. In Morocco, leishmaniasis is a public health problem, with 8,862 cases notified in 2010 alone. It is widely distributed from the mountains of the Riff to the peri-arid foothills of the Anti-Atlas [3,4]. Three Leishmania species are responsible for cutaneous leishmaniasis (CL) in Morocco: L. major, L. tropica and L. infantum. CL due to L. tropica is endemic in the semi-arid regions and the sand fly species Phlebotomus (Paraphlebotomus) sergenti Parrot is considered to be the vector [5,6]. This form of CL was described for the first time in Tanant, a rural locality in Azilal province in High Atlas [7]. Thereafter, a large rural focus has been further identified in the centre and the south of Morocco [6]. In the mid 1990s, CL caused by L. tropica was then reported in Taza province, in northern Morocco [8] and early in 2000, outbreaks occurred in emerging foci in the central and northern parts of the country, and L. tropica was also even found in areas previously only known as L. major foci [9,10]. CL due to L. tropica is basically anthroponotic, although the parasite was also isolated from dogs in some leishmaniasis foci in Morocco, suggesting that L. tropica could also be zoonotic at least in some areas [5,11,12]. Blood meal analysis of haematophagous insects has been informative in elucidating their natural host-feeding preference or host-feeding pattern for identifying potential reservoirs [13]. For sand flies, the blood meal sources have been identified initially by serological techniques, like ELISA [14,15], agarose gel diffusion [16], counter immunoelectrophoresis [17], precipitin test [18,19] and a more laborious histological technique [20]. Although all these methods have been useful in identifying vertebrate hosts for many haematophagous insects, they lack sensitivity and are time-consuming. In the last decade, several PCR-based molecular approaches have been developed to enhance the specificity of insect blood meal identification [21,22]. Analysis of blood meal sources of the engorged vectors has been greatly facilitated by the increasing availability of animal genomic database at the genus and species levels (NCBI, EBI…). Downstream applications of PCR based on primers designed from multiple alignments of the mitochondrial Cytochrome b gene of different vertebrate species have identified avian and mammalian blood sources in various species of sand flies. The sequencing of Cyt b gene, as well as the Cyt b-PCR-RFLP and the Cyt b-PCR combined with reverse line blot have been used to identify the blood meals sources of sand flies [23–25]. The vector incrimination is classically based on the dissection of freshly caught individual sand flies and subsequent culture of parasites present in the gut [26]. This method needs dissecting expertise and is time consuming, since the Leishmania infection rate in sand flies is usually very low [27]. In recent years, several molecular techniques have been developed to identify Leishmania infection in infected phlebotomine sand flies. Because of excellent sensitivity and specificity, these methods have the potential to challenge the time-consuming isoenzyme technique for timely characterization of Leishmania strains. They also play a prominent role in epidemiological prospective surveillance activities which are critical for planning targeted control measures of leishmaniasis [15,28,29]. In the present study, we used molecular tools (ITS1-PCR-RFLP) to detect and identify Leishmania species within naturally infected sand flies collected from three rural localities of Azilal province in Morocco. We also applied nested-PCR of Leishmania ITS-rDNA genes for sequencing analyses to characterize L. tropica. The blood meals of engorged flies were also analyzed to determine the putative animal reservoirs of L. tropica, using a PCR technique based on the Cyt b gene of vertebrate mitochondrial DNA (mtDNA).
The sand flies were collected in three neighboring rural areas in the province of Azilal, High Atlas of Morocco. The three regions are: Ait Makhlouf II: 31°01’23”N, 6°58’53”O and Guimi: 32°00’12”N, 6°55’03”O, which are located in Beni Hassan sector, and Agmeroul: 31°58’38”N, 6°51’12”O in Tabia sector. These localities are at different altitude from 679 m to 840 m, separated from one another at a distance of 5 to 11 km (Fig. 1). Vegetation is sparse and mainly dominated by cactus, jujube plant, olive, and almond trees. Agriculture remains the primary source of income and is mainly based on the production of wheat, almonds and olives. Livestock include chicken, sheep, cattle and horses and most homes have at least one dog. Houses are mainly of the traditional type, constructed with adobe. The villages are surrounded by mountains, with some groundwater sources used for irrigation. Phlebotomine sand flies were caught during three consecutive nights from June to October 2011, using CDC light traps placed inside the houses and in domestic animal shelters, in the three rural sites. The traps placed about 1. 5 meter above the ground were set before sunset and collected the next morning. The collected sand flies were then placed in 1. 5 mL Eppendorf tubes, transferred to the lab in dry ice and kept at -80°C. All the sand flies were washed and dissected. The head and genitalia were used for morphological identification using Moroccan morphologic key [4]. To distinguish P. perniciosus and P. longicuspis a very closely related species, we used the morphological criteria updated by Benabdennbi et al. (1999) [30] and Guernaoui et al. (2005) [31] based on the morphological features of genitalia and the number of coxite hairs. P. perniciosus is characterized by 2 aedeagi forms: (i) copulatory valve with apex bifid; and (ii) copulatory valve with curved apex, the number of coxite hairs being up to 14. P. longicuspis has copulatory valve ending with a single and long point, and 19 or more coxite hairs. The remainder of the body of engorged and unengorged females after dissection was stored sterilely in 1. 5 mL microtube at -20°C until use. Genomic DNAs were extracted from unengorged and engorged female P. sergenti sand fly by using the Kit PureLink™ Genomic DNA Mini Kit and the QIAamp® DNA Blood Mini Kit, respectively, according to the manufacturer’s instruction. For P. sergenti females, the leishmanial ribosomal internal transcribed spacer 1 (ITS1) region was amplified, using the primers LITSR (5′-CTGGATCATTTTCCGATG-3′) and L5. 8S (5′-TGATACCACTTATCGCACTT-3′), following the protocol described by Schonian et al. (2003). This PCR was used to amplify a 300 to 350 bp fragment; then the ITS1-PCR products were digested with restriction endonuclease HaeIII enzyme for Leishmania species identification [32]. The restriction profiles were analyzed by electrophoresis on 3% agarose gel containing ethidium bromide, and compared with Leishmania reference strains: L. major (MHOM/SU/73/5ASKH), L. tropica (MHOM/SU/74/K27) and L. infantum (MHOM/TN/80/IPT1). A 100 bp DNA size marker was used. Negative controls (containing water, without DNA) were added to each PCR run. Further identification of Leishmania parasites was done using nested PCR amplification of ITS1-5. 8S rRNA gene-ITS2. The primers IR1 and IR2 were used for the first PCR stage and ITS1F and ITS2R4 for the second stage, as previously described by Parvizi and Ready [33]. Negative controls were included to each PCR run. All PCR products were analyzed by electrophoresis on 1% agarose gel containing ethidium bromide. We used the standard DNA fragment 100 bp ladder as a size marker. The nested-PCR products (~460 bp including primers) were sequenced, using the primers ITS1F and ITS2R4. The sequences obtained were edited using Seq Scape, BioEdit softwares and aligned using ClustalW in MEGA6 software. A phylogenetic tree was constructed by using the Maximum Likelihood method based on the Tamura-Nei model. Initial tree (s) for the heuristic search were obtained automatically by applying the Maximum Parsimony method. Bootstrap replicates were performed to estimate the node reliability, and values were obtained from 1,000 randomly selected samples of the aligned sequence data. Sequences were compared with entries retrieved from GenBank. Evolutionary analyses were conducted in MEGA6. Morphologically identified P. sergenti females were confirmed by sequence analyses of PCR-amplified mitochondrial DNA fragment encompassing Cyt b region using PhleF (5’-AAT AAA TTA GGA GTA ATT GC-3’) and PhleR (5’-GCC TCG AWT TCG WTT ATG ATA AAT T-3’) primers [34]. A 500 bp fragment was PCR-amplified using the following reagents: 1X PCR buffer, 2. 5 mM MgCl2,0. 4 mM primers (F and R), 0. 2 mM dNTPs, 4 μL DNA, made up with distilled water to a final volume of 25 μL for each reaction. The PCR runs were each initiated as follows: 94°C for 12 min, followed by 5 cycles of 94°C for 30 s; 52°C for 30 s, 72°C for 1 min. The PCR reaction was then continued for 30 cycles under the same conditions, except for the annealing step (54°C for 1 min) and a final extension at 72°C for 10 min. The PCR products were sequenced, employing the same primers used for the PCR. Sequences were processed and aligned, using the multiple alignment programs Seq Scape and BioEdit. A phylogenetic tree was constructed by using the Neighbor-Joining method, in agreement with Kimura 2-parameter model, at uniform rate for transitions and transversions. Bootstrap replicates were performed to estimate the node reliability, and values were obtained from 1,000 randomly selected samples of the aligned sequence data. Sequences were compared with entries retrieved from GenBank. The phylogenetic tree was constructed by using the MEGA6 software. A region of Cyt b gene from host mtDNA was amplified for blood meal source identification of the blood-fed female specimens. The assay is based on specific amplification and sequencing of the blood meal–derived the region of the host mtDNA Cyt b gene. Modified vertebrate-universal specific primers L14841 (5’-CCATCCAACATCTCAGCATGATGAAA-3’) and H15149 (5’-CCCCTCAGAATGATATTTGTCCTCA-3’) were used to amplify a 359 bp segment of the Cyt b gene [35,36]. PCR reaction mixture contained 1x PCR Rxn Buffer, 3 mM MgCl2,300 μM deoxynucleotide triphosphates, 0. 4 μM of each primer, 1 unit of Taq DNA polymerase, 3 μL of DNA template solution and distilled water in a final volume of 25 μL. Samples were incubated at 95°C for 5 min; followed by 36 cycles each at 95°C for 30 s, 60°C for 50 s, and 72°C for 40 s; and a final extension at 72°C for 5 min. A negative control was included for each run. Amplified products were analyzed by electrophoresis on 1% agarose gels stained with 2 mg/mL ethidium bromide and visualized under ultraviolet (UV) light. The amplified products were sequenced on both strands, using a Big Dye Terminator v3. 1. Sequences were edited and aligned using Seq-Scape and BioEdit softwares. The sequences obtained were aligned with those deposited in the GenBank database through the BLAST program for the identification of blood meal vertebrate sources. Sequences of a given pairwise alignment with the lowest E-value were selected as the most likely species of blood meal vertebrate source. In order to characterize the species diversity of sand fly populations, several parameters were used: (i) species richness (S) corresponding to the number of species present in the studied habitat; (ii) relative abundance (pi) reflecting the proportion of individuals belonging to the species i. These parameters were used to calculate different ecological indexes [37,38]. Simpson’s diversity index (D) that quantifies the biodiversity of a habitat. It takes into account the number of species present as well as the abundance of each species. D=1∑i=1Spi2 Shannon index (H’) that quantifies species diversity: H′=−∑i=1spilnpi Pielou index (J’): to calculate the equitability of the distribution of species. J′ = H′/ln (S) where H’ is the value of Shannon index and S is the species richness.
A total of 4,407 sand flies (2,007 females and 2,400 males) were collected in the three rural areas of Azilal province. Morphological analysis identified 10 Phlebotomus species of three subgenera as follows: P. sergenti, P. kazeruni and P. chabaudi of the subgenus Paraphlebotomus; P. perniciosus, P. longicuspis, P. perfiliewi, P. ariasi, P. langeroni and P. chadlii of the subgenus Larroussius and P. papatasi of the subgenus Phlebotomus. Three Sergentomyia species were also identified: Sergentomyia minuta, S. antennata and S. fallax. The predominant species was P. sergenti; with a relative abundance of 47. 32% of the sand flies belonging to the genus Phlebotomus and 45. 75% of the total sample. P. perniciosus and P. longicuspis are the next most abundant species, representing 26. 82% and 14. 23% of the sample respectively (Table 1). July had the highest sand fly abundance, whereas the greatest specific richness (S) was found in September, with the occurrence of eleven species against ten, nine and five species recorded in June, July and October, respectively. The three biodiversity indices were maximal in September (Table 2). All the 39 L. tropica infected female P. sergenti that were identified morphologically, were validated by molecular characteristics based on the sequencing of Cyt b fragment. The alignment of the Cyt b sequences obtained confirmed that all samples corresponded to P. sergenti. All sequences were similar to P. sergenti Cyt b sequences: DQ840350, JN036763, DQ840392 and JN036764. Phylogenetic analyses grouped them in a distinct clade with other P. sergenti retrieved from GenBank (Fig. 4). One hundred and twelve full or partially blood-fed females (11. 6% of the total P. sergenti females collected) were tested for their blood feeding preferences. The Cyt b gene sequencing revealed that 64 P. sergenti fed on a variety of vertebrate hosts, including humans, chickens (Gallus gallus), rodents (Myomyscus brockmani, Mesocritus brandti, Dipodomys ordii: 83–91% of identity), birds, cattle (Bos Taurus), rabbit (Oryctolagus cuniculus), bat (Molossus molossus 89% of identity) and monkey (Macaca nigra 95% of identity) (Table 4).
In the present study, we demonstrated a high infection rate of L. tropica in a single Phlebotomus species, that predominantly, but not exclusively, feeds on humans. This confirms the status of P. sergenti as the primary vector of L. tropica in this focus. Our results highlighted the high diversity of L. tropica in Morocco. A zoonotic transmission of L. tropica in this region is further suggested from the blood meal analysis of the vector. The potential for an animal reservoir host for this species awaits further exploration.
In Morocco three Leishmania species have been reported to cause cutaneous leishmaniasis: Leishmania major, Leishmania tropica and less frequently Leishmania infantum. Amongst these clinically important Leishmania species, Leishmania tropica is considered as a public health problem by the Ministry of Health in Morocco and other endemic countries. Phlebotomus sergenti is the known vector, which is thought to take blood meals mainly from humans, since they appear to be the sole reservoir, considering anthroponosis of the cutaneous leishmaniasis caused by L. tropica in many endemic areas. In the present study, we investigated by molecular tools the presence of Leishmania in field caught Phlebotomus, as well as the heterogeneity of Leishmania tropica in a cutaneous leishmaniasis focus in High Atlas of Morocco. Our results showed a high infection rate of Phlebotomus sergenti, which may be a consequence of high level of the parasite circulating in this region; and underlined the important genetic heterogeneity of Leishmania tropica in Morocco. Analysis of the blood meals of the vectors showed that Phlebotomus sergenti fed on a variety of vertebrates, including wild animals, such as rodent, monkey and bat. Whether these animals play any role in the maintenance of Leishmania tropica in this focus awaits further investigation.
Phlebotomus sergenti in a Cutaneous Leishmaniasis Focus in Azilal Province (High Atlas, Morocco): Molecular Detection and Genotyping of Leishmania tropica, and Feeding Behavior