First Report of Two Jaculus Rodents as Potential Reservoir Hosts of Leishmania Parasites in Tunisia

This study shows, for the first time, natural Leishmania infection among Jaculus spp. in an endemic region of Tataouine, South Tunisia. To better characterize the transmission cycles in this complex focus of mixed transmission, Leishmania detection and species identification were performed by direct examination, internal transcribed spacer-1 (ITS1)-PCR-restriction fragment length polymorphism (RFLP), and sequencing of Jaculus (J.) jaculus (Linnaeus, 1758) and J. hirtipes (Lichtenstein, 1823) rodent species, which are frequently encountered in this area. Leishmania parasites were observed in 19 (41.3%) smears, while DNA parasites were detected in 28 (60.9%) Jaculus spp. spleens; among them, 12 (54.5%) were from 22 J. jaculus individuals and 16 (66.7%) were from 24 J. hirtipes individuals. Leishmania parasites were confirmed as Leishmania (L.) killicki (syn. L. tropica) in two J. hirtipes individuals (4.3%) and L. major (n = 24; 52.2%) in 10 J. jaculus and 14 J. hirtipes individuals. This finding represents the first evidence of natural infection with Leishmania parasites in rodents belonging to the Jaculus genus, providing the rationale to consider them as potential reservoir hosts of Old World Leishmania parasites in Tunisia and North Africa.


Introduction
Leishmaniasis is caused by protozoan parasites of the genus Leishmania. Among the three main clinical manifestations of this disease (cutaneous, visceral, and mucocutaneous), the cutaneous form is the most abundant [1]. In Tunisia, human cutaneous leishmaniasis (CL) is widely distributed and is prevalent in central and southern Tunisia. The annual 2 of 14 incidence of CL ranges from 1 to 10,000 cases according to climatic factors and the importance of herd immunity [2]. Furthermore, three epidemiological scenarios of cutaneous leishmaniasis have been described: sporadic (caused by Leishmania (L.) infantum Nicolle, 1908), chronic (caused by L. tropica Wright, 1903), and zoonotic (caused by L. major Yakimoff & Shokkor, 1914). The latter is predominant in the central and southern parts of the country [3]. In these areas, zoonotic cutaneous leishmaniasis (ZCL) circulates in an endemoepidemic mode [4] and is recognized as a serious public health problem [2,5].
Leishmaniasis caused by L. major and L. infantum are commonly zoonotic, with some species of rodents belonging to the Muridae family and domestic dogs as reservoir hosts, respectively [6,7]. In contrast, L. killicki (syn. L. tropica), the causative agent of chronic CL in urban endemic microfoci, known to be anthroponotic, has been strongly suspected to be zoonotic, especially with evidence of Ctenodactylus (C.) gundi (Rothmann, 1776) infection [8,9]. Several studies have shown that Leishmania parasites can infect small mammals in Tunisia, including the least weasel [7] and hedgehog [10]. Furthermore, the geographical distribution of CL cases caused by L. killicki in the region of Tataouine suggested the absence of the classically conventional reservoir host, C. gundi, while other small mammals and rodent species frequently encountered in the environment (Regional Directorate of Primary Care, Tataouine) could play a role as reservoir hosts, which remain to be identified. Rodent species belonging to the Jaculus genus were suspicious due to their wide distribution in central and southern Tunisia [11]. Moreover, these rodents were frequently faced by the regional health workers in Tataouine (Regional Directorate of Primary Care, Tataouine). This clinical and epidemiological context justified this study, which aimed to investigate rodent species belonging to the Jaculus genus in an old emerging focus of CL, that is, the Tataouine Governorate in southeastern Tunisia, to explore them as novel reservoir hosts of Leishmania parasites in the region.
This will better characterize transmission cycles and provide stronger evidence for intervention policies.

