Glacial refugia of limber pine (Pinus flexilis James) inferred from the population structure of mitochondrial DNA

Plant mitochondrial introns as genetic markers - conservation and variation

Plant genomes are comprised of nuclear, plastid and mitochondrial components characterized by different patterns of inheritance and evolution. Genetic markers from the three genomes provide complementary tools for investigations of inheritance, genetic relationships and phenotypic contributions. Plant mitochondrial genomes are challenging for universal marker development because they are highly variable in terms of size, gene order and intergenic sequences and highly conserved with respect to protein-coding sequences. PCR amplification of introns with primers that anneal to conserved, flanking exons is effective for the development of polymorphic nuclear genome markers. The potential for plant mitochondrial intron polymorphisms to distinguish between congeneric species or intraspecific varieties has not been systematically investigated and is possibly constrained by requirements for intron secondary structure and interactions with co-evolved organelle intron splicing factors. To explore the potential for broadly applicable plant mitochondrial intron markers, PCR primer sets based upon conserved sequences flanking 11 introns common to seven angiosperm species were tested across a range of plant orders. PCR-amplified introns were screened for indel polymorphisms among a group of cross-compatible Citrus species and relatives; two Raphanus sativus mitotypes; representatives of the two Phaseolus vulgaris gene pools; and congeneric pairs of Cynodon, Cenchrus, Solanum, and Vaccinium species. All introns were successfully amplified from each plant entry. Length polymorphisms distinguishable by gel electrophoresis were common among genera but infrequent within genera. Sequencing of three introns amplified from 16 entries identified additional short indel polymorphisms and nucleotide substitutions that separated Citrus, Cynodon, Cenchrus and Vaccinium congeners, but failed to distinguish Solanum congeners or representatives of the Phaseolus vulgaris major gene pools. The ability of primer sets to amplify a wider range of plant species' introns and the presence of intron polymorphisms that distinguish congeners was confirmed by in silico analysis. While mitochondrial intron variation is limited in comparison to nuclear introns, these exon-based primer sets provide robust tools for the amplification of mitochondrial introns across a wide range of plant species wherein useful polymorphisms can be identified.

Hybrid breakdown is elevated near the historical cores of a species' range

New species form when they become reproductively isolated. A classic model of speciation posits that derived mutations appear in isolated populations and reduce fitness when combined in hybrids. While these Bateson-Dobzhansky-Muller incompatibilities are known to accumulate as populations diverge over time, they may also reflect the amount of standing genetic variation within populations. We analysed the fitness of F₂ hybrids in crosses between 24 populations of a plant species (Campanula americana) with broad variation in standing genetic variation and genetic differentiation driven by post-glacial range expansions. Hybrid breakdown varied substantially and was strongest between populations near the historical cores of the species range where within-population genetic diversity was high. Nearly half of the variation in hybrid breakdown was predicted by the combined effects of standing genetic variation within populations, their pairwise genetic differentiation and differences in the climates they inhabit. Hybrid breakdown was enhanced between populations inhabiting distinct climates, likely reflecting local adaptation. Results support that the mutations causing hybrid breakdown, the raw material for speciation, are more common in long-inhabited areas of the species range. Genetic diversity harboured in refugial areas is thus an important source of incompatibilities critical to the speciation process.

Adaptive evolution in a conifer hybrid zone is driven by a mosaic of recently introgressed and background genetic variants

Extant conifer species may be susceptible to rapid environmental change owing to their long generation times, but could also be resilient due to high levels of standing genetic diversity. Hybridisation between closely related species can increase genetic diversity and generate novel allelic combinations capable of fuelling adaptive evolution. Our study unravelled the genetic architecture of adaptive evolution in a conifer hybrid zone formed between Pinus strobiformis and P. flexilis. Using a multifaceted approach emphasising the spatial and environmental patterns of linkage disequilibrium and ancestry enrichment, we identified recently introgressed and background genetic variants to be driving adaptive evolution along different environmental gradients. Specifically, recently introgressed variants from P. flexilis were favoured along freeze-related environmental gradients, while background variants were favoured along water availability-related gradients. We posit that such mosaics of allelic variants within conifer hybrid zones will confer upon them greater resilience to ongoing and future environmental change and can be a key resource for conservation efforts.

Genotyping-by-sequencing and ecological niche modeling illuminate phylogeography, admixture, and Pleistocene range dynamics in quaking aspen (Populus tremuloides)

