Within- and between-population comparisons suggest independently acting selection maintaining parallel clines in Scots pine (Pinus sylvestris).

Parallel clines in traits related to adaptation in a species can be due to independent selection on a pair of traits, or due to selection in one trait resulting in a parallel cline in a correlated trait. To distinguish between the mechanisms giving rise to parallel adaptive population divergence of multiple traits along an environmental gradient we need to study variation, correlations, and selective forces within individual populations along the gradient. In many tree species, budset timing (BST) forms a latitudinal cline, and parallel clinal variation is also found in other seedling traits, such as first-year height (FYH) and fall frost injury (FFI). In this study, we set up a common garden experiment with open pollinated progeny from natural populations of Scots pine (Pinus sylvestris), with one large sample from single population (500 families) and smaller samples from across a latitudinal gradient. BST, FYH and induced FFI were first measured in a greenhouse. The seedlings were then planted in the field, where survival and height were measured at the age of 9 years as fitness proxies. We compared between- and within-population variation and genetic correlations of these three seedling traits, and estimated selection gradients at the family level in our main population, taking into account the potential effects of seed weight. Between-population genetic correlations between seedling traits were high (0.76-0.95). Within-population genetic correlations in the main population were lower (0.14-0.35), as in other populations (0.10-0.39). Within population, extensive adaptive variation persists in the seedling traits, in line with rather weak selection gradients, yet maintaining the clines. Although our sampling does not cover the whole cline equally, the results suggest that the individual clines in these traits are maintained by largely independently acting selection, which results in fewer constraints in adaptation under changing climate.

Taming the massive genome of Scots pine with PiSy50k, a new genotyping array for conifer research.

Pinus sylvestris (Scots pine) is the most widespread coniferous tree in the boreal forests of Eurasia, with major economic and ecological importance. However, its large and repetitive genome presents a challenge for conducting genome-wide analyses such as association studies, genetic mapping and genomic selection. We present a new 50K single-nucleotide polymorphism (SNP) genotyping array for Scots pine research, breeding and other applications. To select the SNP set, we first genotyped 480 Scots pine samples on a 407 540 SNP screening array and identified 47 712 high-quality SNPs for the final array (called 'PiSy50k'). Here, we provide details of the design and testing, as well as allele frequency estimates from the discovery panel, functional annotation, tissue-specific expression patterns and expression level information for the SNPs or corresponding genes, when available. We validated the performance of the PiSy50k array using samples from Finland and Scotland. Overall, 39 678 (83.2%) SNPs showed low error rates (mean = 0.9%). Relatedness estimates based on array genotypes were consistent with the expected pedigrees, and the level of Mendelian error was negligible. In addition, array genotypes successfully discriminate between Scots pine populations of Finnish and Scottish origins. The PiSy50k SNP array will be a valuable tool for a wide variety of future genetic studies and forestry applications.

Selection patterns on early-life phenotypic traits in Pinus sylvestris are associated with precipitation and temperature along a climatic gradient in Europe.

Understanding the dynamics of selection is key to predicting the response of tree species to new environmental conditions in the current context of climate change. However, selection patterns acting on early recruitment stages and their climatic drivers remain largely unknown in most tree species, despite being a critical period of their life cycle. We measured phenotypic selection on Pinus sylvestris seed mass, emergence time and early growth rate over 2 yr in four common garden experiments established along the latitudinal gradient of the species in Europe. Significant phenotypic plasticity and among-population genetic variation were found for all measured phenotypic traits. Heat and drought negatively affected fitness in the southern sites, but heavy rainfalls also decreased early survival in middle latitudes. Climate-driven directional selection was found for higher seed mass and earlier emergence time, while the form of selection on seedling growth rates differed among sites and populations. Evidence of adaptive and maladaptive phenotypic plasticity was found for emergence time and early growth rate, respectively. Seed mass, emergence time and early growth rate have an adaptive role in the early stages of P. sylvestris and climate strongly influences the patterns of selection on these fitness-related traits.

Author Correction: The GenTree Dendroecological Collection, tree-ring and wood density data from seven tree species across Europe.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

275 years of forestry meets genomics in Pinus sylvestris.

