Genomic basis of European ash tree resistance to ash dieback fungus

Transcriptional time-course analysis during ash dieback infection revealed different responses in tolerant and susceptible Fraxinus excelsior genotypes

Hymenoscyphus fraxineus, the causal agent of Ash Dieback (ADB), has been introduced to eastern Europe in the 1990s from where it spread causing decline in European ash populations. However, the genetic basis of the molecular response in tolerant and susceptible ash trees to this disease is still largely unknown. We performed RNA-sequencing to study the transcriptomic response to the disease in four ash genotypes (ADB-tolerant FAR3 and FS36, and ADB-susceptible UW1 and UW2), during a time-course of 7, 14, 21, and 28 days post-inoculation, including mock-inoculated trees as control samples for each sampling time point. The analysis yielded 395 and 500 Differentially Expressed Genes (DEGs) along the response for ADB-tolerant FAR3 and FS36, respectively, while ADB-susceptible UW1 and UW2 revealed 194 and 571 DEGs, respectively, with most DEGs found exclusively in just one of the genotypes. DEGs shared between tolerant genotypes FAR3 and FS36, included genes involved in the production of phytoalexins and other secondary metabolites with roles in plant defense. Moreover, we identified an earlier expression of genes involved in both pattern- and effector-triggered immunity (PTI and ETI) in ADB-tolerant genotypes, while in ADB-susceptible genotypes both responses were delayed (late response). Overall, these results revealed different transcriptomic expression patterns not only between ADB-tolerant and ADB-susceptible genotypes, but also within these two groups. This hints to individual responses in the natural tolerance to ADB, possibly revealing diversified strategies across ash genotypes.

Multiple, Single Trait GWAS and Supervised Machine Learning Reveal the Genetic Architecture of Fraxinus excelsior Tolerance to Ash Dieback in Europe

Common ash (Fraxinus excelsior) is under intensive attack from the invasive alien pathogenic fungus Hymenoscyphus fraxineus, causing ash dieback at epidemic levels throughout Europe. Previous studies have found significant genetic variation among genotypes in ash dieback susceptibility and that host phenology, such as autumn yellowing, is correlated with susceptibility of ash trees to H. fraxineus; however, the genomic basis of ash dieback tolerance in F. excelsior requires further investigation. Here, we integrate quantitative genetics based on multiple replicates and genome-wide association analyses with machine learning to reveal the genetic architecture of ash dieback tolerance and of phenological traits in F. excelsior populations in six European countries (Austria, Denmark, Germany, Ireland, Lithuania, Sweden). Based on phenotypic data of 486 F. excelsior replicated genotypes we observed negative genotypic correlations between crown damage caused by ash dieback and intensity of autumn leaf yellowing within multiple sampling sites. Our results suggest that the examined traits are polygenic and using genomic prediction models, with ranked single nucleotide polymorphisms (SNPs) based on GWAS associations as input, a large proportion of the variation was predicted by unlinked SNPs. Based on 100 unlinked SNPs, we can predict 55% of the variation in disease tolerance among genotypes (as phenotyped in genetic trials), increasing to a maximum of 63% when predicted from 9155 SNPs. In autumn leaf yellowing, 52% of variation is predicted by 100 unlinked SNPs, reaching a peak of 72% using 3740 SNPs. Based on feature permutations within genomic prediction models, a total of eight nonsynonymous SNPs linked to ash dieback crown damage and autumn leaf yellowing (three and five SNPs, respectively) were identified, these were located within genes related to plant defence (pattern triggered immunity, pathogen detection) and phenology (regulation of flowering and seed maturation, auxin transport). We did not find an overlap between genes associated with crown damage level and autumn leaf yellowing. Hence, our results shed light on the difference in the genomic basis of ADB tolerance and autumn leaf yellowing despite these two traits being correlated in quantitative genetic analysis. Overall, our methods show the applicability of genomic prediction models when combined with GWAS to reveal the genomic architecture of polygenic disease tolerance enabling the identification of ash dieback tolerant trees for breeding or conservation purposes.

β-Aminobutyric acid promotes stress tolerance, physiological adjustments, as well as broad epigenetic changes at DNA and RNA nucleobases in field elms (Ulmus minor)

BACKGROUND

β-Aminobutyric acid (BABA) has been successfully used to prime stress resistance in numerous plant species; however, its effectiveness in forest trees has been poorly explored thus far. This study aimed to investigate the influence of BABA on morphological, physiological, and epigenetic parameters in field elms under various growth conditions. Epigenetic changes were assessed in both DNA and RNA through the use of reversed-phase ultra-performance liquid chromatography (UPLC) coupled with sensitive mass spectrometry.

RESULTS

The presented results confirm the influence of BABA on the development, physiology, and stress tolerance in field elms. However, the most important findings are related to the broad epigenetic changes promoted by this amino acid, which involve both DNA and RNA. Our findings confirm, for the first time, that BABA influences not only well-known epigenetic markers in plants, such as 5-methylcytosine, but also several other non-canonical nucleobases, such as 5-hydroxymethyluracil, 5-formylcytosine, 5-hydroxymethylcytosine, N6-methyladenine, uracil (in DNA) and thymine (in RNA). The significant effect on the levels of N6-methyladenine, the main bacterial epigenetic marker, is particularly noteworthy. In this case, the question arises as to whether this effect is due to epigenetic changes in the microbiome, the plant genome, or both.

CONCLUSIONS

