Next Article in Journal
Establishment of a Novel Baculovirus–Silkworm Expression System
Previous Article in Journal
COVID-19 and C. auris: A Case-Control Study from a Tertiary Care Center in Lebanon
Previous Article in Special Issue
The Parasite Load of Trypanosoma cruzi Modulates Feeding and Defecation Patterns of the Chagas Disease Vector Triatoma infestans
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Microorganisms
Volume 10
Issue 5
10.3390/microorganisms10051012
microorganisms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Diversity, Abundance and Leishmania infantum Infection Rate of Phlebotomine Sandflies in an Area with Low Incidence of Visceral Leishmaniasis in Northern Tunisia

by
Marwa Weslati
1,
Jamila Ghrab
1,2,
Meriem Benabid
1,
Olfa Souissi
1,
Karim Aoun
1 and
Aïda Bouratbine
1,*
1
«Laboratoire de Recherche Parasitoses Médicales, Biotechnologies et Biomolécules, LR 16IPT06 and LR 20IPT06», Institut Pasteur de Tunis, Université Tunis El-Manar, Tunis 1002, Tunisia
2
Institut Supérieur des Sciences et Technologies de l’Environnement de Borj Cedria, Université de Carthage, Carthage 1054, Tunisia
*
Author to whom correspondence should be addressed.
Microorganisms 2022, 10(5), 1012; https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms10051012
Submission received: 29 March 2022 / Revised: 29 April 2022 / Accepted: 6 May 2022 / Published: 11 May 2022
(This article belongs to the Special Issue Vector-Parasite Relationships)
Browse Figures
Review Reports Versions Notes

Abstract

:
We report the study of sandfly Leishmania infection in an area of low incidence of visceral leishmaniasis in Tunisia. Sandflies were collected monthly using CDC light-traps set in houses and animal shelters during May–November 2016 and 2017. All males were identified at the species level. A sample of 878 females including all gravid specimens was subjected to kDNA qPCR for Leishmania detection and parasite load estimation. Leishmania species were determined by ITS1 PCR sequencing, and species identification of infected sandflies was performed by DNA barcoding. Phlebotomus perfiliewi and P. perniciosus were the dominant species during the two-year period. However, comparison of their relative abundances showed that P. perniciosus was more abundant during peaks of 2017 with longer activity duration. Real-time kDNA PCR did not detect Leishmania infection in 2016, although it identified four positive specimens (0.7%) in 2017. All four infected specimens were identified as P. perniciosus. ITS1 PCR sequencing allowed L. infantum identification in one kDNA qPCR-positive specimen. This was a P. perniciosus gravid female with a high parasite load caught during the long-lasting peak of 2017. This work highlights the usefulness of multi-seasonal studies of sandfly dynamics and kDNA qPCR in screening Leishmania infection and determining L. infantum vectors in hypo-endemic foci of human leishmaniasis.

1. Introduction

Throughout the Mediterranean basin, Leishmania (L.) infantum is the causative agent of visceral leishmaniasis (VL) and also causes cases of cutaneous leishmaniasis [1]. Domestic dogs are the main reservoir hosts and a dozen of Phlebotomus species are involved in transmission [2,3,4]. Among potential vector species, three, namely, Phlebotomus ariasi, P. perniciosus, and P. perfiliewi, all members of the Larroussius subgenus, have been incriminated as vectors according to conventional criteria in the western Mediterranean region [5,6,7]. In this Mediterranean region, seasonal dynamics of competent vectors showed a regular bell-shaped density curve with a broad central peak encompassing the period July–September and falling between the beginning of May and the end of October [8].
In Tunisia, which is localized in the Southern part of western Mediterranean basin, L. infantum is known to be involved in canine leishmaniasis and both cutaneous and VL in humans [9]. The cutaneous form prevails in the Northern part of the country while VL is distributed in the North and the Center [10,11]. In this large VL endemic area, the spatial study of disease incidence by district identified a hotspot region where VL was homogeneously endemic at a medium/high magnitude. This VL hotspot encompassed the known historic northern VL foci located in the semi-arid zone and supported the presence of a spatial environmental factor underlying the level of endemicity in the area [12]. The entomological survey showed that P. perniciosus was the almost exclusive Larroussius species in the hot spot foci [13]. It was very abundant, especially inside houses, with long seasonal activity and a high L. infantum infection rate [13]. On the other hand, the same spatial study showed that VL incidence rates in districts outside the hotspot were highly variable and not spatially related to each other; the high and low incidence rates in these districts were most likely the result of local factors [12]. This was the case of the emerging VL foci located in the arid bioclimatic zone of Central Tunisia that evidenced a high mobilization of water resources recent years [12,14]. In these emerging foci of high VL incidence, several entomological surveys have been carried out showing L. infantum infection of several Laroussius species, notably P. perniciosus, P. perfiliewi, and P. longicuspis [15,16,17], but also showed L. infantum infection of Phlebotomus papatasi [15,17] and Sergentomyia minuta [15]. These studies have also highlighted the role of P. perfiliewi in these particular foci [15,16]. To our knowledge, no information is available on L. infantum infection rate of phlebotomine sandflies in Northern VL endemic areas of low incidence where P. perniciosus and P. perfiliewi are known to be sympatric species [18].
Quantitative polymerase chain reaction (qPCR) targeting minicircle kinetoplast DNA (kDNA) represents a suitable tool for the diagnosis of Leishmania in vectors [19]. In fact, kDNA qPCR is a sensitive method that allows a detection limit of one flagellate per dissected gut, which legitimates its use as a screening method [20]. Furthermore, kDNA qPCR can determine parasite load [20], with a high burden being correlated with strong evidence of Leishmania transmission [21]. On the other hand, most sandfly species require a blood meal to develop their eggs, which facilitates contact with infected and uninfected hosts and transmission of Leishmania [22]. Targeting gravid female specimens for infection study is therefore a good method to ensure the presence of sandflies with at least one ingested blood meal and increase the chance to find Leishmania [23].
The objective of this work was to study the diversity, abundance, and L. infantum infection rate of phlebotomine sandflies in an area with low incidence of visceral leishmaniasis in northern Tunisia and to discuss the usefulness of quantitative kDNA real-time PCR and molecular investigation of gravid female specimens in screening Leishmania infection and determining L. infantum vectors.

2. Materials and Methods

2.1. Study Area

