Analysis of microsatellite variation in Pinus radiata reveals effects of genetic drift but no recent bottlenecks

Most conifer species occur in large continuous populations, but radiata pine, Pinus radiata, occurs only in five disjunctive natural populations in California and Mexico. The Mexican island populations were presumably colonized from the mainland millions of years ago. According to Axelrod (1981), the mainland populations are relicts of an earlier much wider distribution, reduced some 8,000 years ago, whereas according to Millar (1997, 2000), the patchy metapopulation-like structure is typical of the long-term population demography of the species. We used 19 highly polymorphic microsatellite loci to describe population structure and to search for signs of the dynamics of population demography over space and time. Frequencies of null alleles at microsatellite loci were estimated using an approach based on the probability of identity by descent. Microsatellite genetic diversities were high in all populations [expected heterozygosity (H(e)) = 0.68-0.77], but the island populations had significantly lower estimates. Variation between loci in genetic differentiation (F(ST)) was high, but no locus deviated statistically significantly from the rest at an experiment wide level of 0.05. Thus, all loci were included in subsequent analysis. The average differentiation was measured as F(ST) = 0.14 (SD 0.012), comparable with earlier allozyme results. The island populations were more diverged from the other populations and from an inferred common ancestral gene pool than the mainland ones. All populations showed a deficiency of expected heterozygosity given the number of alleles, the mainland populations more so than the island ones. The results thus do not support a recent important contraction in the mainland range of radiata pine.     

Phylogeography of the Korean pine (Pinus koraiensis) in northeast Asia: inferences from organelle gene sequences

Range-wide genetic variation of Korean pine (Pinus koraiensis) was assessed using maternally inherited mtDNA and paternally inherited cpDNA for 16 natural populations throughout northeast Asia in order to study its phylogeographical history during the Quaternary. The cpDNA variation indicated that there was no difference between populations on the Asian continent and those in the Japanese archipelago. In contrast, the mtDNA variation indicated that there was significant difference between the populations from the two regions, with each region having a different lineage. The continental populations exhibited no diversity in the mtDNA examined despite the species’ current extensive range and large populations. Conversely, while the Korean pine is rare in Japan, the Japanese populations exhibited greater levels of mtDNA diversity (H _T?=?0.502). The higher mtDNA diversity and evidence from numerous Korean pine macrofossil remains dated to the Pleistocene and recovered various sites in Japan suggest that the Japanese archipelago once served as a refugium to a much larger Korean pine population with a more extensive range than is the case today. The presence of the single mtDNA haplotype across the Asian continent suggests that the present widespread populations could have expanded from a single refugium population after the last glacial periods.     

Geographic pattern of genetic variation in Pinus resinosa: area of greatest diversity is not the origin of postglacial populations

Genetic diversity is low in natural populations of red pine, Pinus resinosa, a species that has a vast range across north-eastern North America. In this study, we examined 10 chloroplast microsatellite or simple sequence repeats (cpSSR) loci in 136 individuals from 10 widespread populations. Substantial variation for the cpSSR loci was observed in the study populations. The contrast with red pine's lack of variation for other types of loci is likely to be due to the higher mutation rates typical of SSR loci. The amount of variation is lower than that generally found for cpSSR loci in other pine species. In addition, the variation exhibits a striking geographical pattern. Most of the genetic diversity is among populations, with little within populations, indicating substantial isolation of and genetic drift within many populations in the southern half of the species distribution. The greatest diversity now occurs in the north-eastern part of New England, which is especially intriguing because this entire area was glaciated. Thus the centre of diversity cannot be the origin of postglacial populations, rather it is likely caused by admixture, most probably because of influences from two separate refugia. Furthermore, the pattern indicates that the spread of red pine since the last glaciation is rather more complex than usually described, and it likely includes more than one refugia, complex migration routes, and postglacial-retreat isolation and genetic drift among shrinking populations in regions of the present southern range.     

Glacial vicariance in Eurasia: mitochondrial DNA evidence from Scots pine for a complex heritage involving genetically distinct refugia at mid-northern latitudes and in Asia Minor

Background  

