Department of Biology, University of Oulu, Finland
![]() |
Abstract |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Introduction |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
In this study, we examine nucleotide variation in Scots pine (Pinus sylvestris L.). It is a long-lived predominantly outcrossing perennial with continuous distribution that ranges from Scotland to eastern Siberia and from northern Finland to Turkey (fig. 1 ). It has high variability at enzyme loci, RFLPs, and microsatellites (Gullberg et al. 1985
; Muona and Harju 1989
; Müller-Stark, Baradat, and Bergmann 1992; Karvonen and Savolainen 1993
; Goncharenko, Silin, and Padutov 1994
; Karhu et al. 1996
). In 10 allozyme studies, the average heterozygosity values per population ranged from 0.25 to 0.39 for Scots pine (Müller-Stark, Baradat, and Bergmann 1992), among the highest for the European tree species compared. These values are high compared with many other groups, such as Drosophila allozyme heterozygosity with an average of 0.14. The differentiation between populations at marker loci is low within Scandinavia (FST = 0.02) (Karhu et al. 1996
) and even between Scandinavia and the eastern part of the range (Wang, Szmidt, and Lindgren 1991
). There is much recombination in the genome of pines: the genetic genome size of Pinus taeda is 1,7001,800 cM, likely to be similar to that in many pines (Remington et al. 1999
). No linkage disequilibrium was found between enzyme loci either in pollen grains or female gametes (Muona and Szmidt 1985
). Scots pine has very high early and late inbreeding depression (Koski 1971
; Kärkkäinen, Koski, and Savolainen 1996). In a partially selfing species, most deleterious alleles should be rapidly eliminated. The most reasonable way to account for high inbreeding depression in a partially selfing species is a high genome-wide mutation rate to deleterious alleles (Lande, Schemske, and Schultz 1994
), and Scots pine seems to conform to this model (Koelewijn, Koski, and Savolainen 1999
).
|
However, earlier data do not allow us to make a well-justified prediction on the range of linkage disequilibrium in the genome. The level of disequilibrium between neutral genes is governed by the product 4Nec, where c is the recombination rate between loci (Hill and Robertson 1968
). Close linkage or inbreeding will restrict effective recombination. The larger this product, the less disequilibrium we expect between neutral loci. The current population size is very large and there is much recombination between genes, but we might still observe effects of historical bottlenecks between closely linked nucleotide sites. The current population structure of Scots pine in Europe appears to have been formed mainly by the postglacial migration of the species from refugia located in southwestern, southern, and southeastern parts of its present range. This expansion took place about 7,0009,000 years ago (Hyvärinen 1987
). Scots pine mitochondrial DNA variants provide further evidence for multiple origins after the last glaciation from the south: through western Europe via France, Germany, and Denmark and from the northeast via Finland (Sinclair, Morman, and Ennos 1999
; Soranzo et al. 2000
). Are we able to detect traces of the glacial history at the nucleotide level, or has ample gene flow because of pollen migration (Koski 1970
) and abundant recombination already eliminated these traces? Are the populations at mutation drift equilibrium, or do they still show traces of the possible glacial bottlenecks, in disequilibria between closely linked nucleotide sites or in frequency spectra that are predicted in expanding populations (Harpending 1994
).
We studied nucleotide polymorphism in the gene encoding phenylalanine ammonia-lyase (Pal, EC 4.3.1.5), a key enzyme in the secondary metabolism of higher plants catalyzing the nonoxidative conversion of l-phenylalanine into transcinnamate, the initial substrate of the phenylpropanoid pathway in wood formation. In Scots pine, Pal is believed to be related to ozone tolerance (Rosemann, Heller, and Sandermann Jr. 1991
), defense against pathogens (Lange et al. 1994
), and metabolism of exogenous compounds (Laukkanen and Sarjala 1997
). cDNA sequences of pal have been reported for a number of angiosperms (e.g., Yamada et al. 1992
). In angiosperms this gene usually consists of two exons and one intron. In gymnosperms, pal was first reported from loblolly pine, P. taeda, as a single-copy gene with no introns (Whetten and Sederoff 1992
). Later, Butland, Chow, and Ellis (1998)
sequenced a 366-bplong conserved region of pal from Pinus banksiana and discovered that pal is actually a multigene family of at least 810 loci.
Our results first show that Scots pine has rather low nucleotide diversity at synonymous sites. Following this finding, we examined whether other aspects of the nucleotide data are consistent with the predictions for a species with a large effective population size.
![]() |
Materials and Methods |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
DNA Isolation, Amplification, and Sequencing
Total genomic DNA from the megagametophytes was purified either by CTAB method (Doyle and Doyle 1990
) or by FastDNA® Kit (BIO101). On the basis of the P. taeda pal-gene nucleotide sequence (one exon of 2,265 bp, GenBank U39792), we designed PCR-primers that would amplify the pal1 gene in two overlapping parts in P. sylvestris. To amplify the other P. sylvestris genes used in this study, we designed primers on the basis of P. sylvestris and P. taeda gene sequences from GenBank (table 1
). Short fragments were directly sequenced, first using AmpliTaq Gold (PE Applied Biosystems) for the PCR. The PCR products were separated in 1% agarose gels (0.5 TBE), from which they were purified by centrifugation through aquarium filter fibers. Direct sequencing was done from both strands with BigDye-kit (PE Applied Biosystems). For cloning of longer parts of the genes, PCR was done with Dynazyme EXT (with proofreading property, Finnzyme). The amplification products were first cloned with the TOPO TA cloning kit (version J, Invitrogen) and then sequenced (BigDye, PE Applied Biosystems). The PCR-program for amplification of the pal1 gene reactions with AmpliTaq Gold-enzyme was: 94°C 12 min; 94°C 1.10 min, 53°C 1 min, 72°C 1.50 min, 30 cycles; 72°C 10 min, and for Dynazyme-enzyme: 94°C 4 min; 94°C 1.10 min, 53°C 1 min, 72°C 1 min, 30 cycles; 72°C 10 min. Three cloned products were sequenced in order to avoid errors made by the DNA polymerase. If nucleotide differences between the three clones were detected, then the majority nucleotide was taken to be the correct one in that position. In both sequencing methods, we used several sequencing primers designed on the basis of the P. taeda sequence and located every 300400 bp. Newly determined P. sylvestris sequences are deposited in the GenBank database and the accession numbers are listed in table 1
.
|
![]() |
Results |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Among the sequenced 2,045 sites in 20 alleles, there were 12 segregating sites (table 2
). Nine of the polymorphisms were synonymous and three were nonsynonymous. A transition of C to T at site 243 resulted in a change of alanine to valine in two individuals of the Kolari population. In the Kirov population, a transversion of T to G at site 762 resulted in a change of valine into glycine in one individual, and a transversion of A to T at site 915 alters histidine to leucine. These amino acid changes are considered to be conservative for protein function (Graur and Li 2000
, p. 20).
|
|
All polymorphic sites had only two alternative nucleotides, and four of the sites were singletons. As the populations showed so little divergence, we combined them into one set for Tajima's test. A significant negative Tajima's test could suggest either purifying selection or recent population expansion. However, the value of -0.56 was not statistically significantly different from zero.
Differentiation Between Populations
The genetic differentiation between the populations as estimated with AMOVA (table 4
) was very low. If the population from Spain was excluded, the value of FST was only 0.017, but if it was included in the calculations, the FST increased to 0.110. Statistically significant pairwise differentiation was found only between the most distant populations from Puebla de Lillo and Kirov (FST = 0.310).
|
![]() |
Discussion |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
|
Our results on rather low diversity could be because of several factors that we can examine. First, the effective population size (at least for this locus) could be lower than that assumed. Second, the mutation rate could be lower than that suggested by the inbreeding depression data. Third, the synonymous sites may in fact not be fully neutral, but under directional selection for major codons, e.g., for reasons of translational efficiency (see Akashi 1997
). Further, the low variability could be due to selective forces at other loci that reduce levels of variation, such as hitchhiking (Maynard Smith and Haigh 1974
; Kaplan, Hudson, and Langley 1989
), background selection (Charlesworth, Morgan, and Charlesworth 1993
), or to pseudohitchhiking due to alternating direction of selection in large populations (Gillespie 2000
). The high generation time and short time after glaciation bottlenecks may also mean that nucleotide polymorphism is not at equilibrium, even if our small data set did not show evidence of this (see below).
Large Effective Population Size and Mutation Drift Equilibrium
Many aspects of the data are consistent with a large effective population size. The differentiation between populations separated by thousands of kilometers was very low (e.g., Kolari in northern Finland and Kirov in Russia). Even the currently isolated Lillo was not much diverged. Further, there was no evidence of linkage disequilibrium in the data set, even between the very closely linked polymorphic sites. Linkage disequilibrium at neutral sites is governed by 1/(4Nec), and as c between closely linked sites is very small, Ne must be large to account for the lack of disequilibrium.
We also examined the data for signs of departure from demographic equilibrium, i.e., of postglacial expansion from several refugia, as explained earlier. A recent expansion could show up as an excess of singletons, which can be measured by Tajima's D. The value of -0.56 for the nuclear loci did not provide evidence for recent expansion. However, there are still local traces of possible refugia. The small Spanish populations may be remnants of a larger distribution. The observed lower level of DNA polymorphism in the Spanish population than in the Russian and Finnish populations agree with the previous data obtained with allozymes. The Spanish populations of Scots pine are slightly less variable than the populations from the continuous range of the distribution, with He of 0.325 and 0.363, respectively (Prus-Glowacki and Stephan 1994
). In all, the data on the Scots pine main range are consistent with a large effective population size and mutation drift equilibrium.
High Mutation Rate in Pines?
Because P. sylvestris does not have higher nucleotide diversity compared with other species, we need to examine again the hypothesis of a high mutation rate in pines. As the mutation rate governs the rate of neutral substitution between species, the synonymous divergence between species (Ks) provides an estimate of the mutation rate. Thus, µ = Ks/2T, where T is the time since divergence.
The best estimate of divergence time, based on fossil evidence, for P. sylvestris and P. taeda is 120 Myr (Millar 1998
), and the rate of synonymous substitutions is 0.15 x 10-9 per year for pal1 and about 0.05 x 10-9 per year for Adh (table 6
). The pal1 and the Adh genes have similar rates of neutral nucleotide substitutions, and both appear to evolve much more slowly in pines than in the monocots or dicots. Lu, Szmidt, and Wang (1998)
have observed a low number of substitutions between Gingko biloba and Larix sibirica at the nuclear coxI locus (based on their data, Ks = 0.016). This Ks estimate suggests an even lower mutation rate in those lineages than in pines because Gingko and Larix have diverged before the pines (about 200 MYA). A similar low rate (0.24 x 10-9) has also been found in the chloroplast gene rbcL of the chestnut tree (family Fagaceae), another group of long-lived woody perennials (Frascaria et al. 1993
). As shown in table 6
, the estimates of angiosperm mutation rates are higher by an order of magnitude. Note that even if we use a generation time of 25 years for pine, the per generation mutation rate estimates are not higher than in other plants in either pal or Adh. Lande, Schemske, and Schultz (1994)
have suggested that mutation rates in pines must be much higher than those measured in annual plants to account for the levels of inbreeding depression that pines maintain. It seems that the observed high inbreeding depression in Scots pine fits the model by Lande, Schemske, and Schultz (1994)
quite well (Koelewijn, Koski, and Savolainen 1999
). Mutation rates per generation at nuclear genes leading to chlorophyll deficiency have also been estimated to be about 10 times as high in pines and other woody plants as in annuals (Kärkkäinen, Koski, and Savolainen 1996; Klekowski 1998
). Thus, the genome-wide per generation mutation rate to deleterious alleles gives a very different picture from the neutral mutation rate per site per year. This difference merits further study.
|
Table 5 shows that Dr. melanogaster, like Scots pine, also has high enzyme variation compared with other species, but low nucleotide diversity compared with Dr. simulans. Enough comparisons are not yet available to account for these differences. We cannot exclude the possibility that the high variability at allozyme loci is partly maintained by selection, and thus not directly comparable with the patterns at presumably neutral sites.
Are Synonymous Sites Neutral: Codon Usage in P. sylvestris
Mutations occurring in the third positions of codons should be selectively neutral because they do not cause amino acid changes (Kimura 1968, 1983
). However, unequal use of synonymous codons points to selection at these sites also. Selection may favor some synonymous codons over the others, especially in highly expressed genes. There is a negative correlation between synonymous substitution rate and codon usage bias in many species (e.g., Sharp and Li 1986, 1989
). Selection for alternative codons is expected to be rather weak, and thus efficient only in large populations. Drosophila simulans has high nucleotide diversity and codon bias, and Dr. melanogaster has lower diversity and less codon bias (see table 5
). This is presumably because of the larger effective population size of Dr. simulans (Akashi 1997
). A further point of evidence for the large effective population size of Scots pine could be found by examining codon bias. As Scots pine has very large populations, even weak selection could be efficient. Table 5
shows a comparison of codon bias in Scots pine and P. radiata, a species with very limited population size and less genetic diversity than Scots pine (Moran, Bell, and Elridge 1988
). There is no evidence for increased codon bias in Scots pine relative to P. radiata, and thus at least not for the kind of selection at third positions because of major codon bias in large populations.
Selection at Linked Loci
Selection at linked loci has proved to have a major influence on patterns of diversity, especially in areas of low recombination (Aguadé, Miyashita, and Langley 1989
). This could be because of selective sweeps (Kaplan, Hudson, and Langley 1989
), background selection (Charlesworth, Morgan, and Charlesworth 1993
), or alternating selection pressures (pseudohitchhiking) (Gillespie 2000
). These selective explanations require that there be linkage disequilibrium between the neutral site and the target of selection, and then make additional assumptions. We do not yet know anything about the distribution of the rate of recombination in the P. sylvestris genome. Scots pine has high levels of deleterious mutations (as required by background selection), selection at many loci is efficient in large populations (as required by hitchhiking), and the direction of selection at these many loci can vary (as assumed in the pseudohitchhiking model). However, the similarity of the nucleotide diversity at pal1 gene to the 11 other loci and the low level of disequilibrium within the gene suggest that disequilibrium will also be low between the locus itself and adjacent areas. Thus, selection at linked loci is not the main cause for the low diversity. More genomic areas need to be studied.
![]() |
Acknowledgements |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Footnotes |
---|
Present address: Osteoporosis Research Center, Omaha, Nebraska, USA.
Keywords: nucleotide polymorphism
SNP
Pinus sylvestris
codon bias
linkage disequilibrium
Address for correspondence and reprints: Outi Savolainen, Department of Biology, University of Oulu, P.O. Box 3000, FIN-90014 Oulu, Finland. outi.savolainen{at}oulu.fi
.
![]() |
References |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Abbot R. J., M. F. Gomes, 1989 Population genetic structure and outcrossing rate of Arabidopsis thaliana (L.) Heynh Heredity 62:411-418[ISI]
Aguadé M., N. T. Miyashita, C. H. Langley, 1989 Reduced variation in the yellow-achaete-acute region in natural populations of Drosophila melanogaster Genetics 122:607-616
Akashi H., 1997 Codon bias evolution in Drosophila. Population genetics of mutation-selection drift Gene 205:269-278[ISI][Medline]
Butland S. L., M. L. Chow, B. E. Ellis, 1998 A diverse family of phenylalanine ammonia-lyase genes expressed in pine trees and cell cultures Plant Mol. Biol 37:15-24[ISI][Medline]
Cargill M., D. Altshuler, J. Ireland, et al. (17 co-authors) 1999 Characterization of single-nucleotide polymorphisms in coding regions of human genes Nat. Genet 22:231-238[ISI][Medline]
Charlesworth B., M. T. Morgan, D. Charlesworth, 1993 The effect of deleterious mutations on neutral molecular variation Genetics 134:1289-1303
Charlesworth D., Z. Yang, 1998 Allozyme diversity in Leavenworthia populations with different inbreeding levels Heredity 81:453-461[ISI][Medline]
Doebley J. F., M. M. Goodman, 1984 Isoenzymatic variation in Zea (Gramineae) Syst. Bot 9:203-218[ISI]
Doyle J. J., J. L. Doyle, 1990 Isolation of plant DNA from fresh tissue BRL Focus 12:13-15
Excoffier L., P. E. Smouse, J. M. Quattro, 1992 Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data Genetics 131:479-491
Eyre-Walker A., R. B. Gaut, H. Hilton, D. L. Feldman, B. S. Gaut, 1998 Investigation of the bottleneck leading to the domestication of maize Proc. Natl. Acad. Sci. USA 95:4441-4446
Filatov D. A., D. Charlesworth, 1999 DNA polymorphism, haplotype structure and balancing selection in the Leavenworthia PgiC locus Genetics 153:1423-1434
Frascaria N., L. Maggia, M. Michaud, J. Bousquet, 1993 The rbcL gene sequence from chestnut indicates a slow rate of evolution in the Fagaceae Genome 36:668-671[ISI][Medline]
Gaut B. S., M. T. Clegg, 1993 Nucleotide polymorphism in the Adh1 locus of pearl millet (Pennisetum glaucum) (Poaceae) Genetics 135:1091-1097
Gaut B. S., B. R. Morton, B. C. McCaig, M. T. Clegg, 1996 Substitution rate comparisons between grasses and palms: synonymous rate differences at the nuclear gene Adh parallel rate differences at the plastid gene rbcL Proc. Natl. Acad. Sci. USA 93:10274-10279
Gillespie J. H., 2000 Genetic drift in an infinite population: the pseudohitchhiking model Genetics 155:909-919
Goncharenko G. G., A. E. Silin, V. E. Padutov, 1994 Allozyme variation in natural populations of Eurasian pines: III. Population structure, diversity, differentiation and gene flow in central and isolated populations of Pinus sylvestris L. in Eastern Europe and Siberia Silvae Genet 43:119-132[ISI]
Graur D., W.-H. Li, 2000 Fundamentals of Molecular Evolution. Second edition Sinauer Associates, Inc., Sunderland
Gullberg U., R. Yazdani, D. Rudin, N. Ryman, 1985 Allozyme variation in Scotch pine (Pinus sylvestris) in Sweden Silvae Genet 34:193-201[ISI]
Hamrick J. L., M. J. W. Godt, 1996 Effects of life history traits on genetic diversity in plant species Phil. Trans. R. Soc. Lond. B. Biol. Sci 351:1291-1298[ISI]
Harpending H. C., 1994 Signature of ancient population growth in a low-resolution mitochondrial DNA mismatch distribution Hum. Biol 66:591-600[ISI][Medline]
Hill W. G., A. Robertson, 1968 Linkage disequilibrium in finite populations Theor. Appl. Genet 38:226-231
Hudson R. R., 1987 Estimating the recombination parameter of a finite population model without selection Genet. Res 50:245-250[ISI][Medline]
Huttley G. A., M. L. Durbin, D. E. Glover, M. T. Clegg, 1997 Nucleotide polymorphism in the chalcone synthase-A locus and evolution of the chalcone synthase multigene family of common morning glory Ipomoea purpurea Mol. Ecol 6:549-558[ISI]
Hyvärinen H., 1987 History of forests in northern Europe since the last glaciation. Ann. Acad. Sci. Fenn. Ser. A III Geol. Geogr 145:7-18
Innan H., F. Tajima, R. Terauchi, N. T. Miyashita, 1996 Intragenic recombination in the Adh locus of the wild plant Arabidopsis thaliana Genetics 143:1761-1770
Kaplan N. L., R. R. Hudson, C. H. Langley, 1989 The "hitchhiking effect" revisited Genetics 123:887-899
Karhu A., P. Hurme, M. Karjalainen, P. Karvonen, K. Kärkkäinen, D. Neale, O. Savolainen, 1996 Do molecular markers reflect patterns of differentiation in adaptive traits of conifers? Theor. Appl. Genet 93:215-221[ISI]
Kärkkäinen K., V. Koski, D. Savolainen, 1996 Geographical variation in inbreeding depression of Scots pine Evolution 50:111-119[ISI]
Karvonen P., O. Savolainen, 1993 Variation and inheritance of ribosomal DNA in Pinus sylvestris L. (Scots pine) Heredity 71:614-622[ISI]
Kawabe A., H. Innan, R. Terauchi, N. T. Miyashita, 1997 Nucleotide polymorphism in the acidic chitinase locus (ChiA) region of the wild plant Arabidopsis thaliana Mol. Biol. Evol 14:1303-1315[Abstract]
Kawabe A., N. T. Miyashita, 1999 DNA variation in the basic chitinase locus (ChiB) region of the wild plant Arabidopsis thaliana Genetics 153:1445-1453
Kawabe A., N. T. Miyashita, R. Terauchi, 1997 Phylogenetic relationship among the section Stenophora in the genus Dioscorea based on the analysis of nucleotide sequence variation in the phosphoglucose isomerase (Pgi) locus Genes Genet. Syst 72:253-262[ISI]
Kimura M., 1968 Genetic variability maintained in a finite population due to mutational production of neutral and nearly neutral isoalleles Genet. Res 11:247-269[ISI][Medline]
Kimura M., 1983 The neutral theory of molecular evolution Cambridge University Press, Cambridge, U.K.
Klekowski E., 1998 Mutation rates in mangroves and other plants Genetica 102/103:325-331
Koch M. A., B. Haubold, T. Mitchell-Olds, 2000 Comparative evolutionary analysis of chalcone synthase and alcohol dehydrogenase loci in Arabidopsis, Arabis, and related genera (Brassicaceae) Mol. Biol. Evol 17:1483-1498
Koelewijn H. P., V. Koski, O. Savolainen, 1999 Magnitude and timing of inbreeding depression in Scots pine (Pinus sylvestris L.) Evolution 53:758-768[ISI]
Koski V., 1970 A study of pollen dispersal a mechanism of gene flow in conifers Commun. Inst. For. Fenn 70:1-78
Koski V., 1971 Embryonic lethals of Picea abies and Pinus sylvestris Commun. Inst. For. Fenn 75:1-30
Lande R., D. W. Schemske, S. T. Schultz, 1994 High inbreeding depression, selective interference among loci, and the threshold selfing rate for purging recessive lethal mutations Evolution 48:965-978[ISI]
Lange B. M., M. Trost, W. Heller, C. Langebartels, H. Sandermann Jr., 1994 Elicitor-induced formation of free and cell-wall-bound stilbenes in cell-suspension cultures of Scots pine (Pinus sylvestris L.) Planta 194:143-148[ISI]
Laukkanen H., T. Sarjala, 1997 Effect of exogenous polyamines on Scots pine callus in vitro J. Plant Physiol 150:167-172[ISI]
Liu F., D. Charlesworth, M. Kreitman, 1999 The effect of mating system differences on nucleotide diversity at the phosphoglucose isomerase locus in the plant genus Leavenworthia Genetics 151:343-357
Liu F., L. Zhang, D. Charlesworth, 1998 Genetic diversity in Leavenworthia populations with different inbreeding levels Proc. R. Soc. Lond. Biol. Sci 265:293-301[ISI][Medline]
Lu M. Z., A. E. Szmidt, X.-R. Wang, 1998 RNA editing in gymnosperms and its impact on the evolution of the mitochondrial coxI gene Plant Mol. Biol 37:225-234[ISI][Medline]
Maynard Smith J., J. Haigh, 1974 The hitch-hiking effect of a favourable gene Genet. Res 23:23-35[ISI][Medline]
Millar C. I., 1998 Early evolution of pines Pp. 6991 in D. M. Richardson, ed. Ecology and biogeography of Pinus. Cambridge University Press, Cambridge, U.K
Miyashita N. T., A. Kawabe, H. Innan, R. Terauchi, 1998 Intra- and interspecific DNA variation and codon bias of the alcohol dehydrogenase (Adh) locus in Arabis and Arabidopsis species Mol. Biol. Evol 15:1420-1429
Moniz de Sá M., G. Drouin, 1996 Phylogeny and substitution rates in angiosperm actin genes Mol. Biol. Evol 13:1198-1212[Abstract]
Moran G. F., J. C. Bell, K. Elridge, 1988 The genetic structure and conservation of five natural populations of Pinus radiata Can. J. For. Res 18:506-514[ISI]
Moriyama E. N., J. R. Powell, 1996 Intraspecific nuclear DNA variation in Drosophila Mol. Biol. Evol 13:261-277[Abstract]
Müller-Starck G., P. H. Baradat, F. Bergmann, 1992 Genetic variation within European tree species Pp. 2347 in W. T. Adams, S. H. Straus, D. L. Copes, and A. R. Griffin, eds. Population genetics of forest trees. Kluwer Academic Publishers, Dordrecht, The Netherlands
Muona O., A. Harju, 1989 Effective population sizes, genetic variability, and mating system in natural stands and seed orchards of Pinus sylvestris Silvae Genet 38:221-228[ISI]
Muona O., A. Harju, K. Kärkkäinen, 1988 Background pollination in Pinus sylvestris seed orchards Scand. J. For. Res 4:513-520
Muona O., A. E. Szmidt, 1985 A multilocus study of natural populations of Pinus sylvestris Pp. 226240 in H. R. Gregorius, ed. Population genetics in forestry. Lecture notes in biomathematics, Vol. 60. Springer Verlag, Berlin
Nakamura Y., T. Gojobori, T. Ikemura, 2000 Codon usage tabulated from the international DNA sequence databases: status for the year 2000 Nucleic Acids Res 28:292
Nei M., 1987 Molecular evolutionary genetics Columbia University Press, New York
Nei M., W.-H. Li, 1979 Mathematical model for studying genetic variation in terms of restriction endonucleases Proc. Natl. Acad. Sci. USA 76:5269-5273[Abstract]
Nevo E., 1978 Genetic variation in natural populations: patterns and theory Theor. Popul. Biol 13:121-177[ISI][Medline]
Perry D. J., G. R. Furnier, 1996 Pinus banksiana has at least seven expressed alcohol dehydrogenase genes in two linked groups Proc. Natl. Acad. Sci. USA 93:13020-13023
Prus-Glowacki W., B. R. Stephan, 1994 Genetic variation of Pinus sylvestris from Spain in relation to other European populations Silvae Genet 43:7-14[ISI]
Remington D. L., R. W. Whetten, B. H. Liu, D. M. O' Malley, 1999 Construction of an AFLP genetic map with nearly complete genome coverage in Pinus taeda Theor. Appl. Genet 98:1279-1292[ISI][Medline]
Rosemann D., W. Heller, H. Sandermann Jr., 1991 Biochemical plant responses to ozone: II. Induction of stilbene biosynthesis in Scots pine (Pinus sylvestris L.) seedlings Plant Physiol 97:1280-1286[ISI]
Rozas J., Rozas R., 1999 DnaSP version 3: an integrated program for molecular population genetics and molecular evolution analysis Bioinformatics 15:174-175
Savolainen O., C. H. Langley, B. P. Lazzaro, H. Freville, 2000 Contrasting patterns of nucleotide polymorphism at the alcohol dehydrogenase locus in the outcrossing Arabidopsis lyrata and the selfing Arabidopsis thaliana Mol. Biol. Evol 17:645-655
Schneider S., D. Roessli, L. Excoffier, 2000 Arlequin Version 2.000: a software for population genetics analysis
Sewell M. M., B. K. Sherman, D. B. Neale, 1999 A consensus map for loblolly pine (Pinus taeda L.) I. Construction and integration of individual linkage maps from two outbred three-generation pedigrees Genetics 151:321-330
Sharp P. M., W.-H. Li, 1986 An evolutionary perspective on synonymous codon usage in unicellular organisms J. Mol. Evol 24:28-38[ISI][Medline]
Sharp P. M., W.-H. Li, 1989 On the rate of DNA sequence evolution in Drosophila J. Mol. Evol 28:398-402[ISI][Medline]
Sinclair W. T., J. D. Morman, R. A. Ennos, 1999 The postglacial history of Scots pine (Pinus sylvestris L.) in Western Europe: evidence from mitochondrial DNA variation Mol. Ecol 8:83-88[ISI]
Small R. L., J. A. Ryburn, J. F. Wendel, 1999 Low levels of nucleotide diversity at homologous Adh loci in allotetraploid cotton (Gossypium L.) Mol. Biol. Evol 16:491-501[Abstract]
Soranzo N., R. Alia, J. Provan, W. Powell, 2000 Patterns of variation at a mitochondrial sequence-tagged-site locus provides new insights into the postglacial history of European Pinus sylvestris populations Mol. Ecol 9:1205-1211[ISI][Medline]
Stephan W., C. H. Langley, 1998 DNA polymorphism in Lycopersicon and crossing-over per physical length Genetics 150:1585-1593
Terauchi R., T. Terachi, N. T. Miyashita, 1997 DNA polymorphism in the Pgi locus of a wild yam, Dioscorea tokoro Genetics 147:1899-1914
von Treuren R., H. Kuittinen, K. Kärkkäinen, E. Baena-Gonzalez, O. Savolainen, 1997 Evolution of microsatellites in Arabis petrea and Arabis lyrata, outcrossing relatives of Arabidopsis thaliana Mol. Biol. Evol 14:220-229[Abstract]
Wang X. R., A. E. Szmidt, D. Lindgren, 1991 Allozyme differentiation among populations of Pinus sylvestris (L.) from Sweden and China Hereditas 114:219-226[ISI]
Watterson G. A., 1975 On the number of segregating sites in genetical models without recombination Theor. Popul. Biol 7:256-276[ISI][Medline]
Whetten R. W., R. R. Sederoff, 1992 Phenylalanine ammonia-lyase from loblolly pine: purification of the enzyme and isolation of complementary DNA clones Plant Physiol 98:380-386[ISI]
Wright F., 1990 The "effective number of codons" used in a gene Gene 87:23-29[ISI][Medline]
Yamada T., Y. Tanaka, P. Sriprasertsak, H. Kato, T. Hashimoto, S. Kawamata, Y. Ichinose, H. Kato, T. Shiraishi, H. Oku, 1992 Phenylalanine ammonia-lyase genes from Pisum sativum: structure, organ-specific expression and regulation by fungal elicitor and suppressor Plant Cell Physiol 33:715-725[ISI]