DNA barcoding of Sri Lankan phlebotomine sand flies using cytochrome c oxidase subunit I reveals the presence of cryptic species

Impact of climate and land use on the temporal variability of sand fly density in Sri Lanka: A 2-year longitudinal study

BACKGROUND

Leishmaniasis has emerged as an escalating public health problem in Sri Lanka, with reported cases increasing nearly three folds over past decade, from 1,367 in 2014 to 3714 cases in 2023. Phlebotominae sand flies are the vectors of leishmaniasis. Their density is known to be influenced by context-specific climatic and land use patterns. Thus, we aimed to investigate how these factors drive sand fly density across Sri Lanka.

METHODOLOGY/PRINCIPAL FINDINGS

We analysed monthly collections of sand flies (n = 38,594) and weather data from ten sentinel sites representing three main geo-climatic zones across Sri Lanka, over 24 months. Site-specific land use data was also recorded. The influence of climate and land use patterns on sand fly density across the sentinel sites were estimated using distributed lag non-linear models and machine learning. We found that climate played a major role on sand fly density compared to land use structure. Increase in rainfall and relative humidity at real time, and ambient temperature and soil temperature with a 2-month lag were associated with a statistically significant increase in sand fly density. The maximum relative risk (RR) observed was 3.76 (95% CI: 1.58-8.96) for rainfall at 120 mm/month, 2.14 (95% CI: 1.04-4.38) for relative humidity at 82% (both at real time). The maximum RR was 2.81 (95% CI: 1.09-7.35) for ambient temperature at 34.5°C, and 11.6 (95% CI, 4.38-30.76) for soil temperature (both at a 2-month lag). The real-time increase in ambient temperature, sunshine hours, and evaporation rate, however, reduced sand fly density homogeneously in all study settings. The high density of chena and coconut plantations, together with low density of dense forests, homesteads, and low human footprint values, positively influenced sand fly density.

CONCLUSIONS/SIGNIFICANCE

The findings improve our understanding of the dynamic influence of environment on sand fly densities and spread of leishmaniasis. This knowledge lays a foundation for forecasting of sand fly densities and designing targeted interventions for mitigating the growing burden of leishmaniasis among the most vulnerable populations, particularly in an era of changing climate.

Impact of climatic factors on temporal variability of sand fly abundance in Sri Lanka: Longitudinal study (2018 to 2020) with two-stage hierarchical analysis

Background

Phlebotomine sand flies serve as vectors for leishmaniasis, a major health concern, but a neglected tropical disease. The risk of vector activity is governed by climatic factors that vary in different geographic zones in the country. Thus, we aimed to quantify the effect of climatic variables on sand fly vector activity in ten sentinel sites across Sri Lanka.

Methods

Mean rainfall, ambient temperature, relative humidity, wind speed, soil temperature, evaporation, sunshine hours, and vector densities were recorded at monthly intervals in each location from March 2018 to February 2020. The association between weather variables and sand fly densities was analysed using a two-staged hierarchical procedure; Distributed Lag Non-Linear (DLNM) modelling framework and the DLNM method implemented in the R package dlnm (version number 2.4.6).

Results

Moderate rainfall values up to 120 mm per month and increasing RH up to 82 at lag of 0 months along with increasing soil temperature and evaporation rate at lag of 2 months were associated with statistically significant increase in the sand fly activity. These associations were heterogeneous across study settings. Whereas increasing ambient and soil temperature, sunshine hours, evaporation rate appeared to reduce the sand fly activity homogeneously at lag of 0 month in all the study settings.

Conclusions

The abundance of sand fly vectors varied in relation to selected climatic variables, either in real-time or with a time lag. This information can be utilized for predicting sand fly densities and for the development of effective strategies to prevent leishmaniasis transmission in specific settings.

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.

A DNA barcode library for katydids, cave crickets, and leaf-rolling crickets (Tettigoniidae, Rhaphidophoridae and Gryllacrididae) from Zhejiang Province, China

