Higher fine-scale genetic structure in peripheral than in core populations of a long-lived and mixed-mating conifer - eastern white cedar (Thuja occidentalis L.)

Background Fine-scale or spatial genetic structure (SGS) is one of the key genetic characteristics of plant populations. Several evolutionary and ecological processes and population characteristics influence the level of SGS within plant populations. Higher fine-scale genetic structure may be expected in peripheral than core populations of long-lived forest trees, owing to the differences in the magnitude of operating evolutionary and ecological forces such as gene flow, genetic drift, effective population size and founder effects. We addressed this question using eastern white cedar (Thuja occidentalis) as a model species for declining to endangered long-lived tree species with mixed-mating system. Results We determined the SGS in two core and two peripheral populations of eastern white cedar from its Maritime Canadian eastern range using six nuclear microsatellite DNA markers. Significant SGS ranging from 15 m to 75 m distance classes was observed in the four studied populations. An analysis of combined four populations revealed significant positive SGS up to the 45 m distance class. The mean positive significant SGS observed in the peripheral populations was up to six times (up to 90 m) of that observed in the core populations (15 m). Spatial autocorrelation coefficients and correlograms of single and sub-sets of populations were statistically significant. The extent of within-population SGS was significantly negatively correlated with all genetic diversity parameters. Significant heterogeneity of within-population SGS was observed for 0-15 m and 61-90 m between core and peripheral populations. Average Sp, and gene flow distances were higher in peripheral (Sp = 0.023, σg = 135 m) than in core (Sp = 0.014, σg = 109 m) populations. However, the mean neighborhood size was higher in the core (Nb = 82) than in the peripheral (Nb = 48) populations. Conclusion Eastern white cedar populations have significant fine-scale genetic structure at short distances. Peripheral populations have several-folds higher within-population fine-scale genetic structure than core populations. Anthropogenic disturbances and population fragmentation presumably have significant effects on fine-scale genetic structure in eastern white cedar. Core populations have higher neighborhood size than peripheral populations, whereas gene flow distances are higher in peripheral than in core populations. The results of our study contribute to the knowledge of poorly-understood spatial genetic structure of core versus peripheral populations in plants. As well, the information is of significance for conservation of genetic resources of eastern white cedar and perhaps of other long-lived forest trees with mixed-mating system.


