Genomic prediction of resistance to Hymenoscyphus fraxineus in common ash (Fraxinus excelsior L.) populations.

The increase in introduced insect pests and pathogens due to anthropogenic environmental changes has become a major concern for tree species worldwide. Common ash (Fraxinus excelsior L.) is one of such species facing a significant threat from the invasive fungal pathogen Hymenoscyphus fraxineus. Some studies have indicated that the susceptibility of ash to the pathogen is genetically determined, providing some hope for accelerated breeding programs that are aimed at increasing the resistance of ash populations. To address this challenge, we used a genomic selection strategy to identify potential genetic markers that are associated with resistance to the pathogen causing ash dieback. Through genome-wide association studies (GWAS) of 300 common ash individuals from 30 populations across Poland (ddRAD, dataset A), we identified six significant SNP loci with a p-value ≤1 × 10-4 associated with health status. To further evaluate the effectiveness of GWAS markers in predicting health status, we considered two genomic prediction scenarios. Firstly, we conducted cross-validation on dataset A. Secondly, we trained markers on dataset A and tested them on dataset B, which involved whole-genome sequencing of 20 individuals from two populations. Genomic prediction analysis revealed that the top SNPs identified via GWAS exhibited notably higher prediction accuracies compared to randomly selected SNPs, particularly with a larger number of SNPs. Cross-validation analyses using dataset A showcased high genomic prediction accuracy, predicting tree health status with over 90% accuracy across the top SNP sets ranging from 500 to 10,000 SNPs from the GWAS datasets. However, no significant results emerged for health status when the model trained on dataset A was tested on dataset B. Our findings illuminate potential genetic markers associated with resistance to ash dieback, offering support for future breeding programs in Poland aimed at combating ash dieback and bolstering conservation efforts for this invaluable tree species.

Genetic resources of common ash (Fraxinus excelsior L.) in Poland.

BACKGROUND: Knowledge of genetic structure and the factors that shape it has an impact on forest management practices. European ash (Fraxinus excelsior L.) has declined dramatically throughout its range as a result of a disease caused by the fungus Hymenoscyphus fraxineus. Despite the need for conservation and restoration of the species, genetic data required to guide these efforts at the country level are scarce. Thereofore, we studied the chloroplast and nuclear genetic diversity of 26 natural common ash populations (1269 trees) in Poland. RESULTS: Chloroplast polymorphisms grouped the populations into two geographically structured phylogenetic lineages ascribed to different glacial refugia (the Balkans and the Eastern Alps). However, the populations demonstrated high genetic diversity (mean AR = 12.35; mean Ho = 0.769; mean He = 0.542) but low differentiation based on nuclear microsatellites (FST = 0.045). Significant spatial genetic structure, consistent with models of isolation by distance, was detected in 14 out of 23 populations. Estimated effective population size was moderate-to-high, with a harmonic mean of 57.5 individuals per population. CONCLUSIONS: Genetic diversity was not homogeneously distributed among populations within phylogenetic gene pools, indicating that ash populations are not equal as potential sources of reproductive material. Genetic differences among populations could be related to their histories, including founder effects or gene flow between evolutionary lineages (admixture). Our results suggest that ash stands across Poland could be treated as two main management units (seed zones). Therefore, despite the homogenizing effect of pollen gene flow known for this species, the genetic structure should be taken into account in the management of the genetic resources of the common ash. Although ash dieback poses an additional challenge for the management of genetic resources, efforts should be directed towards protecting populations with high genetic diversity within defined phylogenetic units, as they may be an important source of adaptive variation for future stands.

Relatives reproduce in synchrony: kinship and individual condition shape intraspecific variation in masting phenotype.

Masting (synchronous and interannually variable seed production) is frequently called a reproductive strategy; yet it is unclear whether the reproductive behaviour of individuals has a heritable component. To address this, we used 22 years of annual fruit production data from 110 Sorbus aucuparia L. trees to examine the contributions of genetic factors to the reproductive phenotype of individuals, while controlling for environmental variation. Trees sharing close genetic relationships and experiencing similar habitat conditions exhibited similar levels of reproductive synchrony. Trees of comparable sizes displayed similar levels of year-to-year variation in fruiting, with relatedness contributing to this variation. External factors, such as shading, influenced the time intervals between years with abundant fruit production. The effects of genetic relatedness on the synchrony of reproduction among trees and on interannual variation provide long-awaited evidence that the masting phenotype is heritable, and can respond to natural selection.

