Am. J. Bot.
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (19)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Clark, C. M.
Right arrow Articles by O'Malley, D. M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Clark, C. M.
Right arrow Articles by O'Malley, D. M.
Agricola
Right arrow Articles by Clark, C. M.
Right arrow Articles by O'Malley, D. M.
(American Journal of Botany. 2000;87:774-782.)
© 2000 Botanical Society of America, Inc.

Genetic discontinuity revealed by chloroplast microsatellites in eastern North American Abies (Pinaceae)1

Catherine M. Clark4,2, Thomas R. Wentworth3 and David M. O'Malley2

2 Department of Forestry, Forest Biotechnology Group, Box 7247, North Carolina State University, Raleigh, North Carolina 27695 USA; 3 Department of Botany, Box 7612, North Carolina State University, Raleigh, North Carolina 27695 USA

Received for publication March 30, 1999. Accepted for publication September 10, 1999.

ABSTRACT

Development of conservation strategies for Fraser fir (Abies fraseri) in the southern Appalachian Mountains depends in part on recognition of the extent to which Fraser fir is genetically distinct from the closely related balsam (A. balsamea) and intermediate (A. balsamea var. phanerolepis) fir. These sibling species have exhibited intergrading, clinal variation in morphological, chemical, and genetic characteristics in prior research. Chloroplast microsatellite markers were polymerase chain reaction amplified from genomic DNA samples of 78 individuals representing the geographic ranges of Fraser, balsam, and intermediate fir. Gene diversity levels at two loci ranged among taxa from 0.65 to 0.84. Allele frequencies demonstrated significant differentiation among taxa, with RST values of 0.36 and 0.10. Haplotype diversity and DSH were highest for balsam fir and lowest for intermediate fir. A haplotype network analysis based on allele size distribution for the two loci revealed two distinct clusters of haplotypes and population-specific haplotypes. Ninety-two percent of the haplotypes in one cluster were from balsam fir and intermediate fir, and 84% of the haplotypes in the other cluster were from Fraser fir and intermediate fir. The genetic differentiation of chloroplast DNA markers provides justification for the recognition of Fraser fir as a distinct Management Unit (MU) for conservation purposes, regardless of its taxonomic classification.

Key Words: Abies • chloroplast DNA • chloroplast haplotypes • chloroplast microsatellites • chloroplast SSRs • conservation • Fraser fir • Pinaceae

Eastern North American Abies (Mill.), Fraser fir [Abies fraseri (Pursh) Poir.], balsam fir [Abies balsamea (L.) Mill.], and a variety of balsam, intermediate fir (Abies balsamea var. phanerolepis Fern.), exhibit minimal reproductive isolation (Hawley and DeHayes, 1985 ) and are considered closely related sibling species in the Section Balsameae (Flora of North America Editorial Committee, 1993 ). Classification of these three taxa is based largely on cone bract length to subtending scale ratio, number of leaf hypodermal cells, and geographic location (Fig. 1) (Myers and Bormann, 1963 ; Robinson and Thor, 1969 ; Liu, 1971 ; Thor and Barnett, 1974 ). In recent times, populations in the southern limits of growth are found only at higher elevations along the Appalachian Mountains, due to a lack of suitable habitat. The resulting gap of several hundred kilometres in the distribution of fir contributes to the recognition of two species, balsam and Fraser fir. Based on the low level of character divergence, some authorities consider eastern North American fir a single species with three varieties (Thor and Barnett, 1974 ; Jacobs, Werth, and Guttman, 1984 ). Taxonomic boundaries are obscured by continuous variation in foliage, seed, and cone morphology (Robinson and Thor, 1969 ), and clinal variation in terpene concentrations (Thor and Barnett, 1974 ). Prior evaluation of allozyme variation among the three taxa showed nearly identical complements of alleles at 13 polymorphic loci (Jacobs, Werth, and Guttman, 1984 ).



View larger version (29K):
[in this window]
[in a new window]
 
Fig. 1. Geographic ranges of Fraser (Abies fraseri), balsam (A. balsamea), and intermediate (A.b. var. phanerolepis) fir. Locations of individuals included in the range-wide and regional surveys are indicated by symbols ()

 
Efforts to focus attention on the conservation of species or populations depend in part on our ability to recognize genetic and evolutionary distinctiveness. Fraser fir has been recognized as a V2 (vulnerable) taxon (Oldfield, Lusty, and MacKinven, 1998 ). Declines of 44–91% of mature, reproductive Fraser fir in spatially isolated, high-elevation ecosystems are in part attributed to the presence of an introduced insect pest, the balsam woolly adelgid (Dull et al., 1988 ), although environmental factors are also implicated. In situ and ex situ gene conservation efforts will be facilitated by techniques that estimate diversity and genetic discontinuity for closely related groups.

Molecular genetic markers may better reveal genetic discontinuities and distinctiveness among or between taxa with slight morphological or reproductive differentiation. At the level of nuclear genes, interbreeding and genetic exchange obscure the evidence of past demographic patterns. Characteristics of allozyme diversity have been examined for many representatives of the Pinaceae, which are typically outcrossing, wind-pollinated species. Conifers often display minimal levels of differentiation among populations and high levels of genetic diversity at loci in the nuclear genome (Hamrick and Godt, 1990, 1996 ). Within the Pinaceae opposite uniparental inheritance of mitochondrial and chloroplast genomes provides opportunities to assess discontinuity within an effectively haploid genome. The alleles of loci in an organelle genome can be collectively considered as a haplotype since they are nonrecombining. In these organelles, effective population size is approximately one-half that of nuclear genes (McCauley, 1995 ) and, consequently, both the time to allele fixation within populations and response time to stochastic change are reduced. In plants, mitochondrial genomes have not typically been useful for phylogenetic analysis due to a high rate of sequence reorganization (Sederoff et al., 1981 ; Wu, Krutovskii, and Strauss, 1998 ). Mitochondrial haplotype diversity related to sequence rearrangement has, however, shown utility in population differentiation of pine and fir taxa (Strauss, Hong, and Hipkins, 1993 ; Tsumura and Suyama, 1998 ). In chloroplasts, the conservation and homology of sequence make it possible to compare genes across the plant kingdom (Clegg and Zurawski, 1992 ) and examine phylogenetic relationships in taxa that have diverged for hundreds of thousands to millions of years. Conversely, until hypervariable regions were located, the chloroplast genome had shown minimal utility in examining closely related taxa or intraspecific variation (McCauley, 1995 ).