At the last glacial maximum, Fennoscandia was covered by an ice sheet while the tundra occupied most of the rest of northern Eurasia. More or less disjunct refugial populations of plants were dispersed in southern Europe, often trapped between mountain ranges and seas. Genetic and paleobotanical evidences indicate that these populations have contributed much to Holocene recolonization of more northern latitudes. Less supportive evidence has been found for the existence of glacial populations located closer to the ice margin. Scots pine (Pinus sylvestris L.) is a nordic conifer with a wide natural range covering much of Eurasia. Fractures in its extant genetic structure might be indicative of glacial vicariance and how different refugia contributed to the current distribution at the continental level. The population structure of Scots pine was investigated on much of its Eurasian natural range using maternally inherited mitochondrial DNA polymorphisms.     

Microsatellite analysis reveals genetically distinct populations of red pine (Pinus resinosa, Pinaceae)

Red pine (Pinus resinosa Ait.) is an ecologically and economically important forest tree species of northeastern North America and is considered one of the most genetically depauperate conifer species in the region. We have isolated and characterized 13 nuclear microsatellite loci by screening a partial genomic library with di-, tri-, and tetranucleotide repeat oligonucleotide probes. In an analysis of over 500 individuals representing 17 red pine populations from Manitoba through Newfoundland, five polymorphic microsatellite loci with an average of nine alleles per locus were identified. The mean expected and observed heterozygosity values were 0.508 and 0.185, respectively. Significant departures from Hardy-Weinberg equilibrium with excess homozygosity indicating high levels of inbreeding were evident in all populations studied. The population differentiation was high with 28-35% of genetic variation partitioned among populations. The genetic distance analysis showed that three northeastern (two Newfoundland and one New Brunswick) populations are genetically distinct from the remaining populations. The coalescence-based analysis suggests that "northeastern" and "main" populations likely became isolated during the most recent Pleistocene glacial period, and severe population bottlenecks may have led to the evolution of a highly selfing mating system in red pine.     

Regional patterns of genetic diversity in Pinus flexilis (Pinaceae) reveal complex species history

Pinus flexilis (limber pine) is patchily distributed within its large geographic range; it is mainly restricted to high elevations in the Rocky Mountains and the Basin and Range region of western North America. We examined patterns of allozyme diversity in 30 populations from throughout the species' range. Overall genetic diversity (H(e) = 0.186) was high compared with that of most other pine species but was similar to that of other pines widespread in western North America. The proportion of genetic diversity occurring among populations (G(ST) = 0.101) was also high relative to that for other pines. Observed heterozygosity was less than expected in 28 of the 30 populations. When populations were grouped by region, there were notable differences. Those in the Basin and Range region had more genetic diversity within populations, a higher proportion of genetic diversity among populations, and higher levels of inbreeding within populations than populations from either the Northern or Utah Rocky Mountain regions. Patterns of genetic diversity in P. flexilis have likely resulted from a complex distribution of Pleistocene populations and subsequent gene flow via pollen and seed dispersal.     

Landscape scale genetic effects of habitat fragmentation on a high gene flow species: Speyeria idalia (Nymphalidae)

Detection of the genetic effects of recent habitat fragmentation in natural populations can be a difficult task, especially for high gene flow species. Previous analyses of mitochondrial DNA data from across the current range of Speyeria idalia indicated that the species exhibited high levels of gene flow among populations, with the exception of an isolated population in the eastern portion of its range. However, some populations are found on isolated habitat patches, which were recently separated from one another by large expanses of uninhabitable terrain, in the form of row crop agriculture. The goal of this study was to compare levels of genetic differentiation and diversity among populations found in relatively continuous habitat to populations in both recently and historically isolated habitat. Four microsatellite loci were used to genotype over 300 individuals from five populations in continuous habitat, five populations in recently fragmented habitat, and one historically isolated population. Results from the historically isolated population were concordant with previous analyses and suggest significant differentiation. Also, microsatellite data were consistent with the genetic effects of habitat fragmentation for the recently isolated populations, in the form of increased differentiation and decreased genetic diversity when compared to nonfragmented populations. These results suggest that given the appropriate control populations, microsatellite markers can be used to detect the effects of recent habitat fragmentation in natural populations, even at a large geographical scale in high gene flow species.     

High genetic diversity and distinct origin of recently fragmented Scots pine (<Emphasis Type="Italic">Pinus sylvestris</Emphasis> L.) populations along the Carpathians and the Pannonian Basin