Background
Recently, in addition to core, peripheral populations are also gaining importance for gene conservation, because under global climate change conditions they may possess genotypes of future adaptive potential. Plant populations tend to migrate further from their range periphery to ensure their survival by developing adaptability in a changing environmental condition due to global warming [1]. Genetic studies involving both core and peripheral populations are still limited [2]. Additionally, the available studies mainly focused on withinand among-population genetic diversity [2][3][4], and majority of them reported higher genetic diversity and lower population differentiation in core than in peripheral populations. Although still limited, there are increasing reports on SGS in fragmented versus continuous populations e.g., [5][6][7]. However, information on the comparison of SGS between core and peripheral populations of plant species is almost non-existent. The only published report on this aspect by Gapare & Aitken [8] reported significantly higher SGS in peripheral and random distribution of genotypes in core populations of predominantly outcrossing conifer, Sitka spruce (Picea sitchensis). The level of SGS within plant populations depends on seed and pollen dispersal, mating system, mode of reproduction, size of the founder population, species longevity, adult population density, effective population size, and micro environmental selection [9][10][11][12][13]. Anthropogenic disturbance and fragmentation of populations can also influence SGS by restricting gene flow and increasing local drift due to severe modification of forest environment and ecosystem, and density of populations [14,15]. Information on SGS levels in core and peripheral populations can unravel differential effect of the key evolutionary and ecological genetic processes operating in natural plant populations, such as gene dispersal, natural selection and genetic drift. Analysis of SGS is also important for understanding the population dynamics of a species as its response to microenvironmental variation and global climate change [16,17]. Furthermore, knowledge of SGS in natural populations is also important for developing scientifically-sound sampling strategies for genetic resource conservation and breeding purposes.
Higher SGS may be expected in peripheral than in core populations of long-lived plants, such as forest trees, owing to the differences in founding population, seed and pollen dispersal, adult plant density, inbreeding levels, effective population size (N e ), and disturbance and selection regimes. Restricted gene flow due to limited seed and pollen dispersal can promote mating among relatives which can result in isolation by distance within a population [18]. Therefore, a local genetic subdivision can develop SGS within populations. Selfing and other forms of inbreeding can limit gene dispersal, consequently influencing the spatial distribution of genetic variation within populations [16]. In contrast, generally, outcrossing decreases population subdivision by reducing mating among relatives and increasing N e [10]. Therefore, higher level of SGS can be expected in inbreeding species as compared to outcrossing species. Vegetative clonal reproduction can significantly affect the SGS through the development of spatial clusters of the members of a clone [19]. Long-lived plants are expected to have low levels of SGS within their populations because of decreased population divergence due to increased N e and decreased genetic drift over time [10]. SGS has been found to be inversely correlated with adult density within populations [20]. Additionally, the levels of SGS were found to be positively correlated with the level of fragmentation of plant populations e.g., [6,14,21].
Conifers are generally long-lived plants that are ecologically and economically important globally. They constitute a major component of the northern boreal and temperate forests, especially dominating the boreal forest which comprises one of the largest biomes in the world http://www.borealforest.org/index.php?category=-world_boreal_forest&page=overview. Most of the conifer tree species are characterized by predominantly (~90%) outcrossing mating system. Studies in conifers reported either no or weak levels of SGS within populations [21][22][23][24][25][26][27][28][29][30][31][32]. This is consistent with the conifers' characteristics of predominant outcrossing, wind-pollination and seed dispersal by wind. Except for two studies [24,32], all of the SGS studies in conifers are reported for Pinus or Picea species. Furthermore, with the exception of two [8,28], all other studies on SGS in conifers are based on only one or two populations. All of the SGS studies in conifers have been reported for predominantly outcrossing (≥90%) species, and there is no SGS information on plant species that have mixed-mating system with high selfing levels.
Eastern (or northern) white cedar (Thuja occidentalis L.) provides an ideal model species to study SGS in core versus peripheral populations of declining and longlived plant species with mixed mating system. Eastern white cedar is a wind-pollinated, monoecious, evergreen tree species belonging to the Cypress family (Cupressaceae). It reproduces both sexually and vegetatively. Vegetative reproduction occurs mainly in swampy sites by layering [33]. Eastern white cedar is an ecologically and economically important tree species, highly suitable for many wood products. Since its foliage is a favorite winter food for deer, it is very important for wildlife conservation. It can reach an age of 400 years [34]. populations [34,35]. Historical records indicate that the range and number of populations of eastern white cedar have diminished due to logging, and clearing for agriculture over the past 400 years in Nova Scotia and PEI [33,35]. Presently, there are only 32 identified populations of eastern white cedar in Nova Scotia [34], which are highly fragmented. Due to its rarity in Nova Scotia, eastern white cedar has been designated as a Red (endangered) species in this province. In Prince Edward Island, occurrence of eastern white cedar is limited to only a few populations, mainly found in the western part of the Island (West of Summerside). Eastern white cedar is common in New Brunswick, where it forms its best development [35]. Thus, New Brunswick represents the core or central range of this species. However, here also eastern white cedar has been listed as a species of concern. Like other cedar species (Thuja plicata, T. orientalis, [36][37][38]), eastern white cedar also possesses relatively low genetic diversity [39][40][41]. It is characterized by a mixed-mating system with high levels (mean 54%) of selfing rates [41,42]. Although a few genetic diversity studies have been conducted in the genus Thuja [36][37][38][39][40][41]43], to our best knowledge, there is no information available on SGS in any Thuja species. The only SGS study in the Cupressaceae family was done for one population of an outcrossing species, Cryptomeria japonica [32], which reported significant SGS at short distance classes (up to 60 m) based on nuclear microsatellites, but non-significant SGS for chloroplast microsatellites.
The aim of the present study was to determine withinpopulation fine-scale genetic structure of eastern white cedar growing in core and peripheral populations. Our hypotheses were: (i) since eastern white cedar is a mixed-mating long-lived tree with high levels of selfing rate, and has capability to reproduce vegetatively, this species is expected to have some levels of SGS within population; (ii) the peripheral populations of eastern white cedar are expected to have higher levels of SGS than the core populations due to their lower population density; potentially restricted gene dispersal, smaller N e , and higher anthropogenic fragmentation as compared to the core populations [9][10][11][12][13].

