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

Mitochondrial COI and Cytb gene as valid molecular identification marker of sandfly species (Diptera: Psychodidae) in China

The accurate identification of sandfly species is crucial because some species transmit medically significant diseases, including leishmaniasis, bartonellosis and sandfly fever. However, due to the high similarity of the external morphology in sandfly species, identification can only be performed using internal morphological characteristics after dissection, which is time consuming and requires highly experienced staff. Thus, the introduction of suitable molecular markers may solve these identification problems. This study screened suitable DNA barcodes to identify common sandfly species in China. The phlebotomine sandflies were collected from Sichuan, Henan and Hainan Provinces from 2014 to 2016. The species were identified by the morphological characteristics of the pharyngeal armature and spermatheca. The genomic DNA of sandfly was extracted individually, and mitochondrial DNA (mtDNA) cytochrome C oxidase subunit I (COI) and cytochrome B (Cytb) as well as the 18S subunit of ribosomal DNA (rDNA) were amplified using polymerase chain reaction (PCR). Additionally, intraspecific and interspecific differences (p-distance) were calculated to evaluate the feasibility of the three gene fragments as a DNA barcode. The phylogeny trees of all sandfly species in this study were constructed using neighbor joining (NJ) method. Six species were identified by the morphological features, belonging to Phlebotomus and Sergentomyia, as Ph. chinensis s. l., Ph. stantoni, Se. bailyi, Se. iyengari, Se. squamirostris, and Se. squamipleuris. Analysis based on three gene fragments revealed some degree of intraspecific polymorphism among these sandfly species in China. The largest intraspecific variation occurred in Ph. chinensis s. l. (mtDNA COI, p-distance = 0.042; mtDNA Cytb, p-distance = 0.071), but the 18S rDNA fragment showed a small variation (p-distance = 0.005). The ranges of interspecific p-distances for mtDNA COI and mtDNA Cytb were 0.138 - 0.231 and 0.128 - 0.274, respectively. However, the interspecific p-distances of 18S rDNA are relatively low ranging from 0.003 to 0.055. Both mitochondrial COI and Cytb gene fragments are valid molecular identification markers in theses sandfly species. The topological structure of phylogeny trees based on mtDNA COI, mtDNA Cytb and 18S rDNA genes were all consistent with morphological classification. And we also found there were significant intraspecies differences within Ph. chinensis s. l. (0.006-0.071) and Se. bailyi (0.002-0.032) based on mtDNA Cytb gene fragment. Sequence alignment data suggested that Ph. chinensis s. l. from Sichuan should be Ph. sichuanensis, and the sandfly specimen collected from Henan was Ph. chinensis s. s.. There could be cryptic species in Se. bailyi from China.

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.

Comment on Xu et al. Isolation and Identification of a Novel Phlebovirus, Hedi Virus, from Sandflies Collected in China. Viruses 2021, 13, 772

The article from Xu et al. [...].

Macrogeographic genetic structure of Lutzomyia longipalpis complex populations using Next Generation Sequencing

Lutzomyia longipalpis is the main vector of Leishmania infantum, the causative agent of visceral leishmaniasis in the Neotropical realm. Its taxonomic status has been widely discussed once it encompasses a complex of species. The knowledge about the genetic structure of insect vector populations helps the elucidation of components and interactions of the disease ecoepidemiology. Thus, the objective of this study was to genotypically analyze populations of the Lu. longipalpis complex from a macrogeographic perspective using Next Generation Sequencing. Polymorphism analysis of three molecular markers was used to access the levels of population genetic structure among nine different populations of sand flies. Illumina Amplicon Sequencing Protocol^® was used to identify possible polymorphic sites. The library was sequenced on paired-end Illumina MiSeq platform. Significant macrogeographical population differentiation was observed among Lu. longipalpis populations via PCA and DAPC analyses. Our results revealed that populations of Lu. longipalpis from the nine municipalities were grouped into three clusters. In addition, it was observed that the levels of Lu. longipalpis population structure could be associated with distance isolation. This new sequencing method allowed us to study different molecular markers after a single sequencing run, and to evaluate population and inter-species differences on a macrogeographic scale.

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.

Development of polymorphic EST microsatellite markers for the sand fly, Phlebotomus papatasi (Diptera: Psychodidae)

Background

Phlebotomus papatasi is a widely distributed sand fly species in different tropical and sub-tropical regions including the Middle East and North Africa. It is considered an important vector that transmits Leishmania major parasites, the causative agents of cutaneous leishmaniasis. The development of microsatellite markers for this sand fly vector is of high interest to understand its population structure and to monitor its geographic dispersal.

Results

Fourteen polymorphic microsatellite markers were developed with simple di-, tri- and tetra-nucleotide repeats. The F statistics calculated for the 14 markers revealed high genetic diversity; expected heterozygosity (He) ranged from 0.407 to 0.767, while observed heterozygosity (Ho) was lower and ranged from 0.083 to 0.514. The number of alleles sampled fall in the range of 9–29. Three out of 14 markers deviated from Hardy-Weinberg expectations, no significant linkage disequilibrium was detected and high values for inbreeding coefficient (F_IS) were likely due to inbreeding.