Populus tremuloides is the widest-ranging tree species in North America and an ecologically important component of mesic forest ecosystems displaced by the Pleistocene glaciations. Using phylogeographic analyses of genome-wide SNPs (34,796 SNPs, 183 individuals) and ecological niche modeling, we inferred population structure, ploidy levels, admixture, and Pleistocene range dynamics of P. tremuloides, and tested several historical biogeographical hypotheses. We found three genetic lineages located mainly in coastal-Cascades (cluster 1), east-slope Cascades-Sierra Nevadas-Northern Rockies (cluster 2), and U.S. Rocky Mountains through southern Canadian (cluster 3) regions of the P. tremuloides range, with tree graph relationships of the form ((cluster 1, cluster 2), cluster 3). Populations consisted mainly of diploids (86%) but also small numbers of triploids (12%) and tetraploids (1%), and ploidy did not adversely affect our genetic inferences. The main vector of admixture was from cluster 3 into cluster 2, with the admixture zone trending northwest through the Rocky Mountains along a recognized phenotypic cline (Utah to Idaho). Clusters 1 and 2 provided strong support for the "stable-edge hypothesis" that unglaciated southwestern populations persisted in situ since the last glaciation. By contrast, despite a lack of clinal genetic variation, cluster 3 exhibited "trailing-edge" dynamics from niche suitability predictions signifying complete northward postglacial expansion. Results were also consistent with the "inland dispersal hypothesis" predicting postglacial assembly of Pacific Northwestern forest ecosystems, but rejected the hypothesis that Pacific-coastal populations were colonized during outburst flooding from glacial Lake Missoula. Overall, congruent patterns between our phylogeographic and ecological niche modeling results and fossil pollen data demonstrate complex mixtures of stable-edge, refugial locations, and postglacial expansion within P. tremuloides. These findings confirm and refine previous genetic studies, while strongly supporting a distinct Pacific-coastal genetic lineage of quaking aspen.

Spatial genetic structure of Rocky Mountain bighorn sheep (Ovis canadensis canadensis) at the northern limit of their native range

The Canadian Rocky Mountains are one of the few places on Earth where the spatial genetic structure of wide-ranging species has been relatively unaffected by anthropogenic disturbance. We characterized the spatial genetic structure of Rocky Mountain bighorn sheep (Ovis canadensis canadensis Shaw, 1804) in the northern portion of their range. Using microsatellites from 1495 individuals and mitochondrial DNA sequences from 188 individuals, we examined both broad- and fine-scale spatial genetic structure, assessed sex-biased gene flow within the northern portion of the species range, and identified geographic patterns of genetic diversity. We found that broad-scale spatial genetic structure was consistent with barriers to movement created by major river valleys. The fine-scale spatial genetic structure was characterized by a strong isolation-by-distance pattern, and analysis of neighborhood size using spatial autocorrelation indicated gene flow frequently occurred over distances of up to 100 km. However, analysis of sex-specific spatial autocorrelation and analysis of mitochondrial haplotype distributions failed to detect any evidence of sex-biased gene flow. Finally, our analyses reveal decreasing genetic diversity with increasing latitude, consistent with patterns of post-glacial recolonization of the Rocky Mountains.

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.

Low offspring survival in mountain pine beetle infesting the resistant Great Basin bristlecone pine supports the preference-performance hypothesis

The preference-performance hypothesis states that ovipositing phytophagous insects will select host plants that are well-suited for their offspring and avoid host plants that do not support offspring performance (survival, development and fitness). The mountain pine beetle (Dendroctonus ponderosae), a native insect herbivore in western North America, can successfully attack and reproduce in most species of Pinus throughout its native range. However, mountain pine beetles avoid attacking Great Basin bristlecone pine (Pinus longaeva), despite recent climate-driven increases in mountain pine beetle populations at the high elevations where Great Basin bristlecone pine grows. Low preference for a potential host plant species may not persist if the plant supports favorable insect offspring performance, and Great Basin bristlecone pine suitability for mountain pine beetle offspring performance is unclear. We infested cut bolts of Great Basin bristlecone pine and two susceptible host tree species, limber (P. flexilis) and lodgepole (P. contorta) pines with adult mountain pine beetles and compared offspring performance. To investigate the potential for variation in offspring performance among mountain pine beetles from different areas, we tested beetles from geographically-separated populations within and outside the current range of Great Basin bristlecone pine. Although mountain pine beetles constructed galleries and laid viable eggs in all three tree species, extremely few offspring emerged from Great Basin bristlecone pine, regardless of the beetle population. Our observed low offspring performance in Great Basin bristlecone pine corresponds with previously documented low mountain pine beetle attack preference. A low preference-low performance relationship suggests that Great Basin bristlecone pine resistance to mountain pine beetle is likely to be retained through climate-driven high-elevation mountain pine beetle outbreaks.

Altitudinal gradients, biogeographic history and microhabitat adaptation affect fine-scale spatial genetic structure in African and Neotropical populations of an ancient tropical tree species

