Phenotypic selection on growth rhythm in whitebark pine under climatic conditions warmer than seed origins

Growth rhythm that is well synchronized with seasonal changes in local climatic conditions is understood to enhance fitness; however, rapid ongoing climate change threatens to disrupt this synchrony. To evaluate phenotypic selection on growth rhythm under expected warmer and drier future climate, seedlings from 49 populations of whitebark pine (Pinus albicaulis Engelm.) were grown and measured over more than 10 years in two common garden field experiments on sites that approximate the projected future climate of the seed origins. Selection on growth rhythm was assessed by relating individual plant fitness to timing and rate of shoot elongation. Differential survival clearly evidenced selection on growth rhythm. We detected directional and stabilizing selection that varied in magnitude between experimental sites and among years. The observed phenotypic selection supports the interpretation of clinal variation among populations within tree species as reflecting adaptive variation in response to past natural selection mediated by climate. To the extent that growth rhythm is heritable, results of the present study suggest evolution of whitebark pine toward a more distinct timing of shoot elongation and generally more rapid elongation in the immediate next generation under ongoing climate change in environments similar to the study sites.     

Genomewide scan for adaptive differentiation along altitudinal gradient in the Andrew's toad Bufo andrewsi

Recent studies of humans, dogs and rodents have started to discover the genetic underpinnings of high altitude adaptations, yet amphibians have received little attention in this respect. To identify possible signatures of adaptation to altitude, we performed a genome scan of 15 557 single nucleotide polymorphisms (SNPs) obtained with restriction site‐associated DNA sequencing of pooled samples from 11 populations of Andrew's toad (Bufo andrewsi) from the edge of the Tibetan Plateau, spanning an altitudinal gradient from 1690 to 2768 m.a.s.l. We discovered significant geographic differentiation among all sites, with an average F_ST = 0.023 across all SNPs. Apart from clear patterns of isolation by distance, we discovered numerous outlier SNPs showing strong associations with variation in altitude (1394 SNPs), average annual temperature (1859 SNPs) or both (1051 SNPs). Levels and patterns of genetic differentiation in these SNPs were consistent with the hypothesis that they have been subject to directional selection and reflect adaptation to altitudinal variation among the study sites. Genes with footprints of selection were significantly enriched in binding and metabolic processes. Several genes potentially related to high altitude adaptation were identified, although the identity and functional significance of most genomic targets of selection remain unknown. In general, the results provide genomic support for results of earlier common garden and low coverage genetic studies that have uncovered substantial adaptive differentiation along altitudinal and latitudinal gradients in amphibians.     

Climate‐associated genetic variation in Fagus sylvatica and potential responses to climate change in the French Alps

Local adaptation patterns have been found in many plants and animals, highlighting the genetic heterogeneity of species along their range of distribution. In the next decades, global warming is predicted to induce a change in the selective pressures that drive this adaptive variation, forcing a reshuffling of the underlying adaptive allele distributions. For species with low dispersion capacity and long generation time such as trees, the rapidity of the change could impede the migration of beneficial alleles and lower their capacity to track the changing environment. Identifying the main selective pressures driving the adaptive genetic variation is thus necessary when investigating species capacity to respond to global warming. In this study, we investigate the adaptive landscape of Fagus sylvatica along a gradient of populations in the French Alps. Using a double‐digest restriction‐site‐associated DNA (ddRAD) sequencing approach, we identified 7,000 SNPs from 570 individuals across 36 different sites. A redundancy analysis (RDA)‐derived method allowed us to identify several SNPs that were strongly associated with climatic gradients; moreover, we defined the primary selective gradients along the natural populations of F. sylvatica in the Alps. Strong effects of elevation and humidity, which contrast north‐western and south‐eastern site, were found and were believed to be important drivers of genetic adaptation. Finally, simulations of future genetic landscapes that used these findings allowed identifying populations at risk for F. sylvatica in the Alps, which could be helpful for future management plans.     

The population genomic signature of environmental association and gene flow in an ecologically divergent tree species Metrosideros polymorpha (Myrtaceae)

