Glacial vicariance in the Pacific Northwest: evidence from a lodgepole pine mitochondrial DNA minisatellite for multiple genetically distinct and widely separated refugia

Characterizing the elusive Vancouver Island wolverine, Gulo gulo vancouverensis, using historical DNA

The wolverine (Gulo gulo) is a Holarctic species found in North America primarily across the boreal forest, the subarctic, and along the Pacific coast, including Vancouver Island (VI), British Columbia. While wolverines on VI are rare and possibly extirpated, they have been previously described as a unique subspecies, G. g. vancouverensis, distinct from G. g. luscus from the mainland of North America. However, the validity of the VI subspecies is contentious, with conflicting results from studies of skull morphology. Here, we used molecular analyses to characterize the genetic diversity of the VI population and resolve this taxonomic debate to assist with conservation priorities. Historical DNA of VI wolverines was obtained from museum specimens, amplified at 16 nuclear microsatellite loci, and sequenced at the mitochondrial D-loop control region to compare with wolverines from mainland British Columbia. The VI population had lower allelic richness and was fixed for a single common mtDNA haplotype. Bayesian and non-Bayesian assignments using microsatellites generally revealed admixture across populations, implying allele frequencies between the VI and mainland populations were not significantly different. Hence, both types of genetic markers showed little evolutionary divergence between VI and the mainland population. Combined, these results do not provide evidence of significant genetic distinction for VI wolverines, nor support the subspecific classification. Immediate conservation efforts should focus on estimating population size, while future conservation planning can assume VI wolverines likely are not a unique genetic population and there remains the potential for natural recolonization of wolverines to VI.

Molecular phylogeny and ancestral biogeographic reconstruction of Platanthera subgenus Limnorchis (Orchidaceae) using target capture methods

Platanthera is one of the largest genera of temperate orchids in the Holarctic and exemplifies a lineage that has adaptively radiated into diverse habitats within North America, Asia, Europe, North Africa, Borneo, and Sarawak. Major centers of diversity in this genus are North America and eastern Asia. Despite its diversity, a thorough phylogenetic hypothesis for the genus is lacking because no studies have yet sampled taxa exhaustively or developed a robust molecular toolkit. While there is strong evidence that suggests monophyly of subgenus Limnorchis, most taxa in this group have not been included in a phylogenetic analysis. In this study, we developed a new toolkit for Platanthera consisting of genomic information from 617 low-copy nuclear loci. Using a targeted enrichment approach, we collected high-throughput sequence data in 23 accessions of nine of the 12 diploid species of subgenus Limnorchis and outgroup species across Platanthera. A maximum likelihood analysis resolved a strongly supported monophyletic clade for subgenus Limnorchis. Ancestral biogeographic reconstruction indicated that subgenus Limnorchis originated in western North America ca. 3-4.5 Mya from an ancestor that was widespread in western North America and eastern Asia and subsequently diversified in western North America, followed by dispersal of some species to eastern North America. Our results indicate complex biogeographic connections between Asia and North America, and therefore it suggests that Platanthera is a suitable system to test biogeographic hypotheses over time and space in the Holarctic. Our results are also expected to facilitate further study of diversification and biogeographic spread across Platanthera and lay the groundwork for understanding independent origins, biogeography, and morphological diversification of polyploid species within subgenus Limnorchis.

Signatures of Post-Glacial Genetic Isolation and Human-Driven Migration in the Dothistroma Needle Blight Pathogen in Western Canada

Many current tree improvement programs are incorporating assisted gene flow strategies to match reforestation efforts with future climates. This is the case for the lodgepole pine (Pinus contorta var. latifolia), the most extensively planted tree in western Canada. Knowledge of the structure and origin of pathogen populations associated with this tree would help improve the breeding effort. Recent outbreaks of the Dothistroma needle blight (DNB) pathogen Dothistroma septosporum on lodgepole pine in British Columbia and its discovery in Alberta plantations raised questions about the diversity and population structure of this pathogen in western Canada. Using genotyping-by-sequencing on 119 D. septosporum isolates from 16 natural pine populations and plantations from this area, we identified four genetic lineages, all distinct from the other DNB lineages from outside of North America. Modeling of the population history indicated that these lineages diverged between 31.4 and 7.2 thousand years ago, coinciding with the last glacial maximum and the postglacial recolonization of lodgepole pine in western North America. The lineage found in the Kispiox Valley from British Columbia, where an unprecedented DNB epidemic occurred in the 1990s, was close to demographic equilibrium and displayed a high level of haplotypic diversity. Two lineages found in Alberta and Prince George (British Columbia) showed departure from random mating and contemporary gene flow, likely resulting from pine breeding activities and material exchanges in these areas. The increased movement of planting material could have some major consequences by facilitating secondary contact between genetically isolated DNB lineages, possibly resulting in new epidemics.

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.

Cryptic and extensive hybridization between ancient lineages of American crows