The analysis of fine-scale spatial genetic structure (FSGS) within populations can provide insights into eco-evolutionary processes. Restricted dispersal and locally occurring genetic drift are the primary causes for FSGS at equilibrium, as described in the isolation by distance (IBD) model. Beyond IBD expectations, spatial, environmental or historical factors can affect FSGS. We examined FSGS in seven African and Neotropical populations of the late-successional rain forest tree Symphonia globulifera L. f. (Clusiaceae) to discriminate the influence of drift-dispersal vs. landscape/ecological features and historical processes on FSGS. We used spatial principal component analysis and Bayesian clustering to assess spatial genetic heterogeneity at SSRs and examined its association with plastid DNA and habitat features. African populations (from Cameroon and São Tomé) displayed a stronger FSGS than Neotropical populations at both marker types (mean Sp = 0.025 vs. Sp = 0.008 at SSRs) and had a stronger spatial genetic heterogeneity. All three African populations occurred in pronounced altitudinal gradients, possibly restricting animal-mediated seed dispersal. Cyto-nuclear disequilibria in Cameroonian populations also suggested a legacy of biogeographic history to explain these genetic patterns. Conversely, Neotropical populations exhibited a weaker FSGS, which may reflect more efficient wide-ranging seed dispersal by Neotropical bats and other dispersers. The population from French Guiana displayed an association of plastid haplotypes with two morphotypes characterized by differential habitat preferences. Our results highlight the importance of the microenvironment for eco-evolutionary processes within persistent tropical tree populations.

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.

Linking carbon and water relations to drought-induced mortality in Pinus flexilis seedlings

Survival of tree seedlings at high elevations has been shown to be limited by thermal constraints on carbon balance, but it is unknown if carbon relations also limit seedling survival at lower elevations, where water relations may be more important. We measured and modeled carbon fluxes and water relations in first-year Pinus flexilis seedlings in garden plots just beyond the warm edge of their natural range, and compared these with dry-mass gain and survival across two summers. We hypothesized that mortality in these seedlings would be associated with declines in water relations, more so than with carbon-balance limitations. Rather than gradual declines in survivorship across growing seasons, we observed sharp, large-scale mortality episodes that occurred once volumetric soil-moisture content dropped below 10%. By this point, seedling water potentials had decreased below -5 MPa, seedling hydraulic conductivity had decreased by 90% and seedling hydraulic resistance had increased by >900%. Additionally, non-structural carbohydrates accumulated in aboveground tissues at the end of both summers, suggesting impairments in phloem-transport from needles to roots. This resulted in low carbohydrate concentrations in roots, which likely impaired root growth and water uptake at the time of critically low soil moisture. While photosynthesis and respiration on a leaf area basis remained high until critical hydraulic thresholds were exceeded, modeled seedling gross primary productivity declined steadily throughout the summers. At the time of mortality, modeled productivity was insufficient to support seedling biomass-gain rates, metabolism and secondary costs. Thus the large-scale mortality events that we observed near the end of each summer were most directly linked with acute, episodic declines in plant hydraulic function that were linked with important changes in whole-seedling carbon relations.

Molecular markers reveal limited population genetic structure in a North American corvid, Clark's nutcracker (Nucifraga columbiana)

The genetic impact of barriers and Pleistocene glaciations on high latitude resident species has not been widely investigated. The Clark’s nutcracker is an endemic North American corvid closely associated with Pinus-dominated forests. The nutcracker’s encompasses known barriers to dispersal for other species, and glaciated and unglaciated areas. Clark’s nutcrackers also irruptively disperse long distances in search of pine seed crops, creating the potential for gene flow among populations. Using the highly variable mitochondrial DNA control region, seven microsatellite loci, and species distribution modeling, we examined the effects of glaciations and dispersal barriers on population genetic patterns and population structure of nutcrackers. We sequenced 900 bp of mitochondrial control region for 169 individuals from 15 populations and analysed seven polymorphic microsatellite loci for 13 populations across the Clark’s nutcracker range. We used species distribution modeling and a range of phylogeographic analyses to examine evolutionary history. Clark’s nutcracker populations are not highly differentiated throughout their range, suggesting high levels of gene flow among populations, though we did find some evidence of isolation by distance and peripheral isolation. Our analyses suggested expansion from a single refugium after the last glacial maximum, but patterns of genetic diversity and paleodistribution modeling of suitable habitat were inconclusive as to the location of this refugium. Potential barriers to dispersal (e.g. mountain ranges) do not appear to restrict gene flow in Clark’s nutcracker, and postglacial expansion likely occurred quickly from a single refugium located south of the ice sheets.

Deep genetic divergence between disjunct Refugia in the Arctic-Alpine King's Crown, Rhodiola integrifolia (Crassulaceae)

