DIRECT AND INDIRECT ESTIMATES OF SEED VERSUS POLLEN MOVEMENT WITHIN A POPULATION OF PONDEROSA PINE

Capabilities and limitations of using DNA metabarcoding to study plant-pollinator interactions

Many pollinator populations are experiencing declines, emphasizing the need for a better understanding of the complex relationship between bees and flowering plants. Using DNA metabarcoding to describe plant-pollinator interactions eliminates many challenges associated with traditional methods and has the potential to reveal a more comprehensive understanding of foraging behavior and pollinator life history. Here we use DNA metabarcoding of ITS2 and rbcL gene regions to identify plant species present in pollen loads of 404 bees from three habitats in eastern Oregon. Our specific objectives were to 1) determine whether plant species identified using DNA metabarcoding are consistent with plant species identified using observations, 2) compare characterizations of diet breadth derived from foraging observations to those based on plant species assignments obtained using DNA metabarcoding, and 3) compare plant species assignments produced by DNA metabarcoding using a "regional" reference database to those produced using a "local" database. At the three locations, 31-86% of foraging observations were consistent with DNA metabarcoding data, 8-50% of diet breadth characterizations based on observations differed from those based on DNA metabarcoding data, and 22-25% of plant species detected using the regional database were not known to occur in the study area in question. Plant-pollinator networks produced from DNA metabarcoding data had higher sampling completeness and significantly lower specialization than networks based on observations. Here, we examine some strengths and limitations of using DNA metabarcoding to identify plant species present in bee pollen loads, make ecological inferences about foraging behavior, and provide guidance for future research.

Population genetic structure and connectivity of a riparian selfing herb Caulokaempferia coenobialis at a fine-scale geographic level in subtropical monsoon forest

BACKGROUND

Rivers and streams facilitate movement of individuals and their genes across the landscape and are generally recognized as dispersal corridors for riparian plants. Nevertheless, some authors have reported directly contrasting results, which may be attributed to a complex mixture of factors, such as the mating system and dispersal mechanisms of propagules (seed and pollen), that make it difficult to predict the genetic diversity and population structure of riparian species. Here, we investigated a riparian self-fertilizing herb Caulokaempferia coenobialis, which does not use anemochory or zoochory for seed dispersal; such studies could contribute to an improved understanding of the effect of rivers or streams on population genetic diversity and structure in riparian plants. Using polymorphic ISSR and cpDNA loci, we studied the effect at a microgeographic scale of different stream systems (a linear stream, a dendritic stream, and complex transverse hydrological system) in subtropical monsoon forest on the genetic structure and connectivity of C. coenobialis populations across Dinghu Mountain (DH) and Nankun Mountain (NK).

RESULTS

The results indicate that the most recent haplotypes (DH: H7, H8; NK: h6, h7, h11, h12) are not shared among local populations of C. coenobialis within each stream system. Furthermore, downstream local populations do not accumulate genetic diversity, whether in the linear streamside local populations across DH (H: 0.091 vs 0.136) or the dendritic streamside local populations across NK (H: 0.079 vs 0.112, 0.110). Our results show that the connectivity of local C. coenobialis populations across DH and NK can be attributed to historical gene flows, resulting in a lack of spatial genetic structure, despite self-fertilization. Selfing C. coenobialis can maintain high genetic diversity (H = 0.251; I = 0.382) through genetic differentiation (GST = 0.5915; FST = 0.663), which is intensified by local adaptation and neutral mutation and/or genetic drift in local populations at a microgeographic scale.

CONCLUSION

We suggest that streams are not acting as corridors for dispersal of C. coenobialis, and conservation strategies for maintaining genetic diversity of selfing species should be focused on the protection of all habitat types, especially isolated fragments in ecosystem processes.

Genetic variation and genetic structure within metapopulations of two closely related selfing and outcrossing Zingiber species (Zingiberaceae)

