Key triggers of adaptive genetic variability of sessile oak [Q. petraea (Matt.) Liebl.] from the Balkan refugia: outlier detection and association of SNP loci from ddRAD-seq data.

Knowledge on the genetic composition of Quercus petraea in south-eastern Europe is limited despite the species' significant role in the re-colonisation of Europe during the Holocene, and the diverse climate and physical geography of the region. Therefore, it is imperative to conduct research on adaptation in sessile oak to better understand its ecological significance in the region. While large sets of SNPs have been developed for the species, there is a continued need for smaller sets of SNPs that are highly informative about the possible adaptation to this varied landscape. By using double digest restriction site associated DNA sequencing data from our previous study, we mapped RAD-seq loci to the Quercus robur reference genome and identified a set of SNPs putatively related to drought stress-response. A total of 179 individuals from eighteen natural populations at sites covering heterogeneous climatic conditions in the southeastern natural distribution range of Q. petraea were genotyped. The detected highly polymorphic variant sites revealed three genetic clusters with a generally low level of genetic differentiation and balanced diversity among them but showed a north-southeast gradient. Selection tests showed nine outlier SNPs positioned in different functional regions. Genotype-environment association analysis of these markers yielded a total of 53 significant associations, explaining 2.4-16.6% of the total genetic variation. Our work exemplifies that adaptation to drought may be under natural selection in the examined Q. petraea populations.

SNP marker development in Pinus sylvestris L. in stress-responsive genes characterized from Pinus cembra L. transcriptomes.

Massively parallel sequencing of cDNA is an efficient route for generating sequence collections that represent expressed genes under different environmental control. The analysis of their sequence helps in developing molecular markers, such as SNPs, which represent a useful tool in detecting adaptive signals in populations. In this study novel PCR markers, based on stress responsive genes, were designed from the transcriptome of the haploxylon Swiss stone pine (Pinus cembra L.) and tested for SNPs in the diploxylon Scots pine (Pinus sylvestris L.). 84 primers were tested on P. sylvestris DNA samples originating from three different types of habitat. After sequencing and BLAST search of the amplified products, parts of 19 different candidate genes were analysed by considering the polymorphic sites, insertions/deletions as well as synonymous and non-synonymous SNPs. In a total of 3735 sites no indels, eight synonymous and 11 non-synonymous SNPs were found. By providing de novo molecular markers developed in P. cembra and tested for transferability in Scots pine, our results give support for the use of de novo markers targeting conserved regions across different pines. The SNPs detected may have important applications in further studies of adaptive genetic variation, providing tools to study relevant genes important in the long-term adaptation of pine species.

Geographic isolation and climatic variability contribute to genetic differentiation in fragmented populations of the long-lived subalpine conifer Pinus cembra L. in the western Alps.

BACKGROUND: Genetic processes shape the modern-day distribution of genetic variation within and between populations and can provide important insights into the underlying mechanisms of evolution. The resulting genetic variation is often unequally partitioned within species' distribution range and especially large differences can manifest at the range limit, where population fragmentation and isolation play a crucial role in species survival. Despite several molecular studies investigating the genetic diversity and differentiation of European Alpine mountain forests, the climatic and demographic constrains which influence the genetic processes are often unknown. Here, we apply non-coding microsatellite markers to evaluate the sporadic peripheral and continuous populations of cembra pine (Pinus cembra L.), a long-lived conifer species that inhabits the subalpine treeline ecotone in the western Alps to investigate how the genetic processes contribute to the modern-day spatial distribution. Moreover, we corroborate our findings with paleoecological records, micro and macro-remains, to infer the species' possible glacial refugia and expansion scenarios. RESULTS: Four genetically distinct groups were identified, with Bayesian and FST based approaches, across the range of the species, situated in the northern, inner and south-western Alps. We found that genetic differentiation is substantially higher in marginal populations than at the center of the range, and marginal stands are characterized by geographic and genetic isolation due to spatial segregation and restricted gene flow. Moreover, multiple matrix regression approaches revealed effects of climatic heterogeneity in species' spatial genetic pattern. Also, population stability tests indicated that all populations had experienced a severe historical bottleneck, no heterozygosity excess was detected, suggesting that more recently population sizes have remained relatively stable. CONCLUSIONS: Our study demonstrated that cembra pine might have survived in multiple glacial refugia and subsequently recolonized the Alps by different routes. Modern-day marginal populations, at the edge of the species' range, could maintain stable sizes over long periods without inbreeding depression and preserve high amounts of genetic variation. Moreover, our analyses indicate that climatic variability has played a major role in shaping differentiation, in addition to past historical events such as migration and demographic changes.