High genetic diversity at the extreme range edge: nucleotide variation at nuclear loci in Scots pine (Pinus sylvestris L.) in Scotland

Scots pine - panmixia and the elusive signal of genetic adaptation

Scots pine is the foundation species of diverse forested ecosystems across Eurasia and displays remarkable ecological breadth, occurring in environments ranging from temperate rainforests to arid tundra margins. Such expansive distributions can be favored by various demographic and adaptive processes and the interactions between them. To understand the impact of neutral and selective forces on genetic structure in Scots pine, we conducted range-wide population genetic analyses on 2321 trees from 202 populations using genotyping-by-sequencing, reconstructed the recent demography of the species and examined signals of genetic adaptation. We found a high and uniform genetic diversity across the entire range (global FST 0.048), no increased genetic load in expanding populations and minor impact of the last glacial maximum on historical population sizes. Genetic-environmental associations identified only a handful of single-nucleotide polymorphisms significantly linked to environmental gradients. The results suggest that extensive gene flow is predominantly responsible for the observed genetic patterns in Scots pine. The apparent missing signal of genetic adaptation is likely attributed to the intricate genetic architecture controlling adaptation to multi-dimensional environments. The panmixia metapopulation of Scots pine offers a good study system for further exploration into how genetic adaptation and plasticity evolve under gene flow and changing environment.

Molecular characterisation of Pinus sylvestris (L.) in Ireland at the western limit of the species distribution

BACKGROUND

Scots pine (Pinus sylvestris L.) underwent significant population declines across much of northwest Europe during the mid-to-late Holocene and was thought to have become extirpated in Ireland from about 400 AD. However, most extant populations are plantations reintroduced from Scotland. Others are naturalised therefrom and one in Western Ireland is a putative relict. In this paper, Scots pine in Ireland are genetically described for the first time.

RESULTS

Using two mitochondrial (mtDNA) loci, eight chloroplast (cpSSR) and 18 nuclear (nSSR) loci, the genetic composition and diversity of 19 Irish Scots pine populations is described and compared to other European populations. All trees sampled in Ireland were fixed for mitotype a, which is the most common across northwest Europe. By contrast, cpSSR (HCP = 0.967) and nSSR (He = 0.540) variation was high, and comparable with estimates for other regions across the species range. Differentiation at both sets of loci were similarly low (cpSSR FST = 0.019; nSSR FST = 0.018), but populations from continental Europe were significantly differentiated from all Irish populations based on nSSR variation.

CONCLUSIONS

All Irish Scots pine are likely part of a common Irish-Scottish gene pool which diverged from continental Scots pine following post-glacial recolonisation. A high genetic diversity and an absence of evidence of inbreeding suggests the regional decline of Scots pine did not critically reduce allelic variation. The post-glacial relationship between Irish and Scottish pine is discussed, and a suggestion from recent palaeoecological work that reintroduced Scots pine be managed as a native species is now further supported by genetic data.

Invasion success of a Lessepsian symbiont-bearing foraminifera linked to high dispersal ability, preadaptation and suppression of sexual reproduction

Among the most successful Lessepsian invaders is the symbiont-bearing benthic foraminifera Amphistegina lobifera. In its newly conquered habitat, this prolific calcifier and ecosystem engineer is exposed to environmental conditions that exceed the range of its native habitat. To disentangle which processes facilitated the invasion success of A. lobifera into the Mediterranean Sea we analyzed a ~ 1400 bp sequence fragment covering the SSU and ITS gene markers to compare the populations from its native regions and along the invasion gradient. The genetic variability was studied at four levels: intra-genomic, population, regional and geographical. We observed that the invasion is not associated with genetic differentiation, but the invasive populations show a distinct suppression of intra-genomic variability among the multiple copies of the rRNA gene. A reduced genetic diversity compared to the Indopacific is observed already in the Red Sea populations and their high dispersal potential into the Mediterranean appears consistent with a bridgehead effect resulting from the postglacial expansion from the Indian Ocean into the Red Sea. We conclude that the genetic structure of the invasive populations reflects two processes: high dispersal ability of the Red Sea source population pre-adapted to Mediterranean conditions and a likely suppression of sexual reproduction in the invader. This discovery provides a new perspective on the cost of invasion in marine protists: The success of the invasive A. lobifera in the Mediterranean Sea comes at the cost of abandonment of sexual reproduction.

Genetic diversity analysis and fingerprint construction of Korean pine (Pinus koraiensis) clonal seed orchard

