Genetic diversity and population structure in the Leishmania guyanensis vector Lutzomyia anduzei (Diptera, Psychodidae) from the Brazilian Amazon

Population structure analysis of Phlebotomus papatasi populations using transcriptome microsatellites: possible implications for leishmaniasis control and vaccine development

BACKGROUND

Phlebotomus papatasi is considered the primary vector of Leishmania major parasites that cause zoonotic cutaneous leishmaniasis (ZCL) in the Middle East and North Africa. Phlebotomus papatasi populations have been studied extensively, revealing the existence of different genetic populations and subpopulations over its large distribution range. Genetic diversity and population structure analysis using transcriptome microsatellite markers is important to uncover the vector distribution dynamics, essential for controlling ZCL in endemic areas.

METHODS

In this study, we investigated the level of genetic variation using expressed sequence tag-derived simple sequence repeats (EST-SSRs) among field and colony P. papatasi samples collected from 25 different locations in 11 countries. A total of 302 P. papatasi sand fly individuals were analyzed, including at least 10 flies from each region.

RESULTS

The analysis revealed a high-level population structure expressed by five distinct populations A through E, with moderate genetic differentiation among all populations. These genetic differences in expressed genes may enable P. papatasi to adapt to different environmental conditions along its distribution range and likely affect dispersal.

CONCLUSIONS

Elucidating the population structuring of P. papatasi is essential to L. major containment efforts in endemic countries. Moreover, the level of genetic variation among these populations may improve our understanding of Leishmania-sand fly interactions and contribute to the efforts of vaccine development based on P. papatasi salivary proteins.

Phylogeographic inference for Bichromomyia flaviscutellata sensu stricto (Diptera: Psychodidae: Phebotominae) from the Brazilian Amazon rainforest, based on the 3' region of the COI gene

Bichromomyia flaviscutellata (Mangabeira, 1942) sensu stricto (Diptera: Psychodidae) has been recognized as the main vector of Leishmania amazonensis in the Brazilian Amazon. For this reason, it is of paramount importance to understand the distribution of genetic diversity of populations of this vector, particularly the genetic structure and gene flow, for its management and control efforts. This study investigated the phylogeographic structure of five B. flaviscutellata s.s. populations from the central Brazilian Amazon region by analyzing 1,141 bp fragment of the 3' region of the COI gene. A total of 85 specimens of B. flaviscutellata s.s. were sequenced from Manaus (14), Rio Preto da Eva (10), Pitinga (14), Novo Airão (21), and Autazes (26); all in the state of Amazonas. The dataset yielded 59 haplotypes, most of them connected to each other in the main network. There were high levels of intrapopulation genetic variability (h = 0.945 ± 0.035 - 0.978 ± 0.054). The genetic distance values among populations varied from moderate (0.0873) to very high (0.3535), and all comparisons were significant, as well as the hierarchical analysis (ΦST = 0.2145). In contrast, these comparisons revealed a high number of shared sites (Ss = 6-34) and no difference in fixed sites (Sf = 0) among populations indicating absence of historical isolation. The Mantel test indicated that 67.92% (r = 0.6792; P = 0.06) of the genetic structure observed in B. flaviscutellata s.s. cannot be explained by the isolation-by-distance (IBD) model. This genetic structure, weakly explained by the IBD, may be due mainly by the forest habitat fragmentation and the low dispersal (flight) capacity of sand flies. Both factors could lead to population fragmentation and isolation, which promote genetic differentiation. Taken together, these findings suggest that the genetic structure observed in the studied populations of B. flaviscutellata s.s. is likely generated by microevolutionary processes acting at the population level at the present time and, therefore, evolutionary lineages were not recognized among the populations analyzed.

Genetic variability highlights the invasion route of the Lutzomyia longipalpis complex, the main vector of Visceral Leishmaniasis in Uruguay

