Water availability drives signatures of local adaptation in whitebark pine (Pinus albicaulis Engelm.) across fine spatial scales of the Lake Tahoe Basin,USA

Patterns of local adaptation at fine spatial scales are central to understanding how evolution proceeds, and are essential to the effective management of economically and ecologically important forest tree species. Here, we employ single and multilocus analyses of genetic data (n = 116 231 SNPs) to describe signatures of fine‐scale adaptation within eight whitebark pine (Pinus albicaulis Engelm.) populations across the local extent of the environmentally heterogeneous Lake Tahoe Basin, USA. We show that despite highly shared genetic variation (F_ST = 0.0069), there is strong evidence for adaptation to the rain shadow experienced across the eastern Sierra Nevada. Specifically, we build upon evidence from a common garden study and find that allele frequencies of loci associated with four phenotypes (mean = 236 SNPs), 18 environmental variables (mean = 99 SNPs), and those detected through genetic differentiation (n = 110 SNPs) exhibit significantly higher signals of selection (covariance of allele frequencies) than could be expected to arise, given the data. We also provide evidence that this covariance tracks environmental measures related to soil water availability through subtle allele frequency shifts across populations. Our results replicate empirical support for theoretical expectations of local adaptation for populations exhibiting strong gene flow and high selective pressures and suggest that ongoing adaptation of many P. albicaulis populations within the Lake Tahoe Basin will not be constrained by the lack of genetic variation. Even so, some populations exhibit low levels of heritability for the traits presumed to be related to fitness. These instances could be used to prioritize management to maintain adaptive potential. Overall, we suggest that established practices regarding whitebark pine conservation be maintained, with the additional context of fine‐scale adaptation.     

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.     

Time‐specific and pleiotropic quantitative trait loci coordinately modulate stem growth in Populus

In perennial woody plants, the coordinated increase of stem height and diameter during juvenile growth improves competitiveness (i.e. access to light); however, the factors underlying variation in stem growth remain unknown in trees. Here, we used linkage‐linkage disequilibrium (linkage‐LD) mapping to decipher the genetic architecture underlying three growth traits during juvenile stem growth. We used two Populus populations: a linkage mapping population comprising a full‐sib family of 1,200 progeny and an association mapping panel comprising 435 unrelated individuals from nearly the entire natural range of Populus tomentosa. We mapped 311 quantitative trait loci (QTL) for three growth traits at 12 timepoints to 42 regions in 17 linkage groups. Of these, 28 regions encompassing 233 QTL were annotated as 27 segmental homology regions (SHRs). Using SNPs identified by whole‐genome re‐sequencing of the 435‐member association mapping panel, we identified significant SNPs (P ≤ 9.4 × 10^?7) within 27 SHRs that affect stem growth at nine timepoints with diverse additive and dominance patterns, and these SNPs exhibited complex allelic epistasis over the juvenile growth period. Nineteen genes linked to potential causative alleles that have time‐specific or pleiotropic effects, and mostly overlapped with significant signatures of selection within SHRs between climatic regions represented by the association mapping panel. Five genes with potential time‐specific effects showed species‐specific temporal expression profiles during the juvenile stages of stem growth in five representative Populus species. Our observations revealed the importance of considering temporal genetic basis of complex traits, which will facilitate the molecular design of tree ideotypes.     

Climatic adaptation and ecological divergence between two closely related pine species in Southeast China