Habitat fragmentation strongly affects the genetic diversity of plant populations, and this has always attracted much research interest. Although numerous studies have investigated the effects of habitat fragmentation on the genetic diversity of plant populations, fewer studies have compared species with contrasting breeding systems while accounting for phylogenetic distance. Here, we compare the levels of genetic diversity and differentiation within and among subpopulations in metapopulations (at fine-scale level) of two closely related Zingiber species, selfing Zingiber corallinum and outcrossing Zingiber nudicarpum. Comparisons of the genetic structure of species from unrelated taxa may be confounded by the effects of correlated ecological traits or/and phylogeny. Thus, we possibly reveal the differences in genetic diversity and spatial distribution of genetic variation within metapopulations that relate to mating systems. Compared to outcrossing Z. nudicarpum, the subpopulation genetic diversity in selfing Z. corallinum was significantly lower, but the metapopulation genetic diversity was not different. Most genetic variation resided among subpopulations in selfing Z. corallinum metapopulations, while a significant portion of variation resided either within or among subpopulations in outcrossing Z. nudicarpum, depending on whether the degree of subpopulation isolation surpasses the dispersal ability of pollen and seed. A stronger spatial genetic structure appeared within subpopulations of selfing Z. corallinum potentially due to restricted pollen flow and seed dispersal. In contrast, a weaker genetic structure was apparent in subpopulations of outcrossing Z. nudicarpum most likely caused by extensive pollen movement. Our study shows that high genetic variation can be maintained within metapopulations of selfing Zingiber species, due to increased genetic differentiation intensified primarily by the stochastic force of genetic drift among subpopulations. Therefore, maintenance of natural variability among subpopulations in fragmented areas is key to conserve the full range of genetic diversity of selfing Zingiber species. For outcrossing Zingiber species, maintenance of large populations is an important factor to enhance genetic diversity. Compared to outcrossing Z. nudicarpum, the subpopulation genetic diversity in selfing Z. corallinum was significantly lower, but the metapopulation genetic diversity did not differ. Most genetic variation resided among subpopulations in selfing Z. corallinum metapopulations, while a significant portion of variation resided either within or among subpopulations in outcrossing Z. nudicarpum, depending on whether the degree of subpopulation isolation surpasses the dispersal ability of pollen and seed. Our study shows that selfing Z. corallinum could maintain high genetic diversity through differentiation intensified primarily by the stochastic force of genetic drift among subpopulations at fine-scale level, but not local adaptation.

Landscape Genetics of Plants: Challenges and Opportunities

Dispersal is one of the most important but least understood processes in plant ecology and evolutionary biology. Dispersal of seeds maintains and establishes populations, and pollen and seed dispersal are responsible for gene flow within and among populations. Traditional views of dispersal and gene flow assume models that are governed solely by geographic distance and do not account for variation in dispersal vector behavior in response to heterogenous landscapes. Landscape genetics integrates population genetics with Geographic Information Systems (GIS) to evaluate the effects of landscape features on gene flow patterns (effective dispersal). Surprisingly, relatively few landscape genetic studies have been conducted on plants. Plants present advantages because their populations are stationary, allowing more reliable estimates of the effects of landscape features on effective dispersal rates. On the other hand, plant dispersal is intrinsically complex because it depends on the habitat preferences of the plant and its pollen and seed dispersal vectors. We discuss strategies to assess the separate contributions of pollen and seed movement to effective dispersal and to delineate the effects of plant habitat quality from those of landscape features that affect vector behavior. Preliminary analyses of seed dispersal for three species indicate that isolation by landscape resistance is a better predictor of the rates and patterns of dispersal than geographic distance. Rates of effective dispersal are lower in areas of high plant habitat quality, which may be due to the effects of the shape of the dispersal kernel or to movement behaviors of biotic vectors. Landscape genetic studies in plants have the potential to provide novel insights into the process of gene flow among populations and to improve our understanding of the behavior of biotic and abiotic dispersal vectors in response to heterogeneous landscapes.

Genetic diversity and relationship between cultivated, weedy and wild rye species as revealed by chloroplast and mitochondrial DNA non-coding regions analysis

The genus Secale is small but very diverse. Despite the high economic importance, phylogenetic relationships of rye species have not been fully determined, and they are extremely important for the process of breeding of new cultivars that can be enriched with functional traits derived from wild rye species. The study analyzed the degree of relationship of 35 accessions of the genus Secale, representing 13 most often distinguished species and subspecies, originating from various seed collections in the world, based on the analysis of non-coding regions of the chloroplast (cpDNA) and mitochondrial genome (mtDNA), widely used in phylogenetic and population plant studies, because of a higher rate of evolution than the coding regions. There was no clear genetic structure between different species and subspecies, which may indicated the introgression between these taxa. The obtained data confirmed that S. vavilovii was very similar to S. cereale, which confirmed the assumption that they might share a common ancestor. The results also confirmed the divergence of S. sylvestre from other species and subspecies of rye. Areas that may be useful molecular markers in studies on closely related species of the genus Secale were also indicated.

Genetic Diversity among Age Classes of a Pinus sylvestris (L.) Population from the Białowieża Primeval Forest, Poland

The Białowieża Primeval Forest (BPF) is Europe’s last primeval forest and an irreplaceable area for biodiversity conservation due to its size, protection status, and substantially undisturbed nature. There is no other forest in Europe with such a large surface representing highly-advanced natural succession. This article reports on the first analysis of the genetic variability and demographic structure of a self-renewed Pinus sylvestris population located in BPF, using both chloroplast and mitochondrial DNA markers. The analysis of molecular variance (AMOVA) for chloroplast simple sequence repeats (cpSSRs) revealed a significant genetic differentiation among age classes that accounted for about 2% of the total variance, comparable to those reported among different populations of Scots pine. None of the 117 detected chloroplast haplotypes were common to all age classes. Haplotype diversity ranged from 0.370 to 0.415 for cpSSRs and from 0.320 to 0.455 for mitochondrial markers. The genetic variation of the studied age classes—represented by mitochondrial markers—strongly depicts the maternal genetic structure, indicating limited seed dispersal. Temporal genetic substructuring is maintained within a self-renewed population of Scots pine from the BPF.