Genomewide markers enable us to study genetic differentiation within a species and the factors underlying it at a much higher resolution than before, which advances our understanding of adaptation in organisms. We investigated genomic divergence in Metrosideros polymorpha, a woody species that occupies a wide range of ecological habitats across the Hawaiian Islands and shows remarkable phenotypic variation. Using 1659 single nucleotide polymorphism (SNP) markers annotated with the genome assembly, we examined the population genetic structure and demographic history of nine populations across five elevations and two ages of substrates on Mauna Loa, the island of Hawaii. The nine populations were differentiated into two genetic clusters distributed on the lower and higher elevations and were largely admixed on the middle elevation. Demographic modelling revealed that the two genetic clusters have been maintained in the face of gene flow, and the effective population size of the high‐altitude cluster was much smaller. A F_ST‐based outlier search among the 1659 SNPs revealed that 34 SNPs (2.05%) were likely to be under divergent selection and the allele frequencies of 21 of them were associated with environmental changes along elevations, such as temperature and precipitation. This study shows a genomic mosaic of M. polymorpha, in which contrasting divergence patterns were found. While most genomic polymorphisms were shared among populations, a small fraction of the genome was significantly differentiated between populations in diverse environments and could be responsible for the dramatic adaptation to a wide range of environments.     

netview p: a network visualization tool to unravel complex population structure using genome‐wide SNPs

Network‐based approaches are emerging as valuable tools for the analysis of complex genetic structure in wild and captive populations. netview p combines data quality control with the construction of population networks through mutual k‐nearest neighbours thresholds applied to genome‐wide SNPs. The program is cross‐platform compatible, open‐source and efficiently operates on data ranging from hundreds to hundreds of thousands of SNPs. The pipeline was used for the analysis of pedigree data from simulated (n = 750, SNPs = 1279) and captive silver‐lipped pearl oysters (n = 415, SNPs = 1107), wild populations of the European hake from the Atlantic and Mediterranean (n = 834, SNPs = 380) and grey wolves from North America (n = 239, SNPs = 78 255). The population networks effectively visualize large‐ and fine‐scale genetic structure within and between populations, including family‐level structure and relationships. netview p comprises a network‐based addition to other population analysis tools and provides user‐friendly access to a complex network analysis pipeline through implementation in python .     

Genetic signatures of natural selection in response to air pollution in red spruce (Picea rubens,Pinaceae)

One of the most important drivers of local adaptation for forest trees is climate. Coupled to these patterns, however, are human‐induced disturbances through habitat modification and pollution. The confounded effects of climate and disturbance have rarely been investigated with regard to selective pressure on forest trees. Here, we have developed and used a population genetic approach to search for signals of selection within a set of 36 candidate genes chosen for their putative effects on adaptation to climate and human‐induced air pollution within five populations of red spruce (Picea rubens Sarg.), distributed across its natural range and air pollution gradient in eastern North America. Specifically, we used F_ST outlier and environmental correlation analyses to highlight a set of seven single nucleotide polymorphisms (SNPs) that were overly correlated with climate and levels of sulphate pollution after correcting for the confounding effects of population history. Use of three age cohorts within each population allowed the effects of climate and pollution to be separated temporally, as climate‐related SNPs (n = 7) showed the strongest signals in the oldest cohort, while pollution‐related SNPs (n = 3) showed the strongest signals in the youngest cohorts. These results highlight the usefulness of population genetic scans for the identification of putatively nonneutral evolution within genomes of nonmodel forest tree species, but also highlight the need for the development and application of robust methodologies to deal with the inherent multivariate nature of the genetic and ecological data used in these types of analyses.     

EFFECTS OF AVIAN SEED DISPERSAL ON THE GENETIC STRUCTURE OF WHITEBARK PINE POPULATIONS

We used allozyme analysis to examine family structure, the spatial patterning of related individuals, in two populations of whitebark pine (Pinus albicaulis), a subalpine conifer that commonly displays a multistem form. The individual stems within clumps are genetically distinct individuals, having arisen from separate seeds. Individuals within a clump are genetically more similar than individuals in different clumps, but individuals in neighboring clumps do not appear to be more similar than individuals in distant clumps. This family structure appears to be a direct result of the seed-caching behavior of Clark's nutcrackers (Nucifraga columbiana), the primary dispersal agent for whitebark pine seeds.     