Climate is one of the most important drivers for adaptive evolution in forest trees. Climatic selection contributes greatly to local adaptation and intraspecific differentiation, but this kind of selection could also have promoted interspecific divergence through ecological speciation. To test this hypothesis, we examined intra‐ and interspecific genetic variation at 25 climate‐related candidate genes and 12 reference loci in two closely related pine species, Pinus massoniana Lamb. and Pinus hwangshanensis Hisa, using population genetic and landscape genetic approaches. These two species occur in Southeast China but have contrasting ecological preferences in terms of several environmental variables, notably altitude, although hybrids form where their distributions overlap. One or more robust tests detected signals of recent and/or ancient selection at two‐thirds (17) of the 25 candidate genes, at varying evolutionary timescales, but only three of the 12 reference loci. The signals of recent selection were species specific, but signals of ancient selection were mostly shared by the two species likely because of the shared evolutionary history. F_ST outlier analysis identified six SNPs in five climate‐related candidate genes under divergent selection between the two species. In addition, a total of 24 candidate SNPs representing nine candidate genes showed significant correlation with altitudinal divergence in the two species based on the covariance matrix of population history derived from reference SNPs. Genetic differentiation between these two species was higher at the candidate genes than at the reference loci. Moreover, analysis using the isolation‐with‐migration model indicated that gene flow between the species has been more restricted for climate‐related candidate genes than the reference loci, in both directions. Taken together, our results suggest that species‐specific and divergent climatic selection at the candidate genes might have counteracted interspecific gene flow and played a key role in the ecological divergence of these two closely related pine species.     

Exploring genomic variation associated with drought stress in Picea mariana populations

Predicted increases in drought and heat stress will likely induce shifts in species bioclimatic envelopes. Genetic variants adapted to water limitation may prove pivotal for species response under scenarios of increasing drought. In this study, we aimed to explore this hypothesis by investigating genetic variation in 16 populations of black spruce (Picea mariana) in relation to climate variables in Alaska. A total of 520 single nucleotide polymorphisms (SNPs) were genotyped for 158 trees sampled from areas of contrasting climate regimes. We used multivariate and univariate genotype‐by‐environment approaches along with available gene annotations to investigate the relationship between climate and genetic variation among sampled populations. Nine SNPs were identified as having a significant association with climate, of which five were related to drought stress response. Outlier SNPs with respect to the overall environment were significantly overrepresented for several biological functions relevant for coping with variable hydric regimes, including osmotic stress response. This genomic imprint is consistent with local adaptation of black spruce to drought stress. These results suggest that natural selection acting on standing variation prompts local adaptation in forest stands facing water limitation. Improved understanding of possible adaptive responses could inform our projections about future forest dynamics and help prioritize populations that harbor valuable genetic diversity for conservation.     

Targeted high‐throughput growth hormone 1 gene sequencing reveals high within‐breed genetic diversity in South African goats

This study assessed the genetic diversity in the growth hormone 1 gene (GH1) within and between South African goat breeds. Polymerase chain reaction‐targeted gene amplification together with Illumina MiSeq next‐generation sequencing (NGS) was used to generate the full length (2.54 kb) of the growth hormone 1 gene and screen for SNPs in the South African Boer (SAB) (n = 17), Tankwa (n = 15) and South African village (n = 35) goat populations. A range of 27–58 SNPs per population were observed. Mutations resulting in amino acid changes were observed at exons 2 and 5. Higher within‐breed diversity of 97.37% was observed within the population category consisting of SA village ecotypes and the Tankwa goats. Highest pairwise F_ST values ranging from 0.148 to 0.356 were observed between the SAB and both the South African village and Tankwa feral goat populations. Phylogenetic analysis indicated nine genetic clusters, which reflected close relationships between the South African populations and the other international breeds with the exception of the Italian Sarda breeds. Results imply greater potential for within‐population selection programs, particularly with SA village goats.     

Two missense mutations in exon 9 of caprine PRLR gene were associated with litter size

Guanzhong (n = 321) and Boer (n = 191) goat breeds were used to detect single nucleotide polymorphisms (SNPs) in the coding regions of the prolactin receptor (PRLR) gene by DNA sequencing and PCR‐RFLP. Two SNPs (c.1457G>A and c.1645G>A) were identified that caused amino acid variations p.Ser485Asn and p.Val548Met respectively. Statistical results indicated that the c.1457G>A and c.1645G>A SNPs were significantly associated with litter size in Boer and Guanzhong goat breeds. Further analysis revealed that combined genotype C4 (GGGG) and haplotype G‐G were better than the others for litter size in both goat breeds. These results might contribute to goat genetic resources and breeding.     

