Phylogeography of the Lutzomyia gomezi (Diptera: Phlebotominae) on the Panama Isthmus

Insights into the Genetic Diversity of Leishmania (Viannia) panamensis in Panama, Inferred via Multilocus Sequence Typing (MLST)

Leishmaniasis is a disease caused by parasites of the genus Leishmania and transmitted by sand fly vectors. Tegumentary leishmaniasis is the most prevalent clinical outcome in Latin America, afflicting people from 18 countries. In Panama, the annual incidence rate of leishmaniasis is as high as 3000 cases, representing a major public health problem. In endemic regions, L. panamensis is responsible for almost eighty percent of human cases that present different clinical outcomes. These differences in disease outcomes could be the result of the local interplay between L. panamensis variants and human hosts with different genetic backgrounds. The genetic diversity of L. panamensis in Panama has only been partially explored, and the variability reported for this species is based on few studies restricted to small populations and/or with poor resolutive markers at low taxonomic levels. Accordingly, in this study, we explored the genetic diversity of sixty-nine L. panamensis isolates from different endemic regions of Panama, using an MLST approach based on four housekeeping genes (Aconitase, ALAT, GPI and HSP70). Two to seven haplotypes per locus were identified, and regional differences in the genetic diversity of L. panamensis were observed. A genotype analysis evidenced the circulation of thirteen L. panamensis genotypes, a fact that might have important implications for the local control of the disease.

Similarity but not equivalence: Ecological niche comparison between sandflies from the Pleistocene and future scenarios in Central and South America

Sandfly species (Diptera: Psychodidae) are suspected or proven vectors of Leishmania spp. in the American region. Understanding niche conservatism (NC) in insect vectors allows an understanding of constraints on adaptive responses, and thus implications for disease ecology. Therefore, in this study, the authors evaluated NC in three vector species of leishmaniasis (Lutzomyia gomezi, Psathyromyia shannoni and Pintomyia ovallesi) in Central and South America. For this, the authors performed niche identity and similarity testing through paired comparisons in ENMTools and niche overlap in Niche Analyst. The authors found that species niches were more similar to each other than if the points had been randomly extracted, and they also found extensive similarity between Pa. shannoni and Lu. gomezi niches and in Pa. shannoni niches over different timescales. The authors suggest Pa. shannoni as a priority species due to fundamental niche similarity with phylogenetically related species and also its extensive evolutionary history and ecological plasticity that could affect the emergence and resurgence of leishmaniasis in areas endemic by this vector.

FLEBOTOMÍNEOS (PSYCHODIDAE: PHLEBOTOMINAE) EN LA ZONA PERIURBANA DE MONTERÍA (CÓRDOBA- COLOMBIA)

El objetivo de este trabajo fue determinar la diversidad de flebotomíneos (Diptera: Psychodidae) en la zona peri-urbana de Montería, (Córdoba-Colombia). Los especímenes se colectaron en cuatro localidades alrededor de Montería con trampas CDC y aspiradores bucales durante febrero y junio de 2018. Los individuos se fragmentaron y aclararon con lactofenol para su posterior identificación mediante claves taxonómicas especializadas. Se estimó la diversidad con el índice de Shannon-Wiener, las localidades se compararon usando la distancia de Bray-Curtis y se cuantificó la riqueza con los índices de Chao 2 y Jacknife 2. Se identificaron seis especies: Micropygomyia cayennensis, Lutzomyia gomezi, Pintomyia rangeliana, Micropygomyia micropyga, Psathyromyia shannoni y Pintomyia evansi. Los sitios de muestreo fueron similares en cuanto a diversidad. Se concluye que la zona peri-urbana de Montería está habitada por especies generalistas de dípteros de la familia Psychodidae que constituyen un riesgo epidemiológico. La riqueza en el área se ve afectada por las características propias de las especies encontradas.

Population genetics and phylogeography of Tabanus bromius (Diptera: Tabanidae)

Background

Tabanus bromius (Diptera: Tabanidae) is one of the most notable Tabanidae species of veterinary and medical importance distributed throughout the Palearctic region. In this study, we investigate the genetic diversity and the phylogeographic structure of T. bromius sampled from Turkey, Croatia, and Iran.