The plant phenotype is the result of complex interactions between the plant's DNA, the microbiome, and the environment. We propose that different types of epigenetic changes in the plant and microbiome may play important roles in the largely unknown memory process that enables plants to adapt faster to changing environmental conditions.

Conifers Concentrate Large Numbers of NLR Immune Receptor Genes on One Chromosome

Nucleotide-binding domain and leucine-rich repeat (NLR) immune receptor genes form a major line of defense in plants, acting in both pathogen recognition and resistance machinery activation. NLRs are reported to form large gene clusters in limber pine (Pinus flexilis), but it is unknown how widespread this genomic architecture may be among the extant species of conifers (Pinophyta). We used comparative genomic analyses to assess patterns in the abundance, diversity, and genomic distribution of NLR genes. Chromosome-level whole genome assemblies and high-density linkage maps in the Pinaceae, Cupressaceae, Taxaceae, and other gymnosperms were scanned for NLR genes using existing and customized pipelines. The discovered genes were mapped across chromosomes and linkage groups and analyzed phylogenetically for evolutionary history. Conifer genomes are characterized by dense clusters of NLR genes, highly localized on one chromosome. These clusters are rich in TNL-encoding genes, which seem to have formed through multiple tandem duplication events. In contrast to angiosperms and nonconiferous gymnosperms, genomic clustering of NLR genes is ubiquitous in conifers. NLR-dense genomic regions are likely to influence a large part of the plant's resistance, informing our understanding of adaptation to biotic stress and the development of genetic resources through breeding.

A glimmer of hope - ash genotypes with increased resistance to ash dieback pathogen show cross-resistance to emerald ash borer

Plants rely on cross-resistance traits to defend against multiple, phylogenetically distinct enemies. These traits are often the result of long co-evolutionary histories. Biological invasions can force naïve plants to cope with novel, coincident pests, and pathogens. For example, European ash (Fraxinus excelsior) is substantially threatened by the emerald ash borer (EAB), Agrilus planipennis, a wood-boring beetle, and the ash dieback (ADB) pathogen, Hymenoscyphus fraxineus. Yet, plant cross-resistance traits against novel enemies are poorly explored and it is unknown whether naïve ash trees can defend against novel enemy complexes via cross-resistance mechanisms. To gain mechanistic insights, we quantified EAB performance on grafted replicates of ash genotypes varying in ADB resistance and characterized ash phloem chemistry with targeted and untargeted metabolomics. Emerald ash borer performed better on ADB-susceptible than on ADB-resistant genotypes. Moreover, changes in EAB performance aligned with differences in phloem chemical profiles between ADB-susceptible and ADB-resistant genotypes. We show that intraspecific variation in phloem chemistry in European ash can confer increased cross-resistance to invasive antagonists from different taxonomic kingdoms. Our study suggests that promotion of ADB-resistant ash genotypes may simultaneously help to control the ADB disease and reduce EAB-caused ash losses, which may be critical for the long-term stability of this keystone tree species.

A first draft genome of holm oak (Quercus ilex subsp. ballota), the most representative species of the Mediterranean forest and the Spanish agrosylvopastoral ecosystem "dehesa"

The holm oak (Quercus ilex subsp. ballota) is the most representative species of the Mediterranean Basin and the agrosylvopastoral Spanish "dehesa" ecosystem. Being part of our life, culture, and subsistence since ancient times, it has significant environmental and economic importance. More recently, there has been a renewed interest in using the Q. ilex acorn as a functional food due to its nutritional and nutraceutical properties. However, the holm oak and its related ecosystems are threatened by different factors, with oak decline syndrome and climate change being the most worrying in the short and medium term. Breeding programs informed by the selection of elite genotypes seem to be the most plausible biotechnological solution to rescue populations under threat. To achieve this and other downstream analyses, we need a high-quality and well-annotated Q. ilex reference genome. Here, we introduce the first draft genome assembly of Q. ilex using long-read sequencing (PacBio). The assembled nuclear haploid genome had 530 contigs totaling 842.2 Mbp (N50 = 3.3 Mbp), of which 448.7 Mb (53%) were repetitive sequences. We annotated 39,443 protein-coding genes of which 94.80% were complete and single-copy genes. Phylogenetic analyses showed no evidence of a recent whole-genome duplication, and high synteny of the 12 chromosomes between Q. ilex and Quercus lobata and between Q. ilex and Quercus robur. The chloroplast genome size was 142.3 Kbp with 149 protein-coding genes successfully annotated. This first draft should allow for the validation of omics data as well as the identification and functional annotation of genes related to phenotypes of interest such as those associated with resilience against oak decline syndrome and climate change and higher acorn productivity and nutraceutical value.

Exploring the crop epigenome: a comparison of DNA methylation profiling techniques

Epigenetic modifications play a vital role in the preservation of genome integrity and in the regulation of gene expression. DNA methylation, one of the key mechanisms of epigenetic control, impacts growth, development, stress response and adaptability of all organisms, including plants. The detection of DNA methylation marks is crucial for understanding the mechanisms underlying these processes and for developing strategies to improve productivity and stress resistance of crop plants. There are different methods for detecting plant DNA methylation, such as bisulfite sequencing, methylation-sensitive amplified polymorphism, genome-wide DNA methylation analysis, methylated DNA immunoprecipitation sequencing, reduced representation bisulfite sequencing, MS and immuno-based techniques. These profiling approaches vary in many aspects, including DNA input, resolution, genomic region coverage, and bioinformatics analysis. Selecting an appropriate methylation screening approach requires an understanding of all these techniques. This review provides an overview of DNA methylation profiling methods in crop plants, along with comparisons of the efficacy of these techniques between model and crop plants. The strengths and limitations of each methodological approach are outlined, and the importance of considering both technical and biological factors are highlighted. Additionally, methods for modulating DNA methylation in model and crop species are presented. Overall, this review will assist scientists in making informed decisions when selecting an appropriate DNA methylation profiling method.