The recent development of microsatellite, or SSR (Simple Sequence Repeat), markers in chloroplast DNA (cpDNA) with high levels of variation (Powell et al., 1995a, b ; Vendramin et al., 1996 ; Vendramin and Ziegenhagen, 1997 ) has provided greater resolution in differentiation among closely related taxa. Paternal inheritance of chloroplast genomes in most conifers (Neale, Wheeler, and Allard, 1986 ; Neale and Sederoff, 1989 ) makes chloroplast microsatellites particularly effective markers for studying mating systems, gene flow via both pollen and seeds, and uniparental lineage. Predominant paternal inheritance of chloroplast DNA has been demonstrated in European Abies with previously characterized (Vendramin et al., 1996 ) Pinus thunbergii primers at two highly variable microsatellite loci (Vendramin and Ziegenhagen, 1997 ; Ziegenhagen et al., 1998 ; Vendramin et al., 1999 ). Intraspecific diversity and population structure in heterologous amplification with these primers have been demonstrated in many pine species as well (Morgante and Olivieri, 1993 ; Powell et al., 1995a ; Cato and Richardson, 1996 ; Vendramin et al., 1996 ; Morgante, Felice, and Vendramin, 1998 ). Chloroplast microsatellite loci have also revealed population structure in species considered genetically depauperate (e.g., P. leucodermis, Powell et al., 1995a ; P. resinosa, Echt et al., 1998 ). In this study, we assessed genetic discontinuity and variability of chloroplast microsatellite loci in eastern North American Abies using fluorescent DNA fragment analysis. The geographic distribution of chloroplast haplotypes was evaluated, as well as criteria for recognizing distinct populations in eastern North American fir.

METHODS AND MATERIALS

DNA extraction and PCR amplification
Genomic DNA was extracted from fir needle tissue following the method of Doyle and Doyle (1990) . Three range-wide samples (26 from each taxon) of Fraser, balsam, and intermediate fir were analyzed (Table 1). Intermediate fir were classified primarily by their geographic locations described in earlier studies (Robinson and Thor, 1969 ; Thor and Barnett, 1974 ; Jacobs, Werth, and Guttman, 1984 ). Oligonucleotide primers (Vendramin et al., 1996 ) were used to PCR (polymerase chain reaction) amplify genomic DNA in 10-µL reactions (Jang, Genetic Analysis Bulletin number 30, LI-COR Inc., Lincoln, Nebraska, USA) containing: 20 ng genomic DNA, 10x PCR buffer (Boehringer Mannheim, Indianapolis, Indiana, USA) [50 mmol/L KCl, 10 mmol/L Tris-HCl (pH 8.3)], IRD-41 labeled forward strand primer (0.8 pmol) (LI-COR Inc., Lincoln, Nebraska), reverse strand primer (1.6 pmol), 0.1 mmol/L each dATP, dTTP, dCTP, and dGTP, MgCl2 (2 mmol/L), 1 unit Taq polymerase (Boehringer Mannheim, Indianapolis, Indiana, USA) and sterile deionized H2O. A "touchdown" program on an MJ Research (Watertown, Massachusetts, USA) PTC-100 Programmable Thermal Controller was used for amplification. After an initial denaturation at 94°C x 3 min, eight cycles with descending annealing temperatures (from 56°C x 30 s at 1.0°C per cycle) were followed by a constant annealing temperature of 55°C for 28 cycles. Extension was at 72°C x 30 s for all cycles, and a final 30-min extension time at 72°C was used to ensure that the terminal transferase reaction by the polymerase was carried out to completion. Then 5 µL of stop buffer (95% formamide, 0.1 mmol/L EDTA, and 0.1% bromophenol blue) were added to each well, and reactions were denatured at 94°C for 3 min prior to electrophoresis.


View this table:
[in this window]
[in a new window]
 
Table 1. Location and sources of Fraser (Abies fraseri), balsam (A. balsamea), and intermediate (A. b. var. phanerolepis) fir in the range-wide chloroplast microsatellite survey

 
PCR products were run on 0.25 mm 8% Long Ranger denaturing (7.5 mol/L urea) polyacrylamide gels on a LI-COR DNA 4000L automated sequencer at 50°C with constant current (35 mA) and power (70 W) and 16-bit image data were collected. A typical gel (Fig. 2) consisted of five molecular mass markers (LI-COR STR molecular mass marker), 42 sample lanes, and a water control (Clark, 1998 ). A primer producing a monomorphic band with a size similar to that of the polymorphic alleles was included in all reactions as an internal standard in each lane. DNA from Abies nordmanniana (J. C. Raulston Arboretum, North Carolina State University, Raleigh, North Carolina, Accession number 950295) was also included in two lanes to standardize fragment sizing within and between gels. DNA fragment sizes were analyzed with the software program RFLPScan (Scanalytics, Billerica, Massachusetts, USA). Repeatability of allele assignments was assessed in two independent PCR amplifications with the samples randomly rearranged on plates. PCR products from Pt 30204 were cloned and sequenced using a pGEM T vector (Promega, Madison, Wisconsin, USA) and dRhodamine dye labeled terminator reactions (GenBank accession numbers GBAN-AF 183871–AF 183875; the prefix GBAN- has been added to link the online version of American Journal of Botany but is not part of the actual accession number). PCR products at Pt 71936 from BESS (Base Excision Sequence Scanning) (Hawkins and Hoffman, 1997 ) reactions were visualized using infrared labeled primers on a LI-COR DNA 4000L automated sequencer.