Signatures of local adaptation to current and future climate in phenology-related genes in natural populations of Quercus robur.

BACKGROUND: Local adaptation is a key evolutionary process that enhances the growth of plants in their native habitat compared to non-native habitats, resulting in patterns of adaptive genetic variation across the entire geographic range of the species. The study of population adaptation to local environments and predicting their response to future climate change is important because of climate change. RESULTS: Here, we explored the genetic diversity of candidate genes associated with bud burst in pedunculate oak individuals sampled from 6 populations in Poland. Single nucleotide polymorphism (SNP) diversity was assessed in 720 candidate genes using the sequence capture technique, yielding 18,799 SNPs. Using landscape genomic approaches, we identified 8 FST outliers and 781 unique SNPs in 389 genes associated with geography, climate, and phenotypic variables (individual/family spring and autumn phenology, family diameter at breast height (DBH), height, and survival) that are potentially involved in local adaptation. Then, using a nonlinear multivariate model, Gradient Forests, we identified vulnerable areas of the pedunculate oak distribution in Poland that are at risk from climate change. CONCLUSIONS: The model revealed that pedunculate oak populations in the eastern part of the analyzed geographical region are the most sensitive to climate change. Our results might offer an initial evaluation of a potential management strategy for preserving the genetic diversity of pedunculate oak.

A novel synthesis of two decades of microsatellite studies on European beech reveals decreasing genetic diversity from glacial refugia.

Genetic diversity influences the evolutionary potential of forest trees under changing environmental conditions, thus indirectly the ecosystem services that forests provide. European beech (Fagus sylvatica L.) is a dominant European forest tree species that increasingly suffers from climate change-related die-back. Here, we conducted a systematic literature review of neutral genetic diversity in European beech and created a meta-data set of expected heterozygosity (He) from all past studies providing nuclear microsatellite data. We propose a novel approach, based on population genetic theory and a min-max scaling to make past studies comparable. Using a new microsatellite data set with unprecedented geographic coverage and various re-sampling schemes to mimic common sampling biases, we show the potential and limitations of the scaling approach. The scaled meta-dataset reveals the expected trend of decreasing genetic diversity from glacial refugia across the species range and also supports the hypothesis that different lineages met and admixed north of the European mountain ranges. As a result, we present a map of genetic diversity across the range of European beech which could help to identify seed source populations harboring greater diversity and guide sampling strategies for future genome-wide and functional investigations of genetic variation. Our approach illustrates how to combine information from several nuclear microsatellite data sets to describe patterns of genetic diversity extending beyond the geographic scale or mean number of loci used in each individual study, and thus is a proof-of-concept for synthesizing knowledge from existing studies also in other species. Supplementary Information: The online version contains supplementary material available at 10.1007/s11295-022-01577-4.

Phylogeny of Maleae (Rosaceae) Based on Complete Chloroplast Genomes Supports the Distinction of Aria, Chamaemespilus and Torminalis as Separate Genera, Different from Sorbus sp.

Several genera formerly contained within the genus Sorbus L. sensu lato have been proposed as separate taxa, including Aria, Chamaemespilus and Torminalis. However, molecular evidence for such distinctions are rather scarce. We assembled the complete chloroplast genome of Sorbus aucuparia, another representative of Sorbus s.s., and performed detailed comparisons with the available genomes of Aria edulis, Chamaemespilus alpina and Torminalis glaberrima. Additionally, using 110 complete chloroplast genomes of the Maleae representatives, we constructed the phylogenetic tree of the tribe using Maximum Likelihood methods. The chloroplast genome of S. aucuparia was found to be similar to other species within Maleae. The phylogenetic tree of the Maleae tribe indicated that A. edulis, C. alpina and T. glaberrima formed a concise group belonging to a different clade (related to Malus) than the one including Sorbus s.s. (related to Pyrus). However, Aria and Chamaemespilus appeared to be more closely related to each other than to Torminalis. Our results provide additional support for considering Aria, Chamaemespilus and Torminalis as separate genera different from Sorbus s.s.