Applying molecular and genetic methods to trees and their fungal communities

Forests provide invaluable economic, ecological, and social services. At the same time, they are exposed to several threats, such as fragmentation, changing climatic conditions, or increasingly destructive pests and pathogens. Trees, the inherent species of forests, cannot be viewed as isolated organisms. Manifold (micro)organisms are associated with trees playing a pivotal role in forest ecosystems. Of these organisms, fungi may have the greatest impact on the life of trees. A multitude of molecular and genetic methods are now available to investigate tree species and their associated organisms. Due to their smaller genome sizes compared to tree species, whole genomes of different fungi are routinely compared. Such studies have only recently started in forest tree species. Here, we summarize the application of molecular and genetic methods in forest conservation genetics, tree breeding, and association genetics as well as for the investigation of fungal communities and their interrelated ecological functions. These techniques provide valuable insights into the molecular basis of adaptive traits, the impacts of forest management, and changing environmental conditions on tree species and fungal communities and can enhance tree-breeding cycles due to reduced time for field testing. It becomes clear that there are multifaceted interactions among microbial species as well as between these organisms and trees. We demonstrate the versatility of the different approaches based on case studies on trees and fungi. KEY POINTS: • Current knowledge of genetic methods applied to forest trees and associated fungi. • Genomic methods are essential in conservation, breeding, management, and research. • Important role of phytobiomes for trees and their ecosystems.

Genomic prediction in the wild: A case study in Soay sheep

Genomic prediction, the technique whereby an individual's genetic component of their phenotype is estimated from its genome, has revolutionised animal and plant breeding and medical genetics. However, despite being first introduced nearly two decades ago, it has hardly been adopted by the evolutionary genetics community studying wild organisms. Here, genomic prediction is performed on eight traits in a wild population of Soay sheep. The population has been the focus of a >30 year evolutionary ecology study and there is already considerable understanding of the genetic architecture of the focal Mendelian and quantitative traits. We show that the accuracy of genomic prediction is high for all traits, but especially those with loci of large effect segregating. Five different methods are compared, and the two methods that can accommodate zero-effect and large-effect loci in the same model tend to perform best. If the accuracy of genomic prediction is similar in other wild populations, then there is a real opportunity for pedigree-free molecular quantitative genetics research to be enabled in many more wild populations; currently the literature is dominated by studies that have required decades of field data collection to generate sufficiently deep pedigrees. Finally, some of the potential applications of genomic prediction in wild populations are discussed.

The Unusual Value of Long-Term Studies of Individuals: The Example of the Isle of Rum Red Deer Project

Long-term studies of individuals enable incisive investigations of questions across ecology and evolution. Here, we illustrate this claim by reference to our long-term study of red deer on the Isle of Rum, Scotland. This project has established many of the characteristics of social organization, selection, and population ecology typical of large, polygynous, seasonally breeding mammals, with wider implications for our understanding of sexual selection and the evolution of sex differences, as well as for their population dynamics and population management. As molecular genetic techniques have developed, the project has pivoted to investigate evolutionary genetic questions, also breaking new ground in this field. With ongoing advances in genomics and statistical approaches and the development of increasingly sophisticated ways to assay new phenotypic traits, the questions that long-term studies such as the red deer study can answer become both broader and ever more sophisticated. They also offer powerful means of understanding the effects of ongoing climate change on wild populations.

Genome-Wide SNP Markers Accelerate Perennial Forest Tree Breeding Rate for Disease Resistance through Marker-Assisted and Genome-Wide Selection

The ecological and economic importance of forest trees is evident and their survival is necessary to provide the raw materials needed for wood and paper industries, to preserve the diversity of associated animal and plant species, to protect water and soil, and to regulate climate. Forest trees are threatened by anthropogenic factors and biotic and abiotic stresses. Various diseases, including those caused by fungal pathogens, are one of the main threats to forest trees that lead to their dieback. Genomics and transcriptomics studies using next-generation sequencing (NGS) methods can help reveal the architecture of resistance to various diseases and exploit natural genetic diversity to select elite genotypes with high resistance to diseases. In the last two decades, QTL mapping studies led to the identification of QTLs related to disease resistance traits and gene families and transcription factors involved in them, including NB-LRR, WRKY, bZIP and MYB. On the other hand, due to the limitation of recombination events in traditional QTL mapping in families derived from bi-parental crosses, genome-wide association studies (GWAS) that are based on linkage disequilibrium (LD) in unstructured populations overcame these limitations and were able to narrow down QTLs to single genes through genotyping of many individuals using high-throughput markers. Association and QTL mapping studies, by identifying markers closely linked to the target trait, are the prerequisite for marker-assisted selection (MAS) and reduce the breeding period in perennial forest trees. The genomic selection (GS) method uses the information on all markers across the whole genome, regardless of their significance for development of a predictive model for the performance of individuals in relation to a specific trait. GS studies also increase gain per unit of time and dramatically increase the speed of breeding programs. This review article is focused on the progress achieved in the field of dissecting forest tree disease resistance architecture through GWAS and QTL mapping studies. Finally, the merit of methods such as GS in accelerating forest tree breeding programs is also discussed.

