Assessing the adaptative potential to temperature and precipitation along a steep environmental gradient in populations of European beech

Climate change poses a significant threat to European beech. These concerns highlight the need to assess the adaptive potential of European beech populations to climate change. Landscape genomics, also known as environmental association analysis, is a powerful tool for identifying gene loci that contribute to local adaptation to environmental pressures. Genotypic data was collected from [~]100 adult beech trees per stand in five locations in the South-Eastern Romanian Carpathians along an altitudinal gradient associated with precipitation and temperature. In total, 53 environmental variables, comprising frost frequency change, temperature and precipitation, were extracted from the climatology data base CHELSA. Based on these variables the Ellenberg-Quotient (EQ) was calculated. We performed environmental association analysis using LFMM (latent factor mixed models) to identify Single Nucleotide Polymorphism (SNP) markers associated with environmental variables and with the principal components calculated based on these. We identified 446 SNP markers significantly associated with the first principal component (PC). These were overlapping with the SNP markers significantly associated with all environmental variables except precipitation accumulated during the growing season. The first PC was correlated with all temperature-based variables and elevation at |r| [~]0.989 to [~]0.997 and with all precipitation-based and Ellenberg-Quotient variables at |r| [~]0.945 to 0.950, except precipitation accumulated during the growing season. A high peak region on chromosome 2 from [~]4.56 to [~]16.27 Mb appeared in all results. This region was [~]3.47 Mb downstream from a region for local adaptation identified by Lazic et al. (2024). In this peak, 273 markers located in the coding region of 22 genes were found. Ten out these 22 were described based on a literature review. Among these ten genes, two may be involved in local adaptation based on our literature review. These two genes are polygalacturonase QRT3-like and NRT1/PTR_FAMILY 5.4-like. The gene polygalacturonase QRT3-like plays a role in pollen development in Arabidopsis thaliana L. and Brassica rapa L. We observed at the corresponding SNP markers, a correlation of the minor allele frequency and temperature-based environmental variables.