Historical evolutionary events highly affect the modern-day genetic structure of natural populations. Scots pine (Pinus sylvestris L.), as a dominant tree species of the Eurasian taiga communities following the glacial cycles of the Pleistocene, has survived in small, scattered populations at the range limits of its south-eastern European distribution. In this study, we examined genetic relationships, genetic divergence and demographic history of peripheral populations from central-eastern Europe, the Carpathian Mountains and the Pannonian Basin. Four hundred twenty-one individuals from 20 populations were sampled and characterized with both nuclear and chloroplast simple sequence repeat (SSR) markers. Standard population genetic indices, the degree of genetic differentiation and spatial genetic structure were analysed. Our results revealed that peripheral Scots pine populations retained high genetic diversity despite the recently ongoing fragmentation and isolation of the persisting relict populations. Analysis of molecular variance (AMOVA) showed 7% among-population genetic differentiation, and there was no isolation by distance among the island-like occurrences. Genetic discontinuities with strong barriers (99–100% bootstrap support) were identified in the Carpathians. Based on both marker types, populations of the Western Carpathians were delimited from those inhabiting the Eastern Carpathians, and two main genetic lineages were traced that most probably originate from two main refugia. One refugium presumably existed in the region of the Eastern Alps with the Hungarian Plain, while the other was probably found in the Eastern Carpathians. These findings are supported by recent palynological records. The strongest genetic structure was revealed within the Romanian Carpathians on the basis of both marker types. With only some exceptions, no signs of recent bottlenecks or inbreeding were detected. However, Carpathian natural populations of Scots pine are highly fragmented and have a small census size, though they have not yet been affected by genetic erosion induced by isolation.     

Limited emigration from an outbreak of a forest pest insect

Population density and individual dispersal behaviour affect species' distribution dynamics. Population densities vary over time, and some species occasionally increase to very high numbers, for example during outbreaks. In such situations, populations are expected to expand into new areas as a result of density-dependent dispersal which sometimes even results in range expansion. A local population of the northern pine processionary moth Thaumetopoea pinivora has recently reached outbreak densities at the edge of its northern range at the southern tip of Gotland Island in the Baltic Sea. We first investigated whether the outbreak had resulted in establishment of populations in suitable habitats on Gotland Island outside the outbreak area. Six small populations were found that could potentially have originated from the outbreak area. However, data from 12 microsatellite markers strongly suggest that these populations did not originate from the recent outbreak. Genetic variability was not reduced in these small, isolated populations, and there were several unique alleles, indicating instead a different population history and that there has been no recent range expansion. In addition, there was apparent genetic isolation by geographic distance, implying that despite the high density of the outbreak population, significant gene flow has not occurred.     

High‐density SNP assay development for genetic analysis in maritime pine (Pinus pinaster)

Maritime pine provides essential ecosystem services in the south‐western Mediterranean basin, where it covers around 4 million ha. Its scattered distribution over a range of environmental conditions makes it an ideal forest tree species for studies of local adaptation and evolutionary responses to climatic change. Highly multiplexed single nucleotide polymorphism (SNP) genotyping arrays are increasingly used to study genetic variation in living organisms and for practical applications in plant and animal breeding and genetic resource conservation. We developed a 9k Illumina Infinium SNP array and genotyped maritime pine trees from (i) a three‐generation inbred (F2) pedigree, (ii) the French breeding population and (iii) natural populations from Portugal and the French Atlantic coast. A large proportion of the exploitable SNPs (2052/8410, i.e. 24.4%) segregated in the mapping population and could be mapped, providing the densest ever gene‐based linkage map for this species. Based on 5016 SNPs, natural and breeding populations from the French gene pool exhibited similar level of genetic diversity. Population genetics and structure analyses based on 3981 SNP markers common to the Portuguese and French gene pools revealed high levels of differentiation, leading to the identification of a set of highly differentiated SNPs that could be used for seed provenance certification. Finally, we discuss how the validated SNPs could facilitate the identification of ecologically and economically relevant genes in this species, improving our understanding of the demography and selective forces shaping its natural genetic diversity, and providing support for new breeding strategies.     

Chloroplast DNA Transgresses Species Boundaries and Evolves at Variable Rates in the California Closed-Cone Pines (Pinus radiata, P. muricata, and P. attenuata)