Methods

For this purpose, a 686-base-pair (bp) fragment of mitochondrial DNA cytochrome oxidase I gene (COI) and 1339 bp of the nuclear DNA internal transcribed spacer (ITS) were sequenced from 247 individuals representing 15 populations.

Results

The sequences generated 169 COI haplotypes and 90 ITS alleles. A higher haplotype/allele diversity (h = 0.9909 for the COI gene and Ad = 0.8193 for the ITS region) compared to a low nucleotide diversity (π = 0.020605 for COI gene and π = 0.013667 for the ITS region), present for a high number of singleton and private haplotypes/alleles imply population expansion in the past. The results of phylogenetic analysis led to the uncovering of geographically significant groupings of lineages with regard to the entrance of the species into Anatolia and the location of major geographic barriers. According to current data, the species appears to have entered Turkey from Caucasia and Iran. A molecular clock applied to the COI data suggests that T. bromius diverged from the outgroup species nearly 8.83 million years ago, around the end of the Miocene era.

Conclusions

The results of this study indicate remarkable genetic diversity across the studied range of the species. High haplotype/allele versus low nucleotide diversity and demographic analyses implied that the T. bromius populations have undergone a series of expansions and retreats in the past. Our current findings suggest that T. bromius split from outgroups around the Late Miocene. Subsequent diversification events during the climatic and environmental fluctuation times of the Late Pliocene and Early Pleistocene periods also significantly influenced the species, resulting in the formation of some major genetic lineages. The phylogenetic analyses indicate that T. bromius most likely entered Turkey from the Caucasus region and Iran.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-021-04970-5.

Gut Microbiota Dynamics in Natural Populations of Pintomyia evansi under Experimental Infection with Leishmania infantum

Pintomyia evansi is recognized by its vectorial competence in the transmission of parasites that cause fatal visceral leishmaniasis in rural and urban environments of the Caribbean coast of Colombia. The effect on and the variation of the gut microbiota in female P. evansi infected with Leishmania infantum were evaluated under experimental conditions using 16S rRNA Illumina MiSeq sequencing. In the coinfection assay with L. infantum, 96.8% of the midgut microbial population was composed mainly of Proteobacteria (71.0%), followed by Cyanobacteria (20.4%), Actinobacteria (2.7%), and Firmicutes (2.7%). In insect controls (uninfected with L. infantum) that were treated or not with antibiotics, Ralstonia was reported to have high relative abundance (55.1–64.8%), in contrast to guts with a high load of infection from L. infantum (23.4–35.9%). ASVs that moderately increased in guts infected with Leishmania were Bacillus and Aeromonas. Kruskal–Wallis nonparametric variance statistical inference showed statistically significant intergroup differences in the guts of P. evansi infected and uninfected with L. infantum (p < 0.05), suggesting that some individuals of the microbiota could induce or restrict Leishmania infection. This assay also showed a negative effect of the antibiotic treatment and L. infantum infection on the gut microbiota diversity. Endosymbionts, such as Microsporidia infections (<2%), were more often associated with guts without Leishmania infection, whereas Arsenophonus was only found in guts with a high load of Leishmania infection and treated with antibiotics. Finally, this is the first report that showed the potential role of intestinal microbiota in natural populations of P. evansi in susceptibility to L. infantum infection.

Leishmania spp. Infection Rate and Feeding Patterns of Sand Flies (Diptera: Psychodidae) from a Hyperendemic Cutaneous Leishmaniasis Community in Panamá

