Genetic structure of Mediterranean populations of the sandfly Phlebotomus papatasi by mitochondrial cytochrome b haplotype analysis

Genetic diversity and phylogeography of Phlebotomus argentipes (Diptera: Psychodidae, Phlebotominae), using COI and ND4 mitochondrial gene sequences

BACKGROUND

Phlebotomus argentipes complex is the primary vector for cutaneous leishmaniasis, a burgeoning health concern in contemporary Sri Lanka, where effective vector control is important for proper disease management. Understanding the genetic diversity of the P. argentipes population in Sri Lanka is vital before implementing a successful vector control program. Various studies have indicated that genetic divergence, caused by genetic drift or selection, can significantly influence the vector capacity of arthropod species. To devise innovative control strategies for P. argentipes, exploring genetic diversity and phylogeography can offer valuable insights into vector competence, key genetic trait transfer, and impact on disease epidemiology. The primary objective is to analyze the genetic diversity and phylogeography of the P. argentipes complex in Sri Lanka, based on two mitochondrial genomic regions in modern representatives of P. argentipes populations.

METHODOLOGY

A total of 159 P. argentipes specimens were collected from five endemic areas of cutaneous leishmaniasis and identified morphologically. Two mitochondrial regions (Cytochrome c oxidase subunit I (COI) and NADH dehydrogenase subunit 4 (ND4) were amplified using the total DNA and subsequently sequenced. Partial sequences of those mitochondrial genes were utilized to analyze genetic diversity indices and to explore phylogenetic and phylogeographic relationships.

PRINCIPAL FINDINGS

Among five sampling locations, the highest genetic diversity for COI and ND4 was observed in Hambantota (Hd-0.749, π-0.00417) and Medirigiriya (Hd-0.977, π-0.01055), respectively. Phylogeographic analyses conducted using COI sequences and GenBank retrieved sequences demonstrated a significant divergence of P. argentipes haplotypes found in Sri Lanka. Results revealed that they have evolved from the Indian ancestral haplotype due to historical- geographical connections of the Indian subcontinent with Sri Lanka.

CONCLUSIONS

Utilizing high-mutation-rate mitochondrial genes, such as ND4, can enhance the accuracy of genetic variability analysis in P. argentipes populations in Sri Lanka. The phylogeographical analysis of COI gene markers in this study provides insights into the historical geographical relationship between India and P. argentipes in Sri Lanka. Both COI and ND4 genes exhibited consistent genetic homogeneity in P. argentipes in Sri Lanka, suggesting minimal impact on gene flow. This homogeneity also implies the potential for horizontal gene transfer across populations, facilitating the transmission of genes associated with traits like insecticide resistance. This dynamic undermines disease control efforts reliant on vector control strategies.

Genomic analysis of two phlebotomine sand fly vectors of Leishmania from the New and Old World

Phlebotomine sand flies are of global significance as important vectors of human disease, transmitting bacterial, viral, and protozoan pathogens, including the kinetoplastid parasites of the genus Leishmania, the causative agents of devastating diseases collectively termed leishmaniasis. More than 40 pathogenic Leishmania species are transmitted to humans by approximately 35 sand fly species in 98 countries with hundreds of millions of people at risk around the world. No approved efficacious vaccine exists for leishmaniasis and available therapeutic drugs are either toxic and/or expensive, or the parasites are becoming resistant to the more recently developed drugs. Therefore, sand fly and/or reservoir control are currently the most effective strategies to break transmission. To better understand the biology of sand flies, including the mechanisms involved in their vectorial capacity, insecticide resistance, and population structures we sequenced the genomes of two geographically widespread and important sand fly vector species: Phlebotomus papatasi, a vector of Leishmania parasites that cause cutaneous leishmaniasis, (distributed in Europe, the Middle East and North Africa) and Lutzomyia longipalpis, a vector of Leishmania parasites that cause visceral leishmaniasis (distributed across Central and South America). We categorized and curated genes involved in processes important to their roles as disease vectors, including chemosensation, blood feeding, circadian rhythm, immunity, and detoxification, as well as mobile genetic elements. We also defined gene orthology and observed micro-synteny among the genomes. Finally, we present the genetic diversity and population structure of these species in their respective geographical areas. These genomes will be a foundation on which to base future efforts to prevent vector-borne transmission of Leishmania parasites.

Population Genetic Structure of Phlebotomus sergenti (Diptera: Psychodidae) Collected in Four Regions of Morocco Based on the Analysis of Cyt b and EF-1α Genes

Phlebotomus (Ph.) sergenti is the main vector of Leishmania (L.) tropica (Trypanosomatida: Trypanosomatidae), the causative agent of anthroponotic cutaneous leishmaniasis in Morocco. This species has an extended geographical distribution, wider than that of the parasite. The main objective of our study was to analyze the genetic diversity of Ph. sergenti collected in four foci in Morocco: Taza, Foum Jemâa, El Hanchane, and Ouarzazate. We studied a set of diversity and population structure indices by sequencing two markers; nuclear EF-1α and mitochondrial Cyt b from 175 individual sand flies. Our results showed a considerable degree of intraspecific polymorphism with a high number of haplotypes identified in both genes. Many polymorphic sites detected in the Cyt b sequences (SCyt b = 45) indicate that it is the most polymorphic marker showing a distinct distribution of haplotypes according to their geographical origin, whereas the EF-1α marker showed no geographical isolation. Analysis by Tajima's D and Fu's Fs tests revealed a possible recent expansion of the populations, especially with the EF-1α marker, showing significant values in Taza and Ouarzazate sequences. The present study revealed significant genetic diversity within Ph. sergenti populations in Morocco. The results warrant further research using a combination of more than two markers including mitochondrial and non-mitochondrial markers, which may provide more information to clarify the genetic status of Ph. sergenti.

Molecular taxonomy of phlebotomine sand flies (Diptera, Psychodidae) with emphasis on DNA barcoding: A review

The taxonomy and systematics of sand flies (Diptera, Psychodidae, Phlebotominae) are one of the pillars of research aimed to identifying vector populations and the agents transmitted by these insects. Traditionally, the use of morphological traits has been the main line of evidence for the definition of species, but the use of DNA sequences is useful as an integrative approach for their delimitation. Here, we discuss the current status of the molecular taxonomy of sand flies, including their most sequenced molecular markers and the main results. Only about 37% of all sand fly species have been processed for any molecular marker and are publicly available in the NCBI GenBank or BOLD Systems databases. The genera Phlebotomus, Nyssomyia, Psathyromyia and Psychodopygus are well-sampled, accounting for more than 56% of their sequenced species. However, less than 34% of the species of Sergentomyia, Lutzomyia, Trichopygomyia and Trichophoromyia have been sampled, representing a major gap in the knowledge of these groups. The most sequenced molecular markers are those within mtDNA, especially the DNA barcoding fragment of the cytochrome c oxidase subunit I (coi) gene, which has shown promising results in detecting cryptic diversity within species. Few sequences of conserved genes have been generated, which hampers higher-level phylogenetic inferences. We argue that sand fly species should be sequenced for at least the coi DNA barcoding marker, but multiple markers with different mutation rates should be assessed, whenever possible, to generate multilocus analysis.

Taxonomy, Population Structure and Genetic Diversity of Iranian Leishmania Strains of Cutaneous and Visceral Leishmaniasis

BACKGROUND

Despite the broad distribution of leishmaniasis in Iran, there is a little genetic information about the causative agents and epidemiological status of the disease. Genetic diversity of the parasite is suggested to be one of the factors, which influences the clinical manifestations of the disease. In this study, we investigated the genetic variations, population structure, and evolutionary history of Leishmania species from endemic foci of Iran.

METHODS

Fifty-two isolates from humans, canines, and rodents from different endemic foci of Iran were used to sequence the N-acetyl glucosamine-1-phosphate transferase (Nagt) gene. Phylogenetic and structure analyses were performed to investigate inter- and intra-species diversity of the Leishmania isolates.

RESULTS

In total, 10 haplotypes including L. major (n = 6), L. tropica (n = 2), L. infantum (n = 1) and L. turanica (n = 1) were identified across 52 isolates. Haplotype diversity (Hd) ranged from zero for L. infantum and L. turanica to 0.78 ± 0.136 for L. major. This study identified population structure of Leishmania isolates from different geographical regions of Iran. The results of the phylogenetic tree showed 4 distinct clades for each species of Leishmania. In addition, the highest intraspecies diversity was observed among L. major isolates. No correlation was observed between species and geographic distribution of haplotypes.

CONCLUSIONS

Leishmania isolates were identified at the species level using the Nagt gene, low variation within species indicates conservation of this gene in Leishmania. The results provide knowledge into the evolutionary history of Iranian Leishmania isolates.

Standing genetic variation in laboratory populations of insecticide-susceptible Phlebotomus papatasi and Lutzomyia longipalpis (Diptera: Psychodidae: Phlebotominae) for the evolution of resistance

Insecticides can exert strong selection on insect pest species, including those that vector diseases, and have led to rapid evolution of resistance. Despite such rapid evolution, relatively little is known about standing genetic variation for resistance in insecticide-susceptible populations of many species. To help fill this knowledge gap, we generated genotyping-by-sequencing data from insecticide-susceptible Phlebotomus papatasi and Lutzomyia longipalpis sand flies that survived or died from a sub-diagnostic exposure to either permethrin or malathion using a modified version of the Centers for Disease Control and Prevention bottle bioassay. Multi-locus genome-wide association mapping methods were used to quantify standing genetic variation for insecticide resistance in these populations and to identify specific alleles associated with insecticide survival. For each insecticide treatment, we estimated the proportion of the variation in survival explained by the genetic data (i.e., "chip" heritability) and the number and contribution of individual loci with measurable effects. For all treatments, survival to an insecticide exposure was heritable with a polygenic architecture. Both P. papatasi and L. longipalpis had alleles for survival that resided within many genes throughout their genomes. The implications for resistance conferred by many alleles, as well as inferences made about the utility of laboratory insecticide resistance association studies compared to field observations, are discussed.

First molecular identification of Leishmania major in Phlebotomus papatasi in an outbreak cutaneous leishmaniasis area in Iraq

Cutaneous leishmaniasis (CL) is transmitted by Phlebotomine sand fly vectors, among which Phlebotomus papatasi is prevalent in Western Asia, Northern Africa and Southern Europe, and it is known as a vector for Leishmania major parasite in the world. However, in Iraq, morphological studies showed that P. papatasi is a predominant sand fly species and hypothesised to transmit CL causing Leishmania species including L. major and L. tropica. Few studies have found Leishmania species in sand flies in mixed pools of samples in this country. Accurate identification of sand flies as vectors of Leishmania species is required in Iraq. The current study aims to identify sand fly species, using both morphological and molecular phylogenetic analyses, in a region where CL tends to be endemic. Furthermore, molecular phylogenetic analysis has also used to confirm Leishmania species in the sand fly samples collected in 11 villages between Diyala and Sulaymaniyah Provinces. For the first time, we have found L. major in three individual sand flies, one engorged (with fresh blood meal) and two non-engorged (without visible fresh blood meal) P. papatasi females in an area of CL outbreaks since 2014-till now due to civil wars and internal conflicts happen in the region. Further study should be performed on sand fly population and Leishmania reservoirs in this region.

Phlebotomus papatasi sand fly predicted salivary protein diversity and immune response potential based on in silico prediction in Egypt and Jordan populations

Phlebotomus papatasi sand flies inject their hosts with a myriad of pharmacologically active salivary proteins to assist with blood feeding and to modulate host defenses. In addition, salivary proteins can influence cutaneous leishmaniasis disease outcome, highlighting the potential of the salivary components to be used as a vaccine. Variability of vaccine targets in natural populations influences antigen choice for vaccine development. Therefore, the objective of this study was to investigate the variability in the predicted protein sequences of nine of the most abundantly expressed salivary proteins from field populations, testing the hypothesis that salivary proteins appropriate to target for vaccination strategies will be possible. PpSP12, PpSP14, PpSP28, PpSP29, PpSP30, PpSP32, PpSP36, PpSP42, and PpSP44 mature cDNAs from field collected P. papatasi from three distinct ecotopes in the Middle East and North Africa were amplified, sequenced, and in silico translated to assess the predicted amino acid variability. Two of the predicted sequences, PpSP12 and PpSP14, demonstrated low genetic variability across the three geographic isolated sand fly populations, with conserved multiple predicted MHCII epitope binding sites suggestive of their potential application in vaccination approaches. The other seven predicted salivary proteins revealed greater allelic variation across the same sand fly populations, possibly precluding their use as vaccine targets.

Population Genetics of Phlebotomus papatasi from Endemic and Nonendemic Areas for Zoonotic Cutaneous Leishmaniasis in Morocco, as Revealed by Cytochrome Oxidase Gene Subunit I Sequencing

Zoonotic cutaneous leishmaniasis (ZCL) caused by Leishmania major Yakimoff & Shokhor and transmitted by Phlebotomus papatasi (Scopoli) is a public health concern in Morocco. The disease is endemic mainly in pre-Saharan regions on the southern slope of the High Atlas Mountains. The northern slope of the High Atlas Mountains and the arid plains of central Morocco remain non-endemic and are currently considered high risk for ZCL. Here we investigate and compare the population genetic structure of P. papatasi populations sampled in various habitats in historical foci and non-endemic ZCL areas. A fragment of the mtDNA cytochrome oxidase I (COI) gene was amplified and sequenced in 59 individuals from 10 P. papatasi populations. Haplotype diversity was probed, a median-joining network was generated (F_ST) and molecular variance (AMOVA) were analyzed. Overall, we identified 28 haplotypes with 32 distinct segregating sites, of which seven are parsimony informative. The rate of private haplotypes was high; 20 haplotypes (71.4%) are private ones and exclusive to a single population. The phylogenetic tree and the network reconstructed highlight a genetic structuration of these populations in two well defined groups: Ouarzazate (or endemic areas) and Non-Ouarzazate (or nonendemic areas). These groups are separated by the High Atlas Mountains. Overall, our study highlights differences in terms of population genetics between ZCL endemic and non-endemic areas. To what extent such differences would impact the transmission of L. major by natural P. papatasi population remains to be investigated.

Development of a next generation DNA sequencing-based multi detection assay for detecting and identifying Leishmania parasites, blood sources, plant meals and intestinal microbiome in phlebotomine sand flies

Leishmaniasis is a disease caused by Leishmania parasites transmitted by phlebotomine sand flies (Diptera: Psychodidae). Human infections with different Leishmania species cause characteristic clinical manifestations; cutaneous or visceral leishmaniasis. Here we describe the development and application of a Miseq Next GenerationSequencing (NGS)-based Multi Detection Assay (MDA) designed to characterize metagenomics parameters pertinent to the sand fly vectors which may affect their vectorial capacity for Leishmania. For this purpose, we developed a MDA by which, DNA fragments were amplified through polymerase chain reactions (PCR) and then sequenced by MiSeq/NGS. PCR amplification was achieved using some published and some new primers designed specifically for identifying Leishmania spp. (ITS1), sand fly spp. (cytochrome oxidase I), vertebrate blood (Cytochrome b), plant DNA ribulose-1,5-bisphosphate carboxylase large subunit gene (rbcL), and prokaryotic micobiome (16 s rRNA). This MDA/NGS analysis was performed on two species of wild-caught sand flies that transmit different Leishmania spp. in two ecologically distinct, but geographically neighboring locations. The results were analyzed to identify, quantitate and correlate the measured parameters in order to assess their putative importance in the transmission dynamics of leishmaniasis.

Revision of the species composition and distribution of Turkish sand flies using DNA barcodes

Background

Currently, knowledge regarding the phlebotomine sand fly (Diptera: Psychodidae) fauna of Turkey is restricted to regions with endemic leishmaniasis. However, rapidly changing environmental and social conditions highlight concerns on the possible future expansion of sand fly-borne diseases in Turkey, promoting risk assessment through biosurveillance activities in non-endemic regions. Traditional morphological approaches are complicated by extensive cryptic speciation in sand flies, thus integrated studies utilizing DNA markers are becoming increasingly important for correct sand fly identification. This study contributes to the knowledge of the sand fly fauna in understudied regions of Turkey, and provides an extensive DNA barcode reference library of expertly identified Turkish sand fly species for the first time.

Methods

Fly sampling was conducted at 101 locations from 29 provinces, covering all three biogeographical regions of Turkey. Specimens were morphologically identified using available keys. Cytochrome c oxidase I (cox1) barcode sequences were analyzed both for morphologically distinct species and those specimens with cryptic identity. A taxon identity tree was obtained using Neighbor Joining (NJ) analysis. Species boundaries among closely related taxa evaluated using ABGD, Maximum Likelihood (ML) and haplotype network analyses. Sand fly richness of all three biogeographical regions were compared using nonparametric species richness estimators.

Results

A total of 729 barcode sequences (including representatives of all previously reported subgenera) were obtained from a total of 9642 sand fly specimens collected in Turkey. Specimens belonging to the same species or species complex clustered together in the NJ tree, regardless of their geographical origin. The species delimitation methods revealed the existence of 33 MOTUs, increasing the previously reported 28 recorded sand fly species by 17.8%. The richest sand fly diversity was determined in Anatolia, followed by the Mediterranean, and then the Black Sea regions of the country.

Conclusions

A comprehensive cox1 reference library is provided for the sand fly species of Turkey, including the proposed novel taxa discovered herein. Our results have epidemiological significance exposing extensive distributions of proven and suspected sand fly vectors in Turkey, including those areas currently regarded as non-endemic for sand fly-borne disease.

Electronic supplementary material

The online version of this article (10.1186/s13071-019-3669-3) contains supplementary material, which is available to authorized users.

The changing distribution of Leishmania infantum Nicolle, 1908 and its Mediterranean sandfly vectors in the last 140 kys

The understanding of the effects of past climatic changes on the distribution of vector arthropods can strongly support the understanding of the future potential impact of anthropogenic climatic change on the geographical risk of vector-borne diseases. The zoogeographical patterns of the European sandfly vectors may suffer the continuously changing climate of the last 140 kys. The former range of L. infantum and six Phlebotomus species were modelled for the Last Interglacial, the Last Glacial Maximum and the Mid-Holocene Periods. It was found that the potential distribution of the parasite was much smaller in the Last Glacial Period L. infantum mainly could persist in the western shelves of the Mediterranean Sea. West and East Mediterranean sandfly species inhabited partly distinct refugia. The Apennine Peninsula, Sicily and the Iberian refugium formed a habitat chain along with the coastal areas of the West Mediterranean Basin. There was no direct connection between the Eastern and the Western sandfly refugia in the last 140 kys. The modelled distribution of sandfly taxa for the Middle Holocene Period can explain the relict populations of sandfly taxa in such Central European countries. The former genetic studies strongly confirm the existence of the modelled glacial refugees.

Phylogenetic Studies

Phylogenetics is an important component of the systems biology approach. Knowledge about evolution of the genus Leishmania is essential to understand various aspects of basic biology of these parasites, such as parasite-host or parasite-vector relationships, biogeography, or epidemiology. Here, we present a comprehensive guideline for performing phylogenetic studies based on DNA sequence data, but with principles that can be adapted to protein sequences or other molecular markers. It is presented as a compilation of the most commonly used genetic targets for phylogenetic studies of Leishmania, including their respective primers for amplification and references, as well as details of PCR assays. Guidelines are, then, presented to choose the best targets in relation to the types of samples under study. Finally, and importantly, instructions are given to obtain optimal sequences, alignments, and datasets for the subsequent data analysis and phylogenetic inference. Different bioinformatics methods and software for phylogenetic inference are presented and explained. This chapter aims to provide a compilation of methods and generic guidelines to conduct phylogenetics of Leishmania for nonspecialists.

Genetic structure analysis of Amblyomma mixtum populations in Veracruz State, Mexico

Amblyomma mixtum Koch, 1844 parasitizes livestock, humans, and wildlife in Mexico. However, information on population genetics for this tick species in the country is missing. Tick samples were collected from livestock in ten regions across the state of Veracruz (22°28'N, 17°09'S, 93°36'E, 98°39'W) to analyze the genetic structure of   A. mixtum populations. Ticks were morphologically identified using taxonomic keys. In order to test the intra-specific variability of A. mixtum fragments of the mitochondrial gene 16S-rRNA and cytochrome oxidase subunit 1 (COI) were amplified. Ninety-six sequences were amplified from the 50 specimens' analyzed (96% amplification success). Eleven haplotypes were detected in 16S-rRNA gene and 10 more for COI. Neutrality tests showed negative results in most of the locations analyzed, which is indicative of an excess of recently derived haplotypes. However, these results were not statistically significant. Minimal union network analysis revealed that there is no separation of populations by geography, and that there is an overlap of several haplotypes among diverse populations. Significant genetic differentiation was not detected in the A. mixtum populations sampled in the state of Veracruz, Mexico, this may be due to the frequent movement of livestock hosts. This is the first report on the genetic structure of A. mixtum populations in Mexico.

Genetic diversity, phylogeography and molecular clock of the Lutzomyia longipalpis complex (Diptera: Psychodidae)

BACKGROUND

The Lutzomyia longipalpis complex has a wide but discontinuous distribution in Latin America, extending throughout the Neotropical realm between Mexico and northern Argentina and Uruguay. In the Americas, this sandfly is the main vector of Leishmania infantum, the parasite responsible for Visceral Leishmaniasis (VL). The Lu. longipalpis complex consists of at least four sibling species, however, there is no current consensus on the number of haplogroups, or on their divergence. Particularly in Argentina, there have been few genetic analyses of Lu. longipalpis, despite its southern expansion and recent colonization of urban environments. The aim of this study was to analyze the genetic diversity and structure of Lu. longipalpis from Argentina, and to integrate these data to re-evaluate the phylogeography of the Lu. longipalpis complex using mitochondrial markers at a Latin American scale.

METHODOLOGY/PRINCIPAL FINDINGS

Genetic diversity was estimated from six sites in Argentina, using a fragment of the ND4 and the 3´ extreme of the cyt b genes. Greatest genetic diversity was found in Tartagal, Santo Tomé and San Ignacio. There was high genetic differentiation of Lu. longipalpis in Argentina using both markers: ND4 (FST = 0.452, p < 0.0001), cyt b (FST = 0.201, p < 0.0001). Genetic and spatial Geneland analyses reveal the existence of two primary genetic clusters in Argentina, cluster 1: Tartagal, Santo Tomé, and San Ignacio; cluster 2: Puerto Iguazú, Clorinda, and Corrientes city. Phylogeographic analyses using ND4 and cyt b gene sequences available in GenBank from diverse geographic sites suggest greater divergence than previously reported. At least eight haplogroups (three of these identified in Argentina), each separated by multiple mutational steps using the ND4, are differentiated across the Neotropical realm. The divergence of the Lu. longipalpis complex from its most recent common ancestor (MRCA) was estimated to have occurred 0.70 MYA (95% HPD interval = 0.48-0.99 MYA).

CONCLUSIONS/SIGNIFICANCE

This study provides new evidence supporting two Lu. longipalpis genetic clusters and three of the total eight haplogroups circulating in Argentina. There was a high level of phylogeographic divergence among the eight haplogroups of the Lu. longipalpis complex across the Neotropical realm. These findings suggest the need to analyze vector competence, among other parameters intrinsic to a zoonosis, according to vector haplogroup, and to consider these in the design and surveillance of vector and transmission control strategies.

Population genetics analysis of Phlebotomus papatasi sand flies from Egypt and Jordan based on mitochondrial cytochrome b haplotypes

BACKGROUND

Phlebotomus papatasi sand flies are major vectors of Leishmania major and phlebovirus infection in North Africa and across the Middle East to the Indian subcontinent. Population genetics is a valuable tool in understanding the level of genetic variability present in vector populations, vector competence, and the development of novel control strategies. This study investigated the genetic differentiation between P. papatasi populations in Egypt and Jordan that inhabit distinct ecotopes and compared this structure to P. papatasi populations from a broader geographical range.

METHODS

A 461 base pair (bp) fragment from the mtDNA cytochrome b (cyt b) gene was PCR amplified and sequenced from 116 individual female sand flies from Aswan and North Sinai, Egypt, as well as Swaimeh and Malka, Jordan. Haplotypes were identified and used to generate a median-joining network, F _ST values and isolation-by-distance were also evaluated. Additional sand fly individuals from Afghanistan, Iran, Israel, Jordan, Libya, Tunisia and Turkey were included as well as previously published haplotypes to provide a geographically broad genetic variation analysis.

RESULTS

Thirteen haplotypes displaying nine variant sites were identified from P. papatasi collected in Egypt and Jordan. No private haplotypes were identified from samples in North Sinai, Egypt, two were observed in Aswan, Egypt, four from Swaimeh, Jordan and two in Malka, Jordan. The Jordan populations clustered separately from the Egypt populations and produced more private haplotypes than those from Egypt. Pairwise F _ST values fall in the range 0.024-0.648.

CONCLUSION

The clustering patterns and pairwise F _ST values indicate a strong differentiation between Egyptian and Jordanian populations, although this population structure is not due to isolation-by-distance. Other factors, such as environmental influences and the genetic variability in the circulating Le. major parasites, could possibly contribute to this heterogeneity. The present study aligns with previous reports in that pockets of genetic differentiation exists between populations of this widely dispersed species but, overall, the species remains relatively homogeneous.

Variability of Cutaneous Leishmaniasis Lesions Is Not Associated with Genetic Diversity of Leishmania tropica in Khyber Pakhtunkhwa Province of Pakistan

Leishmania tropica is the causative agent of cutaneous leishmaniasis in Pakistan. Here, intraspecific diversity of L. tropica from northern Pakistan was investigated using multilocus microsatellite typing. Fourteen polymorphic microsatellite markers were typed in 34 recently collected L. tropica isolates from Pakistan along with 158 archival strains of diverse Afro-Eurasian origins. Previously published profiles for 145 strains of L. tropica originating from different regions of Africa, Central Asia, Iran, and Middle East were included for comparison. Six consistently well-supported genetic groups were resolved: 1) Asia, 2) Morroco A, 3) Namibia and Kenya A, 4) Kenya B/Tunisia and Galilee, 5) Morocco B, and 6) Middle East. Strains from northern Pakistan were assigned to Asian cluster except for three that were placed in a geographically distant genetic group; Morocco A. Lesion variability among these Pakistani strains was not associated with specific L. tropica genetic profile. Pakistani strains showed little genetic differentiation from strains of Iraq, Afghanistan, and Syria (F_ST = 0.00-0.06); displayed evidence of modest genetic flow with India (F_ST = 0.14). Furthermore, genetic structuring within these isolates was not geographically defined. Pak-Afghan cluster was in significant linkage disequilibrium (I_A = 1.43), had low genetic diversity, and displayed comparatively higher heterozygosity (F_IS = -0.62). Patterns of genetic diversity observed suggest dominance of a minimally diverse clonal lineage within northern Pakistan. This is surprising as a wide clinical spectrum was observed in patients, suggesting the importance of host and other factors. Further genotyping studies of L. tropica isolates displaying different clinical phenotypes are required to validate this potentially important observation.

Mitochondrial genomes of two phlebotomine sand flies, Phlebotomus chinensis and Phlebotomus papatasi (Diptera: Nematocera), the first representatives from the family Psychodidae

BACKGROUND

Leishmaniasis is a worldwide but neglected disease of humans and animal transmitted by sand flies, vectors that also transmit other important diseases. Mitochondrial genomes contain abundant information for population genetic and phylogenetic studies, important in disease management. However, the available mitochondrial sequences of these crucial vectors are limited, emphasizing the need for developing more mitochondrial genetic markers.

METHODS

The complete mitochondrial genome of Phlebotomus chinensis was amplified in eight fragments and sequenced using primer walking. The mitochondrial genome of Phlebotomus papatasi was reconstructed from whole-genome sequencing data available on Genbank. The phylogenetic relationship of 24 selected representatives of Diptera was deduced from codon positions 1 and 2 for 13 protein coding genes, using Bayesian inference (BI) and maximum likelihood (ML) methods.

RESULTS

We provide the first Phlebotomus (P. chinensis and P. papatasi) mitochondrial genomes. Both genomes contain 13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes, and an A + T-rich region. The gene order of Phlebotomus mitochondrial genomes is identical with the ancestral gene order of insect. Phylogenetic analyses demonstrated that Psychodidae and Tanyderidae are sister taxa. Potential markers for population genetic study of Phlebotomus species were also revealed.