Most species and therefore most hybrid zones have historically been defined using phenotypic characters. However, both speciation and hybridization can occur with negligible morphological differentiation. Recently developed genomic tools provide the means to better understand cryptic speciation and hybridization. The Northwestern Crow (Corvus caurinus) and American Crow (Corvus brachyrhynchos) are continuously distributed sister taxa that lack reliable traditional characters for identification. In this first population genomic study of Northwestern and American crows, we use genomic SNPs (nuDNA) and mtDNA to investigate the degree of genetic differentiation between these crows and the extent to which they may hybridize. Our results indicate that American and Northwestern crows have distinct evolutionary histories, supported by two nuDNA ancestry clusters and two 1.1%-divergent mtDNA clades dating to the late Pleistocene, when glacial advances may have isolated crow populations in separate refugia. We document extensive hybridization, with geographic overlap of mtDNA clades and admixture of nuDNA across >900 km of western Washington and western British Columbia. This broad hybrid zone consists of late-generation hybrids and backcrosses, but not recent (e.g., F1) hybrids. Nuclear DNA and mtDNA clines had concordant widths and were both centred in southwestern British Columbia, farther north than previously postulated. Overall, our results suggest a history of reticulate evolution in American and Northwestern crows, perhaps due to recurring neutral expansion(s) from Pleistocene glacial refugia followed by lineage fusion(s). However, we do not rule out a contributing role for more recent potential drivers of hybridization, such as expansion into human-modified habitats.

Living on the edge: comparative phylogeography and phylogenetics of Oreohelix land snails at their range edge in Western Canada

BACKGROUND

The biodiversity and distributions of terrestrial snails at local and regional scales are influenced by their low vagility and microhabitat specificity. The accessibility of large-bodied species and their characteristically high levels of genetic polymorphism make them excellent ecological and evolutionary models for studies on the phylogeography, phylogenetics, and conservation of organisms in fragmented populations. This study aims to elucidate the biodiversity, systematics, and distributions of genetic lineages within the genus Oreohelix at the northern and western periphery of their range.

RESULTS

We found four mitochondrial clades, three of which are putative subspecies of Oreohelix subrudis. One clade was geographically widespread, occurring within numerous sites in Cypress Hills and in the Rocky Mountains, a second was geographically restricted to the Rocky Mountains in Alberta, and a third was restricted to the Cypress Hills region. A fourth clade was the small-bodied species, O. cooperi. ITS2 sequence and screening data revealed three genetic clusters, of which one was O. cooperi. Cluster 1 contained most individuals in COI clade X and some from clade B and cluster 2 was predominantly made up of individuals from COI clades B and B' and a few from clade X. ITS2 alleles were shared in a narrow contact zone between two COI clades, suggestive of hybridization between the two.

CONCLUSIONS

A sky island known as Cypress Hills, in southeastern Alberta, Canada, is a biodiversity hotspot for terrestrial land snails in the genus Oreohelix. The observed phylogeographic patterns likely reflect reproductive isolation during the Last Glacial Maximum, followed by secondary contact due to passive, long-range dispersal resulting from low vagility, local adaptation, and complex glacial history.

Purifying selection does not drive signatures of convergent local adaptation of lodgepole pine and interior spruce

BACKGROUND

Lodgepole pine (Pinus contorta) and interior spruce (Picea glauca, Picea engelmannii, and their hybrids) are distantly related conifer species. Previous studies identified 47 genes containing variants associated with environmental variables in both species, providing evidence of convergent local adaptation. However, if the intensity of purifying selection varies with the environment, clines in nucleotide diversity could evolve through linked (background) selection that would yield allele frequency-environment signatures resembling local adaptation. If similar geographic patterns in the strength of purifying selection occur in these species, this could result in the convergent signatures of local adaptation, especially if the landscape of recombination is conserved. In the present study, we investigated whether spatially/environmentally varying purifying selection could give rise to the convergent signatures of local adaptation that had previously reported.

RESULTS

We analyzed 86 lodgepole pine and 50 interior spruce natural populations spanning heterogeneous environments in western Canada where previous analyses had found signatures of convergent local adaptation. We estimated nucleotide diversity and Tajima's D for each gene within each population and calculated the strength of correlations between nucleotide diversity and environmental variables. Overall, these estimates in the genes with previously identified convergent local adaptation signatures had no similar pattern between pine and spruce. Clines in nucleotide diversity along environmental variables were found for interior spruce, but not for lodgepole pine. In spruce, genes with convergent adaption signatures showed a higher strength of correlations than genes without convergent adaption signatures, but there was no such disparity in pine, which suggests the pattern in spruce may have arisen due to a combination of selection and hybridization.

CONCLUSIONS

The results rule out purifying/background selection as a driver of convergent local adaption signatures in lodgepole pine and interior spruce.

Northern forest tree populations are physiologically maladapted to drought

Northern forests at the leading edge of their distributions may not show increased primary productivity under climate warming, being limited by climatic extremes such as drought. Looking beyond tree growth to underlying physiological mechanisms is fundamental for accurate predictions of forest responses to climate warming and drought stress. Within a 32-year genetic field trial, we analyze relative contributions of xylem plasticity and inferred stomatal response to drought tolerance in regional populations of a widespread conifer. Genetic adaptation leads to varying responses under drought. Trailing-edge tree populations produce fewer tracheids with thicker cell walls, characteristic of drought-tolerance. Stomatal response explains the moderate drought tolerance of tree populations in central areas of the species range. Growth loss of the northern population is linked to low stomatal responsiveness combined with the production of tracheids with thinner cell walls. Forests of the western boreal may therefore lack physiological adaptations necessary to tolerate drier conditions.

Reproductive isolation and environmental adaptation shape the phylogeography of mountain pine beetle (Dendroctonus ponderosae)