American cutaneous leishmaniasis (ACL) is a common and important vector-borne parasitic zoonosis in Panamá. Here, we study Leishmania spp. infection rates and blood-feeding patterns among common sand flies in Trinidad de Las Minas, a rural community with hyperendemic ACL transmission, and where a deltamethrin fogging trial was performed. Sand flies were collected from April 2010 to June 2011 with light traps installed inside and in the peridomicile of 24 houses. We restricted our analysis to the most abundant species at the study site: Lutzomyia trapidoi, Lutzomyia gomezi, Lutzomyia panamensis, Lutzomyia triramula, and Lutzomyia dysponeta. We detected Leishmania spp. infection in sand flies by polymerase chain reaction (PCR) amplification of the internal transcribed spacer region 1 (ITS-1) in pooled females (1-10 females per pool). Host species of engorged sand flies were identified using a cytochrome b PCR. From 455 sand fly pools analyzed, 255 pools were positive for Leishmania spp., with an estimated infection rate (confidence interval) of 0.096 [0.080-0.115] before the deltamethrin fogging which slightly, but not significantly (P > 0.05), increased to 0.116 [0.098-0.136] after the deltamethrin fogging. Blood meal analysis suggested that pigs, goats, and birds were the most common sand fly blood sources, followed by humans and domestic dogs. DNA sequencing from a subsample of ITS-1 positive pools suggests that Leishmania panamensis, Leishmania naiffi, and other Leishmania spp. were the parasite species infecting the most common vectors at the study site. Our data confirm an association between sand fly species, humans, domestic dogs, and pigs and Leishmania spp. parasites in rural Panamá.

Genetic structure of the Chagas disease vector Triatoma infestans (Hemiptera: Reduviidae) based on nuclear and mitochondrial DNA sequences

In the present study, we analysed the genetic structure of Triatoma infestans populations with a phylogeographical approach using sequences of the mitochondrial NADH dehydrogenase subunit 5 (ND5) and the nuclear elongation factor-1 alpha (EF-1α) genes of bugs obtained from Argentina and Bolivia. Spatially circumscribed haplogroups were distinguished from the ND5 gene sequences, one distributed exclusively to the south of the studied area and, in agreement with the results from the EF-1α gene, one haplogroup limited to Bolivia and another to Morajú located in the Chaco region of Argentina. In both the ND5 and EF-1α networks, the most widespread haplogroup or allele group showed a star-like topology, which is compatible with a recent demographic expansion. The asymmetric historical gene flow detected from a population of the Chaco region towards Bolivia and the spatiotemporal phylogeographical reconstruction of lineage dispersal would support the hypothesis that postulates the Chaco biogeographical region as the area of origin for the species. However, additional studies with a broader sampling in the Andean region are needed to define with certainty whether the origin of T. infestans is Chacoan or Andean.

Molecular validation of anthropophilic Phlebotominae sandflies (Diptera: Psychodidae) in Central Panama

Six Phlebotominae sand fly species are incriminated as biological vectors of human pathogens in Panama, but molecular corroboration is still needed. We aim at confirming the identity of Phlebotominae species documented as anthropophilic in Panama. Adult sandflies were collected from August 2010 to February 2012 in Central Panama using CDC light traps. Species confirmation was accomplished through molecular barcodes and allied sequences from GenBank. A total of 53,366 sand fly specimens representing 18 species were collected. Five species were validated molecularly as single phylogenetic clusters, but Psychodopygus thula depicted two genetically divergent lineages, which may be indicative of cryptic speciation.

Genetic diversity, phylogeography and molecular clock of the Lutzomyia longipalpis complex (Diptera: Psychodidae)

BACKGROUND

The Lutzomyia longipalpis complex has a wide but discontinuous distribution in Latin America, extending throughout the Neotropical realm between Mexico and northern Argentina and Uruguay. In the Americas, this sandfly is the main vector of Leishmania infantum, the parasite responsible for Visceral Leishmaniasis (VL). The Lu. longipalpis complex consists of at least four sibling species, however, there is no current consensus on the number of haplogroups, or on their divergence. Particularly in Argentina, there have been few genetic analyses of Lu. longipalpis, despite its southern expansion and recent colonization of urban environments. The aim of this study was to analyze the genetic diversity and structure of Lu. longipalpis from Argentina, and to integrate these data to re-evaluate the phylogeography of the Lu. longipalpis complex using mitochondrial markers at a Latin American scale.

METHODOLOGY/PRINCIPAL FINDINGS