Exploring Climate Niches of Ponderosa Pine (Pinus ponderosa Douglas ex Lawson) Haplotypes in the Western United States: Implications for Evolutionary History and Conservation

Ponderosa pine (Pinus ponderosa Douglas ex Lawson) occupies montane environments throughout western North America, where it is both an ecologically and economically important tree species. A recent study using mitochondrial DNA analysis demonstrated substantial genetic variation among ponderosa pine populations in the western U.S., identifying 10 haplotypes with unique evolutionary lineages that generally correspond spatially with distributions of the Pacific (P. p. var. ponderosa) and Rocky Mountain (P. p. var. scopulorum) varieties. To elucidate the role of climate in shaping the phylogeographic history of ponderosa pine, we used nonparametric multiplicative regression to develop predictive climate niche models for two varieties and 10 haplotypes and to hindcast potential distribution of the varieties during the last glacial maximum (LGM), ~22,000 yr BP. Our climate niche models performed well for the varieties, but haplotype models were constrained in some cases by small datasets and unmeasured microclimate influences. The models suggest strong relationships between genetic lineages and climate. Particularly evident was the role of seasonal precipitation balance in most models, with winter- and summer-dominated precipitation regimes strongly associated with P. p. vars. ponderosa and scopulorum, respectively. Indeed, where present-day climate niches overlap between the varieties, introgression of two haplotypes also occurs along a steep clinal divide in western Montana. Reconstructed climate niches for the LGM suggest potentially suitable climate existed for the Pacific variety in the California Floristic province, the Great Basin, and Arizona highlands, while suitable climate for the Rocky Mountain variety may have existed across the southwestern interior highlands. These findings underscore potentially unique phylogeographic origins of modern ponderosa pine evolutionary lineages, including potential adaptations to Pleistocene climates associated with discrete temporary glacial refugia. Our predictive climate niche models may inform strategies for further genetic research (e.g., sampling design) and conservation that promotes haplotype compatibility with projected changes in future climate.

Mitochondrial DNA haplotype distribution patterns in Pinus ponderosa (Pinaceae): range-wide evolutionary history and implications for conservation

PREMISE OF THE STUDY

Ponderosa pine (Pinus ponderosa Douglas ex P. Lawson & C. Lawson) exhibits complicated patterns of morphological and genetic variation across its range in western North America. This study aims to clarify P. ponderosa evolutionary history and phylogeography using a highly polymorphic mitochondrial DNA marker, with results offering insights into how geographical and climatological processes drove the modern evolutionary structure of tree species in the region.

METHODS

We amplified the mtDNA nad1 second intron minisatellite region for 3,100 trees representing 104 populations, and sequenced all length variants. We estimated population-level haplotypic diversity and determined diversity partitioning among varieties, races and populations. After aligning sequences of minisatellite repeat motifs, we evaluated evolutionary relationships among haplotypes.

KEY RESULTS

The geographical structuring of the 10 haplotypes corresponded with division between Pacific and Rocky Mountain varieties. Pacific haplotypes clustered with high bootstrap support, and appear to have descended from Rocky Mountain haplotypes. A greater proportion of diversity was partitioned between Rocky Mountain races than between Pacific races. Areas of highest haplotypic diversity were the southern Sierra Nevada mountain range in California, northwestern California, and southern Nevada.

CONCLUSIONS

Pinus ponderosa haplotype distribution patterns suggest a complex phylogeographic history not revealed by other genetic and morphological data, or by the sparse paleoecological record. The results appear consistent with long-term divergence between the Pacific and Rocky Mountain varieties, along with more recent divergences not well-associated with race. Pleistocene refugia may have existed in areas of high haplotypic diversity, as well as the Great Basin, Southwestern United States/northern Mexico, and the High Plains.

European invasion of North American Pinus strobus at large and fine scales: high genetic diversity and fine-scale genetic clustering over time in the adventive range

BACKGROUND

North American Pinus strobus is a highly invasive tree species in Central Europe. Using ten polymorphic microsatellite loci we compared various aspects of the large-scale genetic diversity of individuals from 30 sites in the native distribution range with those from 30 sites in the European adventive distribution range. To investigate the ascertained pattern of genetic diversity of this intercontinental comparison further, we surveyed fine-scale genetic diversity patterns and changes over time within four highly invasive populations in the adventive range.

RESULTS

Our data show that at the large scale the genetic diversity found within the relatively small adventive range in Central Europe, surprisingly, equals the diversity found within the sampled area in the native range, which is about thirty times larger. Bayesian assignment grouped individuals into two genetic clusters separating North American native populations from the European, non-native populations, without any strong genetic structure shown over either range. In the case of the fine scale, our comparison of genetic diversity parameters among the localities and age classes yielded no evidence of genetic diversity increase over time. We found that SGS differed across age classes within the populations under study. Old trees in general completely lacked any SGS, which increased over time and reached its maximum in the sapling stage.