CONCLUSION

The generated mitochondrial genomes of P. chinensis and P. papatasi represent a useful resource for comparative genomic studies and provide valuable future markers for the population genetic study of these important Leishmania vectors. Our results also preliminary demonstrate the phylogenetic placement of Psychodidae based on their mitochondrial genomes.

Mitochondrial DNA Intraspecific Variability in Sergentomyia minuta (Diptera: Psychodidae)

Recently, there has been growing interest in analysis of the geographical variation between populations of different Phlebotomus spp. and American sand flies by comparing the sequences of various genes. However, little is known about the genetic structure of the genus Sergentomyia França & Parrot. No study has been carried out on Sergentomyia minuta Rondani. Most authors recognize this as a species with a high degree of morphological polymorphism, and some suspect that there are two subspecies: Se. minuta minuta Rondani in Europe, having about 40 horizontal cibarial teeth (sticks aligned along a straight line in the cibarial cavity), and Se. minuta parroti Adler & Theodor in North Africa, having about 70 cibarial teeth. Here we analyzed phylogeographic patterns using cytochrome b (Cytb) and cytochrome C oxidase I mtDNA for 29 populations from 10 countries: Algeria, Cyprus, France (continental and Corsica), Greece (continental and Crete), Malta, Montenegro, Morocco, Portugal (continental and Atlantic Savage Islands), Spain, and Tunisia. We analyzed intra- and interpopulation patterns of genetic diversity. Our results from Bayesian inference showed a complex genetic structure of Se. minuta with four haplogroups including many different haplotypes. One haplogroup includes all the specimens from North Africa. A second haplogroup includes a few specimens from the south of France, Spain, and one from Portugal. The third includes many specimens from southern France, all the specimens from Corsica, one from Spain, and all specimen from Portugal except one. A fourth branch includes specimens from the Balkans, Malta, Crete, Cyprus, and curiously some from the Atlantic Savage Islands; settlement of the latter population remains unexplained. However, our results suggest that the settlement of the Mediterranean basin could have occurred at the same time for Se. minuta and both Phlebotomus perniciosus Newstead and Phlebotomus ariasi Tonnoir. The spatial distribution of haplotypes was congruent with phylogenetic findings.

