Single-locus DNA barcoding and species delimitation of the sandfly subgenus Evandromyia (Aldamyia)

Hidden diversity in anthropophilic sand flies of the Monticola Series (Diptera, Psychodidae)

The Monticola series comprises two anthropophilic and widely distributed species in Brazil: Pintomyia (Pifanomyia) monticola (Costa Lima, 1932) and Pintomyia (Pifanomyia) misionensis (Castro, 1959). They mainly occur in the Atlantic Rainforest, and it is known that Pi. monticola comprises at least two well-structured genetic lineages regarding a fragment of the cytochrome c oxidase subunit I (COI) gene. Here, we aim to elucidate the taxonomic status of this group using integrative taxonomy tools. Collections were performed in nine localities of four Brazilian states, and COI fragments were sequenced and merged with publicly available data. Several single-locus species delimitation algorithms, genetic distance metrics, phylogenetic trees, and haplotype networks were used to uncover cryptic diversity and population structure within Pi. monticola and Pi. misionensis. The resulting genetic clusters were then tested for morphological differences through linear and geometric morphometry of several characters. We analyzed 152 COI sequences, comprising 48 haplotypes. The maximum intraspecific p distances were 8.21% (mean 4.17%) and 9.12% (mean 4.4%) for Pi. monticola and Pi. misionensis, respectively, while interspecific ones ranged from 10.94 to 14.09% (mean 12.33%). Phylogenetic gene trees showed well-supported clades for both species, with clear structuring patterns within them. Species-delimitation algorithms split our dataset into at least three putative species for each taxon. Moreover, population structure analysis showed a strong correlation between Atlantic Forest areas of endemism as sources of molecular variation in Pi. monticola. Morphometric analyses were significant for wing shape variation and some linear measurements (mainly of the head) when comparing specimens of different genetic clusters for both taxa. These results indicate strong genetic structuring of Monticola series species, confirmed by morphometry, indicating two possible cryptic species complexes.

Molecular and morphometric study of Brazilian populations of Psychodopygus davisi

In this study, we analysed the molecular and morphometric differences of several populations of the putative sand fly vector Psychodopygus davisi (Root, 1934) (Diptera, Psychodidae, Phlebotominae) in Brazil. We amplified the 658 base pair fragments of the DNA barcoding region-cytochrome c oxidase subunit 1 (COI) gene-for 57 specimens of P. davisi and three specimens of Psychodopygus claustrei (Abonnenc, Léger & Fauran, 1979). We merged our data with public sequences of the same species available from GenBank. Then, the combined dataset-87 sequences and 20 localities-was analysed using population structure analysis and different species delimitation approaches. Geometric morphometry of wings was performed for 155 specimens of P. davisi populations from the North, Midwest and Southeast Brazilian regions, analysing the differences in centroid sizes and canonical variates. Molecular analysis indicated high intraspecific genetic distance values for P. davisi (maximum p distance = 5.52%). All algorithms identified P. davisi and P. claustrei as distinct molecular taxonomic units, despite the low interspecific distance (p distance to the nearest neighbour = 4.79%). P. davisi sequences were split into four genetic clusters by population structure analysis and at least five genetic lineages using intermediate scenarios of the species delimitation algorithms. The species validation analysis of BPP strongly supported the five-species model in our dataset. We found high genetic diversity in this taxon, which is in agreement with its wide geographic distribution in Brazil. Furthermore, the wing analysis showed that specimens from the Southeast Region of Brazil are different from those in the North and the Midwest. The evolutionary patterns of P. davisi populations in Brazil suggest the presence of candidate species, which need to be validated in future studies using a more comprehensive approach with both genomic data and morphological characters.

New sand fly (Diptera, Psychodidae) records and COI DNA barcodes in the state of Maranhão, Eastern Amazon, Brazil