Pathogenic strategies and immune mechanisms to necrotrophs: Differences and similarities to biotrophs and hemibiotrophs

Pathogenesis in plant diseases is complex comprising diverse pathogen virulence and plant immune mechanisms. These pathogens cause damaging plant diseases by deploying specialized and generic virulence strategies that are countered by intricate resistance mechanisms. The significant challenges that necrotrophs pose to crop production are predicted to increase with climate change. Immunity to biotrophs and hemibiotrophs is dominated by intracellular receptors that recognize specific effectors and activate resistance. These mechanisms play only minor roles in resistance to necrotrophs. Pathogen- or host-derived conserved pattern molecules trigger immune responses that broadly contribute to plant immunity. However, certain pathogen or host-derived immune elicitors are enriched by the virulence activities of necrotrophs. Different plant hormones modulate systemic resistance and cell death that have differential impacts on resistance to pathogens of different lifestyles. Knowledge of mechanisms that contribute to resistance to necrotrophs has expanded. Besides toxins and cell wall degrading enzymes that dominate the pathogenesis of necrotrophs, other effectors with subtle contributions are being identified.

Multiomics Molecular Research into the Recalcitrant and Orphan Quercus ilex Tree Species: Why, What for, and How

The holm oak (Quercus ilex L.) is the dominant tree species of the Mediterranean forest and the Spanish agrosilvopastoral ecosystem, "dehesa." It has been, since the prehistoric period, an important part of the Iberian population from a social, cultural, and religious point of view, providing an ample variety of goods and services, and forming the basis of the economy in rural areas. Currently, there is renewed interest in its use for dietary diversification and sustainable food production. It is part of cultural richness, both economically (tangible) and environmentally (intangible), and must be preserved for future generations. However, a worrisome degradation of the species and associated ecosystems is occurring, observed in an increase in tree decline and mortality, which requires urgent action. Breeding programs based on the selection of elite genotypes by molecular markers is the only plausible biotechnological approach. To this end, the authors' group started, in 2004, a research line aimed at characterizing the molecular biology of Q. ilex. It has been a challenging task due to its biological characteristics (long life cycle, allogamous, high phenotypic variability) and recalcitrant nature. The biology of this species has been characterized following the central dogma of molecular biology using the omics cascade. Molecular responses to biotic and abiotic stresses, as well as seed maturation and germination, are the two main objectives of our research. The contributions of the group to the knowledge of the species at the level of DNA-based markers, genomics, epigenomics, transcriptomics, proteomics, and metabolomics are discussed here. Moreover, data are compared with those reported for Quercus spp. All omics data generated, and the genome of Q. ilex available, will be integrated with morphological and physiological data in the systems biology direction. Thus, we will propose possible molecular markers related to resilient and productive genotypes to be used in reforestation programs. In addition, possible markers related to the nutritional value of acorn and derivate products, as well as bioactive compounds (peptides and phenolics) and allergens, will be suggested. Subsequently, the selected molecular markers will be validated by both genome-wide association and functional genomic analyses.

Molecular and Chemical Screening for Inherent Disease Resistance Factors of Norway Spruce (Picea abies) Clones Against Conifer Stem Rot Pathogen Heterobasidion parviporum

Root and stem rot of conifer trees caused by Heterobasidion annosum species complex leads to huge economic losses in Europe, yet not much is known about the molecular and chemical basis for host resistance. To identify inherent chemical or molecular markers in clones found to be either resistant or susceptible, we sampled needle tissues of all the clones before pathogen inoculation. We conducted a short-term resistance screening by using the pathogen H. parviporum to inoculate 70 Norway spruce clones. Based on lesion size, subsets of highly susceptible and resistant clones were further analyzed. Terpene detection and RNA sequencing were performed to explore inherent variations in genotypes differing in resistance to pathogenic challenge at chemical and transcriptional levels. A negative correlation emerged between resistance and growth. Terpene profiles of resistant clones showed higher content of monoterpenes and sesquiterpenes, with concomitant increased transcript abundance of genes involved in the terpenoid pathway. A set of upregulated genes relevant to flavonoid biosynthesis was observed in resistant genotypes, whereas higher transcripts of lignin biosynthetic genes were prevalent in susceptible clones. Genes involved in flavonoid and lignin biosynthesis as well as terpene content may have a role in facilitating resistance of Norway spruce against H. parviporum. Our results provide strong support on the feasibility of sampling needle tissues before pathogen inoculation, and the approach could be of value for large-scale screening of novel biomarkers for durable resistance. The additional insights could form a basis for further research on resistance screening in this pathosystem.

European-wide forest monitoring substantiate the neccessity for a joint conservation strategy to rescue European ash species (Fraxinus spp.)

European ash (Fraxinus excelsior) and narrow-leafed ash (F. angustifolia) are keystone forest tree species with a broad ecological amplitude and significant economic importance. Besides global warming both species are currently under significant threat by an invasive fungal pathogen that has been spreading progressively throughout the continent for almost three decades. Ash dieback caused by the ascomycete Hymenoscyphus fraxineus is capable of damaging ash trees of all age classes and often ultimately leads to the death of a tree after years of progressively developing crown defoliation. While studies at national and regional level already suggested rapid decline of ash populations as a result of ash dieback, a comprehensive survey at European level with harmonized crown assessment data across countries could shed more light into the population decline from a pan-European perspective and could also pave the way for a new conservation strategy beyond national boarders. Here we present data from the ICP Forests Level I crown condition monitoring from 27 countries resulting in > 36,000 observations. We found a substantial increase in defoliation and mortality over time indicating that crown defoliation has almost doubled during the last three decades. Hotspots of mortality are currently situated in southern Scandinavia and north-eastern Europe. Overall survival probability after nearly 30 years of infection has already reached a critical value of 0.51, but with large differences among regions (0.20–0.86). Both a Cox proportional hazard model as well as an Aalen additive regression model strongly suggest that survival of ash is significantly lower in locations with excessive water regime and which experienced more extreme precipitation events during the last two decades. Our results underpin the necessity for fast governmental action and joint rescue efforts beyond national borders since overall mean defoliation will likely reach 50% as early as 2030 as suggested by time series forecasting.