Phylogeography and population genetics of pine butterflies: Sky islands increase genetic divergence

The sky islands of southeastern Arizona (AZ) mark a major transition zone between tropical and temperate biota and are considered a neglected biodiversity hotspot. Dispersal ability and host plant specificity are thought to impact the history and diversity of insect populations across the sky islands. We aimed to investigate the population structure and phylogeography of two pine‐feeding pierid butterflies, the pine white (Neophasia menapia) and the Mexican pine white (Neophasia terlooii), restricted to these “islands” at this transition zone. Given their dependence on pines as the larval hosts, we hypothesized that habitat connectivity affects population structure and is at least in part responsible for their allopatry. We sampled DNA from freshly collected butterflies from 17 sites in the sky islands and adjacent high‐elevation habitats and sequenced these samples using ddRADSeq. Up to 15,399 SNPs were discovered and analyzed in population genetic and phylogenetic contexts with Stacks and pyRAD pipelines. Low genetic differentiation in N. menapia suggests that it is panmictic. Conversely, there is strong evidence for population structure within N. terlooii. Each sky island likely contains a population of N. terlooii, and clustering is hierarchical, with populations on proximal mountains being more related to each other. The N. menapia habitat, which is largely contiguous, facilitates panmixia, while the N. terlooii habitat, restricted to the higher elevations on each sky island, creates distinct population structure. Phylogenetic results corroborate those from population genetic analyses. The historical climate‐driven fluxes in forest habitat connectivity have implications for understanding the biodiversity of fragmented habitats.     

Whitebark pine facilitation at treeline: potential interactions for disruption by an invasive pathogen

In stressful environments, facilitation often aids plant establishment, but invasive plant pathogens may potentially disrupt these interactions. In many treeline communities in the northern Rocky Mountains of the U.S. and Canada, Pinus albicaulis, a stress‐tolerant pine, initiates tree islands at higher frequencies than other conifers – that is, leads to leeward tree establishment more frequently. The facilitation provided by a solitary (isolated) P. albicaulis leading to tree island initiation may be important for different life‐history stages for leeward conifers, but it is not known which life‐history stages are influenced and protection provided. However, P. albicaulis mortality from the non‐native pathogen Cronartium ribicola potentially disrupts these facilitative interactions, reducing tree island initiation. In two Rocky Mountain eastern slope study areas, we experimentally examined fundamental plant–plant interactions which might facilitate tree island formation: the protection offered by P. albicaulis to leeward seed and seedling life‐history stages, and to leeward krummholz conifers. In the latter case, we simulated mortality from C. ribicola for windward P. albicaulis to determine whether loss of P. albicaulis from C. ribicola impacts leeward conifers. Relative to other common solitary conifers at treeline, solitary P. albicaulis had higher abundance. More seeds germinated in leeward rock microsites than in conifer or exposed microsites, but the odds of cotyledon seedling survival during the growing season were highest in P. albicaulis microsites. Planted seedling survival was low among all microsites examined. Simulating death of windward P. albicaulis by C. ribicola reduced shoot growth of leeward trees. Loss of P. albicaulis to exotic disease may limit facilitation interactions and conifer community development at treeline and potentially impede upward movement as climate warms.     

The influence of genetics,defensive chemistry and the fungal microbiome on disease outcome in whitebark pine trees

The invasive fungal pathogen Cronartium ribicola infects and kills whitebark pine (Pinus albicaulis) throughout western North America. Whitebark pine has been proposed for listing under the Endangered Species Act in the USA, and the loss of this species is predicted to have severe impacts on ecosystem composition and function in high‐elevation forests. Numerous fungal endophytes live inside whitebark pine tissues and may influence the severity of C. ribicola infection, either directly by inhibition of pathogen growth or indirectly by the induction of chemical defensive pathways in the tree. Terpenes, a form of chemical defence in pine trees, can also influence disease. In this study, we characterized fungal endophyte communities in whitebark pine seedlings before and after experimental inoculation with C. ribicola, monitored disease progression and compared fungal community composition in susceptible vs. resistant seedlings in a common garden. We analysed the terpene composition of these same seedlings. Seed family identity or maternal genetics influenced both terpenes and endophyte communities. Terpene and endophyte composition correlated with disease severity, and terpene concentrations differed in resistant vs. susceptible seedlings. These results suggest that the resistance to C. ribicola observed in natural whitebark pine populations is caused by the combined effects of genetics, endophytes and terpenes within needle tissue, in which initial interactions between microbes and hosts take place. Tree genotype, terpene and microbiome combinations associated with healthy trees could help to predict or reduce disease severity and improve outcomes of future tree breeding programmes.     