Chromosomal rearrangement can be an important mechanism driving population differentiation and incipient speciation. In the mountain pine beetle (MPB, Dendroctonus ponderosae), deletions on the Y chromosome that are polymorphic among populations are associated with reproductive incompatibility. Here, we used RAD sequencing across the entire MPB range in western North America to reveal the extent of the phylogeographic differences between Y haplotypes compared to autosomal and X-linked loci. Clustering and geneflow analyses revealed three distinct Y haplogroups geographically positioned within and on either side of the Great Basin Desert. Despite close geographic proximity between populations on the boundaries of each Y haplogroup, there was extremely low Y haplogroup mixing among populations, and gene flow on the autosomes was reduced across Y haplogroup boundaries. These results are consistent with a previous study suggesting that independent degradation of a recently evolved neo-Y chromosome in previously isolated populations causes male sterility or inviability among Y haplotype lineages. Phylogeographic results supported historic contraction of MPB into three separate Pleistocene glacial refugia followed by postglacial range expansion and secondary contact. Distinct sets of SNPs were statistically associated with environmental data among the most genetically distinct sets of geographic populations. This finding suggests that the process of adaptation to local climatic conditions is influenced by population genetic structure, with evidence for largely independent evolution in the most genetically isolated Y haplogroup.

Multilocus genetic analyses and spatial modeling reveal complex population structure and history in a widespread resident North American passerine (Perisoreus canadensis)

An increasing body of studies of widely distributed, high latitude species shows a variety of refugial locations and population genetic patterns. We examined the effects of glaciations and dispersal barriers on the population genetic patterns of a widely distributed, high latitude, resident corvid, the gray jay (Perisoreus canadensis), using the highly variable mitochondrial DNA (mtDNA) control region and microsatellite markers combined with species distribution modeling. We sequenced 914 bp of mtDNA control region for 375 individuals from 37 populations and screened seven loci for 402 individuals from 27 populations across the gray jay range. We used species distribution modeling and a range of phylogeographic analyses (haplotype diversity, Φ_ST, SAMOVA, F_ST, Bayesian clustering analyses) to examine evolutionary history and population genetic structure. MtDNA and microsatellite markers revealed significant genetic differentiation among populations with high concordance between markers. Paleodistribution models supported at least five potential areas of suitable gray jay habitat during the last glacial maximum and revealed distributions similar to the gray jay's contemporary during the last interglacial. Colonization from and prolonged isolation in multiple refugia is evident. Historical climatic fluctuations, the presence of multiple dispersal barriers, and highly restricted gene flow appear to be responsible for strong genetic diversification and differentiation in gray jays.

Pest risk assessment of Atropellis spp. for the EU territory

Following a request from the European Commission, the EFSA Plant Health (PLH) Panel performed a risk assessment for Atropellis spp. in the EU focusing on the risk of entry, the host range and the potential impacts. Atropellis is a fungal pathogen of several Pinus spp. in North America. The pathogen has not been reported from Europe and is a quarantine pest regulated in Annex IIAI of Council Directive 2000/29/EC on plants (other than fruit and seeds), isolated bark and wood of Pinus. The main pathways of entry considered were Pinus plants, wood and isolated bark. Given the ban of importing Pinus plants from outside Europe into the EU and the lack of information on EU imports of isolated Pinus bark, only the wood pathway was assessed quantitatively. The conclusion of the assessment of entry for scenario A0 (current regulatory situation) is that the risk of entry of Atropellis spp. is close to zero. This conclusion is expected to apply also in the case of removing the specific Atropellis regulations, because of the remaining generic Pinus requirements, as well as in a scenario with additional risk reduction options. The uncertainty associated with this assessment is relatively limited, given that all the quartiles of the estimated distribution of the number of potential founder populations are close to zero. For the North American Pinus spp. known to be susceptible and widely planted in the EU (mainly P. contorta and P. strobus), the damage observed in North America (loss of wood quality, stem deformations, mortality in young stands, environmental consequences) is expected in the EU to a similar (or higher) degree, should the pathogen be introduced. Similar impacts are expected on the European Pinus spp. known to be host of Atropellis spp. These include widespread and locally abundant species such as P. nigra, P. sylvestris and P. pinaster. There are, however, large uncertainties associated with this impact assessment due to the unknown susceptibility of several other Pinus spp. present in Europe. There is a need for research on the susceptibility to Atropellis spp. of those European Pinus spp.

Post-glacial biogeography of trembling aspen inferred from habitat models and genetic variance in quantitative traits

Using species distribution models and information on genetic structure and within-population variance observed in a series of common garden trials, we reconstructed a historical biogeography of trembling aspen in North America. We used an ensemble classifier modelling approach (RandomForest) to reconstruct palaeoclimatic habitat for the periods 21,000, 14,000, 11,000 and 6,000 years before present. Genetic structure and diversity in quantitative traits was evaluated in common garden trials with 43 aspen collections ranging from Minnesota to northern British Columbia. Our main goals were to examine potential recolonisation routes for aspen from southwestern, eastern and Beringian glacial refugia. We further examined if any refugium had stable habitat conditions where aspen clones may have survived multiple glaciations. Our palaeoclimatic habitat reconstructions indicate that aspen may have recolonised boreal Canada and Alaska from refugia in the eastern United States, with separate southwestern refugia for the Rocky Mountain regions. This is further supported by a southeast to northwest gradient of decreasing genetic variance in quantitative traits, a likely result of repeated founder effects. Stable habitat where aspen clones may have survived multiple glaciations was predicted in Mexico and the eastern United States, but not in the west where some of the largest aspen clones have been documented.

Defense traits in the long-lived Great Basin bristlecone pine and resistance to the native herbivore mountain pine beetle