Korean pine is a native tree species in Northeast China. In order to meet the needs of germplasm resource evaluation and molecular marker-assisted breeding of Korean pine, we collected Korean pine clones from 7 populations in Northeast China, analyzed the genetic diversity and genetic structure by SSR molecular marker technology and clustered them to revealed the inter- and intrapopulation differentiation characteristics of each clone. The fingerprint profiles of 161 Korean pine clones were also constructed. 77 alleles were detected for 11 markers, and 18 genotypes were identified on average for each marker. The PIC of the different markers ranged from 0.155-0.855, and the combination of PI and PIsibs for the 11 markers was 3.1 × 10^-8 and 1.14 × 10^-3, respectively. MANOVA showed that genetic variation existed mainly within populations, accounting for 98% of the total variation. The level of genetic differentiation among populations was low, with an average Nm between populations of 11.036. Genetic diversity is lower in the Lushuihe population and higher in the Tieli population. The 161 Korean pine clones were divided into 4 or 7 populations, and the 7 populations were not clearly distinguished from each other, with only the Lushuihe population showing partial differentiation. There is no significant correlation between the genetic distance of Korean pine populations and the geographical distance of their superior tree sources. This result can provide recommendations for future Korean pine breeding programs. The combination of 11 markers could completely distinguish 161 clones and establish the fingerprint. Genetic diversity of Korean pine clones from the 7 populations was abundant, and the genetic distances of individuals and populations were evenly dispersed. The fingerprint map can be used for the identification of Korean pine clones.

Imprints of selection in peripheral and ecologically marginal central-eastern European Scots pine populations

Knowledge of the molecular mechanisms underlying the stress response in plants is essential to understand evolutionary processes that result in long-term persistence of populations. Populations inhabiting marginal ecological conditions at the distribution range periphery may have preserved imprints of natural selection that have shaped functional genetic variation of the species. Our aim was to evaluate the extent of selection processes in the extremely fragmented, peripheral and isolated populations of Scots pine in central-eastern Europe. Autochthonous populations of the Carpathian Mts. and the Pannonian Basin were sampled and drought stress-related candidate genes were re-sequenced. Neutrality tests and outlier detection approaches were applied to infer the effect and direction of selection. Populations retained high genetic diversity by preserving a high number of alleles and haplotypes, many of them being population specific. Neutrality tests and outlier detection highlighted nucleotide positions that are under divergent selection and may be involved in local adaptation. The detected genetic pattern confirms that natural selection has played an important role in shaping modern-day genetic variation in marginal Scots pine populations, allowing for the long-term persistence of populations. Selection detected at functional regions possibly acts to maintain diversity and counteract the effect of genetic erosion.

Invasion, isolation and evolution shape population genetic structure in Campanula rotundifolia

The distribution and genetic structure of most plant species in Britain and Ireland bear the imprint of the last ice age. These patterns were largely shaped by random processes during recolonization but, in angiosperms, whole-genome duplication may also have been important. We investigate the distribution of cytotypes of Campanula rotundifolia, considering DNA variation, postglacial colonization, environmental partitioning and reproductive barriers. Cytotypes and genome size variation from across the species' range were determined by flow cytometry and genetic variation was assessed using cpDNA markers. A common garden study examined growth and flowering phenology of tetraploid, pentaploid and hexaploid cytotypes and simulated a contact zone for investigation of reproductive barriers. Irish populations were entirely hexaploid. In Britain, hexaploids occurred mostly in western coastal populations which were allopatric with tetraploids, and in occasional sympatric inland populations. Chloroplast markers resolved distinct genetic groups, related to cytotype and geographically segregated; allopatric hexaploids were distinct from tetraploids, whereas sympatric hexaploids were not. Genome downsizing occurred between cytotypes. Progeny of open-pollinated clones from the contact zone showed that maternal tetraploids rarely produced progeny of other cytotypes, whereas the progeny of maternal hexaploids varied, with frequent pentaploids and aneuploids. The presence of distinctive hexaploid chloroplast types in Ireland, Scottish islands and western mainland Britain indicates that its establishment preceded separation of these land masses by sea-level rise c. 16 000 years BP. This group did not originate from British tetraploids and probably diverged before postglacial invasion from mainland Europe. The combination of cytotype, molecular, contact zone and common garden data shows an overall pattern reflecting postglacial colonization events, now maintained by geographic separation, together with more recent occasional local in situ polyploidisation. Reproductive barriers favour the persistence of the tetraploid to the detriment of the hexaploid.