Heterogeneous genetic invasions of three insecticide resistance mutations in Indo‐Pacific populations of Aedes aegypti (L.)

Nations throughout the Indo‐Pacific region use pyrethroid insecticides to control Aedes aegypti, the mosquito vector of dengue, often without knowledge of pyrethroid resistance status of the pest or origin of resistance. Two mutations (V1016G + F1534C) in the sodium channel gene (Vssc) of Ae. aegypti modify ion channel function and cause target‐site resistance to pyrethroid insecticides, with a third mutation (S989P) having a potential additive effect. Of 27 possible genotypes involving these mutations, some allelic combinations are never seen whereas others predominate. Here, five allelic combinations common in Ae. aegypti from the Indo‐Pacific region are described and their geographical distributions investigated using genome‐wide SNP markers. We tested the hypothesis that resistance allele combinations evolved de novo in populations versus the alternative that dispersal of Ae. aegypti between populations facilitated genetic invasions of allele combinations. We used latent factor mixed‐models to detect SNPs throughout the genome that showed structuring in line with resistance allele combinations and compared variation at SNPs within the Vssc gene with genome‐wide variation. Mixed‐models detected an array of SNPs linked to resistance allele combinations, all located within or in close proximity to the Vssc gene. Variation at SNPs within the Vssc gene was structured by resistance profile, whereas genome‐wide SNPs were structured by population. These results demonstrate that alleles near to resistance mutations have been transferred between populations via linked selection. This indicates that genetic invasions have contributed to the widespread occurrence of Vssc allele combinations in Ae. aegypti in the Indo‐Pacific region, pointing to undocumented mosquito invasions between countries.     

Evidence of genomic adaptation to climate in Eucalyptus microcarpa: Implications for adaptive potential to projected climate change

Understanding whether populations can adapt in situ or whether interventions are required is of key importance for biodiversity management under climate change. Landscape genomics is becoming an increasingly important and powerful tool for rapid assessments of climate adaptation, especially in long‐lived species such as trees. We investigated climate adaptation in Eucalyptus microcarpa using the DArTseq genomic approach. A combination of F_ST outlier and environmental association analyses were performed using >4200 genomewide single nucleotide polymorphisms (SNPs) from 26 populations spanning climate gradients in southeastern Australia. Eighty‐one SNPs were identified as putatively adaptive, based on significance in F_ST outlier tests and significant associations with one or more climate variables related to temperature (70/81), aridity (37/81) or precipitation (35/81). Adaptive SNPs were located on all 11 chromosomes, with no particular region associated with individual climate variables. Climate adaptation appeared to be characterized by subtle shifts in allele frequencies, with no consistent fixed differences identified. Based on these associations, we predict adaptation under projected changes in climate will include a suite of shifts in allele frequencies. Whether this can occur sufficiently rapidly through natural selection within populations, or would benefit from assisted gene migration, requires further evaluation. In some populations, the absence or predicted increases to near fixation of particular adaptive alleles hint at potential limits to adaptive capacity. Together, these results reinforce the importance of standing genetic variation at the geographic level for maintaining species’ evolutionary potential.     

Cunningham's skinks show low genetic connectivity and signatures of divergent selection across its distribution