Paraphyly of the subgenus Anaphlebotomus and creation of Madaphlebotomus subg. nov. (Phlebotominae: Phlebotomus)

The systematic position of the Malagasy Phlebotomus (Diptera: Psychodidae) species was assessed in molecular phylogenetic studies. Three molecular markers were sequenced: cytochrome b of the mitochondrial DNA; ITS2, and the D8 domain of the ribosomal DNA. The following species were studied: Phlebotomus (Anaphlebotomus) berentiensis, Phlebotomus (Anaphlebotomus) fertei, Phlebotomus (Anaphlebotomus) fontenillei, Phlebotomus (Anaphlebotomus) vaomalalae and Phlebotomus (Anaphlebotomus) vincenti from Madagascar; Phlebotomus (Anaphlebotomus) stantoni from Asia, and Phlebotomus (Anaphlebotomus) rodhaini from Africa. The following outgroups were selected: Phlebotomus (Euphlebotomus) argentipes, Phlebotomus (Euphlebotomus) barguesae, Phlebotomus (Larroussius) perfiliewi s.l. and Phlebotomus (Adlerius) simici. Each marker analysed by maximum parsimony and maximum likelihood supports the monophyly of the Malagasy Phlebotomus spp. Consequently, we create a new subgenus for these species: Madaphlebotomus subg. nov. This molecular individualization is reinforced by the originality of their spermathecae and by the fact that their geographical distribution is limited to Madagascar, and considers the high level of endemism on this island.