Barcode libraries are generally assembled with two main objectives in mind: specimen identification and species discovery/delimitation. In this study, the standard COI barcode region was sequenced from 681 specimens belonging to katydids (Tettigoniidae), cave crickets (Rhaphidophoridae), and leaf-rolling crickets (Gryllacrididae) from Zhejiang Province, China. Of these, four COI-5P sequences were excluded from subsequent analyses because they were likely NUMTs (nuclear mitochondrial pseudogenes). The final dataset consisted of 677 barcode sequences representing 90 putative species-level taxa. Automated cluster delineation using the Barcode of Life Data System (BOLD) revealed 118 BINs (Barcodes Index Numbers). Among these 90 species-level taxa, 68 corresponded with morphospecies, while the remaining 22 were identified based on reverse taxonomy using BIN assignment. Thirteen of these morphospecies were represented by a single barcode (so-called singletons), and each of 19 morphospecies were split into more than one BIN. The consensus delimitation scheme yielded 55 Molecular Operational Taxonomic Units (MOTUs). Only four morphospecies (I _max > DNN) failed to be recovered as monophyletic clades (i.e., Elimaeaterminalis, Phyllomimusklapperichi, Sinochloraszechwanensis and Xizicushowardi), so it is speculated that these may be species complexes. Therefore, the diversity of katydids, cave crickets, and leaf-rolling crickets in Zhejiang Province is probably slightly higher than what current taxonomy would suggest.

Multiple species delimitation approaches with COI barcodes poorly fit each other and morphospecies – An integrative taxonomy case of Sri Lankan Sericini chafers (Coleoptera: Scarabaeidae)

DNA taxonomy including barcoding and metabarcoding is widely used to explore the diversity in biodiversity hotspots. In most of these hotspot areas, chafers are represented by a multitude of species, which are well defined by the complex shape of male genitalia. Here, we explore how well COI barcode data reflect morphological species entities and thus their usability for accelerated species inventorization. We conducted dedicated field surveys in Sri Lanka to collect the species‐rich and highly endemic Sericini chafers (Coleoptera: Scarabaeidae). Congruence among results of a series of protocols for de novo species delimitation and with morphology‐based species identifications was investigated. Different delimitation methods, such as the Poisson tree processes (PTP) model, Statistical Parsimony Analysis (TCS), Automatic Barcode Gap Discovery (ABGD), Assemble Species by Automatic Partitioning (ASAP), and Barcode Index Number (BIN) assignments, resulted in different numbers of molecular operational taxonomic units (MOTUs). All methods showed both over‐splitting and lumping of morphologically identified species. Only 18 of the observed 45 morphospecies perfectly matched MOTUs from all methods. The congruence of delimitation between MOTUs and morphospecies expressed by the match ratio was low, ranging from 0.57 to 0.67. TCS and multirate PTP (mPTP) showed the highest match ratio, while (BIN) assignment resulted in the lowest match ratio and most splitting events. mPTP lumped more species than any other method. Principal coordinate analysis (PCoA) on a match ratio‐based distance matrix revealed incongruent outcomes of multiple DNA delimitation methods, although applied to the same data. Our results confirm that COI barcode data alone are unlikely to correctly delimit all species, in particular, when using only a single delimitation approach. We encourage the integration of various approaches and data, particularly morphology, to validate species boundaries.

Genetic diversity and population structure of Phlebotomus argentipes: Vector of Leishmania donovani in Sri Lanka

