Phenotypic differentiation in love song traits among sibling species of the Lutzomyia longipalpis complex in Brazil
© Vigoder et al. 2015
Received: 28 November 2014
Accepted: 15 May 2015
Published: 28 May 2015
Brazilian populations of Lutzomyia longipalpis may constitute a complex of cryptic species, and this report investigates the distribution and number of potential sibling species. One of the main differences observed among Brazilian populations is the type of acoustic signal produced by males during copulation. These copulation song differences seem to be evolving faster than neutral molecular markers and have been suggested to contribute to insemination failure observed in crosses between these sibling species. In previous studies, two main types of copulation songs were found, burst-type and pulse-type. The latter type can, in turn, be further subdivided into five different patterns.
We recorded male song from 13 new populations of the L. longipalpis complex from Brazil and compared the songs with 12 already available.
Out of these 25 populations, 16 produce burst-type and 9 produce pulse-type songs. We performed a principal component analysis in these two main groups separately and an additional discriminant analysis in the pulse-type group. The pulse-type populations showed a clear separation between the five known patterns with a high correspondence of individuals to their correct group, confirming the differentiation between them. The distinctiveness of the burst-type subgroups was much lower than that observed among the pulse-type groups and no clear population structure was observed. This suggests that the burst-type populations represent a single species.
Overall, our results are consistent with the existence in Brazil of at least six species of the L. longipalpis complex, one with a wide distribution comprising all the populations with burst-type songs, and five more closely related allopatric siblings with different pulse-type song patterns and more restricted distribution ranges.
Understanding speciation is one of the central questions in evolutionary biology. Many authors consider sexual selection to be one of the main causes of speciation, creating reproductive barriers that can prevent gene flow [1, 2]. Acoustic communication has been implicated in sexual selection and can act as a recognition signal in many animals, from insects to primates, with sibling species showing distinct songs [3–8]. Hence song variation could represent an important phenotype for understanding patterns of speciation in species complexes.
Lutzomyia longipalpis (Diptera: Psychodidae) males produce a song to females during copulation that may contribute to reproductive success [9–11]. This sand fly is the main vector of American visceral leishmaniasis and it constitutes a species complex [12–14]. However, the distribution and number of sibling species is still unclear, particularly in Brazil [11, 15, 16]. Among the evidences for the existence of cryptic species in Brazil is the reproductive isolation observed in laboratory crosses between some populations [12, 17, 18] and the fact that they have differences in phenotypic characters, such as male copulation songs and sex pheromones, and moderate to high levels of genetic divergence in some molecular markers such as microsatellites and several protein coding genes, including some associated with sexual behaviour in Drosophila [11, 13, 16, 18–30].
The reproductive isolation observed between Brazilian populations of the L. longipalpis complex is caused by insemination failure during copulation in crosses between only certain populations [13, 17, 18]. This suggests that after copulation has started, there is signalling that is important for insemination to occur. Because males sing during copulation and the insemination failure is observed in crosses between populations producing different types of copulation songs, this acoustic signal might have a role in species recognition, acting as a reproductive barrier reducing gene flow and potentially preventing gamete wastage [10, 11].
The different cryptic species of the L. longipalpis complex can be separated into two main groups according to the types of copulation song which males produce: pulse-type and burst-type [10, 11]. The Burst-type song is composed of trains with highly polycyclic pulses (“bursts”) modulated in frequency and amplitude. The pulse-type group on the other hand, has previously been shown to include five different song patterns designated as P1, P2, P3, P4 and P5 [10, 11].
Here we test the variability of song amongst L. longipalpis using a more comprehensive geographic sampling analysing 13 new Brazilian populations and comparing them with 12 populations analysed by Araki et al. . We examine the number of phenotypic song types found within each of the pulse- and burst-type song patterns, specifically asking if the pulse-type group represents multiple species and the burst-type only one. Overall, our results are consistent with the existence of at least six species of the L. longipalpis complex in Brazil, one with a wide distribution comprising all the populations with burst-type songs despite some level of geographic structuring, and five more closely related allopatric siblings with different pulse-type song patterns and more restricted distribution ranges.
Sand flies were collected using CDC light traps and the L. longipalpis individuals were identified according to Young & Duncan . Samples were obtained from the following localities: Afonso Claudio (September 2009) (20°04'S, 41°08'W) (Espírito Santo State); Ipanema (March 2010) (19°48′S, 41°42′ W), Nova Porteirinha (December 2007) (15°48′S, 43°18′W) and Lassance (January 2009) (17°53′S, 44°34′W) (Minas Gerais State); Pirenópolis (December 2007) (15°51′S, 48°57′W) (Goiás State); Aracajú (October 2010) (10°54′S, 37°4′W) (Sergipe State); Itamaracá (September 2007) (7°45′S, 34°51′W) and Passira (December 2007) (7°56′S, 35°35′W) (Pernambuco State); Barcarena (December 2007) (1°31′S, 48°37′W), Cametá (February 2010) (2°15′S, 49°30′W) and Camará (June 2010) (2°25′S, 54°43′W) (Pará State); and Palmas (July 2006 and April 2010) (10°10′S, 48°19′W) (Tocantins State).
Lutzomyia longipalpis males show a polymorphism in the number of abdominal spots, which are sex pheromone glands. Males can have either one pair in the fourth tergite (called 1S) or two pairs in the third and fourth tergites (called 2S). The second pair in the third tergite can also be smaller than the pair in the fourth and this is sometimes called the intermediate form [13, 17]. Although the spot polymorphism usually has no taxonomical value , in localities where sympatric species producing either pulse-type or burst-type song are found, males can also be differentiated by the number of spots and in these cases individuals with an intermediate spot phenotype are rare [10, 11, 20, 21, 24]. Of the new samples analysed in the present work we observed 1S and 2S males in sympatry only in the locality of Palmas, therefore this sample was subdivided into Palmas 1S and Palmas 2S.
The recordings were made in accordance with Souza et al.  using a male and a female in each trial. Insects were placed in a small chamber inside an INSECTAVOX , which contains the microphone. Trials were also filmed using a Sony Hi8CCD-TRV65 video camera. Both sound and video were recorded using a Panasonic DMR-ES10 DVD recorder. All trials were performed at 25 °C ± 1 °C and if no copulation occurred in 4–5 min the couple was replaced. Most recordings were performed using either wild caught sand flies or their F1 progeny.
After recording, songs were digitized using CED 1401 and analysed using Spike 2 software (version 4.08), both from Cambridge Electronic Design, United Kingdom. The parameters measured were: inter-pulse or inter-burst interval (IPI/IBI), number of pulses or bursts per train (NP/NB), train length (TL) and the carrier frequency of the pulse train (Freq). In the pulse-type populations (see Results) we also measured the cycles per pulse (CPP) and the proportion of pulses with alternated amplitudes (AmpAlt), which is the proportion of pulses that have amplitudes either higher or lower than both adjacent pulses. The alternated amplitude analysis was also performed in the populations previously described by Souza et al.  and Araki et al. . Note that in the previous analysis of AmpAlt carried out in Araki et al.  for a couple of populations, this song parameter was computed based on a simple visual inspection of the song traces and small differences in pulse amplitude were therefore ignored. In the current analysis this was carried out automatically based on exact amplitude values.
Statistical analysis was performed using SPSS version 19 software. All statistical analyses were done using the new populations reported here combined with data from the populations analysed by Souza et al.  and Araki et al. .
Males from all the new samples produced song during copulation, i.e., after they grabbed the female with their genital clamp. The copulation songs showed the same characteristics described previously for L. longipalpis, with primary and secondary songs [10, 11]. The analysis here concentrates on the primary song as it shows more variation among populations with a clear and consistent pattern and is produced in every copulation, while the secondary song is a long train of song with an uneven distribution of polycyclic pulses not produced by every male (see Araki et al. ). Also even when it is produced, the secondary song is similar in all populations.
The P1 songs are composed of trains of similar pulses with usually two or three cycles per pulse. P2 songs differ from P1 mainly by the presence of interspersed polycyclic pulses between nearly monocyclic pulses. The P3 songs are characterized by an almost perfect alternation of high and low amplitude pulses. P4 also presents a modulation in pulse amplitude but this occurs throughout the train with a more gradual amplitude oscillation or with alternating sequences of 2 or more high and low amplitude pulses. Finally, P5 is characterized by having a very short inter-pulse interval (~35 ms while other Pulse songs have > 50 ms) with some very polycyclic pulses in the end of the train (Fig. 1).
Mean (±SE) values of all parameter analysed in the pulse-type populations and their respective pattern
Loadings of the principal components analysis in the pulse-type populations
The five groups had a perfect 100 % correct assignment of individuals to each original group, confirming that the five song types are quantitatively completely distinct.
Mean (±SE) values of all parameter analysed in the burst-type populations
Loadings of the principal components analysis in the burst-type populations
The occurrence of cryptic species is common among insects and very often they can be separated by phenotypic traits involved in mating . In L. longipalpis sensu lato, the male song produced during copulation clearly suggests the existence of cryptic species within this taxon. The different song patterns observed in Brazilian populations suggests that song is evolving faster than other phenotypic traits, probably because it is under sexual selection [2, 3, 5, 34]. Characteristics under sexual selection tend to evolve faster than other traits making them good markers to differentiate closely related species [35, 36]. Acoustic communication has often been implicated in sexual selection  and this is likely that it has this role in the L. longipalpis species complex since the songs are produced during copulation and some studies suggest the existence of a mechanism of recognition during mating that is important for insemination success [13, 18].
Until now, 9 populations have been found to produce five distinct patterns of pulse-type songs (Fig. 7). These five patterns or subtypes not only have qualitative differences, but the discriminant analysis also demonstrates quantitative variation that allows high discrimination among them (Fig. 4). These five groups also show geographical separation and no overlap between their distributions have yet been observed (Fig. 7).
Statistical analysis showed clear differentiation among the five groups and the discriminant function analysis found that 100 % of the individuals could be assigned to the correct group based on song pattern. P2 song type was found so far in a single locality (Lapinha Cava) in south-eastern Brazil (Fig. 7) and the P1 song is present only in the north-eastern state of Bahia (Fig. 7). P5 was also found in a single locality (Estrela de Alagoas) in the Northeast region, in sympatry with a burst-type population in (Fig. 7).
Consistent with the available molecular data [11, 25], the song analysis also suggests that P3 and P4 are the two closest related groups among the pulse-type groups, since they both have pulse amplitude variation being differentiated mainly by the type of this variation seen, either an almost perfect alternation between high and low amplitude pulses (P3) or a more continuous variation throughout the train (P4) (Fig. 2) and also having more similar values in most parameters (Table 1). Such amplitude variation is extremely unusual in dipteran song, but presumably could be detected by females as song is always produced in copula. The currently known distribution places P3 populations in the far northern part of north-eastern Brazil while P4 is found in the more central part of Brazil (Fig. 7).
Further investigation is necessary to confirm that all five different pulse-type patterns represent distinct species and if the songs are associated with reproductive isolation among them, especially as they are allopatric in nature. The divergence in song patterns observed among the different pulse-type L. longipalpis populations is similar to the variation seen among some closely related Drosophila species [37, 38], so it is quite possible that they represent distinct species. In addition, populations with different song types also produce different pheromones [11, 23] and crossing experiments have found reproductive isolation between at least two of these populations: Jacobina (P1) and Lapinha (P2) . Moreover pulse-type populations with different patterns tend to have a considerable level of genetic divergence [11, 16, 17, 22].
Some small differences were also observed amongst the burst-type populations, but not close to the same extent as the Pulse-type populations. This group has a much wider geographical distribution in Brazil, ranging from the Southeast region up to the North region, crossing several different ecosystems (Fig. 7). Despite the similar song pattern, the statistical analysis of the copulation song of this group showed a difference between Jaíba 2S and a pair of other population. This difference could be an indication of the beginning of a separation of the Jaíba 2S population from the other Burst-type populations. However, the fact that the other 7 populations did not have significant differences suggests that even if it were true, the separation would be very recent and is qualitatively much less significant than that seen among pulse populations.
It is important to note that crossing experiments using different populations with burst-type songs showed they have normal insemination rates . Finally, all burst-type populations that have had their sex pheromone analysed produced the same compound (cembrene) . They have an overall smaller mean Fst (a widely used measure of genetic differentiation) [39, 40] of ~0.16 in the per gene in comparison among them [11, 13, 23] in contrast to the mean Fst observed among pulse-type populations, ~0.26 . All this strongly suggests that the burst-type group constitute a single species.
In Palmas, as previously observed for the localities of Sobral, Jaíba and Estrela [10, 11], males of two sympatric species can be distinguished by the spot phenotype, 1S or 2S. Palmas 1S males produce pulse-type songs (P4) while Palmas 2S males produce burst-type songs. The occurrence of sympatric species in many localities is one of the strongest pieces of evidence that L. longipalpis is a species complex in Brazil. It also suggests the absence of very large ecological differences between the sibling species, reinforcing the idea that sexual selection is acting as the primary force of speciation in the L. longipalpis complex .
A striking aspect of the geographic distribution of the L. longipalpis complex is the large distribution of the burst-type species coupled with a low level of phenotypic variation and the restricted distribution and high level of variation among the pulse-type populations (Fig. 7). What could cause these contrasting patterns of differentiation between the two song-type groups? Multilocus coalescent analysis estimates a recent separation between the burst- and pulse types around 0.5 mya , so if the differences between pulse-types have evolved since then, the evolution of these song types has been extremely rapid. The biogeographic history of the region is poorly understood, but it is possible that after an initial separation between burst- and pulse populations the ecological and climatic changes during the Quaternary could have caused vicariance of the pulse-type populations in refugia with more stable climatic zones , whereas the burst-type was confined to a single refugium with a more rapid expansion following climatic amelioration. If the pulse-types had a more complex history of secondary contacts between developing forms , there would have been more potential for character displacement to influence song and pheromone divergence. More biological reasons for greater divergence could include greater dispersal and gene flow among the burst-type or higher effective population sizes. If sexual selection is stronger in the pulse-type due to mating system variation then mating signals may diverge more quickly, leading to greater levels of sexual isolation among species. Little is known about mating rates or variance in mating success in natural populations of these species, but some of these intriguing hypotheses could be investigated by analyses of broader patterns of genomic divergence between the forms .
Our results show clear qualitative and quantitative variation among the male song patterns found in Brazilian populations of the L. longipalpis complex. These results confirm that the acoustic signals are a very good marker to differentiate these potential cryptic species. Six song patterns were observed, five among pulse-type populations and another one among burst-type populations. The statistical analysis shows that they can be easily differentiated suggesting that populations producing each pattern belong to a different cryptic species of the L. longipalpis complex. Further studies will be important to better understand the mechanism by which the song may be contributing to reproductive isolation.
We would like to thank Gabriel E. M. Ferreira and Fernando José Silva for the help with field collections. We dedicate this work to the memory of Dr Alexandre A. Peixoto.
This work was supported by grants from the Howard Hughes Medical Institute, Conselho Nacional de Desenvolvimento Científico e Tecnológico, Fundação de Amparo à Pesquisa do Estado de São Paulo Coordenação de Aperfeiçoamento de Pessoal de Nível Superior and Fundação Oswaldo Cruz.
- Kraaijeveld K, Kraaijeveld-Smit FJ, Maan ME. Sexual selection and speciation: the comparative evidence revisited. Biol Rev Camb Philos Soc. 2011;86:367–77.PubMedView ArticleGoogle Scholar
- Rodriguez RL, Boughman JW, Gray DA, Hebets EA, Höbel G, Symes LB. Diversification under sexual selection: the relative roles of mate preference strength and the degree of divergence in mate preferences. Ecol Lett. 2013;16:964–74.PubMed CentralPubMedView ArticleGoogle Scholar
- Gleason JM, Ritchie MG. Evolution of courtship song and reproductive isolation in the Drosophila willistoni species complex: Do sexual signals diverge the most quickly? Evolution. 1998;52:1493–500.View ArticleGoogle Scholar
- Ritchie MG, Halsey EJ, Gleason JM. Drosophila song as a species-specific mating signal and the behavioural of Kyriacou & Hall cycles in Drosophila melanogaster song. Anim Behav. 1999;58:649–57.PubMedView ArticleGoogle Scholar
- Mendelson TC, Shaw KL. Close-range acoustic signaling and mate choice in Hawaiian crickets (Gryllidae: Laupala). Behav Ecol Sociobiol. 2006;59:770–6.View ArticleGoogle Scholar
- Reinhold K. Variation in acoustic signalling traits exhibits footprints of sexual selection. Evolution. 2011;65:738–45.PubMedView ArticleGoogle Scholar
- Thinh VN, Hallam C, Roos C, Hammerschmidt K. Concordance between vocal and genetic diversity in crested gibbons. BMC Evol Biol. 2011;11:36.PubMed CentralPubMedView ArticleGoogle Scholar
- de Mello VF, Ritchie MG, Gibson G, Peixoto AA. Acoustic communication in insect disease vectors. Mem Inst Oswaldo Cruz. 2013;108 Suppl 1:26–33.Google Scholar
- de Souza NA, Ward RD, Hamilton JGC, Kyriacou CP, Peixoto AA. Copulation songs in three siblings of Lutzomyia longipalpis (Diptera: Psychodidae). T Roy Soc Trop Med H. 2002;96:102–3.View ArticleGoogle Scholar
- Souza NA, Vigoder FM, Araki AS, Ward RD, Kyriacou CP, Peixoto AA. Analysis of the copulatory courtship songs of Lutzomyia longipalpis in six populations from Brazil. J Med Entomol. 2004;41:906–13.PubMedView ArticleGoogle Scholar
- Araki AS, Vigoder FM, Bauzer LG, Ferreira GE, 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, e365.PubMed CentralPubMedView ArticleGoogle Scholar
- Mangabeira O. Sobre a sistemática e biologia dos flebótomos do Ceará. Rev Bras Malariol Doencas Trop. 1969;21:3–26.Google Scholar
- Ward RD, Phillips A, Burnet B, Marcondes CB. The Lutzomyia longipalpis complex: reproduction and distribution. In: Service MW, editor. Biosystematics of Haematophagous Insects. Oxford: Oxford University Press; 1988. p. 258–69.Google Scholar
- 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:839–47.PubMedGoogle Scholar
- Arrivillaga J, Mutebi JP, Pinango H, Norris D, Alexandre B, Feliciangeli MD, et al. 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:615–27.PubMedView ArticleGoogle Scholar
- Watts PC, Hamilton JG, Ward RD, Noyes HA, Souza NA, Kemp SJ, et al. 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:734–43.PubMedGoogle Scholar
- Ward RD, Ribeiro AL, Ready PR, Murtagh A. Reproductive isolation between different forms of Lutzomyia longipalpis (Lutz & Neiva), (Diptera: Psychodidae), the vector of Leishmania donovani chagasi Cunha & Chagas and its significance to kala-azar distribution in South America. Mem Inst Oswaldo Cruz. 1983;8:269–80.Google Scholar
- Souza NA, Andrade-Coelho CA, Vigoder FM, Ward RD, Peixoto AA. Reproductive isolation between sympatric and allopatric Brazilian populations of Lutzomyia longipalpis sl (Diptera: Psychodidae). Mem Inst Oswaldo Cruz. 2008;103:216–9.PubMedView ArticleGoogle Scholar
- Bauzer LGSR, Souza NA, Ward RSJ, Kyriacou CP, Peixoto AA. The period gene and genetic differentiation between three Brazilian populations of Lutzomyia longipalpis. Insect Mol Biol. 2002;11:315–23.PubMedView ArticleGoogle Scholar
- Bauzer LGSR, Gesto JSM, Souza NA, Ward RD, Hamilton JG, Kyriacou CP, et al. Molecular divergence in the period gene between two putative sympatric species of Lutzomyia longipalpis complex. Mol Biol Evol. 2002;19:1624–7.PubMedView ArticleGoogle Scholar
- Maingon RDC, Ward RD, Hamilton JGC, Noyes HA, Souza N, Kemp SJ, et al. Genetic identification of two sibling species of Lutzomyia longipalpis (Diptera: Psychodidae) that produce distinct male sex pheromones in Sobral, Ceará State, Brazil. Mol Ecol. 2003;12:1879–94.PubMedView ArticleGoogle Scholar
- Bottecchia M, Oliveira SG, Bauzer LGSR, Souza NA, Ward RD, Garner KJ, et al. Genetic divergence in the cacophony IVS6 intron among five Brazilian populations of Lutzomyia longipalpis. J Mol Evol. 2004;58:754–61.PubMedView ArticleGoogle Scholar
- 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:480–8.PubMedView ArticleGoogle Scholar
- Lins RM, 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 lovesong production. Mem Inst Oswaldo Cruz. 2008;103:736–40.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- Araki AS, Ferreira GEM, Mazzoni CJ, Souza NA, Machado RC, Bruno RV, et al. Multilocus analysis of divergence and introgression in sympatric and allopatric sibling species of the Lutzomyia longipalpis complex in Brazil. PLoS Negl Trop Dis. 2013;7, e2495.PubMed CentralPubMedView ArticleGoogle Scholar
- Santos M, Ribolla P, Alonso D, Andrade-Filho J, Casaril A, Ferreira A, et al. Genetic structure of Lutzomyia longipalpis populations in Mato Grosso do Sul, Brazil, based on microsatellite markers. PLoS One. 2013;8, e74268.PubMed CentralPubMedView ArticleGoogle Scholar
- 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–12.PubMedView ArticleGoogle Scholar
- 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:12–7.PubMedView ArticleGoogle Scholar
- Kyriacou CP, Hall JC. Circadian rhythm mutations in Drosophila melanogaster affect short-term fluctuations in the male’s courtship song. PNAS. 1980;77:6729–33.PubMed CentralPubMedView 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). Mem Amer Ent Inst. 1994;54:1–881.Google Scholar
- Gorczyca M, Hall JC. The INSECTAVOX, an integrated device for recording and amplifying courtship songs of Drosophila. Dros Inf Serv. 1987;66:157–60.Google Scholar
- Bickford D, Lohman DJ, Sodhi NS, Ng PK, Meier R, Winker K, et al. Cryptic species as a window on diversity and conservation. Trends Ecol Evol. 2007;22:148–55.PubMedView ArticleGoogle Scholar
- Ritchie MG. Sexual Selection and Speciation. Ann Rev Ecol Evol Syst. 2007;38:79–102.View ArticleGoogle Scholar
- Mullen SP, Mendelson TC, Schal C, Shaw KL. Rapid evolution of cuticular hydrocarbons in a species radiation of acoustically diverse Hawaiian crickets (Gryllidae: Trigonidiinae: Laupala). Evolution. 2007;61:223–31.PubMedView ArticleGoogle Scholar
- Safran RJ, Vitousek MN, Hauber ME, Ghalambor CK. Sexual selection: a dynamic state of affairs. Response to the comments of Cornwallis and Uller in the article: towards an evolutionary ecology of sexual traits. Trends Ecol Evol. 2010;25:429–30.PubMedView ArticleGoogle Scholar
- Ritchie MG, Gleason JM. Rapid evolution of courtship song pattern in Drosophila willistoni sibling species. J Evol Biol. 1995;8:463–79.View ArticleGoogle Scholar
- Costa CTA, Sene FM. Characterization of courtship sounds of species of the subgroup fasciola (Diptera, Drosophilidae, Drosophila repleta group): interspecific and interpopulational analyses. Braz J Biol. 2002;62:573–83.PubMedView ArticleGoogle Scholar
- Weir BS, Clark Cockerham C. Estimating F-Statistics for the analysis of population structure. Evolution. 1984;38(6):1358.View ArticleGoogle Scholar
- Holsinger KE, Weir BS. Genetics in geographically structured populations: defining, estimating and interpreting FST. Nat Rev Genet. 2009;10(9):639–50.PubMedView ArticleGoogle Scholar
- Willis KJ, Whittaker RJ. The Refugial Debate. Science. 2000;287:1406–7.PubMedView ArticleGoogle Scholar
- Hewitt GM. Speciation, hybrid zones and phylogeography - or seeing genes in space and time. Mol Ecol. 2001;10:537–49.PubMedView ArticleGoogle Scholar
- Hey J. Isolation with migration models for more than two populations. Mol Biol Evol. 2010;27:905–2.PubMed CentralPubMedView ArticleGoogle 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.