Research | Open | Published:
Phylogeography of the Lutzomyia gomezi (Diptera: Phlebotominae) on the Panama Isthmus
Parasites & Vectorsvolume 7, Article number: 9 (2014)
Lutzomyia gomezi (Nitzulescu, 1931) is one of the main Leishmania (Vianna) panamensis vectors in Panama, and despite its medical significance, there are no population genetic studies regarding this species. In this study, we used the sequences of the mitochondrial gene cytochrome b/start of NADH1 and the nuclear elongation gene α-1 in order to analyze genetic variation and phylogeographic structure of the Lu. gomezi populations.
A total of 86 Lu. gomezi individuals were captured in 38 locations where cutaneous leishmaniasis occurred. DNA was extracted with phenol/chloroform methods and amplification of genes was performed using PCR primers for mitochondrial and nuclear markers.
We found a total of 37 and 26 haplotypes of mitochondrial and nuclear genes, high haplotype diversity (h) for all three populations were detected with both molecular markers. Nucleotide diversity (π) was estimated to be high for all three populations with the mitochondrial marker, which was opposite to the estimate with the nuclear marker. In the AMOVA Φst recorded moderate (mitochondrial) and small (nuclear) population structure with statistical significance among populations. The analysis of the fixation index (Fst) used to measure the differentiation of populations showed that with the exception of the population located in the region of Bocas del Toro, the other populations presented with minor genetic differentiation. The median-Joining network of the mitochondrial marker reveled three clusters and recorded four haplotypes exclusively of localities sampled from Western Panama, demonstrating strong divergence. We found demographic population expansion with Fu´s Fs neutrality test. In the analysis mismatch distribution was observed as a bimodal curve.
Lu. gomezi is a species with higher genetic pool or variability and mild population structure, due to possible capacity migration and local adaptation to environmental changes or colonization potential. Thus, knowledge of the genetic population and evolutionary history is useful to understand the implications of different population genetic structures for cutaneous leishmaniasis epidemiology.
Cutaneous leishmaniasis is the most common form of leishmaniasis reported in the Republic of Panama; its clinical manifestations range from minor lesions to severe skin ulcers [1–3]. In Panama it was first recorded in 1910 , an overall of 15 cases occurred during 1910–1944; but there has been a sharp increase since 2000 . According to the Epidemiology Department of the Health Ministry of Panama, a total of 26,163 cases of cutaneous leishmaniasis occurred during 2000–2010 and the distribution of occurrence of CL was in regions of the provinces of Bocas de Toro (28%), Panama West (20%), Coclé (17%), Colón (11%), Panama East (5%), Veraguas (3%). The disease is mainly caused by Leishmania (Vianna) panamensis and Lutzomyia ylephiletor, Lu. sanguinaria, Lu. panamensis, Lu. trapidoi, Lu. gomezi were identified as the transmission vectors [6–8].
In Panama, recent research has show that Lutzomyia (Lutzomyia) gomezi (Nitzulescu, 1931) is the most abundant species with wide geographical distribution and their abundance has been associated with cases of clinical CL acquired in households of rural communities . Alongside Lu. gomezi other Lutzomyia species vectors have been caught in the same place and associated with a focus of cutaneous leishmaniasis infection in other American countries, incriminating this species with high potential in the transmission of the disease [2, 10–14].
Lutzomyia (Lutzomia) gomezi in Central America is present from Mexico to Panama and Trinidad Tobago (Caribbean), also recently reported in Guatemala . In South America, it has been reported in Colombia, Venezuela, Ecuador, Peru, French Guiana and Brazil, specifically in Amapá, Acre, Pará, Mata Grosso, Goiás, Bahia, Maranhão and Rondônia [15–17].
In general, its natural environment is a humid and dark place, such as nests, rock crevices, animal burrows and tree bark in the tropical rainforest . However, it is also reported in forest gaps and the canopies were light availability produces humid change [2, 18]. The foraging and blood feeding behavior occurs at twilight and at night from 18:00–20:00 hours . Nevertheless, the deforestation and loss of natural habitat have caused its adaptation to peridomestic areas, perhaps changing the hourly activity and feeding on a large variety of domestic animals [19, 20]. In the focus of leishmaniasis, this species is predominant near peridomestic or outdoor areas more than indoors households, which is difficult to control [2, 19, 21–23].
Despite of the potential significance of Lu. gomezi as a Leishmania vector, few researches have targeted genetic aspects or interrelationship of host-vector species. For instance,  detected a natural infection of Lu. gomezi with the Le. braziliensis in Venezuela, in order to establish a methods for determine the circulation of Leishmania parasites in leishmaniasis endemics areas. On the other hand,  analyzed the changes to the primary and secondary structures of tRNAser of the species; Lu. trinidadensis (Oswaldoi group), Lu. (Psychodopygus) panamensis, Lu. (Micropygomyia) cayennensis cayennensis, Lu. dubitans (Migonei group), Lu. (Lutzomyia) gomezi, Lu. rangeliana (ungrouped) and Lu. evansi (Verrucarum group) for taxonomic purposes, considering that morphological identification can be difficult. Also,  detected with ITS-1 a pool of Lu. gomezi infected with Le. naiffi in Panama, the first report for the country which prompted several hypotheses on the introduction of this parasite into this country.
Attempts to understand the role that arthropod vectors play in disease dynamics and pathogen transmission of leishmaniasis, several studies over genetic population of sandfly vectors has been performance in the Latin America, Iran, Turkey, Palestine, Israel and Egypt [27–30]. Many of these studies are focused in two principal species, Lu. longipalpis and Phlebotomus papatasi to determinate their genetic variation, structures and differentiation of populations [27, 28, 31, 32]. Thereby, molecular evidence for divergence of vectors has been found and an assessment of the impact on leishmaniasis epidemiology. For instance, the taxonomic status of the Lu. longipalpis complex has been fundamental to the understanding of leishmaniasis epidemiology. Symptoms observed for transmission of Leishmania infantum chagasi by Lu. longipalpis in Brazil and Colombia results in visceral infections, whereas the transmission of the same parasite by Lu. longipalpis in Costa Rica results in non-ulcerative lesions [33, 34].
Thus, it becomes important to analyze the population genetics and demographic history of this vector species to assess the genetic pool, gene flow and colonization potential, consequently known as the population changes and genetic variation affecting the vector competence and resistance to insecticides [35, 36]. On the other hand, the population genetics analysis can provide information about migration capacity, reproduction, and adaptation to the conditions of the new habitat, favoring the emergence of vector diseases . In addition to this, the determination of cryptic vector species and their ability to transmit pathogens are even more relevant to understand their implications in the epidemiology of vector diseases and to suggest appropriate and effective prevention and control programs without any environmental risks .
Due to the lack of population genetics information of Lu. gomezi species in Panama and the Americas as a whole, the goal of this research was to evaluate and compare the intra- and inter-population variability of Lu. gomezi from localities with a high incidence of cutaneous leishmaniasis in Panama. Also to identify the barriers that may influence gene flow among Lu. gomezi populations. Moreover, to infer the historical processes that defines actual geographic distributions of the species. These results may contribute on the knowledge of leishmaniasis epidemiology and improve the development of focal or large-scale programs control for leishmaniasis vectors in Panama.
The study was performed in the Republic of Panama between the coordinates 7º11′-9°39′N and 77º10′-83º03′W. The most prevalent climatic regime in Panama is tropical humid, with the dry season (January-March) presenting an average temperature of 31.5ºC and relative humidity of 75%. In other months (April-December), the average temperature is 27ºC and the relative humidity averages 90%. The vegetation in Panama varies according to climate zones and consists mainly of tropical humid forest or savannahs resulting from agricultural activity.
Sandfly collections and identification
Lu. gomezi samples were collected from the thirty-eight localities in Panama, during the dry season (January-April) and rainy season (May-June) of 2010. Thirty-four sites were chosen, representing incidence with cases of cutaneous leishmaniasis in humans during the 2006–2009 period, according to data from the Department of Epidemiology of the Ministry of Health of Panama; four sampling sites were in Panamanian tropical forests. The geographic distribution of the samples locations are indicated in Figure 1.
The characteristic of the thirty-four collecting sites according to the landscape was fragmented environment (secondary vegetation and areas used for agricultural activities) or rural (an area without a large concentration of people). The Coiba National Park, San Lorenzo Park, Oil Pipeline Road Park, and the Serrania Filo del Tallo Park represent natural reserves (forests) protected by the National Environment Authority of Panama (ANAM). The samples localities are listed in Table 1.
Specimens of Lu. gomezi used in the study were caught using CDC light traps and octenol solution utilized to attract hematophagous insects . The traps were positioned in the fragmented and rural locations, near peridomestic areas, with the presence of domestic animal (dogs, chickens, pigs, and cats) and ornamental vegetation. In the forest, traps were placed next of stone crevices and tree barks.
A total of nine CDC light traps were exposed at a height of 1.5 m in each locality and were installed at approximately 50 m intervals. The collection period was of 12 h (18:00 h-06:00 h) for two consecutive days at each point. The specimens were sacrificed with chloroform and stored in 95% ethanol at −20°C. Adult Lu. gomezi were identified by morphology and genital structure using the identification keys of .
DNA extraction and amplification
The DNA of 86 individual sandflies was isolated by a standard phenol/chloroform technique, precipitated with ethanol and resuspended in Tris-EDTA (TE) buffer according to , with some changes and stored at −20°C until used.
Two mitochondrial and nuclear genes were analyzed to estimate the structure and geographic genetics within and among populations of Lu. gomezi. The region of mitochondrial genes analyzed was final Cytb gene, inter-genetic region IGS-1, tRNA-Ser, inter-genetic region IGS-2, and start of the NADH1 (CB3-NIN), and amplified using primers CB3-PDR: 5′CAYATTCAACCWGAATGATA3′/N1N-PDR:5′GGCAYWTTGCCTCGAWTTCGWTATGA3′ . Nuclear region elongation factor alpha-1 (EF α-1) coding an intron region were amplified with primers EF-F03 5′CCTGGACATCGTGATTTCAT3′/EF-R04 5′AGTGCTTCGTGGTGTAT(C/T)TC3′ .
PCR amplification of genes were performed in a volume of 25 μl containing 1× buffer, 1.5 mM MgCl2++, 0.2 mM of each dNTPs, 0.4 μM of each primer, 1.5 units of Taq polymerase and approximately 2 μl of DNA (≈50 ng/μl) template. The details of thermal cycling conditions were outlined by  for CB3-N1N and  for EF α-1. A negative control was included in each PCR reaction and visualization of PCR product using 1.5% agarose gel with a 100 bp ladder.
Sequencing and alignment
PCR products were purified and sequenced directly by the Macrogen INC., Seoul, Korea Sequencing Service in both directions with the same primers used in the amplifications. The sequences were compared to the nucleotide BLAST (Basic Local Alignment Search Tool) tool available in GenBank NCBI (National Center for Biotechnology Information) in order to verify the similarities between the sequences of Lu. gomezi stored in the database. The CodonCode Aligner software was used to edit the nucleotide sequences and identify the heterozygous as double peaks on the chromatograms of the nuclear region. The sequences were aligned with the Muscle option included in the MEGA 5.05 software . The polymorphic sites, single substitution (singleton) and average base frequencies were computed with MEGA 5.05 software.
To carry out an analysis of population genetics (genetic diversity, structure and differentiation) to specimens collected in different localities, we clustered and defined samples in western, central and eastern populations from Panama (Table 1). Thus, the specimens collected in Bocas del Toro, Veraguas and Coiba represented the western populations, those collected from Coclé, Colon, and west of the Panama canal were central populations, whereas specimens from east of the Panama canal, Darién and San Blas corresponding to the eastern population (Figure 1).
We emphasized that heterozygotes detected in the EF α-1 were treated with recombination model in PHASE option, incorporated in the DnaSP v.5. To estimate the genetic diversity we used the following indexes: haplotypes number, haplotype diversity (h), and nucleotide diversity (π) for each population. DnaSP v.5  and Arlequin 3.11  software was used for this purpose.
Genetic structure and differentiation
The population structure was determined performing an Analysis of Molecular Variance (AMOVA) [45, 46]. The Arlequin 3.11  software analyzed the variance components and significance levels p = 0.01 using non-parametric permutations of 1,000 times. To estimate levels of genetic differentiation among the populations a pairwise comparison test was performed. Non-parametric permutations of 1,000 times and a significance p = 0.01 was used to estimate Fst-value  in Arlequin 3.11  software.
A haplotype network for both genes was inferred through the median-joining network method , using the Network v.4.5 software. In order to obtain the values of posteriori probabilities among the clades, a Bayesian analysis was calculated using MrBayes 3.1  program. MrModelTest  software was used to choose the best model of nucleotide substitution according to the AIC (Akaike Information Criterion). The model obtained by MrModeltest for CB3-N1N Bayesian inference analysis was HKY + I, with the proportion of invariable places (I) = 0.7594. For EF α-1, the model obtained was K80 + I, with a proportion of invariable places (I) = 0.6725. 10 million generations of Markov Chain Monte Carlo (MCMC) and burn-in of 5 million were used.
Demographic history of population
In addition to this, Tajima’s D  was used for estimate deviations from neutrality mutation or selective neutrality of these populations, P-value 0.01 was considered statistically significant. Also Fu’s Fs  were calculated for neutrality of mutation against population growth, hitchhiking and background selection, P-value 0.02 as significant. Both values Tajima’s D and Fu’s Fs were calculated using the software Arlequin 3.11  with permutations of 1,000 times.
The mismatch distribution was calculated for CB3-N1N and EF α-1 as supplementary test to evidence of demographic expansion of the population, through sudden expansion model . In Arlequin 3.11  we computed the sum of square deviations (SSD) and Harpending´s raggedness index to test the goodness-of-fit of the model to the data. The time since population expansion was inferred using CB3-NIN region by means of the equation t = τ /2 μ, where tau (τ) was obtained on the mismatch distribution outputs and a mutation rate (μ) of 1.1-1.2% per million years in Heliconius erato.
Overall, 86 sequences isolated of CB3-NIN and EF α-1 of Lu. gomezi individuals from 38 selected localities from Panama Isthmus were analyzed. Of these, an alignment of 501 nucleotides was obtained of CB3-N1N with thirty-nine polymorphic sites and nineteen single substitutions (singleton). The average base frequencies were A = 37%, C = 15%, G = 10%, and T = 38% recording a sequence rich in A-T (75%). A fragmented size of EF α-1 was 597 bp, in this case 166 sequences were obtained by clone generating with the recombination model algorithm of DnaSP v.5 considering the heterozygous individuals. Of those sequences, thirty-three polymorphic sites, six singletons and a larger proportion of G-C (53%): 23%T, 27%C, 23%A, 27%G, were detected.
For the analysis of genetic diversity indexes both CB3-NIN and EF α-1 region demonstrated high haplotype diversity (h) for all three populations [Table 2]. While the nucleotide diversity (π) for the mitochondrial region was estimated to be high for all three populations, nuclear nucleotide diversity was opposite to the estimate with mitochondrial, showing a low value [Table 2]. The haplotype and the landscape features per collection site are showed in [Table 1].
Genetic structure and differentiation
In the AMOVA Φst estimated between the three populations, CB3-N1N region recorded moderate and EF α-1 small population structure with statistical significance among populations. The total genetic variation found within the populations was ≈ 86% (CB3-NIN) and ≈ 95% (EF α-1) as shown in Table 3. Pairwise comparison analysis of the mitochondrial region revealed strong genetic differentiation with highly significance between Western Panama and Central Panama, whereas the Fst-value between Western and Eastern Panama, showed a moderate and significant differentiation. However, when populations were analyzed with nuclear region it revealed weak differentiation among populations of Western Panama and Eastern Panama. Also a moderate level of genetic differentiation was observed among the populations a Central Panama and Eastern Panama [Table 4].
The median-Joining network reconstructed with CB3-NIN region linked a total of 37 different haplotypes with three defined clusters (A, B, C). The lack of hierarchic structure is observed in the cluster “A”, this clade consists of 28 haplotypes, the most frequent, H10, was found in eleven localities sampled in Central Panama and six in Eastern Panama. It also included Altos de Piedra and Coiba haplotypes from Western Panama. Cluster “B” is constituted by four haplotypes. Of these, H15 was collected from six sites in Central Panama and nine from Eastern Panama, with H23 being the most frequent [Figure 2A]. In cluster “C” we recorded four haplotypes exclusively of localities sampled from Western Panama, demonstrating strong divergence support by Bayesian analysis [complementary data]. Thus, the network is characterized by few mutational steps [Figure 2A].
On the other hand, analysis linked a total of 26 haplotypes of the EF α-1 region demonstrating a lack of genetic structure in the Lu. gomezi populations. The patterns of a star-like shape as similarly obtained with the mitochondrial region, showed haplotype H1 as the most frequent in the center and several haplotypes surrounding it with one and three mutational steps [Figure 2B].
The Bayesian phylogenetic tree for the CB3-NIN region showed three clusters with a posteriori probability of 0.99-1.0. Other minor clusters were supported by low a posteriori probability values. The cluster with a higher a posteriori probability corroborated the divergences obtained in the haplotype network, inferred by the median-joining method [Figure 2A and B]. The Bayesian inference of EF α-1 region did not show similitude either, this analysis performed for mitochondrial and nuclear genes are supplementary data attached as Additional files 1 and 2.
We found 50% of the haplotypes in the peridomestic areas of fragmented environments, 40% in rural environments, and 10% in forest environments. Although this information is additional and not the main focus of the research, the predominance of the haplotypes in these sites is important to highlight due to the anthropogenic impact on the landscape of Panama [Table 1].
Demographic history population
The results of Tajima D test (CB3-NIN) were negative but not significant (−0.998; p > 0.01) that indicated neutrality deviation due to a demographic expansion or not neutral selection. Nevertheless, Fu´s Fs values were significantly negative (−19.68; p < 0.02) detecting that a population expanded in the past. The assessment of demographic change with EF α-1 region using both Tajima D and Fu´s Fs was significantly negative (−2.08; p < 0.01 and −13.65; p < 0.02) respectively.
Considering that Fu´s Fs are sensitive estimators for detecting demographic population expansion, we included the analysis of mismatch distribution for characterization of expansion detected with both genes. A bimodal curve was observed with CB3-NIN region and smoother curves with EF α-1 region. SSD value (0.011; p > 0.05) for CB3-NIN did not reject the sudden expansion model, while SSD value (0.049; p > 0.05) for EF α-1 region was observed. The Harpending´s raggedness index determined a good fit of data (Additional file 3). The time since expansion of the population was estimated to be approximately 4,223,825 years ago corresponding to the Pliocene Epoch.
The use of the mitochondrial and nuclear marker in ecology has provided enough information to allow contrast of the past population processes of species . The high content of A-T (74%) of the mitochondrial gene analyzed in Lu. gomezi is similar to that detected in other sandfly species: Phlebotomus papatasi and Lu. evansi[30, 55]. This suggests a high rate of substitutions as a result of natural selection, favorable or unfavorable to populations . A higher proportion of G-C (53%) observed in the nuclear gene in Lu. gomezi was comparable to the proportion verified in Anopheles gambiae (45%) and Culex pipiens (63%), characteristically of intron region [57, 58].
The high haplotype diversity detected among populations indicated a large polymorphism or gene pool of the Lu. gomezi. The existence of this polymorphism is important for the survival and adaptability, considering that the natural habitat of Phlebotominae sandflies is destroyed because of anthropogenic factors . However, we observed discordance in the level of nucleotide diversity between CB3-NINI and EF α-1 sequences. This may be because the mitochondrial gene is highly variable due to an elevated mutation rate and inheritence maternally . While the EF α-1 is an intron region of low copy number of protein coding highly conserved .
The result obtained from AMOVA showed mild genetic structure but greater genetic variation within populations. The explanation is that gene flow occurs between the different populations, probably caused by the flying capacity of the species increasing their migratory potential. Similar results were found from the analysis of population genetics of Lu. longipalpis, Lu. intermedia and Lu. shannoni over a distance between 8 to 15 km, which suggested that individuals from sandfly populations were capable of flying more than 57 m, as previously recorded for sandfly species [2, 60–64]. Thus, we hypothesized that Lu. gomezi is able to disperse over long distances, which could be gradually helped by the winds (adults) or passive movement of eggs and larvae. Studies across distant localities using other molecular markers are necessary to define vector capacity to disperse.
Meanwhile, pairwise Fst-value demonstrated significant genetic differentiation with the CB3-NIN region of the populations from Western Panama and Central Panama. We think that this genetic differentiation observed is influenced by the central mountain range with altitudes of 2,468 to 3,475 m. The mountain range represents a geographic barrier restricting altitudinal distribution of Lu. gomezi, estimated between 175–300 meters high . Similarly, strong genetic differentiation was found between Anopheles albimanus populations from Costa Rica and Panama, also suggesting that the Central American mountain range is a physical barrier that limited the gene pool of malaria vector mosquitoes [66, 67].
On the other side, slight genetic differences were detected from CB3-NIN and EF α-1 regions, among the sandfly populations from Eastern Panama and Central Panama (separated by the Panama Canal basin) could be explained because of the variety of climate and vegetation in Panama. The climate of Panama is influenced by the Inter-tropical Convergence Zone, oscillating one year to other and from one region to another. This variation allowed three types of tropical climate: Afi (very humid), Ami (humid), Awi (dry), according to Köppen Climate Classification . The sampling locations from Veraguas, Coclé, Colon, and Panama (west of the province of Panama) registered precipitations and temperatures between 3,001-3,500 mm and 22°C-24°C, respectively, while in the region of Eastern Panama, the amount of rain and the temperature reported varied between 1,000-1,800 mm and 26°C-27°C, respectively.
As a result, the type of tropical climate of Panama provides enough conditions to enable vegetation heterogeneity and refuge for a great diversity of species . Lu. gomezi samples from Central Panama were collected from borderline locations of tropical forest vegetation of low lands. While the Eastern Panama populations were from areas of vegetation in regeneration with evident anthropogenic impact. Thus, both climate and vegetation defined different ecotypes that can restrict gene flow due to environmental gradient. Climate changes were also used by  to justify the genetic differences found among the Brazilian populations of Lu. whitmani.
The star-like shape of network haplotypes, both for CB3-NIN and EF α-1, reflected a rapid and recent expansion of Lu. gomezi populations. The haplotypes of CB3-NIN (H-10, H-15, H-23) and haplotype EF α-1 (H-1) are older alleles or ancestral haplotypes dispersed in the Isthmus. The network patterns in clusters A and B for the CB3-NIN region and network EF α-1 region could be explained due to the sandflies migration with their hosts, mainly mammals.
Several mammal species (Canidae, Procyonidae, Sciuridae, Didelphidae, Dasypodidade) that currently are sandfly hosts established themselves in Panama during the north–south migration that took place during the closure of the Isthmus between the Miocene-Pliocene epochs [71–73]. We suggest this event as causal of the dispersal of Lu. gomezi because fossil records from the Miocene showed the presence of sandflies in fragments of mammal hair and in the microstructure of bird plumages in the Dominican Republic, suggesting their relationship with several hosts and possibly their transportation mechanism .
The H6 of CB3-NIN region from Coiba National Park (Veraguas) probably derived from the H11 detected in the populations of Central Panama, which indicates a possible migration of Lu. gomezi from the continent. Coiba was a continental island separated from South America in the Eocene, laying next to the Isthmus of Panama due to the decrease in the sea level in the last glaciation. This could have enabled the migration of wild animals and later, when the sea level increased, Coiba was isolated [75, 76]. These factors could have isolated this haplotype, and it is also evidence that wild hosts could help disperse these vectors.
The cluster B clearly showed the divergence in a west population (Bocas del Toro) and this strengthened significantly with the Fst-value [Table 4], we suppose that the population could have become geographically isolated. However other observations are necessary to establish phenotypic differences and corroborate this assumption.
The fragmentation of the forests on Panama is reported since pre-Columbian times [77, 78]. The prevalence of 10 (CB3-NIN) and 1 (EF α-1) ancestral haplotypes all types of environment and the percent of haplotypes found in the fragment and rural environments could indicate the displacement of Lu. gomezi because of disturbed areas. However, it is necessary to establish the consequences of fragmentation forest in the population genetics by increasing sampling in local forests.
The neutrality tests Fu´s Fs results were significantly negative, and mismatch distribution statistic and network haplotype support the assumption of population expansion in the past. Moreover, high haplotype diversity but low nucleotide diversity in EF α- region also suggests population expansion.
Is it the fact that the closure of the Isthmus of Panama affected by climate change in the world was the main factor that enabled the establishment of several species in America, allowing for faunal exchanges between North and South . An alternative explanation of the population expansion is that conditions later to closure of the isthmus may have been favorable to its adaptation to the suitable habitat and different ecological conditions in Panama territory. According to our estimated date the Pliocene climate was much cooler, allowing the expansion of plant and animal species into new habitats . The colonization and adaptation into new habitat can be observed in the star-like shape of network haplotype.
In conclusion, Lu. gomezi is a species with a higher genetic pool or variability and mild population structure, due to possible capacity migration and local adaptation to environmental changes or colonization potential. However the existence of geographic barriers, such as the central mountain range separates a subpopulation. The establishment of Lu. gomezi in the Isthmus of Panama could be interpreted by the exchange of mammals in the American continent, vegetation and climate conditions of the epoch. Their relevance as a carrier of cutaneous leishmaniasis in Panama is important for health schemes, especially if we consider that sandflies have a wide-spread geographic distribution in Panama, and this distribution is possibly associated to the landscape changes caused by deforestation . In this study, mitochondrial and nuclear analysis has provided important information to allow assembles about genetic population and evolutionary history useful to understand the implications of different population genetic structures for cutaneous leishmaniasis epidemiology. Thereby, proposing new perspectives for the control of leishmaniasis vectors.
Cytb gene, inter-genetic region IGS-1
Inter-genetic region IGS-2, and start of the NADH1
- EF α-1:
Elongation factor alpha-1.
Christensen HA, Herrer A, Telford SR: Enzootic cutaneous leishmaniasis in eastern Panama. II. Entomological investigations. Ann Trop Med Parasitol. 1972, 66: 55-66.
Christensen H, Fairchild GB, Herrer A, Johnson C, Young D, de Vásquez A: The ecology of cutaneous, leishmaniasis in the Republic of Panama. J Med Entomol. 1983, 20: 463-484.
Christensen H, Johnson C, Vasquez AM: Leishmaniasis cutánea en Panamá: Un breve resumen. Rev Med Panama. 1984, 9: 182-1987.
Darling ST: Oriental sore in Panama. Proc Canal Zone Med Assn. 1910, 2: 7-20.
Miranda A, Carrasco R, Paz H, Pascale JM, Samudio F, Saldaña A, Santamaría G, Mendoza Y, Calzada JE: Molecular epidemiology of American tegumentary leishmaniasis in Panama. Am J Trop Med Hyg. 2009, 81: 565-571. 10.4269/ajtmh.2009.08-0265.
Herrer A, Telford S, Christensen HA: Enzootic cutaneous leishmaniasis in eastern Panama I: investigation of the infection among forest mammals. Ann Trop Med Parasitol. 1971, 65: 349-358.
Tesh RB, Chaniotis BN, Aronson MD, Johnson KM: Natural host preferences of Panamanian Phlebotominae sandflies as determined by precipitin test. Am J Trop Med Hyg. 1971, 20: 150-156.
Herrer A, Christensen HA: Epidemiological patterns of cutaneous leishmaniasis in Panama. Am J Trop Med Hyg. 1973, 25: 54-58.
Saldaña A, Chaves LF, Rigg CA, Wald C, Smucker JE, Calzada JE: Clinical cutaneous leishmaniasis rates are associated with household Lutzomyia gomezi, Lu. Panamensis, and Lu. trapidoi Abundance in Trinidad de Las Minas, Western Panama. Am J Trop Med Hyg. 2013, 88 (3): 579-
Jaramillo C, Travi BL, Montoya J: Vector competence of some neotropical sandflies for the Leishmania (Viannia) braziliensis complex. Med Vet Entomol. 1994, 8: 1-7. 10.1111/j.1365-2915.1994.tb00375.x.
Jiménez AE, Rojas JC, Vargas F, Herrero MV: Temporal and spatial variation of phlebotomine (Diptera: Psychodidae) community diversity in a cutaneous leishmaniasis endemic area of Costa Rica. J Med Entomol. 2000, 37: 216-221. 10.1603/0022-2585-37.2.216.
Bejarano EE, Uribe S, Rojas W, Vélez ID: Phlebotomine sand flies (Diptera: Psychodidae) associated with the appearance of urban leishmaniasis in the city of Sincelejo, Colombia. Mem Inst Oswaldo Cruz. 2002, 97: 645-647. 10.1590/S0074-02762002000500010.
Santamaría E, Ponce N, Zipa Y, Ferro C: Presencia en el peridomestic de vectores infectados con Leishmania (Viannia) panamensis en dos focos endémicos en el occidente de Boyacá, Piedemonte del valle del Magdalena medio, Colombia. Biomedica. 2006, 26: 82-94.
Vivero RJ, Contreras MA, Cadena H, Acosta LA, Mondragon K, Vélez A, Uribe SI, Vélez ID: Sand flies Phlebotominae (Diptera: Psychodidae) in five Central American countries. International Symposium on Phlebotominae Sandflies ISOPS VII: 25–30 april. 2011, Kudasy-Turkey, 95-
Martins AV, Williams P, Falcão AL: American sand flies (Diptera: Psychodidae: Phlebotominae). 1978, Rio de Janeiro: Acad Bras Ciênc, 195-
Young DG, Duncan MA: Guide to the identification and geographic distribution of Lutzomyia sandflies in Mexico, the West Indies, Central and South America (Diptera: Psychodidae). Mem Amer Ent Inst. 1994, 54: 1-881.
Galati EAB: Classificação de Phlebotominae. Flebotomíneos do Brasil. Edited by: Rangel ER, Lainson R. 2003, Rio de Janeiro: Editora Fiocruz Brazil, 23-52.
Chaniotis BN, Correa MA, Tesh RB, Jonhson KM: Daily and seasonal man-biting activity of Phlebotominae sandflies in Panama. J Med Entomol. 1971, 8: 415-420.
Travi BL, Alder GH, Lozano M, Cadena H, Montoya-Lerma J: Impact of habitat degradation on Phlebotominae (Diptera: Psychodidae) of tropical dry forest in northern Colombia. J Med Entomol. 2002, 39: 451-456. 10.1603/0022-2585-39.3.451.
Valderrama A, Tavares MG, Andrade Filho DJ: Anthropogenic influence on the distribution, abundance and diversity of sandfly species (Diptera: Phlebotominae: Psychodidae), vectors of cutaneous leishmaniasis in Panama. Mem Inst Oswaldo Cruz. 2011, 106: 1024-1031. 10.1590/S0074-02762011000800021.
Feliciangeli MD, Rabinovich J: Abundance of Lutzomyia ovallesi but not Lu. gomezi (Diptera: Psychodidae) correlated with cutaneous leishmaniasis incidence in north-central Venezuela. Med Vet Entomol. 1988, 12: 121-131.
Duque P, Vélez ID, Morales M, Sierra D: Sand flies fauna involved in the transmission of cutaneous leishmaniasis in Afro-Colombian and Amerindian communities of Choco, Pacific Coast of Colombia. Neotrop Entomol. 2004, 33: 255-264.
Cortés LA, Fernández JJ: Especies de Lutzomyia en un foco urbano de leishmaniasis visceral y cutánea en El Carmen de Bolívar, Bolívar, Colombia. Biomedica. 2008, 28: 433-440.
Jorquera A, González R, Marchán-Marcano E, Oviedo M, Matos M: Multiplex-PCR for detection of natural Leishmania infection in Lutzomyia spp captured in an endemic region for cutaneous leishmaniasis in state of Sucre, Venezuela. Mem Inst Oswaldo Cruz. 2005, 100: 43-46.
Vivero RJ, Contreras-Gutiérrez MA, Bejarano EE: Análisis de la estructura primaria y secundaria del ARN de transferencia mitocondrial para serina en siete especies de Lutzomyia. Biomedica. 2007, 27: 429-438.
Azpurua J, de la Cruz D, Valderrama A, Windsor D: Lutzomyia sandfly diversity and rates of infection by Wolbachia and an exotic Leishmania species on Barro Colorado Island, Panama. PLoS Negl Trop Dis. 2010, 4: e627-10.1371/journal.pntd.0000627.
Ready PD, Day JC, de Souza AA, Rangel EF, Davies CR: Mitochondrial DNA characterization of populations of Lutzomyia whitmani (Diptera: Psychodidae) incriminated in the peridomestic and silvatic transmission of Leishmania species. Bull Entomol Res. 1997, 87: 187-195. 10.1017/S0007485300027346.
Soto SI, Lehmann T, Rowton ED, Velez BID, Porter CH: Speciation and population structure in the morphospecies Lutzomyia longipalpis (Lutz & Neiva) as derived from the mitochondrial ND4 gene. Mol Phylogen Evol. 2001, 18: 84-93. 10.1006/mpev.2000.0863.
Hodgkinson VH, Birungi J, Quintana M, Deitze R, Munstermann LE: Mitochondrial Cytocrome B variation in population of the visceral leishmanisis vector Lutzomyia longipalpis across eastern Brazil. Am J Trop Med Hyg. 2003, 69: 386-392.
Hamarsheh O, Presber W, Abdeen Z, Sawalha S, Al-Lahem A, Schoenian G: Genetic structure of mediterranean populations of the sandfly Phlebotomus papatasi by mitochondrial cytochrome b haplotype analysis. Med Vet Entomol. 2007, 21: 270-277. 10.1111/j.1365-2915.2007.00695.x.
Esseghir S, Ready PD, Killick-Kendrick R, Ben-Ismail R: Mitochondrial haplotypes and phylogeography of Phlebotomus vectors of Leishmania major. Insect Mol Biol. 1997, 6: 211-225. 10.1046/j.1365-2583.1997.00175.x.
Depaquit J, Lienard E, Verzeaux-Griffon A, Ferte A, Bounamous H, Gantier JC, Hanafi HA, Jacobson RL, Maroli M, Moin-Vaziri V, Muller F, Ozbel Y, Svobodova M, Volf P, Leger N: Molecular homogeneity in diverse geographical populations of Phlebotomus papatasi (Diptera, Psychodidae) inferred from ND4 mtDNA and ITS2 rDNA Epidemiological consequences. Infect Genet Evol. 2008, 8: 159-170. 10.1016/j.meegid.2007.12.001.
Zeledón R, Hidalgo H, Víquez A, Urbina A: Atypical cutaneous leishmaniasis in a semiarid region of north-west Costa Rica. Trans R Soc Trop Med Hyg. 1989, 83: 786-10.1016/0035-9203(89)90328-3.
Warburg A, Saraiva E, Lanzaro GC, Titus R, Neva F: Saliva of Lutzomyia longipalpis sibling species differs in its composition and propensity to enhance leishmaniasis. Phil Trans Roy Soc Lond. 1994, 345: 223-230.
Lourenço-de-Oliveira R, Vazeille M, Fillipis AMB, Failloux A: Aedes aegypti in Brazil: genetically differentiated populations with high susceptibility to dengue and yellow fever viruses. Tran R Soc Trop Med Hyg. 2004, 98: 43-54. 10.1016/S0035-9203(03)00006-3.
Black WC, Tabachnick WJ: Population genetics in vector biology. The Biology of Disease Vectors. Edited by: Beaty BJ, Marquardt WC. 2004, Colorado: University Press of Colorado, Niwot, 417-437.
Reisen WK, Lothrop HD, Presser SB, Hardy JL, Gordon EW: Landscape ecology of arboviruses in southeastern California: temporal and spatial patterns of enzootic activity in Imperial Valley, 1991–94. J Med Entomol. 1997, 34: 179-188.
Sudia WD, Chamberlain RW: Battery operated light trap, an improved model. Mosq News. 1962, 22: 126-129.
Michalsky EM, Fortes-Dias CL, Pimenta PFP, Secundino NFC, Dias ES: Assessment of PCR in the detection of Leishmania spp in experimentally infected individual phlebotomine sandflies (Diptera: Psychodidae: Plebotominae). Rev Inst Med Trop Sao Paulo. 2002, 44: 255-259.
Testa JM, Montoya-Lerna J, Cadena H, Oviedo M, Ready PD: Molecular identification of vectors of Leishmania in Colombia: mitochondrial introgression in the Lutzomyia townsendi series. Acta Trop. 2002, 84: 205-218. 10.1016/S0001-706X(02)00187-0.
Parvizi P, Assmar M: Nuclear elongation factor 1-α gene, A molecular marker for Iranian sandfly identification. Iranian J Publ Health. 2007, 36: 25-37.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28: 2731-2739. 10.1093/molbev/msr121.
Librado P, Rozas J: DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009, 25: 1451-1452. 10.1093/bioinformatics/btp187.
Excoffier L, Laval G, Schneider S: Arlequin ver. 3.0: An integrated software package for population genetics data analysis. Evol Bioinform Online. 2005, 1: 47-50.
Weir BS, Cockerham C: Estimating F-statistics for the analysis of population structure. Evolution. 1984, 38: 1358-1370. 10.2307/2408641.
Excoffier L, Smouse PE, Quattro JM: Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data. Genetics. 1992, 131: 479-491.
Bandelt HJ, Forster P, Röhl A: Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999, 16: 37-48. 10.1093/oxfordjournals.molbev.a026036.
Huelsenbecr JP, Ronquist , MrBayes F: Bayesian inference of phylogenetics tree. Bioinformatics. 2001, 17: 754-755. 10.1093/bioinformatics/17.8.754.
Nylander JAA, Ronquist F, Huelsenbeck JP, Nieves-Aldrey JL: Bayesian phylogenetic analysis of combined data. System Biol. 2004, 53: 47-67. 10.1080/10635150490264699.
Tajima F: Evolutionary relationship of DNA sequences in finite populations. Genetics. 1983, 105: 437-460.
Fu YX: Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics. 1997, 147: 915-925.
Rogers AR, Harpending HC: Population growth makes waves in the distribution of pairwise genetic differences. Mol Biol Evol. 1992, 9: 552-569.
Brower AVZ: Phylogeny of Heliconius butterflies inferred from mitochondrial DNA sequences (Lepidoptera: Nymphalidae). Mol Phylogenet Evol. 1994, 3: 159-174. 10.1006/mpev.1994.1018.
Galtier N, Nabholz B, Glémin S, Hurst GD: Mitochondrial DNA as a marker of molecular diversity: a reappraisal. Mol Ecol. 2009, 18: 4541-4550. 10.1111/j.1365-294X.2009.04380.x.
Bejarano EE, Rojas W, Uribe S, Vélez ID: Genetic analysis of a recently detected urban population of Lutzomyia evansi (Diptera: Psychodidae) in Colombia. Rev Soc Entomol Argent. 2009, 68: 135-141.
Hedrick PW: Sex differences in mutation, recombination, selection, gene flow, and genetic drift. Evolution. 2007, 61 (12): 2750-2771. 10.1111/j.1558-5646.2007.00250.x.
Calcagnotto D: Taxa de evolução e o relógio molecular. Biologia Molecular e Evolução. Edited by: Matioli SR. 2001, Ribeirão Preto São Paulo: Holos Editora, 52-63.
Nahum LA: Evolução dos genomas. Biologia Molecular e Evolução. Edited by: Matioli SR. 2011, Ribeirão Preto SP: Holos Editora, 82-96.
Simon C, Frati F, Beckenbach A, Crespi B, Liu H, Flook P: Evolution, weighting and phylogenetic utility of mitochondrial gene sequences and a compilation of conserved polymerase chain reaction primers. Ann Entomol Soc Am. 1994, 87: 651-701.
Alexander JB, Young DG: Dispersal of phlebotomine sand flies (Diptera: Psychodidae) in a Colombian focus of Leishmania (Viannia) braziliensis. Mem Inst Oswaldo Cruz. 1992, 87: 397-403. 10.1590/S0074-02761992000300010.
Morrison AC, Ferro C, Morales A, Tesh RB, Wilson ML: Dispersal of the sand fly Lutzomyia longipalpis (Diptera: Psychodidae) at an endemic focus of visceral leishmaniasis in Colombia. J Med Entomol. 1993, 30: 427-435.
Márquez LM, Lampo M, Rinaldi M, Lau P: Gene flow natural and domestic populations of Lutzomyia longipalpis (Diptera: Psychodidae) in a restricted focus of American visceral leishmaniasis in Venezuela. J Med Entomol. 2001, 38: 12-16. 10.1603/0022-2585-38.1.12.
De Queiroz BV, Coutinho-Abreu IV, Sonoda IV, Melo MA, de Andrade PP, de Castro JA, Rebelo JM, Carvalho SM, Ramalho-Ortigao M: Genetic structure of natural populations of the sand fly Lutzomyia longipalpis (Diptera: Psychodidae) from the Brazilian northeastern region. Acta Trop. 2006, 98 (1): 15-24. 10.1016/j.actatropica.2006.01.007.
Rocha LS, Falqueto A, dos Santos CB, Ferreira AL, da Graça CC, Grimaldi G, Cupolillo E: Survey of natural infection by Leishmania in sand fly species collected in southeastern Brazil. Trans R Soc Trop Med Hyg. 2010, 104: 461-466. 10.1016/j.trstmh.2010.02.005.
Añez N, Cazorla D, Nieves E, Chataing B, Castro M, de Yarbuh AL: Epidemiología de la leishmaniasis tegumentaria en Mérida, Venezuela. I. Diversidad y dispersión de especies flebotominas en tres pisos altitudinales y su posible rol en la transmisión de la enfermedad. Mem Inst Oswaldo Cruz. 1988, 83: 455-463. 10.1590/S0074-02761988000400011.
Molina-Cruz A, de Mérida AM, Mills K, Rodríguez F, Schoua C, Yurrita MM, Molina E, Palmieri M, Black WC: Gene Flow Among Anopheles albimanus populations in Central América, South América and the Caribbean Assessed by microsatellites and mitochondrial DNA. Am J Trop Med Hyg. 2004, 71: 350-359.
Loaiza JR, Bermingham E, Scott ME, Rovira JR, Conn JE: Species composition and distribution of adult Anopheles (Diptera: Culicidae) in Panama. J Med Entomol. 2008, 45: 841-851. 10.1603/0022-2585(2008)45[841:SCADOA]2.0.CO;2.
Condit R, Hubbell SP, Foster RB: Changes in a tropical forest with a shifting climate: results from a 50 ha permanent census plot in Panama. J Trop Ecol. 1996, 12: 231-256. 10.1017/S0266467400009433.
Ready PD, de Souza AA, Macario Rebelo JM, Day JC, Silveira FT, Campbell-Lendrum D, Davies CR, Costa JML: Phylogenetic species and domesticity of Lutzomyia whitmani at the southeast boundary of Amazonian Brazil. Trans R Soc Trop Med Hyg. 1998, 92: 159-160. 10.1016/S0035-9203(98)90726-X.
Coates AG, Obando JA: Geological evolution of the Central American Isthmus. Evolution and environment in tropical America. Edited by: Jackson JBC, Budd AF, Coates AG. 1996, Chicago: University of Chicago Press, 21-56.
Cox CB: Plate tectonics, seaways and climate in the historical biogeography of mammals. Mem Inst Oswaldo Cruz. 2000, 95: 509-516.
Webb SD: The great American biotic interchange: patterns and processes. Ann Mo Bot Gard. 2006, 93: 245-257. 10.3417/0026-6493(2006)93[245:TGABIP]2.0.CO;2.
Peñalver E, Grimaldi D: Assemblages of mammalian hair and blood-feeding midges (Insecta: Diptera: Psychodidae: Phlebotominae) in Miocene amber. Trans R Soc Edin: Earth Sci. 2005, 96: 177-195.
Heckadon MS: Panamá: Puente biológico. Panama. 2001, Panamá: Smithsonian Tropical Research Institute, 233-
Hoernle KA, Bogaard PVD, Werner R, Lissinna B, Hauff G: Missing history (16–71 Ma) of the Galapagos hotspot: Implications for the tectonic and biological evolution of the Americas. Geology. 2002, 30: 795-798. 10.1130/0091-7613(2002)030<0795:MHMOTG>2.0.CO;2.
Linares OF, Sheets PD, Rosenthal EJ: Prehistoric agriculture in tropical highlands. Science. 1975, 187: 137-145. 10.1126/science.187.4172.137.
Piperno DR, Jones J: Paleoecological and archaeological implications of a late pleistocene/early holocene record of vegetation and climate from the Pacific coastal plain of Panama. Quaternary Res. 2003, 59: 79-87. 10.1016/S0033-5894(02)00021-2.
Marshall LG: Land Mammals and the Great American Interchange. Am Sci. 1988, 76 (4): 380-388.
Ravelo AC, Andreasen DH, Lyle M, Lyle AO, Wara MW: Regional climate shifts caused by gradual global cooling in the Pliocene epoch. Nature. 2004, 429: 263-267. 10.1038/nature02567.
We are grateful to Maykon Cristiano Passos, for helping with laboratory and analyses data; Jorge Edergarm Santos, Gustavo Martins and Juliana Lopes Fietto for comments in previous drafts of this manuscript; furthermore we thank the anonymous reviewers and editor for their comments.
This research was supported by the Secretaria Nacional de Ciencias y Tecnología Innovadora (SENACYT) Grant Col09-008.
The authors declare that they have no competing interests.
Conceived and designed the experiments: AVC, MGT. Performed the experiments: AVC. Analyzed the data: AVC. Contributed reagents/materials/analysis tools/specimens identifications: MGT, JDAF. Wrote the paper: AVC, JDAF and MGT. All authors read and approved the final version of the manuscript.