Phlebotomus papatasi SP15: mRNA expression variability and amino acid sequence polymorphisms of field populations

Background

The Phlebotomus papatasi salivary protein PpSP15 was shown to protect mice against Leishmania major, suggesting that incorporation of salivary molecules in multi-component vaccines may be a viable strategy for anti-Leishmania vaccines.

Methods

Here, we investigated PpSP15 predicted amino acid sequence variability and mRNA profile of P. papatasi field populations from the Middle East. In addition, predicted MHC class II T-cell epitopes were obtained and compared to areas of amino acid sequence variability within the secreted protein.

Results

The analysis of PpSP15 expression from field populations revealed significant intra- and interpopulation variation.. In spite of the variability detected for P. papatasi populations, common epitopes for MHC class II binding are still present and may potentially be used to boost the response against Le. major infections.

Conclusions

Conserved epitopes of PpSP15 could potentially be used in the development of a salivary gland antigen-based vaccine.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-015-0914-2) contains supplementary material, which is available to authorized users.

Genetic divergence in populations of Lutzomyia ayacuchensis, a vector of Andean-type cutaneous leishmaniasis, in Ecuador and Peru

Haplotype and gene network analyses were performed on mitochondrial cytochrome oxidase I and cytochrome b gene sequences of Lutzomyia (Lu.) ayacuchensis populations from Andean areas of Ecuador and southern Peru where the sand fly species transmit Leishmania (Leishmania) mexicana and Leishmania (Viannia) peruviana, respectively, and populations from the northern Peruvian Andes, for which transmission of Leishmania by Lu. ayacuchensis has not been reported. The haplotype analyses showed higher intrapopulation genetic divergence in northern Peruvian Andes populations and less divergence in the southern Peru and Ecuador populations, suggesting that a population bottleneck occurred in the latter populations, but not in former ones. Importantly, both haplotype and phylogenetic analyses showed that populations from Ecuador consisted of clearly distinct clusters from southern Peru, and the two populations were separated from those of northern Peru.

Population genetics of Leishmania (Leishmania) major DNA isolated from cutaneous leishmaniasis patients in Pakistan based on multilocus microsatellite typing

Background

Cutaneous leishmaniasis (CL) is a major and fast increasing public health problem, both among the local Pakistani populations and the Afghan refugees in camps. Leishmania (Leishmania) major is one of the etiological agents responsible for CL in Pakistan. Genetic variability and population structure have been investigated for 66 DNA samples of L. (L.) major isolated from skin biopsy of CL patients.

Methods

Multilocus microsatellite typing (MLMT), employing 10 independent genetic markers specific to L. (L.) major, was used to investigate the genetic polymorphisms and population structures of Pakistani L. (L.) major DNA isolated from CL human cases. Their microsatellite profiles were compared to those of 130 previously typed strains of L. (L.) major from various geographical localities.

Results

All the markers were polymorphic and fifty-one MLMT profiles were recognized among the 66 L. (L.) major DNA samples. The data displayed significant microsatellite polymorphisms with rare allelic heterozygosities. A Bayesian model-based approach and phylogenetic analysis inferred two L. (L.) major populations in Pakistan. Thirty-four samples belonged to one population and the remaining 32 L. (L.) major samples grouped together into another population. The two Pakistani L. (L.) major populations formed separate clusters, which differ genetically from the populations of L. (L.) major from Central Asia, Iran, Middle East and Africa.

Conclusions

The considerable genetic variability of L. (L.) major might be related to the existence of different species of sand fly and/or rodent reservoir host in Sindh province, Pakistan. A comprehensive study of the epidemiology of CL including the situation or spreading of reservoirs and sand fly vectors in these foci is, therefore, warranted.

Genetic diversity of Leishmania donovani/infantum complex in China through microsatellite analysis

The Leishmania strains from different epidemic areas in China were assessed for their genetic relationship. Twenty-nine strains of Leishmania infantum isolated from 1950 to 2001 were subjected to multilocus microsatellite typing (MLMT) using 14 highly polymorphic microsatellite markers. Twenty-two MLMT profiles were recognized among the 29 L. infantum strains, which differed from one another in 13 loci. Bayesian model-based and distance-based analysis of the data inferred two main populations in China. Sixteen strains belonged to one population, which also comprised previously characterized strains of L. infantum non-MON1 and Leishmania donovani. The parasites within this population are assignable to a distinct cluster that is clearly separable from the populations of L. donovani elsewhere, i.e. India, Sri Lanka and East Africa, and L. infantum non-MON1 from Europe. The remaining 13 Chinese strains grouped together with strains of L. infantum MON1 into another population, but formed a separate cluster which genetically differs from the populations of L. infantum MON1 from Europe, the Middle East, Central Asia and North Africa. The existence of distinct groups of L. infantum MON1 and non-MON1/L. donovani suggests that the extant parasites in China may have been restricted there, but not recently introduced from elsewhere.

