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.

Transcriptomic Analysis Reveals Novel Regulators of the Scots Pine Stilbene Pathway.

Stilbenes accumulate in Scots pine heartwood where they have important roles in protecting wood from decaying fungi. They are also part of active defense responses, and their production is induced by different (a)biotic stressors. The specific transcriptional regulators as well as the enzyme responsible for activating the stilbene precursor cinnamate in the pathway are still unknown. UV-C radiation was the first discovered artificial stress activator of the pathway. Here, we describe a large-scale transcriptomic analysis of pine needles in response to UV-C and treatment with translational inhibitors, both activating the transcription of stilbene pathway genes. We used the data to identify putative candidates for the missing CoA ligase and for pathway regulators. We further showed that the pathway is transcriptionally activated by phosphatase inhibitor, ethylene and jasmonate treatments, as in grapevine, and that the stilbene synthase promoter retains its inducibility in some of the tested conditions in Arabidopsis, a species that normally does not synthesize stilbenes. Shared features between gymnosperm and angiosperm regulation and partially retained inducibility in Arabidopsis suggest that pathway regulation occurs not only via ancient stress-response pathway(s) but also via species-specific regulators. Understanding which genes control the biosynthesis of stilbenes in Scots pine aids breeding of more resistant trees.

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.

A transcriptomic view to wounding response in young Scots pine stems.

We studied the stress response of five-year-old Scots pine xylem to mechanical wounding using RNA sequencing. In general, we observed a bimodal response in pine xylem after wounding. Transcripts associated with water deficit stress, defence, and cell wall modification were induced at the earliest time point of three hours; at the same time, growth-related processes were down-regulated. A second temporal wave was triggered either at the middle and/or at the late time points (one and four days). Secondary metabolism, such as stilbene and lignan biosynthesis started one day after wounding. Scots pine synthesises the stilbenes pinosylvin and its monomethyl ether both as constitutive and induced defence compounds. Stilbene biosynthesis is induced by wounding, pathogens and UV stress, but is also developmentally regulated when heartwood is formed. Comparison of wounding responses to heartwood formation shows that many induced processes (in addition to stilbene biosynthesis) are similar and relate to defence or desiccation stress, but often specific transcripts are up-regulated in the developmental and wounding induced contexts. Pine resin biosynthesis was not induced in response to wounding, at least not during the first four days.

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.

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.

Developmental Changes in Scots Pine Transcriptome during Heartwood Formation.

Scots pine (Pinus sylvestris L.) wood is desired in woodworking industries due to its favorable timber characteristics and natural durability that is contributed by heartwood extractives. It has been discussed whether the Scots pine heartwood extractives (mainly stilbenes and resin acids) are synthesized in the cells of the transition zone between sapwood and heartwood, or if they are transported from the sapwood. Timing of heartwood formation during the yearly cycle has also not been unambiguously defined. We measured steady-state mRNA levels in Scots pine transition zone and sapwood using RNA sequencing. Year-round expression profiles of selected transcripts were further investigated by quantitative RT-PCR. Differentially accumulating transcripts suggest that, of the Scots pine heartwood extractives, stilbenes are synthesized in situ in the transition zone and gain their carbon-skeletons from Suc and triglycerides. Resin acids, on the other hand, are synthesized early in the spring mainly in the sapwood, meaning that they must be transported to the heartwood transition zone. Heartwood formation is marked by programmed cell death that occurs during the summer months in the transition zone.

Association of FLOWERING LOCUS T/TERMINAL FLOWER 1-like gene FTL2 expression with growth rhythm in Scots pine (Pinus sylvestris).

Understanding the genetic basis of the timing of bud set, an important trait in conifers, is relevant for adaptation and forestry practice. In common garden experiments, both Scots pine (Pinus sylvestris) and Norway spruce (Picea abies) show a latitudinal cline in the trait. We compared the regulation of their bud set biology by examining the expression of PsFTL2, a Pinus sylvestris homolog to PaFTL2, a FLOWERING LOCUS T/TERMINAL FLOWER 1 (FT/TFL1)-like gene, the expression levels of which have been found previously to be associated with the timing of bud set in Norway spruce. In a common garden study, we analyzed the relationship of bud phenology under natural and artificial photoperiods and the expression of PsFTL2 in a set of Scots pine populations from different latitudes. The expression of PsFTL2 increased in the needles preceding bud set and decreased during bud burst. In the northernmost population, even short night periods were efficient to trigger this expression, which also increased earlier under all photoperiodic regimes compared with the southern populations. Despite the different biology, with few limitations, the two conifers that diverged 140 million yr ago probably share an association of FTL2 with bud set, pointing to a common mechanism for the timing of growth cessation in conifers.