Despite the strength of climatic variability at high latitudes and upper elevations, we still do not fully understand how plants in North America that are distributed between Arctic and alpine areas responded to the environmental changes of the Quaternary. To address this question, we set out to resolve the evolutionary history of the King's Crown, Rhodiola integrifolia using multi-locus population genetic and phylogenetic analyses in combination with ecological niche modeling. Our population genetic analyses of multiple anonymous nuclear loci revealed two major clades within R. integrifolia that diverged from each other ~ 700 kya: one occurring in Beringia to the north (including members of subspecies leedyi and part of subspecies integrifolia), and the other restricted to the Southern Rocky Mountain refugium in the south (including individuals of subspecies neomexicana and part of subspecies integrifolia). Ecological niche models corroborate our hypothesized locations of refugial areas inferred from our phylogeographic analyses and revealed some environmental differences between the regions inhabited by its two subclades. Our study underscores the role of geographic isolation in promoting genetic divergence and the evolution of endemic subspecies in R. integrifolia. Furthermore, our phylogenetic analyses of the plastid spacer region trnL-F demonstrate that among the native North American species, R. integrifolia and R. rhodantha are more closely related to one another than either is to R. rosea. An understanding of these historic processes lies at the heart of making informed management decisions regarding this and other Arctic-alpine species of concern in this increasingly threatened biome.

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.

Colonizing the American continent: Systematics of the genus Arabis in North America (Brassicaceae)

The circumscription of the genus Arabis underwent many and drastic changes within the past. Using DNA sequence information from the nuclear ribosomal RNA and parts of the plastid genome (trnL-trnLF), as well as a critical evaluation of herbarium material from East Asia and North America, we circumscribe the various Arabis taxa of North America. The American and East Asian Arabis species are closely related and, contrary to what was previously believed, they are not closely related to the Eurasian A. hirsuta. Using cpDNA, we found five North American lineages of Arabis with distinct distribution patterns, of which only the purple/red-flowered lineage consists of proven diploids that evolved directly from East Asian progenitors. All other four lineages evolved via ancient hybridization either on the Asian continent prior to migration to North America or showed significant evidence for hybridization and reticulation while diversifying on the American continent. We also provide the first evidence for the systematic circumscription of East Asian Arabis taxa, which together with the North American taxa, form one clade distantly related to European A. ciliata and Eurasian A. hirsuta. The findings also represent the first record of A. pycnocarpa for the floras of China, Japan, and Russian Far East.

Finding a (pine) needle in a haystack: chloroplast genome sequence divergence in rare and widespread pines

Critical to conservation efforts and other investigations at low taxonomic levels, DNA sequence data offer important insights into the distinctiveness, biogeographic partitioning and evolutionary histories of species. The resolving power of DNA sequences is often limited by insufficient variability at the intraspecific level. This is particularly true of studies involving plant organelles, as the conservative mutation rate of chloroplasts and mitochondria makes it difficult to detect polymorphisms necessary to track genealogical relationships among individuals, populations and closely related taxa, through space and time. Massively parallel sequencing (MPS) makes it possible to acquire entire organelle genome sequences to identify cryptic variation that would be difficult to detect otherwise. We are using MPS to evaluate intraspecific chloroplast-level divergence across biogeographic boundaries in narrowly endemic and widespread species of Pinus. We focus on one of the world's rarest pines - Torrey pine (Pinus torreyana) - due to its conservation interest and because it provides a marked contrast to more widespread pine species. Detailed analysis of nearly 90% ( approximately 105 000 bp each) of these chloroplast genomes shows that mainland and island populations of Torrey pine differ at five sites in their plastome, with the differences fixed between populations. This is an exceptionally low level of divergence (1 polymorphism/ approximately 21 kb), yet it is comparable to intraspecific divergence present in widespread pine species and species complexes. Population-level organelle genome sequencing offers new vistas into the timing and magnitude of divergence within species, and is certain to provide greater insight into pollen dispersal, migration patterns and evolutionary dynamics in plants.

SNP, ARMS and SSH authentication of medicinal Dendrobium officinale KIMURA et MIGO and application for identification of Fengdou drugs

Dried stems of Dendrobium officinale have been used as crude drugs in traditional Chinese medicine (TCM) with good tonic efficacy. Sequences of chloroplast, nuclear and mitochondria genes and the method of genomic DNA (gDNA) suppression subtraction hybridization (SSH) were used to authenticate different populations during the process of good agriculture practice (GAP) and crude drug quality control. Six populations could be authenticated successfully by nine single sucleotide polymorphism (SNP) sites and six pairs of diagnostic primers for amplification refractory mutation system (ARMS) were also designed to identify six populations on the basis of single nucleotide polymorphism (SNPs). The remainder two populations (JSR, GGL) with the same sequences could be authenticated by SSH. One population-specific fragment was obtained by SSH and a pair of specific primers (SSH-JB01, SSH-JB02) on the specific sequence was designed to authenticate GGL population from the other populations tested. As the resultants were population-specific, the botanic origins of fifty "Fengdou" drug samples from markets could be classified. It is evident that the combined methods provide a high throughput and reliable approach for identification of D. officinale plants and "Fengdou" drugs.

Mitochondrial and chloroplast phylogeography ofPicea crassifoliaKom. (Pinaceae) in the Qinghai-Tibetan Plateau and adjacent highlands

