Decoupled mitochondrial and chloroplast DNA population structure reveals Holocene collapse and population isolation in a threatened Mexican-endemic conifer

Chihuahua spruce (Picea chihuahuana Martínez) is a montane subtropical conifer endemic to the Sierra Madre Occidental in northwestern México. Range-wide variation was investigated using maternally inherited mitochondrial (mtDNA) and paternally inherited chloroplast (cpDNA) DNA markers. Among the 16 mtDNA regions analysed, only two mitotypes were detected, while the study of six cpDNA microsatellite markers revealed eight different chlorotypes. The average cpDNA diversity (H = 0.415) was low but much higher than that for mtDNA (H = 0). The distribution of mitotypes revealed two clear nonoverlapping areas (G(ST) = N(ST) = 1), one including northern populations and the second one including the southern and central stands, suggesting that these two regions may represent different ancestral populations. The cpDNA markers showed lower population differentiation (G(ST) = 0.362; R(ST) = 0.230), implying that the two ancestral populations continued to exchange pollen after their initial geographic separation. A lack of a phylogeographic structure was revealed by different spatial analyses of cpDNA (G(ST) > R(ST); and samova), and reduced cpDNA gene flow was noted among populations (Nm = 0.873). Some stands deviated significantly from the mutation-drift equilibrium, suggesting recent bottlenecks. Altogether, these various trends are consistent with the hypothesis of a population collapse during the Holocene warming and suggest that most of the modern P. chihuahuana populations are now effectively isolated with their genetic diversity essentially modelled by genetic drift. The conservation efforts should focus on most southern populations and on the northern and central stands exhibiting high levels of genetic diversity. Additional mtDNA sequence analysis confirmed that P. martinezii (Patterson) is not conspecific with P. chihuahuana, and thus deserves separate conservation efforts.     

Phylogeographic structure of jack pine (Pinus banksiana; Pinaceae) supports the existence of a coastal glacial refugium in northeastern North America

? Premise of the study: The genetic structure of jack pine (Pinus banksiana Lamb.), a North American boreal conifer with large longitudinal distribution, was investigated to test for the possible existence of a genetically distinct lineage in the Maritimes region in northeastern North America, which could be indicative of a mid-latitude coastal refuge during the last glaciation. ? Methods: One maternally inherited mitochondrial DNA (mtDNA) minisatellite marker and four paternally inherited chloroplast DNA (cpDNA) microsatellite markers were used to assess the range-wide geographical structure of jack pine populations with particular focus on northeastern North America. ? Key results: The populations from the Maritimes region presented a unique mtDNA background characterized by very low diversity and the preponderance of a distinctive mitotype. The distribution of cpDNA diversity was not spatially structured, though three chlorotypes were restricted to the east. ? Conclusions: MtDNA data suggest that populations from the Maritimes region derive from a genetically depauperated north-coastal refugium. Contrastingly, the much higher geographical uniformity observed for cpDNA variation indicates that gene flow by pollen had been much more effective than seed gene flow at homogenizing population structure.     

A selective sweep in the chloroplast DNA of dioecious silene (section Elisanthe)

Gene flow occurs predominantly via pollen in angiosperms, leading to stronger population subdivision for maternally inherited markers, relative to paternally or biparentally inherited genes. In contrast to this trend, population subdivision within Silene latifolia and S. dioica, as well as subdivision between the two species, is substantially lower in maternally inherited chloroplast genes compared to paternally inherited Y-linked genes. A significant frequency spectrum bias toward rare polymorphisms and a significant loss of polymorphism in chloroplast genes compared to Y-linked and autosomal genes suggest that intra- and inter-specific subdivision in the chloroplast DNA may have been eroded by a selective sweep that has crossed the S. latifolia and S. dioica species boundary.     

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