Phylogeography of the Lutzomyia gomezi (Diptera: Phlebotominae) on the Panama Isthmus

BACKGROUND

Lutzomyia gomezi (Nitzulescu, 1931) is one of the main Leishmania (Vianna) panamensis vectors in Panama, and despite its medical significance, there are no population genetic studies regarding this species. In this study, we used the sequences of the mitochondrial gene cytochrome b/start of NADH1 and the nuclear elongation gene α-1 in order to analyze genetic variation and phylogeographic structure of the Lu. gomezi populations.

METHODS

A total of 86 Lu. gomezi individuals were captured in 38 locations where cutaneous leishmaniasis occurred. DNA was extracted with phenol/chloroform methods and amplification of genes was performed using PCR primers for mitochondrial and nuclear markers.

RESULTS

We found a total of 37 and 26 haplotypes of mitochondrial and nuclear genes, high haplotype diversity (h) for all three populations were detected with both molecular markers. Nucleotide diversity (π) was estimated to be high for all three populations with the mitochondrial marker, which was opposite to the estimate with the nuclear marker. In the AMOVA Φst recorded moderate (mitochondrial) and small (nuclear) population structure with statistical significance among populations. The analysis of the fixation index (Fst) used to measure the differentiation of populations showed that with the exception of the population located in the region of Bocas del Toro, the other populations presented with minor genetic differentiation. The median-Joining network of the mitochondrial marker reveled three clusters and recorded four haplotypes exclusively of localities sampled from Western Panama, demonstrating strong divergence. We found demographic population expansion with Fu´s Fs neutrality test. In the analysis mismatch distribution was observed as a bimodal curve.

CONCLUSION

Lu. gomezi is a species with higher genetic pool or variability and mild population structure, due to possible capacity migration and local adaptation to environmental changes or colonization potential. Thus, knowledge of the genetic population and evolutionary history is useful to understand the implications of different population genetic structures for cutaneous leishmaniasis epidemiology.

Molecular and morphological systematics of the sandfly Sergentomyia (Sintonius) clydei Sinton, 1928 and questions about its record in the Seychelles

In the Phlebotomine sandflies, a few molecular studies related on the genus Sergentomyia have been published. The present study explored the genetic variability within Sergentomyia (Sintonius) clydei (Diptera, Psychodidae). The sampling included 15 populations originating from 12 countries. A morphological approach was coupled to the sequencing of two molecular markers (cytochrome b mtDNA and cacophony nuclear DNA). The most variable morphological characters resided in the cibarium of the females, especially (i) the pigment patch pattern and (ii) the number of cibarial teeth and denticles in the armature. However this morphological approach was unable to individualize any population within S. clydei. The NJ trees based on both molecular markers individualized the specimens from the Aldabra group of islands in the Seychelles. Surprisingly, cyt b variability was not compatible with the known data about the complete submersion of Aldabra occurring relatively recently some 125,000 years ago. The settlement of these islands by S. clydei from continental Africa, the Middle East or Asia, and the value of mtDNA markers are discussed.

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.

Genetic diversity of the mitochondrial cytochrome b gene in Lutzomyia spp., with special reference to Lutzomyia peruensis, a main vector of Leishmania (Viannia) peruviana in the Peruvian Andes

The genetic divergence caused by genetic drift and/or selection is suggested to affect the vectorial capacity and insecticide susceptibility of sand flies, as well as other arthropods. In the present study, cytochrome b (cyt b) gene sequences were determined in 13 species circulating in Peru to establish a basis for analysis of the genetic structure, and the intraspecific genetic diversity was assessed in the Lutzomyia (Lu.) peruensis, a main vector species of Leishmania (Viannia) peruviana in Peruvian Andes. Analysis of intraspecific genetic diversity in the cyt b gene sequences from 36 Lu. peruensis identified 3 highly polymorphic sites in the middle region of the gene. Haplotype and gene network analyses were performed on the cyt b gene sequences of 130 Lu. peruensis in 9 Andean areas from 3 Departments (Ancash, Lima and La Libertad). The results showed that the populations of La Libertad were highly polymorphic and that their haplotypes were distinct from those of Ancash and Lima, where dominant haplotypes were observed, suggesting that a population bottleneck may have occurred in Ancash and Lima, but not in La Libertad. The present study indicated that the middle region of the cyt b gene is useful for the analysis of genetic structure in sand fly populations.

Genetic differentiation between sandfly populations of Phlebotomus chinensis and Phlebotomus sichuanensis (Diptera: Psychodidae) in China inferred by microsatellites

Background

Phlebotomus chinensis is a primary vector of visceral leishmaniasis; it occurs in various biotopes with a large geographical distribution, ranging from Yangtze River to northeast China. Phlebotomus sichuanensis, a species closely related to P. chinensis in high altitude regions, has a long term disputation on its taxonomic status. Both species occur in the current epidemic regions and are responsible for the transmission of leishmaniasis. Population genetic analysis will help to understand the population structure and infer the relationship for morphologically indistinguishable cryptic species. In this study, microsatellite markers were used for studying the genetic differentiation between P. chinensis and P. sichuanensis.

Methods

Sandflies were collected in 6 representative localities in China in 2005-2009. Ten microsatellite loci were used to estimate population genetic diversity. The intra-population genetic diversity, genetic differentiation and effective population size were estimated.

Results

All 10 microsatellite loci were highly polymorphic across populations, with high allelic richness and heterozygosity. Hardy-Weinberg disequilibrium was found in 23 out of 60 (38.33%) comparisons associated with heterozygote deficits, which was likely caused by the presence of null allele and the Wahlund effect. Bayesian clustering analysis revealed three clusters. The cluster I included almost all specimens in the sample SCD collected at high altitude habitats in Sichuan. The other two clusters were shared by the remaining 5 populations, SCJ in Sichuan, GSZ in Gansu, SXL and SXX in Shaanxi and HNS in Henan. The diversity among these 5 populations was low (F_ST = -0.003-0.090) and no isolation by distance was detected. AMOVA analysis suggested that the variations were largely derived from individuals within populations and among individuals. Consistently, the analysis of ribosomal DNA second internal transcribed spacer (ITS2) sequence uncovered three types of variants, which corresponded with the three gene pools revealed by microsatellites.

Conclusions

The data suggested that the SCD population carried a distinct gene pool, which was differentiated from the other populations. The high altitude ecological habitats, distinctive ITS2 and herein divergence inferred by microsatellite loci support the species status of P. sichuanensis. The P. chinensis populations did not have a significant divergence from each another.

Genetic structure and divergence in populations of Lutzomyia cruciata, a phlebotomine sand fly (Diptera: Psychodidae) vector of Leishmania mexicana in southeastern Mexico