Mountain pine beetle (MPB, Dendroctonus ponderosae) is a significant mortality agent of Pinus, and climate-driven range expansion is occurring. Pinus defenses in recently invaded areas, including high elevations, are predicted to be lower than in areas with longer term MPB presence. MPB was recently observed in high-elevation forests of the Great Basin (GB) region, North America. Defense and susceptibility in two long-lived species, GB bristlecone pine (Pinus longaeva) and foxtail pine (P. balfouriana), are unclear, although they are sympatric with a common MPB host, limber pine (P. flexilis). We surveyed stands with sympatric GB bristlecone-limber pine and foxtail-limber pine to determine relative MPB attack susceptibility and constitutive defenses. MPB-caused mortality was extensive in limber, low in foxtail and absent in GB bristlecone pine. Defense traits, including constitutive monoterpenes, resin ducts and wood density, were higher in GB bristlecone and foxtail than in limber pine. GB bristlecone and foxtail pines have relatively high levels of constitutive defenses which make them less vulnerable to climate-driven MPB range expansion relative to other high-elevation pines. Long-term selective herbivore pressure and exaptation of traits for tree longevity are potential explanations, highlighting the complexity of predicting plant-insect interactions under climate change.

Phylogeographic structure in long-tailed voles (Rodentia: Arvicolinae) belies the complex Pleistocene history of isolation, divergence, and recolonization of Northwest North America's fauna

Quaternary climate fluctuations restructured biodiversity across North American high latitudes through repeated episodes of range contraction, population isolation and divergence, and subsequent expansion. Identifying how species responded to changing environmental conditions not only allows us to explore the mode and tempo of evolution in northern taxa, but also provides a basis for forecasting future biotic response across the highly variable topography of western North America. Using a multilocus approach under a Bayesian coalescent framework, we investigated the phylogeography of a wide‐ranging mammal, the long‐tailed vole, Microtus longicaudus. We focused on populations along the North Pacific Coast to refine our understanding of diversification by exploring the potentially compounding roles of multiple glacial refugia and more recent fragmentation of an extensive coastal archipelago. Through a combination of genetic data and species distribution models (SDMs), we found that historical climate variability influenced contemporary genetic structure, with multiple isolated locations of persistence (refugia) producing multiple divergent lineages (Beringian or northern, southeast Alaska or coastal, and southern or continental) during glacial advances. These vole lineages all occur along the North Pacific Coast where the confluence of numerous independent lineages in other species has produced overlapping zones of secondary contact, collectively a suture zone. Finally, we detected high levels of neoendemism due to complex island geography that developed in the last 10,000 years with the rising sea levels of the Holocene.

Glacial refugia and modern genetic diversity of 22 western North American tree species

North American tree species, subspecies and genetic varieties have primarily evolved in a landscape of extensive continental ice and restricted temperate climate environments. Here, we reconstruct the refugial history of western North American trees since the last glacial maximum using species distribution models, validated against 3571 palaeoecological records. We investigate how modern subspecies structure and genetic diversity corresponds to modelled glacial refugia, based on a meta-analysis of allelic richness and expected heterozygosity for 473 populations of 22 tree species. We find that species with strong genetic differentiation into subspecies had widespread and large glacial refugia, whereas species with restricted refugia show no differentiation among populations and little genetic diversity, despite being common over a wide range of environments today. In addition, a strong relationship between allelic richness and the size of modelled glacial refugia (r(2) = 0.55) suggest that population bottlenecks during glacial periods had a pronounced effect on the presence of rare alleles.

Contrasting patterns of genetic diversity across the ranges of Pinus monticola and P. strobus: a comparison between eastern and western North American postglacial colonization histories

• Premises of the study: Understanding the influence of recent glacial and postglacial periods on species' distributions is key for predicting the effects of future environmental changes. We investigated the influence of two physiographic landscapes on population structure and postglacial colonization of two white pine species of contrasting habitats: P. monticola, which occurs in the highly mountainous region of western North America, and P. strobus, which occurs in a much less mountainous area in eastern North America.•

METHODS

To characterize the patterns of genetic diversity and population structure across the ranges of both species, 158 and 153 single nucleotide polymorphism (SNP) markers derived from expressed genes were genotyped on range-wide samples of 61 P. monticola and 133 P. strobus populations, respectively.•

KEY RESULTS

In P. monticola, a steep latitudinal decrease in genetic diversity likely resulted from postglacial colonization involving rare long-distance dispersal (LDD) events. In contrast, no geographic patterns of diversity were detected in P. strobus, suggesting recolonization via a gradually advancing front or frequent LDD events. For each species, structure analyses identified two distinct southern and northern genetic groups that likely originated from two different glacial lineages. At a finer scale, and for the two species, smaller subgroups were detected that could be remnants of cryptic refugia.•

CONCLUSION

During postglacial colonization, the western and eastern North American landscapes had different impacts on genetic signatures in P. monticola compared with P. strobus. We discuss the importance of our findings for conservation programs and predictions of species' response to climate change.

Beringian sub-refugia revealed in blackfish (Dallia): implications for understanding the effects of Pleistocene glaciations on Beringian taxa and other Arctic aquatic fauna

Background