CONCLUSIONS

Based on (1) the absence of difference in genetic diversity between the native and adventive ranges, together with the lack of structure in the native range, and (2) the lack of any evidence of any temporal increase in genetic diversity at four highly invasive populations in the adventive range, we conclude that population amalgamation probably first happened in the native range, prior to introduction. In such case, there would have been no need for multiple introductions from previously isolated populations, but only several introductions from genetically diverse populations.

Fine-scale spatial genetic structure and gene dispersal in Silene latifolia

Plants are sessile organisms, often characterized by limited dispersal. Seeds and pollen are the critical stages for gene flow. Here we investigate spatial genetic structure, gene dispersal and the relative contribution of pollen vs seed in the movement of genes in a stable metapopulation of the white campion Silene latifolia within its native range. This short-lived perennial plant is dioecious, has gravity-dispersed seeds and moth-mediated pollination. Direct measures of pollen dispersal suggested that large populations receive more pollen than small isolated populations and that most gene flow occurs within tens of meters. However, these studies were performed in the newly colonized range (North America) where the specialist pollinator is absent. In the native range (Europe), gene dispersal could fall on a different spatial scale. We genotyped 258 individuals from large and small (15) subpopulations along a 60 km, elongated metapopulation in Europe using six highly variable microsatellite markers, two X-linked and four autosomal. We found substantial genetic differentiation among subpopulations (global F(ST)=0.11) and a general pattern of isolation by distance over the whole sampled area. Spatial autocorrelation revealed high relatedness among neighboring individuals over hundreds of meters. Estimates of gene dispersal revealed gene flow at the scale of tens of meters (5-30 m), similar to the newly colonized range. Contrary to expectations, estimates of dispersal based on X and autosomal markers showed very similar ranges, suggesting similar levels of pollen and seed dispersal. This may be explained by stochastic events of extensive seed dispersal in this area and limited pollen dispersal.

Effect of inbreeding depression on outcrossing rates among populations of a tropical pine

Inbreeding depression is common among plants and may distort mating system estimates. Mating system studies traditionally ignore this effect, nonetheless an assessment of inbreeding depression that may have occurred before progeny evaluation could be necessary. In the neotropical Pinus chiapensis inbreeding depression was evaluated using regression analysis relating progeny F-values with seed germinability, the mating system was analysed in three populations with contrasting size, using isozymes, obtained a corrected outcrossing rate. Selfing decreased seed viability by 19%, relative to an outcrossed plant. Multilocus outcrossing rates, t(m), varied widely among populations. In the two smallest populations t(m) congruent with 1. Therefore, inbreeding depression did not affect the estimates, but overestimated t(m) by 10% in the third population, which has a true mixed mating system (selfing was the major source of inbreeding), and an unusually low t(m) for pines (t(m) = 0.54, uncorrected, t(m) = 0.49, corrected). Inbreeding depression may be an uneven source of bias for outcrossing estimates even at the infraspecific level. Accuracy [corrected] but not precision [corrected] may be gained by including inbreeding depression in outcrossing estimates. Therefore, caution should be taken when comparing t(m) among species or even populations within the same species.

Fat-tailed gene flow in the dioecious canopy tree species Fraxinus mandshurica var. japonica revealed by microsatellites

Pollen flow, seed dispersal and individual reproductive success can be simultaneously estimated from the genotypes of adults and offspring using stochastic models. Using four polymorphic microsatellite loci, gene flow of the wind-pollinated and wind-seed-dispersed dioecious tree species, Fraxinus mandshurica var. japonica, was quantified in a riparian forest, in northern Japan. In a 10.5-ha plot, 74 female adults, 76 male adults and 292 current-year seedlings were mapped and genotyped, together with 200 seeds. To estimate dispersal kernels of pollen and seeds, we applied normal, exponential power, Weibull, bivariate t-distribution kernels, and two-component models consisting of two normal distribution functions, one with a small and one with a large variance. A two-component pollen flow model with a small contribution (26.1%) from short-distance dispersal (sigma = 7.2 m), and the rest from long-distance flow (sigma = 209.9 m), was chosen for the best-fitting model. The average distance that integrated pollen flows inside and outside the study plot was estimated to be 196.8 m. Tree size and flowering intensity affected reproduction, and there appeared to be critical values that distinguished reproductively successful and unsuccessful adults. In contrast, the gene flow model that estimated both pollen and seed dispersal from established seedlings resulted in extensive seed dispersal, and the expected spatial genetic structures did not satisfactorily fit with the observations, even for the selected model. Our results advanced small-scale individual-based parentage analysis for quantifying fat-tailed gene flow in wind-mediated species, but also clarified its limitations and suggested future possibilities for gene flow studies.

A mitochondrial DNA minisatellite reveals the postglacial history of jack pine (Pinus banksiana), a broad-range North American conifer