The study was conducted in the governorate of Siliana, North-West Tunisia, and specifically in the district of Siliana-North (Figure 1). Average of VL incidence rates observed in this district during 1996–2006 was one of the lowest of the VL endemic area [12]. Furthermore, no VL cases originating from this district have been reported by the Regional Health Directorate during the period 2008–2017, with VL being a notifiable disease in Tunisia.
Weather conditions of the region are known to be those of a Mediterranean semi-arid bioclimate, marked by a high degree of continentality and characterized by typical vegetation series of Pinus halepensis and Tetraclinus articulata. However, climatic parameters are subjected to annual variations, affecting annual degree of aridity and arid season duration.
Monthly rainfall and mean temperature recorded by the National Meteorological Office from January 2016 to December 2017 at the meteorological station of Ramlia located in Siliana-North district (Figure 1) were used to establish the ombro-thermic diagram (Figure 1) and allowed for the computation of the annual aridity index (Ia) by the following formula:
I a = P a T + 10  
where Pa (mm) corresponds to the annual precipitation and T (°C) corresponds to the annual average of temperature. The higher the value of this index, the lower is the aridity of the period. These parameters indicated that the climate in the study region was globally more arid in 2017 than in 2016. In fact, the annual aridity index dropped from 15.4 (in 2016) to 8.6 (in 2017), and the arid season was longer in 2017 than in 2016 (Figure 1).
The sampling station (36°10′ N/9°28′ E at 708 m above sea level) was located at Djema, a rural locality surrounded by a large cereal field and characterized by the presence of a variety of domestic animals (goats, sheep, chickens, dogs, donkeys) living in semi-open shelters or free in yards next to houses (Figure 2). The study was conducted on private land. The land’s owner gave permission to conduct the study on the site.

2.2. Sandfly Collection, Habitat, and Physiological State of Female Specimens

Sandflies were collected using CDC miniature light traps equipped with a fine net cage (John W. Hock Company, Florida, USA) during two successive entomological seasons (June–November 2016 and May–November 2017). CDC traps were set monthly during two consecutive nights (with replacement of traps after every night) in 3 habitats: one CDC indoor (ID), one outdoor (OD), and two in animal shelters (AS). They were set operating 1–2 h before sunset until 1–2 h after sunrise.
All collected sandfly specimens were stored in ethanol 70% then rapidly sorted by sex. All male specimens were conserved by habitat and date of capture for species identification. To identify physiological state of females, they were checked under stereomicroscope. Females were categorized according to their abdomen status and on the basis of blood digestion described by Dolmatova and Demina [22] into (i) unfed (non-distended abdomen without visible blood); (ii) fed (presence of blood meal in the abdomen); (iii) gravid (with abdomen full of eggs). For fed females, blood meal was classified on fresh (red fresh blood filling the abdomen) and digested (dark red to brown blood clot in the proximal part of the abdomen). All female specimens were stored individually for study of Leishmania infection.

2.3. Sandfly Species Identification and Study of Leishmania Infection

Sandfly male specimens were morphologically identified. Briefly, head and genitalia were removed and mounted in Marc-André solution (40 g of chloral hydrate, 30 mL glacial acetic acid, and 30 mL of distelled water) observed under a microscope and identified at the species level according to morphological characteristics of aedeagus, stylopodite, and coxopodite, as previously described [24].
DNA was extracted separately from each female individual. Briefly, the whole body of female specimens were washed twice with distilled water, and then DNA was extracted using the method described by Ready et al. (1991) [25] with minor modifications. DNA was eluted in 20 µL TE and stored at −20 °C until use.

2.3.1. Molecular Screening of Leishmania Infection by kDNA qPCR

Leishmania infection was first screened by pools of 5 sandflies’ DNA. Then, DNAs from sandflies grouped into each positive pool were analyzed individually to determine the positive specimens in each pool and to quantify the parasite load of each specimen. kDNA qPCR was performed as described by Mary et al. [26]. Briefly, qPCR was conducted in a final volume of 25 µL by using a TaqMan universal master mixture (Applied Biosystems, Foster City, CA, USA) containing 100 mM of each primer (5′-CTTTTCTGGTCCTCCGGGTAGG-3′) and (5′-CCACCCGGCCCTATTTTACACCAA-3′), 50 mM of probe (FAM 5′-TTTTCGCAGAACGCCCCTACCCGC-3′), and 1 µL of DNA extract. The DNA was amplified in an Applied Biosystems® (Applied Biosystems, Foster City, CA, USA) for 40 cycles at 95 and 60 °C. A dilution series of Leishmania DNA extracted from 5 × 106  L. infantum promastigotes was used (from 50 to 0.005 parasites/µL) to generate a standard curve for Leishmania quantification. Real-time PCR was considered positive for a qPCR threshold corresponding to 0.01 parasites/PCR reaction for the pools and 0.05 parasites/PCR reaction for individual specimens. Taking into account the (i) pooling dilution, (ii) the amount of biological sample (1 μL of sample DNA), and (iii) the elution volume of the extracted DNA (20 μL), the selected qPCR threshold corresponded to one parasite DNA per sandfly specimen.

2.3.2. Identification of Leishmania Species by ITS1 PCR Sequencing

All positive kDNA qPCR specimens were systematically amplified by ITS1 PCR as previously described [27] using LITSR (5′-CTGGATCATTTTCCGATG-3′) and L5.8S (5′-TGATACCACTTATCGCACTTA-3′) primers. Analysis on a 2% agarose gel was used to verify the amplified product size. ITS1 PCR products were purified by using ExoSAP (ThermoScientific, EU), and Leishmania species identification was performed by DNA sequencing using an ABI Prism® Big Dye™ Terminator, Cycle Sequencing Ready Reaction Kit, and AB1 3130 sequencing system (ABI, PE Applied Biosystems), with the same primers used for PCR. The sequence (LiDj033F4) was deposited in GenBank under the accession number OL804107.

2.3.3. Identification of Sandfly Species by DNA Barcoding

All positive kDNA qPCR female specimens as well as P. perniciosus and P. perfiliewi specimens collected and microscopically identified during our study were molecularly typed on the basis of DNA barcoding. The mitochondrial cytochrome c oxidase gene subunit 1 (COI) was amplified using primers LCO 1490 (5′ GGTCAACAAATCATAAAGATATTGG-3′) and HCO 2198 (5′TAAACTTCAGGGTGACCAAAAAATCA-3′) as previously described [28]. Analysis on a 1% agarose gel was used to verify the amplified product size (about 650 bp). PCR products were purified and cycle-sequenced as described previously. The generated sequences, Dj030F1, Dj032B1, Dj033F4, DJF01, DjM02, Dj033E10 and DjM01 were deposited in GenBank under the accession numbers OL814952, OL814953, OL814954, OL814955, OL819881, OL821654 and OL821655 respectively. MEGA software version 7.0.26 (Kumar, Stecher, and Tamura 2016) was used to conduct multiple sequence alignments and phylogenic analysis.

2.4. Data Analysis

2.4.1. Species Diversity

To compare either species diversity in a given habitat over time or species diversity in the different habitats, two diversity indexes widely used in ecology, and their corresponding equitability indexes were used [29,30]. Shannon’s index (H’) takes into account rare species and was used to study the temporal changes in a giving habitat. The Simpson’s index (D) gives more weight to dominant species than to rare ones, and it was used to make a comparison between habitats.

2.4.2. Statistical Analysis