View larger version (50K):
[in this window]
[in a new window]
 
Fig. 2. Allelic diversity in Fraser (Abies fraseri) and balsam (A. balsamea) fir with Pinus thunbergii primers: (a) Pt 30204, (b) Pt 71936. A monomorphic band at 131 bp (Pt 15169) is included on both gels as an internal standard. Lanes 1, 12, 24, 35, and 48 are LI-COR STR molecular weight markers. Lanes 13 and 36 are A. nordmanniana (A. nord.) samples included as internal standards. Lanes 2-11 and 14-23 are Fraser fir. Lanes 25-34 and 37-46 are balsam fir. Lane 47 is a negative control

 
Data analysis
Allele and haplotype frequencies within each taxon were estimated directly by counting. Gene diversity (He) was calculated with the formula He = 1 - ({Sigma}i=1n pi2) where n equals the number of alleles and pi equals the frequency of the ith allele in the population. Departure from homogeneity of allele frequencies among taxa was tested by {chi}2 analysis (e.g., Neter et al., 1996 ). Haplotype diversity was calculated in the same manner as gene diversity, with n and pi referring to haplotypes.

Population structure was assessed at individual loci with RST (Slatkin, 1995 ) as an index of genetic structure found within and among populations. RST, based on a stepwise mutation model (SMM), was calculated with the formula RST = (S - Sw)/S, where Sw is the sum over all loci of twice the weighted mean of the within population variances, and S is the sum over all loci of twice the variance of the combined populations (Microsat version 1.5b; Minch et al., 1995/1996 ).