Rodent Trapping and Identification
Rodents were collected between May 2016 and June 2018 using butterfly nets at two study sites: BniMhira and Guermessa, belonging to the Tataouine Governorate, South Tunisia [12].
Some of the captured rodents (n = 41) were identified by amplification of the mitochondrial complete cytochrome b (Cytb) gene, as previously described [12].
To identify all captured rodents consistently (n = 46), we also used polymerase chain reaction (PCR) amplification of the mitochondrial partial Cytb gene after DNA extraction from their spleens using the QIAamp ® DNA Mini Kit (QIAgen, Hilden, Germany) according to the manufacturer's instructions. Modified vertebrate-universal specific primers (cytB1-F and cytB2-R) were used to amplify a 370 bp segment [13]. PCR reactions were conducted in a 50 µL reaction volume, including 0.2 mM dNTPs, 2.5 U of Taq DNA polymerase (GE Healthcare, Chicago, IL, USA), and 60 pmol of each primer (Forward 5 -CCATCCAACATYTCADCATGATGAAA-3 and Reverse 5 -GCHCCTCAGAATGATATTTGKCCTCA-3 ). Reactions were performed in an Applied Biosystems GeneAmp PCR System 2700 (Foster City, CA, USA) using the following cycling conditions: an initial denaturation step at 94 • C for 5 min; 30 cycles of 94 • C for 1 min, 54 • C for 1 min, and 72 • C for 1 min; and a final elongation step at 72 • C for 10 min. The amplicons were electrophoresed on 1% agarose gels and visualized under UV light. PCR products were purified using a QIAquick PCR purification kit (QIAGEN, Hilden, Germany) according to the manufacturer's instructions. Amplicon sequencing was performed using a conventional Big Dye Terminator Cycle Sequencing Ready Reaction Kit (Perkin Elmer, Applied Biosystems, Foster City, CA, USA) and an ABI373 automated DNA sequencer (Thermo Fisher Scientific, Waltham, MA, USA). The obtained sequences were edited using ChromasPro software (Technelysium Pty Ltd., South Brisbane, Australia). For homology searches, sequences were submitted to the National Center for Biotech- Before autopsy, each rodent was examined for cutaneous lesions in different parts of the body and skin. Clinical manifestations, including hyperpigmentation of the higher edge of the ear, depilation, small nodules, infiltration, and/or partial or total organ destruction, were recorded as previously described [7].

Parasitological Diagnosis
Extracts from two macerated ears of each caught rodent were subjected to microscopic examination (Direct Exam) and parasite culture in Novy-MacNeal-Nicolle (NNN) medium [7].
ITS1-PCR products were directly purified and sequenced in both directions using a conventional Big Dye Terminator Cycle Sequencing Ready Reaction Kit (Perkin Elmer, Applied Biosystems, Foster City, CA, USA) and an ABI373 Automated DNA Sequencer (Thermo Fisher Scientific, Waltham, MA, USA).

Collection of Publicly Available Sequences and Phylogenetic Analyses
For our study, we searched for partial and complete ITS1 sequences belonging to Leishmania parasites in the NCBI database using the following keywords: Leishmania, internal transcribed spacer 1, ITS1, Tunisia, host, and rodents. Three reference sequences for Leishmania spp. (accession numbers: KF152937, MK474651.1, and MK474641.1 for L. major, L. killicki, and L. infantum, respectively) were used to identify the sequences.
Moreover, two datasets were collected: the first included all sequences originating from Tunisia (Supplementary Table S2) and the second regrouped all sequences detected from rodent hosts worldwide (Supplementary Table S3). The retrieved sequences were aligned using the MEGA 6 software (Pennsylvania State University, State College, PA, USA) [16] to study the genetic variability of the ITS1 fragment from the two datasets. The output was edited and regions of unambiguously aligned sequences were retained for the final analysis. Accession numbers of sequences used in the present phylogenetic study are provided in Supplementary Tables S2 and S3.
Maximum likelihood (ML) phylogenetic analysis "SMS: Smart Model Selection in PhyML" (Montpellier University, Montpellier, France) [17] was run to determine the most suitable model of DNA evolution for the two haplotype sets of sequences using hierarchical log-likelihood-ratio tests and the Akaike information criterion. The best evolutionary model for both datasets was the GTR +G model with estimated gamma distribution shape parameters of 0.915 and 0.497 based on the Tunisian sequences (Supplementary Table S2) and those detected among rodents (Supplementary Table S3), respectively.
Furthermore, identified haplotypes from the two sets of sequences, and their population genetic analysis, including (i) the number of haplotypes (H), (ii) the haplotype diversity (Hd), (iii) the nucleotide diversity (π), (iv) the number of segregating sites (S), (v) the total number of mutations (η), (vi) the average number of nucleotide differences (k), and (vii) the Tajima's D test for neutrality, were realized and calculated using the DnaSP (Julio Rozas & Universitat de Barcelona, Barcelona, Spain) [18].