Establishing corridors of connecting habitat has become a mainstay conservation strategy to maintain gene flow and facilitate climate‐driven range shifts. Yet, little attention has been given to ascertaining the extent to which corridors will benefit philopatric species, which might exhibit localized adaptation. Measures of genetic connectivity and adaptive genetic variation across species’ ranges can help fill this knowledge gap. Here, we characterized the spatial genetic structure of Cunningham's skink (Egernia cunninghami), a philopatric species distributed along Australia's Great Dividing Range, and assessed evidence of localized adaptation. Analysis of 4,274 SNPs from 94 individuals sampled at four localities spanning 500 km and 4° of latitude revealed strong genetic structuring at neutral loci (mean F_ST ± SD = 0.603 ± 0.237) among the localities. Putatively neutral SNPs and those under divergent selection yielded contrasting spatial patterns, with the latter identifying two genetically distinct clusters. Given low genetic connectivity of the four localities, we suggest that the natural movement rate of this species is insufficient to keep pace with spatial shifts to its climate envelope, irrespective of habitat availability. In addition, our finding of localized adaptation highlights the risk of outbreeding depression should the translocation of individuals be adopted as a conservation management strategy.     

The challenge of separating signatures of local adaptation from those of isolation by distance and colonization history: The case of two white pines

Accurately detecting signatures of local adaptation using genetic‐environment associations (GEAs) requires controlling for neutral patterns of population structure to reduce the risk of false positives. However, a high degree of collinearity between climatic gradients and neutral population structure can greatly reduce power, and the performance of GEA methods in such case is rarely evaluated in empirical studies. In this study, we attempted to disentangle the effects of local adaptation and isolation by environment (IBE) from those of isolation by distance (IBD) and isolation by colonization from glacial refugia (IBC) using range‐wide samples in two white pine species. For this, SNPs from 168 genes, including 52 candidate genes for growth and phenology, were genotyped in 133 and 61 populations of Pinus strobus and P. monticola, respectively. For P. strobus and using all 153 SNPs, climate (IBE) did not significantly explained among‐population variation when controlling for IBD and IBC in redundancy analyses (RDAs). However, 26 SNPs were significantly associated with climate in single‐locus GEA analyses (Bayenv2 and LFMM), suggesting that local adaptation took place in the presence of high gene flow. For P. monticola, we found no evidence of IBE using RDAs and weaker signatures of local adaptation using GEA and F_ST outlier tests, consistent with adaptation via phenotypic plasticity. In both species, the majority of the explained among‐population variation (69 to 96%) could not be partitioned between the effects of IBE, IBD, and IBC. GEA methods can account differently for this confounded variation, and this could explain the small overlap of SNPs detected between Bayenv2 and LFMM. Our study illustrates the inherent difficulty of taking into account neutral structure in natural populations and the importance of sampling designs that maximize climatic variation, while minimizing collinearity between climatic gradients and neutral structure.     

Population structure and local selection yield high genomic variation in Mimulus guttatus

Across western North America, Mimulus guttatus exists as many local populations adapted to site‐specific environmental challenges. Gene flow between locally adapted populations will affect genetic diversity both within demes and across the larger metapopulation. Here, we analyse 34 whole‐genome sequences from the intensively studied Iron Mountain population (IM) in conjunction with sequences from 22 Mimulus individuals sampled from across western North America. Three striking features of these data address hypotheses about migration and selection in a locally adapted population. First, we find very high levels of intrapopulation polymorphism (synonymous π = 0.033). Variation outside of genes is likely even higher but difficult to estimate because excessive divergence reduces the efficiency of read mapping. Second, IM exhibits a significantly positive genomewide average for Tajima's D. This indicates allele frequencies are typically more intermediate than expected from neutrality, opposite the pattern observed in many other species. Third, IM exhibits a distinctive haplotype structure with a genomewide excess of positive associations between rarer alleles at linked loci. This suggests an important effect of gene flow from other Mimulus populations, although a residual effect of population founding might also contribute. The combination of multiple analyses, including a novel tree‐based analytic method, illustrates how the balance of local selection, limited dispersal and metapopulation dynamics manifests across the genome. The overall genomic pattern of sequence diversity suggests successful gene flow of divergent immigrant genotypes into IM. However, many loci show patterns indicative of local adaptation, particularly at SNPs associated with chromosomal inversions.     

Development and preliminary evaluation of a genomewide single nucleotide polymorphisms resource generated by RAD‐seq for the small yellow croaker (Larimichthys polyactis)

