Scots pine - panmixia and the elusive signal of genetic adaptation

Scots pine is the foundation species of diverse forested ecosystems across Eurasia and displays remarkable ecological breadth, occurring in environments ranging from temperate rainforests to arid tundra margins. Such expansive distributions can be favored by various demographic and adaptive processes and the interactions between them. To understand the impact of neutral and selective forces on genetic structure in Scots pine, we conducted range-wide population genetic analyses on 2321 trees from 202 populations using genotyping-by-sequencing, reconstructed the recent demography of the species and examined signals of genetic adaptation. We found a high and uniform genetic diversity across the entire range (global FST 0.048), no increased genetic load in expanding populations and minor impact of the last glacial maximum on historical population sizes. Genetic-environmental associations identified only a handful of single-nucleotide polymorphisms significantly linked to environmental gradients. The results suggest that extensive gene flow is predominantly responsible for the observed genetic patterns in Scots pine. The apparent missing signal of genetic adaptation is likely attributed to the intricate genetic architecture controlling adaptation to multi-dimensional environments. The panmixia metapopulation of Scots pine offers a good study system for further exploration into how genetic adaptation and plasticity evolve under gene flow and changing environment.

The GenTree Platform: growth traits and tree-level environmental data in 12 European forest tree species

BACKGROUND

Progress in the field of evolutionary forest ecology has been hampered by the huge challenge of phenotyping trees across their ranges in their natural environments, and the limitation in high-resolution environmental information.

FINDINGS

The GenTree Platform contains phenotypic and environmental data from 4,959 trees from 12 ecologically and economically important European forest tree species: Abies alba Mill. (silver fir), Betula pendula Roth. (silver birch), Fagus sylvatica L. (European beech), Picea abies (L.) H. Karst (Norway spruce), Pinus cembra L. (Swiss stone pine), Pinus halepensis Mill. (Aleppo pine), Pinus nigra Arnold (European black pine), Pinus pinaster Aiton (maritime pine), Pinus sylvestris L. (Scots pine), Populus nigra L. (European black poplar), Taxus baccata L. (English yew), and Quercus petraea (Matt.) Liebl. (sessile oak). Phenotypic (height, diameter at breast height, crown size, bark thickness, biomass, straightness, forking, branch angle, fructification), regeneration, environmental in situ measurements (soil depth, vegetation cover, competition indices), and environmental modeling data extracted by using bilinear interpolation accounting for surrounding conditions of each tree (precipitation, temperature, insolation, drought indices) were obtained from trees in 194 sites covering the species' geographic ranges and reflecting local environmental gradients.

CONCLUSION

The GenTree Platform is a new resource for investigating ecological and evolutionary processes in forest trees. The coherent phenotyping and environmental characterization across 12 species in their European ranges allow for a wide range of analyses from forest ecologists, conservationists, and macro-ecologists. Also, the data here presented can be linked to the GenTree Dendroecological collection, the GenTree Leaf Trait collection, and the GenTree Genomic collection presented elsewhere, which together build the largest evolutionary forest ecology data collection available.

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.

The extent and meaning of hybridization and introgression between Siberian spruce (Picea obovata) and Norway spruce (Picea abies): cryptic refugia as stepping stones to the west?

Boreal species were repeatedly exposed to ice ages and went through cycles of contraction and expansion while sister species alternated periods of contact and isolation. The resulting genetic structure is consequently complex, and demographic inferences are intrinsically challenging. The range of Norway spruce (Picea abies) and Siberian spruce (Picea obovata) covers most of northern Eurasia; yet their geographical limits and histories remain poorly understood. To delineate the hybrid zone between the two species and reconstruct their joint demographic history, we analysed variation at nuclear SSR and mitochondrial DNA in 102 and 88 populations, respectively. The dynamics of the hybrid zone was analysed with approximate Bayesian computation (ABC) followed by posterior predictive structure plot reconstruction and the presence of barriers across the range tested with estimated effective migration surfaces. To estimate the divergence time between the two species, nuclear sequences from two well-separated populations of each species were analysed with ABC. Two main barriers divide the range of the two species: one corresponds to the hybrid zone between them, and the other separates the southern and northern domains of Norway spruce. The hybrid zone is centred on the Urals, but the genetic impact of Siberian spruce extends further west. The joint distribution of mitochondrial and nuclear variation indicates an introgression of mitochondrial DNA from Norway spruce into Siberian spruce. Overall, our data reveal a demographic history where the two species interacted frequently and where migrants originating from the Urals and the West Siberian Plain recolonized northern Russia and Scandinavia using scattered refugial populations of Norway spruce as stepping stones towards the west.

Absence of geographical structure of chloroplast DNA variation in sallow, Salix caprea L

In the present study, we have used PCR-RFLP markers to investigate the chloroplast DNA variation in 24 European populations of Salix caprea L. A subset of these populations has also been analysed with chloroplast microsatellites. The main feature of both markers is the absence of a clear geographic structure (G(ST(PCR-RFLP))=0.090, G(ST(microsatellites))=-0.017) and high levels of variation within populations. This lack of phylogeographic structure in S. caprea is suggested to be the consequence of the joint action of several factors: (i) presence of intermediate latitude refugia with large population sizes during the last glacial maximum, (ii) a high speed of recolonisation and dispersal ability, (iii) a high mutation rate and (iv) extensive hybridisation with other willow species. In addition to the S. caprea samples, a limited number of individuals from several other Salix species were also analysed with PCR-RFLP: S. cinerea, S. aurita, S. purpurea, S. atrocinerea, S. appendiculata, S. elaeagnos, S. fragilis and S. alba. Many of the haplotypes found in Salix caprea were also detected in S. cinerea, S. aurita, S. purpurea, S. atrocinerea and/ or S. appendiculata but not in S. alba, S. elaeagnos or S. fragilis. Our data suggest that hybridisation and gene flow have occurred within these two groups but not between them.

Genetic mapping of sex-linked markers in Salix viminalis L

A total of 88 selective primer combinations were screened using bulked males and females sampled from four families of Salix viminalis. A total of more than 1000 polymorphic fragments was obtained, of which only four cosegregated with sex. These four sex-linked markers were subsequently scored in individuals that were used for bulked sample preparation in additional individuals of the same families, and in individuals in other families. A pair of primers that amplified the sex-linked fragments was constructed from one of the sex-linked amplified fragment length polymorphism (AFLP) fragments. In hybridization of Southern blot filters with the sex-linked DNA fragments, the band was present in females and absent in males, but the opposite pattern of band segregation (a band found in males and no band in females) was never observed in either the AFLP or RFLP experiments. Two of the sex-linked markers were placed on a linkage map. They both map at the same location in a linkage group comprising other markers not segregating with sex. Our data suggest that a single locus governs the sex determination and that nonrecombining sex chromosomes are absent in S. viminalis. A close association was found between skewed sex ratio and segregation distortion at this locus.