An analysis of molecular variance (AMOVA) (Schneider et al., 1997 ) using FST as the measure of genetic distance was carried out to assess population structure of haplotypes within and among Fraser, balsam, and intermediate fir. Empirical significance of the AMOVA values was calculated with >25 000 permutations. Departure from random association of the two alleles in haplotypes was measured by exact tests for linkage disequilibrium (Fisher's exact test with 3200 permutations) using Genetic Data Analysis software (Lewis and Zaykin, 1998 ).

An estimate of genetic distance among individuals within populations was calculated with the mean pairwise haplotype distance measure D2SH, also based on a stepwise mutation approach. It is the squared sum of the pairwise absolute differences in allele sizes for both loci, averaged over all possible pairs of haplotypes and was computed with the formula:

where n equals the number of individuals, L is the number of polymorphic loci, aik and ajk are the allele sizes of the ith and jth individuals at the kth locus (Echt et al., 1998 ; Vendramin et al., 1998 ).

A chloroplast haplotype network, based on similar representations of parsimonious relationships in chloroplast microsatellites, or nonrecombining portions of the human Y chromosome, was constructed for the three taxa (Cooper et al., 1996 ; Doyle et al., 1998 ; Morgante, Felice, and Vendramin, 1998 ). Haplotypes were plotted directly with allele sizes found at Pt 30204 (vertical axis) and Pt 71936 (horizontal axis) and then connected to their nearest neighbors.

RESULTS

The chloroplast genomes of Fraser, balsam, and intermediate fir contain 12 sites that can be PCR amplified using the P. thunbergii microsatellite primers (Table 2). Of the 12 primer pairs, all amplified consistently under standard conditions except PT 9383, which showed poor and inconsistent amplification. All PCR reactions produced a single, robust band per primer pair; thus, there was no evidence for heteroplasmy. Fragments that differed in size by 1 bp could be reliably distinguished and assigned to allele categories for the same samples run on different gels. The mean absolute difference in estimated fragment size of gels was <0.6 bp (Pt 71936: X = 0.20 bp, SD = 0.31; Pt 30204: X = 0.5 bp, SD = 0.36). PCR products of Pt 30204 from five individuals representing the three taxa were sequenced and the presence of a (dA)n (dG)n site similar to that found in P. thunbergii was confirmed. Repeat-length variation both at this site and in a (dT)n site within the same fragment was observed. There was 91% conserved identity between the P. thunbergii and the Abies sequences.


View this table:
[in this window]
[in a new window]
 
Table 2. Pinus thunbergii (PT) primers that PCR amplify in Fraser (Abies fraseri), balsam (A. balsamea), and intermediate (A. b. var. phanerolepis) fir under standard conditions. Alleles are those found in the 78 individuals from the range-wide survey

 
The gene diversity levels for two polymorphic loci, Pt 30204 and Pt 71936, ranged from 0.65 to 0.84 across the fir taxa (Table 3). The gene diversity level at Pt 63718, the third polymorphic locus, was 0.09 for all trees. Pt 63718 was not included in subsequent analyses with the two hypervariable sites. There was significant gene frequency differentiation among the fir taxa. {chi}2 tests for departure from homogeneity in allele frequency at Pt 30204 showed highly significant differentiation between Fraser and intermediate, Fraser and balsam, and balsam and intermediate fir (P <= 0.001). At Pt 71936, there was a highly significant difference in allele frequency between Fraser and intermediate fir (P <= 0.001), and a significant difference (P <= 0.05) between balsam and Fraser, and balsam and intermediate fir. A high degree of population subdivision was shown in RST for Pt 30204 (0.36) and a much lower value for Pt 71936 (0.10) (Table 4).


View this table:
[in this window]
[in a new window]
 
Table 3. Gene diversity in the range-wide samples of Fraser (Abies fraseri), balsam (A. balsamea), and intermediate (A. b. var. phanerolepis) fir (26 samples per taxon)

 

View this table:
[in this window]
[in a new window]
 
Table 4. Distribution of chloroplast microsatellite diversity in Fraser (Abies fraseri), balsam (A. balsamea), and intermediate (A. b. var. phanerolepis) fir at Pt 30204 and Pt 71936 in the range-wide survey

 
The high level of allelic diversity observed for the single locus analysis was also reflected in the values obtained for haplotype diversity and D2SH, a distance-based measure of within-population haplotypic variation (Table 5). There were 38 unique haplotypes in the sample of 78 individuals representing the three taxa. Balsam fir had the highest diversity with 20 haplotypes, Fraser fir with 16, and intermediate fir with nine (Table 5). Seven haplotypes occurred in two taxa, four occurred in both balsam and Fraser fir, and three were shared between balsam and intermediate fir. None of the shared haplotypes occurred in all three taxa (Table 6), and haplotype frequencies exhibited minimal overlap between taxa. In balsam fir, only two haplotypes occurred with a frequency >0.05. Six haplotypes in Fraser fir and five haplotypes in intermediate fir were seen with a frequency >0.05. Haplotype diversity was highest in balsam fir in the range-wide sample and lowest in intermediate fir, though all levels were >0.75. Values of D2SH were >10 for each of the three taxa. Balsam fir showed a much higher distance between individuals within the taxon than either Fraser or intermediate fir, reflecting a higher level of genetic variation in the balsam fir. Intermediate fir (P <= 0.001) and balsam fir (P <= 0.055) showed departure from random association of alleles at the two loci. Highly significant haplotypic differentiation was detected in an AMOVA analysis of haplotypic differentiation (FST = 0.06, P < 0.001, IAM, infinite alleles model).


View this table:
[in this window]
[in a new window]
 
Table 5. Haplotype diversity in Fraser (Abies fraseri), balsam (A. balsamea), and intermediate (A. b. var. phanerolepis) fir in the range-wide survey (26 samples per taxon)

 

View this table:
[in this window]
[in a new window]
 
Table 6. Chloroplast haplotype frequencies for Fraser (Abies fraseri), balsam (A. balsamea), and intermediate (A. b. var. phanerolepis) fir in the range-wide survey at Pt 30204 and Pt 71936. (The most abundant haplotype in each taxon and haplotypes shared between taxa are indicated in boldface type.)

 
Two distinct clusters of haplotypes among Fraser, balsam, and intermediate fir were apparent in the chloroplast haplotype network (Fig. 3). In the two clusters seen in this network, 70 out of the total 78 trees, or ~90% of individuals, were connected to other haplotypes by one mutation or fewer. Four mutations separated the two groups of interconnected individuals, and there were nine mutations separating the most common haplotype in Fraser fir from the most common haplotype in balsam fir. In the upper cluster there were 38 individuals connected to each other by one or fewer mutational changes in allele size. This group included ~62% of the balsam fir samples and 73% of the intermediate. In the lower cluster, ~77% of Fraser fir and 23% of intermediate fir were represented. The two most common Fraser fir haplotypes were also found at low frequency in the upper cluster, and the most common balsam fir haplotype was found at low frequency in the lower cluster. Three individuals were located at least five steps away from the majority of the trees. These individuals were found in the extreme eastern and western portions of the range of balsam fir, Alberta and Newfoundland, Canada. Intermediate fir was found predominantly (18 out of 26 individuals) in the network containing the majority of the balsam fir samples.



View larger version (74K):
[in this window]
[in a new window]
 
Fig. 3. Chloroplast haplotype network showing stepwise differences in allele size among Fraser (Abies fraseri), balsam (A. balsamea), and intermediate (A.b. var. phanerolepis) fir. Haplotypes are separated by the sum of the absolute value of the difference in size of the two polymorphic repeat units. Lines joining haplotypes indicate a single nucleotide difference unless otherwise indicated. B = balsam, F = Fraser, and I = intermediate. Allele sizes (bp) for Pt 71936 on the horizontal axis and Pt 30204 on the vertical axis define the haplotypes

 
DISCUSSION

Chloroplast microsatellites at two highly variable loci provided clear evidence for genetic divergence in eastern North American Abies. The high level of population subdivision at Pt 30204 (RST > 0.3) suggests that recent gene flow between the geographically separated fir taxa is limited. If gene flow via pollen were extensive between stands, little population subdivision would be expected in paternally inherited DNA. Strong among-population subdivision at chloroplast microsatellite loci has also been demonstrated in P. halepensis (RST > 0.2) and P. brutia (RST > 0.3) at nine chloroplast microsatellite loci (Bucci et al., 1998 ). More extensive population subdivision is typically found in maternally inherited mitochondrial genes in conifers, due to limited dispersal via seeds (e.g., Latta et al., 1998 ; Strauss, Hong, and Hipkins, 1993 ). Chloroplast DNA is maternally inherited in angiosperm species and spatial clustering of haplotypes is often pronounced related to limited seed dispersal. Petit et al. (1993) demonstrated that genetic markers for cpDNA in oak have highly significant geographic variation and high genetic differentiation (GST = 0.895). Petit et al. (1997) showed large patches (hundreds of square kilometers) of oak forest in France that had a single cpDNA marker haplotype. Other analyses of Quercus robur and Q. petrea in Europe have also shown nonrandom patterns of cpDNA variation that likely arose through the "resampling" process involved in northward migration of tree populations during the Holocene warming 7000 yr ago (e.g., Dumolin-Lapegue et al., 1997 ; Johnk and Siegismund, 1997 ; Ferris et al., 1998 ).

Studies of chloroplast microsatellites in European Abies found high levels of diversity at the loci used in this study (Vendramin and Ziegenhagen, 1997 ; Ziegenhagen et al., 1998 ; Vendramin et al., 1999 ). An examination of sequence variation in European Abies (Vendramin and Ziegenhagen, 1997 ) showed a highly conserved sequence at Pt 71936. Length variation at Pt 30204 in European Abies was found in two different mononucleotide units within the DNA fragment, and in insertions and deletions in the flanking sequence. In A. alba, Pt 30204 appeared to be a "mutational hotspot" (Vendramin et al., 1999 ). DNA sequence of Pt 30204 in eastern North American Abies also had two mononucleotide repeat units and insertions/deletions in the flanking sequence. Mutational processes at this locus may be different than those found in Pt 71936. Currently, little is known about the processes or rates of mutation at chloroplast microsatellite loci. Predominant paternal inheritance of Abies cpDNA was demonstrated in a study of European Abies species (Vendramin and Ziegenhagen, 1997 ).

Fraser, balsam, and intermediate fir individuals within populations showed high levels of genetic variation at the level of haplotypes reflected by the elevated values found for D2SH. Balsam fir, which had the largest range, had the largest value of D2SH, 26.17. Haplotypes from Alberta, which were distinct from the balsam haplotype cluster in the network analysis, contributed to the large value. Intermediate fir, which occurs in small, isolated stands, had the lowest value, 13.90. In A. alba, values of D2SH were also high, 28.08 (Vendramin 1999 , personal communication, Instituto Miglioramento Genetico Piante Forestali). Values of mean D2SH for chloroplast microsatellite loci in P. pinaster averaged 0.51, and in P. resinosa ranged from 0.05 to 0.44 (Echt et al., 1998 ; Vendramin et al., 1998 ). Within populations of P. halepensis, D2SH ranged from 0.04 to 14.41 and showed evidence of two major branches within the taxon (Morgante, Felice, and Vendramin, 1998 ).

A network analysis has the potential to reveal some of the haplotypic structure that results from ancestral relationships and the nonrandom distribution of mutations among lineages. Ancestry is not taken into account by analyses based on pairwise haplotypic distance measures. The construction of the haplotype networks (Goldstein et al., 1995 ; Slatkin, 1995 ) assumes that microsatellites add or subtract repeat units with equal probability and, in most cases, one repeat unit at a time (i.e., the stepwise mutation model, SMM). Therefore, alleles of a similar size could have a history of fewer intervening mutations and contain information about time since divergence from a common ancestor. However, identical states of microsatellite alleles can arise independently, and the potential for homoplasy must be considered when assessing relationships at these sites (Doyle et al., 1998 ). Fast-evolving microsatellite loci on the human Y-chromosome have provided nonrecombining haplotypes that have been useful for interpreting the expansion and migration of human populations during the past 200 000 yr (e.g., Wilson and Balding, 1998 ). The mutation rate for chloroplast microsatellite loci could be low, ~10-5 (Provan et al., in press ), which would suggest that the discontinuity in fir chloroplast haplotypes could predate the most recent glacial episode. Further research is needed to evaluate a coalescent-with-stepwise-mutation model for the nonrecombining fir microsatellite haplotypes as an explanation for the two groups revealed by the fir haplotype network analysis. The current distribution of fir haplotypes could also be related to environmental gradients and that possibility needs to be investigated as well.

The geographic distribution of the two groups of haplotypes in eastern North American fir showed a nonrandom pattern and provides evidence for the divergence of the taxa by revealing two evolutionary lineages. The number of population-specific haplotypes (balsam 65%, Fraser 75%, intermediate 55%) could make these markers useful for discrimination among the taxa. A study of 714 A. alba individuals from 17 populations revealed 90 haplotypes with evidence for spatial distribution of specific haplotypes (Vendramin et al., 1999 ). Thirteen of the 17 populations in the study of A. alba were characterized by at least one population-specific haplotype. Population-specific haplotypes were found in 88% of the populations analyzed in the P. halepensis species complex (Bucci et al., 1998 ). The majority of haplotypes in this study, with small sample sizes, were detected with a frequency <0.05. Further study and larger sample sizes are needed to characterize the distribution of low-frequency haplotypes within populations and the presence of population-specific haplotypes.

Some authorities believe that the divergence in eastern North American Abies taxa began after the end of the most recent glacial phase (Oosting and Billings, 1951 ; Myers and Bormann, 1963 ). Eastern North American Abies were present as a continuous population 8000 yr ago along the Appalachian Mountains from Georgia to Newfoundland, but the once uninterrupted range was broken into discontinuous populations with warming during the early and mid-Holocene 7000 yr ago (Delcourt and Delcourt, 1987 ). Glacial phases account for 90% of the 2 million yr period of the Quaternary. Interglacial periods were relatively brief, lasting 10 000 to 15 000 yr with temperate climates similar to that of the present (Davis, 1983 ). The last interglacial period prior to the present began 125 000 yr ago and lasted about 15 000 yr. This period was followed by brief cycles of warming and cooling until the beginning of the most recent glacial period 70 000 yr ago. This last interval culminated in a glacial maximum 18 000 to 20 000 yr ago, and much of the current range of balsam fir was under ice (Delcourt and Delcourt, 1987, 1991 ). During this interval, many tree species most likely persisted in small, isolated population centers in refugial areas (Davis, 1983 ). Evidence suggests that tree species adapt very slowly to changing climate and are forced to migrate to survive (Bradshaw and McNeilly, 1991 ). In recent times, migration and population history of Abies taxa may have contributed a larger part to current patterns of diversity than has natural selection.

Regardless of the taxonomic identity of Fraser, balsam, and intermediate fir, we can recognize evolutionary divergence in cpDNA haplotypes among populations. The almost reciprocal monophyly present in haplotypes between balsam and Fraser fir provides justification for recognizing these taxa as independent Management Units according to the definition by Moritz (1994) : "Populations with significant divergence of allele frequencies at nuclear or mitochondrial loci, regardless of the phylogenetic distinctiveness of the alleles." The naturally fragmented population centers make it critically important to understand the distribution of genetic diversity within these three taxa for the development of effective conservation strategies for Fraser fir. The highly variable repeat sequences in chloroplast genomes could reveal population history and phylogeographic structure for plant taxa similar to that revealed for animal taxa in studies of mitochondrial DNA sequence variation (Avise, 1994 ; Avise and Hamrick, 1996 ).

Summary
In this paper, we have shown a strong discontinuity in chloroplast microsatellite haplotypes of eastern North American Abies taxa. The level of divergence found for two chloroplast microsatellite loci was greater than previously shown in allozymes, terpenes, and morphological characteristics (Robinson and Thor, 1969 ; Zavarin and Snajberk, 1972 ; Jacobs, Werth, and Guttman, 1984 ). These markers exhibit a high degree of among-population differentiation and significant difference in allele frequency. Population-specific haplotypes were found in all of the three taxa and could provide information useful for conservation planning, although larger sample sizes are necessary to evaluate haplotype distribution. The high level of variability found at these loci could provide a method to characterize pollen dispersal among and between populations. The large between-taxon variation also suggests geographical structure that contrasts with earlier intergrading variation. These markers show a potential, as demonstrated in larger studies of European Abies populations, to delineate species boundaries and populations with a high value for gene conservation in closely related taxa.

FOOTNOTES

1 The authors thank Giovanni Vendramin for comments on an earlier version of this manuscript; J. Frampton, NCSU, G. Hawley, UVM, L. De Verno, CFS, New Brunswick, L. Corriveau and D. Stepnisky, Weyerhaeuser Canada, J. Brown, OSU, N. Dhir, Alberta Environmental Protection, F. C. Yeh, Univ. Alberta, T. Blount, Shenandoah National Park, and K. Garrard for their assistance in obtaining samples for this study; and the Forest Biotechnology Group and the Department of Botany at North Carolina State University, and a grant-in-aid from Sigma Xi for providing funding. This work represents a portion of a Master's thesis submitted to the Department of Botany at North Carolina State University. Back

4 Author for correspondence, current address: North Carolina State University, Forest Biotechnology Group, Campus Box 7247, Raleigh, North Carolina 27695 USA. (e-mail: catherine_clark{at}ncsu.edu ). Back

LITERATURE CITED

Avise, J. C. 1994 Molecular markers, natural history and evolution. Chapman and Hall, New York, New York, USA.

———, and J. L. Hamrick [eds.]. 1996 Conservation genetics, case histories from nature. Chapman and Hall, New York, New York, USA.

Bradshaw, A. D., and T. McNeilly. 1991 Evolutionary response to global climate change. Annals of Botany 67 (Supplement1): 5–14.

Bucci, G., M. Anzidei, A. Madaghiele, and G. G. Vendramin. 1998 Detection of haplotypic variation and natural hybridization in halepensis-complex pine species using chloroplast simple sequence repeat (SSR) markers. Molecular Ecology 7: 1633–1643.[CrossRef]

Cato, S. A., and T. E. Richardson. 1996 Inter- and intraspecific polymorphism at chloroplast SSR loci and the inheritance of plastids in Pinus radiata D. Don. Theoretical and Applied Genetics 93: 587–592.[CrossRef][Web of Science]

Clark, C. M. 1998 Chloroplast microsatellite differentiation in eastern North American Abies Mill. Master's thesis, North Carolina State University, Raleigh, North Carolina, USA.

Clegg, M. T., and G. Zurawski. 1992 Chloroplast DNA and the study of plant phylogeny: present status and future prospects. In P. E. Soltis, D. E. Soltis, and J. J. Doyle [eds.], Molecular systematics of plants, 1–13. Chapman and Hall, New York, New York, USA.

Cooper, G., W. Amos, D. Hoffman, and D. C. Rubinsztein. 1996 Network analysis of human Y microsatellite haplotypes. Human Molecular Genetics 5: 1759–1766.[Abstract/Free Full Text]

Davis, M. B. 1983 Quaternary history of deciduous forests of eastern North America and Europe. Annals of the Missouri Botanical Garden 70: 550–563.

Delcourt, H. R., and P. A. Delcourt. 1991 Quaternary ecology: a paleoecological perspective. Chapman and Hall, London, UK.

Delcourt, P. A., and H. R. Delcourt. 1987 Long-term forest dynamics of the temperate zone. Springer-Verlag, New York, New York, USA.

Doyle, J. J., and J. I. Doyle. 1990 Isolation of plant DNA from fresh tissue. Focus 12: 13–15.

———, M. Morgante, S. V. Tingey, and W. Powell. 1998 Size homoplasy in chloroplast microsatellites of wild perennial relatives of soybean (Glycine Subgenus Glycine). Molecular Biology and Evolution 15: 215–218.[Web of Science][Medline]

Dull, C. W., J. D. Ward, H. D. Brown, G. W. Ryan, W. H. Clerke, and R. J. Uhler. 1988 Evaluation of Spruce and Fir Mortality in the Southern Appalachian Mountains, United States Department of Agriculture, Forest Service, Southern Region.

Dumolin-Lapegue, S., B. Demesure, S. Fineschi, V. LeCorre, and R. J. Petit. 1997 Phylogeographic structure of white oaks throughout the European continent. Genetics 146: 1475–1487.[Abstract]

Echt, C., L. DeVerno, M. Anzidei, and G. Vendramin. 1998 Chloroplast microsatellites reveal population genetic diversity in red pine, Pinus resinosa Ait. Molecular Ecology 7: 307–316.

Ferris, C., R. A. King, R. Vainola, and G. M. Hewitt. 1998 Chloroplast DNA recognizes three refugial sources of European oaks and suggests independent eastern and western immigrations to Finland. Heredity 80: 584–593.

Flora of North America Editorial Committee. 1993 Flora of North America, vol. 2, Pteridophytes and gymnosperms. Oxford University Press, New York, New York, USA.

Garza, J. C., M. Slatkin, and N. B. Freimer. 1995 Microsatellite allele frequencies in humans and chimpanzees, with implications for constraints on allele size. Molecular Biology and Evolution 124: 594–603.

Goldstein, D. B., A. R. Linares, L. L. Cavalli-Sforza, and M. W. Feldman. 1995 An evaluation of genetic distances for use with microsatellite loci. Genetics 139: 463–471.[Abstract]

Hamrick, J. L., and M. J. W. Godt. 1990 Allozyme diversity in plant species. In A. H. D. Brown, M. T. Clegg, A. L. Kahler, and B. S. Weir [eds.], Plant population genetics, breeding, and genetic resources, 43–82. Sinauer, Sunderland, Massachusetts, USA.

———, and ———. 1996 Conservation genetics of endemic plant species. In J. C. Avise and J. L. Hamrick [eds.], conservation genetics: case histories from nature, 281–304. Chapman and Hall, New York, New York, USA.

Hawley, G. J., and D. H. DeHayes. 1985 Hybridization among several North American firs. I. Crossability. Canadian Journal of Forest Research 15(1): 42–49.

Hawkins, G. A., and L. M. Hoffman. 1997 Base excision sequence scanning. Nature Biotechnology 15: 803–804.[CrossRef][Web of Science][Medline]

Jacobs, B. F., C. R. Werth, and S. I. Guttman. 1984 Genetic relationships in Abies (fir) of eastern United States: an electrophoretic study. Canadian Journal of Botany 62: 609–616.

Johnk, N., and H. R. Siegismund. 1997 Population structure and post-glacial migration routes of Quercus robur and Quercus petraea in Denmark, based on chloroplast DNA analysis. Scandinavian Journal of Forest Research 12: 130–137.[Web of Science]

Latta, R. G., Y. B. Linhart, D. Fleck, and M. Elliot. 1998 Direct and indirect estimates of seed versus pollen movement within a population of ponderosa pine. Evolution 52: 61–67.[CrossRef][Web of Science]

Lewis, P. O., and D. Zaykin. 1998 Genetic data analysis software. Computer program for the analysis of allelic data. Version 1.0. http://alleyn.eeb.ucon.edu/gda/.

Liu, T. 1971 A monograph of the genus Abies. Department of Forestry, National Taiwan University, Tapei, Taiwan, China.

McCauley, D. E. 1995 The use of chloroplast DNA polymorphism in studies of gene flow in plants. Trends in Ecology and Evolution 10(5): 198–202.

Minch, E., A. Ruiz-Linares, D. Goldstein, M. Feldman, and L. L. Cavalli-Sforza. 1995/1996 Microsat Software, Stanford, California, USA.

Morgante, M., N. Felice, and G. G. Vendramin. 1998 Analysis of hypervariable chloroplast microsatellites in Pinus halepensis reveals a dramatic genetic bottleneck. In A. Karp, P. G. Isaac, and D. S. Ingram [eds.], Molecular tools for screening biodiversity, 407–412. Chapman and Hall, London, UK.

———, and A. M. Olivieri. 1993 PCR—amplified microsatellites as markers in plant genetics. Plant Journal 3: 175–182.[CrossRef][Web of Science][Medline]

Moritz, C. 1994 Defining ‘Evolutionarily Significant Units' for conservation. Trends in Ecology and Evolution 9(10): 373–375.

Myers, O., and F. H. Bormann. 1963 Phenotypic variation in Abies balsamea in response to altitudinal and geographic gradients. Ecology 44: 429–436.[CrossRef][Web of Science]

Neale, D. B., and R. R. Sederoff. 1989 Paternal inheritance of chloroplast DNA and maternal inheritance of mitochondrial DNA in loblolly pine. Theoretical and Applied Genetics 77: 212–216.[Web of Science]

———, N. C. Wheeler, and R. W. Allard. 1986 Paternal inheritance of chloroplast DNA in Douglas fir. Canadian Journal of Forest Research 16: 1152–1154.[CrossRef]

Neter, J., M. H. Kutner, C. J. Nachtsheim, and W. Wasserman. 1996 Applied Linear Statistical Models. Irwin, Chicago.

Oldfield, S., C. Lusty, and A. MacKinven [eds.]. 1998 The world list of threatened trees. World Conservation Press, Cambridge, UK.

Oosting, H. J., and W. E. Billings. 1951 A comparison of virgin spruce-fir forest in the northern and southern Appalachian system. Ecology 32: 84–103.[CrossRef][Web of Science]

Petit, R. J., A. Kremer, and D. B. Wagner. 1993 Geographic structure of chloroplast DNA polymorphisms in European oaks. Theoretical and Applied Genetics 87: 122–128.[Web of Science]

———, E. Pineau, B. Demesure, R. Bacilieri, A. Ducousso, and A. Kreme. 1997 Chloroplast DNA footprints of postglacial colonization by oaks. Proceedings of the National Academy of Sciences, USA 94: 9996–10 001.[Abstract/Free Full Text]

Powell, W., M. Morgante, R. McDevitt, G. G. Vendramin, and J. A. Rafalski. 1995a Polymorphic simple sequence repeat regions in chloroplast genomes: applications to the population genetics of pines. Proceedings of the National Academy of Sciences, USA 92: 7759–7763.[Abstract/Free Full Text]

———, ———, C. Andre, J. W. McNicol, G. C. Machray, J. J. Doyle, S. V. Tingey, and J. A. Rafalski. 1995b Hypervariable microsatellites provide a general source of polymorphic DNA markers for the chloroplast genome. Current Biology 5: 1023–1029.[CrossRef][Web of Science][Medline]

Primmer, C. R., A. P. MØller, and H. Ellegren. 1996 A wide-range survey of cross-species microsatellite amplification in birds. Molecular Ecology 5: 365–378.[CrossRef][Medline]

Provan, J., N. Soranzo, N. Wilson, D. Goldstein, and W. Powell. In press A low mutation rate for chloroplast microsatellites. Genetics.

Robinson, J. F., and E. Thor. 1969 Natural variation in Abies of the Southern Appalachians. Forest Science 15: 238–245.[Web of Science]

Schneider, S., J-M. Kueffer, D. Roessli, and L. Excoffier. 1997 Arlequin: a software for population genetic data analysis. Genetics and Biometry Lab, Department of Anthropology, University of Geneva, Geneva, Switzerland.

Sederoff, R. R., C. S. Levings III, D. H. Timothy, and W. W. L. Hu. 1981 Evolution of DNA sequence organization in mitochondrial genomes of Zea. Proceedings of the National Academy of Sciences, USA 78: 5953–5957.[Abstract/Free Full Text]

Slatkin, M. 1995 A measure of population subdivision based on microsatellite allele frequencies. Genetics 139: 457–462.[Web of Science][Medline]

Steinkellner, C., C. Lexer, E. Turetschek, and J. Glössl. 1997 Conservation of (GA)n microsatellite loci between Quercus species. Molecular Ecology 6: 1189–1194.[CrossRef]

Strauss, S. H., Y-P. Hong, and V. D. Hipkins. 1993 High levels of population differentiation for mitochondrial DNA haplotypes in Pinus radiata, muricata, and attenuata. Theoretical and Applied Genetics 86: 605–611.

Thor, E., and P. E. Barnett. 1974 Taxonomy of Abies in the Southern Appalachians: variation in balsam monoterpenes and wood properties. Forest Science 20: 32–40.

Tsumura, Y., and Y. Suyama. 1998 Differentiation of mitochondrial DNA polymorphisms in populations of five Japanese Abies species. Evolution 52: 1031–1042.[CrossRef][Web of Science]

Vendramin, G. G., M. Anzidei, A. Madaghiele, and G. Bucci. 1998 Distribution of genetic diversity in Pinus pinaster Ait. as revealed by chloroplast microsatellites. Theoretical and Applied Genetics 97: 456–463.[CrossRef][Web of Science]

———, B. Degen, R. J. Petit, M. Anzidei, A. Madaghiele, and B. Ziegenhagen. 1999 High level of variation at Abies alba chloroplast microsatellite loci in Europe. Molecular Ecology 8: 1117–1126.[CrossRef][Medline]

———, L. Lelli, P. Rossi, and M. Morgante. 1996 A set of primers for the amplification of 20 chloroplast microsatellites in Pineaceae. Molecular Ecology 5: 595–598.[CrossRef][Medline]

———, and B. Ziegenhagen. 1997 Characterization and inheritance of polymorphic plastid microsatellites in Abies. Genome 40: 857–864.[Medline]

Wilson, I. J., and D. J. Balding. 1998 Geneological inference from microsatellite data. Genetics 150: 499–510.[Abstract/Free Full Text]

Wu, J. Y., K. V. Krutovskii, and S. H. Strauss. 1998 Abundant mitochondrial genome diversity, population differentiation and convergent evolution in pines. Genetics 150: 1605–1614.[Abstract/Free Full Text]

Zavarin, E., and K. Snajberk. 1972 Geographical variability of monoterpenes from Abies balsamea and A. fraseri. Phytochemistry 11: 1407–1421.

Ziegenhagen, B., F. Scholz, A. Madaghiele, and G. Vendramin. 1998 Chloroplast microsatellites as markers for paternity analysis in Abies alba. Canadian Journal of Forest Research 28: 317–321.




This article has been cited by other articles:


Home page
ANN BOT (LOND)Home page
A. C. S. Ramos, J. P. Lemos-Filho, R. A. Ribeiro, F. R. Santos, and M. B. Lovato
Phylogeography of the Tree Hymenaea stigonocarpa (Fabaceae: Caesalpinioideae) and the Influence of Quaternary Climate Changes in the Brazilian Cerrado
Ann. Bot., November 1, 2007; 100(6): 1219 - 1228.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Bot.Home page
A. Terrab, O. Paun, S. Talavera, K. Tremetsberger, M. Arista, and T. F. Stuessy
Genetic diversity and population structure in natural populations of Moroccan Atlas cedar (Cedrus atlantica; Pinaceae) determined with cpSSR markers
Am. J. Botany, September 1, 2006; 93(9): 1274 - 1280.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (19)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Clark, C. M.
Right arrow Articles by O'Malley, D. M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Clark, C. M.
Right arrow Articles by O'Malley, D. M.
Agricola
Right arrow Articles by Clark, C. M.
Right arrow Articles by O'Malley, D. M.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS