- Research article
- Open Access
Multilocus dataset reveals demographic histories of two peat mosses in Europe
© Szövényi et al; licensee BioMed Central Ltd. 2007
- Received: 22 December 2006
- Accepted: 22 August 2007
- Published: 22 August 2007
Revealing the past and present demographic history of populations is of high importance to evaluate the conservation status of species. Demographic data can be obtained by direct monitoring or by analysing data of historical and recent collections. Although these methods provide the most detailed information they are very time consuming. Another alternative way is to make use of the information accumulated in the species' DNA over its history. Recent development of the coalescent theory makes it possible to reconstruct the demographic history of species using nucleotide polymorphism data. To separate the effect of natural selection and demography, multilocus analysis is needed because these two forces can produce similar patterns of polymorphisms. In this study we investigated the amount and pattern of sequence variability of a Europe wide sample set of two peat moss species (Sphagnum fimbriatum and S. squarrosum) with similar distributions and mating systems but presumably contrasting historical demographies using 3 regions of the nuclear genome (appr. 3000 bps). We aimed to draw inferences concerning demographic, and phylogeographic histories of the species.
All three nuclear regions supported the presence of an Atlantic and Non-Atlantic clade of S. fimbriatum suggesting glacial survival of the species along the Atlantic coast of Europe. Contrarily, S. squarrosum haplotypes showed three clades but no geographic structure at all. Maximum likelihood, mismatch and Bayesian analyses supported a severe historical bottleneck and a relatively recent demographic expansion of the Non-Atlantic clade of S. fimbriatum, whereas size of S. squarrosum populations has probably decreased in the past. Species wide molecular diversity of the two species was nearly the same with an excess of replacement mutations in S. fimbriatum. Similar levels of molecular diversity, contrasting phylogeographic patterns and excess of replacement mutations in S. fimbriatum compared to S. squarrosum mirror unexpected differences in the demography and population history of the species.
This study represents the first detailed European wide phylodemographic investigation on bryophytes and shows how pattern of nucleotide polymorphism can reveal unexpected differences in the population history of haploid plants with seemingly similar characteristics.
- Last Glacial Maximum
- Historical Demography
- Mismatch Distribution
- Demographic Expansion
- Ancient Population
Geographic distribution of species is continuously shaped by several extrinsic factors causing repeated historical range and demographic expansions/contractions [1, 2]. Using historical observations, short-term range expansions and contractions, together with associated demographic changes, have been documented in several, mostly invasive species in Europe and also world-wide [3–6]. From a conservation biological point of view, it is of high importance to reconstruct the historical demography of populations, which can help to work out conservation strategies. In addition, correlations between demographic changes and climatic oscillations may be useful to predict species responses to the accelerated climate change . If a series of temporarily spaced data is available, demographic changes over several decades are relatively easy to reconstruct. However, comparison of historical and recent observations is challenging, among others because of the difficulties emerging when comparing data gathered by different sampling strategies .
Methods using the current level and distribution of genetic variability to estimate the historical demography of populations provide an additional tool to test demographic hypotheses suggested by simple historical observations [9–12]. Historical demography considerably influences both qualitative and quantitative properties of sequence-level polymorphisms at neutrally evolving sites [13, 14]. Consequently, based on the current pattern of sequence-level polymorphisms, it is possible to distinguish among populations, which currently underwent a demographic expansion and those, which are in a declining phase. Populations after a demographic expansion are described by rapid coalescent events following long branches on the coalescent tree . In the mismatch distribution, long branches lead to an excess of low frequency polymorphisms [15, 16].
Unfortunately, selective forces might lead to the same pattern of molecular variation as would be expected in an expanding population. Investigating noncoding, neutrally evolving sequence portions does not overcome the problem of selection because they might be influenced by hitchhiking effects of linked loci . Although selective forces and demography might produce the same pattern of polymorphism, selection affects the genome locally, whereas demography extends its effect on the whole genome. Consequently, investigating multiple, unlinked loci can help to separate the influence of selection and demography on the pattern of nucleotide polymorphism [18, 19].
In the present study, we aimed to reconstruct the demographic histories of two peat moss species based on nucleotide polymorphism data. Peat mosses are an ancient plant group with a worldwide distribution  and represent a model system for population genetic research in bryophytes [21–24]. In addition, phylogenetically closely related species pairs with similar distribution patterns and life history characteristics provide an outstanding opportunity to test the influence of species-specific demography on nucleotide polymorphism.
To concentrate on the effect of demography and minimize the influence of other possible factors on the pattern and amount of sequence-level variability, two phylogenetically closely related species with similar distributions, ecological requirements and mating systems have been selected. Both species are monoecious (produce antheridia and archegonia on the same gametophore) and distributed from the Iberian Peninsula to Svalbard across Europe. They frequently occur together at the edge of mires under Alnus or Salix stands and can tolerate higher nutrient concentrations [20, 25]. Although they share similar characteristics, recent investigations show that S. fimbriatum is more successful in colonizing open soil surfaces and produces more rapidly sporophytes after establishment than S. squarrosum . In addition, analyses of past and current distributions of the species show that S. fimbriatum is presumably experiencing a rapid expansion at least in some parts of Europe, whereas S. squarrosum shows no evidence of a recent population demographic change .
Previous investigations on a limited number of accessions using chloroplast sequences, suggested contrasting phylogeographic patterns in the monoecious S. fimbriatum and S. squarrosum in Europe . Although chloroplast markers allowed drawing of phylogeographic hypotheses, very low resolution of the dataset hampered any further statistical inference. Moreover, the chloroplast genome consists of completely linked loci and thus may be strongly influenced by selective forces. Consequently, it is impossible to separate the effect of demography and natural selection on nucleotide polymorphisms using only chloroplast sequence data.
Therefore, this study aims to investigate the historical demography of the two species in detail using more detailed sampling, and sequence data of three unlinked nuclear regions covering appr. 3000 bps of the genome. We intended to differentiate between selective and demographic forces and were particularly interested in the following questions: 1. Does phylogeographic structure depend on the loci investigated? 2. Does multilocus analysis support assumed population expansion of S. fimbriatum? 3. What kind of demographic histories could lead to the current amount and pattern of molecular variability?
Geographic distribution of haplotypes
Molecular polymorphism, tests of neutrality and recombination
Nucleotide polymorphisms of the regions investigated. Number of sequences, aligned length, nucleotide diversity estimates and neutrality tests of the nuclear regions investigated. Tajima's D was calculated taking nucleotide substitutions or nucleotide substitutions and gaps into account.
Aligned lenght in bps
Tajima's D (with gaps/without gaps)
Fay and Wu's H
Total data set
Total data set
The RAPDa region has been excluded from the multilocus HKA test because it showed no point mutations in S. squarrosum and only the microsatellite repeat was variable. No significant deviation from the neutral model was detected using the multilocus HKA test (p = 0.48). Despite of the excess of non-synonymous substitutions in S. fimbriatum compared to S. squarrosum (Table 1), the McDonald-Kreitman test (GapC gene) was not significant either (p = 0.50). In S. squarrosum, Tajima's D and Fay and Wu's H statistics were never significantly different from the neutral expectations.
Using all European accessions of S. fimbriatum, Tajima's D was significant for all three nuclear regions investigated when including gaps (Table 1). However, it turned out to be non-significant when excluding gaps, except for the GapC gene. Fay and Wu's H was never significant. Analysing accessions of the Non-Atlantic clade of S. fimbriatum, Tajima's D was significantly negative except the RAPDa region. The latter showed a significant negative value as well when gaps were included. The Atlantic clade showed no significant deviations from neutrality.
No signs of recombination were found in RAPDa and ITS data sets of both species. GapC sequences of S. fimbriatum showed one incompatible pair and an Rm value of 1. The two accessions causing the incompatibility were removed from the data set in further analyses. No signs of recombination were found in GapC sequences of S. squarrosum.
Test of population growth using Fu's Fs statistic. Calculations were done including only nucleotide substitutions or nucleotide substitutions and gaps in the analysis.
Fu's Fs (with gaps/without gaps)
Total data set
Test of population growth using mismatch distribution analysis. Analysis was conducted using all European accessions of both species. Note that estimations are given only for the GapC gene.
Pure demographic expansion
Test of population growth using maximum likelihood estimation. Analysis was conducted using all European accessions of both species. Note that estimations are given only for the GapC gene.
Maximum likelihood estimates
Fitting the data to the IM model
In contrast to S. fimbriatum, current thetas in S. squarrosum almost gave the same point estimates. Point estimate of the ancestral neutral mutation parameter was an order of magnitude higher than that of the current populations, but the posterior distribution was relatively flat. It converged to zero at higher theta values but did not reach zero even after extending the prior of theta considerably. Scaled migration parameters showed non zero values but were flat. Posterior probability distribution of the splitting parameter was flat and uninformative (not shown).
Results of the present analysis confirm the divergence of an Atlantic and a Non-Atlantic lineage in S. fimbriatum and almost no geographic population structure in S. squarrosum . Moreover, our current results evidently support a genome-wide historical event, which influenced all regions in the same way regardless of genomic position.
Since there is no reliable calibration point in the history of the genus, it is not possible to date the split in S. fimbriatum. In absence of dating, at least two hypotheses are plausible. The first one assumes that Quaternary cold periods were long enough to allow accumulation of the observed number of mutations among the Atlantic and Non-Atlantic groups of populations. In addition, it presumes that current geographic position of populations roughly correspond to that of the historical ones. The Atlantic group might have survived along the Atlantic coast of Spain, France and South England, but location of refugia of the Non-Atlantic group remains ambiguous. Results suggest that Eastern Iberian populations represent recent colonization events. Using reconstructions of Gajewski et al.  based on peat moss spores we hypothesize survival of the Non-Atlantic lineage along the border of Poland and Slovakia where extensive Sphagnum dominated habitats occurred in about 19 ky before present.
The second hypothesis assumes intercontinental migration between Europe and America. North American and European populations of several Sphagnum species appear to be closely related [30–33]. Among others, long-range dispersal as a plausible explanation for the observed pattern has been suggested [26, 34, 35]. Assuming effective intercontinental migration, colonization from a divergent American population might account for the difference between the two lineages. To test this hypothesis we analysed additional samples from different parts of North America. Most of the North American accessions either shared haplotypes with the Non-Atlantic group or differed by only one substitution and represented derived tip haplotypes. Therefore recent colonization from America does not account for the divergence of the Atlantic and Non-Atlantic clades.
In addition, several issues indicate that the Quaternary period has played an important role in the observed genetic structure of the species. First, the patterns detected in this study were very similar to those found in other bryophyte species [36–38] and seed plants [39, 40]. The second reason comes from fossil data and climatic reconstructions. Occurrence of relic Sphagnum species , fossil remains of peat moss spores or leaves [41, 42] and climatic reconstructions also support survival of populations along the Atlantic coast of Europe [29, 43]. Consequently, our results and the above-mentioned observations support the separation of the Atlantic and Non-Atlantic lineages of S. fimbriatum by Quaternary glacial cycles and reject recent colonization from another continent as the cause of the split.
Selection influences each gene differently, whereas demographic processes affect the individual as a whole and thus influence all parts of the genome uniformly regardless of genomic position . Hence, in an expanding population an excess of low frequency variants is expected for all loci. This leads to negative Tajima's D and Fu's Fs statistics [44, 45]. These statistics all analyse the total or relative number of singleton mutations compared to the neutral coalescent. Their power depends on the total number of segregating sites in the data set and the time elapsed since the demographic event [45, 46]. Therefore it is not surprising that some of the statistics were not significant when excluding gaps as segregating sites from the analyses. This shows that substitution rates are relatively low compared to the time elapsed since the demographic event. Increasing the mutation rate with including gaps as mutations in the calculations improves power of the tests and extends the time frame in which those are applicable. Using these assumptions, all test statistics for all regions of S. fimbriatum were significantly different from the neutral expectations, whereas none of them provided significant deviations in S. squarrosum. This confirms a general excess of low frequency variants, in S. fimbriatum especially in the Non-Atlantic clade. Since the regions used are unlinked, we favour a population demographic expansion against a genome wide selection hypothesis in S. fimbriatum. Sequencing errors are negligible because both species were analysed using the same protocols and polymorphisms have been rechecked on the original chromatograms.
Mismatch distribution analysis revealed further details about the demographic history of the two species. Mismatch distribution of spatial expansion only slightly differs from that of pure demographic expansion, the latter showing a higher frequency of identical sequences . In case of S. fimbriatum, data support the spatial expansion model, which is in concordance with the hypotheses about rapid spatial expansion of the species after the Quaternary glaciations. Data from S. squarrosum fit both scenarios, however, considerably better match the latter. This suggests that S. squarrosum probably went through spatial expansion as well. Tau value of S. squarrosum was approximately an order of magnitude higher than that of S. fimbriatum, although they were not significantly different. This suggests that expansion of S. fimbriatum is more recent than that of S. squarrosum. Assuming that S. fimbriatum expanded immediately after the Last Glacial Maximum (LGM), the estimate for S. squarrosum suggests expansion related to earlier geological or climatic events than the LGM. However, results of mismatch analyses need to be interpreted with caution. These analyses do not use genealogical information included in the data and testing the fit between observed and expected distributions and estimating approximate confidence intervals represent a complex problem as well. Estimation and confidence intervals of tau and θ are skewed and biased, especially in data sets with low resolution . Maximum likelihood (LAMARC) and Bayesian analyses (IM) also supported rapid population expansion in S. fimbriatum. Results of both analyses showed no significant evidence for population growth in S. squarrosum. The IM analysis even supported a population decline since the separation of the two lineages of S. squarrosum.
Information on the nucleotide polymorphism of three independent regions gives the opportunity to reconstruct the historical demography of the two species in details. Populations of S. fimbriatum and S. squarrosum have been separated by the Quaternary glaciations, which led to the differentiation of two and three genetically well diverged European lineages, respectively. Fixed replacement mutations among the Atlantic and Non-Atlantic groups of S. fimbriatum indicate that both lineages went through a severe bottleneck, which led to the fixation of slightly deleterious mutations. It is worth noting that populations of the Atlantic group retained relatively divergent haplotypes (Figure 1) despite of their scattered occurrence and small extent (occurrences are usually not larger than 5 × 5 m). Hence, these populations very likely represent the leftovers of a formerly more widely distributed genetic lineage. Fitting the IM model to the data set suggests that a very high proportion of the ancient population of S. fimbriatum founded the Atlantic clade and the size of the ancient population was very small. These results show that the Atlantic clade is very likely a direct descendant of the ancestral population, which retained ancestral haplotypes and went through a severe bottleneck. IM and maximum likelihood estimations show that the Quaternary glaciations shaped population size of the two species differently. Ancient population size of S. squarrosum is considerably larger than that of S. fimbriatum indicating a more severe bottleneck in the latter. Larger effective population size of the ancient population of S. squarrosum also explains the lack of replacement mutations in coding regions as well as the lack of geographic structure of haplotypes.
Pattern and amount of molecular variability
Overall diversity values for the ITS region were usually an order of magnitude higher in S. fimbriatum than in S. squarrosum, however all values were very low. Moreover, in S. fimbriatum, a considerable portion of the diversity was due to the genetic differentiation between Atlantic and Non-Atlantic clades. In contrast, GapC sequences of the species showed nearly the same nucleotide diversity values. This is an unexpected result because the pioneer characteristic of S. fimbriatum (Sundberg et al. 2006), even when combined with a parallel reduction of the historical population size of S. squarrosum, implies lower level of molecular polymorphism in the first species. Peat moss populations are likely to function as a metapopulation [49, 50]. Theory shows that species-wide genetic diversity is rapidly lost in metapopulations of pioneer species if among population migration rates are lower than extinction/recolonization rates [51, 52]. Although pioneer characteristic reduces the level of molecular diversity, in range expansion, due to the phenomenon of surfing mutations, molecular diversity is increased by a higher proportion of retained neutral mutations . Consequently, it seems likely that range expansion contributed significantly to the molecular variability of S. fimbriatum.
Although the magnitude of molecular polymorphism was similar in both species, its pattern differed considerably showing an excess of polymorphism in coding regions of S. fimbriatum compared to S. squarrosum. Recent increase of effective population size results in a significant excess of slightly deleterious replacement mutations . Weakened selection due to reduced effective population size can lead to the accumulation of slightly deleterious mutations as well and thus to elevated nucleotide diversity in selfing plants . The McDonald-Kreitman test showed no significant deviations from neutral expectations for the GapC gene, however, exons of this gene contained replacement mutations in S. fimbriatum whereas in S. squarrosum all mutations were restricted to noncoding regions. Replacement mutations in S. fimbriatum (Lys-Arg, Ile-Val, Lys-Gln) caused no change in hydrophoby and almost no change in pI values and thus are likely to be only slightly deleterious. Consequently, elevated number of replacement mutations in S. fimbriatum also supports a rapid population demographic expansion in Europe.
Both species are easily dispersed by spores to longer distances [26, 27], however, nuclear markers resolved distinct Atlantic and Non-Atlantic clades only in S. fimbriatum. All markers indicate either restricted migration or establishment potential measured on the time scale of the markers' mutation rates. S. fimbriatum produces considerable amounts of spores, which can germinate on a wider range of substrates (especially on substrates with low phosphate availability) compared to other Sphagnum species [26, 56]. Based on these characteristics of the species, the Non-Atlantic clade of S. fimbriatum might have rapidly recolonized available soil surfaces after the last glaciation. In absence of further space, the Atlantic clade might have been unable to follow this lineage and remained restricted to the Atlantic coast of Europe. Similar patterns of leading edge colonization have been observed in several tree species in Europe [57, 58].
Overall nucleotide diversity of the three regions was low in both species compared to estimations in seed plants [reviewed in ], which might be a general trend in bryophytes. Due to the haploid gametophores, no sheltering of recessive alleles exists in bryophytes and the haploid genome of each individual is directly exposed to selection [30, 59]. This, coupled with the presence of deleterious mutations, considerably reduces the amount of genetic diversity. Mutation rate also influences genetic variability. Selection against high mutation rates should exist in bryophytes; otherwise the presence of slightly deleterious mutations will lead to a severe mutational load . Our results show that slightly deleterious mutations are present but are removed very efficiently. The theory of selection directly acting on the gametophores and low mutation rates are supported by the data set presented here. Assuming the same mutation rates in both species, selection is so efficient, that no replacement mutations were found in the GapC gene of S. squarrosum at all. Interestingly, in a previous study McDaniel and Shaw  investigated the level of nucleotide polymorphism in the moss Ceratodon purpureus which turned out to be an order of magnitude higher than in our species. As this investigation relied on a worldwide sample and used multi-copy genes, high level of polymorphism may reflect the complex interplay of various forces, such as evolution of gene families, selective sweeps and phylogeographic history. Hence, the two studies are hardly comparable and further investigations on the extent of population level sequence polymorphism are needed to test theoretical predictions. In contrast to our results, high levels of isozyme diversity were found within populations of several bryophytes [30, 60]. This has partly been interpreted as a result of local adaptation . Although local adaptation might be important in several species, our results do not support this hypothesis, since only one of the three replacement mutations in the GapC gene resulted in charge change. Hence these mutations probably do not influence enzyme functions significantly.
Analysing nucleotide polymorphisms of multiple loci combined with the application of the coalescent theory provides an effective way to draw inferences about historical demography of species. Furthermore, investigating polymorphisms of several independent loci makes it possible to separate the effects of selective and demographic forces. Our study evidently shows that species with similar ecology, mating systems and current distributions may have very different demographic and phylogeographic histories. Such cryptic differences are probably frequent among bryophytes and must be taken into account when assessing the vulnerability of species and planning conservation strategies. Due to their haploid gametophores, bryophytes are especially appropriate for studies on sequence level polymorphisms. Future research using nucleotide polymorphism data will certainly discover undetected demographic, phylogeographic and evolutionary processes in this group of haploid plants.
Sampling and DNA extraction
S. fimbriatum and S. squarrosum plants were collected by the authors, provided by colleagues or obtained from herbaria [see Additional file 1]. Plants were dried and stored in silica gel at 4°C. DNA was extracted using the DNAeasy plant mini kit (Qiagen, Switzerland) following the manufacturers protocol. Three regions of the nuclear genome were used for sequencing. The region spanning the 5th–9'th exons of the GapC (glyceraldehyde-3-phosphate-dehydrogenase C) gene (appr. 1500 bps) was amplified using protocols and primers described in Szövényi et al. . A genomic region with unknown function [RAPDa sensu ] was also used. New primers were designed to improve amplification and sequencing of the fragment. As a third locus the ITS1-5.8S-ITS2 region (hereafter ITS) were amplified and sequenced. Yield of PCR reactions were improved by applying a semi-nested PCR, where 1 μl of the 10 × diluted first round product was used as template in the second reaction [see Additional file 2]. The PCR reaction contained the same amount of components as described in Szövényi et al.  except that final DMSO concentration was set to 1% v/v in the PCR mastermix to amplify the ITS region. PCR reactions were run on Biometra (T1 or Tgradient, Whatman Biometra) and Techne ptc 412 (Barloworld Scientific Ltd.) thermocyclers. The PCR program was the following: 94°C denaturation for 4 min, then 94°C for 1 min, different annealing temperatures for 1 min and 72°C extension for 1–1.5 min [see Additional file 2] with 35 cycles with a final extension step at 72°C for 7 min. The products were checked on 0.8% agarose gels and cleaned using the GFX PCR and gel band purification kit (Amersham Biosciences, Switzerland). Approximately 10 ng product was sequenced in 10 μl cycle sequencing reaction with the Big dye v3.1 in an ABI prism 3100 (Applied Biosystems) genetic analyzer using either the original and/or different internal PCR primers [see Additional file 2]. Sequences were contiged and corrected if necessary with the Sequencher 4.5 software (Gene Code Corporation).
Sequence alignment, molecular diversity and recombination event estimates
Sequences were aligned using ClustalW , checked by eye and adjusted if needed. All polymorphic sites were rechecked on the original electropherograms and corrected to avoid false base callings. Haplotype networks were reconstructed using the TCS software . To describe the infra specific molecular diversity of the regions Watterson's (θw)  and Tajima's (θπ)  moment estimators of the neutral mutation parameter were calculated. In case of the GapC gene, levels of polymorphism were calculated for synonymous, nonsynonymous and silent sites as well.
All regions were tested for recombination using the four gamete test and the minimum number of recombination events during the history of the sample was calculated [Rm, ]. For several following analyses (cf. below) sequences with incompatible nucleotide pairs were cut into recombination free blocks. All calculations were done using dnaSP v.4.10.4 .
Testing deviation from neutrality
To test whether the sequenced regions evolved under neutral expectations, four different methods were applied. The multilocus HKA test was conducted using the HKA software . 10 000 coalescent simulations were run to generate the null distribution of χ2 values. Neutral evolution of the GapC gene was also tested with the McDonald-Kreitman test . To test the neutraliy of each locus separately, Tajima's D  and Fay and Wu's H statistics  were calculated for each of them. For each locus the latter two test statistics were calculated considering nucleotide substitutions alone and recoding gaps as well. Gaps were coded as single and independent events. Due to the low variability and non-overlapping occurrence of gaps in the data set, alignment, assessment of homology and subsequent coding of gaps was straightforward.
Analysis of population growth
As a powerful descriptive statistic of population growth, Fu's Fs , was calculated using frequency distribution of alleles for each region separately. Fu's Fs statistic compares the number of alleles expected under the neutral coalescent in a constant sized population based on the neutral mutation parameter and the size of the sample. In an expanding population large negative values are expected due to the excess of rare mutations. Populations after a demographic expansion show a star-like phylogeny and a unimodal mismatch distribution , where the peak corresponds to the time of the expansion. The method of Schneider and Excoffier  implemented in ARLEQUIN 3.0  was used to fit observed and expected mismatch distributions and to estimate parameters of a pure demographic and a spatial expansion model. The two models are expected to produce different mismatch distributions when the number of migrants between populations is not large. Otherwise the two hypotheses can not be distinguished. In this study both models were considered to complement genealogical methods that fail to detect signs of a spatial expansion.
As an alternative, maximum likelihood estimation of the exponential growth rate and historical θ was conducted. The software LAMARC 2.0  was used to get estimates of growth rates for all European accessions of S. squarrosum and S. fimbriatum and for the Non-Atlantic clade of S. fimbriatum. 30 short chains (20 000 steps each) and 2 long chains (200 000 steps) with a sampling interval of 20 steps and a burn in period of 1000 and 10 000 respectively were applied.
The IM software [76, 77] was used to fit the data to the Isolation with Migration model (IM). This model assumes two populations connected by migration, which are derived by the splitting of an ancient population. The IM model estimates several parameters simultaneously and therefore describes processes in natural populations in their complexity. Populations were defined based on our a priory knowledge of chloroplast lineages using the parsimony networks of each species . Runs were conducted for the combined dataset including all three nuclear regions sequenced. Initial runs with wide priors were used to delimit the plausible range of priors. Since indels of all three regions were informative, analyses have been conducted both with and without indels coded. Including gaps as additional characters did not change the multilocus parameter estimates but narrowed their confidence intervals, however, only results without coding gaps are shown. The infinite site mutation model of sequence evolution was applied and loci were cut into non-recombining blocks. In case of the RAPDa sequences of S. squarrosum, only the microsatellite region was used. Each run was introduced by a 100 000 steps long burn-in period. Multilocus analyses were run using 5–7 parallel chains under a linear heating scheme with a heating value of 0.05 – 0.005. Convergence of parameters and mixing of chains were followed by visual inspection of parameter trend lines and checking of ESS values. Analyses were run until the lowest ESS value reached minimum 200. In all analyses population demographic changes were allowed.
We thank Liudmila Abramova, Baiba Bambe, Sam Bosanquet, Jürgen Brandt, Philippe de Zuttere, Christian Dolnik, Andrea Fodor, Anna Ganeva, Irina Goldberg, Pierre Goubet, Alessia Guggisberg, Patxi Heras, Heike Hofmann, Adam Hölzer, Michael Ignatov, Luca Miserere, Zoltán Illyés, Marta Infante, Nele Ingerpuu, Beatriz Itten, Monika Koperski, Robert Krisai, Ildikó László, Andrej Martinčič, Péter Ódor, Beáta Papp, Claes Persson, Ron Porley, Markus Preussing, José Pujos, Sigurd Såstad, Ana Seneca, Blanka and Jonathan Shaw, Rudolf Soltés, Sture Westerberg, Sebastian Sundberg, Erzsébet Szurdoki, Zoltán Tóth, Edwin Urmi, Alain Vanderpoorten, Kai Vellak, Risto Virtanen and Harald Zechmeister for collecting material and/or for field assistance. Critical comments of two anonymous reviewers on an earlier version of the manuscript are greatly acknowledged. This research was financially supported by a Young Academics Grant of the University of Zurich (Forschungskredit der Universität Zürich), by a Travel Grant of the Swiss National Academy of Science and by a grant of the Hungarian Scientific Research Found (OTKA, NI 68218).
- Hewitt GM: Post-glacial re-colonization of European biota. Biol J Linn Soc. 1999, 68: 87-112. 10.1006/bijl.1999.0332.View ArticleGoogle Scholar
- Hewitt GM: The genetic legacy of the Quaternary ice ages. Nature. 2000, 405: 907-913. 10.1038/35016000.View ArticlePubMedGoogle Scholar
- Fraser SM, Lawton JH: Host range expansion by British moths onto introduced conifers. Ecol Entomol. 1994, 19: 127-137.View ArticleGoogle Scholar
- Losos JB: Dynamics of range expansion by three introduced species of Anolis lizards on Bermuda. J Herpetol. 1996, 30: 204-210. 10.2307/1565511.View ArticleGoogle Scholar
- Sakai AK, Allendorf FW, Holt JS, Lodge DM, Molofsky J, With KA, Baughman S, Cabin RJ, Cohen JE, Ellstrand NC, McCauley DE, O'Neil P, Parker IM, Thompson JN, Weller SG: The population biology of invasive species. Annu Rev Ecol Syst. 2001, 32: 305-332. 10.1146/annurev.ecolsys.32.081501.114037.View ArticleGoogle Scholar
- Radosevich SR, Stubbs MM, Ghersa CM: Plant invasions--process and patterns. Weed Sci. 2003, 51: 254-259. 10.1614/0043-1745(2003)051[0254:PIPAP]2.0.CO;2.View ArticleGoogle Scholar
- Weber E: Current and potential ranges of three exotic goldenrods (Solidago) in Europe. Conserv Biol. 2001, 15: 122-128. 10.1046/j.1523-1739.2001.99424.x.View ArticleGoogle Scholar
- Hedenäs L, Bisang I, Tehler A, Hamnede M, Jaederfelt K, Odelvik G: A herbarium-based method for estimates of temporal frequency changes: mosses in Sweden. Biol Conserv. 2002, 105: 321-331. 10.1016/S0006-3207(01)00212-9.View ArticleGoogle Scholar
- Beaumont MA: Detecting population expansion and decline using microsatellites. Genetics. 1999, 153: 2013-2029.PubMed CentralPubMedGoogle Scholar
- Storz JF, Beaumont MA, Alberts SC: Genetic evidence for long-term population decline in a savannah-dwelling primate: inferences from a hierarchical Bayesian model. Mol. Biol. Evol. 2002, 19: 1981-1990.View ArticlePubMedGoogle Scholar
- Kuhner MK, Yamato J, Beerli P: LAMARC v 2.0. 2005, University of Washington, [http://evolution.gs.washington.edu/lamarc.html]Google Scholar
- Goossens B, Chikhi L, Ancrenaz M, Lackman-Ancrenaz I, Anaud P, Bruford MW: Genetic signature of anthropogenic population collapse in orangutans. PLoS Biology. 2006, 4: e25-10.1371/journal.pbio.0040025. doi:10.1371/journal.pbio.0040025PubMed CentralView ArticlePubMedGoogle Scholar
- Rosenberg NA, Nordborg M: Genealogical trees, coalescent theory and the analysis of genetic polymorphisms. Nat Rev Genet. 2002, 3: 380-390. 10.1038/nrg795.View ArticlePubMedGoogle Scholar
- Hein J, Schierup MH, Wiuf C: Gene Genealogies, Variation and Evolution. A Primer in Coalescent Theory. 2004, Oxford: Oxford University PressGoogle Scholar
- Slatkin M, Hudson RR: Pairwise comparisons of mitochondrial DNA sequences in stable and exponentially growing populations. Genetics. 1991, 129: 555-562.PubMed CentralPubMedGoogle Scholar
- Rodgers AR, Harpending H: Population growth makes waves in the distribution of pairwise genetic differences. Mol. Biol. Evol. 1992, 9 (3): 552-569.Google Scholar
- Maynard Smith J, Haigh J: The hitchhiking effect of a favourable gene. Genet Res Camb. 1974, 219: 1114-1116.Google Scholar
- Wright SI, Gaut BS: Molecular population genetics and the search for adaptive evolution in plants. Mol. Biol. Evol. 2004, 22: 506-519. 10.1093/molbev/msi035.View ArticlePubMedGoogle Scholar
- Nielsen R: Molecular signatures of natural selection. Annu Rev Genet. 2005, 39: 197-218. 10.1146/annurev.genet.39.073003.112420.View ArticlePubMedGoogle Scholar
- Daniels RE, Eddy A: Handbook of European Sphagna. 1985, Huntingdon: Institute of Terrestrial EcologyGoogle Scholar
- Cronberg N: Clonal structure and fertility in a sympatric population of the peat mosses Sphagnum rubellum and Sphagnum capillifolium. Can J Bot. 1996, 74: 1375-1385.View ArticleGoogle Scholar
- Cronberg N, Natcheva R: Hybridization between the peat mosses, Sphagnum capillifolium and S. quinquefarium (Sphagnaceae, Bryophyta) as inferred by morphological characters and isozyme markers. Plant Syst Evol. 2002, 234: 53-70. 10.1007/s00606-002-0196-y.View ArticleGoogle Scholar
- Shaw AJ, Melosik I, Cox CJ, Boles SB: Divergent and reticulate evolution in closely related species of Sphagnum section Subsecunda. Bryologist. 2005, 108: 363-376. 10.1639/0007-2745(2005)108[0363:DAREIC]2.0.CO;2.View ArticleGoogle Scholar
- Flatberg KI, Thingsgaard K, Såstad SM: Interploidal gene flow and introgression in bryophytes. Sphagnum girgensohnii × S. russowii, a case of spontaneous neotriploidy. J Bryol. 2006, 28: 27-37. 10.1179/174328206X90459.View ArticleGoogle Scholar
- Smith AJE: The Moss Flora of Britain and Ireland. 2004, Cambridge: Cambridge Univeristy Press, SecondView ArticleGoogle Scholar
- Sundberg S, Hansson J, Rydin H: Colonisation of Sphagnum on land uplift islands in the Baltic Sea: time, area, distance and life history. J Biogeogr. 2006, 33: 1479-1491. 10.1111/j.1365-2699.2006.01520.x.View ArticleGoogle Scholar
- Szurdoki E, Ódor P: Distribution and expansion of Sphagnum fimbriatum Wils. in Hungary. Lindbergia. 2005, 29: 136-142.Google Scholar
- Szövényi P, Hock Zs, Urmi E, Schneller JJ: Contrasting phylogeographic patterns in Sphagnum fimbriatum and S. squarrosum (Bryophyta, Sphagnopsida) in Europe. New Phytol. 2006, 172: 784-794. 10.1111/j.1469-8137.2006.01870.x.View ArticlePubMedGoogle Scholar
- Gajewski K, Viau A, Sawada M, Atkinson D, Wilson S: Sphagnum peatland distribution in North America and Eurasia during the past 21,000 years. Global Biogeochem Cyc. 2001, 15: 297-310. 10.1029/2000GB001286.View ArticleGoogle Scholar
- Wyatt R, Odrzykoski I, Stoneburner J: Isozyme evidence of reticulate evolution in mosses: Plagiomnium medium is anallopolyploid of P. ellipticum × P. insigne. Syst Bot. 1992, 17: 532-550. 10.2307/2419725.View ArticleGoogle Scholar
- Shaw AJ, Schneider RE: Genetic biogeography of the rare „copper moss“ Mielichhoferia elongata (Bryaceae). Am J Bot. 1995, 82: 8-17. 10.2307/2445780.View ArticleGoogle Scholar
- Såstad SM, Flatberg IK, Hanssen L: Origin, taxonomy and population structure of the allopolyploid peat moss Sphagnum majus. Plant Syst Evol. 2000, 225: 73-84. 10.1007/BF00985459.View ArticleGoogle Scholar
- Shaw AJ, Werner O, Ros RM: Intercontinental mediterranean disjunct mosses: morphological and molecular patterns. Am J Bot. 2003, 90: 540-550.View ArticlePubMedGoogle Scholar
- Van Zanten BO: Experimental studies on trans-oceanic long-range dispersal of moss spores in the southern hemisphere. J Hattori Bot Lab. 1978, 44: 455-482.Google Scholar
- Cronberg N: Colonization dynamics of the clonal moss Hylocomium splendens on islands in a Baltic land uplift area: reproduction, genet distribution and genetic variation. J Ecol. 2002, 90: 925-935. 10.1046/j.1365-2745.2002.00723.x.View ArticleGoogle Scholar
- Cronberg N: Population structure and interspecific differentiation of the peat moss sister species Sphagnum rubellum and S. capillifolium (Sphagnaceae) in northern Europe. Plant Syst and Evol. 1998, 209: 139-158. 10.1007/BF00985226.View ArticleGoogle Scholar
- Van der Velde M, Bijlsma R: Phylogeography of five Polytrichum species within Europe. Biol J Linn Soc. 2003, 78: 203-213. 10.1046/j.1095-8312.2003.00151.x.View ArticleGoogle Scholar
- Hedderson TA, Nowell TL: Phylogeography of Homalothecium sericeum (Hedw.) Br Eur; toward a reconstruction of glacial survival and postglacial migration. J Bryol. 2006, 28: 283-292. 10.1179/174328206X136287.View ArticleGoogle Scholar
- Rendell S, Ennos RA: Chloroplast diversity in Calluna vulgaris (heather) populations in Europe. Mol Ecol. 2002, 11: 69-78. 10.1046/j.0962-1083.2001.01413.x.View ArticlePubMedGoogle Scholar
- Hofmann MH, Glass AS, Tomiuk J, Schmuths H, Fritsch RM, Bachman K: Analysis of molecular data of Arabidopsis thaliana (L.) Heyn. (Brassicaceae) with geographical information systems (GIS). Mol Ecol. 2003, 12: 1007-1019. 10.1046/j.1365-294X.2003.01791.x.View ArticleGoogle Scholar
- Forment P, Jovet-Ast S: Sphagnum imbricatum Russ., subfossil en France dans une tourbe des environs de Calais. Rev bryol lichénol. 1950, 19: 125-126.Google Scholar
- Infante M, Heras P: Sphagnum imbricatum Hornsch. ex Russ., semifosil, en la turbera de Saldropo-Brazar (Ceanuri, Vizcaya). Est Inst Alavés Nat. 1987, 2: 201-203.Google Scholar
- VanAndel TH, Tzedakis PC: Paleolithic landscapes of Europe and environs, 150,000–25,000 years ago: an overview. Quaternary Sci Rev. 1996, 15: 481-500. 10.1016/0277-3791(96)00028-5.View ArticleGoogle Scholar
- Tajima F: Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989, 123: 585-595.PubMed CentralPubMedGoogle Scholar
- Ramos-Onsins SE, Rosas J: Statistical properties of new neutrality tests against population growth. Mol Biol Evol. 2002, 19: 2092-2100.View ArticlePubMedGoogle Scholar
- Sano A, Tachida H: Gene genealogy and properties of test statistics of neutrality under population growth. Genetics. 2005, 169: 1687-1697. 10.1534/genetics.104.032797.PubMed CentralView ArticlePubMedGoogle Scholar
- 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. 10.1046/j.1365-294X.2003.02004.x.View ArticlePubMedGoogle Scholar
- Schneider S, Excoffier L: Estimation of past demographic parameters from the distribution of pairwise differences when the mutation rates vary among sites: application to human mitochondrial DNA. Genetics. 1999, 152: 1079-1089.PubMed CentralPubMedGoogle Scholar
- Söderström L, Herben T: Dynamics of bryophyte metapopulations. Advances in Bryology. Edited by: Longton RE. 1997, Stuttgart: Science Publishers, 6: 205-241.Google Scholar
- Wakeley J, Alicar N: Gene genealogies in a metapopulation. Genetics. 2001, 159: 893-905.PubMed CentralPubMedGoogle Scholar
- Panell JR, Charlesworth B: Neutral genetic diversity in a metapopulation with recurrent local extinction and recolonization. Evolution. 2001, 53: 664-676. 10.2307/2640708.View ArticleGoogle Scholar
- Ingvarsson PK: A metapopulation perspective on genetic diversity and differentiation in partially self-fertilizing plants. Evolution. 2002, 56: 2368-2373.View ArticlePubMedGoogle Scholar
- Klopfstein S, Currat M, Excoffier L: The fate of mutations surfing on the wave of a range expansion. Mol Biol Evol. 2006, 23: 482-490. 10.1093/molbev/msj057.View ArticlePubMedGoogle Scholar
- Eyre-Walker A: Changing effective population size and the McDonald-Kreitman test. Genetics. 2002, 162: 2017-2024.PubMed CentralPubMedGoogle Scholar
- Charlesworth D, Wright SI: Breeding system and genome evolution. Curr Opin Genet Dev. 2001, 11: 685-690. 10.1016/S0959-437X(00)00254-9.View ArticlePubMedGoogle Scholar
- Sundberg S, Rydin H: Habitat requirements for establishment of Sphagnum from spores. J Ecol. 2002, 90: 268-278. 10.1046/j.1365-2745.2001.00653.x.View ArticleGoogle Scholar
- Demesure B, Comps B, Petit RJ: Chloroplast DNA phylogeography of the common beech (Fagus sylvatica L) in Europe. Evolution. 1995, 50: 2515-2520. 10.2307/2410719.View ArticleGoogle Scholar
- Heuertz M, Fineschi S, Anzidei M, Pastorelli R, Salvini D, Paule L, Frascaria-Lacoste N, Hardy OJ, Vekemans X, Vendramin GG: Chloroplast DNA variation and postglacial recolonization of common ash (Fraxinus excelsior L.) in Europe. Mol Ecol. 2004, 13: 3437-3452. 10.1111/j.1365-294X.2004.02333.x.View ArticlePubMedGoogle Scholar
- Wyatt R, Stoneburner A, Odrzykoski IJ: Bryophyte isozymes: systematic and evolutionary implications. Isozymes in Plant Biology. Edited by: Soltis DE, Soltis P. 1989, Portland: Dioscorides Press, 221-240.View ArticleGoogle Scholar
- Stenøien HK, Såstad SM: Genetic variability in bryophytes: does mating system really matter?. J Bryol. 2001, 23: 313-319.View ArticleGoogle Scholar
- McDaniel SF, Shaw AJ: Selective sweeps and intercontinental migration in the cosmopolitan moss Ceratodon purpureus (Hedw) Brid. Mol Ecol. 2005, 14: 1121-1132. 10.1111/j.1365-294X.2005.02484.x.View ArticlePubMedGoogle Scholar
- Stenøien HK, Såstad SM: Genetic structure in three haploid peat mosses (Sphagnum). Heredity. 1999, 82: 391-400. 10.1038/sj.hdy.6884940.View ArticlePubMedGoogle Scholar
- Shaw AJ: Population ecology, population genetics, and microevolution. Bryophyte Biology. Edited by: Shaw AJ, Goffinet B. 2000, Cambridge: Cambridge University Press, 369-403.View ArticleGoogle Scholar
- Szövényi P, Hock Zs, Urmi E, Schneller JJ: New primers for amplifying the GapC gene of bryophytes and its utility in infraspecific phylogenies in the genus Sphagnum. Lindbergia. 2006, 31: 32-39.Google Scholar
- Shaw AJ, Cox CJ, Boles SB: Polarity of peatmoss evolution: who says mosses have no roots?. Am J Bot. 2003, 90: 1777-1787.View ArticlePubMedGoogle Scholar
- Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.PubMed CentralView ArticlePubMedGoogle Scholar
- Clement M, Posada D, Crandall K: TCS: a computer program to estimate gene genealogies. Mol Ecol. 2000, 9: 1657-1660. 10.1046/j.1365-294x.2000.01020.x.View ArticlePubMedGoogle Scholar
- Hudson RR: Gene genealogies and the coalescent process. Oxford Surveys in Evolutionary Biology. Edited by: Harvey PH, Partridge L. 1990, New York: Oxford University Press, 7: 1-44.Google Scholar
- Hudson RR, Kaplan NL: Statistical properties of the number of recombination events in the history of a sample of DNA sequences. Genetics. 1985, 111: 147-164.PubMed CentralPubMedGoogle Scholar
- Rozas J, Sanches-Delbarrio JC, Messeguer X, Rozas R: DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics. 2003, 19: 2496-2497. 10.1093/bioinformatics/btg359.View ArticlePubMedGoogle Scholar
- HKA – a program for conducting the HKA test of natural selection. [http://lifesci.rutgers.edu/~heylab/HeylabSoftware.htm#HKA]
- McDonald JH, Kreitman M: Adaptive protein evolution at the ADH locus in Drosophila. Nature. 1991, 351: 652-654. 10.1038/351652a0.View ArticlePubMedGoogle Scholar
- Fay JC, Wu CI: Hitchhiking under positive Darwinian selection. Genetics. 2000, 155: 1405-1413.PubMed CentralPubMedGoogle Scholar
- Fu YX: Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics. 1997, 147: 915-925.PubMed CentralPubMedGoogle Scholar
- Excoffier L, Laval G, Schneider S: Arlequin ver. 3.0: an integrated software package for population genetics data analysis. Evol Bioinf. 2005, 1: 47-50.Google Scholar
- Hey J, Nielsen R: Multilocus methods for estimating population sizes, migration rates and divergence time, with applications to the divergence of Drosophila pseudoobscura and D. persimilis. Genetics. 2004, 167: 747-760. 10.1534/genetics.103.024182.PubMed CentralView ArticlePubMedGoogle Scholar
- Hey J: On the number of new world founders: a population genetic portrait of the peopling of the Americas. PLoS Biol. 2005, 3: e193-10.1371/journal.pbio.0030193.PubMed CentralView ArticlePubMedGoogle Scholar
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.