Genetic diversity was estimated from six sites in Argentina, using a fragment of the ND4 and the 3´ extreme of the cyt b genes. Greatest genetic diversity was found in Tartagal, Santo Tomé and San Ignacio. There was high genetic differentiation of Lu. longipalpis in Argentina using both markers: ND4 (FST = 0.452, p < 0.0001), cyt b (FST = 0.201, p < 0.0001). Genetic and spatial Geneland analyses reveal the existence of two primary genetic clusters in Argentina, cluster 1: Tartagal, Santo Tomé, and San Ignacio; cluster 2: Puerto Iguazú, Clorinda, and Corrientes city. Phylogeographic analyses using ND4 and cyt b gene sequences available in GenBank from diverse geographic sites suggest greater divergence than previously reported. At least eight haplogroups (three of these identified in Argentina), each separated by multiple mutational steps using the ND4, are differentiated across the Neotropical realm. The divergence of the Lu. longipalpis complex from its most recent common ancestor (MRCA) was estimated to have occurred 0.70 MYA (95% HPD interval = 0.48-0.99 MYA).

CONCLUSIONS/SIGNIFICANCE

This study provides new evidence supporting two Lu. longipalpis genetic clusters and three of the total eight haplogroups circulating in Argentina. There was a high level of phylogeographic divergence among the eight haplogroups of the Lu. longipalpis complex across the Neotropical realm. These findings suggest the need to analyze vector competence, among other parameters intrinsic to a zoonosis, according to vector haplogroup, and to consider these in the design and surveillance of vector and transmission control strategies.

Francisella-Like Endosymbionts and Rickettsia Species in Local and Imported Hyalomma Ticks

Hyalomma ticks (Acari: Ixodidae) are hosts for Francisella-like endosymbionts (FLE) and may serve as vectors of zoonotic disease agents. This study aimed to provide an initial characterization of the interaction between Hyalomma and FLE and to determine the prevalence of pathogenic Rickettsia in these ticks. Hyalomma marginatum, Hyalomma rufipes, Hyalommadromedarii, Hyalommaaegyptium, and Hyalommaexcavatum ticks, identified morphologically and molecularly, were collected from different hosts and locations representing the distribution of the genus Hyalomma in Israel, as well as from migratory birds. A high prevalence of FLE was found in all Hyalomma species (90.6%), as well as efficient maternal transmission of FLE (91.8%), and the localization of FLE in Malpighian tubules, ovaries, and salivary glands in H. marginatum Furthermore, we demonstrated strong cophylogeny between FLE and their host species. Contrary to FLE, the prevalence of Rickettsia ranged from 2.4% to 81.3% and was significantly different between Hyalomma species, with a higher prevalence in ticks collected from migratory birds. Using ompA gene sequences, most of the Rickettsia spp. were similar to Rickettsiaaeschlimannii, while a few were similar to Rickettsiaafricae of the spotted fever group (SFG). Given their zoonotic importance, 249 ticks were tested for Crimean Congo hemorrhagic fever virus infection, and all were negative. The results imply that Hyalomma and FLE have obligatory symbiotic interactions, indicating a potential SFG Rickettsia zoonosis risk. A further understanding of the possible influence of FLE on Hyalomma development, as well as on its infection with Rickettsia pathogens, may lead to novel ways to control tick-borne zoonoses.IMPORTANCE This study shows that Francisella-like endosymbionts were ubiquitous in Hyalomma, were maternally transmitted, and cospeciated with their hosts. These findings imply that the interaction between FLE and Hyalomma is of an obligatory nature. It provides an example of an integrative taxonomy approach to simply differentiate among species infesting the same host and to identify nymphal and larval stages to be used in further studies. In addition, it shows the potential of imported Hyalomma ticks to serve as a vector for spotted fever group rickettsiae. The information gathered in this study can be further implemented in the development of symbiont-based disease control strategies for the benefit of human health.

Spatial Distribution of Sand Fly Vectors and Eco-Epidemiology of Cutaneous Leishmaniasis Transmission in Colombia

Background

Leishmania is transmitted by Phlebotominae insects that maintain the enzootic cycle by circulating between sylvatic and domestic mammals; humans enter the cycles as accidental hosts due to the vector’s search for blood source. In Colombia, leishmaniasis is an endemic disease and 95% of all cases are cutaneous (CL), these cases have been reported in several regions of the country where the intervention of sylvatic areas by the introduction of agriculture seem to have an impact on the rearrangement of new transmission cycles. Our study aimed to update vector species distribution in the country and to analyze the relationship between vectors’ distribution, climate, land use and CL prevalence.

Methods