Populations and sampling
Four natural populations of eastern white cedar were sampled from two Maritime Provinces of Canada. Two of these were core (Edmundston and Perth Andover) populations from New Brunswick and other two were peripheral (Cedar Lake and Rockland) populations from Nova Scotia ( Figure 1; Table 1). We restricted sampling of the core and peripheral populations to the same Maritime region of eastern Canada to avoid other regional and ecological effects confounding the results. The two populations sampled in New Brunswick are located in the central part of the species' range. Despite some harvesting activities, eastern white cedar populations remain intact with contiguous distribution in this province with no barrier to gene flow between populations. On the other hand, the two eastern white cedar populations sampled from Nova Scotia are disjunct, and highly fragmented and disturbed as they have experienced anthropogenic disturbances over a long period of time. It was not possible to find undisturbed pristine natural populations of eastern white cedar in both core and peripheral regions in eastern Canada. Also, fragmented core and contiguous peripheral populations could not be found in our study region. All of the sampled populations are over 100 years old. The populations sampled in Nova Scotia were from almost pure eastern white cedar stands, whereas the populations sampled from New Brunswick are mixed with other tree species, such as black spruce (Picea mariana) and balsam fir (Abies balsamea).
A plot consisting of~100 contiguous eastern white cedar trees was established in each of the four sampled populations irrespective of the plot area ( Figure 2). All of the eastern white cedar individuals were sampled within these plots. We chose this sampling strategy to keep the sample size equal among all four sampled populations, in order to minimize the effects of sample size confounding the comparison of SGS results between the populations. Spatial locations of each of the sampled trees were recorded as their geo-coordinates using a Geographical Positioning System (Garmin GPS MAP 60CSx). The spatial location maps of the sampled trees in each of the four populations are provided in Figure 2.
Small amounts of foliage samples were collected from each of the 394 sampled trees in plastic Ziplock bags containing silica gel. Upon arrival to the lab, the foliage samples were stored at -20°C until processed for DNA extraction.

