Higher fine-scale genetic structure in peripheral than in core populations of a long-lived and mixed-mating conifer - eastern white cedar (Thuja occidentalisL.)
© Pandey and Rajora; licensee BioMed Central Ltd. 2012
Received: 18 October 2011
Accepted: 5 April 2012
Published: 5 April 2012
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.
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.
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.
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 . Genetic studies involving both core and peripheral populations are still limited . Additionally, the available studies mainly focused on within-and among-population genetic diversity [2–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–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  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–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 micro-environmental 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 . 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 . In contrast, generally, outcrossing decreases population subdivision by reducing mating among relatives and increasing N e . 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 . 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 . SGS has been found to be inversely correlated with adult density within populations . 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–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 long-lived 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 . 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 .
Like other cedar species (Thuja plicata, T. orientalis, [36–38]), eastern white cedar also possesses relatively low genetic diversity [39–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–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 , 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 within-population 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–13].
Populations and sampling
Geographical location, sample size and stand density of four sampled populations of eastern white cedar
No. of individuals sampled
Stand density (Trees/ha)
47° 33' 16″
68° 13' 53″
46° 44' 53″
67° 42' 34″
44° 50' 22″
65° 01' 48″
45° 00' 35″
64° 42' 10″
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 . 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 MgCl2, 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 Ne
The standard population genetic diversity parameters (number of alleles, observed and expected heterozygosity) and fixation index (F) were calculated using the program GENALEX . Allelic richness (AR) was estimated using the program FSTAT . Genotypic diversity was determined by calculating observed and expected genotype additivity (GAO and GAE) as the sum of genotypes observed and expected under Hardy-Weinberg equilibrium over the loci . 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 .
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 . 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).
Pattern of spatial distribution of individual trees in the four eastern white cedar populations were estimated using univariate second-order Ripley's K-function for 10, 15, 20 and 30 m distance classes (t). The weighted edge correction  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 .
Spatial Genetic Structure (SGS)
We used Smouse and Peakall's  spatial autocorrelation method to estimate the SGS in the four eastern white cedar populations employing the GENALEX 6.2  program. This method determines spatial distribution of genetic variation based on a frequency-weighted 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  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.  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.  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. . 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.  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 (t2) . SGS heterogeneity over all distance classes (a) in a population and (b) among populations was determined by using the ω test criteria  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.
Spstatistics 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 = -^bF/(1-^F 1 ), where -^bF is the regression slope of the autocorrelogram and ^ F 1 is the mean kinship coefficient (Fij)  between individuals belonging to the first (0-15 m) distance class . Using the Sp statistics, Wright's neighborhood size Nb (Nb = 1/Sp) was also estimated. Nb is defined as Nb = 4πD e σg2, where σg2 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 -^bF and ^ F 1 were carried out using SPAGEDI program . Statistical significance (P ≤ 0.05) of the analyses was performed by permutation tests with 10,000 replications.
Genetic diversity and effective population size
Genetic diversity measures, fixation index and their (SE), and effective population size for four sampled populations of eastern white cedar
Number of pairs of observation per distance class for the four eastern white cedar populations
Distance class (m)
Spatial genetic structure (SGS)
Spatial autocorrelation coefficients (r), single-class squared pair-sample t-test (t2) and multi-class and multi-population (ω) heterogeneity tests of spatial autocorrelation coefficients (r) for individual populations, population pairs and subsets of populations
rand t 2 values and their statistical significance in different distance classes
Overall multi-distance class or multi population correlogram significance (ω-test criterion)
Spatial autocorrelation coefficients (r)
Four pops combined
Single-class squared pair-sample t-test (t2)
ED vs. PA
ED vs. CL
ED vs. RL
PA vs. CL
PA vs. RL
CL vs. RL
Core vs. peripheral
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 over the 0-15 m and 61-90 m distances (up to 90 m) (Table 4).
Spstatistics and dispersal estimates
Spatial genetic structure, Sp statistics and estimates of gene flow distances in four eastern white cedar populations
σg (m) (De/D= 0.1)
σg (m) (De/D= 0.5)
Correlations of SGS parameters with genetic diversity, stand density, and population Ne
There was negative correlation of the significant SGS (spatial autocorrelation coefficient r) with the population genetic diversity measures. The correlations (Pearson correlation coefficient R) of r with genetic diversity measures were as follows: mean number of alleles per locus (A) = -0.85 P = 0.034), observed heterozygosity (HO) = - 0.89 (P = 0.022), expected heterozygosity (HE) = -0.82 (P = 0.044), observed genotypic richness (GAO) = -0.88 (P = 0.024) and expected genotypic richness (GAE) = 08 P = 0.040). The correlation of eastern white cedar stand density with r a -0.69 P = 0.096). The significant SGS r values showed negative high correlation (Pearson R = -0.65)wt Ne, although the correlation was not statistically significant (P = 0.139).
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% ) species Larix larcinia (up to 20 m) . 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 . 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  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  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.  and Trabaud et al.  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 aggregated seedlings after stand disturbance due to fire. These factors may have contributed to the observed spatial pattern of individual trees in eastern white cedar populations. The effect of limited seed dispersal on the significant aggregation of individual trees may be evident from the significant SGS (15-75 m) observed in the four populations. Although we are not aware of any records of fire in the studied populations, past occurrence of fire on the studied populations cannot be ruled out because fire is one of the main agents of natural disturbance in boreal forest .
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 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 . 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.
Spstatistics and gene dispersal distance
In agreement with the SGS results estimated using Smouse and Peakall's  autocorrelation coefficient method, average value for Sp statistic was also higher in peripheral than in core populations of eastern white cedar. Vekemans and Hardy  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 . However, it was lower in eastern white cedar as compared to the average value (Sp = 0.037) of seven mixed mating species . De-Lucas et al.  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  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.
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 poorly-understood spatial genetic structure of core versus peripheral populations in plants.
Total number of alleles
Mean number of alleles per locus
- A E :
Effective number of alleles
Number of private alleles
- b-log (^bF):
Slope of the regression of kinship with ln(dist)
- De :
- F :
- Fij :
Mean kinship coefficient
Geographical positioning system
Observed genotype additivity (richness)
Expected genotype additivity
- Nb :
Wright's neighbourhood size
- Ne :
Effective population size
- r :
Spatial autocorrelation coefficients
- Sp :
Estimate of gene flow distance from Sp
Spatial genetic structure
We thank Tony Duke of the Nova Scotia Department of Natural Resources for drawing our attention to eastern white cedar conservation genetic issues. We also thank Lawrence Benjamin of Nova Scotia Department of Natural Resources, Michael H. McDonald of New Brunswick Department of Natural Resources and Brenda Tse for their help with the field work and Brenda Tse and Sandra Tucker for their laboratory assistance. Our sincere thanks go to Dr. Lisa O'Connell, who conducted field sampling, DNA extraction and genotyping work for two of the four populations sampled in this study. The research was funded by the research grants provided by the Nova Scotia Department of Natural Resources, PEI Model Forest Network Partnerships Ltd., Canada Research Chair Program (CRC950-201869) funds and the Natural Sciences and Engineering Research Council of Canada Discovery Grant (RGPIN 170651) to O.P. Rajora.
- Aitken SN, Yeaman S, Holliday JA, Wang T, Curtis-McLane S: Adaptation, migration or extirpation: climate change outcomes for tree populations. Evol Appl. 2008, 1: 95-111. 10.1111/j.1752-4571.2007.00013.x.PubMedPubMed CentralView ArticleGoogle Scholar
- Eckert CG, Samis KE, Lougheed SC: Genetic variation across species' geographical ranges: the central-marginal hypothesis and beyond. Mol Ecol. 2008, 17: 1170-1188. 10.1111/j.1365-294X.2007.03659.x.PubMedView ArticleGoogle Scholar
- Lesica P, Allendorf FW: When are peripheral populations valuable for conservation?. Conserv Biol. 1995, 9: 753-760. 10.1046/j.1523-1739.1995.09040753.x.View ArticleGoogle Scholar
- Hampe A, Petit RJ: Conserving biodiversity under climate change: the rear edge matters. Ecol Letters. 2005, 8: 461-467. 10.1111/j.1461-0248.2005.00739.x.View ArticleGoogle Scholar
- Williams DA, Wang Y, Borchetta M, Gaines MS: Genetic diversity and spatial structure of a keystone species in fragmented pine rockland habitat. Biol Conserv. 2007, 138: 256-268. 10.1016/j.biocon.2007.04.024.View ArticleGoogle Scholar
- Born C, Hardy OJ, Chevalier MH, Ossari S, Attéke EJW, Hossaert-Mckey M: Small-scale spatial genetic structure in the central African rainforest tree species Aucoumea klaineana: a stepwise approach to infer the impact of limited gene dispersal, population history and habitat fragmentation. Mol Ecol. 2008, 17: 2041-2050. 10.1111/j.1365-294X.2007.03685.x.PubMedView ArticleGoogle Scholar
- De-Lucas AI, Gonza'lez-Marti'nez SC, Vendramin GG, Hidalgo E, Heuertz M: Spatial genetic structure in continuous and fragmented populations of Pinus pinaster Aiton. Mol Ecol. 2009, 18: 4564-4576. 10.1111/j.1365-294X.2009.04372.x.PubMedView ArticleGoogle Scholar
- Gapare WJ, Aitken SN: Strong spatial genetic structure in peripheral but not core populations of Sitka spruce [Picea sitchensis (Bong.) Carr.]. Mol Ecol. 2005, 14: 2659-2667. 10.1111/j.1365-294X.2005.02633.x.PubMedView ArticleGoogle Scholar
- Jain SK: Population structure and the effects of breeding system. Plant Genetic Resources: Today and Tomorrow. Edited by: Frankel O, Hawkes JG. 1975, Cambridge Univ. Press, Cambridge, UK, 15-36.Google Scholar
- Loveless MD, Hamrick JL: Ecological determinants of genetic structure in plant populations. Ann Rev Ecol Syst. 1984, 15: 65-95. 10.1146/annurev.es.15.110184.000433.View ArticleGoogle Scholar
- Hamrick JL, Nason JD: Consequence of dispersal in plants. Population Dynamics in Ecological Space and Time. Edited by: Rhodes OE, Ronald KC, Smith MH. 1996, The University of Chicago Press, Chicago, USA, 203-235.Google Scholar
- Smouse PE, Peakall R: Spatial autocorrelation analysis of individual multiallele and multilocus genetic structure. Heredity. 1999, 82: 561-573. 10.1038/sj.hdy.6885180.PubMedView ArticleGoogle Scholar
- Kalisz S, Nason JD, Hanzawa FM, Tonsor SJ: Spatial population genetic structure in Trillium grandiflorum: the roles of dispersal, mating, history and selection. Evol Biol. 2001, 55: 1560-1568.Google Scholar
- Young AG, Merriam HG: Effects of forest fragmentation on the spatial genetic structure of Acer saccharum Marsh (sugar maple) populations. Heredity. 1994, 72: 201-208. 10.1038/hdy.1994.27.View ArticleGoogle Scholar
- Jump AS, Peñuelas J: Genetic effects of chronic habitat fragmentation in a wind-pollinated tree. Proc Natl Acad Sci. 2006, 103: 8096-8100. 10.1073/pnas.0510127103.PubMedPubMed CentralView ArticleGoogle Scholar
- Epperson BK: Plant dispersal, neighbourhood size and isolation by distance. Mol Ecol. 2007, 16: 3854-3865. 10.1111/j.1365-294X.2007.03434.x.PubMedView ArticleGoogle Scholar
- Troupin D, Nathan R, Vendramin GG: Analysis of spatial genetic structure in an expanding Pinus halepensis population reveals development of fine-scale genetic clustering over time. Mol Ecol. 2006, 15: 3617-3630. 10.1111/j.1365-294X.2006.03047.x.PubMedView ArticleGoogle Scholar
- Wright S: Isolation by distance. Genetics. 1943, 28: 114-138.PubMedPubMed CentralGoogle Scholar
- Reusch TBH, Hukriede W, Stam WT, Olsen JL: Differentiating between clonal growth and limited gene flow using spatial autocorrelation of microsatellites. Heredity. 1999, 83: 120-126. 10.1046/j.1365-2540.1999.00546.x.PubMedView ArticleGoogle Scholar
- Vekemans X, Hardy OJ: New insights from fine-scale spatial genetic structure analyses in plant populations. Mol Ecol. 2004, 13: 921-934. 10.1046/j.1365-294X.2004.02076.x.PubMedView ArticleGoogle Scholar
- Boyle TJB, Liengsiri C, Piewluang C: Genetic structure of black spruce on two contrasting sites. Heredity. 1990, 65: 393-399. 10.1038/hdy.1990.109.View ArticleGoogle Scholar
- Epperson BK, Allard RW: Spatial autocorrelation analysis of the distribution of genotypes within populations of lodgepole pine. Genetics. 1989, 121: 369-377.PubMedPubMed CentralGoogle Scholar
- Knowles P: Spatial genetic structure within two natural stands of black spruce (Picea mariana (Mill.) B. S. P.). Silvae Genet. 1991, 40: 13-19.Google Scholar
- Knowles P, Perry DJ, Foster HA: Spatial genetic structure in two tamarack (Larix laricina (Du Roi) K. Koch) populations with differing establishment histories. Evolution. 1992, 46: 572-576. 10.2307/2409875.View ArticleGoogle Scholar
- Xie CY, Knowles P: Spatial genetic substructure within natural populations of jack pine (Pinus banksiana). Can J Bot. 1991, 69: 547-551. 10.1139/b91-074.View ArticleGoogle Scholar
- Leonardi S, Raddi S, Borghetti M: Spatial autocorrelation of allozyme traits in a Norway spruce (Picea abies) population. Can J For Res. 1996, 26: 63-71. 10.1139/x26-007.View ArticleGoogle Scholar
- Epperson BK, Chung MG: Spatial genetic structure of allozyme polymorphisms within populations of Pinus strobus (Pinaceae). Am J Bot. 2001, 88: 1006-1010. 10.2307/2657081.PubMedView ArticleGoogle Scholar
- Parker KC, Hamrick JL, Parker AJ, Nason JD: Fine-scale genetic structure in Pinus clausa (Pinaceae) populations: effects of disturbance history. Heredity. 2001, 87: 99-113. 10.1046/j.1365-2540.2001.00914.x.PubMedView ArticleGoogle Scholar
- Bucci G, Menozzi P: Spatial autocorrelation and linkage of Mendelian RAPD markers in a population of Picea abies Karst. Mol Ecol. 2002, 11: 305-315. 10.1046/j.1365-294X.2002.01444.x.PubMedView ArticleGoogle Scholar
- González-Martínez SC, Gerber S, Cervera MT, Martínez-Zapater JM, Gil L, Alía R: Seed gene flow and fine-scale structure in a Mediterranean pine (Pinus pinaster Ait.) using nuclear microsatellite markers. Theor Appl Genet. 2002, 104: 1290-1297. 10.1007/s00122-002-0894-4.PubMedView ArticleGoogle Scholar
- Marquardt PE, Epperson BK: Spatial and population genetic structure of microsatellites in white pine. Mol Ecol. 2004, 13: 3305-3315. 10.1111/j.1365-294X.2004.02341.x.PubMedView ArticleGoogle Scholar
- Takahashi T, Tani N, Niiyama K, Yoshida S, Taira H, Tsumura Y: Genetic succession and spatial genetic structure in a natural old growth Cryptomeria japonica forest revealed by nuclear and chloroplast microsatellite markers. For Ecol Manage. 2008, 255: 2820-2828. 10.1016/j.foreco.2008.01.055.View ArticleGoogle Scholar
- Newell RE: Provincial (Nova Scotia) Status report on northern white cedar Thuja occidentalis prepared for the Nova Scotia Department of Natural Resources. 2005Google Scholar
- Jhonston WF: Thuja occidentalis L.: northern white cedar. Silvics of North America. 1. Conifers. Edited by: Burns RM, Honkala BH. 1990, US Department of Agriculture Forest Service Handbook 654, Washington DC, USA, 580-589.Google Scholar
- Bentley PA, Smith EC: A study of northern white cedar and jack pine in Nova Scotia. Proc NS Inst Sci. 1959, 24: 376-398.Google Scholar
- Xie CY, Dancik BP, Yeh FC: The mating system in natural populations of Thuja orientalis. Can J For Res. 1991, 21: 333-339. 10.1139/x91-041.View ArticleGoogle Scholar
- Xie CY, Dancik BP, Yeh FC: Genetic structure of Thuja orientalis. Biochem Syst Ecol. 1992, 20: 433-441. 10.1016/0305-1978(92)90083-P.View ArticleGoogle Scholar
- El-Kassaby YA, Russell J, Ritland K: Mixed mating in an experimental populations of western red cedar, Thuja plicata. J Hered. 1994, 85: 227-231.Google Scholar
- Perry DJ, Knowles P, Yeh FC: Allozyme variation of Thuja occidentalis L. in Northwestern Ontario. Biochem Syst Ecol. 1990, 18: 111-115. 10.1016/0305-1978(90)90045-H.View ArticleGoogle Scholar
- Matthes-Sears U, Mattes-Sears SC, Larson DW: Sources of allozymic variation in Thuja occidentalis in southern Ontario, Canada. Silvae Genet. 1991, 40: 100-105.Google Scholar
- Lamy S, Bouchard A, Simon J-P: Genetic structure, variability and mating system in eastern white cedar (Thuja occidentalis) populations of recent origin in an agricultural landscape in southern Québec. Can J For Res. 1999, 29: 1383-1392.View ArticleGoogle Scholar
- Perry DJ, Knowles P: Evidence of high self-fertilization in natural populations of eastern white cedar (Thuja occidentalis). Can J Bot. 1990, 68: 663-668. 10.1139/b90-086.View ArticleGoogle Scholar
- O'Connell LM, Ritland K, Thompson SL: Patterns of post-glacial colonization by western redcedar (Thuja plicata, Cupressaceae) as revealed by microsatellite markers. Botany. 2008, 86: 194-203. 10.1139/B07-124.View ArticleGoogle Scholar
- O'Connell LM, Ritland CE: Characterization of microsatellite loci in western red cedar (Thuja plicata). Mol Ecol. 2001, 9: 1920-1922.View ArticleGoogle Scholar
- White E, Hunter J, Dubetz C, Brost R, Bratton A, Edes S, Sahota R: Microsatellite markers for individual tree genotyping: application in forest crime prosecutions. J Chem Tech Biotech. 2000, 70: 923-926.View ArticleGoogle Scholar
- Nakao Y, Iwata H, Matsumoto A, Tsumura Y, Tomaru N: Highly polymorphic microsatellite markers in Chamaecyparis obtusa. Can J For Res. 2001, 31: 2248-2251. 10.1139/x01-145.View ArticleGoogle Scholar
- Peakall R, Smouse PE: GENALEX 6: genetic analysis in Excel. Population genetic software for teaching and research. Mol Ecol Notes. 2006, 6: 288-295. 10.1111/j.1471-8286.2005.01155.x.View ArticleGoogle Scholar
- Gaudet J: FSTAT version 1.2: a computer program to calculate F-statistics. J Hered. 1995, 86: 485-486.Google Scholar
- Rajora OP, Rahman MH, Buchert GP, Danick BP: Microsatellite DNA analysis of genetic effects of harvesting in old-growth eastern white pine (Pinus storbus) in Ontario. Mol Ecol. 2000, 9: 339-348. 10.1046/j.1365-294x.2000.00886.x.PubMedView ArticleGoogle Scholar
- Raymond M, Rousset F: GENEPOP version 1.2: population genetics software for the exact tests and ecumenicism. J Hered. 1995, 86: 248-249.Google Scholar
- Beerli P, Felsenstein J: Maximum likelihood estimation of migration rates and population numbers of two populations using a coalescent approach. Genetics. 1999, 152: 763-773.PubMedPubMed CentralGoogle Scholar
- Ripley BD: The second-order analysis of stationary point processes. J Appl Prob. 1976, 13: 255-266. 10.2307/3212829.View ArticleGoogle Scholar
- Ripley BD: Modelling spatial patterns. J R Stat Soc. 1977, B39: 172-212.Google Scholar
- Perry GLW: SpPack: spatial point pattern analysis in Excel using Visual Basic for Applications (VBA). Environ Model Soft. 2004, 19: 559-569. 10.1016/j.envsoft.2003.07.004.View ArticleGoogle Scholar
- Oddou-Muratorio S, Demesure-Musch B, Pelissier R, Gouyon PH: Impacts of gene flow and logging history on the local genetic structure of a scattered tree species, Sorbus torminalis L. Crantz. Mol Ecol. 2004, 13: 3689-3702. 10.1111/j.1365-294X.2004.02373.x.PubMedView ArticleGoogle Scholar
- Smouse PE, Peakall R, Gonzales E: A heterogeneity test for fine-scale genetic structure. Mol Ecol. 2008, 17: 3389-3400. 10.1111/j.1365-294X.2008.03839.x.PubMedView ArticleGoogle Scholar
- Loiselle BA, Sork VL, Nason J, Graham C: Spatial genetic-structure of a tropical understory shrub, Psychotria officinalis (Rubiaceae). Am J Bot. 1995, 82: 1420-1459. 10.2307/2445869.View ArticleGoogle Scholar
- Hardy OJ, Vekemans X: SPAGEDi: a versatile computer program to analyse spatial genetic structure at the individual or population levels. Mol Ecol Notes. 2002, 2: 618-620. 10.1046/j.1471-8286.2002.00305.x.View ArticleGoogle Scholar
- Pandey M, Rajora OP: Genetic diversity and differentiation of core versus peripheral populations of eastern white cedar, Thuja occidentalis L. (Cupressaceae). Am J Botany. 2012, 99: 690-699. 10.3732/ajb.1100116.View ArticleGoogle Scholar
- Knowles P, Furnier GR, Aleksiuk MA, Perry DJ: Significant levels of self-fertilization in natural populations of tamarack. Can J Bot. 1987, 65: 1087-1091. 10.1139/b87-151.View ArticleGoogle Scholar
- Levin DA, Kerster HW: Gene flow in seed plants. In Evolutionary Biology 7. Edited by: Dobzhansky T, Hecht MT, Steere WC. 1974, New York: Plenum Press, 139-220.View ArticleGoogle Scholar
- Schuster WSF, Mitton JB: Paternity and gene dispersal in limber pine (Pinus flexilis James). Heredity. 2000, 84: 348-361. 10.1046/j.1365-2540.2000.00684.x.PubMedView ArticleGoogle Scholar
- Fowells HA: Silvics of forest trees of the United States. U. S. Department of Agriculture. 1965, Agriculture Handbook 271Google Scholar
- De Blois S, Bouchard A: Dynamics of Thuja occidentalis in an agricultural landscape of southern Quebec. J Veg Sci. 1995, 6: 531-542. 10.2307/3236352.View ArticleGoogle Scholar
- Chung MG, Epperson BK: Clonal and spatial genetic structure in Eurya emarginata (Theaceae). Heredity. 2000, 84: 170-177. 10.1046/j.1365-2540.2000.00644.x.PubMedView ArticleGoogle Scholar
- Parker KC, Parker AJ, Beaty M, Fuller MM, Faust TD: Population structure and spatial pattern of two coastal populations of Ocala sand pine (Pinus clausa (Chapm. ex Engelm.) Vasey ex Sarg. var. clausa D. B. Ward). J Torr Bot Soc. 1997, 124: 22-33. 10.2307/2996595.View ArticleGoogle Scholar
- Trabaud L, Michels C, Grosman J: Recovery of burnt Pinus halepensis Mill. forests. II. Pine reconstitution after wildfire. For Ecol Manage. 1985, 113: 67-179.Google Scholar
- McCullough DG, Werner RA, Neumann D: Fire and insects in northern and boreal forest ecosystems of North America. Ann Rev Entomol. 1998, 43: 107-127. 10.1146/annurev.ento.43.1.107.View ArticleGoogle Scholar
- O'Connell L, Mosseler A, Rajora OP: Extensive long-distance pollen dispersal in a fragmented landscape maintains genetic diversity in white spruce. J Heredity. 2007, 98: 640-645. 10.1093/jhered/esm089.View ArticleGoogle 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.