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

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.

Shifts in the ecological niche of Lutzomyia peruensis under climate change scenarios in Peru

The Peruvian Andes presents a climate suitable for many species of sandfly that are known vectors of leishmaniasis or bartonellosis, including Lutzomyia peruensis (Diptera: Psychodidae), among others. In the present study, occurrences data for Lu. peruensis were compiled from several items in the scientific literature from Peru published between 1927 and 2015. Based on these data, ecological niche models were constructed to predict spatial distributions using three algorithms [Support vector machine (SVM), the Genetic Algorithm for Rule-set Prediction (GARP) and Maximum Entropy (MaxEnt)]. In addition, the environmental requirements of Lu. peruensis and three niche characteristics were modelled in the context of future climate change scenarios: (a) potential changes in niche breadth; (b) shifts in the direction and magnitude of niche centroids, and (c) shifts in elevation range. The model identified areas that included environments suitable for Lu. peruensis in most regions of Peru (45.77%) and an average altitude of 3289 m a.s.l. Under climate change scenarios, a decrease in the distribution areas of Lu. peruensis was observed for all representative concentration pathways. However, the centroid of the species' ecological niche showed a northwest direction in all climate change scenarios. The information generated in this study may help health authorities responsible for the supervision of strategies to control leishmaniasis to coordinate, plan and implement appropriate strategies for each area of risk, taking into account the geographic distribution and potential dispersal of Lu. peruensis.

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.

DNA barcoding for identification of sand fly species (Diptera: Psychodidae) from leishmaniasis-endemic areas of Peru

Phlebotomine sand flies are the only proven vectors of leishmaniases, a group of human and animal diseases. Accurate knowledge of sand fly species identification is essential in understanding the epidemiology of leishmaniasis and vector control in endemic areas. Classical identification of sand fly species based on morphological characteristics often remains difficult and requires taxonomic expertise. Here, we generated DNA barcodes of the cytochrome c oxidase subunit 1 (COI) gene using 159 adult specimens morphologically identified to be 19 species of sand flies, belonging to 6 subgenera/species groups circulating in Peru, including the vector species. Neighbor-joining (NJ) analysis based on Kimura 2-Parameter genetic distances formed non-overlapping clusters for all species. The levels of intraspecific genetic divergence ranged from 0 to 5.96%, whereas interspecific genetic divergence among different species ranged from 8.39 to 19.08%. The generated COI barcodes could discriminate between all the sand fly taxa. Besides its success in separating known species, we found that DNA barcoding is useful in revealing population differentiation and cryptic diversity, and thus promises to be a valuable tool for epidemiological studies of leishmaniasis.

Genetic diversity and population structure in the Leishmania guyanensis vector Lutzomyia anduzei (Diptera, Psychodidae) from the Brazilian Amazon

Lutzomyia (Nyssomyia) anduzei has been recognized as a secondary vector of Leishmania guyanensis in the Brazilian Amazon region. Since L. anduzei is anthropophilic, co-exists in areas of high leishmaniasis transmission and has been found infected with L. guyanensis, the understanding of the vector population structure and of the process responsible for it is paramount to the vector management and control efforts. In this study we analyzed 74 and 67 sequences of the COI and Cytb loci, respectively, from mitochondrial DNA, aiming to estimate the intra-population genetic variability and population structure in six L. anduzei samples from the Brazilian Amazon region. For COI, we found 58 haplotypes, low to high (FST=0.0310-0.4128) and significant (P=0.0033) genetic structure, and reduced gene flow among populations. The haplotype network yielded many reticulations that likely resulted from hypervariability in the locus. For Cytb, we observed 27 haplotypes, low to moderate (FST=0.0077-0.1954) and nonsignificant (P>0.05) genetic structure for the majority of comparisons and extensive gene flow among populations, in line with the haplotypes network data. AMOVA analysis indicated that most of the variation occurred within populations (83.41%, 90.94%); nevertheless, there were significant differences (ΦST=0.0906-0.1659; P=0.00098; P=0.00000) among them for both loci. The Mantel test showed that the genetic structure is not associated to an isolation-by-distance (IBD) model in either of both loci. These data suggest that L. anduzei is genetically very diverse. The genetic structure lacking IBD may be due to adaptation to local habitats and the low dispersal capacity of the sandflies, and both could lead to population fragmentation and geographic isolation. These findings have important implications for epidemiology, surveillance and vector control and may be a first step in understanding the evolutionary history of this species.

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.