DNA extraction and genotyping
Total genomic DNA was extracted from the needle samples of individual trees using the DNeasy Plant Mini Kit (Qiagen, USA). Microsatellite DNA markers of the nuclear genome were used for genotyping the sampled trees. A total of 40 microsatellite markers developed for the Cupressaceae species were tested in eastern white cedar. Of these, 37 were developed for Thuja plicata [44,45], [L.M. O'Connell & C.E. Ritland, Unpublished data; C.J. Newton, Unpublished data] and three for Chameacyparis obtusa [46]. Only six (TP1, TP2, TP4, TP6, TP11 and 91 G) of the 40 microsatellites yielded simple and scorable PCR (Polymerase Chain Reaction) products. These six microsatellites were optimized and used for genotyping individual eastern white cedar trees. PCR amplifications were carried out using 10 μl reaction mixtures containing 20 to 30 ng of genomic DNA, 1× PCR buffer, 1.5 mM MgCl 2 , 0.20 mM of each dNTP, 0.50 pmol of each primer and 0.25 U Taq DNA polymerase (Invitrogen, USA). The epGradients Master Cycles (Eppendorf, Germany) were used for PCR amplification of the sample DNA, with the following thermocycling profile: initial denaturation step at 95°C for 5 min, followed by 33 cycles each of 95°C for 45 sec, annealing temperature at 55 to 65°C for 45 sec, 72°C for 45 sec and final extension step of 72°C for 5 min. PCR products were separated on LI-COR4300 genetic analyzers (LI-COR Inc. Lincoln, Nebraska, USA). The genotypes of individual trees were first determined by scoring the LI-COR gels with the SAGA GT program  (LI-COR Inc. Lincoln, Nebraska, USA) and later verified by manual scoring.

Statistical data analysis Genetic diversity and N e
The standard population genetic diversity parameters (number of alleles, observed and expected heterozygosity) and fixation index (F) were calculated using the program GENALEX [47]. Allelic richness (A R ) was estimated using the program FSTAT [48]. Genotypic diversity was determined by calculating observed and expected genotype additivity (GA O and GA E ) as the sum of genotypes observed and expected under Hardy-Weinberg equilibrium over the loci [49]. Deviations of genotypic frequencies from the Hardy-Weinberg expectation were tested using the Markov chain method (parameter values: dememorisation number = 10,000, number of batches = 300, number of iterations per batch = 5,000) in the GENEPOP program [50].
The effective population size (N e ) for each of the four populations was estimated using the maximum likelihood method of coalescent-theory and MIGRATE program [51]. At first, the θ values were estimated for each population using MIGRATE program and assuming an average microsatellite mutation rate of 10 -3 per generation. The N e was calculated as θ/(4×10. -3 ).

Spatial distribution
Pattern of spatial distribution of individual trees in the four eastern white cedar populations were estimated using univariate second-order Ripley's [52]K-function for 10, 15, 20 and 30 m distance classes (t). The weighted edge correction [53] method was applied to correct the edge effect. Monte Carlo procedure was used to test the significance of deviation from a random distribution at 95% upper and lower confidence envelops, generating 10,000 replications. These calculations were performed using the SPPACK program [54].

Spatial Genetic Structure (SGS)
We used Smouse and Peakall's [12] spatial autocorrelation method to estimate the SGS in the four eastern white cedar populations employing the GENALEX 6.2 [47] program. This method determines spatial distribution of genetic variation based on a frequencyweighted average over all alleles and loci and correlation coefficients of genotypes in relation to their spatial distance. It allows for the covariance between alleles (see [12] for detailed description). In order to facilitate the comparison of SGS among the four populations, we used 15 m distance class as representative because three, out of the four populations showed maximum values of Ripley's K at this distance class. Similar approach was used by Oddou-Muratorio et al. [55] to select the distance class for SGS analysis in a tree species, Sorbus torminalis.
In addition to the estimation of multilocus SGS parameters for each of the four populations separately, SGS analyses were also performed for combined subset each of core (ED, and PA), and peripheral (CL and RL) populations, as well by combining all four populations to infer degrees of heterogeneity of SGS within and between subsets based on the method of Smouse et al. [56] using GENALEX 6.2 program. These analyses permit us to formally test whether differences in SGS between the subsets of populations are statistically significant. Since the analysis of SGS using combined multiple populations increases sample sizes, it increases statistical power of the analysis to detect the presence of SGS within populations. Heterogeneity tests for single and multi-populations correlograms were performed as described by Smouse et al. [56]. This method first performs the heterogeneity test of autocorrelation coefficient among individuals lag-by-lag using the pooled within-population correlograms. Then the single-lag tests are used to perform an analogous heterogeneity test of multilag correlograms (see Smouse et al. [56] for detailed description). Single class (t 2 ) and multi-class (ω) test criteria were used to test the heterogeneity of spatial autocorrelation (r) values between single and subsets of populations. Heterogeneity in SGS for single distance classes between any two populations or a subset of populations was determined by performing the squared paired-sample t test (t 2 ) [56]. SGS heterogeneity over all distance classes (a) in a population and (b) among populations was determined by using the ω test criteria [56] under the null hypotheses that (a) there is no significant spatial autocorrelation at any distance class, and that (b) the spatial correlograms from all populations are homogeneous. Statistical significance of all the analyses was performed by permutation test with 9,999 replications at the 95% upper and lower confidence intervals.
The correlations of the SGS parameters (significant spatial autocorrelation coefficient, r) with N e (actual or log transformed), genetic diversity parameters and density of eastern white cedar trees within the sampled populations were determined by calculating the Pearson correlation coefficients (R), using the MINITAB 15 package.

Sp statistics and dispersal estimates
In order to quantify SGS and to compute neighborhood sizes and gene dispersal distances in the four eastern white cedar populations, Sp statistic was estimated as follows: Sp = -^b F /(1-^F 1 ), where -^b F is the regression slope of the autocorrelogram and^F 1 is the mean kinship coefficient (Fij) [57] between individuals belonging to the first (0-15 m) distance class [21]. Using the Sp statistics, Wright's neighborhood size Nb (Nb = 1/Sp) was also estimated. Nb is defined as Nb = 4πD e σ g 2 , where σ g 2 is the axial variance of gene dispersal and D e the effective population density. Estimates of gene dispersal distance (σg) were determined assuming ratios of effective to census density, De/D, of 0.1 and 0.5 [7,21]. Estimates of -^b F and^F 1 were carried out using SPA-GEDI program [58]. Statistical significance (P ≤ 0.05) of the analyses was performed by permutation tests with 10,000 replications.

Genetic diversity and effective population size
In total, 95 alleles were observed at six microsatellite loci in four eastern white cedar populations. The core populations showed higher allelic as well as genotypic genetic diversity than the peripheral populations (Table  2). Also, the core populations possessed higher number of private alleles than the peripheral ones ( Table 2). The two core populations from New Brunswick (ED and PA) showed very similar levels of genetic diversity. All four populations showed significant deviation from the Hardy-Weinberg equilibrium. The N e was 432 and 910, with an average of 671 in the core populations, and 110 and 310, with an average of 210 in the peripheral populations ( Table 2). As the objective of this study was to determine SGS, we do not focus on discussing genetic diversity results here. Results on genetic diversity and population structure of 13 eastern white cedar populations from the same region are presented and discussed elsewhere [59].

Spatial distribution
Results of spatial distribution pattern analysis using Ripley's K-function and the number of pairs of observation for 15 m distance class are presented in Figure 3, and Table 3 (Table 4). An overall significance test (ω) of spatial autocorrelations and correlograms revealed that all single or combinedpopulations were statistically significant ( Table 4).
The maximum values of spatial autocorrelation coefficient (r) were observed for the first (0-15 m) distance class for individual or combined populations ( Table 4). The r values for the first distance class showed great variation among the four populations ranging from 0.030 (PA population) to 0.114 (RL population) ( Figure  4). Thus, r in RL was about four times of that in PA. The highest r for the combined peripheral populations (r = 0.093) was considerably higher than in the intact populations (r = 0.058).
Except for the population pair CL and RL, all pairs revealed significant heterogeneity in SGS between the populations based on the multi-class test criterion (ω) ( Table 4). Significant heterogeneity in SGS was also observed between the core and peripheral populations  Table 4).