Pleistocene climatic instability had profound and diverse effects on the distribution and abundance of Arctic organisms revealed by variation in phylogeographic patterns documented in extant Arctic populations. To better understand the effects of geography and paleoclimate on Beringian freshwater fishes, we examined genetic variability in the genus Dallia (blackfish: Esociformes: Esocidae). The genus Dallia groups between one and three nominal species of small, cold- and hypoxia-tolerant freshwater fishes restricted entirely in distribution to Beringia from the Yukon River basin near Fairbanks, Alaska westward including the Kuskokwim River basin and low-lying areas of Western Alaska to the Amguema River on the north side of the Chukotka Peninsula and Mechigmen Bay on the south side of the Chukotka Peninsula. The genus has a non-continuous distribution divided by the Bering Strait and the Brooks Range. We examined the distribution of genetic variation across this range to determine the number and location of potential sub-refugia within the greater Beringian refugium as well as the roles of the Bering land bridge, Brooks Range, and large rivers within Beringia in shaping the current distribution of populations of Dallia. Our analyses were based on DNA sequence data from two nuclear gene introns (S7 and RAG1) and two mitochondrial genome fragments from nineteen sampling locations. These data were examined under genetic clustering and coalescent frameworks to identify sub-refugia within the greater Beringia refugium and to infer the demographic history of different populations of Dallia.

Results

We identified up to five distinct genetic clusters of Dallia. Four of these genetic clusters are present in Alaska: (1) Arctic Coastal Plain genetic cluster found north of the Brooks Range, (2) interior Alaska genetic cluster placed in upstream locations in the Kuskokwim and Yukon river basins, (3) a genetic cluster found on the Seward Peninsula, and (4) a coastal Alaska genetic cluster encompassing downstream Kuskokwim River and Yukon River basin sample locations and samples from Southwest Alaska not in either of these drainages. The Chukotka samples are assigned to their own genetic cluster (5) similar to the coastal Alaska genetic cluster. The clustering and ordination analyses implemented in Discriminant Analysis of Principal Components (DAPC) and STRUCTURE showed mostly concordant groupings and a high degree of differentiation among groups. The groups of sampling locations identified as genetic clusters correspond to geographic areas divided by likely biogeographic barriers including the Brooks Range and the Bering Strait. Estimates of sequence diversity (θ) are highest in the Yukon River and Kuskokwim River drainages near the Bering Sea. We also infer asymmetric migration rates between genetic clusters. The isolation of Dallia on the Arctic Coastal Plain of Alaska is associated with very low estimated migration rates between the coastal Alaska genetic cluster and the Arctic Coastal Plain genetic cluster.

Conclusions

Our results support a scenario with multiple aquatic sub-refugia in Beringia during the Pleistocene and the preservation of that structure in extant populations of Dallia. An inferred historical presence of Dallia across the Bering land bridge explains the similarities in the genetic composition of Dallia in West Beringia and western coastal Alaska. In contrast, historic and contemporary isolation across the Brooks Range shaped the distinctiveness of present day Arctic Coastal Plain Dallia. Overall this study uncovered a high degree of genetic structuring among populations of Dallia supporting the idea of multiple Beringian sub-refugia during the Pleistocene and which appears to be maintained to the present due to the strictly freshwater nature and low dispersal ability of this genus.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0413-2) contains supplementary material, which is available to authorized users.

Less pollen-mediated gene flow for more signatures of glacial lineages: congruent evidence from balsam fir cpDNA and mtDNA for multiple refugia in eastern and central North America

The phylogeographic structure and postglacial history of balsam fir (Abies balsamea), a transcontinental North American boreal conifer, was inferred using mitochondrial DNA (mtDNA) and chloroplast DNA (cpDNA) markers. Genetic structure among 107 populations (mtDNA data) and 75 populations (cpDNA data) was analyzed using Bayesian and genetic distance approaches. Population differentiation was high for mtDNA (dispersed by seeds only), but also for cpDNA (dispersed by seeds and pollen), indicating that pollen gene flow is more restricted in balsam fir than in other boreal conifers. Low cpDNA gene flow in balsam fir may relate to low pollen production due to the inherent biology of the species and populations being decimated by recurrent spruce budworm epidemics, and/or to low dispersal of pollen grains due to their peculiar structural properties. Accordingly, a phylogeographic structure was detected using both mtDNA and cpDNA markers and population structure analyses supported the existence of at least five genetically distinct glacial lineages in central and eastern North America. Four of these would originate from glacial refugia located south of the Laurentide ice sheet, while the last one would have persisted in the northern Labrador region. As expected due to reduced pollen-mediated gene flow, congruence between the geographic distribution of mtDNA and cpDNA lineages was higher than in other North American conifers. However, concordance was not complete, reflecting that restricted but nonetheless detectable cpDNA gene flow among glacial lineages occurred during the Holocene. As a result, new cpDNA and mtDNA genome combinations indicative of cytoplasmic genome capture were observed.

Intervarietal and intravarietal genetic structure in Douglas-fir: nuclear SSRs bring novel insights into past population demographic processes, phylogeography, and intervarietal hybridization

Douglas-fir (Pseudotsuga menziesii) is one of numerous wide-range forest tree species represented by subspecies/varieties, which hybridize in contact zones. This study examined the genetic structure of this North American conifer and its two hybridizing varieties, coastal and Rocky Mountain, at intervarietal and intravarietal level. The genetic structure was subsequently associated with the Pleistocene refugial history, postglacial migration and intervarietal hybridization/introgression. Thirty-eight populations from the USA and Canada were genotyped for 13 nuclear SSRs and analyzed with simulations and traditional population genetic structuring methods. Eight genetic clusters were identified. The coastal clusters embodied five refugial populations originating from five distinct refugia. Four coastal refugial populations, three from California and one from western Canada, diverged during the Pleistocene (56.9–40.1 ka). The three Rocky Mountain clusters reflected distinct refugial populations of three glacial refugia. For Canada, ice covered during the Last Glacial Maximum, we present the following three findings. (1) One refugial population of each variety was revealed in the north of the distribution range. Additional research including paleodata is required to support and determine whether both northern populations originated from cryptic refugia situated south or north of the ice-covered area. (2) An interplay between intravarietal gene flow of different refugial populations and intervarietal gene flow by hybridization and introgression was identified. (3) The Canadian hybrid zone displayed predominantly introgressants of the Rocky Mountain into the coastal variety. This study provides new insights into the complex Quaternary dynamics of this conifer essential for understanding its evolution (outside and inside the native range), adaptation to future climates and for forest management.