Jack pine (Pinus banksiana Lamb.) is a broadly distributed North American conifer and its current range was covered by the Laurentian ice sheet during the last glacial maximum. To infer about the history and postglacial colonization of this boreal species, range-wide genetic variation was assessed using a new and highly variable minisatellite-like marker of the mitochondrial genome. Among the 543 trees analysed, 14 distinct haplotypes were detected, which corresponded to different repeat numbers of the 32-nucleotide minisatellite-like motif. Several haplotypes were rare with limited distribution, suggesting recent mutation events during the Holocene. At the population level, an average of 2.6 haplotypes and a mean haplotype diversity (H) of 0.328 were estimated. Population subdivision of genetic diversity was quite high with G(ST) and R(ST) values of 0.569 and 0.472, respectively. Spatial analyses identified three relatively homogeneous groups of populations presumably representative of genetically distinct glacial populations, one west and one east of the Appalachian Mountains in the United States and a third one presumably on the unglaciated northeastern coastal area in Canada. These results indicate the significant role of the northern part of the US Appalachian Mountains as a factor of vicariance during the ice age. A fourth distinct group of populations was observed in central Québec where the continental glacier retreated last. It included populations harbouring haplotypes present into the three previous groups, and it had higher level of haplotype diversity per population (H = 0.548) and lower population differentiation (G(ST) = 0.265), which indicates a zone of suture or secondary contact between the migration fronts of the three glacial populations. Introgression from Pinus contorta Dougl. var. latifolia Engelm. was apparent in one western population from Alberta. Altogether, these results indicate that the mitochondrial DNA variation of jack pine is geographically highly structured and it correlates well with large-scale patterns emerging from recent phylogeographical studies of other tree boreal species in North America.

Spatial structure and genetic diversity of two tropical tree species with contrasting breeding systems and different ploidy levels

Analyses of the spatial distribution pattern, spatial genetic structure and of genetic diversity were carried out in two tropical tree species with contrasting breeding systems and different ploidy levels using a 50-ha demographic plot in a lowland dipterocarp forest in Peninsular Malaysia. Shorea leprosula is a diploid and predominantly outcrossed species, whereas S. ovalis ssp. sericea is an autotetraploid species with apomictic mode of reproduction. Genetic diversity parameters estimated for S. leprosula using microsatellite were consistently higher than using allozyme. In comparisons with S. leprosula and other tropical tree species, S. ovalis ssp. sericea also displayed relatively high levels of genetic diversity. This might be explained by the lower pressure of genetic drift due to tetrasomic inheritance, and for autotetraploids each locus can accommodate up to four different alleles and this allows maintenance of more alleles at individual loci. The observed high levels of genetic diversity in S. ovalis ssp. sericea can also be due to a random retention of more heterogeneous individuals in the past, and the apomictic mode of reproduction might be an evolutionary strategy, which allows the species to maintain high levels of genetic diversity. The spatial distribution pattern analyses of both species showed significant levels of aggregation at small and medium but random distribution at the big diameter-class. The decrease in magnitude of spatial aggregation from small- to large-diameter classes might be due to compensatory mortality during recruitment and survival under competitive thinning process. Spatial genetic structure analyses for both species revealed significant spatial genetic structure for short distances in all the three diameter-classes. The magnitude of spatial genetic structure in both species was observed to be decreasing from smaller- to larger-diameter classes. The high spatial genetic structuring observed in S. ovalis ssp. sericea at the small-diameter class is due primarily to limited seed dispersal and apomictic mode of reproduction. The similar observation in S. leprosula, however, can be explained by limited seed and pollen dispersal, which supports further the fact that the species is pollinated by weak fliers, mainly of Thrips and Megalurothrips in the lowland dipterocarp forest.

Estimating seed vs. pollen dispersal from spatial genetic structure in the common ash

Spatial genetic structure was analysed with five highly polymorphic microsatellite loci in a Romanian population of common ash (Fraxinus excelsior L.), a wind-pollinated and wind-dispersed tree species occurring in mixed deciduous forests over almost all of Europe. Contributions of seed and pollen dispersal to total gene flow were investigated by analysing the pattern of decrease in kinship coefficients among pairs of individuals with geographical distance and comparing it with simulation results. Plots of kinship against the logarithm of distance were decomposed into a slope and a shape component. Simulations showed that the slope is informative about the global level of gene flow, in agreement with theoretical expectations, whereas the shape component was correlated with the relative importance of seed vs. pollen dispersal. Hence, our results indicate that insights into the relative contributions of seed and pollen dispersal to overall gene flow can be gained from details of the pattern of spatial genetic structure at biparentally inherited loci. In common ash, the slope provided an estimate of total gene dispersal in terms of Wright's neighbourhood size of Nb = 519 individuals. No precise estimate of seed vs. pollen flow could be obtained from the shape because of the stochasticity inherent to the data, but the parameter combinations that best fitted the data indicated restricted seed flow, sigmas pound 14 m, and moderate pollen flow, 70 m pound sigmap pound 140 m.

