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

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.

Multiple evolutionary lineages for the main vector of Leishmania guyanensis, Lutzomyia umbratilis (Diptera: Psychodidae), in the Brazilian Amazon

Lutzomyia umbratilis is the main vector of Leishmania guyanensis in the Brazilian Amazon and in neighboring countries. Previous biological and molecular investigations have revealed significant differences between L. umbratilis populations from the central Brazilian Amazon region. Here, a phylogeographic survey of L. umbratilis populations collected from nine localities in the Brazilian Amazon was conducted using two mitochondrial genes. Statistical analyses focused on population genetics, phylogenetic relationships and species delimitations. COI genetic diversity was very high, whereas Cytb diversity was moderate. COI genealogical haplotypes, population structure and phylogenetic analyses identified a deep genetic differentiation and three main genetic groups. Cytb showed a shallower genetic structure, two main haplogroups and poorly resolved phylogenetic trees. These findings, allied to absence of isolation by distance, support the hypothesis that the Amazon and Negro Rivers and interfluves are the main evolutionary forces driving L. umbratilis diversification. The main three genetic groups observed represent three evolutionary lineages, possibly species. The first lineage occurs north of the Amazon River and east of Negro River, where Le. guyanensis transmission is intense, implying that L. umbratilis is an important vector there. The second lineage is in the interfluve between north of Amazon River and west of Negro River, an area reported to be free of Le. guyanensis transmission. The third lineage, first recorded in this study, is in the interfluve between south of Amazonas River and west of Madeira River, and its involvement in the transmission of this parasite remains to be elucidated.

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.

Leishmaniasis Beyond East Africa

Climate change is having a substantial impact on our environment and ecosystems and has altered the way humans live, access, and utilize resources with increased risk of zoonotic infectious disease encounters. As global temperatures continue to increase, they impact on public health, migration, food security and land conflict, and as new environments become favorable, exposure to disease carrying vectors. Increased forests or natural habitat clearance for land repurposing, urbanization, road building, and water management are related to an increase in emerging vector borne parasitic diseases. The East African region remains one of the most impacted regions globally for leishmaniasis, a vector borne disease that impacts significantly on the health, wellbeing and livelihoods of affected communities and for which a lack of reporting and control interventions hinder progress toward elimination of this neglected tropical disease. As our world continues to transform, both politically and climatically, it is essential that measures are put in place to improve surveillance and disease management with implementation of control measures, including vector control, especially in low- and middle-income countries that are expected to be most impacted by changes in climate. Only through effective management, now, can we be sufficiently resilient to preventing the inevitable spread of vectors into suitable habitat and expansion of the geographical range of leishmaniasis. This review offers a current perspective on Leishmaniasis as an endemic disease in East Africa and examines the potential of the recent emergence of Leishmania infection in hitherto unaffected regions to become a public health concern if no disease management is achieved.

Studies of host preferences of wild-caught Phlebotomus orientalis and Ph. papatasi vectors of leishmaniasis in Sudan

Introduction

Understanding the feeding behavior and host choice of sand flies provides valuable information on vector-host relationships and elucidates the epidemiological patterns of leishmaniasis transmission. Blood meal analysis studies are essential for estimating the efficiency of pathogen transmission, assessing the relative human disease risk, and assist in identifying the other potential hosts of leishmaniasis. In Sudan and most of East Africa, there are large remaining gaps in knowledge regarding the feeding habits of phlebotomine vectors. The study aimed to identify the blood meal sources and, therefore, the host preferences of the principal vectors Phlebotomus orientalis and Ph. papatasi in leishmaniasis endemic areas of eastern and central Sudan.

Materials and methods

Sand flies were collected from two endemic villages in eastern and central Sudan using CDC light traps and sticky traps. The phlebotomine sand flies were morphologically and then molecularly identified. The source of blood meal of the engorged females was determined using a multiplex PCR methodology and specific primers of cytochrome b gene of mitochondrial DNA for human, goat, cow, and dog. The detection of the Leishmania parasite was done using PCR.

Results