We studied phylogenetic relationships among populations and species in the California closed-cone pines (Pinus radiata D. Don, P. attenuata Lemm., and P. muricata D. Don) via chloroplast DNA restriction site analysis. Data on genetic polymorphism within and among 19 populations in the three species were collected using9 to 20 restriction enzymes and 38 to 384 trees. Because only five clades and extremely low intraclade diversity were found, additional phylogenetic data were collected using a single representative per clade and two outgroup species, P. oocarpa Schiede and P. jeffreyi Loud. In total, 25 restriction enzymes were employed and approximately 2.7 kb surveyed (2.3% of genome). The five clades recognized were Monterey pine, knob-cone pine, and the southern, intermediate, and northern races of bishop pine. On the basis of bootstrapping, both Wagner and Dollo parsimony analyses strongly separated the northern and intermediate races of bishop pine from the southern race; knobcone pine from Monterey and bishop pines; and the closed-cone pines from the two outgroups. Approximate divergence times were estimated for the lineages leading to knob-cone pine and to the intermediate and northern populations of bishop pine. The position of Monterey pine relative to bishop pine within their monophyletic clade was unresolved. Surprisingly, Montery pine and the southern race of bishop pine were much more similar to one another than was the southern race of bishop pine to its conspecific intermediate and northern races. Both the Monterey and southern bishop pine lineages also evolved severalfold more slowly than did the knobcone pine and intermediate-northern bishop pine lineages. These results differ significantly from a recent allozyme study, corroborating previous observations that chloroplast genome phylogeny can depart substantially from that of nuclear genes.     

Genetic differentiation and reduced genetic diversity at the northern range edge of two species with different dispersal modes

Theory predicts that genetic variation should be reduced at range margins, but empirical support is equivocal. Here, we used genotyping‐by‐sequencing technology to investigate genetic variation in central and marginal populations of two species in the marine gastropod genus Crepidula. These two species have different development and dispersal types and might therefore show different spatial patterns of genetic variation. Both allelic richness and the proportion of private alleles were highest in the most central populations of both species, and lower at the margin. The species with low dispersal, Crepidula convexa, showed high degrees of structure throughout the range that conform to the pattern found in previous studies using other molecular markers. The northernmost populations of the high‐dispersing species, Crepidula fornicata, are distinct from more central populations, although this species has been previously observed to have little genetic structure over much of its range. Although genetic diversity was significantly lower at the range margin, the absolute reduction in diversity observed with these genomewide markers was slight, and it is not yet known whether there are functional consequences for the marginal populations.     

Low nucleotide diversity at the pal1 locus in the widely distributed Pinus sylvestris.

Nucleotide polymorphism in Scots pine (Pinus sylvestris) was studied in the gene encoding phenylalanine ammonia-lyase (Pal, EC 4.3.1.5). Scots pine, like many other pine species, has a large current population size. The observed levels of inbreeding depression suggest that Scots pine may have a high mutation rate to deleterious alleles. Many Scots pine markers such as isozymes, RFLPs, and microsatellites are highly variable. These observations suggest that the levels of nucleotide variation should be higher than those in other plant species. A 2,045-bp fragment of the pal1 locus was sequenced from five megagametophytes each from a different individual from each of four populations, from northern and southern Finland, central Russia, and northern Spain. There were 12 segregating sites in the locus. The synonymous site overall nucleotide diversity was only 0.0049. In order to compare pal1 with other pine genes, sequence was obtained from two alleles of 11 other loci (total length 4,606 bp). For these, the synonymous nucleotide diversity was 0.0056. These estimates are lower than those from other plants. This is most likely because of a low mutation rate, as estimated from between-pine species synonymous site divergence. In other respects, Scots pine has the characteristics of a species with a large effective population. There was no linkage disequilibrium even between closely linked sites. This resulted in high haplotype diversity (14 different haplotypes among 20 sequences). This could also give rise to high per locus diversity at the protein level. Divergence between populations in the main range was low, whereas an isolated Spanish population had slightly lower diversity and higher divergence than the remaining populations.     

Variation on reaction norm in lodgepole pine natural populations