275 years of forestry meets genomics in Pinus sylvestris

Pinus sylvestris has a long history of basic and applied research that is relevant for both forestry and evolutionary studies. Its patterns of adaptive variation and role in forest economic and ecological systems have been studied extensively for nearly 275 years, detailed demography for a 100 years and mating system more than 50 years. However, its reference genome sequence is not yet available and genomic studies have been lagging compared to, for example, Pinus taeda and Picea abies, two other economically important conifers. Despite the lack of reference genome, many modern genomic methods are applicable for a more detailed look at its biological characteristics. For example, RNA-seq has revealed a complex transcriptional landscape and targeted DNA sequencing displays an excess of rare variants and geographically homogenously distributed molecular genetic diversity. Current DNA and RNA resources can be used as a reference for gene expression studies, SNP discovery, and further targeted sequencing. In the future, specific consequences of the large genome size, such as functional effects of regulatory open chromatin regions and transposable elements, should be investigated more carefully. For forest breeding and long-term management purposes, genomic data can help in assessing the genetic basis of inbreeding depression and the application of genomic tools for genomic prediction and relatedness estimates. Given the challenges of breeding (long generation time, no easy vegetative propagation) and the economic importance, application of genomic tools has a potential to have a considerable impact. Here, we explore how genomic characteristics of P. sylvestris, such as rare alleles and the low extent of linkage disequilibrium, impact the applicability and power of the tools.

Exome Resequencing Reveals Evolutionary History, Genomic Diversity, and Targets of Selection in the Conifers Pinus taeda and Pinus elliottii

Loblolly pine (Pinus taeda) and slash pine (Pinus elliottii) are ecologically and economically important pine species that dominate many forest ecosystems in the southern United States, but like all conifers, the study of their genetic diversity and demographic history has been hampered by their large genome size. A small number of studies mainly based on candidate-gene sequencing have been reported for P. taeda to date, whereas none are available for P. elliottii. Targeted exome resequencing has recently enabled population genomics studies for conifers, approach used here to assess genomic diversity, signatures of selection, population structure, and demographic history of P. elliottii and P. taeda. Extensive similarities were revealed between these species: both species feature rapid linkage disequilibrium decay and high levels of genetic diversity. Moreover, genome-wide positive correlations for measures of genetic diversity between the species were also observed, likely due to shared structural genomic constraints. Also, positive selection appears to be targeting a common set of genes in both pines. Demographic history differs between both species, with only P. taeda being affected by a dramatic bottleneck during the last glacial period. The ability of P. taeda to recover from a dramatic reduction in population size while still retaining high levels of genetic diversity shows promise for other pines facing environmental stressors associated with climate change, indicating that these too may be able to adapt successfully to new future conditions even after a drastic population size contraction.

Natural variation of DNA methylation and gene expression may determine local adaptations of Scots pine populations

Long-lived conifers are vulnerable to climate change because classical evolutionary processes are slow in developing adaptive responses. Therefore, the capacity of a genotype to adopt different phenotypes is important. Gene expression is the primary mechanism that converts genome-encoded information into phenotypes, and DNA methylation is employed in the epigenetic regulation of gene expression. We investigated variations in global DNA methylation and gene expression between three Scots pine (Pinus sylvestris L.) populations located in northern and southern Finland using mature seeds. Gene expression levels were studied in six DNA methyltransferase (DNMT) genes, which were characterized in this study, and in 19 circadian clock genes regulating adaptive traits. In embryos, expression diversity was found for three DNMT genes, which maintain DNA methylation. The expression of two DNMT genes was strongly correlated with climate variables, which suggests a role for DNA methylation in local adaptation. For adaptation-related genes, expression levels showed between-population variation in 11 genes in megagametophytes and in eight genes in embryos, and many of these genes were linked to climate factors. Altogether, our results suggest that differential DNA methylation and gene expression contribute to local adaptation in Scots pine populations and may enhance the fitness of trees under rapidly changing climatic conditions.

Cryptic genetic variation and adaptation to waterlogging in Caledonian Scots pine, Pinus sylvestris L

Local adaptation occurs as the result of differential selection among populations. Observations made under common environmental conditions may reveal phenotypic differences between populations with an underlying genetic basis; however, exposure to a contrasting novel environment can trigger release of otherwise unobservable (cryptic) genetic variation. We conducted a waterlogging experiment on a common garden trial of Scots pine, Pinus sylvestris (L.), saplings originating from across a steep rainfall gradient in Scotland. A flood treatment was maintained for approximately 1 year; physiological responses were gauged periodically in terms of photochemical capacity as measured via chlorophyll fluorescence. During the treatment, flooded individuals experienced a reduction in photochemical capacity, F v /F m, this reduction being greater for material originating from drier, eastern sites. Phenotypic variance was increased under flooding, and this increase was notably smaller in saplings originating from western sites where precipitation is substantially greater and waterlogging is more common. We conclude that local adaptation has occurred with respect to waterlogging tolerance and that, under the flooding treatment, the greater increase in variability observed in populations originating from drier sites is likely to reflect a relative absence of past selection. In view of a changing climate, we note that comparatively maladapted populations may possess considerable adaptive potential, due to cryptic genetic variation, that should not be overlooked.

Planting exotic relatives has increased the threat posed by Dothistroma septosporum to the Caledonian pine populations of Scotland

To manage emerging forest diseases and prevent their occurrence in the future, it is essential to determine the origin(s) of the pathogens involved and identify the management practices that have ultimately caused disease problems. One such practice is the widespread planting of exotic tree species within the range of related native taxa. This can lead to emerging forest disease both by facilitating introduction of exotic pathogens and by providing susceptible hosts on which epidemics of native pathogens can develop. We used microsatellite markers to determine the origins of the pathogen Dothistroma septosporum responsible for the current outbreak of Dothistroma needle blight (DNB) on native Caledonian Scots pine (Pinus sylvestris) populations in Scotland and evaluated the role played by widespread planting of two exotic pine species in the development of the disease outbreak. We distinguished three races of D. septosporum in Scotland, one of low genetic diversity associated with introduced lodgepole pine (Pinus contorta), one of high diversity probably derived from the DNB epidemic on introduced Corsican pine (Pinus nigra subsp. laricio) in England and a third of intermediate diversity apparently endemic on Caledonian Scots pine. These races differed for both growth rate and exudate production in culture. Planting of exotic pine stands in the UK appears to have facilitated the introduction of two exotic races of D. septosporum into Scotland which now pose a threat to native Caledonian pines both directly and through potential hybridization and introgression with the endemic race. Our results indicate that both removal of exotic species from the vicinity of Caledonian pine populations and restriction of movement of planting material are required to minimize the impact of the current DNB outbreak. They also demonstrate that planting exotic species that are related to native species reduces rather than enhances the resilience of forests to pathogens.

Early phenology and growth trait variation in closely related European pine species

Closely related taxa occupying different environments are valuable systems for studying evolution. In this study, we examined differences in early phenology (bud set, bud burst) and early growth in a common garden trial of closely related pine species: Pinus sylvestris, P. mugo, and P. uncinata. Seeds for the trial were sourced from populations across the ranges of each species in Europe. Over first 4 years of development, clear differences were observed between species, while the most significant intraspecific differentiation was observed among plants from P. sylvestris populations from continental European locations. Trait differences within P. sylvestris were highly correlated with altitude and latitude of the site of origin. Meanwhile, P. mugo populations from the Carpathians had the earliest bud set and bud flush compared to other populations of the species. Overall, populations from the P. mugo complex from heterogeneous mountain environments and P. sylvestris from the Scottish Highlands showed the highest within-population variation for the focal traits. Although the three species have been shown to be genetically highly similar, this study reveals large differences in key adaptive traits both among and within species.

Substantial variation in the timing of pollen production reduces reproductive synchrony between distant populations of Pinus sylvestris L. in Scotland

The ability of a population to genetically adapt to a changing environment is contingent not only on the level of existing genetic variation within that population, but also on the gene flow received from differently adapted populations. Effective pollen-mediated gene flow among plant populations requires synchrony of flowering. Therefore differences in timing of flowering among genetically divergent populations may reduce their ability to adapt to environmental change. To determine whether gene flow among differently adapted populations of native Scots pine (Pinus sylvestris) in Scotland was restricted by differences in their flowering phenology, we measured timing of pollen release among populations spanning a steep environmental gradient over three consecutive seasons (2014-2016). Results showed that, over a distance of 137 km, there were as many as 15.8 days' difference among populations for the predicted timing of peak pollen shedding, with the earliest development in the warmer west of the country. There was much variation between years, with the earliest development and least synchrony in the warmest year (2014) and latest development and greatest synchrony in the coolest year (2015). Timing was negatively correlated with results from a common-garden experiment, indicative of a pattern of countergradient variation. We conclude that the observed differences in reproductive synchrony were sufficient to limit gene flow via pollen between populations of P. sylvestris at opposite ends of the environmental gradient across Scotland. We also hypothesize that continually warming, or asymmetrically warming spring temperatures will decrease reproductive synchrony among pine populations.