Evaluation of reduced subsets of single nucleotide polymorphisms for the prediction of age at puberty in sows

Genomic information could be used efficiently to improve traits that are expensive to measure, sex limited or expressed late in life. This study analyzed the phenotypic variation explained by major SNPs and windows for age at puberty in gilts, an indicator of reproductive longevity. A genome‐wide association study using 56 424 SNPs explained 25.2% of the phenotypic variation in age at puberty in a training set (n = 820). All SNPs from the top 10% of 1‐Mb windows explained 33.5% of the phenotypic variance compared to 47.1% explained by the most informative markers (n = 261). In an evaluation population, consisting of subsequent batches (n = 412), the predictive ability of all SNPs from the major 1‐Mb windows was higher compared to the variance captured by the most informative SNP from each of these windows. The phenotypic variance explained in the evaluation population varied from 12.3% to 36.8% when all SNPs from major windows were used compared to 6.5–23.7% explained by most informative SNPs. The correlation between phenotype and genomic prediction values based on SNP effects estimated in the training population was marginal compared to their effects retrained in the evaluation population for all (0.46–0.81) or most informative SNPs (0.30–0.65) from major windows. An increase in genetic gain of 20.5% could be obtained if genomic selection included both sexes compared to females alone. The pleiotropic role of major genes such as AVPR1A could be exploited in selection of both age at puberty and reproductive longevity.     

The adaptive potential of Populus balsamifera L. to phenology requirements in a warmer global climate

The manner in which organisms adapt to climate change informs a broader understanding of the evolution of biodiversity as well as conservation and mitigation plans. We apply common garden and association mapping approaches to quantify genetic variance and identify loci affecting bud flush and bud set, traits that define a tree's season for height growth, in the boreal forest tree Populus balsamifera L. (balsam poplar). Using data from 478 genotypes grown in each of two common gardens, one near the southern edge and another near the northern edge of P. balsamifera's range, we found that broad‐sense heritability for bud flush and bud set was generally high (H² > 0.5 in most cases), suggesting that abundant genetic variation exists for phenological response to changes in the length of the growing season. To identify the molecular genetic basis of this variation, we genotyped trees for 346 candidate single nucleotide polymorphisms (SNPs) from 27 candidate genes for the CO/FT pathway in poplar. Mixed‐model analyses of variance identified SNPs in 10 genes to be associated with variation in either bud flush or bud set. Multiple SNPs within FRIGIDA were associated with bud flush, whereas multiple SNPs in LEAFY and GIGANTEA 5 were associated with bud set. Although there was strong population structure in stem phenology, the geographic distribution of multilocus association SNP genotypes was widespread except at the most northern populations, indicating that geographic regions may harbour sufficient diversity in functional genes to facilitate adaption to future climatic conditions in many sites.     

Stable carbon isotope analysis reveals widespread drought stress in boreal black spruce forests

Unprecedented rates of climate warming over the past century have resulted in increased forest stress and mortality worldwide. Decreased tree growth in association with increasing temperatures is generally accepted as a signal of temperature‐induced drought stress. However, variations in tree growth alone do not reveal the physiological mechanisms behind recent changes in tree growth. Examining stable carbon isotope composition of tree rings in addition to tree growth can provide a secondary line of evidence for physiological drought stress. In this study, we examined patterns of black spruce growth and carbon isotopic composition in tree rings in response to climate warming and drying in the boreal forest of interior Alaska. We examined trees at three nested scales: landscape, toposequence, and a subsample of trees within the toposequence. At each scale, we studied the potential effects of differences in microclimate and moisture availability by sampling on northern and southern aspects. We found that black spruce radial growth responded negatively to monthly metrics of temperature at all examined scales, and we examined ?¹³C responses on a subsample of trees as representative of the wider region. The negative ?¹³C responses to temperature reveal that black spruce trees are experiencing moisture stress on both northern and southern aspects. Contrary to our expectations, ?¹³C from trees on the northern aspect exhibited the strongest drought signal. Our results highlight the prominence of drought stress in the boreal forest of interior Alaska. We conclude that if temperatures continue to warm, we can expect drought‐induced productivity declines across large regions of the boreal forest, even for trees located in cool and moist landscape positions.     

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.     