Phylogeography of Ramalina menziesii, a widely distributed lichen-forming fungus in western North America

The complex topography and climate history of western North America offer a setting where lineage formation, accumulation and migration have led to elevated inter- and intraspecific biodiversity in many taxa. Here, we study Ramalina menziesii, an epiphytic lichenized fungus with a range encompassing major ecosystems from Baja California to Alaska to explore the predictions of two hypotheses: (i) that the widespread distribution of R. menziesii is due to a single migration episode from a single lineage and (ii) that the widespread distribution is due to the formation and persistence of multiple lineages structured throughout the species' range. To obtain evidence for these predictions, we first construct a phylogenetic tree and identify multiple lineages structured throughout the species' range--some ancient ones that are localized and other more recent lineages that are widely distributed. Second, we use an isolation with migration model to show that sets of ecoregion populations diverged from each other at different times, demonstrating the importance of historical and current barriers to gene flow. Third, we estimated migration rates among ecoregions and find that Baja California populations are relatively isolated, that inland California ecoregion populations do not send out emigrants and that migration out of California coastal and Pacific Northwest populations into inland California ecoregions is high. Such intraspecific geographical patterns of population persistence and dispersal both contribute to the wide range of this genetically diverse lichen fungus and provide insight into the evolutionary processes that enhance species diversity of the California Floristic Province.

Effects of introgression on the genetic population structure of two ecologically and economically important conifer species: lodgepole pine (Pinus contorta var. latifolia) and jack pine (Pinus banksiana)

Forest trees exhibit a remarkable range of adaptations to their environment, but as a result of frequent and long-distance gene flow, populations are often only weakly differentiated. Lodgepole and jack pine hybridize in western Canada, which adds the opportunity for introgression through hybridization to contribute to population structure and (or) adaptive variation. Access to large sample size, high density SNP datasets for these species would improve our ability to resolve population structure, parameterize introgression, and separate the influence of demography from adaptation. To accomplish this, 454 transcriptome reads for lodgepole and jack pine were assembled using Newbler and MIRA, the assemblies mined for SNPs, and 1536 SNPs were selected for typing on lodgepole pine, jack pine, and their hybrids (N = 536). We identified population structure using both Bayesian clustering and discriminate analysis of principle components. Introgressed SNP loci were identified and their influence on observed population structure was assessed. We found that introgressed loci resulted in increased differentiation both within lodgepole and jack pine populations. These findings are timely given the recent mountain pine beetle population expansion in the hybrid zone, and will facilitate future studies of adaptive traits in these ecologically important species.

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.

Large-scale asymmetric introgression of cytoplasmic DNA reveals Holocene range displacement in a North American boreal pine complex

Jack pine (Pinus banksiana) and lodgepole pine (Pinus contorta var. latifolia) are two North American boreal hard pines that hybridize in their zone of contact in western Canada. The main objective of this study was to characterize their patterns of introgression resulting from past and recent gene flow, using cytoplasmic markers having maternal or paternal inheritance. Mitochondrial DNA (mtDNA) and chloroplast DNA (cpDNA) diversity was assessed in allopatric populations of each species and in stands from the current zone of contact containing morphological hybrids. Cluster analyses were used to identify genetic discontinuities among groups of populations. A canonical analysis was also conducted to detect putative associations among cytoplasmic DNA variation, tree morphology, and site ecological features. MtDNA introgression was extensive and asymmetric: it was detected in P. banksiana populations from the hybrid zone and from allopatric areas, but not in P. contorta populations. Very weak cpDNA introgression was observed, and only in P. banksiana populations. The mtDNA introgression pattern indicated that central Canada was first colonized by migrants from a P. contorta glacial population located west of the Rocky Mountains, before being replaced by P. banksiana migrating westward during the Holocene. In contrast, extensive pollen gene flow would have erased the cpDNA traces of this ancient presence of P. contorta. Additional evidence for this process was provided by the results of canonical analysis, which indicated that the current cpDNA background of trees reflected recent pollen gene flow from the surrounding dominant species rather than historical events that took place during the postglacial colonization.

Characterizing the physical and genetic structure of the lodgepole pine × jack pine hybrid zone: mosaic structure and differential introgression

Understanding the physical and genetic structure of hybrid zones can illuminate factors affecting their formation and stability. In north-central Alberta, lodgepole pine (Pinus contorta Dougl. ex Loud. var. latifolia) and jack pine (Pinus banksiana Lamb) form a complex and poorly defined hybrid zone. Better knowledge of this zone is relevant, given the recent host expansion of mountain pine beetle into jack pine. We characterized the zone by genotyping 1998 lodgepole, jack pine, and hybrids from British Columbia, Alberta, Saskatchewan, Ontario, and Minnesota at 11 microsatellites. Using Bayesian algorithms, we calculated genetic ancestry and used this to model the relationship between species occurrence and environment. In addition, we analyzed the ancestry of hybrids to calculate the genetic contribution of lodgepole and jack pine. Finally, we measured the amount of gene flow between the pure species. We found the distribution of the pine classes is explained by environmental variables, and these distributions differ from classic distribution maps. Hybrid ancestry was biased toward lodgepole pine; however, gene flow between the two species was equal. The results of this study suggest that the hybrid zone is complex and influenced by environmental constraints. As a result of this analysis, range limits should be redefined.