Phlebotomus argentipes is the vector of Leishmania donovani which causes the disease leishmaniasis, a neglected tropical disease and a growing health problem in Sri Lanka. A proper understanding of the population genetic structure of sand fly vectors is considered important prior to planning and implementation of a successful vector control program. Thus, the present study was conducted to determine the population genetic structure of sand fly vectors in Sri Lanka. Two mitochondrial genes namely Cytochrome c oxidase subunit 1 (Cox 1) and Cytochrome b (Cytb), and the internal transcribed spacer 2 (ITS2) region from the nuclear ribosomal DNA were used for molecular characterization. Analyses included maximum likelihood method, network analysis and DNA polymorphisms. The outcome revealed unique sequences of all genomic regions studied except the cox 1 gene had a relationship with sand flies isolated previously from Sri Lanka, India and Israel and cytb gene of 4 sand flies that aligned with those isolated earlier from Sri Lanka and 3 from Madagascar. Furthermore, cox 1 gene and ITS 2 region analyses based on F_ST values indicated a possible gene flow between the study sites whereas cytb gene analysis favoured the existence of genetically distinct populations of P. argentipes in each of the study sites. Poor population differentiation of P. argentipes, a possible consequence of a gene flow, is indeed of concern due to the risk imposed by promoting the spread of functionally important phenotypes such as insecticide resistance across the country, making future vector control efforts challenging.

Beyond taxonomy: species complexes in New World phlebotomine sand flies

A species complex (= species group, species series) is an assemblage of species, which are related morphologically and phylogenetically. Recent research has revealed several arthropod vector species that were believed to be a single nominal species actually representing a group of closely related species, which are sometimes morphologically indistinguishable at one or more developmental stages. In some instances, differences in terms of vector competence, capacity, or both have been recorded. It highlights the importance of detecting and studying species complexes to improve our understanding of pathogen transmission patterns, which may be vectored more or less efficiently by different species within the complex. Considering more than 540 species, about one-third of the phlebotomine sand flies in the New World present males and/or females morphologically indistinguishable to one or more species. Remarkably, several of these species may act in transmission of pathogenic agents. In this article, we review recent research on species complexes in phlebotomine sand flies from the Americas. Possible practical implications of recently acquired knowledge and future research needs are also discussed.

Sand fly identification and screening for Leishmania spp. in six provinces of Thailand

Background

Phlebotomine sand flies are vectors of Leishmania spp. At least 27 species of sand flies have been recorded in Thailand. Although human leishmaniasis cases in Thailand are mainly imported, autochthonous leishmaniasis has been increasingly reported in several regions of the country since 1999. Few studies have detected Leishmania infection in wild-caught sand flies, although these studies were carried out only in those areas reporting human leishmaniasis cases. The aim of this study was therefore to identity sand fly species and to investigate Leishmania infection across six provinces of Thailand.

Methods

Species of wild-caught sand flies were initially identified based on morphological characters. However, problems identifying cryptic species complexes necessitated molecular identification using DNA barcoding in parallel with identification based on morphological characters. The wild-caught sand flies were pooled and the DNA isolated prior to the detection of Leishmania infection by a TaqMan real-time PCR assay.

Results

A total of 4498 sand flies (1158 males and 3340 females) were caught by trapping in six provinces in four regions of Thailand. The sand flies were morphologically classified into eight species belonging to three genera (Sergentomyia, Phlebotomus and Idiophlebotomus). Sergentomyia iyengari was found at all collection sites and was the dominant species at most of these, followed in frequency by Sergentomyia barraudi and Phlebotomus stantoni, respectively. DNA barcodes generated from 68 sand flies allowed sorting into 14 distinct species with 25 operational taxonomic units, indicating a higher diversity (by 75%) than that based on morphological identification. Twelve barcoding sequences could not be assigned to any species for which cytochrome c oxidase subunit I sequences are available. All tested sand flies were negative for Leishmania DNA.

Conclusions

Our results confirm the presence of several sand fly species in different provinces of Thailand, highlighting the importance of using DNA barcoding as a tool to study sand fly species diversity. While all female sand flies tested in this study were negative for Leishmania, the circulation of Leishmania spp. in the investigated areas cannot be ruled out.

Graphical abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-021-04856-6.

A Comprehensive Review of Cutaneous Leishmaniasis in Sri Lanka and Identification of Existing Knowledge Gaps

Purpose