Sp statistics and dispersal estimates
Results of Sp statistics and dispersal estimates in the four eastern white cedar populations are presented in Table 5.  Table 5).
Correlations of SGS parameters with genetic diversity, stand density, and population Ne

Spatial genetic structure in eastern white cedar
Our study demonstrates that eastern white cedar has significant within-population fine-scale genetic structure at short distances. This was clearly evident from significant positive SGS at~15 m to~75 m in four studied eastern white cedar populations; suggesting that  individuals with similar alleles or genotypes were more often growing in proximity than expected by chance alone. The extent of the observed significant SGS depended upon the location of the populations (core versus peripheral). These results are discussed later in this section. The extent of SGS for combined four eastern white cedar populations (up to 45 m) is more than two times of that observed in one 20-years old population of another mixed-mating (mean selfing rate of 27% [60]) species Larix larcinia (up to 20 m) [33]. The existence of SGS in Larix larcinia was explained primarily due to the effect of founding population. In contrast, predominantly-outcrossing species showed either weak or no within-population SGS e.g. [17,22,28]. Thus, mixed-mating conifer species tolerating high levels of selfing apparently have higher SGS within their populations than the predominantly outcrossing conifers.
Besides higher levels of selfing or inbreeding, limited gene flow and founding gene pool may be responsible for the observed short-distance SGS in eastern white cedar. Restricted gene flow due to limited pollen and seed dispersal contributes to local grouping of individuals with similar alleles or genotypes [10]. In eastern white cedar, both pollen and seeds are disseminated by wind. We are not aware of any information on the average distance of pollen dispersal in this species. However, in conifers it is well known that pollen is dispersed in enormous amounts over considerable distances [61,62]. Therefore, the restriction in gene flow due to limited pollen flow is less likely to happen in eastern white cedar, especially in the contiguous populations. However, in the peripheral populations, reduced gene flow may have contributed to the observed levels of fine-scale genetic structure as a result of fragmentation. Fowells [63] reported a seeding range of up to 60 m in eastern white cedar under normal conditions. Interestingly, the mean level of SGS observed in eastern white cedar (up to 45 m) falls within this range of the seeding distance. Therefore, restricted gene flow due to limited seed dispersal may be one of the main reasons for the observed short-distance SGS in the species.
Vegetative propagation through natural layering has been reported in eastern white cedar [33,64]. Since the individuals established from vegetative propagation have genetic constitution identical to the source tree and grow in proximity forming clumps, this may generate a significant SGS at short distance intervals. However, we did not find any concrete evidence for vegetative propagation contributing to the observed levels of SGS in any of the four study populations. We detected only four to eight individuals (4 to 8%) per population growing in proximity of each other with the same multilocus genotypes (data not shown). This level of multilocus genotype sharing may be due to vegetative reproduction, selfing, a combination of these factors or just by chance. However, our data does not support vegetative propagation contributing to the observed levels of SGS in eastern white cedar. Although, it is generally expected that vegetative reproduction is one of the major factors that may contribute to the occurrence of significant SGS in vegetatively reproducing species, Chung & Epperson [65] found little effect of vegetative reproduction on SGS in Eurya emarginata. They explained that the restriction in gene flow due to limited dispersal of seed and pollen could be the main reason for the observed SGS in the species.