Spatial genetic structure of the mountain pine beetle (Dendroctonus ponderosae) outbreak in western Canada: historical patterns and contemporary dispersal

Environmental change has a wide range of ecological consequences, including species extinction and range expansion. Many studies have shown that insect species respond rapidly to climatic change. A mountain pine beetle epidemic of record size in North America has led to unprecedented mortality of lodgepole pine, and a significant range expansion to the northeast of its historic range. Our goal was to determine the spatial genetic variation found among outbreak population from which genetic structure, and dispersal patterns may be inferred. Beetles from 49 sampling locations throughout the outbreak area in western Canada were analysed at 13 microsatellite loci. We found significant north-south population structure as evidenced by: (i) Bayesian-based analyses, (ii) north-south genetic relationships and diversity gradients; and (iii) a lack of isolation-by-distance in the northernmost cluster. The north-south structure is proposed to have arisen from the processes of postglacial colonization as well as recent climate-driven changes in population dynamics. Our data support the hypothesis of multiple sources of origin for the outbreak and point to the need for population specific information to improve our understanding and management of outbreaks. The recent range expansion across the Rocky Mountains into the jack/lodgepole hybrid and pure jack pine zones of northern Alberta is consistent with a northern British Columbia origin. We detected no loss of genetic variability in these populations, indicating that the evolutionary potential of mountain pine beetle to adapt has not been reduced by founder events. This study illustrates a rapid range-wide response to the removal of climatic constraints, and the potential for range expansion of a regional population.

Population structure and migration pattern of a conifer pathogen, Grosmannia clavigera, as influenced by its symbiont, the mountain pine beetle

We investigated the population structure of Grosmannia clavigera (Gc), a fungal symbiont of the mountain pine beetle (MPB) that plays a crucial role in the establishment and reproductive success of this pathogen. This insect-fungal complex has destroyed over 16 million ha of lodgepole pine forests in Canada, the largest MPB epidemic in recorded history. During this current epidemic, MPB has expanded its range beyond historically recorded boundaries, both northward and eastward, and has now reached the jack pine of Alberta, potentially threatening the Canadian boreal forest. To better understand the dynamics between the beetle and its fungal symbiont, we sampled 19 populations in western North America and genotyped individuals from these populations with eight microsatellite markers. The fungus displayed high haplotype diversity, with over 250 unique haplotypes observed in 335 single spore isolates. Linkage equilibria in 13 of the 19 populations suggested that the fungus reproduces sexually. Bayesian clustering and distance analyses identified four genetic clusters that corresponded to four major geographical regions, which suggested that the epidemic arose from multiple geographical sources. A genetic cluster north of the Rocky Mountains, where the MPB has recently become established, experienced a population bottleneck, probably as a result of the recent range expansion. The two genetic clusters located north and west of the Rocky Mountains contained many fungal isolates admixed from all populations, possibly due to the massive movement of MPB during the epidemic. The general agreement in north-south differentiation of MPB and G. clavigera populations points to the fungal pathogen's dependence on the movement of its insect vector. In addition, the patterns of diversity and the individual assignment tests of the fungal associate suggest that migration across the Rocky Mountains occurred via a northeastern corridor, in accordance with meteorological patterns and observation of MPB movement data. Our results highlight the potential of this pathogen for both expansion and sexual reproduction, and also identify some possible barriers to gene flow. Understanding the ecological and evolutionary dynamics of this fungus-beetle association is important for the modelling and prediction of MPB epidemics.

Low levels of population genetic structure in Pinus contorta (Pinaceae) across a geographic mosaic of co-evolution

PREMISE OF THE STUDY

Population genetic analyses provide information on the population context in which evolutionary processes operate and are important for understanding the evolution of geographically variable traits. Earlier studies showed that cone structure of lodgepole pine in the Rocky Mountains diverged among populations because of geographic variation in coevolutionary interactions involving mammalian and avian seed predators. Analyses of population genetic variation are needed to determine whether this divergence has arisen despite extensive gene flow and whether populations to the east and west of the Rocky Mountains have evolved convergent phenotypes independently.

METHODS

We investigated genetic structuring across 22 stands of lodgepole pine in the central Rocky Mountains and in isolated peripheral populations that experience different seed predators and exhibit parallel divergence in cone traits using a set of nine simple sequence repeats and 235 AFLP loci.

KEY RESULTS

Our analyses reveal high levels of genetic diversity within and low genetic differentiation among populations. Nonetheless, geographic and genetic distances were correlated, and isolated populations to the east and west of the Rocky Mountains had higher levels of differentiation than did populations in the central part of the range.

CONCLUSIONS

These data indicate not only that adaptive divergence of cone traits across a geographic mosaic of coevolution has occurred despite minimal genetic differentiation, but also that isolated populations to the east and west of the Rocky Mountains have evolved distinctive cones independently and in parallel. The population structure quantified here will inform future research aimed at detecting genetic variants associated with divergent adaptive traits.

Species delimitation and biogeography of two fir species (Abies) in central China: cytoplasmic DNA variation