Statistical analysis was performed using the MedCalc Statistical software (version 11.4.4.0). We compared median sandfly densities using the Mann–Whitney nonparametric test. The Spearman’s non-parametric correlation test was applied to establish the relationship between P. perniciosus and P. perfiliewi densities. The chi-squared test was used for comparison of relative abundances. Fisher’s exact test was applied to a 2 × 2 contingency table as an alternative to the chi-squared test to compare the proportion of qPCR-positive female sandflies in (i) gravid vs. fed specimens or (ii) gravid vs. unfed specimens. The significance level was set at 5% for all tests.

3. Results

3.1. Sandfly Fauna in the Study Site

During the two-year period, a total of 6125 specimens were collected: 2106 were collected during 2016 and 4019 during 2017; 5674 specimens were caught in AS, 243 in ID, and 208 in OD (Table 1). Cumulatively, 3336 females and 2789 males (sex-ratio = 0.84) were caught. Females outnumbered males during the two years of capture (Table 1).
Among the collected females, 2575 (77.2%) were unfed, 735 (22%) were fed, and only 26 (0.8%) were gravid. Physiological state of females and degree of blood digestion is reported in Table 2.

3.2. Species Diversity

On the basis of morphological identification of male specimens, six species were identified, namely, P. (Larroussius) perniciosus (n = 1355, 48.58%), P. (Larroussius) perfiliewi (n = 1281, 45.93%), Sergentomyia (Sergentomyia) antennata (n = 93, 3.33%), P. (Paraphlebotomus) sergenti (n = 48, 1.72%), P. (Phlebotomus) papatasi (n = 9, 0.32%), and P. (Larroussius) longicuspis (n = 3, 0.11%); P. perniciosus and P. perfiliewi corresponded to 94.5% of sandfly male specimens, whereas P. longicuspis, P. papatasi, P. sergenti, and S. antennata were rare species.
In AS, diversity was fairly steady during the two years (H’/H’max = 0.55 in 2016 and H’/H’max = 0.48 in 2017) with the dominance of the two species P. perfiliewi and P. perniciosus during the two-year period (D/Dmax = 0.52 in 2016 and D/Dmax = 0.59 in 2017) (Table 3). Indoor and OD ecological diversity was higher in 2016 than in 2017 (H’/H’max = 0.96 and H’/H’max = 1 in 2016 versus H’/H’max = 0.28 and H’/H’max = 0.21 in 2017). In 2017, P. perniciosus was the only dominant species in ID (pi = 0.9, D/Dmax = 0.23), whereas P. perfiliewi was the only dominant species in OD (pi = 0.95, D/Dmax = 0.15) (Table 3).

3.3. Sandfly Density

Sandfly density was computed, taking into account the number of CDCs set and number of nights of capture, and density was expressed as number of sandfly specimens per CDC.night. Globally, in the whole study site, sandfly density varied from 2 to 250 sandflies/CDC.night in 2016 and from 2 to 286 sandflies/CDC.night in 2017. Although the number of sandfly specimens collected in 2017 was higher than in 2016, no statistical differences were shown between median densities observed in 2016 and 2017.
The recorded densities in the different habitats during each year of collection are reported in Figure 3. In AS, the highest densities were recorded during July 2016 (487 sandflies/CDC.night) and August–September 2017 (535 and 356 sandflies/CDC.night, respectively). In indoor and OD, sandfly densities were usually low (≤56 sandflies/CDC.night). The highest density in OD was observed in August 2017, while the highest density in ID was shown in September 2017 (56 and 46 sandflies/CDC.night, respectively) (Figure 3).

3.4. Seasonal Dynamics

Seasonal density pattern of sandflies during the two-year period is reported in Figure 4. Globally, sandfly activity started in May–June and ended in October. However, in 2017, adults were caught during a longer period than in 2016. On the other hand, in 2017, a sharp bi-modal trend, with a first minor peak in June and a major peak in August–September, was clearly shown. Furthermore, during this second year of collection, the male/female ratio was higher during the early period of collection, whereas females outnumbered males at the end of the season (Figure 4).
Monthly P. perniciosus and P. perfiliewi male densities showed a significant positive correlation during the two-year period (Spearman’s coefficient ρ = 0.83, p = 0.001) and fairly the same seasonal density pattern with one major peak in July 2016 and two peaks during 2017 (a minor one in June and a major one in August–September) was shown (Figure 4). Nevertheless, some differences were observed between the two species. In fact, while the density of P. perfiliewi and its dynamics over time showed little variability during the two years of collection (density around 70 sandflies/CDC.night during a sharp peak in July 2016 or August 2017), the density of P. perniciosus showed more variability. In fact, the density of P. perniciosus, which was much lower than that of P. perfiewi during the 2016 peak (15 sandflies/CDC.night), rose by about five times (reaching around 70 sandflies/CDC.night), and its duration lengthened, resulting in an extended peak encompassing August–September 2017 (Figure 4). Furthermore, study of sandfly species abundance during peaks showed that relative abundance of P. perfiliewi during July 2016 was significantly higher than that of P. perniciosus (76.6% versus 13.6%) (p < 0.0001), whereas relative abundances of P. perniciosus during June, August, and particularly September 2017 were significantly higher than that of P. perfiliewi (58%, 58%, and 77.7% versus 24.7%, 41.8%, and 15.8%, respectively) (p < 0.0001).

3.5. Leishmania Infection

Among the 3336 collected females, 878 specimens (26.3%) were selected for DNA extraction (Table 4). Selection was conducted according to: (i) habitat (ID, OD, AS): almost all female specimens from ID and OD were used for DNA extraction, whereas about 20% of specimens from AS (where female number was much higher than ID and OD) were randomly selected; (ii) month and year of capture: almost all female specimens caught in June, August, September, and October 2016 and May, June, July, and October 2017 were used for DNA extraction, whereas about 20% of specimens collected either during the peak of July 2016 or the peak of August–September 2017 were randomly selected; and (iii) physiological state of females (unfed, fed, and gravid): all gravid females were selected for DNA extraction, whereas a comparable number of fed vs. unfed specimens were used to study infection. For fed females, degree of blood meal digestion was noted (Table 4).
Four specimens were positive by kDNA qPCR; no infection was found among the 311 specimens collected in 2016. The four kDNA qPCR-positive specimens were found among the 567 specimens (0.7%) collected in 2017. All four positive females were caught in September, one from ID and three from AS. Two females were fed, and two were gravid, and the two positive fed females harbored fresh blood (Table 4); infection rate was significantly higher in gravid (8.7%) than either in fed (0.8%) (p < 0.05) or in unfed (0%) (p < 0.01) specimens.
All infected sandfly specimens were identified as P. perniciosus on the basis of the 613 base pairs of mitochondrial cytochrome c oxidase (COI) gene subunit 1. In fact, the four sequences from infected sandflies collected during this study, Dj033E10, Dj033F4, Dj032B1 and Dj030F1 (gene bank accession numbers: OL821654, OL814954, OL814953, OL814952 respectively) showed homology higher than 98% with (i) the Mediterranean haplotypes of P. perniciosus available in GenBank and (ii) nucleotide sequences from P. perniciosus male and female specimens collected and microscopically identified during our study (DjM01 and DjF01, gene bank accession numbers: OL821655, and OL814955).