Spatial distribution pattern of trees
The analysis of spatial distribution pattern using Ripley's K function indicated significant aggregation of trees in all four populations of eastern white cedar. Within population spatial patterns of individual plants can be influenced by several factors such as variation in microsites, level of disturbance, seed dispersal patterns [56,66]. Parker et al. [66] and Trabaud et al. [67] reported a significant aggregation of tress in Pinus clausa and Pinus halepensis populations, respectively, and they explained that the observed aggregation could be due to either limited seed dispersal or recruitment of  on the studied populations cannot be ruled out because fire is one of the main agents of natural disturbance in boreal forest [68].

SGS in core versus peripheral populations
The mean significant SGS in the peripheral populations (up to 90 m) was up to six times of that was observed in the core (15 m) populations. The difference in SGS between the core and peripheral populations was statistically significant (ω = 65.26, P = 0.001). Both Nova Scotia populations (CL and RL) showed the highest within-population SGS, and the lowest genetic diversity, stand density and N e . These populations are highly fragmented and affected by human disturbance for over 400 years ever since the arrival of European settlers in the Province. The several-folds higher SGS levels observed in the two marginal and highly fragmented populations from Nova Scotia is likely due to reduced stand density, gene dispersal, genetic diversity and N e observed in these populations (Tables 1 and 2) as well as increased founder effect and genetic drift due to anthropogenic fragmentation and other disturbances. Thus, the Table 4 Spatial autocorrelation coefficients (r), single-class squared pair-sample t-test (t 2 ) and multi-class and multipopulation (ω) heterogeneity tests of spatial autocorrelation coefficients (r) for individual populations, population pairs and subsets of populations     [20]; Nb Wright's neighbourhood size; σg, estimate of gene flow distance from Sp, assuming different ratios of effective (De) to census density (D); *, significant (P < 0.05) ED, Edmundston; PA, Perth Andover; CL, Cedar Lake; RL, Rockland magnitude of evolutionary forces operating in core versus peripheral populations is probably different. It was not possible to dissect out the effects of anthropogenic fragmentation and disturbances on SGS in the peripheral populations from that of their peripheral location.
Comparisons of SGS between pristine core and pristine peripheral populations would have been ideal. However, in the eastern Maritime Canadian range of eastern white cedar where our study was focused, no peripheral population could be found that is not anthropogenically fragmented and disturbed. Thus, fragmentation and disturbance may also have contributed to higher SGS observed in the peripheral populations. The negative correlations of SGS with stand density, genetic diversity parameters and N e are consistent with the population genetic theories and expectations. The negative correlation between genetic diversity and SGS could be due to the combined effect of founder effects and reduced stand density. The spatial distribution of genotypes was found to be random in the core populations of Sitka spruce, whereas significant positive SGS up to 500 m was observed in the peripheral populations of this species [8]. In eastern white cedar, we observed mean significant positive SGS up to 90 m (combined populations) in the peripheral and up to 15 m in the core populations. Thus, the magnitude of impact of ecological and evolutionary factors on SGS in peripheral populations of Sitka spruce appears to be much higher than in eastern white cedar. Sitka spruce is a predominantly outcrossing species. The extent of differences observed in within-population SGS in core versus peripheral populations between these species may be related to the differences in their mating system. The mixed-mating species with relatively low inherent genetic diversity, such as eastern white cedar, can perhaps tolerate better negative effects (increased inbreeding, decreased genetic diversity and N e ) than the predominantly outcrossing species with relatively high genetic diversity, such as Sitka spruce. However, this aspect remains to be further tested in more details globally.