Contrasting patterns of genetic variation in core and peripheral populations of highly outcrossing and wind pollinated forest tree species

Gene flow tends to have a homogenising effect on a species' background genetic variation over large geographical areas. However, it is usually unknown to what extent the genetic structure of populations is influenced by gene exchange between core and peripheral populations that may represent stands of different evolutionary and demographic history. In this study, we looked at the patterns of population differentiation in Scots pine-a highly outcrossing and wind pollinated conifer species that forms large ecosystems of great ecological and economic importance in Europe and Asia. A set of 13 polymorphic nuclear microsatellite loci was analysed to infer the genetic relationships among 24 populations (676 individuals) from Europe and Asia Minor. The study included specimens from the primary continuous range and from isolated, marginal stands that are considered to be autochthonous populations representative of the species' putative refugial areas. Despite their presumably different histories, a similar level of genetic variation and no evidence of a population bottleneck was found across the populations. Differentiation among populations was relatively low (average FST = 0.035); however, the population structure was not homogenous, which was clearly evident from the allelic frequency spectra and Bayesian assignment analysis. Significant differentiation over short geographical distances was observed between isolated populations within the Iberian and Anatolian Peninsulas (Asia Minor), which contrasted with the absence of genetic differentiation observed between distant populations e.g., between central and northern Europe. The analysed populations were assigned to several groups that corresponded to the geographical regions of their occurrence. These results will be useful in genetics studies in Scots pine that aim to link nucleotide and phenotypic variation across the species distribution range and for development of sustainable breeding and management programs.

Has Scots pine (Pinus sylvestris) co-evolved with Dothistroma septosporum in Scotland? Evidence for spatial heterogeneity in the susceptibility of native provenances

Spatial heterogeneity in pathogen pressure leads to genetic variation in, and evolution of, disease-related traits among host populations. In contrast, hosts are expected to be highly susceptible to exotic pathogens as there has been no evolution of defence responses. Host response to pathogens can therefore be an indicator of a novel or endemic pathosystem. Currently, the most significant threat to native British Scots pine (Pinus sylvestris) forests is Dothistroma needle blight (DNB) caused by the foliar pathogen Dothistroma septosporum which is presumed to be exotic. A progeny-provenance trial of 6-year-old Scots pine, comprising eight native provenances each with four families in six blocks, was translocated in April 2013 to a clear-fell site in Galloway adjacent to a DNB-infected forest. Susceptibility to D. septosporum, measured as DNB severity (estimated percentage nongreen current-year needles), was assessed visually over 2 years (2013-2014 and 2014-2015; two assessments per year). There were highly significant differences in susceptibility among provenances but not among families for each annual assessment. Provenance mean susceptibility to D. septosporum was negatively and significantly associated with water-related variables at site of origin, potentially due to the evolution of low susceptibility in the host in response to high historical pathogen pressure.

Substantial heritable variation for susceptibility to Dothistroma septosporum within populations of native British Scots pine (Pinus sylvestris)

The threat from pests and pathogens to native and commercially planted forest trees is unprecedented and expected to increase under climate change. The degree to which forests respond to threats from pathogens depends on their adaptive capacity, which is determined largely by genetically controlled variation in susceptibility of the individual trees within them and the heritability and evolvability of this trait. The most significant current threat to the economically and ecologically important species Scots pine (Pinus sylvestris) is dothistroma needle blight (DNB), caused by the foliar pathogen Dothistroma septosporum. A progeny-population trial of 4-year-old Scots pine trees, comprising six populations from native Caledonian pinewoods each with three to five families in seven blocks, was artificially inoculated using a single isolate of D. septosporum. Susceptibility to D. septosporum, assessed as the percentage of non-green needles, was measured regularly over a period of 61 days following inoculation, during which plants were maintained in conditions ideal for DNB development (warm; high humidity; high leaf wetness). There were significant differences in susceptibility to D. septosporum among families indicating that variation in this trait is heritable, with high estimates of narrow-sense heritability (0.38-0.75) and evolvability (genetic coefficient of variation, 23.47). It is concluded that native Scots pine populations contain sufficient genetic diversity to evolve lower susceptibility to D. septosporum through natural selection in response to increased prevalence of this pathogen.