The total number of collected female phlebotomine sand flies was 180. Morphological identification revealed the abundance of Ph. orientalis 103 (57.2%), Ph. papatasi 42 (23.3%), Ph. bergeroti 31 (17.2%), Ph. rodhaini 2 (1.1%) and Ph. duboscqi 2 (1.1%) in the study sites. Out of the 180 collected, 31 (17%) were blood-fed flies. Three species were blood-fed and molecularly identified: Ph. papatasi (N = 7, 22.6%), Ph. bergeroti (N = 9, 26%), and Ph. orientalis (N = 15, 48.4%). Blood meal analysis revealed human DNA in two Ph. orientalis (6.4%), hence, the anthropophilic index was 13.3%.

Conclusions

Multiplex PCR protocol described here allowed the identification of blood meal sources of many vertebrate species simultaneously. The results indicate that wild-caught Ph. orientalis are anthropophilic in the study areas. Further studies on larger blood-fed sample size are required to validate the potential applications of this technique in designing, monitoring and evaluating control programs, particularly in investigating the potential non-human hosts of leishmaniasis.

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.

Ecological niche modeling predicting the potential distribution of Leishmania vectors in the Mediterranean basin: impact of climate change

Background

Due to climate change, the geographical distribution of sand flies during the last decades has shifted northward from latitudes below 45°N in southern Europe to latitudes just above 50^○N. Recent studies show that some phlebotomine sand flies were recorded in several parts of Germany and Belgium. In central Europe, some autochthone leishmaniasis cases are being recorded in regions traditionally regarded as leishmaniasis-free. An important challenge is to predict the geographical distribution of leishmaniasis vectors under new climatic conditions. In this study, we attempted to predict the current distribution of six leishmaniasis vectors in the Mediterranean basin and forecast species’ geographical shift under future climate scenarios using an ensemble ecological niche modeling approach. Species records were obtained from scientific surveys published in the research literature between 2006 and 2016. A series of climate metrics describing temperature and precipitation in the study area under two climatic scenarios were obtained from WorldClim database. A consensus model was derived from six varieties of modeling approaches (regression, machine learning and classification techniques) in order to ensure valid prediction of distribution of vectors under different climate scenarios.

Results

Model performance was generally high for the included species with a specificity (true negative rate) ranging from 81.03 to 96.52% (mean = 86.94%) and a sensitivity (true positive rate) ranging from 87.93 to 100% (mean = 96.98%). Our work evidenced the hypothesis of the widespread of Leishmania vectors under climate change scenarios. All of the studied species are prospected to gain new areas that are actually not suitable for vectors’ survival. Phlebotomine sand flies are prospected to invade extra-Mediterranean regions, especially western and central Europe.

Conclusions

Our study confirmed the importance of environmental and climate factors on the distribution of leishmaniasis vectors and demonstrated the performance of ecological niche modeling in the prediction of the geographical spread of vector-borne diseases. Ecological niche modeling should be considered in the future as a valuable tool in addition to experimental laboratory studies for a better understanding of the biology of vector species.

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.

Genetic differentiation over a small spatial scale of the sand fly Lutzomyia vexator (Diptera: Psychodidae)

Background

The geographic scale and degree of genetic differentiation for arthropod vectors that transmit parasites play an important role in the distribution, prevalence and coevolution of pathogens of human and wildlife significance. We determined the genetic diversity and population structure of the sand fly Lutzomyia vexator over spatial scales from 0.56 to 3.79 km at a study region in northern California. The study was provoked by observations of differentiation at fine spatial scales of a lizard malaria parasite vectored by Lu. vexator.

Methods

A microsatellite enrichment/next-generation sequencing protocol was used to identify variable microsatellite loci within the genome of Lu. vexator. Alleles present at these loci were examined in four populations of Lu. vexator in Hopland, CA. Population differentiation was assessed using Fst and D (of Cavalli-Sforza and Edwards), and the program Structure was used to determine the degree of subdivision present. The effective population size for the sand fly populations was also calculated.

Results