A database with geographic information was assembled, and ecological niche modeling was performed to explore the potential distribution of each of the 21 species of medical importance in Colombia, using thirteen bioclimatic variables, three topographic and three principal components derived from NDVI. Binary models for each species were obtained and related to both land use coverage, and a CL prevalence map with available epidemiological data. Finally, maps of species potential distribution were summed to define potential species richness in the country.

Results

In total, 673 single records were obtained with Lutzomyia gomezi, Lutzomyia longipalpis, Psychodopygus panamensis, Psathyromyia shannoni and Pintomyia evansi the species with the highest number of records. Eighteen species had significant models, considering the area under the curve and the jackknife results: L. gomezi and P. panamensis had the widest potential distribution. All sand fly species except for Nyssomyia antunesi are mainly distributed in regions with rates of prevalence between 0.33 to 101.35 cases per 100,000 inhabitants and 76% of collection data points fall into transformed ecosystems.

Discussion

Distribution ranges of sand flies with medical importance in Colombia correspond predominantly to disturbed areas, where the original land coverage is missing therefore increasing the domiciliation potential. We highlight the importance of the use of distribution maps as a tool for the development of strategies for prevention and control of diseases.

Genetic structure analysis of Spirometra erinaceieuropaei isolates from central and southern China

Background

Sparganosis caused by invasion of the plerocercoid larvae (spargana) of Spirometra erinaceieuropaei have increased in recent years in China. However, the population genetic structure regarding this parasite is still unclear. In this study, we used the sequences of two mitochondrial genes cytochrome b (cytb) and cytochrome c oxidase subunit I (cox1) to analyze genetic variation and phylogeographic structure of the S. erinaceieuropaei populations.

Methodology/Principal Findings

A total of 88 S. erinaceieuropaei isolates were collected from naturally infected frogs in 14 geographical locations of China. The complete cytb and cox1 genes of each sample was amplified and sequenced. Total 61 haplotypes were found in these 88 concatenated sequences. Each sampled population and the total population have high haplotype diversity (Hd), accompanied by very low nucleotide diversity (Pi). Phylogenetic analyses of haplotypes revealed two distinct clades (HeN+HuN+GZ-AS clade and GX+HN+GZ-GY clade) corresponding two sub-networks yielded by the median-joining network. Pairwise F_ST values supported great genetic differentiation between S. erinaceieuropaei populations. Both negative Fu’s F_S value of neutrality tests and unimodal curve of mismatch distribution analyses supported demographic population expansion in the HeN+HuN+GZ-AS clade. The BEAST analysis showed that the divergence time between the two clades took place in the early Pleistocene (1.16 Myr), and by Bayesian skyline plot (BSP) an expansion occurred after about 0.3 Myr ago.

Conclusions

S. erinaceieuropaei from central and southern China has significant phylogeographic structure, and climatic oscillations during glacial periods in the Quaternary may affect the demography and diversification of this species.

American Cutaneous Leishmaniasis in Panama: a historical review of entomological studies on anthropophilic Lutzomyia sand fly species

We review existing information on the epidemiology of American Cutaneous Leishmaniasis (ACL) in Panama, with emphasis on the bionomics of anthropophilic Lutzomyia sand fly species. Evidence from Panamanian studies suggests that there are six anthropophilic species in the country: Lutzomyia trapidoi, Lu. panamensis, Lu. gomezi, Lu. ylephiletor, Lu. sanguinaria and Lu. pessoana (Henceforth Lu. carrerai thula). In general, these taxa are abundant, widespread and feed opportunistically on their hosts, which make them potential transmitters of pathogens to a broad range of wildlife, domesticated animals and humans. Furthermore, nearly all man-biting species in Panama (with the exception of Lu. gomezi) expand demographically during the rainy season when transmission is likely higher due to elevated Leishmania infection rates in vector populations. Despite this, data on the distribution and prevalence of ACL suggest little influence of vector density on transmission intensity. Apart from Lu. trapidoi, anthropophilic species seem to be most active in the understory, but vertical stratification, as well as their opportunistic feeding behavior, could vary geographically. This in turn seems related to variation in host species composition and relative abundance across sites that have experienced different degrees of human alteration (e.g., deforestation) in leishmaniasis endemic regions of Panama.