Recent advances in high‐throughput sequencing technologies have offered the possibility to generate genomewide sequence data to delineate previously unidentified genetic structure, obtain more accurate estimates of demographic parameters and to evaluate potential adaptive divergence. Here, we identified 27 556 single nucleotide polymorphisms for the small yellow croaker (Larimichthys polyactis) using restriction‐site‐associated DNA (RAD) sequencing of 24 individuals from two populations. Significant sources of genetic variation were identified, with an average nucleotide diversity (π) of 0.00105 ± 0.000425 across individuals, and long‐term effective population size was thus estimated to range between 26 172 and 261 716. According to the results, no differentiation between the two populations was detected based on the SNP data set of top quality score per contig or neutral loci. However, the two analysed populations were highly differentiated based on SNP data set of both top F_ST value per contig and the outlier SNPs. Moreover, local adaptation was highlighted by an F_ST‐based outlier tests implemented in LOSITAN and a total of 538 potentially locally selected SNPs were identified. blast2go annotation of contigs containing the outlier SNPs yielded hits for 37 (66%) of 56 significant blastx matches. Candidate genes for local adaptation constituted a wide array of biological functions, including cellular response to oxidative stress, actin filament binding, ion transmembrane transport and synapse assembly. The generated SNP resources in this study provided a valuable tool for future population genetics and genomics studies of L. polyactis.     

Inbreeding and conservation genetics in whitebark pine

Whitebark pine (Pinus albicaulisEngelm.) is threatened across its native rangeby an exotic fungal pathogen introduced withinthe last century. Mortality has beenextensive, and projected potential range shiftsbased on impending climate change have revealedfurther pressures to survival and adaptationfor this long-lived, high-elevation conifer. Quantifying genetic variation and the matingsystem of whitebark pine in its northern rangeprovides a basis for effective conservationmeasures. Isozyme analysis of vegetative budtissue revealed high expected heterozygosity(0.262), moderate population differentiation(F_ST = 0.061) and highly significantcorrelations between observed heterozygosityand geographic variables (R² = 0.36,latitude; R² = 0.30 longitude), supportingthe hypothesis that this species recolonizedits current northern range following glacialretreat from several refugia in the Washingtonand Oregon Cascades and in the northernRockies. Mating system analysis based onsimultaneous isozyme analyses of embryo andhaploid megagametophyte tissues foundrelatively high levels of consanguineous matingand selfing for a conifer (t _m =0.73) within populations. Avian seeddistribution by the Clark's nutcracker (Nucifragia columbiana Wilson) appears to bethe overriding factor influencing geneticpatterns: being a mutualistic seed disperser,caches comprised of related seeds develop intoclumped stands with strong family substructure. While it is a critical wildlife habitatcomponent, lack of commercial utilization hasmade in situ adaptation the primaryconservation focus. Encouraging regenerationsuccess and nutcracker caching by maintainingnatural fire regimes will provide anecosystem-based conservation solution; however,in the Rocky Mountains between 52° N and47° N, disease-resistant individualsshould be located and propagated in order toensure long-term survival of the species inhigh pathogen hazard areas.     

Comparison of survey methods for monitoring Clark's Nutcrackers and predicting dispersal of whitebark pine seeds

ABSTRACT Clark's Nutcrackers (Nucifraga columbiana) disperse seeds of whitebark pines (Pinus albicaulis) in western North America by their scatter‐hoarding behavior. Because of declines in whitebark pine, resource managers are seeking an effective means of monitoring nutcracker population trends and the probability of seed dispersal by nutcrackers. We tested the reliability of four survey techniques (standard point counts, playback point counts, line transects, and Breeding Bird Survey routes) for estimating population size by conducting surveys at sites where a portion of the nutcracker population was marked with radio transmitters. The efficacy of distance sampling, based on detection rates from our unadjusted surveys, was also assessed. We conducted counts of whitebark pine cones within stands and related the probability of seed dispersal within stands to cone production and nutcracker abundance. We conducted 70 h of surveys for Clark's Nutcrackers at eight sites from July through November in 2007 and 2009 and estimated cone densities at six of these sites. Detection rates for all survey techniques were low and variable and we detected an average of 5.6 nutcrackers per 30 min of survey time. We also found no difference in detection rates among survey types, although significantly more nutcrackers were detected during surveys conducted during the peak of whitebark pine cone harvest (P < 0.0001). Nutcracker abundance was not correlated with cone density (P= 0.29) and we observed nutcrackers pouching seeds at all sites. Thus, cone density did not provide reliable information on whether seed dispersal was likely to occur. We suggest that alternate methods be considered for monitoring populations and assessing seed dispersal probability because we did not reliably detect nutcrackers using conventional survey techniques and because nutcracker abundance was not correlated with cone density.     