Eight microsatellite markers were characterized and revealed high genetic diversity (uHe = 0.79–0.92, Na = 12–24) and slight but significant differentiation across the fine spatial scale examined (average pairwise D = 0.327; F_ST = 0.0185 (95 % bootstrapped CI: 0.0102–0.0264). Even though the insects are difficult to capture using standard methods, the estimated population size was thousands per local site.

Conclusions

The results argue that Lu. vexator at the study sites are abundant and not highly mobile, which may influence the overall transmission dynamics of the lizard malaria parasite, Plasmodium mexicanum, and other parasites transmitted by this species.

Electronic supplementary material

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

New microsatellite markers for multi-scale genetic studies on Phlebotomus ariasi Tonnoir, vector of Leishmania infantum in the Mediterranean area

The population structure of Phlebotomus ariasi, a proven vector of Leishmania infantum in the Mediterranean area, is still poorly understood. Previously, only two microsatellite loci had been developed to study the population genetics of this species. Herein we use these loci and determined fourteen novel microsatellite loci, useful for the characterization of P. ariasi populations. These loci were tested on three populations of P. ariasi, two from France and one from Portugal. In addition, the usefulness of these markers was also evaluated on seven other sandfly species. We show, that for P. ariasi, 15 of the 16 loci selected were polymorphic, with a mean of 4.25 alleles and an observed heterozygosity of 0.299. Within the P. ariasi population of France, 11 loci were polymorphic, with an average of 2.44 alleles and an observed heterozygosity of 0.2177. The fixation index was moderate among the French populations but high between French and Portuguese populations. In addition, eight loci were also found to be amplifiable in six other Phlebotomus species. These results demonstrate the usefulness of this new set of microsatellite loci for population structure and molecular ecology studies of P. ariasi at various spatial scales, but also of other sandfly species.

Concomitant Infection with Leishmania donovani and L. major in Single Ulcers of Cutaneous Leishmaniasis Patients from Sudan

In Sudan human leishmaniasis occurs in different clinical forms, that is, visceral (VL), cutaneous (CL), mucocutaneous (ML), and post-kala-azar dermal leishmaniasis (PKDL). Clinical samples from 69 Sudanese patients with different clinical manifestations were subjected to a PCR targeting the cytochrome oxidase II (COII) gene for Leishmania species identification. Mixed infections were suspected due to multiple overlapping peaks presented in some sequences of the COII amplicons. Cloning these amplicons and alignment of sequences from randomly selected clones confirmed the presence of two different Leishmania species, L. donovani and L. major, in three out of five CL patients. Findings were further confirmed by cloning the ITS gene. Regarding other samples no significant genetic variations were found in patients with VL (62 patients), PKDL (one patient), or ML (one patient). The sequences clustered in a single homogeneous group within L. donovani genetic group, with the exception of one sequence clustering with L. infantum genetic group. Findings of this study open discussion on the synergetic/antagonistic interaction between divergent Leishmania species both in mammalian and vector hosts, their clinical implications with respect to parasite fitness and response to treatment, and the route of transmission with respect to vector distribution and or adaptation.

Molecular epidemiology and genetic diversity of Entamoeba species in a chelonian collection

Veterinary medicine has focused recently on reptiles, due to the existence of captive collections in zoos and an increase in the acquisition of reptiles as pets. The protozoan parasite, Entamoeba can cause amoebiasis in various animal species and humans. Although amoebiasis disease is remarkably rare in most species of chelonians and crocodiles, these species may serve as Entamoeba species carriers that transmit parasites to susceptible reptile species, such as snakes and lizards, which can become sick and die. In this study, we identified the Entamoeba species in a population of healthy (disease-free) chelonians, and evaluated their diversity through the amplification and sequencing of a small subunit rDNA region. Using this procedure, three Entamoeba species were identified: Entamoeba invadens in 4.76 % of chelonians, Entamoeba moshkovskii in 3.96 % and Entamoeba terrapinae in 50 %. We did not detect mixed Entamoeba infections. Comparative analysis of the amplified region allowed us to determine the intra-species variations. The E. invadens and E. moshkovskii strains isolated in this study did not exhibit marked differences with respect to the sequences reported in GenBank. The analysis of the E. terrapinae isolates revealed three different subgroups (A, B and C). Although subgroups A and C were very similar, subgroup B showed a relatively marked difference with respect to subgroups A and C (Fst = 0.984 and Fst = 1.000, respectively; 10-14 % nucleotide variation, as determined by blast) and with respect to the sequences reported in GenBank. These results suggested that E. terrapinae subgroup B may be either in a process of speciation or belong to a different lineage. However, additional research is necessary to support this statement conclusively.

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.