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

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.

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.

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.

Genetic Structure of Phlebotomus orientalis (Diptera: Psychodidae) in Leishmaniasis Endemic Foci of Sudan

BACKGROUND AND OBJECTIVE

Visceral leishmaniasis (VL) remains a major concern in many parts of Sudan. The disease is transmitted by Phlebotomus orientalis. The objective of this study was to determine genetic structure of Phlebotomus orientalis population from 5 geographical regions in Sudan.

MATERIALS AND METHODS

A total of 194 individual sand flies were collected from 5 geographic regions in Sudan. The field collected sand flies were analyzed by Random Amplified Polymorphic DNA (RAPD) using 30 primers.

RESULTS

Eight hundred and 65 bands from 4 RAPD primers were analyzed for genetic variation. A higher level of intrapopulational variability was detected in populations of P. orientalis from eastern Sudan compared to those populations from central and northern Sudan. Diagnostic bands were detected in populations of P. orientalis central Sudan. Hieratical clustering analysis showed clear clustering into 2 main populations with 1 population subdivided into 4 subpopulations. However, these populations did not show any correlation with their geographical origins. Furthermore, the low genetic differentiation among subpopulations was supported by fixation index (FST) estimated by analysis of molecular variance (AMOVA).

CONCLUSION

It is concluded that the populations of P. orientalis from the selected areas in Sudan have a low genetic differentiation. However, assessment of genetic structure of P. orientalis populations is important for understanding the patterns of transmission of VL in different endemic areas.

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.

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.

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.

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.

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.