The Impact of the Neophyte Tree Fraxinus pennsylvanica [Marshall] on Beetle Diversity under Climate Change

We studied the impact of the neophyte tree Fraxinus pennsylvanica on the diversity of beetles in floodplain forests along the river Elbe in Germany in 2016, 2017 and in 2020, where 80% of all Fraxinus excelsior trees had died following severe droughts. Beetles were collected by insecticidal knock-down from 121 trees (64 F. excelsior and 57 F. pennsylvanica) and identified to 547 species in 15,214 specimens. The trees sampled in 2016 and 2017 showed no signs of drought stress or ash dieback and serve as a reference for the comparison with the 2020 fauna. The data proved that F. excelsior harbours the most diverse beetle community, which differed also significantly in guild composition from F. pennsylvanica. Triggered by extremely dry and long summer seasons, the 2020 ash dieback had profound and forest-wide impacts. Several endangered, red-listed beetle species of Saxonia Anhalt had increased in numbers and became secondary pests on F. excelsior. Diversity decreased whilst numbers of xylobionts increased on all trees, reaching 78% on F. excelsior. Proportions of xylobionts remained constant on F. pennsylvanica. Phytophages were almost absent from all trees, but mycetophages increased on F. pennsylvanica. Our data suggest that as a result of the dieback of F. excelsior the neophyte F. pennsylvanica might become a rescue species for the European Ash fauna, as it provides the second-best habitat. We show how difficult it is to assess the dynamics and the ecological impact of neophytes, especially under conditions similar to those projected by climate change models. The diversity and abundance of canopy arthropods demonstrates their importance in understanding forest functions and maintenance of ecosystem services, illustrating that their consideration is essential for forest adaptation to climate change.

Host−Pathogen Interactions in Leaf Petioles of Common Ash and Manchurian Ash Infected with Hymenoscyphus fraxineus

Some common ash trees (Fraxinus excelsior) show tolerance towards shoot dieback caused by the invasive ascomycete Hymenoscyphus fraxineus. Leaf petioles are considered to serve as a pathogen colonization route to the shoots. We compared four common ash clones with variation in disease tolerance, and included the native host, Manchurian ash (Fraxinus mandshurica), as a reference. Tissue colonization, following rachis inoculation by H. fraxineus, was monitored by histochemical observations and a quantitative polymerase chain reaction (qPCR) assay specific to H. fraxineus. Axial spread of the pathogen towards the petiole base occurred primarily within the phloem and parenchyma, tissues rich in starch in healthy petioles. In inoculated petioles, a high content of phenolics surrounded the hyphae, presumably a host defense response. There was a relationship between field performance and susceptibility to leaf infection in three of the four studied common ash clones, i.e., good field performance was associated with a low petiole colonization level and vice versa. Low susceptibility to leaf infection may counteract leaf-to-shoot spread of the pathogen in common ash, but the limited number of clones studied warrants caution and a larger study. The Manchurian ash clone had the highest petiole colonization level, which may suggest that this native host has evolved additional mechanisms to avoid shoot infection.

Identifying and testing marker-trait associations for growth and phenology in three pine species: Implications for genomic prediction

In tree species, genomic prediction offers the potential to forecast mature trait values in early growth stages, if robust marker-trait associations can be identified. Here we apply a novel multispecies approach using genotypes from a new genotyping array, based on 20,795 single nucleotide polymorphisms (SNPs) from three closely related pine species (Pinus sylvestris, Pinus uncinata and Pinus mugo), to test for associations with growth and phenology data from a common garden study. Predictive models constructed using significantly associated SNPs were then tested and applied to an independent multisite field trial of P. sylvestris and the capability to predict trait values was evaluated. One hundred and eighteen SNPs showed significant associations with the traits in the pine species. Common SNPs (MAF > 0.05) associated with bud set were only found in genes putatively involved in growth and development, whereas those associated with growth and budburst were also located in genes putatively involved in response to environment and, to a lesser extent, reproduction. At one of the two independent sites, the model we developed produced highly significant correlations between predicted values and observed height data (YA, height 2020: r = 0.376, p < 0.001). Predicted values estimated with our budburst model were weakly but positively correlated with duration of budburst at one of the sites (GS, 2015: r = 0.204, p = 0.034; 2018: r = 0.205, p = 0.034-0.037) and negatively associated with budburst timing at the other (YA: r = -0.202, p = 0.046). Genomic prediction resulted in the selection of sets of trees whose mean height was taller than the average for each site. Our results provide tentative support for the capability of prediction models to forecast trait values in trees, while highlighting the need for caution in applying them to trees grown in different environments.

The Impact of Biotic and Abiotic Stress Factors on Development of European Ash Tissue Cultures