Comparative transcriptomics of a complex of four European pine species

Background

Pinus sylvestris, P. mugo, P. uliginosa and P. uncinata are closely related but phenotypically and ecologically very distinct European pine species providing an excellent study system for analysis of the genetic basis of adaptive variation and speciation. For comparative genomic analysis of the species, transcriptome sequence was generated for 17 samples collected across the European distribution range using Illumina paired-end sequencing technology.

Results

De novo transcriptome assembly of a reference sample of P. sylvestris contained 40968 unigenes, of which fewer than 0.5% were identified as putative retrotransposon sequences. Based on gene annotation approaches, 19659 contigs were identified and assigned to unique genes covering a broad range of gene ontology categories. About 80% of the reads from each sample were successfully mapped to the reference transcriptome of P. sylvestris. Single nucleotide polymorphisms were identified in 22041-24096 of the unigenes providing a set of ~220-262 k SNPs identified for each species. Very similar levels of nucleotide polymorphism were observed across species (π=0.0044-0.0053) and highest pairwise nucleotide divergence (0.006) was found between P. mugo and P. sylvestris at a common set of unigenes.

Conclusions

The study provides whole transcriptome sequence and a large set of SNPs to advance population and association genetic studies in pines. Our study demonstrates that transcriptome sequencing can be a very useful approach for development of novel genomic resources in species with large and complex genomes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1401-z) contains supplementary material, which is available to authorized users.

Signatures of volcanism and aridity in the evolution of an insular pine (Pinus canariensis Chr. Sm. Ex DC in Buch)

Oceanic islands of volcanic origin provide useful templates for the study of evolution because they are subjected to recurrent perturbations that generate steep environmental gradients that may promote adaptation. Here we combine population genetic data from nuclear genes with the analysis of environmental variation and phenotypic data from common gardens to disentangle the confounding effects of demography and selection to identify the factors of importance for the evolution of the insular pine P. canariensis. Eight nuclear genes were partially sequenced in a survey covering the entire species range, and phenotypic traits were measured in four common gardens from contrasting environments. The explanatory power of population substrate age and environmental indices were assessed against molecular and phenotypic diversity estimates. In addition, neutral genetic variability (FST) and the genetic differentiation of phenotypic variation (QST) were compared in order to identify the evolutionary forces acting on these traits. Two key factors in the evolution of the species were identified: (1) recurrent volcanic activity has left an imprint in the genetic diversity of the nuclear genes; (2) aridity in southern slopes promotes local adaptation in the driest localities of P. canariensis, despite high levels of gene flow among populations.

Inferring multiple refugia and phylogeographical patterns in Pinus massoniana based on nucleotide sequence variation and DNA fingerprinting

BACKGROUND

Pinus massoniana, an ecologically and economically important conifer, is widespread across central and southern mainland China and Taiwan. In this study, we tested the central-marginal paradigm that predicts that the marginal populations tend to be less polymorphic than the central ones in their genetic composition, and examined a founders' effect in the island population.

METHODOLOGY/PRINCIPAL FINDINGS

We examined the phylogeography and population structuring of the P. massoniana based on nucleotide sequences of cpDNA atpB-rbcL intergenic spacer, intron regions of the AdhC2 locus, and microsatellite fingerprints. SAMOVA analysis of nucleotide sequences indicated that most genetic variants resided among geographical regions. High levels of genetic diversity in the marginal populations in the south region, a pattern seemingly contradicting the central-marginal paradigm, and the fixation of private haplotypes in most populations indicate that multiple refugia may have existed over the glacial maxima. STRUCTURE analyses on microsatellites revealed that genetic structure of mainland populations was mediated with recent genetic exchanges mostly via pollen flow, and that the genetic composition in east region was intermixed between south and west regions, a pattern likely shaped by gene introgression and maintenance of ancestral polymorphisms. As expected, the small island population in Taiwan was genetically differentiated from mainland populations.

CONCLUSIONS/SIGNIFICANCE

The marginal populations in south region possessed divergent gene pools, suggesting that the past glaciations might have low impacts on these populations at low latitudes. Estimates of ancestral population sizes interestingly reflect a recent expansion in mainland from a rather smaller population, a pattern that seemingly agrees with the pollen record.