Quantifying the effects of drought on abrupt growth decreases of major tree species in Switzerland

Drought entails important effects on tree physiology, which may result in short‐ to long‐term radial growth decreases. While the majority of studies have focused on annual drought‐related variability of growth, relatively little is known about sustained growth decreases following drought years. We apply a statistical framework to identify climatic factors that induce abrupt growth decreases and may eventually result in tree mortality. We used tree‐ring data from almost 500 standing dead trees and 200 living trees in eight sites of the Swiss network of strict forest reserves, including four of the most important Central European tree species (Abies alba, Picea abies, Fagus sylvatica and Quercus spp.). First, to assess short‐term growth responses to drought under various climate and site conditions, we calculated correlations and linear mixed‐effects models between ring‐width indices (RWIs) and drought based on the Standardized Precipitation Evapotranspiration Index (SPEI). Second, to quantify drought effects on abrupt growth decreases, we applied distributed lag nonlinear models (DLNMs), which account for both delayed effects and the nonlinear relationship between the SPEI and the occurrence of abrupt growth decreases. Positive correlations between RWIs and the SPEI indicated short‐term growth responses of all species, particularly at arid sites. Results of the DLNMs revealed species‐specific growth responses to drought. For Quercus spp., abrupt growth decreases were more likely to occur several years following severe drought, whereas for P. abies, A. alba, and F. sylvatica abrupt growth decreases started frequently immediately in the drought year. We conclude that the statistical framework allows for quantifying the effects of drought intensity on the probability of abrupt growth decreases, which ultimately contributes to an improved understanding of climate impacts on forest community dynamics.     

Use of intraspecific variation in thermal responses for estimating an elevational cline in the timing of cold hardening in a sub‐boreal conifer

To avoid winter frost damage, evergreen coniferous species develop cold hardiness with suitable phenology for the local climate regime. Along the elevational gradient, a genetic cline in autumn phenology is often recognised among coniferous populations, but further quantification of evolutionary adaptation related to the local environment and its responsible signals generating the phenological variation are poorly understood. We evaluated the timing of cold hardening among populations of Abies sachalinensis, based on time series freezing tests using trees derived from four seed source populations × three planting sites. Furthermore, we constructed a model to estimate the development of hardening from field temperatures and the intraspecific variations occurring during this process. An elevational cline was detected such that high‐elevation populations developed cold hardiness earlier than low‐elevation populations, representing significant genetic control. Because development occurred earlier at high‐elevation planting sites, the genetic trend across elevation overlapped with the environmental trend. Based on the trade‐off between later hardening to lengthen the active growth period and earlier hardening to avoid frost damage, this genetic cline would be adaptive to the local climate. Our modelling approach estimated intraspecific variation in two model components: the threshold temperature, which was the criterion for determining whether the trees accumulated the thermal value, and the chilling requirement for trees to achieve adequate cold hardiness. A higher threshold temperature and a lower chilling requirement could be responsible for the earlier phenology of the high‐elevation population. These thermal responses may be one of the important factors driving the elevation‐dependent adaptation of A. sachalinensis.     

New resources for genetic studies in Populus nigra: genome‐wide SNP discovery and development of a 12k Infinium array