It remains unclear how speciation history might contribute to species-specific variation and affect species delimitation. We examined concordance between cytoplasmic genetic variation and morphological taxonomy in two fir species, Abies chensiensis and A. fargesii, with overlapping distributions in central China. Range-wide genetic variation was investigated using mitochondrial (mt) and plastid (pt) DNA sequences, which contrast in their rates of gene flow. Four mtDNA haplotypes were recovered and showed no obvious species' bias in terms of relative frequency. In contrast, a high level of ptDNA variation was recorded in both species with 3 common ptDNA haplotypes shared between them and 21 rare ptDNA haplotypes specific to one or other species. We argue that the lack of concordance between morphological and molecular variation between the two fir species most likely reflects extensive ancestral polymorphism sharing for both forms of cytoplasmic DNA variation. It is feasible that a relatively fast mutation rate for ptDNA contributed to the production of many species-specific ptDNA haplotypes, which remained rare due to insufficient time passing for their spread and fixation in either species, despite high levels of intraspecific ptDNA gene flow. Our phylogeographic analyses further suggest that polymorphisms in both organelle genomes most likely originated during and following glacial intervals preceding the last glacial maximum, when species distributions became fragmented into several refugia and then expanded in range across central China.

Phylogeography of Pinus subsection Australes in the Caribbean Basin

BACKGROUND AND AIMS

Four species of Pinus subsection Australes occur in the Caribbean Basin: P. caribaea, P. cubensis, P. maestrensis and P. occidentalis. This study analyses the phylogeography of these species to assess possible colonization events from Central America to the islands and subsequent population expansions during glacial periods driven by both drier climate and larger emerged land areas.

METHODS

Allele size data were obtained for plastid microsatellites for 314 individuals from 24 populations, covering the distribution range of subsection Australes in the Caribbean Basin.

KEY RESULTS

In total, 113 plastid haplotypes were identified. The highest genetic diversity was found in populations of P. caribaea. Overall, Caribbean Basin populations fit the isolation by distance model. Significant phylogeographical structure was found (R(ST) = 0·671 > permuted R(ST) = 0·101; P < 0·0001). The haplotype network and a Bayesian analysis of population structure (BAPS) indicated different Central American origins for P. caribaea var. bahamensis and P. caribaea var. caribaea plastids, with Central America populations in northern and south-eastern groups. Sudden expansion times for BAPS clusters were close to three glacial maxima.

CONCLUSIONS

Central America contains ancestral plastid haplotypes. Population expansion has played a major role in the distribution of genetic diversity in P. caribaea var. hondurensis. Two colonization events gave rise to the P. caribaea var. bahamensis and P. caribaea var. caribaea lineages. Plastid variation in the eastern species (P. cubensis, P. maestrensis and P. occidentalis) evolved independently from that in P. caribaea var. caribaea. Incomplete lineage sorting between P. cubensis and P. maestrensis is apparent. Inferred expansion times for P. caribaea var. bahamensis and for the eastern lineages correspond to glacial maxima, whereas those for P. caribaea var. hondurensis correspond to the beginning of the temperature decrease that led to Marine Isotope Stage 8.

Cryptic postzygotic isolation in an eruptive species of bark beetle (Dendroctonus ponderosae)

Studies of postzygotic isolation often involve well-differentiated taxa that show a consistent level of incompatibility, thereby limiting our understanding of the initial stages and development of reproductive barriers. Dendroctonus ponderosae provides an informative system because recent evidence suggests that distant populations produce hybrids with reproductive incompatibilities. Dendroctonus ponderosae shows an isolation-by-distance gene flow pattern allowing us to characterize the evolution of postzygotic isolation (e.g., hybrid inviability, hybrid sterility) by crossing populations along a continuum of geographic/genetic divergence. We found little evidence of hybrid inviability among these crosses. However, crosses between geographically distant populations produced sterile males (consistent with Haldane's rule). This effect was not consistent with the fixation of mutations in an isolation-by-distance pattern, but instead is spatially localized. These reproductive barriers are uncorrelated with a reduction in gene flow suggesting their recent development. Crosses between geographically proximal populations bounding the transition from compatibility to hybrid male sterility showed evidence of unidirectional reduction in hybrid male fecundity. Our study describes significant postzygotic isolation occurring across a narrow and molecularly cryptic geographic zone between the states of Oregon and Idaho. This study provides a view of the early stages of postzygotic isolation in a geographically widespread species.

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.

Gene flow and species delimitation: a case study of two pine species with overlapping distributions in southeast China

Species delimitation detected by molecular markers is complicated by introgression and incomplete lineage sorting between species. Recent modeling suggests that fixed genetic differences between species are highly related to rates of intraspecific gene flow. However, it remains unclear whether such differences are due to high levels of intraspecific gene flow overriding the spread of introgressed alleles or favoring rapid lineage sorting between species. In pines, chloroplast (cp) and mitochondrial (mt) DNAs are normally paternally and maternally inherited, respectively, and thus their relative rates of intraspecific gene flow are expected to be high and low, respectively. In this study, we used two pine species with overlapping geographical distributions in southeast China, P. massoniana and P. hwangshanensis, as a model system to examine the association between organelle gene flow and variation within and between species. We found that cpDNA variation across these two pine species is more species specific than mtDNA variation and almost delimits taxonomic boundaries. The shared mt/cp DNA genetic variation between species shows no bias in regard to parapatric versus allopatric species' distributions. Our results therefore support the hypothesis that high intraspecific gene flow has accelerated cpDNA lineage sorting between these two pine species.