The disjunct distribution of forests in the Qinghai-Tibetan Plateau (QTP) and adjacent Helan Shan and Daqing Shan highlands provides an excellent model to examine vegetation shifts, glacial refugia and gene flow of key species in this complex landscape region in response to past climatic oscillations and human disturbance. In this study, we examined maternally inherited mitochondrial DNA (nad1 intron b/c and nad5 intron 1) and paternally inherited chloroplast DNA (trnC-trnD) sequence variation within a dominant forest species, Picea crassifolia Kom. We recovered nine mitotypes and two chlorotypes in a survey of 442 individuals from 32 populations sampled throughout the species' range. Significant mitochondrial DNA population subdivision was detected (G(ST) = 0.512; N(ST) = 0.679), suggesting low levels of recurrent gene flow through seeds among populations and significant phylogeographical structure (N(ST) > GST, P < 0.05). Plateau haplotypes differed in sequence from those in the adjacent highlands, suggesting a long period of allopatric fragmentation between the species in the two regions and the presence of independent refugia in each region during Quaternary glaciations. On the QTP platform, all but one of the disjunct populations surveyed were fixed for the same mitotype, while most populations at the plateau edge contained more than one haplotype with the mitotype that was fixed in plateau platform populations always present at high frequency. This distribution pattern suggests that present-day disjunct populations on the QTP platform experienced a common recolonization history. The same phylogeographical pattern, however, was not detected for paternally inherited chloroplast DNA haplotypes. Two chlorotypes were distributed throughout the range of the species with little geographical population differentiation (G(ST) = N(ST) = 0.093). This provides evidence for highly efficient pollen-mediated gene flow among isolated forest patches, both within and between the QTP and adjacent highland populations. A lack of isolation to pollen-mediated gene flow between forests on the QTP and adjacent highlands is surprising given that the Tengger Desert has been a geographical barrier between these two regions for approximately the last 1.8 million years.

Phylogeography of the endangeredCathaya argyrophylla(Pinaceae) inferred from sequence variation of mitochondrial and nuclear DNA

Cathaya argyrophylla is an endangered conifer restricted to subtropical mountains of China. To study phylogeographical pattern and demographic history of C. argyrophylla, species-wide genetic variation was investigated using sequences of maternally inherited mtDNA and biparentally inherited nuclear DNA. Of 15 populations sampled from all four distinct regions, only three mitotypes were detected at two loci, without single region having a mixed composition (G(ST) = 1). Average nucleotide diversity (theta(ws) = 0.0024; pi(s) = 0.0029) across eight nuclear loci is significantly lower than those found for other conifers (theta(ws) = 0.003 approximately 0.015; pi(s) = 0.002 approximately 0.012) based on estimates of multiple loci. Because of its highest diversity among the eight nuclear loci and evolving neutrally, one locus (2009) was further used for phylogeographical studies and eight haplotypes resulting from 12 polymorphic sites were obtained from 98 individuals. All the four distinct regions had at least four haplotypes, with the Dalou region (DL) having the highest diversity and the Bamian region (BM) the lowest, paralleling the result of the eight nuclear loci. An AMOVA revealed significant proportion of diversity attributable to differences among regions (13.4%) and among populations within regions (8.9%). F(ST) analysis also indicated significantly high differentiation among populations (F(ST) = 0.22) and between regions (F(ST) = 0.12-0.38). Non-overlapping distribution of mitotypes and high genetic differentiation among the distinct geographical groups suggest the existence of at least four separate glacial refugia. Based on network and mismatch distribution analyses, we do not find evidence of long distance dispersal and population expansion in C. argyrophylla. Ex situ conservation and artificial crossing are recommended for the management of this endangered species.

Contrasting phylogeographical patterns between mainland and island taxa of the Pinus luchuensis complex

Species whose geographical distribution encompasses both mainland and island populations provide an ideal system for examining isolation and genetic divergence. In this study, paternally transmitted chloroplast DNA (cpDNA) and maternally transmitted mitochondrial DNA (mtDNA) were used to estimate population structure and phylogeography of Pinus luchuensis, a species found in eastern China (ssp. hwangshanensis), Taiwan (ssp. taiwanensis), and the Ryukyu Archipelago (ssp. luchuensis). Gene genealogies of both mtDNA and cpDNA reveal two major lineages. Molecular dating indicates that these lineages diverged before the colonization of P. luchuensis subspecies in Taiwan and the Ryukyu Archipelago. Both mtDNA and cpDNA show a lack of correspondence between molecular phylogeny and subspecies designation. Phylogeographical analysis suggests that paraphyly of the subspecies is the result of recent divergence rather than secondary contacts. In spite of the short divergence history of P. luchuensis on islands, the island populations show the same degree of genetic divergence as mainland populations. Low levels of genetic diversity in the mainland ssp. hwangshanensis suggest demographic bottlenecks. In contrast, the high heterogeneity of genetic composition for island populations is likely to be associated with a history of multiple colonization from the mainland. The spatial apportionment of organelle DNA polymorphisms is consistent with a pattern of stepwise colonization on island populations.

Geographic patterns of microsatellite variation in Boechera stricta, a close relative of Arabidopsis