Evidence for gene flow via seed dispersal from crop to wild relatives in Beta vulgaris (Chenopodiaceae): consequences for the release of genetically modified crop species with weedy lineages

Gene flow and introgression from cultivated to wild plant populations have important evolutionary and ecological consequences and require detailed investigations for risk assessments of transgene escape into natural ecosystems. Sugar beets (Beta vulgaris ssp. vulgaris) are of particular concern because: (i) they are cross-compatible with their wild relatives (the sea beet, B. vulgaris ssp. maritima); (ii) crop-to-wild gene flow is likely to occur via weedy lineages resulting from hybridization events and locally infesting fields. Using a chloroplastic marker and a set of nuclear microsatellite loci, the occurrence of crop-to-wild gene flow was investigated in the French sugar beet production area within a 'contact-zone' in between coastal wild populations and sugar beet fields. The results did not reveal large pollen dispersal from weed to wild beets. However, several pieces of evidence clearly show an escape of weedy lineages from fields via seed flow. Since most studies involving the assessment of transgene escape from crops to wild outcrossing relatives generally focused only on pollen dispersal, this last result was unexpected: it points out the key role of a long-lived seed bank and highlights support for transgene escape via man-mediated long-distance dispersal events.

Diverging patterns of mitochondrial and nuclear DNA diversity in subarctic black spruce: imprint of a founder effect associated with postglacial colonization

High-latitude ecotonal populations at the species margins may exhibit altered patterns of genetic diversity, resulting from more or less recent founder events and from bottleneck effects in response to climate oscillations. Patterns of genetic diversity were investigated in nine populations of the conifer black spruce (Picea mariana [Mill.] BSP.) in northwestern Québec, Canada, using seed-dispersed mitochondrial (mt) DNA and nuclear (nc) DNA. mtDNA diversity (mitotypes) was assessed at three loci, and ncDNA diversity was estimated for nine expressed sequence tag polymorphism (ESTP) loci. Sampling included populations from the boreal forest and the southern and northern subzones of the subarctic forest-tundra, a fire-born ecotone. For ncDNA, populations from all three vegetation zones were highly diverse with little population differentiation (thetaN = 0.014); even the northernmost populations showed no loss of rare alleles. Patterns of mitotype diversity were strikingly different: within-population diversity and population differentiation were high for boreal forest populations [expected heterozygosity per locus (HE) = 0.58 and thetaM = 0.529], but all subarctic populations were fixed for a single mitotype (HE = 0). This lack of variation suggests a founder event caused by long-distance seed establishment during postglacial colonization, consistent with palaeoecological data. The estimated movement of seeds alone (effective number of migrants per generation, NmM < 2) was much restricted compared to that estimated from nuclear variants, which including pollen movement (NmN > 17). This could account for the conservation of a founder imprint in the mtDNA of subarctic black spruce. After reduction, presumably in the early Holocene, the diversity in ncDNA would have been replenished rapidly by pollen-mediated gene flow, and maintained subsequently through vegetative layering during the current cooler period covering the last 3000 years.

Mitochondrial haplotype distribution, seed dispersal and patterns of postglacial expansion of ponderosa pine

Maternally inherited mtDNA in a secondary contact zone of ponderosa pine revealed a cline less than 10 km wide - much narrower than previously described. A survey of 76 populations gave no evidence either of intermixing or of a mosaic contact zone. Such sharp contact zones are consistent with diffusive range expansion, rather than long distance colonization. However, evidence for long distance seed dispersal events was found in two populations where haplotypes were observed far from their main area of occurrence. The results suggest a small number of long distance colonists with diffusive dispersal from these centres.

Spatial autocorrelation and linkage of Mendelian RAPD markers in a population of Picea abies Karst

The spatial clustering of single- and di-locus genotypes in a natural, continuous population of Norway spruce was investigated using 69 Mendelian Random Amplified Polymorphic DNA (RAPD) markers that covered about 15 of the species genome, and whose linkage relationships were known. Spatial autocorrelation techniques and randomization tests, applied to both single- and di-locus genotypes, revealed a weak, though significant, spatial structure at the scale 0-200 m (5 of single-locus and 7 of di-locus genotypes). To assess the relative importance of isolation by distance and linkage between markers on their spatial genetic structuring, we grouped joins between sampled trees into equivalence categories expected to show similar, specific patterns of spatial distribution under isolation by distance. Results from both single- and di-locus analyses were consistent with the existence of patches of like homozygotes (about 8 and 11 of loci at the single- and di-locus level, respectively) surrounded by a mix of like heterozygotes. Similar structuring has been predicted by simulation models under isolation by distance and selective neutrality. Overall, linkage between markers accounted for an increase of spatial clumping of di-locus genotypes involving tightly linked loci with recombination fractions up to 0.1, a consequence of limited, stochastic spread of single-locus genotypes in space. Our results support the hypothesis that isolation by distance and linkage have a small, though significant, effect even within continuous forest tree populations. In general, the spatial distribution of multilocus genotypes within populations should be interpreted with caution when linkage relationships among the markers used are unknown.