Pinus sylvestris has a long history of basic and applied research that is relevant for both forestry and evolutionary studies. Its patterns of adaptive variation and role in forest economic and ecological systems have been studied extensively for nearly 275 years, detailed demography for a 100 years and mating system more than 50 years. However, its reference genome sequence is not yet available and genomic studies have been lagging compared to, for example, Pinus taeda and Picea abies, two other economically important conifers. Despite the lack of reference genome, many modern genomic methods are applicable for a more detailed look at its biological characteristics. For example, RNA-seq has revealed a complex transcriptional landscape and targeted DNA sequencing displays an excess of rare variants and geographically homogenously distributed molecular genetic diversity. Current DNA and RNA resources can be used as a reference for gene expression studies, SNP discovery, and further targeted sequencing. In the future, specific consequences of the large genome size, such as functional effects of regulatory open chromatin regions and transposable elements, should be investigated more carefully. For forest breeding and long-term management purposes, genomic data can help in assessing the genetic basis of inbreeding depression and the application of genomic tools for genomic prediction and relatedness estimates. Given the challenges of breeding (long generation time, no easy vegetative propagation) and the economic importance, application of genomic tools has a potential to have a considerable impact. Here, we explore how genomic characteristics of P. sylvestris, such as rare alleles and the low extent of linkage disequilibrium, impact the applicability and power of the tools.

The GenTree Dendroecological Collection, tree-ring and wood density data from seven tree species across Europe.

The dataset presented here was collected by the GenTree project (EU-Horizon 2020), which aims to improve the use of forest genetic resources across Europe by better understanding how trees adapt to their local environment. This dataset of individual tree-core characteristics including ring-width series and whole-core wood density was collected for seven ecologically and economically important European tree species: silver birch (Betula pendula), European beech (Fagus sylvatica), Norway spruce (Picea abies), European black poplar (Populus nigra), maritime pine (Pinus pinaster), Scots pine (Pinus sylvestris), and sessile oak (Quercus petraea). Tree-ring width measurements were obtained from 3600 trees in 142 populations and whole-core wood density was measured for 3098 trees in 125 populations. This dataset covers most of the geographical and climatic range occupied by the selected species. The potential use of it will be highly valuable for assessing ecological and evolutionary responses to environmental conditions as well as for model development and parameterization, to predict adaptability under climate change scenarios.

Utilization of Tissue Ploidy Level Variation in de Novo Transcriptome Assembly of Pinus sylvestris.

Compared to angiosperms, gymnosperms lag behind in the availability of assembled and annotated genomes. Most genomic analyses in gymnosperms, especially conifer tree species, rely on the use of de novo assembled transcriptomes. However, the level of allelic redundancy and transcript fragmentation in these assembled transcriptomes, and their effect on downstream applications have not been fully investigated. Here, we assessed three assembly strategies for short-reads data, including the utility of haploid megagametophyte tissue during de novo assembly as single-allele guides, for six individuals and five different tissues in Pinus sylvestris We then contrasted haploid and diploid tissue genotype calls obtained from the assembled transcriptomes to evaluate the extent of paralog mapping. The use of the haploid tissue during assembly increased its completeness without reducing the number of assembled transcripts. Our results suggest that current strategies that rely on available genomic resources as guidance to minimize allelic redundancy are less effective than the application of strategies that cluster redundant assembled transcripts. The strategy yielding the lowest levels of allelic redundancy among the assembled transcriptomes assessed here was the generation of SuperTranscripts with Lace followed by CD-HIT clustering. However, we still observed some levels of heterozygosity (multiple gene fragments per transcript reflecting allelic redundancy) in this assembled transcriptome on the haploid tissue, indicating that further filtering is required before using these assemblies for downstream applications. We discuss the influence of allelic redundancy when these reference transcriptomes are used to select regions for probe design of exome capture baits and for estimation of population genetic diversity.

Adaptive potential of northernmost tree populations to climate change, with emphasis on Scots pine (Pinus sylvestris L.).

The adaptive potential of the northernmost Pinus sylvestris L. (and other northern tree) populations is considered by examining first the current patterns of quantitative genetic adaptive traits, which show high population differentiation and clines. We then consider the postglacial history of the populations using both paleobiological and genetic data. The current patterns of diversity at nuclear genes suggest that the traces of admixture are mostly visible in mitochondrial DNA variation patterns. There is little evidence of increased diversity due to admixture between an eastern and western colonization lineage, but no signal of reduced diversity (due to sequential bottlenecks) either. Quantitative trait variation in the north is not associated with the colonizing lineages. The current clines arose rapidly and may be based on standing genetic variation. The initial phenotypic response of Scots pine in the north is predicted to be increased survival and growth. The genetic responses are examined based on quantitative genetic predictions of sustained selection response and compared with earlier simulation results that have aimed at more ecological realism. The phenotypic responses of increased growth and survival reduce the opportunity for selection and delay the evolutionary responses. The lengthening of the thermal growing period also causes selection on the critical photoperiod in the different populations. Future studies should aim at including multiple ecological and genetic factors in evaluating potential responses.