<title>Introduction</title> <p>Pleistocene glacial cycles cause repeated range expansions and contractions with important present-day demographic consequences for temporal and boreal populations and species in the Northern Hemisphere ##REF##10879524##[1]##–##UREF##1##[3]##. When populations are separated and isolated in glacial refugia for a sufficient period of time they may evolve pre- or post-zygotic reproductive barriers and become reproductively isolated sister taxa ##UREF##2##[4]##–##REF##17363673##[6]##. Following postglacial range expansions such taxa may, depending on the degree of diversification, form hybrid zones with at least some interbreeding and introgression, and reinforcement processes may complete the speciation event ##UREF##4##[7]##–##UREF##5##[9]##. Assuming a constant mitochondrial molecular clock whereby substitutions take place at a rate of 1% per Myr, genetic distances between sister taxa in birds at higher latitudes translate to divergence times of approximately 0.1–3.0 Myr. As formulated in the “Pleistocene speciation hypothesis”, a common view is that these speciation events took place in isolated refugia over one to several full glacial cycles ##UREF##2##[4]##–##REF##17363673##[6]##.</p> <p>An alternative hypothesis, referred to as the “postglacial expansion hypothesis”, is that speciation takes place during expansion rather than during isolation ##REF##28564874##[10]##–##REF##17725978##[13]##. During rapid postglacial range expansions advancing populations encounter a wide variety of unoccupied habitats with varying selection regimes, which could drive diversification in comparatively short periods of time. Since this scenario does not invoke long-term isolation in refugia, the postglacial expansion hypothesis poses an alternative to traditional models of Pleistocene speciation. In support of the postglacial expansion hypothesis, recent divergence with extensive morphological variation has, for example, been documented in several bird species ##REF##28564874##[10]##–##REF##17725978##[13]##.</p> <p>Here we use molecular genetic data to evaluate postglacial colonization patterns in a bird with only slight phenotypic variation, the great reed warbler (<italic>Acrocephalus arundinaceus</italic>). This is a socially polygynous long-distance migratory passerine bird that breeds in reed marsh habitats in most of the Palaearctic between latitudes 35° and 60° N, east to North-Western Mongolia ##UREF##6##[14]##–##UREF##8##[16]##. Great reed warblers in the eastern parts of the range have slightly paler plumage than western birds and sometimes two subspecies are recognised, <italic>A. a. arundinaceus</italic> in the western parts, and <italic>A. a. zarudnyi</italic> in the eastern parts of the range ##UREF##6##[14]##, ##UREF##9##[17]##. In a previous study, Bensch and Hasselquist ##UREF##10##[18]## looked into the phylogeography of the species in the western part (Western and Central Europe) of the breeding range. Interestingly, they found that great reed warblers had two main mitochodria DNA (mtDNA) haplotype clades (called A and B), where clade A occurred in all populations with declining frequencies towards east, whereas clade B was absent from Western Europe ##UREF##10##[18]##. Based on these patterns it was hypothesised that birds with clade A haplotypes expanded from a glacial refugium in western parts and clade B from a refugium in the eastern parts of the species' breeding range ##UREF##10##[18]##.</p> <p>In the present study we reveal a more complex colonization history of the great reed warbler (i) by using mtDNA sequence data from a total of 15 populations distributed over most of the species' Eurasian breeding range now including also Eastern Europe, the Middle East and Western and Central Asia, and (ii) by using recently developed statistical tools that make it possible to test more sophisticated hypotheses regarding the divergence and expansion, e.g. inferring asymmetric gene flow between populations and time since divergence ##REF##15238526##[19]##, ##REF##17301231##[20]##. We document independent postglacial expansion events from two glacial refugia, where birds from a refugium in Europe (carrying clade A haplotypes) colonised not only Europe, but also Asia, before birds from a refugium in the Middle East (carrying clade B haplotypes) expanded north in Central and Eastern Europe. The consequences of repeated range expansions and contractions in general, and their role in driving diversification and speciation, are discussed.</p>
<title>Materials and Methods</title> <title>Study populations and molecular methods</title> <p>Samples of great reed warblers were obtained from 15 breeding populations throughout the species' breeding range: from Spain in west to Kazakhstan in east; and from Iran in south to Sweden in north (##TAB##0##Table 1##). This means that we have covered a substantial part of the breeding range, except from the most north-eastern and eastern parts of the range (Southern Siberia and North-Western Mongolia; ##UREF##6##[14]##, ##UREF##9##[17]##).</p> <p>For all individuals, except for those in the Netherlands, DNA was extracted from blood (5–50 µL stored in SET-buffer) by a standard phenol-chloroform extraction protocol. The samples from the Netherlands consisted of dried contour feathers plucked from nestlings (different nests) and DNA was extracted by using 5% Chelex 100 (Biorad; for references, see ##UREF##10##[18]##).</p> <p>The mitochondrial control region II in great reed warblers (GenBank: AF111791) was amplified and sequenced with the primers BCML4 (<named-content content-type="gene">5′-TTCACAGATACAAATGCTTGGG- 3′</named-content>) and FTPH3 (<named-content content-type="gene">5′-AAGGCTGGGAGAGTTGTTGA- 3′</named-content>) following the procedures in Bensch and Hasselquist ##UREF##10##[18]##. These two primers give a product of 577 base-pairs (bps) and previous analyses have confirmed that 494 bps of this fragment are from the control region II and the 81 bps flanking the control region (3′) from the tRNA<sup>Phe</sup> gene ##UREF##10##[18]##. Maternal inheritance of haplotypes in families and absence of double base calling ##UREF##10##[18]## confirm the mitochondrial (and not NUMT) origin of this fragment.</p> <p>Polymerase chain reactions (PCRs) were performed in volumes of 25 µL and included 10–50 ng of total genomic DNA, 0.125 mM of each nucleotide, 1.5 mM MgCl<sub>2</sub>, 0.6 µM of each primer and 0.5 U AmpliTaq polymerase. The PCRs were run using the following conditions: 30 s at 94°C, 30 s at 50°C, 30 s at 72°C (35 cycles). Before the cyclic reactions the samples were incubated at 94°C for 2 min, and after completion at 72°C for 10 min. The PCR product was precipitated (NH<sub>4</sub>Ac and ethanol) and then dissolved in 20 µL of water; 2–4 µL was then used for sequencing (BigDye sequencing kit; Applied Biosystems) in an ABI Prism 3100 capillary sequencer (Applied Biosystems).</p> <p>The data set of 93 control region sequences from six populations gathered and analysed in Bensch and Hasselquist ##UREF##10##[18]## was increased to a total number of 281 sequences from 15 populations (##TAB##0##Table 1##).</p> <title>Statistical methods</title> <p>Basic population statistics were calculated using the program D<sc>na</sc>SP 4.10.9 ##REF##14668244##[49]##: haplotype diversity (<italic>H</italic>) and nucleotide diversity (<italic>π</italic>) were calculated according to Nei ##UREF##23##[50]##.</p> <p>Evolutionary relationships between haplotypes were assessed by the Neighbour-Joining method with MEGA 3.1 ##REF##15260895##[51]##. We used Tamura and Nei's ##REF##8336541##[52]## distance measure and a gamma correction parameter alpha of 0.04 to account for among-site variation in evolutionary rate ##REF##21237881##[53]##, as had been estimated previously with P<sc>uzzle</sc>
##UREF##24##[54]## (see ##UREF##10##[18]##). The phylogeny was tested with the Bootstrap procedure (10,000 replications). We also conducted a phylogenetic analysis using the program M<sc>r</sc>B<sc>ayes</sc> 3 ##REF##12912839##[55]##, implementing a GTR+I+G model of molecular evolution as first selected by the program M<sc>r</sc>M<sc>odeltest2</sc>
##UREF##25##[56]##. The Bayesian phylogeny was obtained by using four heated and one cold MCMC chain, which was sampled every 200 generations over 20 million generations generating 100,000 trees. The first 25% of the trees were discarded as the burn-in period, and the remaining 75,000 trees were used to construct a majority consensus tree. The result from the Bayesian analysis was qualitatively similar to that of the Neighbour-Joining method, with a strong support (100%) for the branching of the two major clades (data not shown; cf. ##SUPPL##2##Figure S1##).</p> <p>We calculated the degree of genetic differentiation measured as <italic>F</italic>
<sub>ST</sub> (following ##REF##1427045##[57]##, equation 3) between pairs of populations in D<sc>na</sc>SP. Neighbour-Joining population trees were calculated in MEGA using these pair-wise <italic>F</italic>
<sub>ST</sub>-values. Minimum spanning networks were drawn in N<sc>etwork</sc> 4.2.0.1 (<ext-link ext-link-type="uri" xlink:href="http://www.fluxus-technology.com/">http://www.fluxus-technology.com/</ext-link>).</p> <p>The average number of nucleotide substitutions per site between the clades, <italic>D</italic>
<sub>XY</sub>, and the number of net nucleotide substitutions per site between clades, <italic>D</italic>
<sub>A</sub>, were calculated in D<sc>na</sc>SP <sc>(</sc>following ##UREF##23##[50]##, equation 10.20 and 10.21, respectively). Standard errors of <italic>D</italic>
<sub>XY</sub> and <italic>D</italic>
<sub>A</sub> were calculated in MEGA using the Bootstrap procedure with 1000 replicates. These parameters were used to estimate the time since divergence of the two clades with the expressions <italic>t</italic> = (<italic>D</italic>
<sub>XY</sub>/2<italic>μ</italic>)×<italic>g</italic>, and <italic>t</italic> = (<italic>D</italic>
<sub>A</sub>/2<italic>μ</italic>)×<italic>g</italic>, respectively, where <italic>μ</italic> is the rate of substitution per site Myr<sup>−1</sup> and <italic>g</italic> is the generation time (2 years; ##UREF##10##[18]##). The rate of molecular evolution of the avian control region is uncertain and we use a value of 15% but also provide the results for a range of 10–20% per site per million years (s/s/Myr) ##REF##17725978##[13]##, ##REF##28568852##[21]##–##UREF##12##[24]##.</p> <p>We calculated mismatch distributions in D<sc>na</sc>SP. The observed mismatch distribution was tested against expected values in a stable population (i.e. population with constant population size; following ##REF##1145509##[58]##), and in a growing or declining population (following ##REF##1316531##[59]##, equation 4), with χ<sup>2</sup>-tests (pooling classes to achieve a minimum of 5 expected counts in each class). We estimated the parameter <italic>τ</italic>, i.e. the date of the expansion measured in units of mutational time, which can be used to estimate the time since expansion with the expression <italic>t</italic> = (<italic>τ</italic>/2<italic>μk</italic>)×<italic>g</italic>, where <italic>μ</italic> and <italic>g</italic> are as above and <italic>k</italic> is the sequence length (i.e. 494 bp).</p> <p>We also inferred populations sizes, asymmetric gene flow and time since divergence between pairs of populations with demographic population genetic models implemented in the program IM<sc>a</sc>, Isolation with Migration, analytic, version 7/13/2007 ##REF##15238526##[19]##, ##REF##17301231##[20]##. We did preliminary runs with different settings to determine a suitable multi-dimensional parameter space; we started with default values and increased them gradually to achieve output parameters within the parameter space. In the final runs, we used a value of 100 for the scalars of <italic>θ</italic> of all populations (i.e. population 1, population 2 and the ancestral population); a value of 100 as the maximum migration rate between populations (from population 1 to population 2, and from population 2 to population 1, respectively); and 10 for the maximum time of population splitting. For each pair of populations we used a ‘burn-in’ period of 1×10<sup>6</sup> steps and a sampling period of 1×10<sup>7</sup> steps. The output parameters (<italic>θ</italic>
<sub>1</sub>, <italic>θ</italic>
<sub>2</sub>, <italic>θ</italic>
<sub>A</sub>, <italic>m</italic>
<sub>1</sub>, <italic>m</italic>
<sub>2</sub>, <italic>t</italic>) were expressed in demographic units as follow: effective population size, <italic>N</italic>
<sub>E</sub> = <italic>θ</italic>/4<italic>μkg</italic>; number of migrants per generation, <italic>M</italic> = 2<italic>N</italic>
<sub>E</sub>×(<italic>m</italic>×2<italic>μk</italic>) = <italic>θ</italic>×<italic>m</italic>/2<italic>g</italic>; divergence time in years, <italic>T</italic> = <italic>tg</italic>/2<italic>μk</italic>; where <italic>μ</italic>, <italic>k</italic> and <italic>g</italic> are as above.</p>
<title>Results</title> <title>Haplotype and nucleotide diversity</title> <p>We sequenced 494 bps of the control region II in 281 great reed warblers from 15 Eurasian populations (##TAB##0##Table 1##). This region contained 39 variable sites which defined 58 different haplotypes separated by up to 2.0%. Haplotype and nucleotide diversity for all samples were 0.913 (ranged between 0.286 for Iran and 0.947 for Hungary) and 0.00735 (ranged between 0.00117 for Iran and 0.00912 for the Czech Republic), respectively (##TAB##0##Table 1##).</p> <title>Haplotype relationship and population divergence</title> <p>A Neighbour-Joining analysis of the 281 sequences revealed two major clades (A and B) supported by a bootstrap value of 75% (##SUPPL##2##Figure S1##), and a minimum-spanning network (MSN) was drawn to show the relationship between haplotype within and between clades (##FIG##0##Figure 1##).</p> <p>The average number of nucleotide substitutions per site between the clades, <italic>D</italic>
<sub>XY</sub>, was 0.013 (±0.004 SE) and the number of net nucleotide substitutions per site between clades, <italic>D</italic>
<sub>A</sub>, was 0.010 (±0.004 SE). Assuming a mutation rate of 15% (range 10–20%) per site Myr<sup>−1</sup> for the control region ##REF##17725978##[13]##, ##REF##28568852##[21]##–##UREF##12##[24]##, a constant molecular clock and a generation time of 2 years ##UREF##10##[18]##, we obtain a divergence time for the two clades of 87 (65–141) and 65 (49–98) kyr BP for <italic>D</italic>
<sub>XY</sub> and <italic>D</italic>
<sub>A</sub>, respectively.</p> <p>In ##FIG##1##Figure 2##, we have plotted the frequencies of the two clades in each population. Clade A is more widely dispersed and dominates in most parts of the breeding range except in some populations in the Central Europe where it occurs in intermediate frequencies, and in the Middle East where it is almost absent (##FIG##1##Figure 2##). Clade B is more restricted geographically, and dominates in the Middle East with declining frequencies towards north in South-Eastern and Central Europe (##FIG##1##Figure 2##).</p> <p>The differentiation between pairs of populations measured as <italic>F</italic>
<sub>ST</sub> varied substantially and reached its highest value between the Netherlands and Iran (<italic>F</italic>
<sub>ST</sub> = 0.871; <italic>F</italic>
<sub>ST</sub>-values are given in ##SUPPL##0##Table S1a##). To illustrate the genetic distance between populations we generated a Neighbour-Joining (NJ) population tree from these pair-wise distances (##FIG##2##Figure 3a##). In this tree, three main branches could be seen: one leading to Spain and the Netherlands, another to Kazakhstan and the last one to Turkey and Iran. Thus, the main features of this tree correspond to the geographical distribution of clade A and B, although the two separate branches leading to Spain and the Netherlands, and Kazakhstan, respectively, suggest that there is also detectable population divergence within clade A.</p> <p>To evaluate the latter, we performed separate analyses for sequences belonging to clade A and B. The sample size of each clade was low in some populations and therefore we pooled data from some closely located populations in these analyses. We used eight populations for clade A (##SUPPL##0##Table S1b##) and five populations for clade B (##SUPPL##0##Table S1c##). There were several high <italic>F</italic>
<sub>ST</sub>-values for clade A, with the highest <italic>F</italic>
<sub>ST</sub> of 0.349 for Spain and Kazakhstan (##SUPPL##0##Table S1b##). The highest <italic>F</italic>
<sub>ST</sub>-value for the clade B was 0.238 between Belarus/Ukraine/Russia and Turkey/Iran (##SUPPL##0##Table S1c##). In the corresponding NJ population tree for clade A, there was a long branch separating Kazakhstan from the other populations via Russia, and an additional shorter branch leading to Spain (##FIG##2##Figure 3b##). The NJ population tree for clade B showed overall shorter branches (##FIG##2##Figure 3c##).</p> <p>A closer inspection of the NJ sequence tree reveals that there is a large relatively uniform cluster of 26 sequences from Kazakhstan and Russia in the central part of the tree (indicated with a dotted line in ##SUPPL##2##Figure S1##; shown in grey in ##FIG##0##Figure 1##), which explains the long branch leading to Kazakhstan in the NJ population trees (##FIG##2##Figure 3a,b##).</p> <title>Mismatch distributions and clade specific analyses of gene flow</title> <p>When all haplotypes were tested simultaneously there were two clear peaks in the mismatch distribution at 1–2 bp and 6–7 bp, respectively (##FIG##3##Figure 4a##). This distribution did not correspond to what would be expected in a stable population (χ<sup>2</sup> = 67.0, df = 10, p&lt;0.001) or an expanding population (χ<sup>2</sup> = 48.5, df = 9, p&lt;0.001). This supports the view that there have been expansions from two glacial refugia populations. When the two clades were tested separately, the mismatch distributions for both clades differed significantly from what would be expected in a stable population (clade A: χ<sup>2</sup> = 49.2, df = 6, p&lt;0.001; clade B: χ<sup>2</sup> = 23.5, df = 4, p&lt;0.001), whereas the distributions did not deviate from the expected in an expanding population (clade A: χ<sup>2</sup> = 2.6, df = 5, p = 0.77; clade B: χ<sup>2</sup> = 3.5, df = 3, p = 0.31) (##FIG##3##Figure 4b,c##). The value <italic>τ</italic> was 1.734 for clade A and 1.467 for clade B, which translate to dates of expansion of 23 (18–35) and 20 (15–30) kyr BP, respectively.</p> <p>We inferred the colonisation history of each clade separately by using isolation with migration models implemented in the program IM<sc>a</sc>
##REF##15238526##[19]##, ##REF##17301231##[20]##. To increase sample size in each group, we pooled populations as in the previous analyses. For clade A, the present-day effective population sizes were approximately 50,000, the ancestral population size approximately 10,000 and the time since divergence between pairs of populations approximately 4 kyr (##SUPPL##1##Table S2##). The gene flow varied substantially between populations: there was little gene flow between Spain and neighbouring populations (<italic>M</italic>&lt;8 migrants per generation), strong bi-directional gene flow between almost all central European populations (<italic>M</italic>&gt;100 for several populations), a tendency for an asymmetric gene flow from eastern Europe to Russia (<italic>M</italic> = 8.3 and 5.3, respectively) and pronounced asymmetric gene flow from Russia to Kazakhstan (<italic>M</italic> = 32.6 and 2.2, respectively; ##FIG##4##Figure 5a##; ##SUPPL##1##Table S2##). Clade B had present-day and ancestral population sizes of approximately 10,000 and the time since divergence between pairs of populations was approximately 1.5 kyr (##SUPPL##1##Table S2##). The gene flow between populations was overall much lower for clade B (several <italic>M</italic>≈1; max <italic>M</italic> = 13.7) than for clade A (several <italic>M</italic>&gt;100), and, if anything, there was weak asymmetric gene flow from southern to northern populations (##FIG##4##Figure 5b##; ##SUPPL##1##Table S2##). In a model where all clade A haplotypes were pooled in one population and all B haplotypes in another population, the isolation with migration models produced an estimate of the time since divergence of 18 kyr BP, an ancestral effective population size of 13,000, and present-day effective population sizes of 95,000 and 25,000 for the clade A and B, respectively (##SUPPL##1##Table S2##).</p>
<title>Discussion</title> <title>Evidence for two independent postglacial expansions</title> <p>Our analyses point to the fact that the great reed warbler has expanded along two independent routes from two separated glacial refugia. There are two major clades in the mtDNA sequence tree and a clear geographical pattern in the frequency of these clades still exists among present-day populations. The analyses of mismatch distributions support that two independent expansions have occurred and reject the possibility of one or two large panmictic populations.</p> <p>The most parsimonious interpretation is that clade B was isolated in a refugium located somewhere in the Middle East and this view gains some support from the analyses of gene flow suggesting low levels of asymmetric gene flow towards north. The location of clade A's refugium is more difficult to track down, because the clade is presently widely dispersed and seems to represent a relatively old expansion. However, two sources of information indicate that the refugium was located in the western parts of the breeding range. First, the gene flow analyses show a slight asymmetric gene flow from Eastern Europe to Western Asia and pronounced asymmetric gene flow from Western Asia to Central Asia. Secondly, a closer inspection of the NJ sequence tree, and the minimum-spanning network, reveals that there is a large, relatively uniform cluster of 26 sequences from Western and Central Asia in the central part of the tree, and no indication that these haplotypes have dispersed to Europe. We conclude that clade A's refugium was located somewhere in South-Western Palaearctic, with a colonisation route that went north of the Black sea to Central Asia. There are several potential, well-known refugia in that region including the Iberian peninsula, Italy and the Balkans, as well as North-Western Africa ##REF##10879524##[1]##, ##UREF##13##[25]##.</p> <p>The absence of gene flow between the Middle East and Central Asia probably reflects geographical expansion barriers in form of large mountain ridges located between the Black Sea and the Caspian Sea as well as to the east of the Caspian Sea. The great reed warbler is a long-distant migrant where all individuals migrate to Africa south of the Sahara ##UREF##6##[14]##, ##UREF##14##[26]##. The migratory routes and overwintering strategies of great reed warblers from Central Asia are unknown, but if these populations descend from European ancestors as our data suggest, their migratory habits might be more similar to those of the European populations, than those of populations in the Middle East ##UREF##6##[14]##, ##UREF##10##[18]##.</p> <p>The great reed warbler is highly philopatric in some populations, yet has a relatively high dispersal potential ##REF##28547254##[27]##–##REF##12675819##[29]##, and it is possible that it is capable of rapidly expanding its range during favourable conditions. Rapid expansions have been documented in a few bird species, e.g. the collared dove (<italic>Streptopelia decaocto</italic>) has colonised much of Europe in only 150 years ##UREF##16##[30]##. Our estimates of the divergence time of the two haplotype clades in great reed warblers range from 18 kyr BP with the IMa method to 65–87 kyr BP based the degree of sequence divergence. These values are within the last glaciation period, which reached its maximum approximately 20 kyr BP ##REF##12913123##[31]##, ##REF##15101575##[32]##. Then it took several thousands of years before the new unoccupied habitat became available. Probably the pristine habitat has been available only during the last few thousand years (maybe up to 10 kyr BP), which sets the maximum date for the expansion. Thus, the estimated expansion dates from the mismatch distribution analyses (i.e. 24 and 20 kyr BP for clade A and clade B, respectively) are too high. The divergence dates between present-day populations estimated with the IM<sc>a</sc> method seem more plausible: 4 kyr BP for clade A and 1.5 kyr BP for clade B. Still, these estimates are sensitive to mutation rates and generation times, and assume a constant molecular clock and migration-drift equilibrium, thus should be interpreted cautiously. Although it may be very difficult to date the postglacial expansion in the great reed warbler, our data strongly suggest that the expansion of clade A preceded that of clade B. The level of gene flow is higher between populations for clade A, and clade A has colonised a much wider area and is now present all the way from Western Europe to Central Asia.</p> <title>Consequences of repeated range expansions and contractions</title> <p>Speciation may be the most spectacular consequence of repeated glacial cycles, but there are also other important consequences of glacial population bottlenecking and postglacial expansions. For instance, geographical barriers may result in limits to present-day range distributions and failures to colonize seemingly suitable habitat ##REF##10879524##[1]##, ##REF##3576198##[33]##. Thus, glacial range oscillation is one of several factors that contribute to the evolution of species ranges ##REF##10879524##[1]##, ##REF##3576198##[33]##–##REF##18811273##[35]##. The presence of dispersal barriers to the north and north-eastern of the Middle East seems to provide a plausible explanation for the delayed expansion of clade B in great reed warblers.</p> <p>Postglacial expansions may also leave more subtle traces in form of lower genetic variation in populations at the expanding range limits; referred to as the ‘leading edge hypothesis’ ##REF##10879524##[1]##, ##UREF##17##[36]##. In the present study, the admixture of the two clades makes the interpretation difficult, but it is true that several populations in the expanding range limit have low levels of haplotype and nucleotide diversity (e.g. the Netherlands and Kazakhstan). However, the two supposedly older populations in the Middle East have low diversity and consequently the leading edge hypothesis does not provide an exclusive explanation to mtDNA diversity in great reed warblers.</p> <p>A study of the breeding ecology of the great reed warbler in Sweden has been ongoing for more that two decades ##UREF##7##[15]##, ##UREF##18##[37]##, ##UREF##19##[38]##. An unexpected result from our previous studies is that molecular markers hitchhike with important fitness genes ##UREF##20##[39]##, ##REF##15154561##[40]##, which suggest that there are high levels of linkage disequilibrium in the population ##REF##9650277##[41]##, ##REF##12453232##[42]##. It is well-known that admixture may cause linkage disequilibrium ##UREF##21##[43]##, ##REF##11818140##[44]## and since admixture may follow from rapid populations expansions and contractions this is yet another important consequence of postglacial range oscillations. Whether postglacial admixture has contributed to the high level of linkage disequilibrium in the great reed warbler (##REF##28565235##[45]##, ##REF##15612297##[46]##; B. Hansson, K. Csilléry, unpublished) remains to be evaluated, but the present result is suggestive.</p> <p>The Pleistocene speciation hypothesis suggests that speciation occurs in isolated refugia over one or more full glacial cycles ##UREF##2##[4]##–##REF##17363673##[6]##. The postglacial expansion hypothesis proposes an alternative explanation namely that diversification and speciation takes place during expansion rather than during isolation ##REF##28564874##[10]##–##REF##17725978##[13]##. During postglacial expansions the advancing populations may encounter a variety of unoccupied habitats with varying selection regimes, which could drive diversification. In support of this hypothesis, recent divergence with extensive morphological variation has, for example, been documented in several bird species; redpolls (<italic>Carduelis flammea</italic>; ##REF##28564874##[10]##; but see ##UREF##22##[47]##), yellow wagtails (<italic>Motacilla flava</italic>; ##REF##12859633##[11]##), yellow-rumped warblers (<italic>Dendroica coronata</italic>; ##REF##17181728##[12]##) , willow warblers (<italic>Phylloscopus trochilus</italic>; ##REF##28565443##[48]##) and dark-eyed juncos (<italic>Junco hyemalis</italic>; ##REF##17725978##[13]##).</p> <p>As mentioned above, great reed warblers in the eastern parts of the range (east of the Black Sea; <italic>A. a. zarudnyi</italic>) have paler plumage than western birds (<italic>A. a. arundinaceus</italic>) ##UREF##6##[14]##, ##UREF##9##[17]## and our results suggest that the plumage diversification took place after the last glaciation during the easterly expansion of clade A. This gives support to the postglacial expansion hypothesis and the role of expansions in driving diversification in comparatively short time periods. However, the paler <italic>A. a. zarudnyi</italic> is also recognised in eastern parts of the Middle East ##UREF##6##[14]##, ##UREF##9##[17]##, i.e. in the region where our Iranian population with clade B haplotypes is located. This suggests a more complicated history of the plumage diversification in the great reed warbler. For instance, the plumage variation may have appeared recently and differentiated independently in each expanding clade. Alternatively, the plumage variation may be of old origin and could have survived in both refugia, with the plumage differentiation following during the expansion of each clade. Both these scenarios are in line with the postglacial expansion hypothesis.</p> <p>Our data show that the refugia populations became separated during the last glacial period and that the postglacial expansions have resulted in substantial admixture in Central and Eastern Europe. However, there is no indication of any (strong) pre- or post-zygotic isolation mechanisms between the great reed warbler clades (##REF##28565235##[45]##, ##REF##15612297##[46]##; S. Bensch, B. Hansson, D. Hasselquist, unpublished). This is in line with the Pleistocene speciation hypothesis invoking that much longer periods of time in isolation are needed for speciation to occur ##UREF##2##[4]##–##REF##17363673##[6]##.</p>
<p>Conceived and designed the experiments: BH DH SB. Analyzed the data: BH SB. Wrote the paper: BH. Organized field work and data collection: SB PZ DH BH. Performed lab work and evaluated sequences: SB MT.</p> <p>Pleistocene glacial cycles play a major role in diversification and speciation, although the relative importance of isolation and expansion in driving diversification remains debated. We analysed mitochondrial DNA sequence data from 15 great reed warbler (<italic>Acrocephalus arundinaceus</italic>) populations distributed over the vast Eurasian breeding range of the species, and revealed unexpected postglacial expansion patterns from two glacial refugia. There were 58 different haplotypes forming two major clades, A and B. Clade A dominated in Western Europe with declining frequencies towards Eastern Europe and the Middle East, but showed a surprising increase in frequency in Western and Central Asia. Clade B dominated in the Middle East, with declining frequencies towards north in Central and Eastern Europe and was absent from Western Europe and Central Asia. A parsimonious explanation for these patterns is independent postglacial expansions from two isolated refugia, and mismatch distribution analyses confirmed this suggestion. Gene flow analyses showed that clade A colonised both Europe and Asia from a refugium in Europe, and that clade B expanded much later and colonised parts of Europe from a refugium in the Middle East. Great reed warblers in the eastern parts of the range have slightly paler plumage than western birds (sometimes treated as separate subspecies; <italic>A. a. zarudnyi</italic> and <italic>A. a. arundinaceus</italic>, respectively) and our results suggest that the plumage diversification took place during the easterly expansion of clade A. This supports the postglacial expansion hypothesis proposing that postglacial expansions drive diversification in comparatively short time periods. However, there is no indication of any (strong) reproductive isolation between clades and our data show that the refugia populations became separated during the last glaciation. This is in line with the Pleistocene speciation hypothesis invoking that much longer periods of time in isolation are needed for speciation to occur.</p>
<title>Supporting Information</title>
<p>Samples were kindly provided by R.T. Pinheiro (Spain), J. Graveland (the Netherlands), S. Fischer (Germany), P. Prochazka (the Czech Republic and Turkey), A. Poluda (Ukraine), M. Illieva (Russia), S.M. Amininasab (Iran). M. Svensson and Å. Bensch assisted with data collection in Greece and Latvia; C. Moskat in Hungary; A.K. Olsson, E. Gavrilov and A. Gavrilov in Kazakhstan. This is report no. 150 from Kvismare Bird Observatory.</p>
<fig id="pone-0002794-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002794.g001</object-id><label>Figure 1</label><caption><title>Minimum-spanning network (MSN) for 58 mtDNA control region II haplotypes in great reed warblers from 15 breeding populations.</title><p>Circles are proportional to sample size (281 sequences in total) and each line corresponds to one substitution; with the exception of thicker lines, which correspond to two substitutions. The two major haplotype clades are shown, clade A and B. Grey shading indicates a relatively large uniform cluster of sequences within clade A only found in Russia and Kazakhstan (see text and ##SUPPL##2##Figure S1## for details).</p></caption></fig> <fig id="pone-0002794-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002794.g002</object-id><label>Figure 2</label><caption><title>The frequencies of the two clades in each of the 15 great reed warbler breeding populations (clade A is in white and clade B in black).</title><p>Breeding distribution is indicated (dotted line).</p></caption></fig> <fig id="pone-0002794-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002794.g003</object-id><label>Figure 3</label><caption><title>Neighbour-Joining (NJ) population trees generated from pair-wise population distances (<italic>F</italic>
<sub>ST</sub>-values) for (a) all haplotypes, (b) clade A haplotypes, and (c) clade B haplotypes.</title><p>Because of low samples sizes in some populations, we pooled closely located populations in the clade specific analyses (b and c).</p></caption></fig> <fig id="pone-0002794-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002794.g004</object-id><label>Figure 4</label><caption><title>Mismatch distributions for (a) all haplotypes, (b) clade A haplotypes and (c) clade B haplotypes.</title><p>Shown are the observed distribution (lines), the expected distribution of a population with constant size (dotted) and the expected distribution of an expanding population (dashed).</p></caption></fig> <fig id="pone-0002794-g005" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002794.g005</object-id><label>Figure 5</label><caption><title>Schematic illustration of the gene flow between populations estimated with the IM<sc>a</sc> method for (a) clade A and (b) clade B.</title><p>Size of arrows are proportional to the gene flow (the largest arrow gives M≥100 migrants per generation; see ##SUPPL##1##Table S2## for details). Labels refer to population (SPA–Spain; NET–the Netherlands; GER–Germany; SWE–Sweden; LAT–Latvia; BEL–Belarus; UKR–Ukraine; CZE–the Czech Republic; HUN–Hungary; BUL–Bulgaria; GRE–Greece; RUS–Russia; KAZ–Kazakhstan; TUR–Turkey; IRA–Iran).</p></caption></fig>
<table-wrap id="pone-0002794-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002794.t001</object-id><label>Table 1</label><caption><title>Sampling localities, sampling period, number of sequences (<italic>n</italic>), number of haplotypes (<italic>nH</italic>), haplotype diversity (<italic>H</italic>) and nucleotide diversity (<italic>π</italic>).</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Country</td><td align="left" rowspan="1" colspan="1">Site</td><td align="left" rowspan="1" colspan="1">Coordinates</td><td align="left" rowspan="1" colspan="1">Sampling year</td><td align="left" rowspan="1" colspan="1">
<italic>n</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>nH</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>H</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>π</italic>
</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">Spain</td><td align="left" rowspan="1" colspan="1">Hondo Natural Park, Alicante</td><td align="left" rowspan="1" colspan="1">38°12′N, 0°42′W</td><td align="left" rowspan="1" colspan="1">1996</td><td align="left" rowspan="1" colspan="1">11</td><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1">0.709</td><td align="left" rowspan="1" colspan="1">0.00202</td></tr><tr><td align="left" rowspan="1" colspan="1">The Netherlands</td><td align="left" rowspan="1" colspan="1">Zwarte Meer,Weerribben</td><td align="left" rowspan="1" colspan="1">52°37′N, 5°55′E</td><td align="left" rowspan="1" colspan="1">1995</td><td align="left" rowspan="1" colspan="1">10</td><td align="left" rowspan="1" colspan="1">4</td><td align="left" rowspan="1" colspan="1">0.533</td><td align="left" rowspan="1" colspan="1">0.00205</td></tr><tr><td align="left" rowspan="1" colspan="1">Sweden</td><td align="left" rowspan="1" colspan="1">Kvismaren, Närke</td><td align="left" rowspan="1" colspan="1">59°10′N, 15°25′E</td><td align="left" rowspan="1" colspan="1">1987–1990</td><td align="left" rowspan="1" colspan="1">22</td><td align="left" rowspan="1" colspan="1">7</td><td align="left" rowspan="1" colspan="1">0.649</td><td align="left" rowspan="1" colspan="1">0.00480</td></tr><tr><td align="left" rowspan="1" colspan="1">Latvia</td><td align="left" rowspan="1" colspan="1">Engure/Kanieris, Tukums</td><td align="left" rowspan="1" colspan="1">57°07′N, 23°20′E</td><td align="left" rowspan="1" colspan="1">1992</td><td align="left" rowspan="1" colspan="1">20</td><td align="left" rowspan="1" colspan="1">11</td><td align="left" rowspan="1" colspan="1">0.916</td><td align="left" rowspan="1" colspan="1">0.00654</td></tr><tr><td align="left" rowspan="1" colspan="1">Germany</td><td align="left" rowspan="1" colspan="1">Müggelsee, Berlin</td><td align="left" rowspan="1" colspan="1">52°26′N, 13°39′E</td><td align="left" rowspan="1" colspan="1">1992,1993</td><td align="left" rowspan="1" colspan="1">19</td><td align="left" rowspan="1" colspan="1">12</td><td align="left" rowspan="1" colspan="1">0.918</td><td align="left" rowspan="1" colspan="1">0.00652</td></tr><tr><td align="left" rowspan="1" colspan="1">The Czech Republic</td><td align="left" rowspan="1" colspan="1">Hodonin fishponds, Moravia</td><td align="left" rowspan="1" colspan="1">48°54′N, 17°02′E</td><td align="left" rowspan="1" colspan="1">2006</td><td align="left" rowspan="1" colspan="1">17</td><td align="left" rowspan="1" colspan="1">11</td><td align="left" rowspan="1" colspan="1">0.934</td><td align="left" rowspan="1" colspan="1">0.00912</td></tr><tr><td align="left" rowspan="1" colspan="1">Hungary</td><td align="left" rowspan="1" colspan="1">Apaj Channels, Kiskunlachaza</td><td align="left" rowspan="1" colspan="1">47°07′N, 19°05′E</td><td align="left" rowspan="1" colspan="1">2000</td><td align="left" rowspan="1" colspan="1">20</td><td align="left" rowspan="1" colspan="1">14</td><td align="left" rowspan="1" colspan="1">0.947</td><td align="left" rowspan="1" colspan="1">0.00565</td></tr><tr><td align="left" rowspan="1" colspan="1">Belarus</td><td align="left" rowspan="1" colspan="1">Turov, Zhitkovichi</td><td align="left" rowspan="1" colspan="1">52°01′N, 27°49′E</td><td align="left" rowspan="1" colspan="1">2000</td><td align="left" rowspan="1" colspan="1">17</td><td align="left" rowspan="1" colspan="1">9</td><td align="left" rowspan="1" colspan="1">0.890</td><td align="left" rowspan="1" colspan="1">0.00737</td></tr><tr><td align="left" rowspan="1" colspan="1">Ukraine</td><td align="left" rowspan="1" colspan="1">Usovka, Poltava</td><td align="left" rowspan="1" colspan="1">50°19′N, 32°32′E</td><td align="left" rowspan="1" colspan="1">2000</td><td align="left" rowspan="1" colspan="1">18</td><td align="left" rowspan="1" colspan="1">11</td><td align="left" rowspan="1" colspan="1">0.941</td><td align="left" rowspan="1" colspan="1">0.00799</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Denisyvka, Poltava</td><td align="left" rowspan="1" colspan="1">49°53′N, 32°36′E</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Vilkovo, Odessa</td><td align="left" rowspan="1" colspan="1">45°28′N, 29°35′E</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">Russia</td><td align="left" rowspan="1" colspan="1">Steppe liman, Saratov</td><td align="left" rowspan="1" colspan="1">50°43′N, 46°27′E</td><td align="left" rowspan="1" colspan="1">2006</td><td align="left" rowspan="1" colspan="1">28</td><td align="left" rowspan="1" colspan="1">12</td><td align="left" rowspan="1" colspan="1">0.857</td><td align="left" rowspan="1" colspan="1">0.00524</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Solyanka, Saratov</td><td align="left" rowspan="1" colspan="1">50°49′N, 47°05′E</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Furmanovo, Saratov</td><td align="left" rowspan="1" colspan="1">51°38′N, 49°07′E</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">Kazakhstan</td><td align="left" rowspan="1" colspan="1">Stone Lake, Zhambyl</td><td align="left" rowspan="1" colspan="1">42°49′N, 70°56′E</td><td align="left" rowspan="1" colspan="1">2001</td><td align="left" rowspan="1" colspan="1">35</td><td align="left" rowspan="1" colspan="1">9</td><td align="left" rowspan="1" colspan="1">0.793</td><td align="left" rowspan="1" colspan="1">0.00275</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Lake Balkhash, Almaty</td><td align="left" rowspan="1" colspan="1">45°12′N, 73°59′E</td><td align="left" rowspan="1" colspan="1">2006</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">Greece</td><td align="left" rowspan="1" colspan="1">Limni Mikri Prespa, Florina</td><td align="left" rowspan="1" colspan="1">40°50′N, 21°05′E</td><td align="left" rowspan="1" colspan="1">1990</td><td align="left" rowspan="1" colspan="1">20</td><td align="left" rowspan="1" colspan="1">13</td><td align="left" rowspan="1" colspan="1">0.932</td><td align="left" rowspan="1" colspan="1">0.00906</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Mitrikou Lake, Rodopi</td><td align="left" rowspan="1" colspan="1">40°58′N, 25°17′E</td><td align="left" rowspan="1" colspan="1">2005</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">Bulgaria</td><td align="left" rowspan="1" colspan="1">Kalimok, Tutrakan</td><td align="left" rowspan="1" colspan="1">44°01′N, 26°26′E</td><td align="left" rowspan="1" colspan="1">2005–2006</td><td align="left" rowspan="1" colspan="1">20</td><td align="left" rowspan="1" colspan="1">12</td><td align="left" rowspan="1" colspan="1">0.900</td><td align="left" rowspan="1" colspan="1">0.00836</td></tr><tr><td align="left" rowspan="1" colspan="1">Turkey</td><td align="left" rowspan="1" colspan="1">Mogan Lake, Ankara</td><td align="left" rowspan="1" colspan="1">39°46′N, 32°48′E</td><td align="left" rowspan="1" colspan="1">2005</td><td align="left" rowspan="1" colspan="1">17</td><td align="left" rowspan="1" colspan="1">7</td><td align="left" rowspan="1" colspan="1">0.824</td><td align="left" rowspan="1" colspan="1">0.00374</td></tr><tr><td align="left" rowspan="1" colspan="1">Iran</td><td align="left" rowspan="1" colspan="1">Zarin Kola, Mazandaran</td><td align="left" rowspan="1" colspan="1">36°44′N, 53°00′E</td><td align="left" rowspan="1" colspan="1">2004</td><td align="left" rowspan="1" colspan="1">7</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">0.286</td><td align="left" rowspan="1" colspan="1">0.00117</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>Pooled</italic>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<italic>281</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>58</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>0.913</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>0.00735</italic>
</td></tr></tbody></table></alternatives></table-wrap>
<supplementary-material content-type="local-data" id="pone.0002794.s001"><label>Table S1</label><caption><p>F<sub>ST</sub>-values between populations using data from (a) clade A and B haplotypes; (b) clade A haplotypes; and (c) clade B haplotypes.</p><p>(0.08 MB DOC)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pone.0002794.s002"><label>Table S2</label><caption><p>Estimates of demographic parameters from two-population isolation with migration models implemented in the program IMa ##REF##15238526##[19]##, ##REF##17301231##[20]##.</p><p>(0.08 MB DOC)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pone.0002794.s003"><label>Figure S1</label><caption><p>Neighbour-Joining tree based on 281 mtDNA control region II sequences of great reed warblers from 15 breeding populations. Labels refer to population (SPA-Spain; NET-the Netherlands; GER-Germany; SWE-Sweden; LAT-Latvia; BEL-Belarus; UKR-Ukraine; CZE-the Czech Republic; HUN-Hungary; BUL-Bulgaria; GRE-Greece; RUS-Russia; KAZ-Kazakhstan; TUR-Turkey; IRA-Iran). The two major haplotype clades are shown, clade A and B. The dotted line indicates a relatively large uniform cluster of sequences within clade A only found in Russia and Kazakhstan (see text for details).</p><p>(3.30 MB TIF)</p></caption></supplementary-material>
<fn-group><fn fn-type="COI-statement"><p><bold>Competing Interests: </bold>The authors have declared that no competing interests exist.</p></fn><fn fn-type="financial-disclosure"><p><bold>Funding: </bold>The work was supported by grants from the Swedish Research Council (VR), the Swedish council for Environment, Agricultural Sciences and Spatial Planning (FORMAS), the National Swedish Environment Protection Board, the Carl Tryggers stiftelse, the Royal Swedish Academy of Science (the K. and A. Wallenbergs, the J.A. Ahlstrand, and the Hierta-Retzius Foundations), the Lunds Djurskyddsfond, the Magn Bergvalls Stiftelse, the National Geographic Society and the Crafoordska Stiftelsen.</p></fn></fn-group>
<graphic id="pone-0002794-t001-1" xlink:href="pone.0002794.t001"/> <graphic xlink:href="pone.0002794.g001"/> <graphic xlink:href="pone.0002794.g002"/> <graphic xlink:href="pone.0002794.g003"/> <graphic xlink:href="pone.0002794.g004"/> <graphic xlink:href="pone.0002794.g005"/>
<media xlink:href="pone.0002794.s001.doc"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pone.0002794.s002.doc"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pone.0002794.s003.tif"><caption><p>Click here for additional data file.</p></caption></media>
[{"label": ["2"], "element-citation": ["\n"], "surname": ["Schluter"], "given-names": ["D"], "year": ["2000"], "source": ["The ecology of adaptive radiation"], "publisher-loc": ["Oxford"], "publisher-name": ["Oxford university press"]}, {"label": ["3"], "element-citation": ["\n"], "surname": ["Coyne", "Orr"], "given-names": ["JA", "HA"], "year": ["2004"], "source": ["Speciation"], "publisher-loc": ["Sunderland, MA"], "publisher-name": ["Sinauer Associates, Inc"]}, {"label": ["4"], "element-citation": ["\n"], "surname": ["Klicka", "Zink"], "given-names": ["J", "RM"], "year": ["1997"], "article-title": ["The importance of recent ice ages in speciation: a failed paradigm."], "source": ["Science"], "volume": ["277"], "fpage": ["1666"], "lpage": ["1669"]}, {"label": ["5"], "element-citation": ["\n"], "surname": ["Johnson", "Cicero"], "given-names": ["NK", "C"], "year": ["2004"], "article-title": ["New mitochondrial DNA data affirm the importance of Pleistocene speciation in North American birds."], "source": ["Evolution Int J Org Evolution"], "volume": ["58"], "fpage": ["1122"], "lpage": ["1130"]}, {"label": ["7"], "element-citation": ["\n"], "surname": ["Spencer", "McArdle", "Lambert"], "given-names": ["HG", "BH", "DM"], "year": ["1986"], "article-title": ["A theoretical investigation of speciation by reinforcement."], "source": ["Am Nat"], "volume": ["128"], "fpage": ["241"], "lpage": ["262"]}, {"label": ["9"], "element-citation": ["\n"], "surname": ["S\u00e6tre", "Moum", "Bures", "Kr\u00e1l", "Adamjan", "Moreno"], "given-names": ["G-P", "T", "S", "M", "M", "J"], "year": ["1997"], "article-title": ["A sexually selected character displacement in flycatchers reinforces premating isolation."], "source": ["Nature"], "volume": ["387"], "fpage": ["589"], "lpage": ["592"]}, {"label": ["14"], "element-citation": ["\n"], "surname": ["Cramp"], "given-names": ["S"], "year": ["1992"], "source": ["Handbook of the birds of Europe, the Middle East and North Africa. Vol VI"], "publisher-loc": ["Oxford"], "publisher-name": ["Oxford University Press"]}, {"label": ["15"], "element-citation": ["\n"], "surname": ["Hasselquist"], "given-names": ["D"], "year": ["1998"], "article-title": ["Polygyny in the great reed warbler: a long-term study of factors contributing to male fitness."], "source": ["Ecology"], "volume": ["79"], "fpage": ["2376"], "lpage": ["2390"]}, {"label": ["16"], "element-citation": ["\n"], "surname": ["Hansson", "Bensch", "Hasselquist"], "given-names": ["B", "S", "D"], "year": ["2000"], "article-title": ["Patterns of nest predation contribute to polygyny in the great reed warbler."], "source": ["Ecology"], "volume": ["81"], "fpage": ["319"], "lpage": ["328"]}, {"label": ["17"], "element-citation": ["\n"], "surname": ["del Hoyo", "Elliott", "Christie"], "given-names": ["J", "A", "DA"], "year": ["2006"], "source": ["Handbook of the Birds of the World"], "publisher-name": ["Lynx Edicions"]}, {"label": ["18"], "element-citation": ["\n"], "surname": ["Bensch", "Hasselquist"], "given-names": ["S", "D"], "year": ["1999"], "article-title": ["Phylogeographic population structure of great reed warblers: an analysis of mtDNA control region sequences."], "source": ["Biol J Linn Soc"], "volume": ["66"], "fpage": ["171"], "lpage": ["185"]}, {"label": ["22"], "element-citation": ["\n"], "surname": ["Baker", "Marshall", "Mindell"], "given-names": ["AJ", "HD", "DP"], "year": ["1997"], "article-title": ["Mitochondrial control region sequences as tools for understanding evolution."], "source": ["Avian molecular evolution and systematics"], "publisher-loc": ["San Diego, CA"], "publisher-name": ["Academic Press"], "fpage": ["51"], "lpage": ["82"]}, {"label": ["24"], "element-citation": ["\n"], "surname": ["Perez-Tris", "Bensch", "Carbonell", "Helbig", "Telleria"], "given-names": ["J", "S", "R", "AJ", "JL"], "year": ["2004"], "article-title": ["Historical diversification of migration patterns in a passerine bird."], "source": ["Evolution Int J Org Evolution"], "volume": ["58"], "fpage": ["1819"], "lpage": ["1832"]}, {"label": ["25"], "element-citation": ["\n"], "surname": ["Hewitt"], "given-names": ["GM"], "year": ["1996"], "article-title": ["Some genetic consequences of ice ages, and their role in divergence and speciation."], "source": ["Biol J Linn Soc"], "volume": ["58"], "fpage": ["247"], "lpage": ["276"]}, {"label": ["26"], "element-citation": ["\n"], "surname": ["De Roo", "Deheegher"], "given-names": ["A", "J"], "year": ["1969"], "article-title": ["Ecology of the great reed warbler "], "italic": ["Acrocephalus arundinaceus"], "source": ["Gerfaut"], "volume": ["59"], "fpage": ["260"], "lpage": ["275"]}, {"label": ["28"], "element-citation": ["\n"], "surname": ["Hansson", "Bensch", "Hasselquist"], "given-names": ["B", "S", "D"], "year": ["2002"], "article-title": ["Predictors of natal dispersal in great reed warblers: results from small and large census areas."], "source": ["J Avian Biol"], "volume": ["33"], "fpage": ["311"], "lpage": ["314"]}, {"label": ["30"], "element-citation": ["\n"], "surname": ["Hagemeijer", "Blair"], "given-names": ["WJM", "MJ"], "year": ["1997"], "source": ["The EBCC atlas of European breeding birds"], "publisher-loc": ["London"], "publisher-name": ["Poyser"]}, {"label": ["36"], "element-citation": ["\n"], "surname": ["Ibrahim", "Nichols", "Hewitt"], "given-names": ["KM", "RA", "GM"], "year": ["1996"], "article-title": ["Spatial patterns of genetic variation generated by different forms of dispersal during range expansion."], "source": ["Heredity"], "volume": ["77"], "fpage": ["282"], "lpage": ["291"]}, {"label": ["37"], "element-citation": ["\n"], "surname": ["Bensch", "Hasselquist"], "given-names": ["S", "D"], "year": ["1991"], "article-title": ["Territory infidelity in the polygynous great reed warbler "], "italic": ["Acrocephalus arundinaceus"], "source": ["J Anim Ecol"], "volume": ["60"], "fpage": ["857"], "lpage": ["871"]}, {"label": ["38"], "element-citation": ["\n"], "surname": ["Hansson", "Jack", "Christians", "Pemberton", "\u00c5kesson", "Westerdahl", "Bensch", "Hasselquist"], "given-names": ["B", "L", "JK", "JM", "M", "H", "S", "D"], "year": ["2007"], "article-title": ["No evidence for inbreeding avoidance in a great reed warbler population."], "source": ["Behav Ecol"], "volume": ["18"], "fpage": ["157"], "lpage": ["164"]}, {"label": ["39"], "element-citation": ["\n"], "surname": ["Hansson", "Bensch", "Hasselquist", "\u00c5kesson"], "given-names": ["B", "S", "D", "M"], "year": ["2001"], "article-title": ["Microsatellite diversity predicts recruitment of sibling great reed warblers."], "source": ["Proc R Soc Lond B"], "volume": ["268"], "fpage": ["1287"], "lpage": ["1291"]}, {"label": ["43"], "element-citation": ["\n"], "surname": ["Iles", "Bishop"], "given-names": ["MM", "DT"], "year": ["1998"], "article-title": ["The effect of population structure and mutation rate on linkage disequilibrium."], "source": ["Am J Hum Genet"], "volume": ["63"], "fpage": ["A42"]}, {"label": ["47"], "element-citation": ["\n"], "surname": ["Ottvall", "Bensch", "Walinder", "Lifjeld"], "given-names": ["R", "S", "G", "JT"], "year": ["2002"], "article-title": ["No evidence of genetic differentiation among Lesser and Common Redpolls."], "source": ["Avian Science"], "volume": ["4"], "fpage": ["237"], "lpage": ["244"]}, {"label": ["50"], "element-citation": ["\n"], "surname": ["Nei"], "given-names": ["M"], "year": ["1987"], "source": ["Molecular Evolutionary Genetics"], "publisher-loc": ["New York"], "publisher-name": ["Columbia Univ. Press"]}, {"label": ["54"], "element-citation": ["\n"], "surname": ["Strimmer", "von Haeseler"], "given-names": ["K", "A"], "year": ["1997"], "source": ["PUZZLE 4.0. Maximum likelihood analysis for nucleotide, amino acid, and two-state data"], "publisher-loc": ["Munchen"], "publisher-name": ["Zoologisches Institut, Universit\u00e4t Munchen"]}, {"label": ["56"], "element-citation": ["\n"], "surname": ["Nylander"], "given-names": ["JAA"], "year": ["2004"], "article-title": ["MrModeltest2."], "ext-link": ["http://www.ebc.uu.se/systzoo/staff/nylander.html"]}]
PLoS One. 2008 Jul 30; 3(7):e2794