The sand fly fauna and the usefulness of the DNA barcoding fragment of the cytochrome c oxidase subunit I (COI) gene were accessed in a forest fragment in the municipality of Governador Newton Bello, state of Maranhão, Brazil. We performed entomological collections in three independent campaigns in May and October 2021, and January 2023. Sand flies were morphologically-identified and then DNA barcoded. Sequences were deposited and analyzed in the BOLD System Database, and various species delimitation algorithms, to assess whether DNA sequences merge into taxonomic units in accordance with nominal species. In total, 1,524 sand flies were collected, comprising 32 nominal species. Nyssomyia antunesi was the most abundant species (31.5 %), followed by Psychodopygus davisi (27 %). We reported for the first time in the state of Maranhão, the presence of Lutzomyia evangelistai, Lutzomyia sherlocki, Pressatia equatorialis, and Psathyromyia barrettoi. We amplified and analyzed 67 COI barcodes of 23 species, which were merged with conspecific sequences extracted from GenBank. The maximum intraspecific p distances ranged from 0.0 % to 14.74 %, while the distances to the nearest neighbor varied from 1.67 % to 13.64 %. The phylogenetic gene tree and species delimitation tools clustered sequences into well-supported clades/clusters for each nominal species, except for Pressatia choti/Pr. equatorialis, which have the lowest interspecific genetic distance (1.67 %). We sequenced for the first time COI barcodes of Brumptomyia brumpti, Evandromyia monstruosa, Micropygomyia rorotaensis, Micropygomyia pilosa, Pintomyia christenseni, Pintomyia pacae, Pr. equatorialis, Pa. barrettoi, and Psathyromyia hermanlenti, which will be useful for further molecular identification and classification proposals of Neotropical species. This study updated the current list of the sand fly fauna for the state of Maranhão to 97, and demonstrated that COI barcodes are useful for specific identification.

Trichophoromyia auraensis: evidence for cryptic species and first record in the state of Maranhão, Brazil

Trichophoromyia auraensis (Mangabeira, 1942) (Diptera, Psychodidae, Phlebotominae) has a wide geographic distribution in the western region of the Amazon biome, where it is a putative Leishmania vector. Here, we reported for the first time a population of this species in the Brazilian state of Maranhão, in the eastern Amazon, from which we DNA-barcoded and compared with previously processed specimens from Acre State, in the western Amazon. For this, we analyzed the DNA barcoding fragment (658 bp) of the mitochondrial cytochrome c oxidase subunit I (COI) gene and the nuclear internal transcribed spacer 2 (ITS2) of Trichophoromyia species using phylogenetic gene trees, and species delimitation algorithms. The analyses of COI barcodes showed high values of genetic distance (mean K2P = 5.17) and well-supported clades/MOTUs for the eastern and western populations of T. auraensis, which may indicate a possible complex of cryptic species. The western population of this taxon merged with the close-related sand fly Trichophoromyia velezbernali Posada-López, Galvis and Galati, 2018 from Colombia, which may be associated with the recent speciation history and introgression between these populations. These evidences should be evaluated with a more comprehensive sampling in terms of analyzed populations and molecular markers.

DNA barcodes provide insights into the diversity and biogeography of the non-biting midge Polypedilum (Diptera, Chironomidae) in South America

South America, particularly within its tropical belt, is renowned for its unparalleled high levels of species richness, surpassing other major biomes. Certain neotropical areas harbor fragmented knowledge of insect diversity and face imminent threats from biodiversity loss and climate change. Hence, there is an urgent need for rapid estimation methods to complement slower traditional taxonomic approaches. A variety of algorithms for delimiting species through single-locus DNA barcodes have been developed and applied for rapid species diversity estimates across diverse taxa. However, tree-based and distance-based methods may yield different group assignments, leading to potential overestimation or underestimation of putative species. Here, we investigate the performance of different DNA-based species delimitation approaches to rapidly estimate the diversity of Polypedilum (Chironomidae, Diptera) in South America. Additionally, we test the hypothesis that significant differences exist in the community structure of Polypedilum fauna between South America and its neighboring regions, particularly the Nearctic. Our analysis encompasses a dataset of 1492 specimens from 598 locations worldwide, with a specific focus on South America. Within this region, we analyzed a subset of 247 specimens reported from 37 locations. Using various methods including the Barcode Index Number (BIN), Bayesian Poisson tree processes (bPTP), multi-rate Poisson tree processes (mPTP), single-rate Poisson tree processes (sPTP), and generalized mixed Yule coalescent (sGMYC), we identify molecular operational taxonomic units (MOTUs) ranging from 267 to 520. Our results indicate that the sGMYC method is the most suitable for estimating putative species in our dataset, resulting in the identification of 75 species in the Neotropical region, particularly in South America. Notably, this region exhibited higher species richness in comparison to the Palearctic and Oriental realms. Additionally, our findings suggest potential differences in species composition of Polypedilum fauna between the Neotropical and the adjacent Nearctic realms, highlighting high levels of endemism and species richness in the first. These results support our hypothesis that there are substantial differences exist in species composition between the Polypedilum fauna in South America and the neighboring regions.

