Signals of demographic expansion in Drosophila virilis

  • Patricia M Mirol1, 4Email author,

    Affiliated with

    • Jarkko Routtu2,

      Affiliated with

      • Anneli Hoikkala2 and

        Affiliated with

        • Roger K Butlin3

          Affiliated with

          BMC Evolutionary Biology20088:59

          DOI: 10.1186/1471-2148-8-59

          Received: 12 September 2007

          Accepted: 25 February 2008

          Published: 25 February 2008

          Abstract

          Background

          The pattern of genetic variation within and among populations of a species is strongly affected by its phylogeographic history. Analyses based on putatively neutral markers provide data from which past events, such as population expansions and colonizations, can be inferred. Drosophila virilis is a cosmopolitan species belonging to the virilis group, where divergence times between different phylads go back to the early Miocene. We analysed mitochondrial DNA sequence variation among 35 Drosophila virilis strains covering the species' range in order to detect demographic events that could be used to understand the present characteristics of the species, as well as its differences from other members of the group.

          Results

          Drosophila virilis showed very low nucleotide diversity with haplotypes distributed in a star-like network, consistent with a recent world-wide exponential expansion possibly associated either with domestication or post-glacial colonization. All analyses point towards a rapid population expansion. Coalescence models support this interpretation. The central haplotype in the network, which could be interpreted as ancestral, is widely distributed and gives no information about the geographical origin of the population expansion. The species showed no geographic structure in the distribution of mitochondrial haplotypes, in contrast to results of a recent microsatellite-based analysis.

          Conclusion

          The lack of geographic structure and the star-like topology depicted by the D. virilis haplotypes indicate a pattern of global demographic expansion, probably related to human movements, although this interpretation cannot be distinguished from a selective sweep in the mitochondrial DNA until nuclear sequence data become available. The particular behavioural traits of this species, including weak species-discrimination and intraspecific mate choice exercised by the females, can be understood from this perspective.

          Background

          The development of methods to analyse intraspecific phylogenies has provided very valuable tools to understand how populations have been influenced by historical and contemporary processes [1, 2]. Mitochondrial DNA has proved a very useful tool for reconstructing phylogenies of species and historical demographic events [38]. Its usefulness for the study of closely related taxa and populations within species lies in its very low rate of recombination, maternal inheritance, conserved structure, reduced effective population size and relatively high rates of evolution [6, 7, 9]. However, after an explosion of studies using mitochondrial DNA for phylogeographic inference, there is nowadays rising concern about making inferences based in this single molecule, due to problems related with recombination, effective population size, mutation rates, introgression and neutrality [10, 11]. These problems influence the construction of species-level phylogenies most profoundly [e.g. [12]]. However, they may also confound estimations of demographic history and coalescence times within species. Nevertheless, if these drawbacks are taken into account when examining the evolutionary history of a taxon, mitochondrial DNA surveys can provide efficient means of detecting gene flow, levels of reproductive isolation, species boundaries and historical patterns of population structure [13].

          The pattern of variation of mitochondrial DNA, permits inferences about a species' demographic history, including events such as population expansions and colonizations. The application of coalescent theory [14, 15] for the analysis of population data has allowed inferences about past and present population size. A change in population size over time will be reflected in the DNA sequences in such a way that the analysis of the substitutions can indicate the direction and the timing of the change. It has been shown that demographic expansions lead to star-shaped genealogies [16], an excess of rare mutations [17] and a unimodal mismatch distribution [18]. Although spatial or range expansions also lead to an increase in the global effective size of a species, they lead to the same molecular signature as that of a purely demographic expansion if the migration rate between demes was large [19, 20]. Otherwise, for relatively low levels of gene flow, multimodal mismatch distributions should be expected [19]. During the Pleistocene, global climate oscillations and the associated cycles of glaciation had a profound influence on the population distributions of many organisms [2124]. Populations and species retracted their habitats to limited areas which served as refuges. Range expansions and colonizations followed by demographic expansions often happened after these events, when climate ameliorated, and it is possible to find the signals of those processes in contemporary populations.

          The Drosophila virilis group comprises 13 species and subspecies divided into two clades, the virilis and montana phylads [2530]. Throckmorton [25] suggested that the phylads diverged in the early Miocene, or not later than the Pliocene, when both of them entered the New World by way of Beringia. Divergence times between the phylads have been estimated to be from 7 Mya [31] to 11 Mya [29]. In general, species within the montana phylad have evolved more in terms of chromosomal re-arrangements, and are also more variable regarding the number of inversions segregating within populations than members of the virilis phylad [25], although species in this group show a higher number of fusions.

          Drosophila virilis group species have been used as models in evolutionary studies for decades and their role in this kind of research will become even more important now that the whole genome of D. virilis has been sequenced http://​insects.​eugenes.​org/​species/​data/​dvir/​. Within the virilis phylad, Drosophila virilis is a domestic species of nearly cosmopolitan distribution. The primitive karyotype is found in D. virilis and natural populations are exceptional in having no inversion polymorphism, in contrast to other species of the group. It has been proposed that D. virilis originated somewhere in the ancient deciduous forests of China or in arid regions such as Iran or Afghanistan, and remained isolated from the remaining species of the group until relatively recently [25]. Its ancestral role is supported by its monomorphic karyotype, an extremely high (about 50%) proportion of satellite DNA, and the distribution of various repetitive DNA elements [32].

          Drosophila virilis also differs from other species of the group by the weak species-discrimination and intraspecific mate choice exercised by the females [e.g. [25, 33]]. Although females are able to recognize species-specific characters of the male song, they do not require to hear the song before mating [34]. D. virilis strains originating from different parts of the species distribution have shown significant variation in male courtship song [33]. Rapid expansion into a novel habitat, not occupied by related species, may explain these characteristics. The ability of D. virilis to survive in various kinds of environments can partly be explained by its high thermotolerance [35] and high tolerance for ethanol [36]. Like several other cosmopolitan species, e.g. D. melanogaster, D. virilis possesses intraspecific genetic differences in thermotolerance with an obvious adaptive significance to local thermal conditions [37], which may make it easier for the species to occupy new habitats.

          In this paper, we report an analysis of the demographic history of D. virilis using mitochondrial DNA sequence data, in order to provide an historical framework for evolutionary studies on life history and behavioural traits of the species and comparisons with other species in the group, especially D. montana, the most divergent species within the virilis group.

          Results and Discussion

          We examined a total of 35 lines covering the species range: 9 from Japan, 9 from Eastern Europe and Central Asia, 5 from Western Europe, 2 from the United States and 10 from China (Table 1). These lines have been kept in captivity for widely varying numbers of years: since 1913 in the case of line 15010-1051.8 from California, and only since 2002 for the lines originating from wild flies collected in China.
          Table 1

          Lines of D. virilis used in the study, indicating, when it was available, year of collection and coordinates from which the line originates

          Drosophila virilis

          Line

          Year

          Coordinates

          Matsuyama, Japan

          A11

          1973

           
           

          A12

          1973

           
           

          B15

          1973

           
           

          B31

          1973

           
           

          B42

          1973

           

          Sapporo, Japan

          SBB

          1986

          43° 3'N, 141° 21'E

          Sakata, Japan

          SKT

          1987

           

          Japan

          Jap. Inv.

            

          Mishima, Japan

          W158

            

          Hangzhow, China

          15010-1051.47

          1948

           

          Human, China

          V-Hunan

            

          Zeziping, China

          V-ZZP-01

          2001

           

          Wuwei, China

          V-WW-03

          2002

           
           

          V-WW-05

          2002

           
           

          V-WW-08

          2002

           

          Lanzhou, China

          V-EH-01

          2002

           

          Dunghuang, China

          V-DNH

            

          Nanjing, China

          V-NANJING

            

          Qufu, China

          V-QUFU

            

          Russia

          15010-1051.52

          1976

           

          Baku, Azerbaijan

          1413

          1974

          40° 22'N, 49° 48'E

          Jalta, Ukraine

          1415

          1973

          46° 57'N, 37° 16'E

          Batumi, Georgia

          A

           

          41° 39'N, 41° 39'E

          Batumi, Georgia

          S9

          1970

          41° 39'N, 41° 39'E

          Yerevan, Armenia

          1

           

          40° 12'N, 44° 31'E

          Mzheta, Caucasus

          25

            

          Caucasus

          1411

          1973

          45°N, 50°E

          Tashkent, Uzbekistan

          12

           

          41° 18'N, 69° 16'E

          Holland

          W159

            

          Leeds, UK

          LeedsA

          1995

          53° 47'N, 1° 32'W

           

          LeedsB

          1995

           
           

          1430

          1981

          53°N, 1°E

           

          1433

          1982

          53°N, 1°E

          Pasadena, California, USA

          15010-1051.0

          1913

          34° 8'N, 118° 8'W

          Truckee, California, USA

          15010-1051.8

           

          39° 19'N, 120° 12'W

          The analysis was based on a total of 1358 base pairs of mitochondrial DNA sequence, 670 bp corresponding to the cytochrome oxidase I and 688 bp to the cytochrome oxidase II gene. There was a total of 19 haplotypes and, among them, 24 nucleotide substitutions, of which 20 were transitions and 4 transversions. Nucleotide diversity (per base) was 0.00185 ± 0.00114. The haplotype network (Figure 1), estimated using statistical parsimony, included one common haplotype, represented in six lines and present in all geographical regions except the USA. The second most frequent haplotype was found in four lines from China. The difference between these haplotypes was only two substitutions: in fact, the maximum number of differences between any pair of haplotypes was only 7 substitutions. There was no apparent geographical structure in the relationships among the haplotypes: the network was star-like with low levels of sequence divergence and a high frequency of unique mutations, indicating either a rapid population expansion or that selection has caused the rapid spread of a mitochondrial lineage carrying beneficial mutations. The limited number of lines sampled per region means that rare haplotypes specific to individual regions may have been missed. However, such haplotypes would contribute little to overall measures of differentiation.
          http://static-content.springer.com/image/art%3A10.1186%2F1471-2148-8-59/MediaObjects/12862_2007_Article_628_Fig1_HTML.jpg
          Figure 1

          Network obtained for theD. virilishaplotypes using statistical parsimony. Haplotypes are represented by ellipses, the area of the ellipse is proportional to the frequency of the haplotype, points in the lines connecting circles indicate substitutions. Shading denotes region: white – North America, pale grey – Japan, dark grey – Western Europe, black – Asia and Eastern Europe.

          When analysing demographic scenarios, where migration between populations, recombination and sexual reproduction usually play an important role, population history might not be tree-like. In these cases intraspecific genealogies constitute the appropriate way to make evolutionary inferences. The genealogical approach, which estimates parameters of the genealogical process, provides a coherent framework in which to consider recombination, migration, selection and other processes [38]. This approach takes into account several characteristics of the populations which make them different to species-level analysis, such as low divergence, ancestral nodes extant in the sampled population, multifurcations, reticulation and large sample sizes [1, 2].

          The mismatch distribution was smooth and unimodal (Figure 2a), but it departed from the expectations of the stepwise population expansion model (Harpending's raggedness index = 0.0721, p = 0.11, fit to the stepwise growth model, SSD = 0.0157, p = 0.03). There was an excess of unique mutations compared with that expected for a constant-size population (Table 2). We used the program FLUCTUATE to fit a model of exponential population expansion and estimate θ 0, the present-day scaled population size (Table 2). The site categories and the rate of change were obtained using a maximum likelihood heuristic search in PAUP, using HKY85 as preferred model, following MODELTEST, with 5 rate categories, transition/transversion ratio of 2.3 and shape parameter of 0.105. Model comparisons using GENIE confirm that the exponential growth model was a good fit to the data and this can be seen in the skyline plot (Figure 2b) which indicates a rapid and recent expansion (Table 2).
          http://static-content.springer.com/image/art%3A10.1186%2F1471-2148-8-59/MediaObjects/12862_2007_Article_628_Fig2_HTML.jpg
          Figure 2

          Mismatch distribution (a) and generalised skyline plot (b) among haplotypes ofD. virilis. (a): Expectations from the stepwise growth model, fitted in ARLEQUIN, are superimposed. (b): Observed values (solid line) and fitted values from the best model (broken line) – see Table 2. Smoothing parameter (epsilon) was 4E-5 (maximum likelihood value from option 'maxepsilon' in GENIE). Nu: number of females, Tu: generations.

          Table 2

          Summary statistics of mitochondrial variation and results from the FLUCTUATE and GENIE analyses.

                

          FLUCTUATE results

           

          GENIE results

            

          n

          H

          θS

          θπ

          Tajima's D

          Fu's F

          Θ0 (se)

          g (se)

          AICc Best model (second best model)

          Θ0 (95% CI)

          g (95% CI)

          35

          19

          0.0044 (0.0015)

          0.0019 (0.0011)

          -2.0055 p < 0.01

          -14.311 p < 0.0001

          0.0499 (0.0103)

          3532.6 (358.8)

          Exp 149.83 (Log 148.38)

          0.091 (0.027–0.410)

          3951 (2610–5433)

          Estimate (standard errors): n – number of lines surveyed, H – number of haplotypes, θS – diversity (Waterson 1975), θπ – diversity (Tajima 1989), θ 0 current estimate of 2Nμ, g – scaled population growth parameter. GENIE model abbreviations: Exp – exponential growth, Log – logistic growth.

          The Drosophila virilis lines analysed in this study did not show any geographic structure among mitochondrial haplotypes. There was a very low number of substitutions differentiating haplotypes and they were commonly shared among locations. All analyses point towards a rapid population expansion over a short time-scale, which could be consistent with population growth following the end of the last glaciation and/or a shift into domestic environments. These are the most obvious possible causes of expansion, although they suggest population growth would have initiated less than 10,000 years ago which is more recent than the timing implied by the skyline plot, for example, which suggests continuous population growth for >50,000 years (using a substitution rate of 10-8 per year). However, it is not possible to rule out alternative models, such as recent expansion from a smaller but more stable refugial population and time estimation should be treated with caution given the wide confidence interval on the expansion rate estimate and uncertainty about the substitution rate in D. virilis.

          In natural populations, ancestral haplotypes may persist and be sampled together with their descendants [2]. The coalescence process predicts that high frequency haplotypes are likely to have been present in the population for a long time [2]. Rare haplotypes represent more recent mutations and are more likely to be related to common haplotypes than to other rare variants [39]. Therefore, ancestral haplotypes have greater probability of being central in a network, being common and having a broader geographic distribution [2]. Following these considerations, the central haplotype in the network appears to be ancestral (Figure 1), and because it is widely distributed gives no information on the geographic origin of the population expansion.

          In their study on sequence variation in six X-linked genes of 21 D. virilis strains from different continents, Vieira and Charlesworth [40] found no fixed differences between the Asian strains and strains originating from Europe or North- or South-America, all the variants found outside Asia being also present in Asia, but not vice versa. These data were consistent with either a large population centred in Asia and a smaller migrant population elsewhere, or a large migrant population that went through a bottleneck [40]. Throckmorton [25] has previously suggested that D. virilis originated from an ancestral form in Asia, since the most primitive species of the virilis-repleta section of the genus Drosophila have been observed in Southeast Asia. This hypothesis recently received further support from a molecular phylogeny of the virilis section [41].

          In a recent study, Huttunen et al. [42] analysed 48 microsatellite loci in 30 D. virilis strains, many of which were analysed here for mitochondrial DNA. Although a phylogenetic tree and STRUCTURE analysis showed only moderate clustering of the strains originating from Continental Asia, Europe, America and Japan, an assignment test using a priori information about the geographical origin of the strains gave high posterior probabilities for their correct assignment. Microsatellite variation also showed significant population differentiation, as measured by FST, with evidence for isolation by distance. Variability detected by microsatellites can have a more recent origin than variation measured by mitochondrial DNA. Consequently, the lack of population differentiation in mitochondrial DNA could indicate shared ancestry while microsatellite variation between populations indicates that differentiation is in progress because current gene flow is restricted. Alternatively, the homogeneity of mitochondrial DNA could be the result of a recent selective sweep which did not disturb the pre-existing population structure for nuclear loci. Unfortunately, the use of laboratory strains prohibited testing for a demographic expansion using the microsatellite data.

          There are several papers describing genetic variation in another species of the virilis group, D. americana, which, together with D. lummei, represents dispersal of the virilis group into the Neartic region [4346]. Drosophila americana is represented by two chromosomal forms, D. a. americana and D. a. texana, differing by a fusion between the X and fourth chromosome and the frequency of chromosomal inversions. As in the case of D. virilis, phylogenetic analysis failed to resolve distinct clades between the chromosomal forms and also between geographic regions sampled, with Fst values not significantly different from zero [43, 44]. The absence of differentiation among the geographically distinct populations indicated either that gene flow homogenizes neutral sequence variation, even when different chromosomal rearrangements are involved, or that large populations retain shared ancestral polymorphisms [43]. Furthermore, and opposite to results obtained for D. virilis [42], an analysis of 27 microsatellite loci in 85 individuals from 6 different natural populations of D. americana showed no differentiation between populations [46], which confirms the existence of ongoing gene flow between populations of the species.

          The pattern of geographic variation in haplotypes of D. virilis, and in D. americana, is in sharp contrast with the pattern recently described for another species of the group, D. montana [47]. In this case there was clear differentiation, based on both mitochondrial DNA and microsatellites, between lines from populations in Finland, Canada and USA. Both markers indicated the presence of at least two distinct populations, one in Eurasia and the other one representing the expansion of the species to the New World, with a divergence time between them estimated from 450,000 to 900,000 years ago, within the Pleistocene. Although D. virilis and D. montana belong to the same species group in the genus Drosophila, they represent different phylads and differ in many characteristics, including chromosomal variation, habitat preferences, and courtship behaviour. Differences in biogeographic history, reflected in the pattern of mitochondrial variation in the two species, could have been important in the origin of these characteristics and constitute the basis for the interpretation of their evolution.

          Conclusion

          Mitochondrial DNA sequence variation in Drosophila virilis suggests a worldwide exponential population expansion during the Pleistocene, with extensive migration between demes. Slight differentiation at microsatellite markers results either from bottlenecking during or after expansion, or accumulation of differentiation across current barriers to gene exchange. The alternative possibility, that a selective sweep has homogenised mitochondrial sequences without disturbing pre-existing population structure for nuclear markers, cannot be excluded. However, a recent demographic expansion is consistent with the domestic habitat of D. virilis and helps to explain its lack of chromosomal polymorphism as well as evidence for a history of weak selection on mating behaviour by comparison with related species.

          Methods

          Drosophila stocks and sampling

          In total, 35 Drosophila virilis strains, covering the species' range, were selected for analysis (Table 1). The stocks were collected during a time period covering almost 90 years, from 1913 to flies sampled in China in 2001/2002. A single individual from each strain, either from laboratory stocks or freshly caught, was used to extract DNA and for PCR amplification of the COI and COII mitochondrial genes.

          Amplification and sequencing of mitochondrial DNA

          DNA was extracted from ethanol-preserved flies following a standard protocol [48], where the samples were homogenised in buffer and proteinase k, and DNA was extracted with chloroform-isoamyl alcohol and precipitated with isopropanol. The amplification of mitochondrial DNA was carried out with primers flanking the COII gene in the tRNALYS and tRNALEU [[49]; TL2: 5'-ATGGCAGATTAGTGCAATGG-3', TKN: 5'-GTTTAAGAGACCAGTACTTG-3'], which amplify an 850 bp fragment that includes the 688 bp COII gene, and COI-1460-F: 5'-ATCTATCGCCTAAACTTCAGCC-3' and COI-2195-R: 5'-ACTTCAGGGTGACCAAAAAATC-3' [50, 51] which amplify the complete 670 bp corresponding to the COI gene. PCR reactions were performed in 50 μl volumes including 0.5 μM of each primer, 200 μM dNTPs, 1.5 mM MgCl2 and 1 U Taq polymerase (Bioline) in reaction buffer. Initial denaturation was for 7 minutes at 94°C followed by 35 cycles of 1 minute at 94°C, 1 minute at the annealing temperature (54°C for COI and 56°C for COII) and 1 minute at 72°C, and a final incubation of 5 minutes at 72°C. The products were purified using QIAquick columns (QIAGEN) and sequenced using the forward primer. No ambiguous sites were found in the sequences. Sequences (GenBank accession nos. DQ426800 to DQ426823) were aligned with CLUSTAL-V [52].

          Mitochondrial DNA analysis of population history and phylogeography

          The partition homogeneity test, PHT [53], as implemented in PAUP 4.0, was used to test for incongruence between the COI and COII data sets. There was congruence between the data-sets and therefore the two fragments were combined for all subsequent analyses. ARLEQUIN 2.0 [54] was used to calculate pairwise distances between haplotypes, the mismatch distributions and tests of the standard neutral model for a demographically stable population (Tajima's D [55] and Fu's F [56]).

          The program FLUCTUATE [57] was used to make simultaneous estimates of present day θ and the population growth rate g, assuming an exponential model of growth and using a maximum likelihood approach. The parameters used for the simulations were obtained by running a hierarchy of likelihood-ratio tests in Modeltest 3.0 [58] to choose the model of evolution with the best fit to the data. Skyline plots were constructed using GENIE v. 3.0 [59]. The starting trees were obtained using maximum likelihood with molecular clock enforced. The software requires that the genealogy is estimated under the assumption of a molecular clock [59]. GENIE was also used to calculate the fit to different models of population growth, with fit assessed using the corrected Akaike Information Criterion.

          Networks of haplotypes were constructed based on statistical parsimony using the program TCS 1.06 [60].

          Declarations

          Acknowledgements

          We are grateful to the members of the 'Co-evolved Traits' Research Training Network for their valuable input to the work presented here and to the European Commission for funding the network (HPRN-CT-2002-00266). Special thanks are due to Zhang Wenxia for help in collecting D. virilis flies in China. Laboratory strains of the species were obtained also from Jorge Vieira, Michael Evgenev and Bowling Green stock center.

          Authors’ Affiliations

          (1)
          School of Biology, The University of Leeds
          (2)
          Department of Biological and Environmental Science, University of Jyväskylä
          (3)
          Animal and Plant Sciences, The University of Sheffield
          (4)
          Museo Argentino de Ciencias Naturales, CONICET

          References

          1. Emerson BC, Paradis E, Thëbaud C: Revealing the demographic histories of species using DNA sequences. Trends Ecol Evol 2001, 16: 707–716.View Article
          2. Posada D, Crandall KA: Intraspecific gene genealogies: trees grafting into networks. Trends Ecol Evol 2001, 16: 37–45.View ArticlePubMed
          3. Avise JC, Arnold J, Ball MR, Bermingham E, Lamb T, Neigel JE, Reeb CA, Saunders NC: Intraspecific phylogeography: the mitochondrial DNA bridge between population genetics and systematics. Ann Rev Ecol Syst 1987, 18: 489–522.
          4. Avise JC: Ten unorthodox perspectives on evolution prompted by comparative population genetic findings on mitochondrial DNA. Annu Rev Genet 1991, 25: 45–69.PubMed
          5. Avise JC: Phylogeography: the history and formation of species. Cambridge, MA: Harvard University Press 2000.
          6. Avise JC: Molecular markers, natural history and evolution. 2 Edition Sinauer Associates, Sunderland, Massachusetts 2004.
          7. Moritz C, Dowling TE, Brown WM: Evolution of animal mitochondrial DNA: relevance for population biology and systematics. Ann Rev Ecol Syst 1987, 18: 269–292.View Article
          8. Harrison RG: Animal mitochondrial DNA as a genetic marker in population and evolutionary biology. TREE 1989, 4: 6–11.
          9. Brown WM, George M, Wilson AC: Rapid evolution of animal mitochondrial DNA. Proc Natl Acad Sci USA 1979, 76: 1967–1971.View ArticlePubMed
          10. Shaw KL: Conflict between nuclear and mitochondrial DNA phylogenies of a recent species radiation: what mtDNA reveals and conceals about modes of speciation in Hawaiian crickets. Proc Natl Acad Sci USA 2002, 99: 16122–16127.View ArticlePubMed
          11. Ballard JW, Whitlock MC: The incomplete history of mitochondria. Mol Ecol 2004, 13: 729–744.View ArticlePubMed
          12. Funk DJ, Omland KE: Species-level paraphyly and polyphyly: frequency, causes and consequences, with insights from animal mitochondrial DNA. Ann Rev Ecol Syst 2003, 34: 397–423.View Article
          13. Rubinoff D, Holland BS: Between two extremes: mitochondrial DNA is neither the panacea nor the nemesis of phylogenetic and taxonomic inference. Syst Biol 2005, 54: 952–961.View ArticlePubMed
          14. Kingman JFC: The coalescent. Storch Proc Appl 1982, 13: 235–248.View Article
          15. Kingman JFC: On the genealogy of large populations. J Appl Proba 1982, 19A: 27–43.View Article
          16. Slatkin M, Hudson RR: Pairwise comparisons of mitochondrial DNA sequences in stable and exponentially growing populations. Genetics 1991, 129: 555–562.PubMed
          17. Harpending H, Rogers A: Genetic perspectives on human origins and differentiation. Proc Nat Acad Sci USA 2000, 96: 10597–10602.
          18. Rogers AR, Harpending H: Population growth makes waves in the distribution of pairwise genetic divergences. Mol Biol Evol 1992, 9: 552–569.PubMed
          19. Ray N, Currat M, Excoffier L: Intra-deme molecular diversity in spatially expanding populations. Mol Biol Evol 2003, 20: 76–86.View ArticlePubMed
          20. Excoffier L: Patterns of DNA sequence diversity and genetic structure after a range expansion: lessons from the infinite-island model. Mol Ecol 2004, 13: 853–864.View ArticlePubMed
          21. Hewitt GM: Some genetic consequences of ice ages, and their role in divergence and speciation. Biol J Linn Soc 1996, 58: 247–276.
          22. Taberlet P, Fumagalli L, Wust-Saucy AG, Cosson JF: Comparative phylogeography and postglacial colonization routes in Europe. Mol Ecol 1998, 7: 453–464.View ArticlePubMed
          23. Hewitt GM: The genetic legacy of the Quaternary ice ages. Nature 2000, 405: 907–913.View ArticlePubMed
          24. Hewitt GM: Speciation, hybrid zones and phylogeography or seeing genes in space and time. Mol Ecol 2001, 10: 537–549.View ArticlePubMed
          25. Throckmorton LH: The Genetics and Biology of Drosophila Academic Press: New York 1982., 3b:
          26. Spicer GS: Molecular evolution and phylogeny of the Drosophila virilis species group as inferred by two-dimensional electrophoresis. Journal of Molecular Evolution 1991, 33: 379–394.View ArticlePubMed
          27. Spicer GS: Reevaluation of the phylogeny of the Drosophila virilis . Annual Entomology Society of America 1992, 85: 11–25.
          28. Nurmisky DI, Moriyama EN, Lozovskaya ER, Hartl DL: Molecular phylogeny and genome evolution in the Drosophila virilis species group: duplications of the alcohol dehydrogenase gene. Molecular Biology and Evolution 1996, 13 (1) : 132–149.
          29. Spicer GS, Bell CD: Molecular phylogeny of the Drosophila virilis species group (Diptera: Drosophilidae) inferred from mitochondrial 12S and 16S ribosomal RNA genes. Annual Entomology Society of America 2002, 95: 156–161.View Article
          30. Adrianov BV, Sorokina SY, Gorelova TV, Mitrofanov VG: Mitochondrial DNA polymorphism in natural populations of the Drosophila virilis species group. Russian Journal of genetics 2003, 39 (6) : 762–768.
          31. Ostrega MS: Restriction endonuclease analysis of the relatedness of D. montana and D. virilis lines. Drosophila Inf Serv 1985, 61: 132–133.
          32. Lezin GT, Makarova KW, Velikodvorskaya VV, Zelentsova ES, Kechumyan RR, Kidwell MG, Koonin E, Evgen'ev MB: The structure and evolutionary role of the Penelope element in the Drosophila virilis species group. Mol Biol 2001, 35 (5) : 805–815.View Article
          33. Saarikettu M, Liimatainen JO, Hoikkala A: Intraspecific variation in mating behaviour does not cause sexual isolation between Drosophila virilis strains. Anim Behav 2005, 70: 417–426.View Article
          34. Hoikkala A: The importance of different courtship stimuli in the mating behaviour of European species of the Drosophila virilis group. Ann Zool Fennici 1988, 25: 257–263.
          35. Shilova VY, garbuz DG, Evgen'ev MB, Zatsepina OG: Small heat shock proteins and adaptation of various Drosophila species to hyperthermia. Molecular Biology 2006, 40 (2) : 271–276.View Article
          36. David JR, Kitagawa O: Possible similarities in ethanol tolerance and latitudinal variations between Drosophila virilis and D. melanogaster . Jpn J Genet 1982, 57: 89–95.View Article
          37. Gibert P, Moreteau B, Pétavy G, Karan D, David JR: Chill-coma tolerance, a major climatic adaptation among Drosophila species. Evolution 2001, 55: 1063–1068.View ArticlePubMed
          38. Rosenberg NA, Nordborg M: Genealogical trees, coalescent teory and the analysis of genetic polymorphisms. Nature Reviews in Genetics 2002, 3 (5) : 380–390.View Article
          39. Excoffier L, Langaney A: Origin and differentiation of human mitochondrial DNA. Am J Hum Genet 1989, 44: 73–85.PubMed
          40. Vieira J, Charlesworth B: X chromosome DNA variation in Drosophila virilis . Proc R Soc Lond B 1999, 266 (1431) : 1905–1912.View Article
          41. Wang B, Park J, Watabe H, Gao J, Xiangyu J, Aotsuka T, Chen H, Zhang Y: Molecular phylogeny of the Drosophila virilis section (Diptera: Drosophilidae) based on mitochondrial and nuclear sequences. Mol Phyl Evol 2006, 40 (2) : 484–500. Epub 2006 Mar 28.View Article
          42. Huttunen S, Aspi J, Schlötterer, Routtu J, Hoikkala A: Variation in male courtship song traits in Drosophila virilis : the effects of selection and drift on song divergence at the intraspecific level. Behavior Genetics 2008, 38 (1) : 82–92. Epub 2007 Nov 8.View ArticlePubMed
          43. McAllister B: Sequence differentiation associated with an inversion on the Neo-X chromosome of Drosophila americana . Genetics 2003, 165: 1317–1328.PubMed
          44. Caletka BC, McAllister BF: A genealogical view of chromosomal evolution and species delimitation in the Drosophila virilis species subgroup. Mol Phyl Evol 2004, 33 (3) : 664–670.View Article
          45. McAllister B, Evans AL: Increased nucleotide diversity with transient Y linkage in Drosophila americana . PloS one 2006, 1: e112-.View ArticlePubMed
          46. Schäfer MA, Orsini L, McAllister B, Schlötterer C: Patterns of microsatellite varition through a transition zone of a chromosomal cline in Drosophila americana . Heredity 2006, 97: 291–295.View ArticlePubMed
          47. Mirol PM, Schäfer MA, Orsini L, Routtu J, Schlötterer C, Hoikkala A, Butlin R: Phylogeographic patterns in Drosophila montana . Mol Ecol 2006, 16 (5) : 1085–1097.View Article
          48. Sambrook J, Fritsch EF, Maniatis T: Molecular cloning: a laboratory manual 2 Edition Cold spring Harbor Laboratory Press: New York 1989.
          49. Liu H, Beckenbach AT: Evolution of the mitochondrial cytochrome oxidase II gene among ten orders of insects. Mol Phylog Evol 1992, 41: 31–52.
          50. 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 Entom Soc Am 1994, 87: 651–701.
          51. de Brito RA, Manfrin MH, Sene FM: Mitochondrial DNA phylogeography of Brazilian populations of Drosophila buzzatii . Genet Mol Biol 2002, 25: 161–171.View Article
          52. Higgins D, Bleasby A, Fuchs R: Clustal V: Improved software for multiple sequence alignment. Comp Appl Biosc 1992, 8 (2) : 189–191.PubMed
          53. Farris JS, Kallersjo M, Kluge AG, Bult C: Testing significance of incongruence. Cladistics 1995, 10: 315–319.View Article
          54. Schneider S, Roessli D, Excoffier L: Arlequin ver. 2.000: A software for population genetics data analysis. Genetics and Biometry Laboratory, University of Geneva, Switzerland 2000.
          55. Tajima F: Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics 1989, 123: 585–595.PubMed
          56. Fu YX, Li WH: Statistical tests of neutrality of mutations. Genetics 1993, 133: 693–709.PubMed
          57. Kuhner MK, Yamato J, Felsenstein J: Maximum likelihood estimation of population growth rates based on coalescent. Genetics 1998, 149: 429–434.PubMed
          58. Posada D, Crandall KA: Modeltest: testing the model of DNA substitution. Bioinformatics 1998, 14: 817–818.View ArticlePubMed
          59. Pybus OG, Rambaut A, Harvey PH: An integrated framework for the inference of viral population history from reconstructed genealogies. Genetics 2000, 155: 1429–1437.PubMed
          60. Clement M, Posada D, Crandall KA: TCS: a computer program to estimate gene genealogies. Mol Ecol 2000, 9: 1657–1659.View ArticlePubMed

          Copyright

          © Mirol et al. 2008

          This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://​creativecommons.​org/​licenses/​by/​2.​0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

          Advertisement