Landscape genetic structure of Pinus banksiana: allozyme variation

Many of Wisconsin's jack pine forests originated following fire or agricultural abandonment creating a forest mosaic fragmented by a history of disturbance and past land use. The extent and patterning of genetic diversity at a landscape scale (30 × 30 km) was investigated in 82 natural Pinus banksiana Lamb. (jack pine) stands in Wisconsin using 14 polymorphic allozymes. Most measures of genetic diversity and overall allelic frequencies varied little among these stands, and Reynolds' (coancestry) genetic distances were small (mean = 0.026). Genetic differentiation among stands was limited but significant ([Formula: see text]ST = 0.022). Gene flow (Nm) is extensive, estimated to be slightly more than 11 migrants per generation. Autocorrelation analysis provided evidence for a weak pattern of genetic structure at a few loci with a spatial scale of 8–15 km. However, when all loci were examined together the populations did not present a clear spatial pattern across the landscape, probably because of extensive gene flow among stands.Key words: jack pine, allozymes, genetic diversity, population differentiation, spatial autocorrelation, gene flow.

Fine-scale genetic structure in Pinus clausa (Pinaceae) populations: effects of disturbance history

Spatial autocorrelation analyses of 12 allozyme loci were used to compare genetic structure within populations of two varieties of Pinus clausa. P. clausa var. immuginata populations tend to be uneven-aged, with continuous recruitment in small gaps created by wind damage, whereas P. clausa var. clausa populations are more even-aged, with recruitment postdating periodic canopy fires. Three var. immuginata populations and three matched pairs of var. clausa populations, including both a mature and a nearby recently burned population, were examined. Aggregation of multilocus genotypes at small distances was evident in all young var. clausa populations. Little inbreeding was apparent among juveniles or adults in these populations; their genetic structure is likely to have resulted from limited seed dispersal. Genotypes were not significantly spatially structured in nearby matched mature populations. Genetic structure was less evident in var. immuginata populations. Aggregated genotypes were only apparent in the population where patches included juveniles of similar ages; dense juvenile clumps in the other two var. immuginata populations comprised a variety of ages. Interannual variability in allele frequencies of surviving seedlings may account for the absence of genetic structure in these populations.

Outcrossing and paternity analysis of Pinus densiflora (Japanese red pine) by microsatellite polymorphism

This study employed microsatellite loci to analyse outcrossing rate and pollen dispersal in Japanese red pine (Pinus densiflora) in an isolated stand. The average offspring outcrossing rate for 29 cones was 0.955. Significant differences in outcrossing rates between offspring groups on individual branches that extended in different directions at different heights were not detected. Male parents of 874 offspring collected from the maternal tree were assessed by exclusion using polymorphisms at three microsatellite loci. Paternity analysis indicated that at least 31% of the offspring were fertilized by pollen from trees outside the stand. The average distance of pollen migration within the study stand was 68 m, with a maximum value of 325 m. There was excess mating with nearby P. densiflora trees, of which only a few were predominant pollen donors. In addition, a weakly directional bias in P. densiflora pollination was also detected in the study stand, suggesting that female strobili on a branch of the maternal tree were more easily fertilized by pollen from trees in that direction.

Cytonuclear disequilibrium and genetic drift in a natural population of ponderosa pine

We measured the cytonuclear disequilibrium between 11 nuclear allozyme loci and both mitochondrial and chloroplast DNA haplotypes in a natural population of ponderosa pine (Pinus ponderosa, Laws). Three allozyme loci showed significant associations with mtDNA variation, while two other loci showed significant association with cpDNA. However, the absolute number of individuals involved in any of the associations was small, such that in none of the nuclear-organellar combinations was the difference between observed and expected numbers >11 individuals. Patterns of association were not consistent across loci or organellar genomes, suggesting that they are not the result of mating patterns, which would act uniformly on all loci. This pattern of disequilibria is consistent with the action of genetic drift and with existing knowledge of the structure of this population and thus does not imply the action of other evolutionary processes. The overall magnitude (normalized disequilibrium) of associations was greater for maternally inherited mtDNA than for paternally inherited cpDNA, though this difference was neither large nor significant. Such significant disequilibria involving the paternally inherited organelle indicate that not only are there a limited number of seed parents, but the effective number of pollen parents is also limited.

Chloroplast microsatellites and mitochondrial nad1 intron 2 sequences indicate congruent phylogenetic relationships among Swiss stone pine (Pinus cembra), Siberian stone pine (Pinus sibirica), and Siberian dwarf pine (Pinus pumila)