Sri Lanka is endemic to cutaneous leishmaniasis and reported as the latest focus of leishmaniasis in the Asian subcontinent. Annually, the number of leishmaniasis cases is increasing; therefore, more efficient diagnostic tools, treatment methods and effective prevention measures are indispensable. For this reason, many studies were conducted regarding leishmaniasis infections in Sri Lanka; however, some areas need more attention. Thus, in this review, we comprehensively discussed the studies on leishmaniasis carried out in Sri Lanka.

Methods

Published articles on leishmaniasis in Sri Lanka were searched on PubMed, Google Scholar and ResearchGate databases. Inclusion criteria for the articles were based on keyword searches including ‘Leishmaniasis in Sri Lanka’, ‘Leishmaniasis vector in Sri Lanka’, ‘Sandfly species in Sri Lanka’, ‘Leishmaniasis epidemiology in Sri Lanka’ which are publicly accessible as of 15th July 2019.

Results

In this study, we evaluated and summarized the leishmaniasis reports in Sri Lanka and mainly focused on clinical presentation of leishmaniasis infection, genetic characteristics of Leishmania donovani Sri Lankan strain, geographical distribution and associated environmental factors, immunological aspects of the infection, vector, reservoir host, risk factors, diagnosis and treatment, and prevention and control. Furthermore, we identified the areas where further research is needed to fill the essential knowledge gaps.

Conclusions

Leishmaniasis has become a critically important parasitic infection in Sri Lanka, whereas the significant clinical form is cutaneous leishmaniasis. Prevalence of the leishmaniasis infections is reported from all the districts of the country. Therefore, more studies are essential to be carried out to fill the existing knowledge gaps emphasized in this review.

Alterations and Interchange of Morphometric Characters in Different Life Cycle Stages with Reference to Genomic Variations of Anopheles subpictus (Diptera; Culicidae) Sibling Species Complex in Sri Lanka

The species complex of the mosquito Anopheles subpictus is designated by the sibling species A⁻D, depending on morphological characters of life cycle stages and variations in polytene chromosomes. However, morphological aberrations in the life cycle stages make the identification of sibling species uncertain and imprecise. The objective of the present study is to determine the suitability of morphological variations of sibling species and their genomic variations to identify the sibling species status of an An. subpictus population in Sri Lanka. Life cycle stages of larvae, pupal exuviae, and adults were examined for previously reported distinctive morphological features. Five nuclear and mitochondrial genome regions, including the Internal transcribed spacer 2 (ITS2) region, D3 region, white gene, cytochrome c oxidase I (COI), and Cytochrome b (Cyt-b), were sequenced and analyzed for variations. The eggs changed their distinct sibling morphological characters during metamorphosis (89.33%). The larvae, pupal exuviae, and adult stages showed deviation from their sibling characters by 26.10%, 19.71%, and 15.87%, respectively. However, all the species from the analysis shared two distinct sequence types for all regions, regardless of the morphological variations. In conclusion, the An. subpictus sibling species complex in Sri Lanka is not identifiable using morphological characters due to variations, and the genomic variations are independent from the morphological variations.

Species identification of horse flies (Diptera: Tabanidae) in Thailand using DNA barcoding

Horse flies (Diptera: Tabanidae) are of medical and veterinary importance because they are known to transmit pathogens. Approximately 80 species of horse flies have been reported in Thailand. Monitoring the distribution of horse fly species is important to control the spread of diseases transmitted by them. Currently, the species identification of horse flies is based on their morphology; this requires considerable skills and taxonomic expertise, and it may be difficult to identify morphologically similar species. DNA-based identification methods are increasingly being developed for rapid and accurate identification of various insect species. In this study, we used mitochondrial cytochrome oxidase subunit I (COI) for species identification of horse flies in Thailand. A 658 bp fragment of COI was amplified from 145 adult horse flies belonging to 48 morphologically distinct species and sequenced. Sequence analysis revealed an intraspecific divergence of 0.0%-4.4% and an interspecific divergence of 0.0%-16.2%. Our results showed that COI barcodes were effective in discriminating a majority of horse flies in Thailand on the basis of the barcoding gap and phylogenetic analyses. However, COI barcodes were unable to distinguish among members of the Tabanus striatus complex and some species within the T. ceylonicus group.