Statistical Analysis
Categorical variables are summarized as frequency counts and percentages. Comparisons between group percentages were performed using the chi-squared test or Fisher's exact test. Statistical significance was set at p ≤ 0.05. All statistical analyses were performed using STATA software version 11 (StataCorp LLC, TX, USA).

Rodent Identification
A total of 46 rodents belonging to the Jaculus genus were captured. Phylogenetic analysis of the complete Cytb gene, as previously described [12], identified 21 as belonging to Jaculus (J.) jaculus and 20 to J. hirtipes. A total of 39 short Cytb sequences were successfully amplified, allowing the identification of all rodent populations (Supplementary Table S1). Blasting 360 bp sequences allowed us to identify all the captured rodents as 22 (47.8%) belonging to J. jaculus and 24 (52.2%) to J. hirtipes. No discordance was observed between the two studies. The distribution of these rodents according to capture biotope and gender is illustrated in Table 1.

Leishmania Infection
The process of Leishmania diagnosis and identification among Jaculus rodents is shown in the flowchart in Figure 1.

Clinical Manifestations
Clinical examination of captured rodents revealed the presence of skin lesions as small nodules in the back and tail, and, more specifically, in both ears ( Figure 2

Clinical Manifestations
Clinical examination of captured rodents revealed the presence of skin lesions as small nodules in the back and tail, and, more specifically, in both ears ( Figure 2

Clinical Manifestations
Clinical examination of captured rodents revealed the presence of skin lesions as small nodules in the back and tail, and, more specifically, in both ears ( Figure 2

Leishmania Infection Prevalence
The presence of a confirmed Leishmania infection varied among Jaculus specimens, according to the technique used ( Figure 1).  Figure 3). Specifically, seven (31.8%) J. jaculus and twelve (50%) J. hirtipes showed the presence of amastigote forms. Unfortunately, we did not succeed in growing Leishmania isolates from these rodents on Novy-MacNeal-Nicolle (NNN) medium because of the high rate of contamination.

Leishmania Infection Prevalence
The presence of a confirmed Leishmania infection varied among Jaculus specimens, according to the technique used ( Figure 1).

Parasitological Methods
For all studied rodents, the infection prevalence was 41.3% (19/46) by microscopy ( Figure 3). Specifically, seven (31.8%) J. jaculus and twelve (50%) J. hirtipes showed the presence of amastigote forms. Unfortunately, we did not succeed in growing Leishmania isolates from these rodents on Novy-MacNeal-Nicolle (NNN) medium because of the high rate of contamination.

Combined Methods
Combining these two methods (parasitological and molecular), Leishmania parasites were detected in 37 of the 46 studied Jaculus individuals, reaching a positivity rate of 80.4%. Moreover, there is a non-significant difference between positivity rate of both rodent species (p-value not significant); it was 77.3% (17/22) among J. jaculus and 83.3% (20/24) among J. hirtipes.

Combined Methods
Combining these two methods (parasitological and molecular), Leishmania parasites were detected in 37 of the 46 studied Jaculus individuals, reaching a positivity rate of 80.4%. Moreover, there is a non-significant difference between positivity rate of both rodent species (p-value not significant); it was 77.3% (17/22) among J. jaculus and 83.3% (20/24) among J. hirtipes.
Second, 82% (23/28) of the ITS1-PCR-positive samples were successfully sequenced, and their ML phylogenetic tree was constructed using three reference sequences for each Leishmania spp. (accession numbers of Leishmania ITS1 sequences used as references: KF152937, MK474651.1, and MK474641.1 for L. major, L. killicki, and L. infantum), confirming the results of the RFLP analysis (Figure 4).
Second, 82% (23/28) of the ITS1-PCR-positive samples were successfully sequenced, and their ML phylogenetic tree was constructed using three reference sequences for each Leishmania spp. (accession numbers of Leishmania ITS1 sequences used as references: KF152937, MK474651.1, and MK474641.1 for L. major, L. killicki, and L. infantum), confirming the results of the RFLP analysis (Figure 4).  Only four samples were not double-checked by either method: one of them was not identified by RFLP and found to be L. major by sequencing, and three identified by RFLP as L. major were not successfully sequenced. Finally, two samples were not species-identified either by RFLP or sequencing.
Ultimately, combining the results of both methods, we identified Leishmania spp. in 26 samples: 10 J. jaculus (45.5%) and 14 J. hirtipes (58.3%) were found to be infected by L. major and 2 J. hirtipes (8.3%) were found to be infected by L. killicki.

Leishmania Phylogenetic Analysis
Phylogenetic analysis of all Tunisian Leishmania ITS1 sequences, performed using the ML method, showed three strongly supported clades with high bootstrap values (100) separating Leishmania parasite spp. (Figure 5).
Only four samples were not double-checked by either method: one of them was not identified by RFLP and found to be L. major by sequencing, and three identified by RFLP as L. major were not successfully sequenced. Finally, two samples were not species-identified either by RFLP or sequencing.
Ultimately, combining the results of both methods, we identified Leishmania spp. in 26 samples: 10 J. jaculus (45.5%) and 14 J. hirtipes (58.3%) were found to be infected by L. major and 2 J. hirtipes (8.3%) were found to be infected by L. killicki.

Leishmania Phylogenetic Analysis
Phylogenetic analysis of all Tunisian Leishmania ITS1 sequences, performed using the ML method, showed three strongly supported clades with high bootstrap values (100) separating Leishmania parasite spp. (Figure 5). Interestingly, among each of the three clades, we found three different haplotypes: one dominant haplotype encompassing the majority of the sequences regardless of their location, origin, and hosts, and two other haplotypes encompassing only one sequence each. The sequences detected among Jaculus spp. were grouped among the dominant haplotype for either L. major (Hap 2) or L. tropica (Hap 4).
Indices of genetic diversity among the haplotypes determined for the three species of Leishmania parasites are summarized in Table 2. Tajima's D test, which was used to determine the extent of neutral selection among these Tunisian parasites, generated a positive value (D = 1.2243; p > 0.1).  Interestingly, among each of the three clades, we found three different haplotypes: one dominant haplotype encompassing the majority of the sequences regardless of their location, origin, and hosts, and two other haplotypes encompassing only one sequence each. The sequences detected among Jaculus spp. were grouped among the dominant haplotype for either L. major (Hap 2) or L. tropica (Hap 4).
Indices of genetic diversity among the haplotypes determined for the three species of Leishmania parasites are summarized in Table 2. Tajima's D test, which was used to determine the extent of neutral selection among these Tunisian parasites, generated a positive value (D = 1.2243; p > 0.1). The same topology was also observed in a phylogenetic study of the worldwide ITS1 Leishmania sequences detected in rodent reservoir hosts. These sequences were sufficient to discriminate Leishmania species (Figure 6).
The same topology was also observed in a phylogenetic study of the worldwide ITS1 Leishmania sequences detected in rodent reservoir hosts. These sequences were sufficient to discriminate Leishmania species (Figure 6). The number of haplotypes increased compared to the previous phylogenetic analysis. In fact, L. tropica and L. infantum clades were represented by two and three haplotypes, respectively, with a dominant haplotype, while the remaining haplotypes contained a single sequence, whereas the L. major clade comprised six haplotypes with two dominant haplotypes (Hap 12 and Hap 17). Interestingly, each of the latter haplotypes comprised L. major ITS1 sequences from a single geographic location, namely Tunisia and Iran, and each of the four remaining haplotypes contained a single sequence. Similar to the previous analysis, the sequences described in this study were grouped in the dominant haplotype for both parasite spp. (Hap 17 for L. major clade and Hap 10 for L. tropica clade). The indices of genetic diversity among the haplotypes ( Figure 6) are summarized in Table 2. Tajima's D test, used to determine the extent of neutral selection among parasites detected from different rodent reservoir hosts, generated a negative value (D = -0.6; p > 0.05).

Discussion
To our knowledge, this is the first study to shed light on the importance of Leishmania infection and its clinical manifestations in two rodent spp. belonging to the Jaculus genus. This novel finding could constitute a breakthrough in understanding the transmission cycles of Old World cutaneous leishmaniasis, with implications for control. Confirmed or suspected reservoir hosts have been identified following the discovery of Leishmania parasites in numerous studies [7,9,[21][22][23][24][25][26].
In fact, Jaculus rodents have been found to be infected with many ecto-and endoparasites [27][28][29]. Although Jaculus spp. were found to be susceptible to experimental infection with Leishmania parasites, the infection seems to be chronic and does not impair animal health [30,31], even in Tunisia with L. tropica [32]. However, despite several attempts since 1986, Leishmania parasites have never been isolated and/or detected in these rodents [33][34][35][36][37]. It is worth noting that all of these studies were based on a small sample size of rodents (<10) and were conducted only using parasitological diagnostic methods.
In our study, the Leishmania infection prevalence among Jaculus spp. was 41.3% and 60.9% using parasitological and molecular tools, respectively. These high infection rates The number of haplotypes increased compared to the previous phylogenetic analysis. In fact, L. tropica and L. infantum clades were represented by two and three haplotypes, respectively, with a dominant haplotype, while the remaining haplotypes contained a single sequence, whereas the L. major clade comprised six haplotypes with two dominant haplotypes (Hap 12 and Hap 17). Interestingly, each of the latter haplotypes comprised L. major ITS1 sequences from a single geographic location, namely Tunisia and Iran, and each of the four remaining haplotypes contained a single sequence. Similar to the previous analysis, the sequences described in this study were grouped in the dominant haplotype for both parasite spp. (Hap 17 for L. major clade and Hap 10 for L. tropica clade). The indices of genetic diversity among the haplotypes ( Figure 6) are summarized in Table 2. Tajima's D test, used to determine the extent of neutral selection among parasites detected from different rodent reservoir hosts, generated a negative value (D = −0.6; p > 0.05).

Discussion
To our knowledge, this is the first study to shed light on the importance of Leishmania infection and its clinical manifestations in two rodent spp. belonging to the Jaculus genus. This novel finding could constitute a breakthrough in understanding the transmission cycles of Old World cutaneous leishmaniasis, with implications for control. Confirmed or suspected reservoir hosts have been identified following the discovery of Leishmania parasites in numerous studies [7,9,[21][22][23][24][25][26].
In fact, Jaculus rodents have been found to be infected with many ecto-and endoparasites [27][28][29]. Although Jaculus spp. were found to be susceptible to experimental infection with Leishmania parasites, the infection seems to be chronic and does not impair animal health [30,31], even in Tunisia with L. tropica [32]. However, despite several attempts since 1986, Leishmania parasites have never been isolated and/or detected in these rodents [33][34][35][36][37]. It is worth noting that all of these studies were based on a small sample size of rodents (<10) and were conducted only using parasitological diagnostic methods.
In our study, the Leishmania infection prevalence among Jaculus spp. was 41.3% and 60.9% using parasitological and molecular tools, respectively. These high infection rates fluctuate depending on the rodent spp. Indeed, they increased among J. hirtipes using both methods (50% by direct examination and 66.7% by ITS1-PCR), while they reached only 31.8% by direct examination and 54.4% by ITS1-PCR among tested J. jaculus individuals. Despite these variations, the overall infection prevalence in these rodent spp. (more than 80%) remains the highest compared to previous proportions observed among the potential and/or confirmed rodent reservoir hosts in Tunisia using parasitological, serological, or molecular tools [7,9,[21][22][23]38].
Approximately 60% of the overall infections among these rodents were asymptomatic, similar to previous findings in rodent reservoir hosts [9,39]. This might be explained by the fact that Leishmania parasites can disseminate to the spleen during the early stages of infection in asymptomatic animals [40]. This finding, in addition to the clinical manifestation and location of lesions, is in agreement with previous studies performed on reservoir hosts in Tunisia [7,9,22]. This could indicate a long history of co-evolution between these rodent spp. and the detected Leishmania parasites.
RFLP and sequencing approaches, followed by the construction of a phylogenetic tree (Figure 4), allowed for the identification of Leishmania spp. infecting these rodents. Surprisingly, both rodent Jaculus genera were found to be infected with L. major parasites. However, only J. hirtipes was found to be infected with L. killicki. Moreover, two specimens of J. hirtipes infected with L. killicki were captured in Guermessa, which is known to be an old focus of CL owing to L. killicki in Tunisia (according to the Regional Health Directorate and National Control Program of Cutaneous Leishmaniasis). This concomitant infection with two different Leishmania spp. has already been described in other rodents, such as Nesokia indica (L. major and L. infantum) [41], Psammomys spp. (L. major, L. killicki, and L. infantum) [21]; C. gundi (L. major and L. killicki), ref. [9] and Mus musculus Linnaeus, 1758 (L. tropica and L. infantum) [25].
The Tunisian Leishmania phylogenetic tree ( Figure 5) allowing the separation of the parasite spp. failed to show clear grouping according to the sequences, geographical origin, or related hosts. Moreover, this result confirms the heterogeneity of circulating Leishmania parasites in the country, with a dominant haplotype for each parasite sp., even though two more haplotypes were detected containing a single sequence for each Leishmania spp. (Supplementary Table S2).
Furthermore, the genetic diversity among Tunisian sequences was high, as evidenced by the high haplotype diversity (Hd = 0.707) and nucleotide diversity (π = 0.15), which were comparable to the values determined from kinetoplast DNA sequences of other tested Leishmania complexes [51]. This heterogeneity has been reported using other genetic markers among Tunisian Leishmania spp. in previous studies [52][53][54][55].
In addition, the positive value of Tajima's D test suggests a recent population bottleneck or balancing selection [51,56,57]. The latter seems to be more realistic and can explain the Leishmania parasite evolution in Tunisia, despite the fact that this value was not significant, most likely because of the small number of available Leishmania ITS1 sequences. In addition, genetic variability among Leishmania ITS1 genes at the species level, although not clearly established due to the reduced number of sequences of the determined haplotypes (one sequence for some haplotypes), can be related to the zymodemes of Leishmania parasites circulating in the country. As established for Tunisia, L. major MON-25 and L. killicki MON-8 were the exclusive zymodemes for these dermotropic species, whereas the viscerotropic species were represented by three different zymodemes: MON-24 (the most frequent), MON-1, and MON-80 [2]. In the last decade, many other zymodemes have been detected in Tunisia (L. tropica MON-317) [55] and neighboring countries [2]. This description can explain our results and confirm the presence of other potential zymodemes. These hypotheses should be addressed in future studies.
On the other hand, Leishmania from the phylogenetic tree of rodent hosts ( Figure 6) showed a high level of diversity, with 11 haplotypes detected, and a clear clustering according to the geographical origin in L. major clade (Supplementary Table S3). Similarly, the genetic diversity among rodent host sequences was high, as revealed by the high haplotype diversity (Hd = 0.765) and nucleotide diversity (π = 0.03).
In addition, the negative value of the Tajima's D test suggests, by definition, a population expansion or purifying selection [51]. This could explain the Leishmania parasite evolution among their rodent reservoir hosts, suggesting a purifying geographical link. Enriching databases with more Leishmania ITS1 sequences harbored by rodent reservoir hosts from all over the world could confirm this finding. This is the first description of Leishmania infection in rodents belonging to the Dipodidae family. In fact, previous studies using ITS1 sequences have described Leishmania infection among rodents belonging to two families: seven species belonging to five genera in the Muridea family and only one species from the Ctenodactylidea family ( Figure 6). Confirming that these new rodents are reservoir hosts for cutaneous leishmaniasis might pave the way for more targeted control measures.

Conclusions
The present study suggests the role of the two Jaculus spp. as potential reservoir hosts in the transmission cycles of L. major and L. killicki for both the former and only for J. hirtipes in the latter. Despite the small size of this study, the high proportion of infected specimens, the high density of these rodents, and their wide geographic spread in Tunisia [58] and the Old World [59] demand deeper investigations on larger sample sizes with potential striking consequences for CL transmission and control.
Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/microorganisms10081502/s1, Table S1: Accession numbers of short and long sequences of Cytochrome b used to identify rodents as belonging to J. jaculus and J. hirtipes.; Table S2: Accessions number of Tunisian Leishmania sequences, their hosts, locations and haplotypes.; Table S3: Accessions number of Leishmania sequences detected from rodent all over the world, their hosts, locations and haplotypes. Funding: This work was funded in part by the Tunisian Ministry of Higher Education and Scientific Research for the Laboratory of Transmission, Control and Immunobiology of Infections (LR11IPT02), Institut Pasteur de Tunis, and in the frame of the Centre d'Investigation Clinique (PRF-4; CIC2016IPT02). LeiSHield-MATI project "A multi-disciplinary international effort to identify clinical, molecular and social factors impacting cutaneous leishmaniasis" H2020-MSCA-RISE-2017-778298-LeiSHield-MATI partially supported the study. The funding bodies played no role in the design of the study, sample collection, analysis, interpretation, or writing the manuscript.

Institutional Review Board Statement:
All animal experiments complied with institutional, national, and international guidelines. The study, rodent collection, and handling protocols were approved by the Biomedical Ethics Committee of the Institut Pasteur de Tunis (document reference: 2016/13E/I/CIC).

Informed Consent Statement: Not applicable.
Data Availability Statement: Cytb and ITS1 sequences were submitted to the GenBank database under accession numbers MT340265-MT340304 and MW058066-MW058088, respectively. The dresses are as follows: www.ncbi.nlm.nih.gov/genbank (accessed on 20 May 2022).