The genus Boechera is a widespread North American group with great potential for studies of ecology and evolution: Boechera is closely related to Arabidopsis and exhibits different ecological and reproductive strategies. Boechera stricta (previously Arabis drummondii) is a morphologically and genetically well-defined, perennial crucifer species. Fifteen natural populations of diploid individuals from the Rocky Mountains were analysed using 21 microsatellite loci. In accordance with our expectation for this predominately inbreeding species, a high F IS value (0.89) was observed. Furthermore, populations of B. stricta were highly differentiated, as indicated by F ST = 0.56. Three clusters were identified using structure- the majority of populations belonged to either the Northern or Southern cluster. Together, the north-south partitioning and evenness of genetic variation across the two clusters suggested multiple refugia for this perennial herb in the Rocky Mountains. Pleistocene glaciation, together with the topographically and climatologically heterogeneous cordillera, has profoundly influenced the genetic architecture of B. stricta. Genetic population structure was also influenced by relatively recent genome admixture at two levels: within species (involving individuals from the Northern and Southern clusters) and between species (with the hybridization of B. stricta and Boechera holboellii). This complexity of population structure at presumably neutral microsatellite loci located throughout the genome in B. stricta provides a baseline against which to test whether functional genetic variation is undergoing local adaptive evolution throughout the natural species range.

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.

Mitochondrial genome recombination in the zone of contact between two hybridizing conifers

Variation in mitochondrial DNA was surveyed at four gene loci in and around the zone of contact between two naturally hybridizing conifers, black spruce (Picea mariana) and red spruce (P. rubens) in northeastern North America. Most of the mtDNA diversity of these species was found in populations next to or into the zone of contact, where some individuals bore rare mitotypes intermediate between the common mitotypes observed in the allopatric areas of each species. Sequence analysis and tests for mtDNA recombination point to this phenomenon, rather than to recurrent mutation, as the most tenable hypothesis for the origin of these rare mitotypes. From the 10 mitotypes observed, at least 4 would be the product of recombination between 4 of the 5 putative ancestral mitotypes. Tests for cytonuclear disequilibrium and geographical structure of the putative recombinant mitotypes suggest that mtDNA recombination is not frequent and relatively recent on the geological time scale. mtDNA recombination would have been promoted by transient heteroplasmy due to leakage of paternal mtDNA since the Holocene secondary contact between the two species.

Variation in mitochondrial DNA reveals multiple distant glacial refugia in black spruce (Picea mariana), a transcontinental North American conifer

Range-wide genetic variation of black spruce (Picea mariana) was studied using polymerase chain reaction-random fragment length polymorphism markers of the mitochondrial genome. Four polymorphic mitochondrial DNA (mtDNA) loci were surveyed and two or three alleles were detected at each locus, resulting in 10 multilocus mtDNA types or mitotypes. A significant subdivision of population genetic diversity was detected (GST = 0.671; NST = 0.726), suggesting low levels of gene flow among populations. The distribution of mitotypes was not random (NST > GST; P < 0.05) and revealed four partially overlapping zones, presumably representative of different glacial populations. Comparison of the genetic structure derived from mtDNA markers and the colonization paths previously deduced from the fossil and pollen records allow us to infer at least three southern and one northeastern glacial populations for black spruce. The patterns revealed in this study suggest that black spruce shares its biogeographical history with other forest-associated North American species.

Genetic diversity of nuclear and mitochondrial genomes in Pinus parviflora Sieb. & Zucc. (Pinaceae) populations

Genetic diversities of the nuclear and mitochondrial genomes in Pinus parviflora were studied in 16 populations, which were distributed across most of the species' range in Japan. Six mitochondrial DNA haplotypes were identified among the 16 populations. The intrapopulation diversity of allozymes was similar to that of other endemic woody species (H(S)=0.259). Although P. parviflora is distributed in discrete populations, differentiation between these was very low (G(ST)=0.044). In addition, the extent of genetic differentiation between two varieties (var. pentaphylla and var. parviflora) was extremely low (G(VT)=0.001). Intrapopulation diversity of mitochondrial DNA was also very low (H(S)=0.098), but population differentiation was high (G(ST)=0.863). Moreover, the distribution of haplotypes reflected the taxonomic differences between P. parviflora var. pentaphylla and var. parviflora. The populations of var. pentaphylla and var. parviflora contained different haplotypes. Differing modes of inheritance may account for the differences in nuclear and mitochondrial genetic diversity.

Cross-species amplification of mitochondrial DNA sequence-tagged-site markers in conifers: the nature of polymorphism and variation within and among species in Picea

