Characterization of microsatellite loci in western redcedar (Thuja plicata)

Genetic consequences of fragmentation in "arbor vitae," eastern white cedar (Thuja occidentalis L.), toward the northern limit of its distribution range

We tested the hypothesis that marginal fragmented populations of eastern white cedar (EWC) are genetically isolated due to reduced pollen and gene flow. In accordance with the central-marginal model, we predicted a decrease in population genetic diversity and an increase in differentiation along the latitudinal gradient from the boreal mixed-wood to northern coniferous forest. A total of 24 eastern white cedar populations were sampled along the north-south latitudinal gradient for microsatellite genotyping analysis. Positive F_is values and heterozygote deficiency were observed in populations from the marginal (F_is = 0.244; P_HW = 0.0042) and discontinuous zones (F_is = 0.166; P_HW = 0.0042). However, populations from the continuous zone were in HW equilibrium (F_is = −0.007; P_HW = 0.3625). There were no significant latitudinal effects on gene diversity (H_s), allelic richness (AR), or population differentiation (F_st). Bayesian and NJT (neighbor-joining tree) analyses demonstrated the presence of a population structure that was partly consistent with the geographic origins of the populations. The impact of population fragmentation on the genetic structure of EWC is to create a positive inbreeding coefficient, which was two to three times higher on average than that of a population from the continuous zone. This result indicated a higher occurrence of selfing within fragmented EWC populations coupled with a higher degree of gene exchange among near-neighbor relatives, thereby leading to significant inbreeding. Increased population isolation was apparently not correlated with a detectable effect on genetic diversity. Overall, the fragmented populations of EWC appear well-buffered against effects of inbreeding on genetic erosion.

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.

Genetic diversity and differentiation of core vs. peripheral populations of eastern white cedar, Thuja occidentalis (Cupressaceae)

PREMISE OF THE STUDY

Geographically peripheral (marginal) populations are expected to have lower genetic diversity and higher genetic differentiation than geographically core (central) populations as a result of supposedly lower effective population size (N(e)) and higher genetic drift, founder effect, fragmentation, and isolation in peripheral than in core populations. Here we address this issue for a long-lived plant species, eastern white cedar (Thuja occidentalis).

METHODS

Genetic diversity and population structure of 13 natural populations of eastern white cedar from its Canadian eastern peripheral and core natural ranges in New Brunswick, Nova Scotia, and Prince Edward Island were studied using six nuclear microsatellite DNA markers.

KEY RESULTS

The core populations of eastern white cedar had significantly higher allelic diversity (mean A = 8.83, A(r) = 8.13, A(e) = 4.03) and N(e) (428) than the peripheral populations (A = 6.64, A(r) = 6.15, A(e) = 3.12, N(e) = 198). However, expected heterozygosity was similar in the core (H(e) = 0.64) and peripheral (H(e) = 0.60) populations. Genetic differentiation was significantly higher among the peripheral (F(ST) = 0.089) than among the core (F(ST) = 0.032) populations. No genetic differentiation (F(ST)/Φ(RT) = 0.000) was detected between core and peripheral regions.

CONCLUSIONS

Peripheral populations have significantly lower N(e) and genetic diversity in terms of allelic diversity (richness) and significantly higher genetic differentiation than the core populations of eastern white cedar in its Canadian eastern range. However, core and peripheral populations have similar levels of expected heterozygosity. Implications for conservation of eastern white cedar genetic resources are discussed.

Contrasting microsatellite variation between subalpine and western larch, two closely related species with different distribution patterns

Subalpine larch (Larix lyallii Parl.) and western larch (Larix occidentalis Nutt.) represent two closely related species with contrasting abundance and distribution patterns in Western North America. Genetic diversity at seven informative microsatellite loci was determined for 19 populations of subalpine larch and nine populations of western larch. Contrasting genetic diversity and patterns of population differentiation were observed between the two species. The overall within-population genetic diversity parameters were lower in subalpine larch (A = 3.2; A(P) = 3.6; H(E) = 0.418) than in western larch (A(P) = 5.51; H(E) = 0.580), a pattern that is likely related to historical or demographic factors. No evidence of interspecific hybridization was observed. Significantly more population differentiation (theta = 0.15; R(ST) = 0.07), consistent with more restricted gene flow, was observed for subalpine larch as compared to western larch (theta = 0.05; R(ST) = 0.04). Under the assumption of an infinite allele mutation model, 12 of the 19 subalpine larch populations showed signs of deviation from the mutation-drift equilibrium, which suggests Holocene population bottlenecks and fluctuations in effective population size for this species. None of the western larch populations deviated significantly from the mutation-drift equilibrium. For both species, Mantel's test revealed a significant positive relationship between geographical and genetic distances indicative of isolation by distance. A similar geographical structure was detected in both species, suggesting at least two genetically distinct glacial populations in each species. The various implications for gene conservation are discussed.

Fine-scale estimation of outcrossing in western redcedar with microsatellite assay of bulked DNA

Western redcedar (Thuja plicata, Cupressaceae) is a self-fertile conifer with a mixed mating system and significant variation for outcrossing among populations. In this paper, we conducted a fine-scale study of mating system variation to identify correlates of outcrossing in natural populations. We examined variation for outcrossing within and among individual trees, and describe a new method to estimate outcrossing using bulked DNA samples. Bulking (assaying DNA tissues from several individuals simultaneously) increases the experimental power without increasing the experimental effort. We sampled 80 trees from four natural populations in southwestern British Columbia. From each tree, we sampled from up to six crown positions (three heights and inner vs outer branches). From each position, two samples of three seedlings each were bulked before DNA extractions. Using four microsatellite loci, we obtained outcrossing rates for each tree and for each of the six crown positions. We found individual tree selfing rates to increase with tree height in all four populations, but selfing rates did not differ among crown positions. The higher selfing rate of larger trees is probably due to their greater proportional contribution to local pollen clouds. Individual tree outcrossing rates ranged from 22 to 100% and the population outcrossing rates from 66 to 78%. Missed alleles due to bulking and the estimation method used both cause a downward bias in outcrossing rates, so that these estimates are probably lower than the actual outcrossing rates. Nevertheless, the trends we observed are not affected by systematic biases of estimation.