Complete Chloroplast Genomes of Fagus sylvatica L. Reveal Sequence Conservation in the Inverted Repeat and the Presence of Allelic Variation in NUPTs.

Growing amounts of genomic data and more efficient assembly tools advance organelle genomics at an unprecedented scale. Genomic resources are increasingly used for phylogenetic analyses of many plant species, but are less frequently used to investigate within-species variability and phylogeography. In this study, we investigated genetic diversity of Fagus sylvatica, an important broadleaved tree species of European forests, based on complete chloroplast genomes of 18 individuals sampled widely across the species distribution. Our results confirm the hypothesis of a low cpDNA diversity in European beech. The chloroplast genome size was remarkably stable (158,428 ± 37 bp). The polymorphic markers, 12 microsatellites (SSR), four SNPs and one indel, were found only in the single copy regions, while inverted repeat regions were monomorphic both in terms of length and sequence, suggesting highly efficient suppression of mutation. The within-individual analysis of polymorphisms showed >9k of markers which were proportionally present in gene and non-gene areas. However, an investigation of the frequency of alternate alleles revealed that the source of this diversity originated likely from nuclear-encoded plastome remnants (NUPTs). Phylogeographic and Mantel correlation analysis based on the complete chloroplast genomes exhibited clustering of individuals according to geographic distance in the first distance class, suggesting that the novel markers and in particular the cpSSRs could provide a more detailed picture of beech population structure in Central Europe.

Genomic signatures of natural selection at phenology-related genes in a widely distributed tree species Fagus sylvatica L.

BACKGROUND: Diversity among phenology-related genes is predicted to be a contributing factor in local adaptations seen in widely distributed plant species that grow in climatically variable geographic areas, such as forest trees. European beech (Fagus sylvatica L.) is widespread, and is one of the most important broadleaved tree species in Europe; however, its potential for adaptation to climate change is a matter of uncertainty, and little is known about the molecular basis of climate change-relevant traits like bud burst. RESULTS: We explored single nucleotide polymorphisms (SNP) at candidate genes related to bud burst in beech individuals sampled across 47 populations from Europe. SNP diversity was monitored for 380 candidate genes using a sequence capture approach, providing 2909 unlinked SNP loci. We used two complementary analytical methods to find loci significantly associated with geographic variables, climatic variables (expressed as principal components), or phenotypic variables (spring and autumn phenology, height, survival). Redundancy analysis (RDA) was used to detect candidate markers across two spatial scales (entire study area and within subregions). We revealed 201 candidate SNPs at the broadest scale, 53.2% of which were associated with phenotypic variables. Additive polygenic scores, which provide a measure of the cumulative signal across significant candidate SNPs, were correlated with a climate variable (first principal component, PC1) related to temperature and precipitation availability, and spring phenology. However, different genotype-environment associations were identified within Southeastern Europe as compared to the entire geographic range of European beech. CONCLUSIONS: Environmental conditions play important roles as drivers of genetic diversity of phenology-related genes that could influence local adaptation in European beech. Selection in beech favors genotypes with earlier bud burst under warmer and wetter habitats within its range; however, selection pressures may differ across spatial scales.

A Chromosome-Level Genome Assembly of the European Beech (Fagus sylvatica) Reveals Anomalies for Organelle DNA Integration, Repeat Content and Distribution of SNPs.

The European Beech is the dominant climax tree in most regions of Central Europe and valued for its ecological versatility and hardwood timber. Even though a draft genome has been published recently, higher resolution is required for studying aspects of genome architecture and recombination. Here, we present a chromosome-level assembly of the more than 300 year-old reference individual, Bhaga, from the Kellerwald-Edersee National Park (Germany). Its nuclear genome of 541 Mb was resolved into 12 chromosomes varying in length between 28 and 73 Mb. Multiple nuclear insertions of parts of the chloroplast genome were observed, with one region on chromosome 11 spanning more than 2 Mb which fragments up to 54,784 bp long and covering the whole chloroplast genome were inserted randomly. Unlike in Arabidopsis thaliana, ribosomal cistrons are present in Fagus sylvatica only in four major regions, in line with FISH studies. On most assembled chromosomes, telomeric repeats were found at both ends, while centromeric repeats were found to be scattered throughout the genome apart from their main occurrence per chromosome. The genome-wide distribution of SNPs was evaluated using a second individual from Jamy Nature Reserve (Poland). SNPs, repeat elements and duplicated genes were unevenly distributed in the genomes, with one major anomaly on chromosome 4. The genome presented here adds to the available highly resolved plant genomes and we hope it will serve as a valuable basis for future research on genome architecture and for understanding the past and future of European Beech populations in a changing climate.