We examined the spatial distribution of maternally inherited mitochondrial DNA and paternally inherited chloroplast DNA polymorphisms in a permanently marked stand of ponderosa pine (Pinus ponderosa Laws). Movement of maternally inherited mtDNA occurs only via seed dispersal, and mtDNA haplotypes showed significant patch structure. Moreover, individuals within patches identified by mtDNA haplotypes were related approximately as half-sibs based upon analysis of allozyme genotypes. Thus, seed dispersal is limited within the population, and creates matrilineal clusters in space. By contrast, paternally inherited cpDNA is dispersed by movement of both seed and pollen. Chloroplast DNA polymorphisms showed no evidence of patch structure, but rather a weak (and nonsignificant) trend toward hyperdispersion, suggesting nearly unlimited movement of pollen among trees within this stand. Two of the trees had unique allozyme alleles, which were used to directly measure pollen movement away from those trees. Marked pollen was as likely to disperse across the population as it was to fertilize near neighbors.     

Impacts of seed and pollen flow on population genetic structure for plant genomes with three contrasting modes of inheritance

The classical island and one-dimensional stepping-stone models of population genetic structure developed for animal populations are extended to hermaphrodite plant populations to study the behavior of biparentally inherited nuclear genes and organelle genes with paternal and maternal inheritance. By substituting appropriate values for effective population sizes and migration rates of the genes concerned into the classical models, expressions for genetic differentiation and correlation in gene frequency between populations can be derived. For both models, differentiation for maternally inherited genes at migration-drift equilibrium is greater than that for paternally inherited genes, which in turn is greater than that for biparentally inherited nuclear genes. In the stepping-stone model, the change of genetic correlation with distance is influenced by the mode of inheritance of the gene and the relative values of long- and short-distance migration by seed and pollen. In situations where it is possible to measure simultaneously Fst for genes with all three types of inheritance, estimates of the relative rates of pollen to seed flow can be made for both the short- and long-distance components of migration in the stepping-stone model.     

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.     

HISTORICAL SEPARATION AND PRESENT GENE FLOW THROUGH A ZONE OF SECONDARY CONTACT IN PONDEROSA PINE

I examined the effects of historical division and secondary contact between eastern and western varieties of ponderosa pine (Pinus ponderosa Laws Pinaceae) on extant patterns of genetic variation. Fossil and biogeographic evidence both indicate that the current point of contact between these two varieties represents secondary contact following historical separation during the Wisconsin glaciation. Current gene flow was assessed by observing the degree of introgression of paternally inherited cpDNA and maternally inherited mtDNA polymorphisms. Both seeds and pollen are wind dispersed in ponderosa pine. Introgression was primarily from west to east, the direction of the prevailing wind, for both organelles, but introgression of cpDNA far exceeded that of mtDNA. Thus pollen is the main agent of contemporary gene flow between the two varieties. Neither seeds nor pollen showed enough introgression since secondary contact to have homogenized the two gene pools. However, allozyme differentiation was minimal. This calls into question assumptions of selective neutrality for at least some of the markers. Theory predicts that nuclear markers will show a high locus-to-locus variance of F_ST following historical separation. This prediction is confirmed by the allozyme data for ponderosa pine, and may provide a useful means of identifying historical separations from allele frequency data.     

Regional population dynamics define the local genetic structure in Sorbus torminalis