Primers previously developed to amplify specific non-coding regions of the mitochondrial genome in Angiosperms, and new primers for additional non-coding mtDNA regions, were tested for their ability to direct DNA amplification in 12 conifer taxa and to detect sequence-tagged-site (STS) polymorphisms within and among eight species in Picea. Out of 12 primer pairs, nine were successful at amplifying mtDNA in most of the taxa surveyed. In conifers, indels and substitutions were observed for several loci, allowing them to distinguish between families, genera and, in some cases, between species within genera. In Picea, interspecific polymorphism was detected for four loci, while intraspecific variation was observed for three of the mtDNA regions studied. One of these (SSU rRNA V1 region) exhibited indel polymorphisms, and the two others ( nad1 intron b/c and nad5 intron1) revealed restriction differences after digestion with Sau3AI (PCR-RFLP). A fourth locus, the nad4L- orf25 intergenic region, showed a multibanding pattern for most of the spruce species, suggesting a possible gene duplication. Maternal inheritance, expected for mtDNA in conifers, was observed for all polymorphic markers except the intergenic region nad4L- orf25. Pooling of the variation observed with the remaining three markers resulted in two to six different mtDNA haplotypes within the different species of Picea. Evidence for intra-genomic recombination was observed in at least two taxa. Thus, these mitotypes are likely to be more informative than single-locus haplotypes. They should be particularly useful for the study of biogeography and the dynamics of hybrid zones.

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.

Land barriers and open oceans: effects on gene diversity and population structure in Avicennia germinans L. (Avicenniaceae)

Avicennia germinans L. is a widespread mangrove species occupying the west coast of Africa and the Atlantic and Pacific coasts of the Americas from the Bahamas to Brazil and Baja California to Peru. An amplified fragment length polymorphism (AFLP) molecular analysis was carried out to assess genetic architecture within this species and to evaluate the effects of the Atlantic Ocean and the Central American Isthmus (CAI) on population and regional genetic diversity and differentiation. In total, 349 polymorphic AFLP fragments were identified among 144 individuals from 14 populations from the east Atlantic, west Atlantic and east Pacific. Levels of genetic diversity varied considerably among populations, but were generally higher in populations from the east Atlantic. Regional differentiation between the Pacific coast and Atlantic populations was greater than between east and west Atlantic populations, suggesting that the CAI has had an important influence on population genetic structure in this species. The lower level of divergence of east Atlantic from west Atlantic populations suggests some dispersal across the Atlantic Ocean, although migration rates are probably low; Nm from GST equal to 0.41 and accumulation of private and rare alleles in the east Atlantic. Population differentiation did not appear to follow an isolation by distance model and has probably resulted from complex patterns of population bottlenecks, and founder events due to landscape changes during the Pleistocene, particularly in the west Atlantic. The molecular data provide no support for the treatment of east Atlantic populations as a separate species A. africana.

Mitochondrial DNA diversity, population structure, and gender association in the gynodioecious plant Silene vulgaris

A highly variable mitochondrial DNA (mtDNA) restriction fragment length polymorphism (RFLP) locus is used to assess the population structure of mitochondrial genomes in the gynodioecious plant Silene vulgaris at two spatial scales. Thirteen mtDNA haplotypes were identified within 250 individuals from 18 populations in a 20-km diameter region of western Virginia. The population structure of these mtDNA haplotypes was estimated as thetaST = 0.574 (+/- 0.066 SE) and, surprisingly, genetic differentiation among populations was negatively correlated with geographic distance (Mantel r = -0.246, P < 0.002). Additionally, mtDNA haplotypes were spatially clumped at the scale of meters within one population. Gender in S. vulgaris is determined by an interaction between autosomal male fertility restorers and cytoplasmic male sterility (CMS) factors, and seed fitness is affected by an interaction between gender and population sex ratio; thus, selection acting on gender could influence the distribution of mtDNA RFLP haplotypes. The sex ratio (females:hermaphrodites) varied among mtDNA haplotypes across the entire metapopulation, possibly because the haplotypes were in linkage disequilibrium with different CMS factors. The gender associated with some of the most common haplotypes varied among populations, suggesting that there is also population structure in male fertility restorer genes. In comparison with reports of mtDNA variation from other published studies, we found that S. vulgaris exhibits a large number of mtDNA haplotypes relative to that observed in other species.

DNA from bird-dispersed seed and wind-disseminated pollen provides insights into postglacial colonization and population genetic structure of whitebark pine (Pinus albicaulis)

Uniparentally inherited mitochondrial (mt)DNA and chloroplast (cp)DNA microsatellites (cpSSRs) were used to examine population genetic structure and biogeographic patterns of bird-dispersed seed and wind-disseminated pollen of whitebark pine (Pinus albicaulis Engelm.). Sampling was conducted from 41 populations throughout the range of the species. Analyses provide evidence for an ancestral haplotype and two derived mtDNA haplotypes with distinct regional distributions. An abrupt contact zone between mtDNA haplotypes in the Cascade Range suggests postglacial biogeographic movements. Among three cpSSR loci, 42 haplotypes were detected within 28 cpSSR sample populations that were aggregated into six regions. Analysis of molecular variance (amova) was used to determine the hierarchical genetic structure of cpSSRs. amova and population pairwise comparisons (FST ) of cpSSR, and geographical distribution of mtDNA haplotypes provide insights into historical changes in biogeography. The genetic data suggest that whitebark pine has been intimately tied to climatic change and associated glaciation, which has led to range movements facilitated by seed dispersal by Clark's nutcracker (Nucifraga columbiana Wilson). The two hypotheses proposed to explain the genetic structure are: (i) a northward expansion into Canada and the northern Cascades in the early Holocene; and (ii) historical gene flow between Idaho and the Oregon Cascades when more continuous habitat existed in Central Oregon during the late Pleistocene. Genetic structure and insights gained from historical seed movements provide a basis on which to develop recovery plans for a species that is at risk from multiple threats.