Whole genome resequencing of 51 Populus nigra (L.) individuals from across Western Europe was performed using Illumina platforms. A total number of 1 878 727 SNPs distributed along the P. nigra reference sequence were identified. The SNP calling accuracy was validated with Sanger sequencing. SNPs were selected within 14 previously identified QTL regions, 2916 expressional candidate genes related to rust resistance, wood properties, water‐use efficiency and bud phenology and 1732 genes randomly spread across the genome. Over 10 000 SNPs were selected for the construction of a 12k Infinium Bead‐Chip array dedicated to association mapping. The SNP genotyping assay was performed with 888 P. nigra individuals. The genotyping success rate was 91%. Our high success rate was due to the discovery panel design and the stringent parameters applied for SNP calling and selection. In the same set of P. nigra genotypes, linkage disequilibrium throughout the genome decayed on average within 5–7 kb to half of its maximum value. As an application test, ADMIXTURE analysis was performed with a selection of 600 SNPs spread throughout the genome and 706 individuals collected along 12 river basins. The admixture pattern was consistent with genetic diversity revealed by neutral markers and the geographical distribution of the populations. These newly developed SNP resources and genotyping array provide a valuable tool for population genetic studies and identification of QTLs through natural‐population based genetic association studies in P. nigra.     

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.     

An ecological genetic delineation of local seed‐source provenance for ecological restoration

An increasingly important practical application of the analysis of spatial genetic structure within plant species is to help define the extent of local provenance seed collection zones that minimize negative impacts in ecological restoration programs. Here, we derive seed sourcing guidelines from a novel range‐wide assessment of spatial genetic structure of 24 populations of Banksia menziesii (Proteaceae), a widely distributed Western Australian tree of significance in local ecological restoration programs. An analysis of molecular variance (AMOVA) of 100 amplified fragment length polymorphism (AFLP) markers revealed significant genetic differentiation among populations (Φ_PT = 0.18). Pairwise population genetic dissimilarity was correlated with geographic distance, but not environmental distance derived from 15 climate variables, suggesting overall neutrality of these markers with regard to these climate variables. Nevertheless, Bayesian outlier analysis identified four markers potentially under selection, although these were not correlated with the climate variables. We calculated a global R‐statistic using analysis of similarities (ANOSIM) to test the statistical significance of population differentiation and to infer a threshold seed collection zone distance of ~60 km (all markers) and 100 km (outlier markers) when genetic distance was regressed against geographic distance. Population pairs separated by >60 km were, on average, twice as likely to be significantly genetically differentiated than population pairs separated by <60 km, suggesting that habitat‐matched sites within a 30‐km radius around a restoration site genetically defines a local provenance seed collection zone for B. menziesii. Our approach is a novel probability‐based practical solution for the delineation of a local seed collection zone to minimize negative genetic impacts in ecological restoration.     

Pollen,wind and fire: how to investigate genetic effects of disturbance‐induced change in forest trees

Understanding the consequences of habitat disturbance on mating patterns although pollen and seed dispersal in forest trees has been a long‐standing theme of forest and conservation genetics. Forest ecosystems face global environmental pressures from timber exploitation to genetic pollution and climate change, and it is therefore essential to comprehend how disturbances may alter the dispersal of genes and their establishment in tree populations in order to formulate relevant recommendations for sustainable resource management practices and realistic predictions of potential adaptation to climate change by means of range shift or expansion (Kremer et al. 2012 ). However, obtaining reliable evidence of disturbance‐induced effects on gene dispersal processes from empirical evaluation of forest tree populations is difficult. Indeed, tree species share characteristics such as high longevity, long generation time and large reproductive population size, which may impede the experimenter's ability to assess parameters at the spatial and time scales at which any change may occur (Petit and Hampe 2006 ). It has been suggested that appropriate study designs should encompass comparison of populations before and after disturbance as well as account for demonstrated variation in conspecific density, that is, the spatial distribution of mates, and forest density, including all species and relating to alteration in landscape openness (Bacles & Jump 2011 ). However, more often than not, empirical studies aiming to assess the consequences of habitat disturbance on genetic processes in tree populations assume rather than quantify a change in tree densities in forests under disturbance and generally fail to account for population history, which may lead to inappropriate interpretation of a causal relationship between population genetic structure and habitat disturbance due to effects of unmonitored confounding variables (Gauzere et al. 2013). In this issue, Shohami and Nathan ( 2014 ) take advantage of the distinctive features of the fire‐adapted wind‐pollinated Aleppo pine Pinus halepensis (Fig. 1) to provide an elegant example of best practice. Thanks to long‐term monitoring of the study site, a natural stand in Israel, Shohami and Nathan witnessed the direct impact of habitat disturbance, here taking the shape of fire, on conspecific and forest densities and compared pre‐ and postdisturbance mating patterns estimated from cones of different ages sampled on the same surviving maternal individuals (Fig. 2). This excellent study design is all the more strong that Shohami and Nathan took further analytical steps to account for confounding variables, such as historical population genetic structure and possible interannual variation in wind conditions, thus giving high credibility to their findings of unequivocal fire‐induced alteration of mating patterns in P. halepensis. Most notably, the authors found, at the pollen pool level, a disruption of local genetic structure which, furthermore, they were able to attribute explicitly to enhanced pollen‐mediated gene immigration into the low‐density fire‐disturbed stand. This cleverly designed research provides a model approach to be followed if we are to advance our understanding of disturbance‐induced dispersal and genetic change in forest trees.     