Recent changes in sylvicultural practices in Central Europe have created forests with closed canopies, and tree species preferring open and sunny forests have declined in area and abundance. This led to increased isolation of populations of many rare insect-pollinated, fleshy-fruited species with a naturally scattered distribution. To gain insight into the regional population dynamics of such species, we investigated the consequences of spatial isolation, population size and density on the genetic structure of Sorbus torminalis and simultaneously considered the relationship between fecundity and habitat quality. Genotype data for biparentally (ISSRs) and maternally inherited (cpDNA PCR-RFLPs) molecular markers were generated for 26 Swiss populations of S. torminalis. We applied analyses of molecular variance (amova) to both marker types and separated the relative contributions of pollen and seed dispersal to historical gene flow. amova detected significant differentiation among populations (Phi(ST ISSR) = 0.107; Phi(ST cpDNA) = 0.370) in both marker types. The relative rate of pollen to seed gene flow was low (r = 2.919) and significantly different from equality. Isolation by distance was weak within Eastern and Western Switzerland, although populations were substantially differentiated. Within-population molecular variance was not explained by population size, whereas habitat quality (openness) positively influenced the percentage of fruiting trees and the degree of fruiting per tree, indicating that more open forests enhance sexual reproduction. Our findings of significant genetic differentiation in the absence of clear geographical structuring can be explained by the distinct ecology of S. torminalis and nondirectional colonization events in metapopulation-like dynamics.     

From glacial refugia to modern populations: new assemblages of organelle genomes generated by differential cytoplasmic gene flow in transcontinental black spruce

Assessing species' range-wide cytoplasmic diversity provides valuable insights regarding their dispersal and adaptive potential in a changing environment. Transcontinental chloroplast (cpDNA) and mitochondrial DNA (mtDNA) population structures were compared to identify putative ancestral and new cytoplasmic genome assemblages in black spruce (Picea mariana), a North American boreal conifer. Mean within-population diversity and allelic richness for cpSSR markers were 0.80 and 4.21, respectively, and diminished westward. Population differentiation based on G(ST) was lower for cpDNA than for mtDNA (G(ST) =0.104 and 0.645, respectively) but appeared comparable when estimated using Jost differentiation index (D=0.459 and 0.537, respectively). Further analyses resulted in the delineation of at least three genetically distinct cpDNA lineages partially congruent with those inferred from mtDNA data, which roughly corresponded to western, central and eastern Canada. Additionally, the patterns of variation in Alaska for both cpDNA and mtDNA markers suggested that black spruce survived the last glacial maximum in this northern region. The range-wide comparison of the geographic extent of cytoplasmic DNA lineages revealed that extensive pollen gene flow between ancestral lineages occurred preferentially from west to east during the postglacial expansion of the species, while seed-mediated gene flow remained geographically restricted. This differential gene flow promoted intraspecific cytoplasmic capture that generated new assemblages of cpDNA and mtDNA genomes during the Holocene. Hence, black spruce postglacial colonization unexpectedly resulted in an increase in genetic diversity with possible adaptive consequences.     

Spatial genetic structure among and within populations of Primula sieboldii growing beside separate streams

We investigated the hierarchical genetic structure of SSR (simple sequence repeats) and cpDNA (chloroplast DNA) polymorphisms among and within populations of Primula sieboldii, a heterostylous clonal herb. Seven out of eight populations at the study site, located in a mountainous region of Nagano Prefecture, had each developed alongside a different stream, and the other occurred on a flat area 70 m from the nearest stream. The magnitude of genetic differentiation among streamside populations in maternally inherited cpDNA (Phi = 0.341) was much higher than that in biparentally inherited SSRs (Phi = 0.011). This result suggests that seed dispersal among streams was restricted, and pollen was the primary agent of gene flow among streamside populations. In contrast, genetic differentiation among subpopulations within streams were low at both markers (Phi = 0.053 for cpDNA, Phi = 0.025 for SSR). This low differentiation among subpopulations in cpDNA compared with that among streamside populations suggest that seed dispersal occur along the stream probably during flooding. This hypothesis was supported by the fact that in cpDNA haplotypes, no clear genetic structure was detected within the streamside population, while a significant genetic structure was found within 20 m in the nonstreamside population. Furthermore, within the streamside populations, two pairs of ramets with identical multilocus genotypes for eight SSR loci were distantly (> 50 m) distributed along the same streamside, suggesting dispersal of clonal propagule. Our study showed that the heterogeneity of the landscape can influence gene flow and hence spatial genetic structure.