Sp statistics and gene dispersal distance
In agreement with the SGS results estimated using Smouse and Peakall's [12] autocorrelation coefficient method, average value for Sp statistic was also higher in peripheral than in core populations of eastern white cedar. Vekemans and Hardy [20] reported significantly higher levels of Sp value in populations with low density as compared to the populations with high density. Since the density of trees in the peripheral populations was much lower than in the core ones, the effect of density on the results of Sp cannot be ruled out. Except in the population PA, Sp values observed in the other three populations of eastern white cedar were higher than the average value observed in five wind-dispersed temperate trees (Sp = 0.012) and 18 outcrossing plant species [20]. However, it was lower in eastern white cedar as compared to the average value (Sp = 0.037) of seven mixed mating species [20]. De-Lucas et al. [7] reported higher Sp values in fragmented than in continuous populations of another conifer species, Pinus pinaster, and they explained that the lower effective population size due to few reproductive individuals in fragmented populations may have resulted in restricted gene flow via seed and increased the level of SGS in the populations. This explanation may also be applicable to eastern white cedar because the peripheral populations were highly affected by anthropogenic fragmentation. The average gene flow distance in eastern white cedar was higher in peripheral populations than in the core ones. Higher gene flow distances would be expected in fragmented and peripheral populations than in core and intact populations because of lower tree density in peripheral populations resulting in lower obstruction for pollen movement [20,69], which is the case in peripheral populations of eastern white cedar. In agreement with our results, Vekemans and Hardy [20] also reported higher gene flow distances in low density populations than in populations with high density.

Evolutionary potential and genetic resource conservation implications
The results of genetic diversity, Ne and SGS comparisons in core versus peripheral populations suggest that the evolutionary potential of the studied peripheral eastern white cedar populations may have been adversely affected. The results of our study have significance and implications for conservation and management of genetic resources of eastern white cedar -a declining to endangered species in Atlantic Canada. For ex-situ conservation via seed and vegetative propagule collection, a minimum distance of 75 meters in peripheral and 20 m in core populations is advisable to capture maximum genetic diversity and avoid sampling genetically-related trees. The core populations in New Brunswick should be conserved in-situ. Both of the peripheral Nova Scotia populations have suffered genetic degradation as a result of extensive human disturbance and edge effect. Here, both ex-situ and in-situ genetic resource conservation efforts supplemented with genetic restoration will be required given that eastern white cedar is a red-listed species in Nova Scotia. The results of our study are relevant for within-population sampling of individuals for population genetic studies in this species. The results may also be applicable for conservation and management of genetic resources of other long-lived, declining species with similar mixed-mating characteristics.

Conclusions
In conclusion, our results clearly demonstrate that the studied eastern white cedar populations have significant fine-scale genetic structure at short distances. Peripheral populations have several-folds higher within-population fine-scale genetic structure than core populations. Vegetative propagation does not apparently contribute to SGS in eastern white cedar. Anthropogenic disturbances and fragmentation presumably have significant effects on fine-scale genetic structure in eastern white cedar. The magnitude of evolutionary forces such as founder effect, genetic drift, gene flow and selection, operating in core versus peripheral populations is likely to be different. This may be the first report comparing SGS in intact core versus fragmented peripheral populations of a long-lived, mixed-mating declining plant. The results of our study contribute to the knowledge of poorlyunderstood spatial genetic structure of core versus peripheral populations in plants.