A reference genome of the European beech (Fagus sylvatica L.).

Background: The European beech is arguably the most important climax broad-leaved tree species in Central Europe, widely planted for its valuable wood. Here, we report the 542 Mb draft genome sequence of an up to 300-year-old individual (Bhaga) from an undisturbed stand in the Kellerwald-Edersee National Park in central Germany. Findings: Using a hybrid assembly approach, Illumina reads with short- and long-insert libraries, coupled with long Pacific Biosciences reads, we obtained an assembled genome size of 542 Mb, in line with flow cytometric genome size estimation. The largest scaffold was of 1.15 Mb, the N50 length was 145 kb, and the L50 count was 983. The assembly contained 0.12% of Ns. A Benchmarking with Universal Single-Copy Orthologs (BUSCO) analysis retrieved 94% complete BUSCO genes, well in the range of other high-quality draft genomes of trees. A total of 62,012 protein-coding genes were predicted, assisted by transcriptome sequencing. In addition, we are reporting an efficient method for extracting high-molecular-weight DNA from dormant buds, by which contamination by environmental bacteria and fungi was kept at a minimum. Conclusions: The assembled genome will be a valuable resource and reference for future population genomics studies on the evolution and past climate change adaptation of beech and will be helpful for identifying genes, e.g., involved in drought tolerance, in order to select and breed individuals to adapt forestry to climate change in Europe. A continuously updated genome browser and download page can be accessed from beechgenome.net, which will include future genome versions of the reference individual Bhaga, as new sequencing approaches develop.

Seeing the forest through the trees: comprehensive inference on individual mating patterns in a mixed stand of Quercus robur and Q. petraea.

BACKGROUND AND AIMS: Sexual reproduction is one of the most important moments in a life cycle, determining the genetic composition of individual offspring. Controlled pollination experiments often show high variation in the mating system at the individual level, suggesting a persistence of individual variation in natural populations. Individual variation in mating patterns may have significant adaptive implications for a population and for the entire species. Nevertheless, field data rarely address individual differences in mating patterns, focusing rather on averages. This study aimed to quantify individual variation in the different components of mating patterns. METHODS: Microsatellite data were used from 421 adult trees and 1911 seeds, structured in 72 half-sib families collected in a single mixed stand of Quercus robur and Q. petraea in northern Poland. Using a Bayesian approach, mating patterns were investigated, taking into account pollen dispersal, male fecundity, possible hybridization and heterogeneity in immigrant pollen pools. KEY RESULTS: Pollen dispersal followed a heavy-tailed distribution (283 m on average). In spite of high pollen mobility, immigrant pollen pools showed strong genetic structuring among mothers. At the individual level, immigrant pollen pools showed highly variable divergence rates, revealing that sources of immigrant pollen can vary greatly among particular trees. Within the stand, the distribution of male fecundity appeared highly skewed, with a small number of dominant males, resulting in a ratio of census to effective density of pollen donors of 5·3. Male fecundity was not correlated with tree diameter but showed strong cline-like spatial variation. This pattern can be attributed to environmental variation. Quercus petraea revealed a greater preference (74 %) towards intraspecific mating than Q. robur (36 %), although mating preferences varied among trees. CONCLUSIONS: Mating patterns can reveal great variation among individuals, even within a single even-age stand. The results show that trees can mate assortatively, with little respect for spatial proximity. Such selective mating may be a result of variable combining compatibility among trees due to genetic and/or environmental factors.

NM+: software implementing parentage-based models for estimating gene dispersal and mating patterns in plants.

