Sequence Analysis of the Ribosomal DNA ITS2 Region for Phlebotomus papatasi (Diptera: Psychodidae)

Molecular identification of potential leishmaniasis vector species within the Phlebotomus (Euphlebotomus) argentipes species complex in Sri Lanka

Background

Leishmaniasis is an emerging vector-borne disease in Sri Lanka. Phlebotomus (Euphlebotomus) argentipes sensu lato Annandale and Brunette 1908 is suspected to be a potential vector. Three sibling species have been reported in the species complex based on analysis of morphological data. A study was carried out in different parts of Sri Lanka including cutaneous leishmaniasis prevailing localities to characterise the sibling species of Phlebotomus (Euphlebotomus) argentipes sensu lato and to establish their possible role in Leishmania transmission.

Methods

Sandflies were collected using cattle baited trap nets and mouth aspirator. They were identified based on existing taxonomic keys. Sequences of amplified cytochrome oxidase subunit I (CO I), cytochrome oxidase b (cyt b), internal transcribed spacer 2 (ITS2), 18s and 28s rDNA regions were analysed to confirm the number of sibling species. Vectorial capacity of the sibling species was checked by detecting human and Leishmania DNA.

Results

Sandflies collected using different techniques were processed for identification, parasite detection and molecular characterization. The 18s, 28s rDNA and cytochrome oxidase subunit I (CO I), internal transcribed spacer 2 (ITS2) and cytochrome b oxidase (cytb) sequences confirmed that the species belonged to the Argentipes complex. 18s and 28s sequences did not show any variation among the proposed sibling species. The phylogeny created from mitochondrial CO I and cytochrome b data and from the nuclear ITS2 region supports the existence of only two groups of flies (termed A and B) from Phlebotomus (Euphlebotomus) argentipes complex instead of the previously proposed three. The Leishmania mini-circle kinetoplastid, heat shock protein 70 (hsp70) and internal transcribed spacer I DNA along with human blood were detected from sibling species A only, which has not previously been considered to be a vector.

Conclusions

The taxonomy of the Sri Lankan Argentipes species complex is reassessed based on the molecular data. The existence of two sibling species is proposed; sibling species A has a long sensilla chaetica (> 50% length of the second antennal flagellomere) and sibling species B has a short sensilla cheatica (< 50%). Sibling species A is incriminated as a vector for leishmaniasis in Sri Lanka.

Morphometric and molecular differentiation of Phlebotomus (Phlebotomus) sandflies

The closely related sandfly species of the subgenus Phlebotomus namely, Phlebotomus papatasi (Scopoli, 1786), Phlebotomus duboscqi Neveu-Lemair, 1906 and Phlebotomus bergeroti Parrot, 1934 (Diptera: Psychodidae), are major vectors of Leishmania major (Kinetoplastida: Trypanosomatidae), the causative agent of cutaneous leishmaniasis in the Old World. Although allopatric in most of their distribution, the three species exist sympatrically in many places in central and eastern Sudan. Males of the three species can be distinguished using morphological characters; however, females are much harder to identify, thus complicating epidemiological studies. We carried out a morphometric and a molecular study to determine reliable morphological features and develop a polymerase chain reaction (PCR) assay for distinguishing females of these species. Males and females from each species were collected from sites in Sudan, East Africa and from one site in Mali, West Africa. Males were analysed morphologically and 20 characters and 10 character ratios were used in a stepwise discriminant analysis. This led to the identification of four characters with high discriminant loading scores sufficient for accurate male species identification. Male DNA was then used for the development of a PCR-based species diagnostic based on the second internal transcribed spacer (ITS2) of the ribosomal DNA. A set of four primers was developed to generate fragment sizes that are specific to each species and can reliably identify females as well as hybrid DNA. Both the morphometric and the molecular findings of this study have important applications for studies of the epidemiology of cutaneous leishmaniasis.

Molecular markers for Phlebotomus papatasi (Diptera: Psychodidae) and their usefulness for population genetic analysis

Three molecular typing tools: multilocus microsatellite typing, cytochrome b sequence analysis and internal transcribed spacer 2 (ITS2) sequence analysis, were evaluated for their usefulness in inferring the population structure of Phlebotomus papatasi sand flies. ITS2 sequence analysis did not prove suitable for inferring phylogenetic and population genetic relationships across P. papatasi sand flies. Microsatellite markers showed high resolution in differentiating globally distributed P. papatasi populations, whereas cytochrome b sequence analysis provided insight into the relationships between closely related populations from the Mediterranean. Population structure, differentiation and demographic history among P. papatasi are important for understanding patterns of dispersal in this species and for planning appropriate control measures.

Population structure and geographical subdivision of the Leishmania major vector Phlebotomus papatasi as revealed by microsatellite variation

Multi-locus microsatellite typing (MLMT) has been employed to infer the population structure of Phlebotomus papatasi (Scopoli) (Diptera: Psychodidae) sandflies and assign individuals to populations. Phlebotomus papatasi sandflies were collected from 35 sites in 15 countries. A total of 188 P. papatasi individuals were typed using five microsatellite loci, resulting in 113 different genotypes. Unique microsatellite signatures were observed for some of the populations analysed. Comparable results were obtained when the data were analysed with Bayesian model and distance-based methods. Bayesian statistic-based analyses split the dataset into two distinct genetic clusters, A and B, with further substructuring within each. Population A consisted of five subpopulations representing large numbers of alleles that were correlated with the geographical origins of the sandflies. Cluster B comprised individuals collected in the Middle East and the northern Mediterranean area. The subpopulations B1 and B2 did not, however, show any further correlation to geographical origin. The genetic differentiation between subpopulations was supported by F statistics showing statistically significant (Bonferroni-corrected P < 0.005) values of 0.221 between B2 and B1 and 0.816 between A5 and A4. Identification of the genetic structure of P. papatasi populations is important for understanding the patterns of dispersal of this species and to developing strategies for sandfly control.