Variations on the norm of reaction among ten natural lodgepole pine populations sampled from three lodgepole pine subspecies (Pinus contorta ssp. contorta, ssp. latifolia and ssp. murrayana) were studied by using 20 year heights measured in 57 provenance test sites across interior British Columbia (B.C.). There were significant population by site interactions. Concurrent joint regression and the AMMI model were used to dissect these population by environmental interactions. Joint regression analysis indicated that three populations (from the northwest) had a negative linear regression coefficient with environmental deviation, three (from central and southeast sites) had a positive regression coefficient and four (from the southwest) had a zero regression coefficient. The AMMI model revealed a similar pattern of reaction norm among the ten populations. But the three significant IPCA axes, which captured twice as much of the G × E sum of squares than joint regression, were more effective in separating the ten populations and associating their performance with the climate of test sites and their origin. The variation patterns of reaction norm in lodgepole pine populations demonstrated that adaptation of lodgepole pine natural populations to the various physical environments, at sub-species as well as at population level, was due largely to a balance between selection for high growth potential in less severe environments and selection for high cold hardiness in severe environments. Received: 4 May 2000 / Accepted: 10 November 2000     

Paternally Inherited Chloroplast Polymorphism in Pinus: Estimation of Diversity and Population Subdivision, and Tests of Disequilibrium with a Maternally Inherited Mitochondrial Polymorphism

We have surveyed a chloroplast DNA restriction fragment length polymorphism in 745 individuals, distributed rangewide in eight allopatric natural populations of jack pine (Pinus banksiana Lamb.) and eight allopatric natural populations of lodgepole pine (Pinus contorta Dougl.). The polymorphic region of the chloroplast genome is located near duplicated psbA genes. Fourteen length variants were found in the survey, and these variants distinguished the two species qualitatively. Variant diversities were high in both species (h(es) = 0.43 in jack pine; h(es) = 0.44 in lodgepole pine). Population subdivision was weak within and among lodgepole pine subspecies and in jack pine (i.e., θvalues were less than 0.05). This weak subdivision is compatible with theoretical predictions for paternally inherited markers in wind-pollinated outcrossers, as well as for polymorphisms with high length mutation rates. If these populations are at a drift-migration equilibrium, the chloroplast DNA restriction fragment data and previous mitochondrial frequency data from the same individuals are consistent with gene flow that is differential through seeds and pollen. The new data have permitted the first empirical tests of disequilibrium between maternally and paternally inherited factors. As expected, these tests failed to detect convincing evidence of non-random association between chloroplast and mitochondrial variants.     

The critical role of local refugia in postglacial colonization of Chinese pine: joint inferences from DNA analyses,pollen records,and species distribution modeling

The importance of long‐distance migration from low to high latitudes relative to local spread from northern refugia after the Last Glacial Maximum (LGM) remains a focus of debate for many temperate tree species. We assessed the dynamics of Chinese pine Pinus tabulaeformis, a widespread species endemic to northern China, since the LGM by integrating cytoplasmic DNA data, mapped pollen records and ecological niche modeling. Genetic variation among 544 individuals from 50 populations spanning the entire natural species range revealed eight genetic clusters with distinct geographic distribution, indicating glacial lineages likely originating from multiple local microrefugia. Palynological evidence suggested that the northernmost part of the natural distribution originated from local postglacial spread. Niche modeling indicated high probability of the species being present in the area of the Loess Plateau and coastal areas north of the Yangtze River during the LGM. The three lines of evidence jointly suggest that the species persisted through the last glaciation in the mountains surrounding the Loess Plateau of northern China and that the current distribution of the species originated primarily from the spread of local refugial populations, instead of long‐distance migration. These results cast doubt on the notion that Chinese pine migrated from areas south of the Yangtze River and underscore the importance of northern refugia.     

Allozyme Variation in the Rare Gulf Coast Endemic Polygonella macrophylla Small (Polygonaceae)

Abstract Polygonella macrophylla Small is a rare, perennial, primarily gynodioecious plant species endemic to a narrow zone of coastal sand pine scrub habitat along the Gulf of Mexico from Carrabelle, Florida, to Gulf Shores, Alabama (U.S.A.). The species is comprised of a crimson-red flowered form ("rubra"), known from only two populations at the eastern distributional limit of the species, and a white/pink flowered form ("alba") represented by several populations throughout the remainder of the range of the species. An electrophoretic investigation revealed that commonly used measures of gene diversity are very close to expected values in "rubra" given the combination of ecological attributes displayed by the species; however, gene diversity in "alba" is much lower than expected. Population differentiation is high, with estimates of the coefficient of gene differentiation ( G _ST ) indicating that more than 30% of the diversity is among populations. Lack of gene flow among populations, estimated to be at most 0.6 migrating individuals per generation for the species, is a likely cause of greater population differentiation. High inbreeding, detected both within populations and within the species, is suggestive of high rates of self-fertilization in hermaphroditic plants and/or biparental inbreeding in small populations.     