DNA Barcoding Reveals Hidden Diversity of Sand Flies (Diptera: Psychodidae) at Fine and Broad Spatial Scales in Brazilian Endemic Regions for Leishmaniasis

Sand fly (Diptera: Psychodidae) taxonomy is complex and time-consuming, which hampers epidemiological efforts directed toward controlling leishmaniasis in endemic regions such as northeastern Brazil. Here, we used a fragment of the mitochondrial cytochrome c oxidase I (COI) gene to identify sand fly species in Maranhão State (northeastern Brazil) and to assess cryptic diversity occurring at different spatial scales. For this, we obtained 148 COI sequences of 15 sand fly species (10 genera) from Maranhão (fine spatial scale), and joined them to COI sequences from other Brazilian localities (distant about 2,000 km from Maranhão, broad spatial scale) available in GenBank. We revealed cases of cryptic diversity in sand flies both at fine (Lutzomyia longipalpis (Lutz and Neiva) and Evandromyia termitophila (Martins, Falcão and Silva)) and broad spatial scales (Migonemyia migonei (França), Pressatia choti (Floch and Abonnenc), Psychodopygus davisi (Root), Sciopemyia sordellii (Shannon and Del Ponte), and Bichromomyia flaviscutellata (Mangabeira)). We argue that in the case of Bi. flaviscutellata, the cryptic diversity is associated with a putative new species. Cases in which DNA taxonomy was not as effective as morphological identification possibly involved recent speciation and/or introgressive hybridization, highlighting the need for integrative approaches to identify some sand fly species. Finally, we provide the first barcode sequences for four species (Brumptomyia avellari (Costa Lima), Evandromyia infraspinosa (Mangabeira), Evandromyia evandroi (Costa Lima and Antunes), and Psychodopygus complexus (Mangabeira)), which will be useful for further molecular identification of neotropical species.

Identification of wild-caught phlebotomine sand flies from Crete and Cyprus using DNA barcoding

Background

Phlebotomine sand flies (Diptera: Psychodidae) are vectors of Leishmania spp., protozoan parasites responsible for a group of neglected diseases called leishmaniases. Two sand fly genera, Phlebotomus and Sergentomyia, contain species that are present in the Mediterranean islands of Crete and Cyprus where the visceral (VL), cutaneous (CL) and canine (CanLei) leishmaniases are a public health concern. The risk of transmission of different Leishmania species can be studied in an area by monitoring their vectors. Sand fly species are traditionally identified using morphological characteristics but minute differences between individuals or populations could be overlooked leading to wrong epidemiological predictions. Molecular identification of these important vectors has become, therefore, an essential tool for research tasks concerning their geographical distribution which directly relates to leishmaniasis control efforts. DNA barcoding is a widely used molecular identification method for cataloguing animal species by sequencing a fragment of the mitochondrial gene encoding cytochrome oxidase I.

Results

DNA barcoding was used to identify individuals of five sand fly species (Phlebotomus papatasi, P. similis, P. killicki, Sergentomyia minuta, S. dentata) circulating in the islands of Crete and Cyprus during the years 2011–2014. Phlebotomus papatasi is a known vector of zoonotic CL in the Middle East and it is found in both islands. Phlebotomus similis is the suspected vector of Leishmania tropica in Greece causing anthroponotic CL. Phlebotomus killicki was collected in Cyprus for the first time. Sergentomyia minuta, found to present intraspecific diversity, is discussed for its potential as a Leishmania vector. Molecular identification was consistent with the morphological identification. It successfully identified males and females, which is difficult when using only morphological characters. A phylogenetic tree was constructed based on the barcodes acquired, representing their genetic relationships along with other species from the area studied. All individuals identified were clustered according to their species and subgenus.

Conclusions

Molecular identification of sand flies via DNA barcoding can accurately identify these medically important insects assisting traditional morphological tools, thus helping to assess their implication in Leishmania transmission.