Conclusions

The development of these functional microsatellites enable a high resolution of P. papatasi populations. It opens a path for researchers to perform multi locus-based population genetic structure analyses, and comparative mapping, a part of the efforts to uncover the population dynamics of this vector, which is an important global strategy for understanding the epidemiology and control of leishmaniasis.

Electronic supplementary material

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

Diversity of bacteriome associated with Phlebotomus chinensis (Diptera: Psychodidae) sand flies in two wild populations from China

Sand fly Phlebotomus chinensis is a primary vector of transmission of visceral leishmaniasis in China. The sand flies have adapted to various ecological niches in distinct ecosystems. Characterization of the microbial structure and function will greatly facilitate the understanding of the sand fly ecology, which would provide critical information for developing intervention strategy for sand fly control. In this study we compared the bacterial composition between two populations of Ph. chinensis from Henan and Sichuan, China. The phylotypes were taxonomically assigned to 29 genera of 19 families in 9 classes of 5 phyla. The core bacteria include Pseudomonas and enterobacteria, both are shared in the sand flies in the two regions. Interestingly, the endosymbionts Wolbachia and Rickettsia were detected only in Henan, while the Rickettsiella and Diplorickettsia only in Sichuan. The intracellular bacteria Rickettsia, Rickettsiella and Diplorickettsia were reported for the first time in sand flies. The influence of sex and feeding status on the microbial structure was also detected in the two populations. The findings suggest that the ecological diversity of sand fly in Sichuan and Henan may contribute to shaping the structure of associated microbiota. The structural classification paves the way to function characterization of the sand fly associated microbiome.

A Historical Overview of the Classification, Evolution, and Dispersion of Leishmania Parasites and Sandflies

BACKGROUND

The aim of this study is to describe the major evolutionary historical events among Leishmania, sandflies, and the associated animal reservoirs in detail, in accordance with the geographical evolution of the Earth, which has not been previously discussed on a large scale.

METHODOLOGY AND PRINCIPAL FINDINGS

Leishmania and sandfly classification has always been a controversial matter, and the increasing number of species currently described further complicates this issue. Despite several hypotheses on the origin, evolution, and distribution of Leishmania and sandflies in the Old and New World, no consistent agreement exists regarding dissemination of the actors that play roles in leishmaniasis. For this purpose, we present here three centuries of research on sandflies and Leishmania descriptions, as well as a complete description of Leishmania and sandfly fossils and the emergence date of each Leishmania and sandfly group during different geographical periods, from 550 million years ago until now. We discuss critically the different approaches that were used for Leishmana and sandfly classification and their synonymies, proposing an updated classification for each species of Leishmania and sandfly. We update information on the current distribution and dispersion of different species of Leishmania (53), sandflies (more than 800 at genus or subgenus level), and animal reservoirs in each of the following geographical ecozones: Palearctic, Nearctic, Neotropic, Afrotropical, Oriental, Malagasy, and Australian. We propose an updated list of the potential and proven sandfly vectors for each Leishmania species in the Old and New World. Finally, we address a classical question about digenetic Leishmania evolution: which was the first host, a vertebrate or an invertebrate?

CONCLUSIONS AND SIGNIFICANCE

We propose an updated view of events that have played important roles in the geographical dispersion of sandflies, in relation to both the Leishmania species they transmit and the animal reservoirs of the parasites.

Visceral Leishmaniasis in China: an Endemic Disease under Control

Visceral leishmaniasis (VL) caused by Leishmania spp. is an important vector-borne and largely zoonotic disease. In China, three epidemiological types of VL have been described: anthroponotic VL (AVL), mountain-type zoonotic VL (MT-ZVL), and desert-type ZVL (DT-ZVL). These are transmitted by four different sand fly species: Phlebotomus chinensis, P. longiductus, P. wui, and P. alexandri. In 1951, a detailed survey of VL showed that it was rampant in the vast rural areas west, northwest, and north of the Yangtze River. Control programs were designed and implemented stringently by the government at all administrative levels, resulting in elimination of the disease from most areas of endemicity, except the western and northwestern regions. The control programs consisted of (i) diagnosis and chemotherapy of patients, (ii) identification, isolation, and disposal of infected dogs, and (iii) residual insecticide indoor spraying for vector control. The success of the control programs is attributable to massive and effective mobilization of the general public and health workers to the cause. Nationally, the annual incidence is now very low, i.e., only 0.03/100,000 according to the available 2011 official record. The overwhelming majority of cases are reported from sites of endemicity in the western and northwestern regions. Here, we describe in some depth and breadth the current status of epidemiology, diagnosis, treatment, and prevention of the disease, with particular reference to the control programs. Pertinent information has been assembled from scattered literature of the past decades in different languages that are not readily accessible to the scientific community. The information provided constitutes an integral part of our knowledge on leishmaniasis in the global context and will be of special value to those interested in control programs.

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.