The low dispersal capacity of sand flies could lead to population isolation due to geographic barriers, climate variation, or to population fragmentation associated with specific local habitats due to landscape modification. The phlebotomine sand fly Lutzomyia cruciata has a wide distribution throughout Mexico and is a vector of Leishmania mexicana in the southeast. The aim of this study was to evaluate the genetic diversity, structure, and divergence within and among populations of Lu. cruciata in the state of Chiapas, and to infer the intra-specific phylogeny using the 3' end of the mitochondrial cytochrome b gene. We analyzed 62 sequences from four Lu. cruciata populations and found 26 haplotypes, high genetic differentiation and restricted gene flow among populations (Fst=0.416, Nm=0.701, p<0.001). The highest diversity values were recorded in populations from Loma Bonita and Guadalupe Miramar. Three lineages (100% bootstrap and 7% overall divergence) were identified using a maximum likelihood phylogenetic analysis which showed high genetic divergence (17.2-22.7%). A minimum spanning haplotype network also supported separation into three lineages. Genetic structure and divergence within and among Lu. cruciata populations are hence affected by geographic heterogeneity and evolutionary background. Data obtained in the present study suggest that Lu. cruciata in the state of Chiapas consists of at least three lineages. Such findings may have implications for vector capacity and hence for vector control strategies.

Evidence for genetic differentiation at the microgeographic scale in Phlebotomus papatasi populations from Sudan

Background

Cutaneous Leishmaniasis (CL) is endemic in Sudan. It is caused by Leishmania major parasites and transmitted by Phlebotomus papatasi sandflies. Recently, uncommon clinical manifestations of CL have been reported. Moreover, L. donovani parasites that cause Visceral Leishmaniasis (VL) have been isolated from CL lesions of some patients who contracted the disease in Khartoum State, Central Sudan with no history of travelling to VL endemic sites on south-eastern Sudan. Because different clinical manifestations and the parasite behaviour could be related to genetic differentiation, or even sub-structuring within sandfly vector populations, a population genetic study was conducted on P. papatasi populations collected from different localities in Khartoum State known for their uncommon CL cases and characterized by contrasting environmental conditions.

Methods

A set of seven microsatellite loci was used to investigate the population structure of P. papatasi samples collected from different localities in Khartoum State, Central Sudan. Populations from Kassala State, Eastern Sudan and Egypt were also included in the analyses as outgroups. The level of genetic diversity and genetic differentiation among natural populations of P. papatasi was determined using F_ST statistics and Bayesian assignments.

Results

Genetic analyses revealed significant genetic differentiation (F_ST) between the Sudanese and the Egyptian populations. Within the Sudanese P. papatasi populations, one population from Gerif West, Khartoum State, exhibited significant genetic differentiation from all other populations including those collected as near as 22 km.

Conclusion

The significant genetic differentiation of Gerif West P. papatasi population from other Sudanese populations may have important implication for the epidemiology of leishmaniasis in Khartoum State and needs to be further investigated. Primarily, it could be linked to the unique location of Gerif West which is confined by the River Nile and its tributaries that may act as a natural barrier for gene flow between this site and the other rural sites. The observed high migration rates and lack of genetic differentiation among the other P. papatasi populations could be attributed to the continuous human and cattle movement between these localities.

Expression plasticity of Phlebotomus papatasi salivary gland genes in distinct ecotopes through the sand fly season

BACKGROUND

Sand fly saliva can drive the outcome of Leishmania infection in animal models, and salivary components have been postulated as vaccine candidates against leishmaniasis. In the sand fly Phlebotomus papatasi, natural sugar-sources modulate the activity of proteins involved in meal digestion, and possibly influence vectorial capacity. However, only a handful of studies have assessed the variability of salivary components in sand flies, focusing on the effects of environmental factors in natural habitats. In order to better understand such interactions, we compared the expression profiles of nine P. papatasi salivary gland genes of specimens inhabiting different ecological habitats in Egypt and Jordan and throughout the sand fly season in each habitat.

RESULTS

The majority of investigated genes were up-regulated in specimens from Swaymeh late in the season, when the availability of sugar sources is reduced due to water deprivation. On the other hand, these genes were not up-regulated in specimens collected from Aswan, an irrigated area less susceptible to drought effects.

CONCLUSION

Expression plasticity of genes involved with vectorial capacity in disease vectors may play an important epidemiological role in the establishment of diseases in natural habitats.

Multilocus molecular and phylogenetic analysis of phlebotomine sand flies (Diptera: Psychodidae) from southern Italy

This study reports a combined analysis of mitochondrial and ribosomal DNA target regions of phlebotomine sand flies (Diptera: Psychodidae) from the Mediterranean region. A ∼900 bp long fragment of the mitochondrial DNA encompassing regions within cytb and nd1 gene and the complete ITS2 ribosomal region (∼500 bp) were sequenced and characterized for Phlebotomus perniciosus, Phlebotomus perfiliewi, Phlebotomus neglectus, Phlebotomus papatasi, and Sergentomyia minuta, captured in two sites of southern Italy. From one to eight mitochondrial haplotypes and from one to three ITS2 sequence types were found for the examined specimens according to the different sand fly species. The mean interspecific difference in the mitochondrial sequences was of 16.1%, with an overall intraspecific nucleotide variation from 0.1 to 2.8%. A higher interspecific difference (mean 25.1%) was recorded for the ITS2 sequence, with an overall intraspecific nucleotide variation up to 4.9%. The sequence types alignment of ITS2 region showed that all phlebotomine specimens possessed a split 5.8S rRNA, consisting of a mature 5.8S rRNA and a 2S rRNA separated by a short transcribed spacer. Phylogenetic analysis of the Phlebotomus spp. sequences, herein determined and of those available in GenBank™ were concordant in clustering P. neglectus, P. perfiliewi and P. papatasi with the same species collected from different geographic areas of the Mediterranean basin in four main clades for mtDNA and ITS2, respectively. This study demonstrates the utility of multilocus sequencing, provides a dataset for the molecular identification of the most prevalent phlebotomine sand flies in southern Europe and defines the phylogenetic relationships among species examined.

Intraspecific variability of natural populations of Phlebotomus sergenti, the main vector of Leishmania tropica

An intraspecific study of Phlebotomus sergenti was performed on populations from Turkey, Syria, Israel, and Uzbekistan by four different approaches: geometric morphometrics, RAPD analysis, internal transcribed spacer 2 (ITS2) sequencing (nuclear marker), and cytochrome B sequencing (mitochondrial marker). In RAPD analysis, distinct clades were formed in accordance with the geographical origin of the specimens. There was no distinct grouping according to place of origin within the Turkish samples from various localities in south-eastern Anatolia, which suggests a gene flow between populations separated spatially by the Amanos mountains, a mountain range of a considerable altitude. The results of ITS2 rDNA sequencing complied with the previously published intraspecific division of P. sergenti into two branches, northeastern and southwestern. However, mtDNA haplotypes formed three lineages with specimens from Turkey and Israel, sharing a common clade. A previously postulated hypothesis about a complex of sibling species within P. sergenti is therefore questionable. Cytochrome B seems to be a more discriminative marker for intraspecific variability assessment.

The sandflies (Diptera: Psychodidae, Phlebotominae) in military camps in northern Afghanistan (2007-2009), as identified by morphology and DNA 'barcoding'

As part of a continuous, standardized programme of monitoring the Leishmania vectors in German military camps in northern Afghanistan between 2007 and 2009, a detailed taxonomic analysis of the endemic sandfly fauna, as sampled using light and odour-baited traps, was conducted. Of the 10 sandfly species that were recorded, six may serve as enzootic and/or zooanthroponotic vectors of parasites causing human leishmaniasis. The use of a simple DNA-'barcoding' technique based on the mitochondrial cyt b gene, to identify the collected sandflies to species level, revealed (1) a clear discrimination between the potential vector species, (2) clustering of species within most subgenera, and (3) particularly high heterogeneity within the subgenus Paraphlebotomus (Phlebotomus alexandri being grouped with Ph. papatasi rather than with other Paraphlebotomus species). The data also indicate a high level of genetic heterogeneity within the subgenus Sergentomyia but close similarity between Sergentomyia sintoni and Sergentomyia murgabiensis. The morphological similarity of many medically important sandflies can make species identification difficult, if not impossible. The new DNA-barcoding techniques may provide powerful discriminatory tools in the future.