Different patterns of genetic structure of relict and isolated populations of endangered peat-bog pine (Pinus uliginosa Neumann)

Recent changes in environmental conditions in populations of peat-bog pine (Pinus uliginosa Neumann) caused rapid decline or even extinction of the species in several stands in Central Europe. Conservation strategies forP. uliginosa require information about the evolutionary history and genetic structure of its populations. Using isozymes we assessed the genetic structure ofP. uliginosa from four isolated stands in Poland and compared the results to genetic structures of other closely related pine species including eight populations ofPinus mugo, ten ofPinus sylvestris and one ofPinus uncinata. The level of genetic variability ofP. uliginosa measured by the mean number of alleles per locus and average heterozygosity was similar to others related toP. uliginosa taxa from the reference group but it differs among populations. High genetic similarity was found between two populations ofP. uliginosa from Low Silesian Pinewood. The populations were genetically distinct as compared to other populations includinglocus classicus of the species from the peat bog at Batorów Reserve. Very low genetic distance (D_n = 0.002) and small genetic differentiation (GST = 0.003) were found betweenP. uliginosa andP. mugo in the sympatric populations of the species from Zieleniec peat bog suggesting the ongoing natural hybridisation and genetic contamination of peat-bog pine from this area. Some evidence for skew in allele frequency distribution potentially due to recent bottleneck was found in population from Low Silesian Pinewood. The analysed open pollinated progeny derived from twoP. uliginosa stands from Low Silesian Pine-wood showed the excess of homozygotes as compared to the maternal trees indicating high level of inbreeding (F = 0.105,F = 0.081). The results are discussed in the context of evolution ofP. uliginosa populations, taxonomic relationships between the analysed species and conservation strategies for active protection of peat-bog pine.     

Characteristics of gene pools of three Pinus pumila (Pall.) Regel populations at the range margins

Genetic variation of Siberian dwarf pine Pinus pumila (Pall.) Regel was characterized in three marginal populations in southwestern, southern and eastern parts of the natural species range (Transbaikalia, Primorye, Kamchatka) using isozyme analysis. Analysis involving 16 isozyme loci encoding ten enzyme systems was conducted. Our results confirm that P. pumila is among the most polymorphic species in the family Pinus. Three marginal populations exhibited high genetic variation (P95 = 68.8%, Ho = 0.247, He = 0.291). Populations heterogeneity and significantly high level of divergence in coniferous (F(ST) = 0.050, D(N) = 0.044) reflect their genetic originality. In summary, it was shown that the level of genetic variation characteristic for P. pumila in other parts of the not only is reproduced in the populations examined but even is close to maximum there.     

Geographic pattern of genetic diversity in Pinus resinosa: contact zone between descendants of glacial refugia

Although red pine (Pinus resinosa) generally has low or completely lacks variation for molecular markers, some variation is observed for chloroplast microsatellites (cpSSRs). We sampled and examined 10 cpSSRs for 19 populations. Analysis of these populations plus 10 previously studied populations shows that the geographic distribution of genetic diversity over the range of P. resinosa is markedly nonuniform. Although the pattern exhibits little isolation by distance, there is a region centered in northeastern New England where populations contain much greater chloroplast haplotype diversity than elsewhere. This area is band-shaped, with the longer axis nearly parallel with latitude, and very sharply delineated. The area of high diversity was buried by the Laurentide ice sheet. The geographic pattern indicates that P. resinosa is not at equilibrium, and the species has had a more complex postglacial history than typically purported for forest trees in eastern North America. The results suggest that the area of high diversity is a stable transition zone between descendants of two distinct refugia, one in the southern Appalachians and another near the North Atlantic coastline of the Wisconsinian glacial period. Plausible explanations are given that selection between two lineages, along latitudinal zones, may have maintained the transition zone.