We studied the phylogenetic relationships among the three stone pine species, Pinus cembra, P. sibirica, and P. pumila, using chloroplast microsatellites and mitochondrial nad1 intron 2 sequences. The three chloroplast microsatellite loci combined into a total of 18 haplotypes. Fourteen haplotypes were detected in 15 populations of P. cembra and one population of P. sibirica, five of which were shared between the two species, and the two populations of P. pumila comprised four species-specific haplotypes. Mitochondrial intron sequences confirmed this grouping of species. Sequences of P. cembra and P. sibirica were identical, but P. pumila differed by several nucleotide substitutions and insertions/deletions. A repeat region found in the former two species showed no intraspecific variation. These results indicate a relatively recent evolutionary separation of P. cembra and P. sibirica, despite their currently disjunct distributions. The species-specific chloroplast and mitochondrial markers of P. sibirica and P. pumila should help to trace the hybridization in their overlapping distribution area and to identify fossil remains with respect to the still unresolved postglacial re-colonization history of these two species.

Diversity and genetic structure in populations of Pseudotsuga menziesii (Pinaceae) at chloroplast microsatellite loci

Genetic variation was compared between uniparentally-inherited (chloroplast simple sequence repeats, cpSSRs) vs. biparentally-inherited (isozyme and random amplified polymorphic DNA, RAPD) genetic markers in Douglas-fir (Pseudotsuga mensiezii) from British Columbia. Three-hundred twenty-three individuals from 11 populations were assayed. In Douglas-fir, the cpSSR primer sites were well-conserved relative to Pinus thunbergii (11 of 17 loci clearly amplified), but only 3 loci were appreciably polymorphic. At these cpSSR loci, we found an unexpectedly low level of polymorphism within populations, and no genetic differentiation among populations. By contrast, the nuclear markers showed variation typical of conifers, with significant among-population differentiation. This difference is likely the outcome of both historical factors and high pollen dispersal.

Tandem repeats in plant mitochondrial genomes: application to the analysis of population differentiation in the conifer Norway spruce

Mitochondrial DNA, widely applied in studies of population differentiation in animals, is rarely used in plants because of its slow rate of sequence evolution and its complex genomic organization. We demonstrate the utility of two polymorphic mitochondrial tandem repeats located in the second intron of the nad1 gene of Norway spruce. Most of the size variants showed pronounced population differentiation and a distinct geographical distribution. A GenBank search revealed that mitochondrial tandem repeats occur in a broad range of plant species and may serve as a novel molecular marker for unravelling population processes in plants.

Paternity and gene dispersal in limber pine (Pinus flexilis James)

This study provides empirical information on intrapopulation gene dispersal via pollen, the size of genetic 'neighbourhoods', and interpopulation gene flow in a long-lived conifer, limber pine (Pinus flexilis). We used allozyme loci for a paternity analysis of 518 seeds produced in an isolated population of limber pine located in north-eastern Colorado, U.S.A., separated by 2 km from the nearest conspecific trees and nearly 100 km from populations in the Rocky Mountains. We also used indirect techniques (FST analyses) to estimate gene flow rates among subdivisions of the study population and among five widely separated populations. Within the main study population limber pine exhibited a polymorphism level of 50%, observed heterozygosity of 0.159, and 2.36 alleles per polymorphic locus. Mountain populations were slightly more variable. The main study population showed significant differentiation in allozymes among neighbouring subpopulations. The mean FST was 0.031 and the gene flow rate among subpopulations was estimated as 7.8 migrants per generation. Among widely separated populations the mean FST was 0.035 and the gene flow rate was estimated as 6.9 migrants per generation. The paternity analysis indicated a best estimate of 6.5% pollen immigration (minimum 1.1%) from populations 2 km to 100+ km away. For 4% of the seeds examined, paternity could be ascribed to a single tree in the study population. Fractional paternity and likelihood methods were used to estimate pollen dispersal distances for the remainder of the seeds. Mean pollen dispersal distance was estimated at 140 m using the fractional method, similar to results from the other techniques. This compares with a mean distance of 172 m between potential mates. These results suggest near-panmictic pollen dispersal over this population, which covers about 15 ha. The observed allozyme differences and surprisingly low estimates of among-subpopulation gene flow are ascribed to a probable restriction of gene dispersal by seeds.

Spatial genetic structure of delphinium nuttallianum populations: inferences about gene flow

The spatial genetic structure of a plant population provides a potential record of past gene flow and mating. We used hierarchical F-statistics and spatial autocorrelation to characterize spatial genetic differentiation of allozymes in adult Delphinium nuttallianum plants within and among six natural populations separated from one another by up to 3 km. Previous direct estimates suggested that gene flow is highly localized, averaging << 10 m. Earlier studies of seed-set, pollen-tube growth and progeny fitness suggested that partial reproductive isolation exists between plants growing too close together (<3 m) and too far apart (>100 m). Thus we anticipated substantial genetic differentiation on scales of a few to hundreds of metres. However, we detected little differentiation among the six populations, among replicate study plots within populations, or among subsections of study plots, except at the smallest scale of cm to m. These results suggest that relatively rare long-distance pollen movement has gone undetected and that postpollination selection may further modify genetic structure during the life cycle. Lack of differentiation is not at odds with the observation of partial reproductive isolation, because some loci may respond to spatial variation in selection without this response being evident at marker loci.