Fraxinus excelsior L. is threatened by a variety of environmental factors causing a decline of the species. The most important biotic factors negatively affecting the condition of the F. excelsior population are fungi such as the pathogen Hymenoscyphus fraxineus. Abiotic factors with potentially harmful effect to the F. excelsior population are the accumulation of heavy metals and salinity in soils. Thus, the aim of this study was to investigate the impact of selected biotic and abiotic stress factors to determine which of them pose a threat to European ash. The study was conducted using in vitro techniques based on callus and seedlings regenerated via indirect organogenesis. Tissue cultures exclude the influence of other factors, including the environmental impact on ash extinction. The results confirmed very strong pathogenic potential of H. fraxineus in which after 14 days the callus tissue cells died as the tissue failed to activate its defense mechanisms. Experiments showed the high toxicity of cadmium in concentration of 0.027 mmol/L. Salinity caused the activity of oxidation enzymes to vary among seedlings and calluses in the control suggesting the enzymes play a role in controlling the morphogenetic development of tissue cultures.

Transformation of European Ash (Fraxinus excelsior L.) Callus as a Starting Point for Understanding the Molecular Basis of Ash Dieback

The population of European ash (Fraxinus excelsior L.) is currently facing the risk of collapse, mainly due to ash dieback, a disease caused by a pathogenic fungus, Hymenoscyphus fraxineus. To facilitate studies into the molecular basis of ash dieback and design breeding strategies for a generation of resistant trees, it is necessary to develop tools enabling the study of gene function in F. excelsior. Despite this, a method for the genetic engineering of F. excelsior is still missing. Here, we report the first successful genetic transformation of F. excelsior callus and a selection process enabling the formation of stable transgenic callus lines. The protocol relies on the use of Agrobacterium tumefaciens to transform callus tissue derived from embryos of F. excelsior. In our experiments, we used the β-glucuronidase (GUS) reporter system to demonstrate the transformation of callus cells and performed RT-PCR experiments to confirm the stable expression of the transgene. Since ash dieback threatens the long-term stability of many native F. excelsior populations, we hope that the transformation techniques described in this manuscript will facilitate rapid progress in uncovering the molecular basis of the disease and the validation of gene targets previously proposed to be linked to the resistance of trees to H. fraxineus pathogenicity.

SUMOylation in Phytopathogen Interactions: Balancing Invasion and Resistance

Plants are constantly confronted by a multitude of biotic stresses involving a myriad of pathogens. In crops, pathogen infections result in significant agronomical losses worldwide posing a threat to food security. In order to enter plant tissues and establish a successful infection, phytopathogens have to surpass several physical, and chemical defense barriers. In recent years, post-translational modification (PTM) mechanisms have emerged as key players in plant defense against pathogens. PTMs allow a highly dynamic and rapid response in front of external challenges, increasing the complexity and precision of cellular responses. In this review, we focus on the role of SUMO conjugation (SUMOylation) in plant immunity against fungi, bacteria, and viruses. In plants, SUMO regulates multiple biological processes, ranging from development to responses arising from environmental challenges. During pathogen attack, SUMO not only modulates the activity of plant defense components, but also serves as a target of pathogen effectors, highlighting its broad role in plant immunity. Here, we summarize known pathogenic strategies targeting plant SUMOylation and, the plant SUMO conjugates involved in host-pathogen interactions. We also provide a catalog of candidate SUMO conjugates according to their role in defense responses. Finally, we discuss the complex role of SUMO in plant defense, focusing on key biological and experimental aspects that contribute to some controversial conclusions, and the opportunities for improving agricultural productivity by engineering SUMOylation in crop species.

The utility of genomic prediction models in evolutionary genetics

Variation in complex traits is the result of contributions from many loci of small effect. Based on this principle, genomic prediction methods are used to make predictions of breeding value for an individual using genome-wide molecular markers. In breeding, genomic prediction models have been used in plant and animal breeding for almost two decades to increase rates of genetic improvement and reduce the length of artificial selection experiments. However, evolutionary genomics studies have been slow to incorporate this technique to select individuals for breeding in a conservation context or to learn more about the genetic architecture of traits, the genetic value of missing individuals or microevolution of breeding values. Here, we outline the utility of genomic prediction and provide an overview of the methodology. We highlight opportunities to apply genomic prediction in evolutionary genetics of wild populations and the best practices when using these methods on field-collected phenotypes.

Genomic basis for drought resistance in European beech forests threatened by climate change

In the course of global climate change, Central Europe is experiencing more frequent and prolonged periods of drought. The drought years 2018 and 2019 affected European beeches (Fagus sylvatica L.) differently: even in the same stand, drought-damaged trees neighboured healthy trees, suggesting that the genotype rather than the environment was responsible for this conspicuous pattern. We used this natural experiment to study the genomic basis of drought resistance with Pool-GWAS. Contrasting the extreme phenotypes identified 106 significantly associated single-nucleotide polymorphisms (SNPs) throughout the genome. Most annotated genes with associated SNPs (>70%) were previously implicated in the drought reaction of plants. Non-synonymous substitutions led either to a functional amino acid exchange or premature termination. An SNP assay with 70 loci allowed predicting drought phenotype in 98.6% of a validation sample of 92 trees. Drought resistance in European beech is a moderately polygenic trait that should respond well to natural selection, selective management, and breeding.

Genomic Selection for Forest Tree Improvement: Methods, Achievements and Perspectives