In the Americas, the sandfly Lutzomyia longipalpis is the main vector of the parasitic protozoa Leishmania infantum, the etiological agent of visceral leishmaniasis (VL). The Lu. longipalpis species complex is currently discontinuously distributed across the Neotropical region, from Mexico to the north of Argentina and Uruguay. During its continental spreading, it must have adapted to several biomes and temperature amplitudes, when founder events should have contributed to the high genetic divergence and geographical structure currently observed, reinforcing the speciation process. The first report of Lu. longipalpis in Uruguay was in 2010, calling the attention of Public Health authorities. Five years later, the parasite Le. infantum was recorded and in 2015 the first case of VL in canids was reported. Hitherto seven human cases by VL have been reported in Uruguay. Here, we publish the first DNA sequences from the mitochondrial genes ND4 and CYTB of Lu. longipalpis collected in Uruguay, and we used these molecular markers to investigate their genetic variability and population structure. We described four new ND4 haplotypes in a total of 98 (4/98) and one CYTB in a total of 77 (1/77). As expected, we were able to establish that the Lu. longipalpis collected in two localities (i.e. Salto and Bella Unión) from the north of Uruguay are closely related to the populations from neighbouring countries. We also propose that the possible route for the vector arrival to the region may have been through vegetation and forest corridors of the Uruguay River system, as well as it may have benefited from landscape modifications generated by commercial forestation. The ecological-scale processes shaping Lu. longipalpis populations, the identification of genetically homogeneous groups and the gene flow among them must be carefully investigated by using highly sensible molecular markers (i.e. genome wide SNPs) since it will help to the understanding of VL transmission and contribute to the planification of public policies on its control.

Multiple evolutionary lineages for the main vector of Leishmania guyanensis, Lutzomyia umbratilis (Diptera: Psychodidae), in the Brazilian Amazon

Lutzomyia umbratilis is the main vector of Leishmania guyanensis in the Brazilian Amazon and in neighboring countries. Previous biological and molecular investigations have revealed significant differences between L. umbratilis populations from the central Brazilian Amazon region. Here, a phylogeographic survey of L. umbratilis populations collected from nine localities in the Brazilian Amazon was conducted using two mitochondrial genes. Statistical analyses focused on population genetics, phylogenetic relationships and species delimitations. COI genetic diversity was very high, whereas Cytb diversity was moderate. COI genealogical haplotypes, population structure and phylogenetic analyses identified a deep genetic differentiation and three main genetic groups. Cytb showed a shallower genetic structure, two main haplogroups and poorly resolved phylogenetic trees. These findings, allied to absence of isolation by distance, support the hypothesis that the Amazon and Negro Rivers and interfluves are the main evolutionary forces driving L. umbratilis diversification. The main three genetic groups observed represent three evolutionary lineages, possibly species. The first lineage occurs north of the Amazon River and east of Negro River, where Le. guyanensis transmission is intense, implying that L. umbratilis is an important vector there. The second lineage is in the interfluve between north of Amazon River and west of Negro River, an area reported to be free of Le. guyanensis transmission. The third lineage, first recorded in this study, is in the interfluve between south of Amazonas River and west of Madeira River, and its involvement in the transmission of this parasite remains to be elucidated.

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.

Behavior and population structure of Anopheles darlingi in Colombia

Anopheles darlingi is a widely distributed and important malaria vector in Colombia. Biogeographical and ecological heterogeneity across the Colombian distribution led to the hypothesis of behavioral and genetic differentiation among A. darlingi populations. A total of 2017 A. darlingi specimens were collected during 222 h of sampling. This vector was the most abundant anopheline species in most of the localities sampled. Subdivision between samples collected west and east of the Andes was indicated by 1) mitochondrial COI and nuclear CAD sequences from NW-W and CE-S populations (COI ΦST=0.48761-0.81974, CAD FST=0.11319-0.21321), 2) a COI haplotype network, and 3) SAMOVA. Endo- and exophagy were detected in populations west of the Andes, whereas exophagy was evident in PTG, a locality east of the Andes. Isolation by resistance was significant for COI and explained 26% of the genetic differentiation. We suggest that at a macrogeographic scale, the Andes influence the differentiation of A. darlingi in Colombia and may drive divergence, and, at a microgeographic scale, ecological differences have a significant impact on structure. These data could constitute a baseline for the design of effective vector interventions, locality-specific for the east and similar for panmictic populations west of the Andes.