Influences of the biophysical environment on blister rust and mountain pine beetle,and their interactions,in whitebark pine forests

Aim To understand how the biophysical environment influences patterns of infection by non‐native blister rust (caused by Cronartium ribicola) and mortality caused by native mountain pine beetles (Dendroctonus ponderosae) in whitebark pine (Pinus albicaulis) communities, to determine how these disturbances interact, and to gain insight into how climate change may influence these patterns in the future. Location High‐elevation forests in south‐west Montana, central Idaho, eastern and western Oregon, USA. Methods Stand inventory and dendroecological methods were used to assess stand structure and composition and to reconstruct forest history at sixty 0.1‐ha plots. Patterns of blister rust infection and mountain pine beetle‐caused mortality in whitebark pine trees were examined using nonparametric Kruskal–Wallis ANOVA, Mann–Whitney U‐tests, and Kolmogorov–Smirnov two‐sample tests. Stepwise regression was used to build models of blister rust infection and mountain pine beetle‐related mortality rates based on a suite of biophysical site variables. Results Occurrence of blister rust infections was significantly different among the mountain ranges, with a general gradient of decreasing blister rust occurrence from east to west. Evidence of mountain pine beetle‐caused mortality was identified on 83% of all dead whitebark pine trees and was relatively homogenous across the study area. Blister rust infected trees of all ages and sizes uniformly, while mountain pine beetles infested older, larger trees at all sites. Stepwise regressions explained 64% and 58% of the variance in blister rust infection and beetle‐caused mortality, respectively, indicating that these processes are strongly influenced by the biophysical environment. More open stand structures produced by beetle outbreaks may increase the exposure of surviving whitebark pine trees to blister rust infection. Main conclusions Variability in the patterns of blister rust infection and mountain pine beetle‐caused mortality elucidated the fundamental dynamics of these disturbance agents and suggests that the effects of climate change will be complex in whitebark pine communities and vary across the species’ range. Interactions between blister rust and beetle outbreaks may accelerate declines or facilitate the rise of rust resistance in whitebark pine depending on forest conditions at the time of the outbreak.     

Genetic variation and bidirectional gene flow in the riparian plant Miscanthus lutarioriparius,across its endemic range: implications for adaptive potential

Miscanthus lutarioriparius is an endemic species that grows along the middle and lower reaches of the Yangtze River and is a valuable source of germplasm for the development of second‐generation energy crops. The plant that propagates via seeds, stem nodes, and rhizomes shows high phenotypic variation and strong local adaptation. Here, we examined the magnitude and spatial distribution of genetic variation in M. lutarioriparius across its entire distributional range and tested underlying factors that shaped its genetic variation. Population genetic analyses were conducted on 644 individuals from 25 populations using 16 microsatellite markers. M. lutarioriparius exhibited a high level of genetic variation (H_E = 0.682–0.786; A_r = 4.74–8.06) and a low differentiation (F_ST = 0.063; D_est = 0.153). Of the total genetic variation, 10% was attributed to the differences among populations (df = 24, P < 0.0001), whereas 90% was attributed to the differences among individuals (df = 619, P ≤ 0.0001). Genetic diversity did not differ significantly across longitudes and did not increase in the populations growing downstream of the Yangtze River. However, significant associations were found between genetic differentiation and spatial distance. Six genetic discontinuities were identified, which mostly distributed among downstream populations. We conclude that anthropogenic factors and landscape features both contributed to shaping the pattern of gene flow in M. lutarioriparius, including long‐distance bidirectional dispersal. Our results explain the genetic basis of the high degree of adaptability in M. lutarioriparius and identify potential sources of new germplasm for the domestication of this potential second‐generation energy crop.