The breeding of forest trees is only a few decades old, and is a much more complicated, longer, and expensive endeavor than the breeding of agricultural crops. One breeding cycle for forest trees can take 20–30 years. Recent advances in genomics and molecular biology have revolutionized traditional plant breeding based on visual phenotype assessment: the development of different types of molecular markers has made genotype selection possible. Marker-assisted breeding can significantly accelerate the breeding process, but this method has not been shown to be effective for selection of complex traits on forest trees. This new method of genomic selection is based on the analysis of all effects of quantitative trait loci (QTLs) using a large number of molecular markers distributed throughout the genome, which makes it possible to assess the genomic estimated breeding value (GEBV) of an individual. This approach is expected to be much more efficient for forest tree improvement than traditional breeding. Here, we review the current state of the art in the application of genomic selection in forest tree breeding and discuss different methods of genotyping and phenotyping. We also compare the accuracies of genomic prediction models and highlight the importance of a prior cost-benefit analysis before implementing genomic selection. Perspectives for the further development of this approach in forest breeding are also discussed: expanding the range of species and the list of valuable traits, the application of high-throughput phenotyping methods, and the possibility of using epigenetic variance to improve of forest trees.

Analyzing Ash Leaf-Colonizing Fungal Communities for Their Biological Control of Hymenoscyphus fraxineus

The invasive ascomycete Hymenoscyphus fraxineus has been threatening Fraxinus excelsior populations throughout Europe for over two decades. Since the infection and first colonization by the pathogen occurs in leaves, leaf-colonizing microorganisms have been discussed as a barrier and as possible biocontrol agents against the disease. To identify fungal groups with health-supporting potential, we compared the fungal microbiota of compound leaves from susceptible and tolerant ash trees in four ash stands with high H. fraxineus exposure. The fungal communities were analyzed both culture-independently by ITS2 amplicon sequencing and by the taxonomic classification of 1,704 isolates using matrix-assisted laser desorption/ionization time of flight mass spectrometry (MALDI-TOF MS) or sequencing of the entire ITS region. The fungal community structure did not show significant differences depending on the health status. However, for several OTUs and a MALDI group, a significantly higher abundance was found in tolerant ash trees. Thus, the yeast Papiliotrema flavescens was significantly increased and accounted for 12.3% of the mycobiome of tolerant ashes (OTU0003), and it had also a distinctly higher abundance among the isolates. The filamentous ascomycete Sarocladium strictum was increased 24-fold among the isolates of tolerant trees, but its abundance was comparably low. An in vitro screening for the growth inhibition of the pathogen via cocultivation resulted in 28 yeast-like isolates and 79 filamentous fungi with antagonistic activity. A statistical cocultivation test on two H. fraxineus strains confirmed six of the yeast-like isolates that suppressed H. fraxineus significantly, from 39-50%, two of them through a fungicidal effect. The highest inhibition rates among the yeasts were found for three isolates belonging to Aureobasidium pullulans and P. flavescens. The cocultivation test of the filamentous isolates revealed higher effects compared to the yeasts. Four isolates showed significant inhibition of both H. fraxineus strains with a rate of 72-100%, and five further isolates inhibited only one H. fraxineus strain significantly. The most effective isolates were members of the genus Cladosporium. During the next step, in planta tests will be necessary to verify the efficacy of the antagonistic isolates and to assess their suitability as biocontrol agents.

Modern Strategies to Assess and Breed Forest Tree Adaptation to Changing Climate

Studying the genetics of adaptation to new environments in ecologically and industrially important tree species is currently a major research line in the fields of plant science and genetic improvement for tolerance to abiotic stress. Specifically, exploring the genomic basis of local adaptation is imperative for assessing the conditions under which trees will successfully adapt in situ to global climate change. However, this knowledge has scarcely been used in conservation and forest tree improvement because woody perennials face major research limitations such as their outcrossing reproductive systems, long juvenile phase, and huge genome sizes. Therefore, in this review we discuss predictive genomic approaches that promise increasing adaptive selection accuracy and shortening generation intervals. They may also assist the detection of novel allelic variants from tree germplasm, and disclose the genomic potential of adaptation to different environments. For instance, natural populations of tree species invite using tools from the population genomics field to study the signatures of local adaptation. Conventional genetic markers and whole genome sequencing both help identifying genes and markers that diverge between local populations more than expected under neutrality, and that exhibit unique signatures of diversity indicative of "selective sweeps." Ultimately, these efforts inform the conservation and breeding status capable of pivoting forest health, ecosystem services, and sustainable production. Key long-term perspectives include understanding how trees' phylogeographic history may affect the adaptive relevant genetic variation available for adaptation to environmental change. Encouraging "big data" approaches (machine learning-ML) capable of comprehensively merging heterogeneous genomic and ecological datasets is becoming imperative, too.

Canditate metabolites for ash dieback tolerance in Fraxinus excelsior

Ash dieback, a forest epidemic caused by the invasive fungus Hymenoscyphus fraxineus, threatens ash trees throughout Europe. Within Fraxinus excelsior populations, a small proportion of genotypes show a low susceptibility to the pathogen. We compared the metabolomes from a cohort of low-susceptibility ash genotypes with a cohort of high-susceptibility ash genotypes. This revealed two significantly different chemotypes. A total of 64 candidate metabolites associated with reduced or increased susceptibility in the chemical families secoiridoids, coumarins, flavonoids, phenylethanoids, and lignans. Increased levels of two coumarins, fraxetin and esculetin, were strongly associated with reduced susceptibility to ash dieback. Both coumarins inhibited the growth of H. fraxineus in vitro when supplied at physiological concentrations, thereby validating their role as markers for low susceptibility to ash dieback. Similarly, fungal growth inhibition was observed when the methanolic bark extract of low-susceptibility ash genotypes was supplied. Our findings indicate the presence of constitutive chemical defense barriers against ash dieback in ash.