3.6. Leishmania Infantum Infection and Parasitic Load

Of the four kDNA qPCR-positive P. perniciosus females, only one was positive by ITS1 PCR; DNA sequencing identifying L. infantum (gene bank accession number: OL804107) (Table 5). This was a gravid specimen caught in AS during September 2017. It was highly infected with thousands of parasites (Table 5). The other three specimens were negative by ITS1 PCR, and Leishmania species identification was not possible using this target. All harbored less than 20 parasites (Table 5).

4. Discussion

Sandflies were collected monthly through two entomological seasons by CDC light traps. This sampling method commonly used in VL eco-epidemiological studies allows for the capturing of several species of the Larroussius subgenus that are the potential vectors of L. infantum, the causative agent of VL [18]. Moreover, light traps are useful to collect large numbers of samples, allowing for the measurement of the relative changes in abundance over time [8,31]. They are also useful for determining the temporal activity of species as they capture active sandflies [8,31]. This method, however, may be biased towards capture of unfed female specimens rather than fed and gravid ones. In fact, unfed females are more active in looking for a host to take blood meal, so they are more attracted by CDC traps [32].
Contrary to other studies, females outnumbered males during the two seasons of capture [8]. Furthermore, unfed specimens were the most abundant, whereas gravid females were rarely caught. This seems to be related to method of capture. However, it is important to note that female sandflies were not dissected and visually classified as unfed, fed, and gravid. Thus, a small volume of digested blood, as well as eggs in the first stages of development, may not be revealed and, therefore, numbers of fed and gravid females may be underestimated.
Species identification of male specimens showed that P. perniciosus and P. perfiliewi were the most abundant species in the study site localized in the semi-arid bioclimatic zone. It is well known that these two Larroussius species prevail within a specific bioclimatic continuum from the humid to arid zones of Western Mediterranean countries with inter-specific relationships, with P. perfiliewi and P. perniciosus preferences being for humid and semi-arid zones, respectively [18,33]. In Iran, P. perfiliewi is considered an intermediate species between xerophylious and hygrophylious species that prefer a humid region [34]. In Tunisia, global aridity of the center and the south of the country seems to limit the distribution of this species to the north above the “dorsal” upland [24]. However, P. perfiliewi was reported with high prevalence below the “dorsal” in irrigated zones [35]. This highlights the importance of micro-climate impact on sand fly fauna structure. On the other hand, the sampling station was located in anthropic rural biotope characterized by the presence of many domestic animals such as dogs, sheep, and chicken, which may explain the abundance of Larroussius species in the sample collection [18,31,36,37,38].
To study species diversity in the different habitats (AS, ID, and OD), the Shannon’s and Simpson’s diversity indexes and their correspondent equitability indices were computed. By using these indices, it was possible to follow the variation of the species diversity over time, but also to make comparison between habitats [39]. In AS, equitability indices showed that P. perniciosus and P. perfiliewi were dominant through the two-year period. Furthermore, the high sandfly density observed in AS likely determined the global trend in the study site [13,18]. Indeed, AS seems to constitute ecological niche of these sandfly species [40]. It is usually characterized by biotic and abiotic factors suitable for larval development; adult activity; and female’s rest such as organic matter produced by domestic animals, vertebrate host source of blood meals, suitable brightness, and humidity conditions. Being close to houses, this prolific habitat may maintain and control Phlebotomine abundance ID and OD, where vector contact with humans and dogs is the most frequent. In these two latter habitats, high species diversity was observed in 2016, and a switch toward a quasi-exclusive presence of P. perfiliewi OD and P. perniciosus ID was shown in 2017.
Globally, adult sandfly activity in the studied site started in May–June and extended to October with a large variation in sandfly densities during the season, which was in accordance with phenology studies performed in temperate Mediterranean region [8]. The bimodal trend observed in study site in 2017 likely indicated the occurrence of two sandfly generations at least and was already reported for both P. perniciosus and P. perfiliewi in Tunisia [13,41]. Interestingly, the multi-seasonal entomological survey evidenced two different monthly sandfly density patterns according to the year of collection. In 2016, a sharp peak of abundance was observed in July, while a long-lasting peak was observed during August–September 2017. Moreover, the magnitude of density peak observed in 2016 was mainly related to P. perfiliewi, whereas the long duration of the major peak observed in 2017 was associated with P. perniciosus. These specific species sandfly phenologies according to the year of capture are probably the result of local climatic conditions that could affect the two species differently [42]. In fact, despite that P. perfiliewi and P. perniciosus being usually sympatric species, P. perniciosus prefers more arid regions than P. perfiliewi [33]. Thus, fluctuations in local climatic factors resulting in an increase in the global aridity and a longer arid period such as that shown in the study area in 2017 may have allowed P. perniciosus to last longer at high density levels [37]. This correlation between seasonal variations in the densities with climatic factors, mainly temperature and rainfall, has been abundantly reported [8,34,40,43]. Phlebotomus perniciosus density was positively correlated with temperature and negatively correlated with rainfall in Spain [42].
In the present study, we applied a kDNA qPCR technique to study Leishmania infection in sandfly populations captured in a human dwelling located in Siliana North district, with a focus of human VL with low incidence. No infection was found during the first year of collection, and the qPCR positivity in sandflies was very low (0.7%) during the second year of capture, with Leishmania infections being found only in four P. perniciosus specimens caught at the end of the sandfly activity season. This corroborates reports of other authors in Algeria [44] and highlights the importance of multi-seasonal entomological survey and utility of the qPCR technique in the screening of Leishmania infection in areas of VL low incidence. The ITS1 PCR was, however, positive only in one infected specimen out of four. This could probably be explained by lower sensitivity of ITS1-based PCR than the kDNA-based molecular tool [45,46], with the three negative specimens harboring a very low parasitic load of less than 20 parasites. However, even though we performed a real-time PCR using a TaqMan probe, ensuring more specificity to Leishmania [26], we cannot exclude that kinetoplast DNA of other kinetoplastids has been amplified, with the latter being occasionally found in sandflies, and with sandflies even being vectors of some of them [47,48]. The kDNA PCR is considered to be the most sensitive method for diagnosing leishmaniasis since there are about 10,000 minicircles per parasite. However, these reactions generally either amplify genus- or subgenus-specific conserved regions and do not allow species identification. Recent development of a barcoding-like approach based on high-throughput sequencing of kDNA amplicons and minicircle sequence comparisons may allow in the future a reliable identification of Leishmania species in kDNA qPCR-positive specimens [49,50].
On the other hand, infection rate in gravid females was significantly higher than in fed specimens, which highlights the importance of screening gravid females for epidemiological surveys as suggested by other authors [23]. Targeting gravid female specimens at the end of the sandfly activity season seems particularly relevant in the case of Larroussius species, which are gonotrophically concordant and do not feed again until the blood meal is digested and the eggs are laid, with the probability of detecting infective promastigotes increasing progressively during the second and third ovarian cycles [51].
In this northern area of low incidence of VL, P. perniciosus, the main vector of VL in Western Mediterranean countries, was the only infected species. No infection was detected in P. perfiliewi, although it was as abundant as P. perniciosus in the study site during the two-year period. Interestingly, P. perniciosus infection was detected at a very low level only during the second year of collection when climatic conditions were more favorable to proliferation and long seasonal activity duration of this Larroussius species. The ITS1 PCR was, however, positive only in one infected specimen out of four. Therefore, the identification of L. infantum was only possible on one specimen of P. perniciosus. However, this latter specimen had a very high parasite load, underlining the major involvement of this Larroussius species in the transmission of L. infantum in VL foci of Northern Tunisia with low incidence of the disease [21].