Improving the COI DNA barcoding library for Neotropical phlebotomine sand flies (Diptera: Psychodidae)

Sand fly species are traditionally identified using morphological traits, though this method is hampered by the presence of cryptic species. DNA barcoding is a widely used tool in the case of insects of medical importance, where it is necessary to know quickly which species are present in a transmission area. Here, we assess the usefulness of mitochondrial cytochrome c oxidase subunit I (COI) DNA barcoding as a practical tool for species identification, correct assignment of isomorphic females, and to evaluate the detection of cryptic diversity that occurs in the same species. A fragment of the COI gene was used to generate 156 new barcode sequences for sand flies from different countries of the Neotropical region, mainly Colombia, which had been identified morphologically as 43 species. The sequencing of the COI gene allowed the detection of cryptic diversity within species and correctly associated isomorphic females with males identified by morphology. The maximum intraspecific genetic distances ranged from 0 to 8.32% and 0 to 8.92% using uncorrected p distances and the Kimura 2-parameter (K2P) model, respectively. The minimum interspecific distance (nearest neighbor) for each species ranged from 1.5 to 14.14% and 1.51 to 15.7% using p and K2P distances, respectively. Three species had more than 3% maximum intraspecific distance: Psychodopygus panamensis, Micropygomyia cayennensis cayennensis, and Pintomyia evansi. They also were split into at least two molecular operational taxonomic units (MOTUs) each, using different species delimitation algorithms. Regarding interspecific genetic distances, the species of the genera Nyssomyia and Trichophoromyia generated values lower than 3% (except Nyssomyia ylephiletor and Ny. trapidoi). However, the maximum intraspecific distances did not exceed these values, indicating the presence of a barcode gap despite their proximity. Also, nine sand fly species were DNA barcoded for the first time: Evandromyia georgii, Lutzomyia sherlocki, Ny. ylephiletor, Ny. yuilli pajoti, Psathyromyia punctigeniculata, Sciopemyia preclara, Trichopygomyia triramula, Trichophoromyia howardi, and Th. velezbernali. The COI DNA barcode analysis enabled the correct delimitation of several Neotropical sand fly species from South and Central America and raised questions about the presence of cryptic species for some taxa, which should be further assessed.

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.

Comparing the Efficiency of Single-Locus Species Delimitation Methods within Trochoidea (Gastropoda: Vetigastropoda)

In the context of diminishing global biodiversity, the validity and practicality of species delimitation methods for the identification of many neglected and undescribed biodiverse species have been paid increasing attention. DNA sequence-based species delimitation methods are mainly classified into two categories, namely, distance-based and tree-based methods, and have been widely adopted in many studies. In the present study, we performed three distance-based (ad hoc threshold, ABGD, and ASAP) and four tree-based (sGMYC, mGMYC, PTP, and mPTP) analyses based on Trochoidea COI data and analyzed the discordance between them. Moreover, we also observed the performance of these methods at different taxonomic ranks (the genus, subfamily, and family ranks). The results suggested that the distance-based approach is generally superior to the tree-based approach, with the ASAP method being the most efficient. In terms of phylogenetic methods, the single threshold version performed better than the multiple threshold version of GMYC, and PTP showed higher efficiency than mPTP in delimiting species. Additionally, GMYC was found to be significantly influenced by taxonomic rank, showing poorer efficiency in datasets at the genus level than at higher levels. Finally, our results highlighted that cryptic diversity within Trochoidea (Mollusca: Vetigastropoda) might be underestimated, which provides quantitative evidence for excavating the cryptic lineages of these species.

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.