Transcriptional responses in developing lesions of European common ash (Fraxinus excelsior) reveal genes responding to infection by Hymenoscyphus fraxineus

BACKGROUND

With the expanding ash dieback epidemic that has spread across the European continent, an improved functional understanding of the disease development in afflicted hosts is needed. The study investigated whether differences in necrosis extension between common ash (Fraxinus excelsior) trees with different levels of susceptibility to the fungus Hymenoscyphus fraxineus are associated with, and can be explained by, the differences in gene expression patterns. We inoculated seemingly healthy branches of each of two resistant and susceptible ash genotypes with H. fraxineus grown in a common garden.

RESULTS

Ten months after the inoculation, the length of necrosis on the resistant genotypes were shorter than on the susceptible genotypes. RNA sequencing of bark samples collected at the border of necrotic lesions and from healthy tissues distal to the lesion revealed relatively limited differences in gene expression patterns between susceptible and resistant genotypes. At the necrosis front, only 138 transcripts were differentially expressed between the genotype categories while 1082 were differentially expressed in distal, non-symptomatic tissues. Among these differentially expressed genes, several genes in the mevalonate (MVA) and iridoid pathways were found to be co-regulated, possibly indicating increased fluxes through these pathways in response to H. fraxineus. Comparison of transcriptional responses of symptomatic and non-symptomatic ash in a controlled greenhouse experiment revealed a relatively small set of genes that were differentially and concordantly expressed in both studies. This gene-set included the rate-limiting enzyme in the MVA pathway and a number of transcription factors. Furthermore, several of the concordantly expressed candidate genes show significant similarity to genes encoding players in the abscisic acid- or Jasmonate-signalling pathways.

CONCLUSIONS

A set of candidate genes, concordantly expressed between field and greenhouse experiments, was identified. The candidates are associated with hormone signalling and specialized metabolite biosynthesis pathways indicating the involvement of these pathways in the response of the host to infection by H. fraxineus.

Toward Unifying Global Hotspots of Wild and Domesticated Biodiversity

Global biodiversity hotspots are areas containing high levels of species richness, endemism and threat. Similarly, regions of agriculturally relevant diversity have been identified where many domesticated plants and animals originated, and co-occurred with their wild ancestors and relatives. The agro-biodiversity in these regions has, likewise, often been considered threatened. Biodiversity and agro-biodiversity hotspots partly overlap, but their geographic intricacies have rarely been investigated together. Here we review the history of these two concepts and explore their geographic relationship by analysing global distribution and human use data for all plants, and for major crops and associated wild relatives. We highlight a geographic continuum between agro-biodiversity hotspots that contain high richness in species that are intensively used and well known by humanity (i.e., major crops and most viewed species on Wikipedia) and biodiversity hotspots encompassing species that are less heavily used and documented (i.e., crop wild relatives and species lacking information on Wikipedia). Our contribution highlights the key considerations needed for further developing a unifying concept of agro-biodiversity hotspots that encompasses multiple facets of diversity (including genetic and phylogenetic) and the linkage with overall biodiversity. This integration will ultimately enhance our understanding of the geography of human-plant interactions and help guide the preservation of nature and its contributions to people.

It takes three to tango: the importance of microbes, host plant, and soil management to elucidate manipulation strategies for the plant microbiome

The world’s population is expected to grow to almost 10 billion by 2050, placing unprecedented demands on agriculture and natural resources. The risk in food security is also aggravated by climate change and land degradation, which compromise agricultural productivity. In recent years, our understanding of the role of microbial communities on ecosystem functioning, including plant-associated microbes, has advanced considerably. Yet, translating this knowledge into practical agricultural technologies is challenged by the intrinsic complexity of agroecosystems. Here, we review current strategies for plant microbiome manipulation, classifying them into three main pillars: (i) introducing and engineering microbiomes, (ii) breeding and engineering the host plant, and (iii) selecting agricultural practices that enhance resident soil and plant-associated microbial communities. In each of these areas, we analyze current trends in research, as well as research priorities and future perspectives.

Host Range of Hymenoscyphus fraxineus in Slovak Arboreta

The health of 34 different Fraxinus taxa in association with the pathogenic fungus Hymenoscyphus fraxineus was assessed in four Slovak arboreta. Averaged across all arboreta, nearly one-quarter (24.9%) of all evaluated trees showed ash dieback symptoms. The damage was most serious on the common ash F. excelsior, a native species. The percentage of dead trees did not exceed 2% for all evaluated trees. Generally, ash trees of all ages were affected, though the intensity of the damage varied among the sites. The identity of H. fraxineus was confirmed by conventional PCR targeting the internal transcribed spacer (ITS) sequences of the nuclear ribosomal DNA, as well as the 18S gene/ITS-2 region of the rDNA operon. In Slovakia, the pathogen has expanded its host range from native species not only to their ornamental cultivars, but also to introduced North American (F. cinerea, F. latifolia, F. pennsylvanica, F. quadrangulata) and Asian (F. bungeana, F. chinensis ssp. rhynchophylla, F. mandshurica) ash species. H. fraxineus was also observed on the previous year’s leaf petioles of the native European species F. ornus, considered a weakly susceptible host. In Slovak arboreta, H. fraxineus was found on 23 Fraxinus taxa; 21 of them represent first records for the country. F. bungeana is recorded as a new host species of H. fraxineus.