Author Contributions

Conceptualization, M.W., J.G. and A.B.; methodology, M.W., J.G. and A.B.; formal analysis, M.W., J.G. and A.B.; investigation, M.W., M.B. and O.S.; data curation, M.W., J.G. and A.B.; writing—original draft preparation, M.W., J.G., K.A. and A.B.; writing—review and editing, M.W., J.G., K.A. and A.B.; visualization, M.W.; supervision, J.G. and A.B.; project administration, A.B.; funding acquisition, K.A. and A.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of Higher Education and Scientific Research, Tunisia, in the setting of LR 16-IPT-06 and LR 20-IPT-06 “Medical Parasitology, Biotechnology and Biomolecules”, Pasteur Institute of Tunis, and the APC was funded by Pasteur Institute of Tunis.

Institutional Review Board Statement

Not applicable.

Acknowledgments

We would like to thank Manel Zerzeri and Adel Rhim from the Pasteur Institute of Tunisia and Rafik Garni from the Pasteur Institute of Algeria for their kind contribution to the technical aspect of this work. We also thank the Regional Director of Public Health of Siliana.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pratlong, F.; Lami, P.; Ravel, C.; Balard, Y.; Dereure, J.; Serres, G.; El-Baidouri, F.; Dedet, J.P. Geographical distribution and epidemiological features of Old World Leishmania infantum and Leishmania donovani foci, based on the isoenzyme analysis of 2277 strains. Parasitology 2013, 140, 423–434. [Google Scholar] [CrossRef] [Green Version]
  2. World Health Organization (W.H.O). Control of the Leishmaniases: Report of a Meeting of the WHO Expert Commitee on the Control of Leishmaniases, Geneva, 22–26 March 2010; WHO Technical Report Series; World Health Organization: Geneva, Switzerland, 2010. [Google Scholar]
  3. Gramiccia, M.; Gradoni, L. The current status of zoonotic leishmaniases and approaches to disease control. Int. J. Parasitol. 2005, 35, 1169–1180. [Google Scholar] [CrossRef]
  4. Killick-Kendrick, R. Phlebotomine vectors of the leishmaniases: A review. Med. Vet. Entomol. 1990, 4, 1–24. [Google Scholar] [CrossRef]
  5. Rioux, J.A.; Killick-Kendrick, R.; Leaney, A.J.; Young, C.J.; Turner, D.P.; Lanotte, G.; Bailly, M. Ecology of leishmaniasis in the south of France. 11. Canine leishmaniasis: Successful experimental transmission from dog to dog by the bite of Phlebotomus ariasi Tonnoir, 1921. Ann. Parasitol. Hum. Comp. 1979, 54, 401–407. [Google Scholar] [CrossRef] [Green Version]
  6. Bettini, S.; Gramiccia, M.; Gradoni, L.; Atzeni, M.C. Leishmaniasis in Sardinia: II. Natural infection of Phlebotomus perniciosus Newstead, 1911, by Leishmania infantum Nicolle, 1908, in the province of Cagliari. Trans. R. Soc. Trop. Med. Hyg. 1986, 80, 458–459. [Google Scholar] [CrossRef]
  7. Maroli, M.; Gramiccia, M.; Gradoni, L. Natural infection of Phlebotomus perfiliewi with Leishmania infantum in a cutaneous leishmaniasis focus of the Abruzzi region, Italy. Trans. R. Soc. Trop. Med. Hyg. 1987, 8, 596–598. [Google Scholar] [CrossRef]
  8. Alten, B.; Maia, C.; Afonso, M.O.; Campino, L.; Jiménez, M.; González, E.; Molina, R.; Bañuls, A.L.; Prudhomme, J.; Vergnes, B.; et al. Seasonal dynamics of phlebotomine sandfly species proven vectors of Mediterranean leishmaniasis caused by Leishmania infantum. PLoS Negl. Trop. Dis. 2016, 10, e0004458. [Google Scholar] [CrossRef]
  9. Aoun, K.; Amri, F.; Chouihi, E.; Haouas, N.; Bedoui, K.; Benikhlef, R.; Ghrab, J.; Babba, H.; Chahed, M.K.; Harrat, Z.; et al. Épidémiologie de Leishmania (L.) infantum, L. major et L. killicki en Tunisie: Résultats et analyse de l’identification de 226 isolats humains et canins et revue de la littérature. Bull. Soc. Pathol. Exot. 2008, 101, 323–328. [Google Scholar] [CrossRef]
  10. Aoun, K.; Bouratbine, A. Cutaneous leishmaniasis in North Africa: A review. Parasite 2014, 21, 14. [Google Scholar] [CrossRef] [Green Version]
  11. Aoun, K.; Jeddi, F.; Amri, F.; Ghrab, J.; Bouratbine, A. Current epidemiological data on visceral leishmaniasis in Tunisia. Med. Mal. Infect. 2009, 39, 775–779. [Google Scholar] [CrossRef]
  12. Ben-Ahmed, K.; Aoun, K.; Jeddi, F.; Ghrab, J.; El-Aroui, M.A.; Bouratbine, A. Visceral leishmaniasis in Tunisia: Spatial distribution and association with climatic factors. Am. J. Trop. Med. Hyg. 2009, 81, 40–45. [Google Scholar] [CrossRef]
  13. Benabid, M.; Ghrab, J.; Rhim, A.; Ben-Romdhane, R.; Aoun, K.; Bouratbine, A. Temporal dynamics and Leishmania infantum infection prevalence of Phlebotomus perniciosus (Diptera, Phlebotominae) in highly endemic areas of visceral leishmaniasis in Tunisia. PLoS ONE 2017, 12, 0184700. [Google Scholar] [CrossRef] [Green Version]
  14. Ben-Salah, A.; Ben-Ismail, R.; Amri, F.; Chlif, S.; Ben-Rzig, F.; Kharrat, H.; Hadhri, H.; Hassouna, M.; Dellagi, K. Investigation of the spread of human visceral leishmaniasis in central Tunisia. Trans. R. Soc. Trop. Med. Hyg. 2000, 94, 382–386. [Google Scholar] [CrossRef]
  15. Barhoumi, W.; Fares, W.; Cherni, S.; Derbali, M.; Dachraoui, K.; Chelbi, I.; Ramalho-Ortigao, M.; Beier, J.C.; Zhioua, E. Changes of sand fly populations and Leishmania infantum infection rates in an irrigated village located in arid Central Tunisia. Int. J. Environ. Res. Public Health 2016, 13, 329. [Google Scholar] [CrossRef] [Green Version]
  16. Remadi, L.; Chargui, N.; Jimenez, M.; Molina, R.; Haouas, N.; González, E.; Chaabane-Banaouas, R.; Ben Salah, E.; Haddaji, M.; Chaabouni, Y.; et al. Molecular detection and identification of Leishmania DNA and blood meal analysis in Phlebotomus (Larroussius) species. PLoS Negl. Trop. Dis. 2020, 14, e0008077. [Google Scholar] [CrossRef]
  17. Chaouch, M.; Chaabane, A.; Ayari, C.; Ben-Othman, S.; Sereno, D.; Chemkhi, J.; BenAbderrazak, S. Investigation of natural infection of Phlebotomine (Diptera: Psychodidae) by Leishmania in Tunisian endemic regions. Parasite Epidemiol. Control 2021, 14, e00212. [Google Scholar] [CrossRef]
  18. Ghrab, J.; Rhim, A.; Bach-Hamba, D.; Chahed, M.K.; Aoun, K.; Nouira, S.; Bouratbine, A. Phlebotominae (Diptera: Psychodidae) of human leishmaniosis sites in Tunisia. Parasite 2006, 13, 23–33. [Google Scholar] [CrossRef] [Green Version]
  19. Galluzzi, L.; Ceccarelli, M.; Diotallevi, A.; Menotta, M.; Magnani, M. Real-time PCR applications for diagnosis of leishmaniasis. Parasites Vectors 2018, 11, 273. [Google Scholar] [CrossRef]
  20. Myskova, J.; Votypka, J.; Volf, P. Leishmania in sandflies: Comparison of quantitative polymerase chain reaction with other techniques to determine the intensity of infection. J. Med. Entomol. 2008, 45, 133–138. [Google Scholar] [CrossRef]
  21. Gonzalez, E.; Alvarez, A.; Ruiz, S.; Molina, R.; Jimenez, M. Detection of high Leishmania infantum loads in Phlebotomus perniciosus captured in the leishmaniasis focus of southwestern Madrid region (Spain) by real time PCR. Acta Trop. 2017, 171, 68–73. [Google Scholar] [CrossRef]
  22. Dolmatova, A.V.; Denima, N.A. Les phlébotomes (Phlebotominae) et Maladies qu’ils Transmettent; Documentation Technique; O.R.S.T.O.M.: Paris, France, 1971; p. 169. [Google Scholar]
  23. Torina, A.; Sole, M.; Reale, S.; Vitale, F.; Caracappa, S. Use of phlebotomine sandflies as indicator of Leishmania prevalence in an endemic area. Ann. N. Y. Acad. Sci. 2008, 1149, 355–357. [Google Scholar] [CrossRef]
  24. Croset, H.; Rioux, J.A.; Maistre, M.; Bayar, N. Les Phlébotomes de Tunisie (Diptera, Phlebotomidae). Mise au point systématique, chorologique et éthologique. Ann. Parasitol. Hum. Comp. 1978, 53, 711–749. [Google Scholar] [CrossRef] [Green Version]
  25. Ready, P.D.; Lainson, R.; Shaw, J.J.; Souza, A.A. DNA probes for distinguishing Psychodopygus wellcomei from Psychodopygus complexus (Diptera: Psychodidae). Mem. Inst. Oswaldo Cruz 1991, 86, 41–49. [Google Scholar] [CrossRef] [Green Version]
  26. Mary, C.; Faraut, F.; Lascombe, L.; Dumon, H. Quantification of Leishmania infantum DNA by a real-time PCR assay with high sensitivity. J. Clin. Microbiol. 2004, 42, 5249–5255. [Google Scholar] [CrossRef] [Green Version]
  27. Schönian, G.; Nasereddin, A.; Dinse, N.; Schweynoch, C.; Schallig, H.D.; Presber, W.; Jaffe, C.L. PCR diagnosis and characterization of Leishmania in local and imported clinical samples. Diagn. Microbiol. Infect. Dis. 2003, 47, 349–358. [Google Scholar] [CrossRef]
  28. Kumar, N.P.; Srinivasan, R.; Jambulingam, P. DNA barcoding for identification of sandflies (Diptera: Psychodidae) in India. Mol. Ecol. Res. 2012, 12, 414–420. [Google Scholar] [CrossRef]
  29. De Jong, T.M. A Comparison of three diversity indices based on their components of richness and evenness. Oikos 1975, 26, 222–227. [Google Scholar] [CrossRef]
  30. Mouillot, D.; Leprêtre, A.A. Comparison of species diversity estimators. Rese. Popul. Ecol. 1999, 41, 203–215. [Google Scholar]
  31. Prudhomme, J.; Rahola, N.; Toty, C.; Cassan, C.; Roiz, D.; Vergnes, B.; Thierry, M.; Rioux, J.A.; Alten, B.; Sereno, D.; et al. Ecology and spatiotemporal dynamics of sandflies in the Mediterranean Languedoc region (Roquedur area, Gard, France). Parasites Vectors 2015, 8, 642. [Google Scholar] [CrossRef] [Green Version]
  32. Alexander, B. Sampling methods for phlebotomine sandflies. Med. Vet. Entomol. 2000, 14, 109–122. [Google Scholar] [CrossRef]
  33. Rispail, P.; Dereure, J.; Jarry, D. Risk zones of human Leishmaniases in the Western Mediterranean basin: Correlations between vector sandflies, bioclimatology and phytosociology. Mem. Inst. Oswaldo Cruz 2002, 97, 477–483. [Google Scholar] [CrossRef] [Green Version]
  34. Norouzi, B.; Hanafi-Bojd, A.A.; Moin-Vaziri, V.; Noorallahi, A.; Azari-Hamidian, S. Ecology of sand flies (Diptera: Psychodidae, Phlebotominae) in a new focus of leishmaniasis in northern Iran. Acta Trop. 2020, 212, 105649. [Google Scholar] [CrossRef]
  35. Barhoumi, W.; Qualls, W.A.; Archer, R.S.; Fuller, D.O.; Chelbi, I.; Cherni, S.; Derbali, M.; Arheart, K.L.; Zhioua, E.; Beier, J.C. Irrigation in the arid regions of Tunisia impacts the abundance and apparent density of sand fly vectors of Leishmania infantum. Acta Trop. 2015, 141, 73–78. [Google Scholar] [CrossRef] [Green Version]
  36. Alcover, M.M.; Ballart, C.; Martín-Sánchez, J.; Serra, T.; Castillejo, S.; Portús, M.; Gállego, M. Factors influencing the presence of sandflies in Majorca (Balearic Islands, Spain) with special reference to Phlebotomus pernicious, vector of Leishmania infantum. Parasites Vectors 2014, 7, 421. [Google Scholar] [CrossRef] [Green Version]
  37. Rossi, E.; Bongiorno, G.; Ciolli, E.; Di Muccio, T.; Scalone, A.; Gramiccia, M.; Gradoni, L.; Maroli, M. Seasonal phenology, host-blood feeding preferences and natural Leishmania infection of Phlebotomus perniciosus (Diptera, Psychodidae) in a high-endemic focus of canine leishmaniasis in Rome province, Italy. Acta Trop. 2008, 105, 158–165. [Google Scholar] [CrossRef]
  38. Cazan, C.D.; Păstrav, I.R.; Györke, A.; Oguz, G.; Alten, B.; Mihalca, A.D. Seasonal dynamics of a population of Phlebotomus (Larroussius) perfiliewi Parrot, 1930 (Diptera: Psychodidae) in North-Eastern Romania. Parasitol. Res. 2019, 118, 1371–1384. [Google Scholar] [CrossRef]
  39. Tabbabi, A.; Ghrab, J.; Aoun, K.; Ready, P.D.; Bouratbine, A. Habitats of the sandfly vectors of Leishmania tropica and L. major in a mixed focus of cutaneous leishmaniasis in southeast Tunisia. Acta Trop. 2011, 119, 131–137. [Google Scholar] [CrossRef]
  40. Risueño, J.; Muñoz, C.; Pérez-Cutillas, P.; Goyena, E.; Gonzálvez, M.; Ortuño, M.; Bernal, L.J.; Ortiz, J.; Alten, B.; Berriatua, E. Understanding Phlebotomus perniciosus abundance in south-east Spain: Assessing the role of environmental and anthropic factors. Parasites Vectors 2017, 10, 189. [Google Scholar] [CrossRef] [Green Version]
  41. Croset, H.; Rioux, J.A.; Juminer, B.; Tour, S. Annual fluctuations of populations of Phlebotomus perniciosus Newstead, 1911, Phlebotomus perfiliewi Parrot, 1930 and Sergentomyia minuta parroti (Adler and Theodor, 1927)(Diptera-Psychodidae) in northern Tunisia. Arch. Inst. Pasteur Tunis 1970, 47, 43–56. [Google Scholar]
  42. Ballart, C.; Guerrero, I.; Castells, X.; Barón, S.; Castillejo, S.; Alcover, M.M.; Portús, M.; Gállego, M. Importance of individual analysis of environmental and climatic factors affecting the density of Leishmania vectors living in the same geographical area: The example of Phlebotomus ariasi and P. perniciosus in northeast Spain. Geosp. Health 2014, 8, 389–403. [Google Scholar] [CrossRef] [Green Version]
  43. Boutsini, S.; Athanasiou, L.V.; Spanakos, G.; Ntousi, D.; Dotsika, E.; Bisia, M.; Papadopoulos, E. Phlebotomine sandflies and factors associated with their abundance in the leishmaniasis endemic area of Attiki, Greece. Parasitol. Res. 2018, 117, 107–113. [Google Scholar] [CrossRef]
  44. Bennai, K.; Tahir, D.; Lafri, I.; Bendjaballah-Laliam, A.; Bitam, I.; Parola, P. Molecular detection of Leishmania infantum DNA and host blood meal identification in Phlebotomus in a hypoendemic focus of human leishmaniasis in northern Algeria. PLoS Negl. Trop. Dis. 2018, 12, e0006513. [Google Scholar] [CrossRef]
  45. Marcelino, A.P.; de Souza Filho, J.A.; e Bastos, C.D.V.; Ribeiro, S.R.; Medeiros, F.A.C.; Reis, I.A.; da Rocha, A.C.V.M.; Barbosa, J.R.; Paz, G.F.; Gontijo, C.M.F. Comparative PCR-based diagnosis for the detection of Leishmania infantum in naturally infected dogs. Acta Trop. 2020, 207, 105495. [Google Scholar] [CrossRef]
  46. Bensoussan, E.; Nasereddin, A.; Jonas, F.; Schnur, L.F.; Jaffe, C.L. Comparison of PCR assays for diagnosis of cutaneous leishmaniasis. J. Clin. Microbiol. 2006, 44, 1435–1439. [Google Scholar] [CrossRef] [Green Version]
  47. González, E.; Molina, R.; Jiménez, M. Rabbit trypanosome detection in Phlebotomus perniciosus sand flies from the leishmaniasis outbreak in Madrid, Spain. Acta Trop. 2018, 187, 201–206. [Google Scholar] [CrossRef]
  48. Svobodova, M.; Radrova, J. Phlebotomine Sandflies-Potential Vectors of Avian Trypanosomes. Acta Protozool. 2018, 57, 53–59. [Google Scholar]
  49. Kocher, A.; Valiere, S.; Banuls, A.L.; Murienne, J. High-throughput sequencing of kDNA amplicons for the analysis of Leishmania minicircles and identification of Neotropical species. Parasitology 2018, 145, 585–594. [Google Scholar] [CrossRef]
  50. Bussotti, G.; Benkahla, A.; Jeddi, F.; Souiaï, O.; Aoun, K.; Späth, G.F.; Bouratbine, A. Nuclear and mitochondrial genome sequencing of North-African Leishmania infantum isolates from cured and relapsed visceral leishmaniasis patients reveals variations correlating with geography and phenotype. Microb. Genom. 2020, 6, mgen000444. [Google Scholar] [CrossRef]
  51. Killick-Kendrick, R.; Rioux, J.A. Mark-release-recapture of sandflies fed on leishmanial dogs: The natural life-cycle of Leishmania infantum in Phlebotomus ariasi. Parassitologia 2002, 44, 67–71. [Google Scholar]