Clinal variation at phenology-related genes in spruce: parallel evolution in FTL2 and Gigantea?

Parallel clines in different species, or in different geographical regions of the same species, are an important source of information on the genetic basis of local adaptation. We recently detected latitudinal clines in SNPs frequencies and gene expression of candidate genes for growth cessation in Scandinavian populations of Norway spruce (Picea abies). Here we test whether the same clines are also present in Siberian spruce (P. obovata), a close relative of Norway spruce with a different Quaternary history. We sequenced nine candidate genes and 27 control loci and genotyped 14 SSR loci in six populations of P. obovata located along the Yenisei river from latitude 56°N to latitude 67°N. In contrast to Scandinavian Norway spruce that both departs from the standard neutral model (SNM) and shows a clear population structure, Siberian spruce populations along the Yenisei do not depart from the SNM and are genetically unstructured. Nonetheless, as in Norway spruce, growth cessation is significantly clinal. Polymorphisms in photoperiodic (FTL2) and circadian clock (Gigantea, GI, PRR3) genes also show significant clinal variation and/or evidence of local selection. In GI, one of the variants is the same as in Norway spruce. Finally, a strong cline in gene expression is observed for FTL2, but not for GI. These results, together with recent physiological studies, confirm the key role played by FTL2 and circadian clock genes in the control of growth cessation in spruce species and suggest the presence of parallel adaptation in these two species.

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.

Gene, phenotype and function: GLABROUS1 and resistance to herbivory in natural populations of Arabidopsis lyrata.

The molecular genetic basis of adaptive variation is of fundamental importance for evolutionary dynamics, but is still poorly known. Only in very few cases has the relationship between genetic variation at the molecular level, phenotype and function been established in natural populations. We examined the functional significance and genetic basis of a polymorphism in production of leaf hairs, trichomes, in the perennial herb Arabidopsis lyrata. Earlier studies suggested that trichome production is subject to divergent selection. Here we show that the production of trichomes is correlated with reduced damage from insect herbivores in natural populations, and using statistical methods developed for medical genetics we document an association between loss of trichome production and mutations in the regulatory gene GLABROUS1. Sequence data suggest that independent mutations in this regulatory gene have provided the basis for parallel evolution of reduced resistance to insect herbivores in different populations of A. lyrata and in the closely related Arabidopsis thaliana. The results show that candidate genes identified in model organisms provide a valuable starting point for analysis of the genetic basis of phenotypic variation in natural populations.

Population structure in Arabidopsis lyrata: evidence for divergent selection on trichome production.

Leaf trichomes may serve several biological functions including protection against herbivores, drought, and UV radiation; and their adaptive value can be expected to vary among environments. The perennial, self-incompatible herb Arabidopsis lyrata is polymorphic for trichome production, and occurs in a glabrous and a trichome-producing form. Controlled crosses indicate that the polymorphism is governed by a single gene, with trichome production being dominant. We examined the hypothesis that trichome production is subject to divergent selection (i.e., directional selection favoring different phenotypes in different populations) by comparing patterns of variation at the locus coding for glabrousness and at eight putatively neutral isozyme loci in Swedish populations of A. lyrata. The genetic diversity (He) and allele number at isozyme loci tended to increase with population size and decreased with latitude of origin, whereas genetic diversity at the locus coding for glabrousness did not vary with population size and increased with latitude of origin. The degree of genetic differentiation at the glabrousness locus was much higher than that at isozyme loci. Genetic identity at isozyme loci was negatively related to geographic distance, suggesting isolation by distance. In contrast, there was no significant correlation between genetic identity at the glabrousness locus and at isozyme loci. The results are consistent with the hypothesis that divergent selection contributes to population differentiation in trichome production in A. lyrata.

Genetic basis of trichome production in Arabidopsis lyrata.