Phylogenetic relationships among domesticated and wild species of Cucurbita (Cucurbitaceae) inferred from a mitochondrial gene: Implications for crop plant evolution and areas of origin

We have investigated the phylogenetic relationships among six wild and six domesticated taxa of Cucurbita using as a marker an intron region from the mitochondrial nad1 gene. Our study represents one of the first successful uses of a mtDNA gene in resolving inter- and intraspecific taxonomic relationships in Angiosperms and yields several important insights into the origins of domesticated Cucurbita. First, our data suggest at least six independent domestication events from distinct wild ancestors. Second, Cucurbita argyrosperma likely was domesticated from a wild Mexican gourd, Cucurbita sororia, probably in the same region of southwest Mexico that gave rise to maize. Third, the wild ancestor of Cucurbita moschata is still unknown, but mtDNA data combined with other sources of information suggest that it will probably be found in lowland northern South America. Fourth, Cucurbita andreana is supported as the wild progenitor of Cucurbita maxima, but humid lowland regions of Bolivia in addition to warmer temperate zones in South America from where C. andreana was originally described should possibly be considered as an area of origin for C. maxima. Fifth, our data support other molecular results that indicate two separate domestications in the Cucurbita pepo complex. The potential zone of domestication for one of the domesticated subspecies, C. pepo subsp. ovifera, includes eastern North America and should be extended to northeastern Mexico. The wild ancestor of the other domesticated subspecies, C. pepo subsp. pepo, is undiscovered but is closely related to C. pepo subsp. fraterna and possibly will be found in southern Mexico.

The evolutionary split of Pinaceae from other conifers: evidence from an intron loss and a multigene phylogeny

The second intron in the mitochondrial gene nad1 was surveyed using PCR, DNA sequencing, or Southern hybridization in 323 species (313 genera, 212 families) of seed plants. The intron was absent in all 22 species (22 genera, 8 families) of non-Pinaceae conifers studied, in Welwitschia mirabilis, and in seven angiosperms. Whereas absence of the intron in seven angiosperms and Welwitschia is likely due to seven independent losses when evaluated against the recently published multigene phylogenies, the lack of the intron in all non-Pinaceae conifers can be best explained by a single loss. These data suggest that the non-Pinaceae conifers represent a monophyletic group. We also conducted a phylogenetic analysis of seed plants using a combined data set of the partial exon and intron sequences of nad1 generated from this study and published sequences of mitochondrial cox1 and small subunit (SSU) rDNA, chloroplast rbcL, and nuclear 18S rDNA. The results supported the split of conifers into two groups: Pinaceae and non-Pinaceae conifers. The Gnetales were sister to Pinaceae, in agreement with the conclusion from other recent molecular phylogenetic studies that refute the anthophyte hypothesis.

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.

Did the pleistocene glaciations promote divergence? Tests of explicit refugial models in montane grasshopprers

There is a long-standing debate over whether or not the Pleistocene glaciations promoted speciation. While some models predict that extensive mixing of populations during interglacial expansion would have inhibited divergence, others postulate that divergence among allopatric glacial refuges or founder events during recolonization of previously glaciated areas would have promoted differentiation. Using a combination of traditional and coalescent based population genetic approaches, this study finds that the glaciations did not inhibit divergence among populations of the grasshopper Melanoplus oregonensis. Instead, drift associated with recolonization of previously glaciated areas, as well as divergence among multiple allopatric glacial refugia, have both contributed to differentiation in this montane grasshopper from the 'sky islands' of the northern Rocky Mountains. Significant population structure was detected by phylogenetic and FST analyses, including significant FST values among individual pairs of sky-island populations. In addition to clustering of haplotypes within populations, there is some evidence of regional phylogeographic structure, although none of the 'regional groups' form a monophyletic clade and there is a lack of concordance between the genealogical and geographical positions of some haplotypes. However, coalescent simulations confirm there is significant regional phylogeographic structure that most likely reflects divergence among multiple ancestral refugial populations, and indicate that it is very unlikely that the observed gene tree could have been produced by the fragmentation of a single widespread ancestral population. Thus, rather than inhibiting differentiation, the glaciations appear to have promoted population divergence in M. oregonensis, suggesting that they may have contributed to the radiation of Melanoplus species during the Pleistocene.

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.