Ecological genomics of sand fly salivary gland genes: an overview

Sand fly saliva contains an array of bioactive molecules that facilitate blood feeding and also function as modulators of the vertebrate immune response. Such a complex of biologically active molecules was shown to be both conserved and divergent among sand fly species. It is likely that expression of sand fly salivary molecules could be modulated by environmental factors, both biotic and abiotic, that ultimately dictate the quality, and possibly quantity, of the secreted saliva. Carbohydrates are an integral part of the sand fly diet, and sugar-sources found in natural habitats are potentially involved in defining the profile of sand fly saliva, and may influence vectorial capacity. Saliva can drive the outcome of Leishmania infection in animal models, and salivary molecules are potential targets for development of vaccines to control Leishmania infection. Thus, identifying what environmental factors effectively modulate sand fly saliva in the field is a critical step towards the development of meaningful protection strategies against leishmaniasis that are based on salivary compounds from sand fly vectors.

Genetic structures of sand fly (Diptera: Psychodidae) populations in a leishmaniasis endemic region of Turkey

The object of this study was to determine the genetic structures of three vector species, Phlebotomus tobbi, Phlebotomus papatasi, and Phlebotomus sergenti, in the Cukurova Region of Turkey, an endemic focus of cutaneous leishmaniasis. The genetic diversity indices, neutrality tests and hierarchical analysis of molecular variance (AMOVA) were performed using partial sequences of ITS2 and cytochrome b gene regions. In all species, within population genetic variation was higher than between population variation for ITS2 gene region. Fst values were low and non-significant for P. sergenti, and were higher for P. papatasi and P. tobbi indicating a weak structuring between populations. AMOVA tests suggest any substantial isolation between populations within species. AMOVA analysis of cyt b gene region revealed significant genetic structuring between populations for P. papatasi and P. sergenti. Fst values were relatively high and significant for these species indicating a certain degree of isolation between populations. However, in P. tobbi, any significant population genetic structuring was detected. Tajima's D and Fu's Fs values were negative and significant in all three species might be indicating a demographic expansion.

Distribution of Leishmania major zymodemes in relation to populations of Phlebotomus papatasi sand flies

Phlebotomus papatasi (Scopoli) (Diptera: Psychodidae) is the main vector of Leishmania major Yakimoff & Schokhor (Kinetoplastida: Trypanosomatidae), the causative agent of zoonotic cutaneous leishmaniasis in the Old World. Multilocus enzyme electrophoresis (MLEE) was extensively used to type different L. major stocks allover the world. Multilocus microsatellite typing (MLMT) has been recently used to investigate P. papatasi sand flies at population and subpopulation levels. In this article, the association between geographical distribution of L. major zymodemes and the distribution of populations and subpopulations of L. major vector; P. papatasi are discussed.

Occurrence and genetic variability of Phlebotomus papatasi in an urban area of southern Italy

BACKGROUND

A phlebotomine sand fly was noticed in the second floor of an old building in a highly urbanized area of southern Italy. A short-term entomological survey was carried out in the subsequent weeks to this event, allowing the collection of additional phlebotomine sand flies that were later identified as Phlebotomus papatasi. We assessed the genetic variability among P. papatasi sequences obtained in this study and those available from Italy using a mitochondrial DNA (mtDNA) fragment (from cytochrome b gene to NADH1) and the internal transcribed spacer 2 (ITS2) as genetic markers.

RESULTS

From 9 June to 19 July, eight males and seven females (two blood-fed) of P. papatasi were collected in the old town of Bari (southern Italy). The insects were found near the bed and in the bathroom and potential blood sources (e.g., pigeons and dogs) for them were common in the neighbourhood. Again, five females of P. papatasi collected in Valenzano, another urban area in the province of Bari, were also identified and included in the genetic study. The mtDNA sequences (945 bp) obtained from Bari and Valenzano were identical except for a single transition (T ↔ C) at the 793 nucleotide residue. Pairwise comparison of the last 440 bp of the mtDNA fragment analyzed herein with other sequences of P. papatasi from Italy revealed a nucleotide variation ranging from 0.2 to 1.3%. Three ITS2 sequence types were detected within specimens collected in Valenzano, one of them identical to that from Bari. Pairwise comparison of ITS2 sequences of P. papatasi from Italy revealed a nucleotide variation up to 1.8%.

CONCLUSIONS

This study reports the occurrence of P. papatasi in an urban area of southern Italy and shows a low nucleotide difference among ITS2 and mtDNA sequences of this species available from Italy. The presence of P. papatasi in urban areas might represent a risk for human health, particularly for the potential transmission of sandfly fever viruses.

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.

Identification of geographically distributed sub-populations of Leishmania (Leishmania) major by microsatellite analysis

BACKGROUND

Leishmania (Leishmania) major, one of the agents causing cutaneous leishmaniasis (CL) in humans, is widely distributed in the Old World where different species of wild rodent and phlebotomine sand fly serve as animal reservoir hosts and vectors, respectively. Despite this, strains of L. (L.) major isolated from many different sources over many years have proved to be relatively uniform. To investigate the population structure of the species highly polymorphic microsatellite markers were employed for greater discrimination among it's otherwise closely related strains, an approach applied successfully to other species of Leishmania.

RESULTS

Multilocus Microsatellite Typing (MLMT) based on 10 different microsatellite markers was applied to 106 strains of L. (L.) major from different regions where it is endemic. On applying a Bayesian model-based approach, three main populations were identified, corresponding to three separate geographical regions: Central Asia (CA); the Middle East (ME); and Africa (AF). This was congruent with phylogenetic reconstructions based on genetic distances. Re-analysis separated each of the populations into two sub-populations. The two African sub-populations did not correlate well with strains' geographical origin. Strains falling into the sub-populations CA and ME did mostly group according to their place of isolation although some anomalies were seen, probably, owing to human migration.

CONCLUSION

The model- and distance-based analyses of the microsatellite data exposed three main populations of L. (L.) major, Central Asia, the Middle East and Africa, each of which separated into two sub-populations. This probably correlates with the different species of rodent host.

Isoenzymatic characterization of Phlebotomus papatasi (Diptera: Psychodidae) of the Marrakech area, Morocco

This study reports the genetic characterization of urban and rural populations of Phlebotomus (Phlebotomus) papatasi Scopoli (Diptera: Psychodidae) in Marrakech, Morocco. Using isoenzymatic analysis, four Moroccan populations were compared with other Mediterranean basin populations from Spain, Cyprus, and Syria. Morphological anomalies were noted in the male genitalia of 5.3% of the specimens collected from Marrakech area. Qualitative analysis of zymogram profiles revealed nine polymorphic enzymes (HK, PGM, PGI, 6PGD, MDH1, MDH2, ICD2, FUM, and ACO) and three monomorphic enzymes (ME, ICD1, and alphaGPDH). Genetic distances clearly separated the populations of western Mediterranean countries (Morocco and Spain) from eastern countries (Syria and Cyprus), but they could not be used to differentiate between urban and rural populations in Marrakech area.