NM+ is computer software designed for making inferences on plant gene dispersal and mating patterns by modelling parentage probabilities of offspring based on the spatially explicit neighbourhood model. NM+ requires a sample of mapped and genotyped candidate parents and offspring; however, offspring may optionally be assigned to single maternal parents (forming so-called half-sib progeny arrays). Using maximum likelihood, NM+ estimates a number of parameters, including proportions of offspring due to self-fertilization, pollen immigration from outside of a defined study site, parameters of pollen (and/or seed) dispersal kernels (exponential-power, Weibull, geometric or 2Dt) and selection gradients relating covariates (phenotypic traits) with male (and/or female) reproductive success. NM+ allows for missing data both in parents and in offspring. It accounts for null alleles and their frequencies can optionally be considered as estimable parameters. Data files are formatted in a table-like structure so they can be easily prepared in a spreadsheet software. By default NM+ is for studying plant populations, however, it can be used for any organism as long as data requirements and model assumptions are met. NM+ runs under Windows, but it can be launched under Linux using WINE emulator. NM+ can be downloaded free of charge from http://www.genetyka.ukw.edu.pl/index_pliki/software.htm.

Simultaneous estimation of null alleles and inbreeding coefficients.

Although microsatellites are a very efficient tool for many population genetics applications, they may occasionally produce "null" alleles, which, when present in high proportion, may affect estimates of key parameters such as inbreeding and relatedness coefficients or measures of genetic differentiation. In order to account for the presence of null alleles, it is first necessary to estimate their frequency within studied populations. However, the commonly used null allele frequency estimators are not of general applicability because they can produce upwardly biased estimates when a population under study experiences some inbreeding. In such a case, 2 formerly described approaches, population inbreeding model and individual inbreeding model, can be applied for simultaneous estimation of null allele frequencies and of the inbreeding coefficient. In this study, we demonstrate the properties and utility of these 2 methods and show that they outperform the commonly used approaches in the estimation of null allele frequencies based on genotypic data. The methods are applied to empirical data from a natural population of European beech (Fagus sylvatica L.), and results are briefly discussed. The methods presented in this paper are implemented in the Windows-based user-friendly INEST computer program (available free of charge at http://genetyka.ukw.edu.pl/INEst10_setup.exe).

Presence of triploids among oak species.

BACKGROUND AND AIMS: Polyploids, organisms that have multiple sets of chromosomes, are common in certain plant and animal taxa. However, there are only a few reports of intraspecific ploidy variation within the genus Quercus. The aim of the study was to investigate the suspected ploidy level of two oaks that have unusual microsatellite banding patterns. METHODS: Polyploidy was investigated by using microsatellite analysis, stomata length measurements and nuclear DNA content estimation by flow cytometry. KEY RESULTS: Each putative triploid tree has patterns of microsatellite variation unexpected for diploid genomes, with up to three alleles at some loci, significantly longer stomata and 1.5 times more DNA per nucleus compared with diploids. CONCLUSIONS: To our knowledge, this report contains the first evidence for triploidy in Q. petraea and confirmation of this phenomenon in Q. robur. Regardless of the positive or negative aspects of the presence of triploid oaks in forest stands, it is of value to be able to screen for them. This study demonstrates that nuclear microsatellites and estimation of DNA content by flow cytometry can readily be used for this purpose.

Effective gene dispersal and female reproductive success in Mediterranean maritime pine (Pinus pinaster Aiton).

Understanding population-scale processes that affect allele frequency changes across generations is a long-standing interest in genetic, ecological and evolutionary research. In particular, individual differences in female reproductive success and the spatial scale of gene flow considerably affect evolutionary change and patterns of local selection. In this study, a recently developed maximum-likelihood (ML) method based on established offspring, the Seedling Neighbourhood Model, was applied and exponentially shaped dispersal kernels were fitted to both genetic and ecological data in a widespread Mediterranean pine, Pinus pinaster Aiton. The distribution of female reproductive success in P. pinaster was very skewed (about 10% of trees mothered 50% of offspring) and significant positive female selection gradients for diameter (gamma = 0.7293) and cone crop (gamma = 0.4524) were found. The selective advantage of offspring mothered by bigger trees could be due to better-quality seeds. These seeds may show more resilience to severe summer droughts and microsite variation related to water and nutrient availability. Both approaches, ecological and of parentage, consistently showed a long-distance dispersal component in saplings that was not found in dispersal kernels based on seed shadows, highlighting the importance of Janzen-Connell effects and microenvironmental variation for survival at early stages of establishment in this Mediterranean key forest tree.

Historical and contemporary mating patterns in remnant populations of the forest tree Fraxinus excelsior L.

Genetic variation at microsatellite markers was used to quantify genetic structure and mating behavior in a severely fragmented population of the wind-pollinated, wind-dispersed temperate tree Fraxinus excelsior in a deforested catchment in Scotland. Remnants maintain high levels of genetic diversity, comparable with those reported for continuous populations in southeastern Europe, and show low interpopulation differentiation (E = 0.080), indicating that historical gene exchange has not been limited (Nm = 3.48). We estimated from seeds collected from all trees producing fruits in three of five remnants that F. excelsior is predominantly outcrossing (t(m) = 0.971 +/- 0.028). Use of a neighborhood model approach to describe the relative contribution of local and long-distance pollen dispersal indicates that pollen gene flow into each of the three remnants is extensive (46-95%) and pollen dispersal has two components. The first is very localized and restricted to tens of meters around the mother trees. The second is a long-distance component with dispersal occurring over several kilometers. Effective dispersal distances, accounting for the distance and directionality to mother trees of sampled pollen donors, average 328 m and are greater than values reported for a continuous population. These results suggest that the opening of the landscape facilitates airborne pollen movement and may alleviate the expected detrimental genetic effects of fragmentation.

Parentage versus two-generation analyses for estimating pollen-mediated gene flow in plant populations.

Assessment of contemporary pollen-mediated gene flow in plants is important for various aspects of plant population biology, genetic conservation and breeding. Here, through simulations we compare the two alternative approaches for measuring pollen-mediated gene flow: (i) the NEIGHBORHOOD model--a representative of parentage analyses, and (ii) the recently developed TWOGENER analysis of pollen pool structure. We investigate their properties in estimating the effective number of pollen parents (N(ep)) and the mean pollen dispersal distance (delta). We demonstrate that both methods provide very congruent estimates of N(ep) and delta, when the methods' assumptions considering the shape of pollen dispersal curve and the mating system follow those used in data simulations, although the NEIGHBORHOOD model exhibits generally lower variances of the estimates. The violations of the assumptions, especially increased selfing or long-distance pollen dispersal, affect the two methods to a different degree; however, they are still capable to provide comparable estimates of N(ep). The NEIGHBORHOOD model inherently allows to estimate both self-fertilization and outcrossing due to the long-distance pollen dispersal; however, the TWOGENER method is particularly sensitive to inflated selfing levels, which in turn may confound and suppress the effects of distant pollen movement. As a solution we demonstrate that in case of TWOGENER it is possible to extract the fraction of intraclass correlation that results from outcrossing only, which seems to be very relevant for measuring pollen-mediated gene flow. The two approaches differ in estimation precision and experimental efforts but they seem to be complementary depending on the main research focus and type of a population studied.

Cautions on direct gene flow estimation in plant populations.

Through simulations we have investigated the statistical properties of two of the main approaches for directly estimating pollen gene flow (m) in plant populations: genotypic exclusion and mating models. When the assumptions about accurately known background pollen pool allelic frequencies are met, both methods provide unbiased results with comparable variances across a range of true m values. However, when presumed allelic frequencies differ from actual ones, which is more likely in research practice, both estimators are biased. We demonstrate that the extent and direction of bias largely depend on the difference (measured as genetic distance) between the presumed and actual pollen pools, and on the degree of genetic differentiation between the local population and the actual background pollen sources. However, one feature of the mating model is its ability to estimate pollen gene flow simultaneously with background pollen pool allelic frequencies. We have found that this approach gives nearly unbiased pollen gene flow estimates, and is practical because it eliminates the necessity of providing independent estimates of background pollen pool allelic frequencies. Violations of the mating model assumptions of random mating within local population affect the precision of the estimates only to a limited degree.