Microorganisms 10 01012 g001 550
Figure 1. Geographical situation and ombro-thermic diagram (2016–2017) of the study area. (a) The map of Tunisia (left) shows the governorate of Siliana (shaded grey). The zoom (right) reports the location of Djema, the sampling site (star), and Ramlia, the meteorological station, both located in the district of Siliana-North (shaded grey). (b) Ombro-thermic diagram based on mean temperature and rainfall recorded in Ramlia. Arid period extended from June to August in 2016 and from May to September in 2017.
Figure 1. Geographical situation and ombro-thermic diagram (2016–2017) of the study area. (a) The map of Tunisia (left) shows the governorate of Siliana (shaded grey). The zoom (right) reports the location of Djema, the sampling site (star), and Ramlia, the meteorological station, both located in the district of Siliana-North (shaded grey). (b) Ombro-thermic diagram based on mean temperature and rainfall recorded in Ramlia. Arid period extended from June to August in 2016 and from May to September in 2017.
Microorganisms 10 01012 g001
Microorganisms 10 01012 g002 550
Figure 2. Photographs showing landscape and habitats in the sampling site of Djema. (a) Global view of the site: a rural zone surrounded of hills and cultivated fields. (b) Location of houses in relation to each other. (c) Animal shelter with domestic animals close to houses. (d) Dogs around houses and animal shelters.
Figure 2. Photographs showing landscape and habitats in the sampling site of Djema. (a) Global view of the site: a rural zone surrounded of hills and cultivated fields. (b) Location of houses in relation to each other. (c) Animal shelter with domestic animals close to houses. (d) Dogs around houses and animal shelters.
Microorganisms 10 01012 g002
Microorganisms 10 01012 g003 550
Figure 3. Sandfly densities according to habitat and year of capture. All monthly sandfly densities recorded by habitat and year of capture are plotted. For each group of sandfly collection, the median density is reported by a marker and a horizontal line.
Figure 3. Sandfly densities according to habitat and year of capture. All monthly sandfly densities recorded by habitat and year of capture are plotted. For each group of sandfly collection, the median density is reported by a marker and a horizontal line.
Microorganisms 10 01012 g003
Microorganisms 10 01012 g004 550
Figure 4. Seasonal density pattern of sandflies during the two-year period. (a) Density patterns of male and female sandflies. (b) Density pattern of P. perniciosus and P. perfiliewi male specimens.
Figure 4. Seasonal density pattern of sandflies during the two-year period. (a) Density patterns of male and female sandflies. (b) Density pattern of P. perniciosus and P. perfiliewi male specimens.
Microorganisms 10 01012 g004
Table
Table 1. Number of male and female specimens according to habitat and year of capture.
Table 1. Number of male and female specimens according to habitat and year of capture.
YearHabitatNumber
CDC.night		Number
Males		Number
Females		Number
Sandflies		Sex Ratio
2016ID122237590.59
OD121327400.48
AS24871113620070.77
Total48906120021060.75
2017ID14111781891.42
OD14601041640.58
AS281712195436660.88
Total561883213640190.88
Total 1042789333661250.84
Table
Table 2. Physiological state of female specimens according to habitat and year of capture.
Table 2. Physiological state of female specimens according to habitat and year of capture.
YearHabitatNumber
Females		Number
Unfed (%)		Number Fed (%)Number Gravid (%)
Fresh BMDigested BM
2016ID3730 (81%)6 (16.3%)1 (2.7%)
15
OD2723 (85.2%)4 (14.8%)0 (0%)
13
AS1136888 (78.2%)246 (21.6%)2 (0.2%)
123123
Total1200941 (78.4%)256 (21.3%)3 (0.3%)
125131
2017ID7849 (62.8%)27 (34.6%)2 (2.6%)
1116
OD10486 (82.7%)15 (14.4%)3 (2.9%)
78
AS19541499(76.7%)437 (22.4%)18 (0.9%)
267170
Total21361634 (76.5%)479 (22.4%)23 (1.1%)
285194
Table
Table 3. Diversity of Phlebotominae sandfly fauna according to the sampling period and habitat.
Table 3. Diversity of Phlebotominae sandfly fauna according to the sampling period and habitat.
YearHabitatAbundance (pi)Diversity Indices
P. perfP. pernP. longP. papP. sergS. antSH’H’/H’maxDD/Dmax
2016ID10 (0.45)5 (0.23)-7 (0.32)--31.540.960.671.00
OD6 (0.46)7 (0.54)----21.001.000.541.00
AS665 (0.76)127 (0.15)--40 (0.05)39 (0.04)41.110.550.390.52
2017ID1 (0.01)100 (0.90)-2 (0.02)5 (0.05)3 (0.03)50.650.280.190.23
OD57 (0.95)1 (0.02)---2 (0.03)30.330.210.100.15
AS542 (0.32)1115 (0.65)3 (0.00)-3(0.00)49 (0.03)51.110.480.480.59
P. perf: P. perfiliewi; P. pern: P. perniciosus; P. pap: P. papatasi; P. serg: P. sergenti; P. long: P. longicuspis; pi: relative abundance; S: specific richness; H’: Shannon’s index; H’/H’max: equitability index; D: Simpson’s index; D/Dmax: equitability index.
Table
Table 4. Number of infected specimens according to habitat, physiological state, and month and year of capture.
Table 4. Number of infected specimens according to habitat, physiological state, and month and year of capture.
Number of Specimens in 2016Number of Specimens in 2017
CollectedDNA+ kDNA qPCRCollectedDNA+ kDNA qPCR
(Infection Rate)
HabitatID3732078731 (1.36%)
OD272401041030
AS1136255019543913 (0.76%)
Month of captureMay---770
June2828042420
July11282390770
August55011512750
September121209122234 (0.44%)
October2727017130
November00-000
Physiological
state		Unfed941178016343040
Fed25613004792402 (0.83%)
Gravid33023232 (8.69%)
Blood mealFresh125480285872 (2.29%)
Digested1318201941530
Total female specimens1200311021365674 (0.70%)
Table
Table 5. Parasitic load, ITS1 PCR results, and Leishmania species identification of infected specimens.
Table 5. Parasitic load, ITS1 PCR results, and Leishmania species identification of infected specimens.
CodeHabitatPhysiological StateMonth of CaptureEstimated Number of Parasites per Sandfly Female SpecimenITS1 PCR
Result		Leishmania Species
Dj033E10ASGravidSept 20172 parasitesNegative
Dj033F4ASGravidSept 201738,000 parasitesPositiveL. infatum
Dj032B1ASFedSept 20174 parasitesNegative
Dj030F1IDFedSept 201718 parasitesNegative
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Weslati, M.; Ghrab, J.; Benabid, M.; Souissi, O.; Aoun, K.; Bouratbine, A. Diversity, Abundance and Leishmania infantum Infection Rate of Phlebotomine Sandflies in an Area with Low Incidence of Visceral Leishmaniasis in Northern Tunisia. Microorganisms 2022, 10, 1012. https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms10051012

AMA Style

Weslati M, Ghrab J, Benabid M, Souissi O, Aoun K, Bouratbine A. Diversity, Abundance and Leishmania infantum Infection Rate of Phlebotomine Sandflies in an Area with Low Incidence of Visceral Leishmaniasis in Northern Tunisia. Microorganisms. 2022; 10(5):1012. https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms10051012

Chicago/Turabian Style

Weslati, Marwa, Jamila Ghrab, Meriem Benabid, Olfa Souissi, Karim Aoun, and Aïda Bouratbine. 2022. "Diversity, Abundance and Leishmania infantum Infection Rate of Phlebotomine Sandflies in an Area with Low Incidence of Visceral Leishmaniasis in Northern Tunisia" Microorganisms 10, no. 5: 1012. https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms10051012

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop