Genetic structuring and fixed polymorphisms in the gene period among natural populations of Lutzomyia longipalpis in Brazil
- César Raimundo Lima CostaJr1,
- Moises Thiago de Souza Freitas1,
- Carlos Alberto Santiago FigueirêdoJr1,
- Nádia Consuelo Aragão1,
- Lidiane Gomes da Silva1,
- Carlos Brisola Marcondes2,
- Raimundo Vieira Dias3,
- Tereza Cristina Leal-Balbino4,
- Manuela Barbosa Rodrigues Souza1,
- Marcelo Ramalho-Ortigão5 and
- Valdir de Queiroz Balbino1Email author
© Lima Costa et al.; licensee BioMed Central. 2015
Received: 31 October 2014
Accepted: 5 March 2015
Published: 1 April 2015
Even one hundred years after being originally identified, aspects of the taxonomy of the sand fly Lutzomyia longipalpis, the principal vector of Leishmania infantum in the Americas, remain unresolved for Brazilian populations of this vector. The diversity of morphological, behavioral, biochemical, and ethological characters, as well as the genetic variability detected by molecular markers are indicative of the presence of a complex of species.
In this study, a 525 bp fragment of the period gene was used to evaluate sympatric populations of L. longipalpis. A combination of probabilistic methods such as maximum likelihood and genetic assignment approach to investigate sympatric species of L. longipalpis were applied in three populations of Northeast Brazil.
Fixed polymorphisms in geographically isolated populations of L. longipalpis from two localities in the state of Ceará and one in the state of Pernambuco, Brazil, was identified in a 525 bp fragment of the gene period (per). Our results suggest a direct relationship between the number of spots found in males’ tergites and the genetic variation in cryptic species of L. longipalpis. The fragment used in this study revealed the nature of the ancestral morphotype 1S.
New polymorphisms were identified in the gene per which can be used as a genetic barcode to sympatric taxonomy of L. longipalpis. The per gene fragment confirmed the presence of two siblings species of L. longipalpis in Sobral and showed that these same species are present in two other localities, representing an expansion within the L. longipalpis species complex with regards to the states of Ceará and Pernambuco.
A species complex is generally defined as a group of morphologically similar species that differ with regards to genetic and ethological aspects . Low flight range and geographic isolation between populations are major drivers of allopatric speciation, contributing to constitution of cryptic species in Phlebotomine sand flies [2-4]. The occurrence of cryptic species is a fairly common event in other disease vectors and have been described in Anopheles gambiae, Culex pipiens, Triatoma dimidiata and phlebotomines [5-8].
Lutzomyia longipalpis, the principal vector of Leishmania infantum, the etiologic agent of American visceral leishmaniasis, has a discontinued distribution throughout the Neotropical region with different populations displaying aspects compatible with those of a species complex [9,10]. The existence of cryptic species in the Lutzomyia longipalpis complex was supported by studies using morphological and molecular markers in populations from Central and South America leading to the identification of L. pseudolongipalpis , formally recognized as the first taxon of the L. longipalpis complex (reviewed by ). Later, it was also suggested that L. cruzi should be regarded as a cryptic species within the L. longipalpis complex .
In Brazil, the presence of a Lutzomyia longipalpis complex was first proposed by Mangabeira  and was based on the number of abdominal pale spots, one (1S) or two (2S), visible on abdominal tergites. Analyses of male sexual pheromones and courtship (mating) sounds, as well as microsatellite markers and speciation genes, have provided further evidence of a L. longipalpis species complex .
The gene period (per) controls biological rhythms and plays a central role in eclosion and insect locomotion activity [15,16]. Per has been used in studies of fruit flies (Drosophilidae) population genetics and has been shown to be a useful marker, especially when it comes to identifying cryptic species [17,18]. In sand flies, per has also been utilized in studies of population genetics to identify possible members of the L. longipalpis complex in Brazil. Combined with certain behavioral markers (e.g. male mating songs and the types of sexual pheromones produced), analysis of the variability of per suggested the existence of two major population groups of L. longipalpis in Brazil. One group, with at least five species, is found in most localities where L. longipalpis has been reported. This group displays a single spot (1S) on the abdominal tergites, but different pheromones and patterns of pulsating songs. The second group, with two spots (2S), is represented by a single species present in the North, Northeast and Southeast regions of Brazil. Males in this group produce Burst-type songs and the pheromone cembrene-1 .
The NE region in Brazil is the region with the highest index of visceral leishmaniasis , and studies demonstrated the presence of at least two sister species within L. longipalpis complex [19,21,22]. In Sobral, State of Ceará, these species occur in sympatry and can be separated based either on their abdominal spots 1S or 2S, the different mating songs, the types of pheromones, or the genetic composition of males . Hence, genetic markers related to the abdominal spots may indeed be an important tool for molecular taxonomy of the L. longipalpis complex, especially concerning eco-epidemiological studies to assess potential vectorial capacity between distinct populations.
Here, we investigated the presence of polymorphisms in the gene per and their relatedness to the number of abdominal spots in L. longipalpis males. Sand flies from three localities in the NE of Brazil, separated by distances ranging from 108 km to as much as 457 km from each other, were used for our analyses. Four novel fixed single nucleotide polymorphisms (SNP) were identified, some able to separate between the 1S and 2S flies. In addition, our results point to an ancestral origin of the morphotype 1S. This study contributes to the understanding of the natural history of L. longipalpis populations and provides new insights about the biogeography of this sand fly.
Field collection and identification of phebotomine sand flies
Sand fly trappings were carried out in Sobral (3°41′15″S; 40°21′5″W) and Caririaçu (07°02′31″S; 39°17′02″W) in the State of Ceará, and in Bodocó (07°46′42″S, 39°56′28″W) in the State of Pernambuco, Brazil. All locations included in this study have a BSw‘h’ climate in accordance to the Köppen climate classification , with temperatures ranging from 23°C to 36°C and low annual rainfall (936 mm to 1100 mm). The localities included in this study are fully inserted in Caatinga biome, with vegetation composed mainly of ligneous and herbaceous species with high degree of xerophily .
Sand flies were trapped in the surrounding houses and domestic animal shelters using five CDC-type miniature light traps positioned at approximately 0.6 m from the ground. Sand flies were identified according to Young and Duncan , and L. longipalpis males were separated based on the number of abdominal spots into 1S and 2S.
DNA extraction, PCR and sequencing
Genomic DNA was extracted from each specimen using 100 μl of Chelex® resin (Bio-Rad, Hercules, CA) based on the protocol described by Solano et al.  with modifications. Briefly, each sand fly was homogenized using a hand held homogenizer and pestle in 10% Chelex solution followed by incubation at 56°C for an hour. Sample lysates were incubated at 95°C for 30 min, centrifuged for 6 min at 13,000 × g, and each supernatant was removed and stored at −4°C.
For each DNA sample isolated, a segment of 525 bp of per  was amplified by PCR using the Master Mix kit (Promega, Fitchburg, WI). PCRs were carried out using a PTC-200 thermocycler (MJ Research, Ramsey, MN) as follows: 5 min at 95°C, and 30 cycles of 95°C for 30 seconds, 50°C for 30 seconds, and 72°C for 1 min, with a final extension at 72°C for 10 min. A 5 μl aliquot from each PCR product was separated on 0.5% agarose and the remainder was purified using the Genomic DNA Purification Wizard kit (Promega, Fitchburg).
Bi-directional sequencing reactions were performed on each purified PCR product using the BygDye Terminator v3.1 Matrix Standard (Applied Biosystems, Foster City, CA) and analyzed using an ABI3100 Sequence Analyzer (Applied Biosystems, Foster City). Each sample was sequenced in duplicate and the sequences obtained were assembled and analyzed with the Staden package  based on the values of Phred 40 . The high-quality sequences obtained were deposited in GenBank (accession numbers KF479047-KF479163 and KP013849-KP013902).
Sequence alignment was performed using the program ClustalW (MEGA v.5.1), and the conserved and variables sites (parsimony informative and singletons sites) also were verified using MEGA v.5.1 . The analyses of polymorphisms and the neutrality tests were performed using DnaSP 5.10.1 . For optimal viewing of polymorphisms, parsimony-informative sites were exported to compose sequence logos obtained by Weblogo v. 3.2 .
The genetic structuring of the L. longipalpis populations was verified using sequences alignment throughout STRUCTURE v.2.3 , with model-based approach. This algorithm clusters individuals into populations. The assigned proportion of each individual belonging to each population (membership coefficient Q) was estimated using Bayesian statistics and Markov Chain Monte Carlo simulations (MCMC). MCMC simulations were performed with 100,000 interactions of burn-in period and followed by one million steps. For each value of K (1 to 10), 10 interactions were performed in order to estimate the K values, and the most likely population (or cluster) number was determined by the ΔK analysis . Matched Fst and MN values  were obtained through the Arlequin v.3.1 software , using 1000 random permutations. Population structure was verified by the utilization of AMOVA , based on the evaluation of different hierarchical groups.
For the UPGMA tree, the Fst values were imported into MEGA V. 5.1 software. The jModeltest v.0.1.1 software was used to select the best-fit DNA substitution model for ML analysis based on the Akaike information criterion (AIC) algorithm . The best-fit model for the per sequence analyzed was the TVM + I + G. The maximum likelyhood (ML) tree for per sequences was calculated using PhyML v. 3.0 .
Neutrality tests and intra-population genetic diversity measures for each sample
π ± SE
0.00684 ± 0.00094
0.00885 ± 0.00089
0.01392 ± 0.00108
0.01058 ± 0.00056
0.01648 ± 0.00108
0.01162 ± 0.00089
0.01901 ± 0.01901
Genetic differentiation among samples
AMOVA results for L. longipalpis populations
Source of variation
Percentage of variation
Among populations within groups
FSC (haplotypes/populations within groups)
The composition of the L. longipalpis complex is still a rather controversial topic even though the distribution of this sand fly in Brazil strongly suggests the presence of at least five species as part of this complex [10,19].
Our assessment of polymorphism frequency in the per gene in different populations of L. longipalpis revealed for the first time the presence of two fixed polymorphisms (124 and 424) that can reliably be used to separate 1S and 2S phenotypes. The SNP at position T124C was previously reported in L. longipalpis populations collected in the town of Jaíba, Minas Gerais, located 1,347 km from Sobral and 1,040 km from Caririaçu . The SNP at position C424T was identified in L. longipalpis populations from Bodocó and Caririaçu. Two additional SNPs, at positions C171T and T256C can be used as coadjuvants to separate between the 1S and 2S L. longipalpis morphotypes from Sobral and Bodocó, respectively. Although all polymorphisms are synonymous, the combined SNPs can be used as markers for cryptic species in L. longipalpis.
The importance of the ML method for the detection of recent speciation events also was previously shown [40-42]. The Bayesian and ML approaches for L. longipalpis collected in Bodocó, Caririaçu and Sobral indicate the presence of two genetically separated populations, and that each collected sample belongs to either 1S or 2S morphotype. Our analyses also indicated that these two genetically distinct populations are related to abdominal spots, as previously suggested for sympatric populations of L. longipalpis [21-23].
Analysis the 525 bp fragment of per in sand flies from Sobral showed greater separation between the two morphotypes, even when genetic distance methodology was used (Additional file 1: Figure S1). A Neighbor-joining tree added to the population genetic structure analyses confirmed the presence of cryptic species occurring in sympatry in Sobral, as noted by Bauzer et al. . In that case, the presence of a haplotype shared by distinct morphotypes was observed. In contrast, our analyses using a greater number of polymorphisms revealed patterns of total association between phenotypes and genotypes [22,43].
As the number of shared haplotypes observed was lower in the 525 bp fragment of per in comparison to the 266 bp fragment (15 shared haplotypes versus 26, Additional file 2: Table S1), it was suggestive of the ability of the 525 bp per to detect greater genetic variation in L. longipalpis populations. Analyses by AMOVA of the results obtained with the 525 bp fragment confirmed a greater variation between phenotypically distinct groups characterized by the association of a pattern of abdominal spots and the genetic marker not previously reported for L. longipalpis populations.
In Sobral, L. longipalpis is described as distinct populations of two sympatric species, and commonly separated by the 1S and 2S phenotypes. Genetic markers, copulatory sounds, and sex pheromones have been characterized for each of the L. longipalpis populations found in Sobral and have been used to confirm the separation conferred by the 1S/2S phenotypes [19,23,44,45]. Conversely, L. longipalpis collected in the localities of Bodocó and Caririaçu are, according to our data, two genetically distinct populations, and the results observed, match the pattern of abdominal spots described in the sympatric populations [22,23,43]. Thus, as L. longipalpis from Bodocó and Caririaçu exhibit the same phenotypic and genetic patterns of populations found in Sobral, it is likely that these two populations also share chemical and behavioral characteristics similar to what was previously described for the L. longipalpis found in Sobral.
Our group previously reported on a secondary contact between L. longipalpis populations that were separated by the original course of the São Francisco River . This secondary contact between the 1S populations of the Brazilian SE and NE may have promoted the genetic diversity of 1S. It also further reinforced the hypothesis that 2S actually derives from 1S, in accordance to the maximum likelihood tree.
Understanding the complex population genetics of L. longipalpis and its pattern of distribution is critical in areas with high endemicity for the transmission of visceral leishmaniasis, such as the current situation in the state of Ceará in Brazil. The genetic analyses using the 525 bp fragment of the per gene revealed for the first time a moderate geographical structuring between the L. longipalpis populations, and a significant variability with regards to the 1S and 2S phenotypes. The results presented here also underscore the importance of the abdominal spots for the diagnosis of cryptic species of sympatric populations of the L. longipalpis complex in Brazil, and the use of per as an important barcode marker to populations of L. longipalpis. Further confirmation of the fixed SNPs identified in per and its application as taxonomic markers to differentiate sympatric L. longipalpis populations across Brazil is warranted.
Ethical approval was not required for the current study.
This study was supported by Fundação de Amparo à Ciência e Tecnologia do Estado de Pernambuco (APQ-0651-2.02/10), and Ministério da Saúde, Brazil (TC 300/2013). VQB was supported by CNPq, Grant Number 309124/2012-3. We are also grateful to the Secretaria de Saúde do Município de Sobral for assistance and logistic support.
- Bickford D, Lohman DJ, Sodhi NS, Ng PKL, Meier R, Winker K, et al. Cryptic species as a window on diversity and conservation. Trends Ecol Evol. 2007;22:148–55. doi:10.1016/j.tree.2006.11.004.View ArticlePubMedGoogle Scholar
- Alexander JB. Dispersal of phlebotomine sand flies (Diptera: Psychodidae) in a Colombian coffee plantation. J Med Entomol. 1987;24:552–58.View ArticlePubMedGoogle Scholar
- Morrison AC, Ferro C, Tesh RB. Host preferences of the sand fly Lutzomyia longipalpis at an endemic focus of american visceral leishmaniasis in Colombia. Am J Trop Med Hyg. 1993;49:68–75.PubMedGoogle Scholar
- Cohnstaedt LW, Beati L, Caceres AG, Ferro C, Munstermann LE. Phylogenetics of the phlebotomine sand fly group verrucarum (Diptera: Psychodidae: Lutzomyia). Am J Trop Med Hyg. 2011;84:913–22.View ArticlePubMed CentralPubMedGoogle Scholar
- Anderson TJC, Haubold B, Williams JT, Estrada-Franco JG, Richardson L, Mollinedo R, et al. Microsatellite markers reveal a spectrum of population structures in the malaria parasite Plasmodium falciparum. Mol Biol Evol. 2000;17:1467–82.View ArticlePubMedGoogle Scholar
- Krzywinski J, Besansky NJ. Molecular systematics of Anopheles: from subgenera to subpopulations. Annu Rev Entomol. 2003;48:111–39.View ArticlePubMedGoogle Scholar
- Adler PH, Cheke RA, Post RJ. Evolution, epidemiology, and population genetics of black flies (Diptera: Simuliidae). Infect Genetic Evol. 2010;10:846–65. doi:10.1016/j.meegid.2010.07.003.View ArticleGoogle Scholar
- Monteiro F, Peretolchina T, Lazoski C, Harris K, Dotson E. Phylogeographic pattern and extensive mitochondrial DNA divergence disclose a species complex within the Chagas disease vector Triatoma dimidiata. PLoS One. 2013;8:e70974.View ArticlePubMed CentralPubMedGoogle Scholar
- Lainson R, Rangel EF. Lutzomyia longipalpis and the eco-epidemiology of American visceral leishmaniasis, with particular reference to Brazil. Mem Inst Oswaldo Cruz. 2005;100:811–27.View ArticlePubMedGoogle Scholar
- Bauzer LGSR, Souza NA, Maigon RDC, Peixoto AA. Lutzomyia longipalpisin Brazil: a complex or a single species? A mini-review. Mem Inst Oswaldo Cruz. 2007;102:1–12.View ArticlePubMedGoogle Scholar
- Arrivillaga JC, Feliciangeli MD. Lutzomyia pseudolongipalpis: the first new species within the longipalpis (Diptera: Psychodidae: Phlebotominae) complex from La Rinconada, Curarigua, Lara State, Venezuela. J Med Entomol. 2001;l38:783–90.View ArticleGoogle Scholar
- Soares RP, Turco SJ. Lutzomyia longipalpis (Diptera: Psychodidae: Phlebotominae): a review. An Acad Bras Cienc. 2003;75:301–30.View ArticlePubMedGoogle Scholar
- Vigoder FM, Araki AS, Bauzer LGSR, Souza NA, Brazil RP, Peixoto AA. Lovesongs and period gene polymorphisms indicate Lutzomyia cruzi (Mangabeira, 1938) as a sibling species of the Lutzomyia longipalpis (Lutz and Neiva, 1912) complex. Infec Gen Evol. 2010;10:734–39. doi:10.1016/j.meegid.2010.05.004.View ArticleGoogle Scholar
- Mangabeira O. Sobre a sistemática e biologia dos Phlebotomus do Ceará. Rev Bras Malar Doenç Trop. 1969;21:3–26.Google Scholar
- Hardin PE, Hall JC, Rosbash M. Feedback of the Drosophila period gene product on circadian cycling of its messenger RNA levels. Nature. 1990;343(6258):536–40.View ArticlePubMedGoogle Scholar
- Dunlap JC. Molecular bases for circadian clocks. Cell. 1999;96(2):271–90. doi:10.1016/S0092-8674(00)80566-8.View ArticlePubMedGoogle Scholar
- Rosato E, Peixoto AA, Barbujani G, Costa R, Kyriacou CP. Molecular evolution of the period gene in Drosophila simulans. Genetics. 1994;138:693–707.PubMed CentralPubMedGoogle Scholar
- Sawyer LA, Sandrelli F, Pasetto C, Peixoto AA, Rosato E, Costa R, et al. The period gene Thr-Gly polymorphism in Australian and African Drosophila melanogaster populations: implications for selection. Genetics. 2006;174:465–80.View ArticlePubMed CentralPubMedGoogle Scholar
- Araki AS, Vigoder FM, Bauzer LGSR, Ferreira GEM, Souza NA, Araújo IB, et al. Molecular and behavioral differentiation among brazilian populations of Lutzomyia longipalpis (Diptera: Psychodidae: Phlebotominae). PLoS Negl Trop Dis. 2009;3:1–12.View ArticleGoogle Scholar
- Ministério da Saúde. 2013. http://portalsaude.saude.gov.br/images/pdf/2014/setembro/09/LV-Casos.pdf.
- Lins RM, Souza NA, Brazil RP, Maingon RD, Peixoto AA. Fixed differences in the paralytic gene define two lineages within the Lutzomyia longipalpis complex producing different types of courtship songs. PLoS One. 2012;7:e44323.View ArticlePubMed CentralPubMedGoogle Scholar
- Silva MH, Nascimento MDSB, Leonardo FS, Rebêlo JMM, Pereira SRF. Genetic differentiation in natural populations of Lutzomyia longipalpis (Lutz & Neiva) (Diptera: Psychodidae) with different phenotypic spot patterns on tergites in males. Neotrop Entomol. 2011;40:501–6.View ArticlePubMedGoogle Scholar
- 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:315–23.View ArticlePubMedGoogle Scholar
- Peel MC, Finlayson BL, McMahon TA. Updated world map of the Koppen Geiger climate classification. Hydrol Earth Syst Sci. 2007;11:1633–44.View ArticleGoogle Scholar
- Pennington RT, Prado DE, Pendry CA. Neotropical seasonally dry forests and Quaternary vegetation changes. J Biog. 2000;27(2):261–73.View ArticleGoogle Scholar
- Young DG, Duncan MA. Guide to the identification and geographic distribution of Lutzomyia sand flies in Mexico, the West Indies, Central and South America (Diptera: Psychodidae). Memoirs of the Am Entom Inst. 1994;54:1–881.Google Scholar
- Solano P, Duvallet G, Dumas V, Cuisance D, Cuny G. Microsatellite markers for genetic population studies in Glossina palpalis (Diptera: Glossinidae). Acta Trop. 1997;65:175–80.View ArticlePubMedGoogle Scholar
- Mazzoni CJ, Gomes CA, Souza NA, De Queiroz RG, Justiniano SCB, Ward RD, et al. Molecular evolution of the period gene in sandflies. J Mol Evol. 2002;55:553–62.View ArticlePubMedGoogle Scholar
- Standen R. The Staden sequence analysis package. Mol Biotechnol. 1996;5:233–41.View ArticleGoogle Scholar
- Ewing B, Green P. Base-calling of automated sequencer traces using Phred? II. Error? Probabilities. Genome Res. 1998;8:186–94.View ArticlePubMedGoogle Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA 5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;10:2731–39. doi:10.1093/molbev/msr121.View ArticleGoogle Scholar
- Rozas J, Rozas R. DnaSP version 3: an integrated program for molecular population genetics and molecular evolution analysis. Bioinformatics. 1999;15:174–75.View ArticlePubMedGoogle Scholar
- Crooks GE, Hon G, Onchandia JM, Brenner SE. WebLogo: a sequence logo generator. Genome Res. 2004;14:1188–90.View ArticlePubMed CentralPubMedGoogle Scholar
- Pritchard JK, Stephens P, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155:945–59.PubMed CentralPubMedGoogle Scholar
- Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software structure: a simulation study. Mol Ecol. 2005;14:2611–20.View ArticlePubMedGoogle Scholar
- Rousset F. Genetic differentiation and estimation of gene flow from F statistics under isolation by distance. Genetics. 1997;145:1219–28.PubMed CentralPubMedGoogle Scholar
- Excoffier L, Laval G, Schneider S. Arlequin ver. 3.0: an integrated software package for population genetics data analysis. EvolBioinform Online. 2005;1:47–50.Google Scholar
- Posada D. jModel test: phylogenetic model averaging. Mol Biol Evol. 2008;25:1253–56.View ArticlePubMedGoogle Scholar
- Guindon S, Gascuel O. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003;52:696–704.View ArticlePubMedGoogle Scholar
- Passamonti M, Mantovani B, Scali V. Phylogeny and karyotype evolution of the Iberian Leptyniaattenuata species complex (Insecta: Phasmatodea). Mol Phylogenet Evol. 2004;30:87–96.View ArticlePubMedGoogle Scholar
- Lehr MA, Kilpatrick CW, Wilkerson RC, Conn E. Cryptic species in the Anopheles (Nyssorhynchus) albitarsis (Diptera: Culicidae) complex: incongruence between random amplified polymorphic dna-polymerase chain reaction identification and analysis of mitochondrial DNA COI gene sequences. Ann Entomol Soc Am. 2005;98:908–17.View ArticlePubMed CentralPubMedGoogle Scholar
- Wilkerson RC, Foster PG, Li C, Sallum MAM. Molecular phylogeny of neotropical Anopheles (Nyssorhynchus) albitarsis species complex (Diptera: Culicidae). Ann Entomol Soc Am. 2005;98:918–25.View ArticlePubMed CentralPubMedGoogle Scholar
- Lins RMMA, Souza NA, Peixoto AA. Genetic divergence between two sympatric species of the Lutzomyia longipalpis complex in the paralytic gene, a locus associated with insecticide resistance and love song production. Mem Inst Oswaldo Cruz. 2008;103:736–40.View ArticlePubMedGoogle Scholar
- Hamilton JGC, Maingon R, Alexander B, Ward R, Brazil RP. Analysis of the sex pheromone extract of individual male Lutzomyia longipalpis sandflies from six regions in Brazil. Med Vet Entomol. 2005;19:480–88.View ArticlePubMedGoogle Scholar
- Souza NA, Vigoder FM, Araki AS, Ward RD, Kyriacou CP, Peixoto A. Analysis of the copulatory courtship songs of Lutzomyia longipalpis in six populations from Brazil. J Med Entomol. 2004;41:906–13.View ArticlePubMedGoogle Scholar
- Coutinho-Abreu IV, Sonoda IV, Fonseca JA, Melo MA, Balbino VQ. Lutzomyia longipalpis s.l. in Brazil and the impact of the Sao Francisco River in the speciation of this sand fly vector. Parasit Vectors. 2008;1:16.View ArticlePubMed CentralPubMedGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.