- Open Access
Lutzomyia longipalpis s.l. in Brazil and the impact of the Sao Francisco River in the speciation of this sand fly vector
Parasites & Vectorsvolume 1, Article number: 16 (2008)
Lutzomyia longipalpis s.l. (Diptera: Psychodidae) is the principal vector of Leishmania infantum chagasi in the Americas, and constitutes a complex of species. Various studies have suggested an incipient speciation process based on behavioral isolation driven by the chemotype of male sexual pheromones. It is well known that natural barriers, such as mountains and rivers can directly influence population divergence in several organisms, including insects. In this work we investigated the potential role played by the Sao Francisco River in eastern Brazil in defining the current distribution of Lu. longipalpis s.l. Our studies were based on analyses of polymorphisms of the cytochrome b gene (cyt b) sequences from Lu. longipalpis s.l. available in public databases, and from additional field-caught individuals. Altogether, 9 distinct populations and 89 haplotypes were represented in the analyses. Lu. longipalpis s.l. populations were grouped according to their distribution in regards to the 10°S parallel: north of 10°S (<10°S); and south of 10°S (>10°S). Our results suggest that although no polymorphisms were fixed, moderate genetic divergences were observed between the groups analyzed (i.e., F ST = 0.184; and Nm = 2.22), and were mostly driven by genetic drift. The population divergence time estimated between the sand fly groups was about 0.45 million years (MY), coinciding with the time of the change in the course of the Sao Francisco River, during the Mindel glaciation. Overall, the polymorphisms on the cyt b haplotypes and the current speciation process detected in Lu. longipalpis s.l. with regards to the distribution of male sexual pheromones suggest a role of the Sao Francisco River as a significant geographical barrier in this process.
Lutzomyia longipalpis s.l. (Diptera: Psychodidae) is the vector of Leishmania infantum chagasi, the causative agent of visceral leishmaniasis (VL) in the New World [1–3]. This sand fly species has a wide, though discontinuous distribution ranging from southern Mexico to northern Argentina . The pattern of distribution of this sand fly is directly associated with a notable population divergence due to a reduced gene flow, allowing the appearance of sibling species . In the last few years many investigators suggested the existence of at least four sibling species of Lu. longipalpis s.l. in South and Central Americas, based on the study of various chromosomal and molecular markers [6–11].
Recent studies also have suggested the presence of sibling species in eastern Brazil [12, 13]. These studies relied on the characterization of different male sexual pheromones, the "love songs" resulted from male wing vibration during courtship, as well as microsatellite markers and speciation genes. As many as five different chemotypes of sexual pheromones were characterized from Lu. longipalpis s.l. in Brazil . Three of the chemotypes are pheromones of the homosesquiterpene family, 3-methyl-α-himachalene (3MαH), (S)-9-methyl-germacrene-B (9MGB) and (S)-9-methyl-germacrene-B+ (9MGB+) (the latter two can be differentiated based on the amounts of 9MGB produced). Two other chemotypes are members of the diterpene family: cembrene-1 and cembrene-2. Interestingly, different chemotypes can be associated with different male copulation sounds or love songs generated by wing flapping: 3MαH is associated with the copulation songs Pulse songtype1; 9MGB is associated with either the Pulse songtype2 or the Pulse songtype3; cembrene-1 production is associated with Burstsong . No copulation songs have been associated with the production of pheromones 9MGB+ and cembrene-2 [3, 14, 15].
Lu. longipalpis s.l. in eastern Brazil is comprised of sibling species . However, studies have yet to demonstrate a link between these sibling species and geographic barriers that potentially participate in the speciation process of this sand fly. One such barrier is the Sao Francisco River, the second largest river in Brazil extending 3160 km. The Sao Francisco runs generally north behind the coastal range before turning east to form the border between the state of Bahia and the states of Pernambuco and Alagoas. It enters the Atlantic Ocean between the states of Alagoas and Sergipe (Figure 1A). Until the Mindel glaciation period (0.4 million years ago – MY) the Sao Francisco ran northward towards the State of Piaui (Figure 1B). During the Mindel glaciation, the course of the river was altered to its current location , as shown in Figure 1.
Natural geographical barriers such as rivers have been shown to significantly influence genetic diversification and speciation in distinct taxa [17–21]. Thus, we analyzed the population structuring of Lu. longipalpis s.l. in eastern Brazil based on different haplotypes of the cytochrome b gene (cyt b), and investigated the role of the Sao Francisco River in the speciation of this sand fly. Our results suggest that the natural barrier represented by the Sao Francisco River in Lu. longipalpis s.l. has played a significant role in the development of sibling species of this sand fly. The data are supported by evidence from recent studies on the geological history of the Sao Francisco River . To our knowledge, this work is the first to link the ongoing speciation of Lu. longipalpis s.l. in eastern Brazil with a significant natural barrier.
Lu. longipalpis s.l. used in this study were captured in five separate locations (Figure 1 and Table 1) using CDC light traps. Trapping of sand flies was done as described previously . Adult sand flies were killed by freezing at -20°C and identified as described .
DNA Extraction, Amplification and Sequencing
The DNA samples were prepared as reported previously , and amplification of cyt b was performed using the primers N1N-5'GGCAYWTTGCCTCGAWTTCGWTATGA3' and CB3-5'CAYATTCAACCWGAATGATA3' . PCR reactions were routinely carried out in 50 μl volume containing 10 pmoles of each primer, 1 ng of Lu. longipalpis s.l. DNA, 5 μl of 10× PCR buffer with 2 mM of MgCl2 (Biotools, B&M Labs, S.A. Madrid, Spain), 60 μM of dNTPs and 1.5 units of Taq DNA polymerase (Biotools). PCR reactions were performed as follows: a hot start at 94°C/3 min, followed by five cycles at 94°C/30 sec, 38°C/30 sec and 72°C/90 sec, and 35 cycles at 94°C/30 sec, 42°C/30 sec and 72°C/90 sec. The PCR products were purified with the Nucleospin Extract kit (Macherey-Nagel Inc., Easton, PA) and cycle-sequenced using the BigDye sequencing kit (Applied Biosystems, Foster City, CA), according to the manufacturer instructions. The sequencing reactions were performed using an ABI PRISM® 3100 sequence analyzer (Applied Biosystems). Each sample was sequenced twice and the sequencing quality was assessed using Phred  with a cutoff value of 25. The haplotypes were assembled using CAP3 .
Analysis of Genetic Variability
Analyses of genetic variability were performed in 34 haplotypes identified from Lu. longipalpis s.l. used in this work, including haplotypes described previously [26, 27]. The indexes of genetic variability were estimated using DnaSP 4.10 . The parameters measured were: the number of haplotypes (h); the haplotypic diversity (Hd); the number of segregating sites (S); the pairwise nucleotide diversity (π); the net nucleotide substitution (d A ); and the average number of nucleotide substitutions (d XY ). The θ parameter (2N e μ) was estimated for θπ and θ s . θπ is based on the average pairwise number of differences between sequences  and θ s is estimated from the number of segregation sites in a population . The following tests were performed to determine whether the pattern of polymorphisms were in accordance with the model of neutral evolution: D , D* and F* , and Fs . The Fs test was also used to assess population expansion.
The extent of nucleotide differentiation (F ST )  and the level of gene flow (Nm)  among the populations also were analyzed using DnaSP 4.10 . The F ST values were plotted in a matrix and used to create a distance-based Neighbor-Joining (NJ) tree using MEGA v3.1 . The Kimura-2P (K2P) genetics distance between clades of the NJ tree was also calculated using the MEGA v3.1 software. The genetic distances were used to calculate the time of divergence between clades, based on the rate of mutation of 2.3% sequence divergence per million years . The minimum spanning networks and the mismatch distributions were obtained using Arlequin v2.0 . Raggedness tests were performed to determine whether mismatch distributions fit the model of stable or expanding populations [39–41]. Also, we assessed the significance of the effect of geographic distance on differentiation via the Mantel test . This test was based on the linear regression between the value of the F ST /(1 - F ST )  and the geographical straight-line distance (in kilometers, Km) between each population pair. Geographic distances between capture sites were obtained from the Coordinate Distance Calculator . Results of the Mantel test were a reflection of 1,000 permutations using Mantel v1.18 .
The populations analyzed in this study were combined into two groups, taking into account the sand fly distribution in relation to the 10°S parallel. Thus, individuals that were originally captured north of the 10°S parallel were named <10°S; and those captured south of the 10°S parallel, >10°S
In this work we analyzed phylogeographic patterns of natural populations of Lu. longipalpis s.l.. in the Brazilian NE region. Our studies were based on polymorphisms of a 261 bp DNA fragment from the 3' end of the cyt b gene. The phylogenetic signal shown by this short fragment was higher than those presented by fragments of 318 and 489 nucleotides (unpublished data). Our analyses included twenty nine haplotypes previously described [26, 27] and five novel haplotypes identified from Lu. longipalpis s.l. from the Brazilian NE.
Out of 11 Lu. longipalpis s.l. individuals sequenced, nine distinct haplotypes were identified, including seven described earlier [26, 27]: Ll-8 from Calumbi; Ll-3 (present in two individuals) and Ll-9 from Itamaracá; Ll-15 from João Pessoa; and Ll-9 and Ll-10 from Patos (Table 1). Five novel haplotypes were named Ll-34 and Ll-35 (from Calumbi), Ll-36 (from João Pessoa), and Ll-37 (from Patos). A total of 34 haplotypes (including the ones described by [26, 27]) and their frequencies were analyzed as indicated above. The haplotypes differed up to six mutation steps and the polymorphic sites displayed no more than two variant nucleotides. The relationship between transition and transversion among the haplotypes was 25:1, with only five substitutions being non-synonymous (at position 26, a Leu to Ser change, L26S; then, S53N, P161L, V190M, N254S).
Our data shows that Ll-2 is the most frequent haplotype and the haplotypic diversity is high (0.947) despite the geographically widespread distribution of haplotypes Ll-1, Ll-2 and Ll-7 (Tables 2 and 3). Furthermore, no genetic differences were found between the populations of Juazeiro and Natal, Juazeiro and Patos, Natal and Patos, Natal and Calumbi, and Patos and Calumbi. In contrast, the highest genetic distance was estimated between populations of Itamaraca and Monte Santo (d A = 0.496%; Table 4).
A F ST distance-based Neighbor-Joining (NJ) tree was obtained for all populations investigated in this study. Two clades were identified: (I) including the populations of Natal, Patos, Juazeiro, João Pessoa, and Itamaracá; and (II) including the populations of Calumbi, Pancas, Monte Santo, and Feira de Santana (Figure 2). In light of the low level of genetic structuring and geographic distribution the two clusters seen in Figure 2 were referred to as group <10°S (for populations north of 10°S parallel) and group >10°S (for populations south of 10°S parallel). According to its geographic localization the population of Calumbi was also included in <10°S.
The net nucleotide substitution (d A ) between the groups <10°S and >10°S was 0.217% (Table 6). The estimates of haplotypic diversity within each group were high (Table 3), although some haplotypes (Ll-1, Ll-2, Ll-7 e Ll-8) presented a widespread geographic distribution. The average number of nucleotide differences (π) varied from 0.7% to 1.2% per site, or from 1.85 to 3.18 per nucleotide (Table 3), and no site were fixed at different nucleotides between groups. Within the 55 haplotypes from the group <10°S, 18 (32.7%) were singletons. In contrast, a lower proportion of singletons (20.6%, or 7 out of 34) was observed within the group >10°S.
Intra-population group analysis showed values of genetic structuring (F ST ) and gene flow (Nm) of 0.078 and 5.92 for the group <10°S (Table 6). Equivalent results were also obtained for the group >10°S (F ST = 0.136 and Nm = 3.16, Table 6). The values of F ST among the two groups support moderate genetic divergence and higher gene flow values (F ST = 0.184 and Nm = 2.22, Table 6), whereas the pattern of genetic structuring did not follow the model of isolation by distance (Figure 3).
To assess gene flow between the various populations investigated, we applied the minimum spanning network (MSN) analysis (Figure 4). The MSN constructed with all 34 different haplotypes (Figure 4A) shows that most of the haplotypes shared by groups <10°S and >10°S were positioned at the interior nodes, representing ancestral haplotypes [46, 47]. Thus, the Nm of 2.22 migrates per generation between groups <10°S and >10°S may reflect a greater influence from the shared ancestral polymorphisms than by ongoing gene flow, according to the outline of the MSN. Nevertheless, ongoing gene flow between these populations is still possible, as alternate connections between the terminal haplotypes are observed (Figure 4A).
MSN analysis was separately applied to the population groups <10°S and >10°S (Figure 4B and 4C). In both cases, interior shared haplotypes were also observed for each group. However, the existence of shared haplotypes at the tips of the networks, and the connections between terminal haplotypes of distinct populations points to ongoing gene flow in each case, consistent with panmitic populations. For the group <10°S (Figure 4B) the strong gene flow is supported by the shared haplotype positioned at the tip (Ll-9), and by several connections between terminal haplotypes of distinct populations (Figure 4B). For the group >10°S, despite the presence of shared ancestral haplotypes, terminal haplotypes from three populations were connected (Figure 4C), suggesting ongoing gene flow.
Neutrality tests  performed between the groups <10°S and >10°S resulted in negative values. Although the estimates of θs and θπ are higher than those of π, the D test  and the D* and F* tests  were not significant (Tables 3 and 7). On the other hand, the Fs test  showed significant values for the same population groups. The mismatch distribution for groups <10°S and >10°S was unimodal (not shown), with very small raggedness statistics (0.025 and 0.050, respectively).
The genetic distances (K2P) between groups was 1.2%, between groups <10°S and >10°S. Based on cyt b molecular clock, the divergence time between these groups were estimated to be 0.45 million years ago (MY).
The model of isolation by distance was tested through linear regression between the geographic distances and the values of F ST /(1-F ST ). We found that the pattern of genetic structuring of the populations was in disagreement with the model of isolation by distance (r = 0.347, P = 0.9628; Figure 3).
According to , there are five reproductively isolated populations of L. longipalpis s.l. in Brazil separated on the basis of chemotypes of the sexual pheromones they produce. The sexual pheromone 9MGB (chemotype 1) is produced by Lu. longipalpis s.l from Lapinha (MG); 3MαH (chemotype 2) is produced by flies from Jacobina (BA); cembrene-1 (chemotype 3) is found in flies from Sobral (CE), Santarem (PA), and Estrela de Alagoas and Costa do Sol (AL); cembrene-2 (chemotype 4) by flies from Jaiba (MG); finally, 9MGB+ (chemotype 5) by flies from Sobral (CE) and Montes Claros (MG).
Here, the genetic structuring of ten Lu. longipalpis s.l. populations from eastern Brazil was estimated through the analysis of 26 segregating sites of a partial segment of 261 nucleotides from the mitochondrial gene cyt b. The populations studied were clustered in two groups (as indicated in Materials and Methods), with individuals in group <10°S correspond to the sand flies that produce the chemotype 3 (cembrene-1 and 9MGB+); and those in the group >10°S to chemotypes 1, 2 and 5 (9MGB, 3MαH and 9MGB+).
Intra-population analyses of groups <10°S and >10°S revealed low and moderate genetic differentiation. However, in both cases, the values of Nm (Tables 4) were above the threshold required for differentiation by genetic drift (Nm < 2, ). Also, the presence of several tip connections between haplotypes of distinct populations, and a shared tip haplotype in the MSN for the group <10°S (Figure 4B) suggest that ongoing gene flow is likely responsible for the Nm value observed.
In contrast, inter-population analyses revealed moderate genetic differentiation and effective number of migrants between the groups <10°S and >10°S (Table 6). The Fs test and the smoothness of the mismatch distribution indicated that these groups are in geographic expansion. Since population expansion interferes with the estimates of gene flow [39–41], the Nm values obtained may be overestimated obscuring the actual population structuring. This, along with the prevalence of ancient introgression, rather than ongoing gene flow, indicates that genetic drift is mainly responsible for the moderate genetic differentiation between these groups.
Our results from the cyt b molecular clock suggest that the time of divergence between the Lu. longipalpis s.l. in groups <10°S and >10°S occurred around 0.45 MY. A similar divergence period (0.39 MY) was estimated between the populations of Baturite (located north the 10°S parallel) and Jacobina (located south of the same parallel). That estimation was based on the analysis of the molecular clock of the cytochrome oxidase II gene, according to a genetic distance of K2P = 0.009 , and a rate of divergence of 2.3% per MY . The estimated time of divergence between the two groups investigated here correspond to the Mindel glaciation (0.38–0.41 MY) that took place during the Early to Middle Pleistocene, and was associated with the shift in the Sao Francisco River's course .
The Sao Francisco River is the second largest river in Brazil, separating the States of Pernambuco (PE) and Ceara (CE) from Bahia (BA) and Minas Gerais (MG). On its original course, the river flowed from MG, bearing northward through BA and continuing towards the coast of the State of Piaui (PI) through what is today the Piaui River and the Parnaiba River Basins. Following the uplift of the Parnaiba River Basin and the formation of a transversal geological fault to Northeastern coast, the course of the Sao Francisco River was altered to its current location , as shown in Figure 1A.
The reduced gene flow between sibling species and moderate genetic differentiation between groups <10°S and >10°S are consistent with a distribution pattern likely driven by the Sao Francisco River. Additionally, the estimated time of the divergence of these two groups coincides with the change of the river course. However, this divergence cannot be associated with reduced gene flow as mitochondrial DNA is prone to introgress through incipient species boundaries , and as demonstrated via cyt b gene analysis of Lu. longipalpis s.l. sympatric cryptic species . Taken together, the data suggest that the Sao Francisco River is a significant geographic barrier between populations of Lu. longipalpis s.l., and could also have contributed to the current level of population diversity seen for this sand fly.
This putative role of the Sao Francisco River in the speciation of Lu. longipalpis s.l. may also explain the current distribution of male sexual pheromones chemotypes. A plausible scenario is based on the suggestion that 9MGB represents an ancestral chemotype pheromone and that the diterpenes (i.e., cembrene-1 and 2) are more recent chemotypes . Dispersal of Lu. longipalpis s.l. belonging to the Brazilian clade  is estimated to have occurred in the Plioceno-Pleistocene  when the Sao Francisco River flowed on a South-to-North direction . Thus, the Sao Francisco River would have served as a barrier to the introduction of the 9MGB chemotype in the Brazilian NE coast.
The current chemotype distribution in the Brazilian NE possibly followed a process similar to what is described for the ring species hypothesis . Under this hypothesis Lu. longipalpis s.l. would have reached the coastal areas of the NE by circumventing the River's nascent, along its west margin, in Minas Gerais, South Eastern Brazil (as shown in Figure 1B). This process led to the adaptation to different habitats, and to differentiation into a population producing the cembrane-1 chemotype presently found in the NE of Brazil. This idea is supported by the presence of Lu. longipalpis s.l. in the Brazilian SE producing different pheromone chemotypes, such as 9MGB found in Lapinha Cave, cembrene-2 in found in Jaiba, and 3MαH found in Jacobina [14, 54]. This scenario also supports the suggestion that Lu. longipalpis s.l. sibling species currently found in Sobral were the result of vicariance 1MY ago , and prior to the change in the course of the Sao Francisco River. Thus, the sympatric species in Sobral were isolated by the Sao Francisco River, leading to the development of 9MGB+ on the west and cembrene-1 on the east margins. Following the shift of the Sao Francisco River these species reached the current secondary contact proposed by  (Figure 1A). Our results suggesting greater similarity between Lu. longipalpis s.l. from Juazeiro (located on the southern rim of the Sao Francisco River) and flies belonging to group <10°S in comparison to flies belonging to group >10°S can be explained by a recent break-down of this geographic barrier, either by a recent geographic accident or through human activity in the region. This is consistent with our data indicating that no ongoing mitochondrial DNA haplotype exchange is detected between individuals from Juazeiro and individuals from group >10°S (Figure 4A).
Lu. longipalpis s.l. secreting cembrene-1 also are present in areas other than the Brazilian NE [14, 54, 56]. In two such areas (Marajo Island and Santarem), the presence of cembrene-1 was possibly due to a recent geographic expansion after the change in course of the Sao Francisco River. A study using microsatellite analysis indicated strong similarities between Lu. longipalpis s.l. from Marajo Island and flies from either Sobral or Natal . Accordingly, our data supports a pattern compatible with population expansion from the latter two locations towards the former based on the mismatch distribution. In regards to the cembrene-1 population currently found in Espirito Santo do Pinhal, Sao Paulo , biogeographic and genetic studies are necessary in order to explain its origin.
Based on the data presented here we propose that the Sao Francisco River restricted gene flow between Lu. longipalpis s.l., participating in speciation processes of this sand fly. This paper contributes to our understanding on the expansion of Lu. longipalpis s.l. in Brazil and provides novel clues in regards to several aspects of the divergence of this important vector species.
Ashford RW: Leishmaniasis as emerging and reemerging zoonoses. Int J Parasitol. 2000, 30: 1269-1281. 10.1016/S0020-7519(00)00136-3.
Lainson R, Rangel EF: Lutzomyia longipalpis and the eco-epidemiology of American visceral leishmaniasis, with particular reference to Brazil: a review. Mem Inst Oswaldo Cruz. 2005, 100 (8): 811-827. 10.1590/S0074-02762005000800001.
Bauzer LG, Souza NA, Maingon RD, Peixoto AA: Lutzomyia longipalpis in Brazil: a complex or a single species? A mini-review. Mem Inst Oswaldo Cruz. 2007, 102 (1): 1-12. 10.1590/S0074-02762007000100001.
Young D, Duncan M: Guide for identification and geographic distribution of Lutzomyia sand-flies in Mexico, the West Indies, Central and South America (Diptera: Psychodidae). Mem Am Entomol Inst. 1994, 54: 1-881.
Black IV, Mustermann L: Molecular taxonomy and systematics of arthropod vectors. Biology of Disease Vectors. Edited by: WC M. 2005, Amsterdam , Elsevier Academic Press, 297-310. second
Arrivillaga J, Mutebi JP, Pinango H, Norris D, Alexander B, Feliciangeli MD, Lanzaro GC: The taxonomic status of genetically divergent populations of Lutzomyia longipalpis (Diptera: Psychodidae) based on the distribution of mitochondrial and isozyme variation. J Med Entomol. 2003, 40 (5): 615-627.
Arrivillaga JC, Norris DE, Feliciangeli MD, Lanzaro GC: Phylogeography of the neotropical sand fly Lutzomyia longipalpis inferred from mitochondrial DNA sequences. Infect Genet Evol. 2002, 2 (2): 83-95. 10.1016/S1567-1348(02)00087-4.
Lanzaro GC, Ostrovska K, Herrero MV, Lawyer PG, Warburg A: Lutzomyia longipalpis is a species complex: genetic divergence and interspecific hybrid sterility among three populations. Am J Trop Med Hyg. 1993, 48 (6): 839-847.
Soto SI, Lehmann T, Rowton ED, Velez BI, Porter CH: Speciation and population structure in the morphospecies Lutzomyia longipalpis (Lutz & Neiva) as derived from the mitochondrial ND4 gene. Mol Phylog Evol. 2001, 18 (1): 84-93. 10.1006/mpev.2000.0863.
Yin H, Mutebi JP, Marriott S, Lanzaro GC: Metaphase karyotypes and G-banding in sandflies of the Lutzomyia longipalpis complex. Med Vet Entomol. 1999, 13 (1): 72-77. 10.1046/j.1365-2915.1999.00139.x.
Yin H, Norris DE, Lanzaro GC: Sibling species in the Llutzomyia longipalpis complex differ in levels of mRNA expression for the salivary peptide, maxadilan. Insect Mol Biol. 2000, 9 (3): 309-314. 10.1046/j.1365-2583.2000.00190.x.
Bauzer LG, Souza NA, Ward RD, Kyriacou CP, Peixoto AA: The period gene and genetic differentiation between three Brazilian populations of Lutzomyia longipalpis. Insect Mol Biol. 2002, 11 (4): 315-323. 10.1046/j.1365-2583.2002.00340.x.
Maingon RD, Ward RD, Hamilton JG, Noyes HA, Souza N, Kemp SJ, Watts PC: Genetic identification of two sibling species of Lutzomyia longipalpis (Diptera: Psychodidae) that produce distinct male sex pheromones in Sobral, Ceara State, Brazil. Mol Ecol. 2003, 12 (7): 1879-1894. 10.1046/j.1365-294X.2003.01871.x.
Hamilton JG, Maingon RD, Alexander B, Ward RD, Brazil RP: Analysis of the sex pheromone extract of individual male Lutzomyia longipalpis sandflies from six regions in Brazil. Med Vet Entomol. 2005, 19 (4): 480-488. 10.1111/j.1365-2915.2005.00594.x.
Maingon RD, Ward RD, Hamilton JG, Bauzer LG, Peixoto AA: The Lutzomyia longipalpis species complex: does population sub-structure matter to Leishmania transmission?. Trends Parasitol. 2008, 24 (1): 12-17. 10.1016/j.pt.2007.10.003.
Mabesoone JM: Sedimentary basins of northeast Brazil. 1994, Recife , Federal University of Pernambuco, 310-
Hall JPW, Harvey DJ: The phylogeography of Amazonia revisited: New evidence from riodinid butterflies. Evolution. 2002, 35: 829-843.
Rodriguez-Robles JA, Jezkova T, Leal M: Genetic structuring in the threatened ‘‘Lagartijo del Bosque Seco” (Anolis cooki) from Puerto Rico. Molecular Phylogenetics and Evolution. 2008, 46: 503-514. 10.1016/j.ympev.2007.11.006.
Wilson GM, Den Brussche RA, McBee K, Jonhson LA, Jones CA: Intraspecific phylogeography of red squirrels (Tamiasciurus hudsonicus) in the central rocky mountain region of North America. Genetica. 2005, 125: 141-154. 10.1007/s10709-005-5154-5.
Wood DA, Fisher AN, Reeder TW: Novel patterns of historical isolation, dispersal, and secondary contact across Baja California in the Rosy Boa (Lichanura trivirgata). Molecular Phylogenetics and Evolution. 2008, 46: 484-502. 10.1016/j.ympev.2007.11.014.
Miller MP, Haig SM, Wagner RS: Phylogeography and Spatial Genetic Structure of the Southern Torrent Salamander: Implications for Conservation and Management. Journal of Heredity. 2006, 97 (6): 561-570. 10.1093/jhered/esl038.
de Queiroz Balbino V, 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.
Ready PD, Day JC, Souza AA, Rangel EF, Davies CR: Mitochondrial DNA characterization of populations of Lutzomyia whitmani (Diptera: Psychodidae) incriminated in the peri-domestic and silvatic transmission of Leishmania species in Brazil. Bull Entomol Res. 1997, 87: 187-195.
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 (3): 211-225. 10.1046/j.1365-2583.1997.00175.x.
Huang X, Madan A: CAP3: A DNA sequence assembly program. Genome Res. 1999, 9 (9): 868-877. 10.1101/gr.9.9.868.
Hodgkinson VH, Birungi J, Haghpanah M, Joshi S, Munstermann LE: Rapid identification of mitochondrial cytochrome B haplotypes by single strand conformation polymorphism in Lutzomyia longipalpis (Diptera: Psychodidae) populations. J Med Entomol. 2002, 39 (4): 689-694.
Hodgkinson VH, Birungi J, Quintana M, Dietze R, Munstermann LE: Mitochondrial cytochrome b variation in populations of the visceral leishmaniasis vector Lutzomyia longipalpis across eastern Brazil. Am J Trop Med Hyg. 2003, 69 (4): 386-392.
Rozas J, Sanchez-DelBarrio JC, Messeguer X, Rozas R: DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics. 2003, 19 (18): 2496-2497. 10.1093/bioinformatics/btg359.
Tajima F: Evolutionary relationship of DNA sequences in finite populations. Genetics. 1983, 105 (2): 437-460.
Watterson GA: On the number of segregating sites in genetical models without recombination. Theor Popul Biol. 1975, 7 (2): 256-276. 10.1016/0040-5809(75)90020-9.
Tajima F: Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989, 123 (3): 585-595.
Fu YX, Li WH: Statistical tests of neutrality of mutations. Genetics. 1993, 133 (3): 693-709.
Fu YX: Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics. 1997, 147 (2): 915-925.
Hudson RR, Boos DD, Kaplan NL: A statistical test for detecting geographic subdivision. Mol Biol Evol. 1992, 9 (1): 138-151.
Hudson RR, Slatkin M, Maddison WP: Estimation of levels of gene flow from DNA sequence data. Genetics. 1992, 132 (2): 583-589.
Kumar S, Tamura K, Nei M: MEGA3: Integrated software for Molecular Evolutionary Genetics Analysis and sequence alignment. Brief Bioinform. 2004, 5 (2): 150-163. 10.1093/bib/5.2.150.
Brower AV: Rapid morphological radiation and convergence among races of the butterfly Heliconius erato inferred from patterns of mitochondrial DNA evolution. Proc Natl Acad Sci U S A. 1994 , 91 (14): 6491-6495. 10.1073/pnas.91.14.6491.
Slatkin M, Hudson RR: Pairwise comparisons of mitochondrial DNA sequences in stable and exponentially growing populations. Genetics. 1991, 129 (2): 555-562.
Rogers AR, Harpending H: Population growth makes waves in the distribution of pairwise genetic differences. Mol Biol Evol. 1992, 9 (3): 552-569.
Harpending HC: Signature of ancient population growth in a low-resolution mitochondrial DNA mismatch distribution. Hum Biol. 1994, 66 (4): 591-600.
Mantel N: The detection of disease clustering and a generalized regression approach. Cancer Res. 1967, 27 (2): 209-220.
Rousset F: Genetic differentiation and estimation of gene flow from F-statistics under isolation by distance. Genetics. 1997, 145 (4): 1219-1228.
Mantel-v1.18: . [http://life.bio.sunysb.edu/morph/]
Bandelt HJ, Forster P, Röhl A: Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999, 16: 37-48.
Donnelly MJ, Pinto J, Girod R, Besansky NJ, Lehmann T: Revisiting the role of introgression vs shared ancestral polymorphisms as key processes shaping genetic diversity in the recently separated sibling species of the Anopheles gambiae complex. Heredity. 2004, 92 (2): 61-68. 10.1038/sj.hdy.6800377.
Nei M: Molecular evolutionary genetics. 1987, New York , Columbia University press, 512-
Slatkin M: Gene flow and the geographic structure of natural populations. Science. 1987, 236 (4803): 787-792. 10.1126/science.3576198.
Coyne JA, Orr HA: Speciation. 2004, Sinauer, 545-
Ward RD, Ready P, Watts PC: A molecular analysis of the Lutzomyia longipalpis species complex. 2001, Abstracts Book: 127-
Conn JE, Mirabello L: The biogeography and population genetics of neotropical vector species. Heredity. 2007, 99 (3): 245-256. 10.1038/sj.hdy.6801002.
Mayr E: Systematics and the Origin of Species. 1942, New York , Columbia University Press
Casanova C, Hamilton JGC, Trigo JR, Costa AIP: Identification of Sex Pheromones of Lutzomyia longipalpis (Lutz & Neiva, 1912) populations from the State of Sao Paulo, Brazil. Mem Inst Oswaldo Cruz. 2006, 101 (1): 113-115. 10.1590/S0074-02762006000100023.
Bauzer LG, Gesto JS, Souza NA, Ward RD, Hamilton JG, Kyriacou CP, Peixoto AA: Molecular divergence in the period gene between two putative sympatric species of the Lutzomyia longipalpis complex. Mol Biol Evol. 2002, 19 (9): 1624-1627.
Watts PC, Hamilton JG, Ward RD, Noyes HA, Souza NA, Kemp SJ, Feliciangeli MD, Brazil R, Maingon RD: Male sex pheromones and the phylogeographic structure of the Lutzomyia longipalpis species complex (Diptera: Psychodidae) from Brazil and Venezuela. Am J Trop Med Hyg. 2005, 73 (4): 734-743.
Hamilton JG, Brazil RP, Maingon R: A fourth chemotype of Lutzomyia longipalpis (Diptera: Psychodidae) from Jaibas, Minas Gerais State, Brazil. J Med Entomol. 2004, 41 (6): 1021-1026.
We are grateful to Alexandre Magalhães Martins (PPGG-UFPE) for technical assistance.
The authors declare that they have no competing interests.
IVCA participated in the conception of the study, analyses of the data, and drafting of the manuscript; IVS participated in data analyses and design of the study; JAF and MAM participated in sample collection; VQB and MRO participated in the conception and design of the study and crafting of the manuscript.
Iliano V Coutinho-Abreu, Ivan V Sonoda, Valdir Q Balbino and Marcelo Ramalho-Ortigão contributed equally to this work.