Leaf trichomes may protect plants against herbivorous insects, and may increase tolerance to drought and UV-radiation. The perennial herb Arabidopsis lyrata (Brassicaceae) is polymorphic for trichome production and occurs in a glabrous and trichome-producing form. In addition, there is quantitative variation in trichome density among trichome-producing plants. To examine the genetic basis of glabrousness, we conducted controlled crosses with plants originating from two natural populations in Sweden (one polymorphic for trichome-production, and one consisting of glabrous plants only). In addition, we estimated the heritability of trichome number from parent-offspring regressions for plants originating from the polymorphic population. Crosses between glabrous plants resulted in glabrous offspring only, whereas crosses between glabrous and trichome-producing plants, and crosses between trichome-producing individuals, resulted in either all trichome-producing offspring or both phenotypes. In segregating crosses between trichome-producing plants, the ratio of glabrous:trichome-producing offspring did not deviate significantly from 1:3, while in segregating crosses between glabrous and trichome-producing individuals the ratio did in most cases not deviate from 1:1. Within- and between-population crosses gave similar results. The heritability of trichome number estimated from regression of offspring on mid-parent was high (h2 +/- SE, 0.65 +/- 0.15). The results suggest that glabrousness is inherited in a simple Mendelian fashion, with the allele coding for trichome production being dominant over that for glabrousness. They further indicate that glabrousness is due to a mutation at the same locus in both populations.

GENETIC BASIS OF INBREEDING DEPRESSION IN ARABIS PETRAEA.

Inbreeding depression may be caused by (partially) recessive or overdominant gene action. The relative evolutionary importance of these two modes has been debated; the former mode is emphasized in the "dominance hypothesis," the latter in the "overdominance hypothesis." We analyzed the genetic basis of inbreeding depression in the self-incompatible herb Arabis petraea (L.) Lam.: In the selfed progeny of twelve parental plants, we studied the proportion of chlorophyll-deficient seedlings, the genotypic distributions of marker genes, and associations of marker genotypes with viability and quantitative traits. Early components of fitness were examined by scoring seed size, germination time, and early growth rate and by observing the proportion of chlorophyll-deficient seedlings. Later components of fitness, flowering, and root and aboveground biomass were also measured. Marker genotypes of young seedlings were scored for 11 enzyme loci and three microsatellite markers. We found a high proportion (about 70%) of families with chlorophyll-deficient seedlings, indicating a high mutational load. We found six significant deviations from 1:2:1 ratio at marker loci of 60 tests in seedlings, with three of these significant at the experimentwide level. Deviations from the expected ratio were assumed to be due to linked viability loci. A graphical and a Bayesian method were used to distinguish between the overdominance and dominance hypotheses. Most of the deviant segregation ratios suggested overdominance instead of recessivity of the deleterious allele. Neither the early (seed size, germination time, or early growth trait) nor the late quantitative traits (flowering, and root and aboveground biomass) showed significant linkage to markers at the experimentwide level. Presence of significant associations between markers and early viability, but lack thereof for quantitative traits expressed late, suggests either that there may be relatively low inbreeding depression in later life stages or that individual quantitative trait loci may have smaller effects than loci contributing to early viability.

GEOGRAPHICAL VARIATION IN THE INBREEDING DEPRESSION OF SCOTS PINE.

The magnitude of inbreeding depression caused by recessive mutations in a population is dependent on the mutation rate and on the intensity of selection against the mutations. We studied geographical differences in the level of early inbreeding depression of Scots pine in a common garden experiment. The mean abortion rate of experimentally self-pollinated seeds was significantly lower (75.4%) among trees that originated from northern populations (66-69°N) than among trees from more southern (60-62°N) populations (86.5%). Thus, the number of embryonic lethal equivalents was lower in the northern populations (4.5) than in the southern ones (6.9). The outcrossing rate at the mature seed stage was slightly lower in the northern populations (average 0.93) than in the southern one (0.99). The estimated selfing rate at the zygote stage varied from 0-0.28 in the populations. The reduction in the magnitude of inbreeding depression in the north may have been caused by increased levels of self-fertilization in the northern populations. The proportion of self-fertilized seedlings and adults was very small in all populations (F ≈ 0), indicating high inbreeding depression also in later life stages. The high level of inbreeding depression in the partially selfing Scots pine can be explained by mutation-selection balance only if the mutation rate is high.