The capacity to cope with climate warming declines from temperate to tropical latitudes in two widely distributed Eucalyptus species

As rapid climate warming creates a mismatch between forest trees and their home environment, the ability of trees to cope with warming depends on their capacity to physiologically adjust to higher temperatures. In widespread species, individual trees in cooler home climates are hypothesized to more successfully acclimate to warming than their counterparts in warmer climates that may approach thermal limits. We tested this prediction with a climate‐shift experiment in widely distributed Eucalyptus tereticornis and E. grandis using provenances originating along a ~2500 km latitudinal transect (15.5–38.0°S) in eastern Australia. We grew 21 provenances in conditions approximating summer temperatures at seed origin and warmed temperatures (+3.5 °C) using a series of climate‐controlled glasshouse bays. The effects of +3.5 °C warming strongly depended on home climate. Cool‐origin provenances responded to warming through an increase in photosynthetic capacity and total leaf area, leading to enhanced growth of 20–60%. Warm‐origin provenances, however, responded to warming through a reduction in photosynthetic capacity and total leaf area, leading to reduced growth of approximately 10%. These results suggest that there is predictable intraspecific variation in the capacity of trees to respond to warming; cool‐origin taxa are likely to benefit from warming, while warm‐origin taxa may be negatively affected.     

Minimum sample sizes for population genomics: an empirical study from an Amazonian plant species

High‐throughput DNA sequencing facilitates the analysis of large portions of the genome in nonmodel organisms, ensuring high accuracy of population genetic parameters. However, empirical studies evaluating the appropriate sample size for these kinds of studies are still scarce. In this study, we use double‐digest restriction‐associated DNA sequencing (ddRADseq) to recover thousands of single nucleotide polymorphisms (SNPs) for two physically isolated populations of Amphirrhox longifolia (Violaceae), a nonmodel plant species for which no reference genome is available. We used resampling techniques to construct simulated populations with a random subset of individuals and SNPs to determine how many individuals and biallelic markers should be sampled for accurate estimates of intra‐ and interpopulation genetic diversity. We identified 3646 and 4900 polymorphic SNPs for the two populations of A. longifolia, respectively. Our simulations show that, overall, a sample size greater than eight individuals has little impact on estimates of genetic diversity within A. longifolia populations, when 1000 SNPs or higher are used. Our results also show that even at a very small sample size (i.e. two individuals), accurate estimates of F_ST can be obtained with a large number of SNPs (≥1500). These results highlight the potential of high‐throughput genomic sequencing approaches to address questions related to evolutionary biology in nonmodel organisms. Furthermore, our findings also provide insights into the optimization of sampling strategies in the era of population genomics.