Molecular markers for tolerance of European ash (Fraxinus excelsior) to dieback disease identified using Associative Transcriptomics

Fraxinus excelsior updated long-read genome reveals the importance of MADS-box genes in tolerance mechanisms against ash dieback

Ash dieback caused by the fungus Hymenoscyphus fraxineus has devastated the European ash tree population since it arrived in Europe in 1992. Great effort has been put into breeding programs to increase the genetic diversity of ash trees and find heritable genetic markers associated with resistance, or tolerance mechanisms, to ash dieback. To facilitate identification of molecular markers, we used Oxford Nanopore Technologies combined with Illumina sequencing to obtain an accurate and contiguous ash genome. We used this genome to reanalyze transcriptome data from a Danish ash panel of 182 tree accessions. Using associative transcriptomics, we identified 175 gene expression markers, including 11 genes annotated as dormancy MADS-box transcription factors which are associated with ash bud dormancy, flowering, and senescence. We hypothesize that tolerant trees both break dormancy earlier in the year by increasing the expression of flowering-related SOC1 MADS-box and reducing the expression of SVP-like MADS-box, whilst also accelerating senescence by increasing the expression of JOINTLESS MADS-box genes. DNA methylation differences in the promoters of MADS-box genes between 1 tolerant and 1 susceptible tree indicate potential epigenetic regulation of these traits.

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.

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.

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.

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.

How genomics can help biodiversity conservation

The availability of public genomic resources can greatly assist biodiversity assessment, conservation, and restoration efforts by providing evidence for scientifically informed management decisions. Here we survey the main approaches and applications in biodiversity and conservation genomics, considering practical factors, such as cost, time, prerequisite skills, and current shortcomings of applications. Most approaches perform best in combination with reference genomes from the target species or closely related species. We review case studies to illustrate how reference genomes can facilitate biodiversity research and conservation across the tree of life. We conclude that the time is ripe to view reference genomes as fundamental resources and to integrate their use as a best practice in conservation genomics.

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.

Gold Nanostars Combined with the Searched Antibody for Targeted Oral Squamous Cell Carcinoma Therapy

Oral squamous cell carcinoma (OSCC) is the most common cancer in the oral and maxillofacial region. Due to the special physiological and anatomical position of the oral cavity, the disease often has a significant impact on the chewing, swallowing, language, and breathing functions of patients. In recent years, with the development of medical molecular biology, molecular targeted therapy has received increasing clinical attention and has gradually become a new method for the treatment of malignant tumors. In this research, gold nanostars with a high photothermal effect combined with the searched targeted antibody were used for OSCC therapy. We use the data set in the public database and construct a gene co-expression module by weighted gene co-expression network analysis (WGCNA). It was found that the turquoise module and the midnight blue module had the greatest connection to tumorigenesis. Cytoscape software was used to analyze the important modules, and the top 10 genes of each module were selected; the survival analysis of the top 10 genes was carried out by gene expression profiling interactive analysis (GEPIA), which indicated that these genes (SERPINH1, MMP11, ADAM12, FADS3, SLC36A2, C1QTNF7, SCRG1, and APOBEC2) have statistical significance as key genes that are related to the tumorigenesis of OSCC. Then, the anti-SERPINH1 antibody targeted to SERPINH1 was chosen as the inhibitor and combined with gold nanostars for photothermal assisted targeted therapy. Thus, the searched key genes can be regarded as biomarkers and therapeutic targets for further precise diagnosis.

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.

Transcriptomics of Biostimulation of Plants Under Abiotic Stress

Plant biostimulants are compounds, living microorganisms, or their constituent parts that alter plant development programs. The impact of biostimulants is manifested in several ways: via morphological, physiological, biochemical, epigenomic, proteomic, and transcriptomic changes. For each of these, a response and alteration occur, and these alterations in turn improve metabolic and adaptive performance in the environment. Many studies have been conducted on the effects of different biotic and abiotic stimulants on plants, including many crop species. However, as far as we know, there are no reviews available that describe the impact of biostimulants for a specific field such as transcriptomics, which is the objective of this review. For the commercial registration process of products for agricultural use, it is necessary to distinguish the specific impact of biostimulants from that of other legal categories of products used in agriculture, such as fertilizers and plant hormones. For the chemical or biological classification of biostimulants, the classification is seen as a complex issue, given the great diversity of compounds and organisms that cause biostimulation. However, with an approach focused on the impact on a particular field such as transcriptomics, it is perhaps possible to obtain a criterion that allows biostimulants to be grouped considering their effects on living systems, as well as the overlap of the impact on metabolism, physiology, and morphology occurring between fertilizers, hormones, and biostimulants.

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.

A Brassica napus Reductase Gene Dissected by Associative Transcriptomics Enhances Plant Adaption to Freezing Stress

Cold treatment (vernalization) is required for winter crops such as rapeseed (Brassica napus L.). However, excessive exposure to low temperature (LT) in winter is also a stress for the semi-winter, early-flowering rapeseed varieties widely cultivated in China. Photosynthetic efficiency is one of the key determinants, and thus a good indicator for LT tolerance in plants. So far, the genetic basis underlying photosynthetic efficiency is poorly understood in rapeseed. Here the current study used Associative Transcriptomics to identify genetic loci controlling photosynthetic gas exchange parameters in a diversity panel comprising 123 accessions. A total of 201 significant Single Nucleotide Polymorphisms (SNPs) and 147 Gene Expression Markers (GEMs) were detected, leading to the identification of 22 candidate genes. Of these, Cab026133.1, an ortholog of the Arabidopsis gene AT2G29300.2 encoding a tropinone reductase (BnTR1), was further confirmed to be closely linked to transpiration rate. Ectopic expressing BnTR1 in Arabidopsis plants significantly increased the transpiration rate and enhanced LT tolerance under freezing conditions. Also, a much higher level of alkaloids content was observed in the transgenic Arabidopsis plants, which could help protect against LT stress. Together, the current study showed that AT is an effective approach for dissecting LT tolerance trait in rapeseed and that BnTR1 is a good target gene for the genetic improvement of LT tolerance in plant.

Advanced spectroscopy-based phenotyping offers a potential solution to the ash dieback epidemic

Natural and urban forests worldwide are increasingly threatened by global change resulting from human-mediated factors, including invasions by lethal exotic pathogens. Ash dieback (ADB), incited by the alien invasive fungus Hymenoscyphus fraxineus, has caused large-scale population decline of European ash (Fraxinus excelsior) across Europe, and is threatening to functionally extirpate this tree species. Genetically controlled host resistance is a key element to ensure European ash survival and to restore this keystone species where it has been decimated. We know that a low proportion of the natural population of European ash expresses heritable, quantitative resistance that is stable across environments. To exploit this resource for breeding and restoration efforts, tools that allow for effective and efficient, rapid identification and deployment of superior genotypes are now sorely needed. Here we show that Fourier-transform infrared (FT-IR) spectroscopy of phenolic extracts from uninfected bark tissue, coupled with a model based on soft independent modelling of class analogy (SIMCA), can robustly discriminate between ADB-resistant and susceptible European ash. The model was validated with populations of European ash grown across six European countries. Our work demonstrates that this approach can efficiently advance the effort to save such fundamental forest resource in Europe and elsewhere.

The Endophytic Mycobiome of European Ash and Sycamore Maple Leaves - Geographic Patterns, Host Specificity and Influence of Ash Dieback

The European ash (Fraxinus excelsior) is threatened by the introduced ascomycete Hymenoscyphus fraxineus, the causal agent of ash dieback. Endophytic fungi are known to modulate their host's resistance against pathogens. To understand possible consequences of ash dieback on the endophytic mycobiome, F. excelsior leaves were collected in naturally regenerated forests and the fungal communities analyzed by classic culture and Illumina amplicon sequencing using a newly developed and validated fungal-specific primer. Collections were done in the area infested by ash dieback north of the Alps, and in the disease free area on the south side. Sycamore maple (Acer pseudoplatanus) was additionally collected, as well as the flowering ash (F. ornus), which occurs naturally in the south and shows tolerance to ash dieback. Both cultivation and amplicon sequencing revealed characteristic endophytic fungal communities dominated by several strictly host specific Venturia species. On A. pseudoplatanus, a hitherto undescribed Venturia species was identified. Due to its dominance on F. excelsior, V. fraxini is unlikely to go extinct in case of reduced host densities. A majority of species was not strictly host specific and is therefore likely less affected by ash dieback in the future. Still, shifts in community structure and loss of genetic diversity cannot be excluded. The potentially endangered endophyte Hymenoscyphus albidus was rarely found. In addition to host specificity, species with preferences for leaf laminae or petioles were found. We also detected considerable geographical variation between sampling sites and clear differences between the two sides of the Alps for endophytes of F. excelsior, but not A. pseudoplatanus. Since sycamore maple is not affected by an epidemic, this could point toward an influence of ash dieback on ash communities, although firm conclusions are not possible because of host preferences and climatic differences. Furthermore, the mycobiota of F. excelsior trees with or without dieback symptoms were compared, but no clear differences were detected. Besides methodical refinement, our study provides comprehensive data on the ash mycobiome that we expect to be subject to changes caused by an emerging disease of the host tree.

Survival of European Ash Seedlings Treated with Phosphite after Infection with the Hymenoscyphus fraxineus and Phytophthora Species

The European Fraxinus species are threatened by the alien invasive pathogen Hymenoscyphus fraxineus, which was introduced into Poland in the 1990s and has spread throughout the European continent, causing a large-scale decline of ash. There are no effective treatments to protect ash trees against ash dieback, which is caused by this pathogen, showing high variations in susceptibility at the individual level. Earlier studies have shown that the application of phosphites could improve the health of treated seedlings after artificial inoculation with H. fraxineus. Three-year-old F. excelsior seedlings were inoculated with the following pathogens: a H. fraxineus, Phytophthora species mixture (P. plurivora, P. megasperma, and P. taxon hungarica), in combination with two pathogens and mock-inoculated as the control, and then either watered or treated with ammonium phosphite (Actifos). Results showed significant differences in the survival of seedlings and symptoms of disease development among the treatments. Chlorophyll-a fluorescence parameters indicated a decrease in photosynthetic efficiency in infected plants, suggesting that they were under strong biotic stress, but none of the parameters could be used as a reliable bioindicator for ash decline disease. The application of Actifos enhanced the production of triterpenes (ursolic and oleanolic acid), and decreased the production of phenols (tyrosol) and sterols (β-sitosterol) in seedlings infected with H. fraxineus. Treatment with Actifos caused seedlings to enhance their response to pathogen(s) attack and increase their survival probability.

Genome-wide epigenetic variation among ash trees differing in susceptibility to a fungal disease

Background

European ash trees (Fraxinus excelsior) are currently threatened by ash dieback (ADB) caused by the fungus Hymenoscyphus fraxineus but a small percentage of the population possesses natural low susceptibility. The genome of a European ash tree has recently been sequenced. Here, we present whole genome DNA methylation data for two F. excelsior genotypes with high susceptibility to ADB, and two genotypes with low susceptibility, each clonally replicated. We also include two genotypes of Manchurian ash (F. mandshurica), an ash species which has co-evolved with H. fraxineus and also has low susceptibility to ADB.

Results

In F. excelsior, we find an average methylation level of 76.2% in the CG context, 52.0% in the CHG context, and 13.9% in the CHH context; similar levels to those of tomato. We find higher methylation in transposable elements as opposed to non-mobile elements, and high densities of Non-Differentially Methylation Positions (N-DMPs) in genes with housekeeping functions. Of genes putatively duplicated in whole genome duplication (WGD) events, an average of 25.9% are differentially methylated in at least one cytosine context, potentially indicative of unequal silencing. Variability in methylation patterns exists among clonal replicates, and this is only slightly less than the variability found between different genotypes. Of twenty genes previously found to have expression levels associated with ADB susceptibility, we find only two of these have differential methylation between high and low susceptibility F. excelsior trees. In addition, we identify 1683 significant Differentially Methylated Regions (DMRs) (q-value< 0.001) between the high and low susceptibility genotypes of F. excelsior trees, of which 665 remain significant when F. mandshurica samples are added to the low susceptibility group.

Conclusions

We find a higher frequency of differentially methylated WGD-derived gene duplicates in ash than other plant species previously studied. We also identify a set of genes with differential methylation between genotypes and species with high versus low susceptibility to ADB. This provides valuable foundational data for future work on the role that epigenetics may play in gene dosage compensation and susceptibility to ADB in ash.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4874-8) contains supplementary material, which is available to authorized users.

Genotyping-by-sequencing through transcriptomics: implementation in a range of crop species with varying reproductive habits and ploidy levels

The application of genomics in crops has the ability to significantly improve genetic gain for agriculture. Many marker-dense tools have been developed, but few have seen broad adoption in plant genomics due to issues of significant variations of genome size, levels of ploidy, single nucleotide polymorphism (SNP) frequency and reproductive habit. When combined with limited breeding activities, small research communities and scant sequence resources, the suitability of popular systems is often suboptimal and routinely fails to effectively balance cost-effectiveness and sample throughput. Genotyping-by-sequencing (GBS) encompasses a range of protocols including resequencing of the transcriptome. This study describes a skim GBS-transcriptomics (GBS-t) approach developed to be broadly applicable, cost-effective and high-throughput while still assaying a significant number of SNP loci. A range of crop species with differing levels of ploidy and degree of inbreeding/outbreeding were chosen, including perennial ryegrass, a diploid outbreeding forage grass; phalaris, a putative segmental allotetraploid outbreeding forage grass; lentil, a diploid inbreeding grain legume; and canola, an allotetraploid partially outbreeding oilseed. GBS-t was validated as a simple and largely automated, cost-effective method which generates sufficient SNPs (from 89 738 to 231 977) with acceptable levels of missing data and even genome coverage from c. 3 million sequence reads per sample. GBS-t is therefore a broadly applicable system suitable for many crops, offering advantages over other systems. The correct choice of subsequent sequence analysis software is important, and the bioinformatics process should be iterative and tailored to the specific challenges posed by ploidy variation and extent of heterozygosity.

Validation of an updated Associative Transcriptomics platform for the polyploid crop species Brassica napus by dissection of the genetic architecture of erucic acid and tocopherol isoform variation in seeds

An updated platform was developed to underpin association genetics studies in the polyploid crop species Brassica napus (oilseed rape). Based on 1.92 × 10¹² bases of leaf mRNAseq data, functional genotypes, comprising 355 536 single-nucleotide polymorphism markers and transcript abundance were scored across a genetic diversity panel of 383 accessions using a transcriptome reference comprising 116 098 ordered coding DNA sequence (CDS) gene models. The use of the platform for Associative Transcriptomics was first tested by analysing the genetic architecture of variation in seed erucic acid content, as high-erucic rapeseed oil is highly valued for a variety of applications in industry. Known loci were identified, along with a previously undetected minor-effect locus. The platform was then used to analyse variation for the relative proportions of tocopherol (vitamin E) forms in seeds, and the validity of the most significant markers was assessed using a take-one-out approach. Furthermore, the analysis implicated expression variation of the gene Bo2g050970.1, an orthologue of VTE4 (which encodes a γ-tocopherol methyl transferase converting γ-tocopherol into α-tocopherol) associated with the observed trait variation. The establishment of the first full-scale Associative Transcriptomics platform for B. napus enables rapid progress to be made towards an understanding of the genetic architecture of trait variation in this important species, and provides an exemplar for other crops.

Ash leaf metabolomes reveal differences between trees tolerant and susceptible to ash dieback disease

European common ash, Fraxinus excelsior, is currently threatened by Ash dieback (ADB) caused by the fungus, Hymenoscyphus fraxineus. To detect and identify metabolites that may be products of pathways important in contributing to resistance against H. fraxineus, we performed untargeted metabolomic profiling on leaves from five high-susceptibility and five low-susceptibility F. excelsior individuals identified during Danish field trials. We describe in this study, two datasets. The first is untargeted LC-MS metabolomics raw data from ash leaves with high-susceptibility and low-susceptibility to ADB in positive and negative mode. These data allow the application of peak picking, alignment, gap-filling and retention-time correlation analyses to be performed in alternative ways. The second, a processed dataset containing abundances of aligned features across all samples enables further mining of the data. Here we illustrate the utility of this dataset which has previously been used to identify putative iridoid glycosides, well known anti-herbivory terpenoid derivatives, and show differential abundance in tolerant and susceptible ash samples.

A first assessment of Fraxinus excelsior (common ash) susceptibility to Hymenoscyphus fraxineus (ash dieback) throughout the British Isles

Ash dieback (ADB), caused by Hymenoscyphus fraxineus, has severely damaged a large proportion of ash trees (Fraxinus excelsior) in continental Europe. We have little damage data for the British Isles where the disease was found only five years ago in the Southeast, and is still spreading. A large-scale screening trial to evaluate ADB damage to provenances of F. excelsior sourced from throughout the British Isles was planted in 2013 in the southeast of England. In 2016, we scored trees by their level of ADB damage observed in field at the two worst affected (based on assessments in 2015) of the 14 sites. Significant differences were found in average ADB damage among planting sites and seed source provenances. Trees from certain provenances in Scotland were the least damaged by ADB, whereas trees from Wales and Southeast England were the most badly damaged in both trial sites. Thus the levels of ADB damage currently seen in ash populations in Southeast England may not be an accurate predictor of the damage expected in future throughout the British Isles. Given all provenances contained some healthy trees, a breeding programme to produce genetically variable native ash tree populations with lower ADB susceptibility may be feasible.

Gene flow of common ash (Fraxinus excelsior L.) in a fragmented landscape

Gene flow dynamics of common ash (Fraxinus excelsior L.) is affected by several human activities in Central Europe, including habitat fragmentation, agroforestry expansion, controlled and uncontrolled transfer of reproductive material, and a recently introduced emerging infectious disease, ash dieback, caused by Hymenoscyphus fraxineus. Habitat fragmentation may alter genetic connectivity and effective population size, leading to loss of genetic diversity and increased inbreeding in ash populations. Gene flow from cultivated trees in landscapes close to their native counterparts may also influence the adaptability of future generations. The devastating effects of ash dieback have already been observed in both natural and managed populations in continental Europe. However, potential long-term effects of genetic bottlenecks depend on gene flow across fragmented landscapes. For this reason, we studied the genetic connectivity of ash trees in an isolated forest patch of a fragmented landscape in Rösenbeck, Germany. We applied two approaches to parentage analysis to estimate gene flow patterns at the study site. We specifically investigated the presence of background pollination at the landscape level and the degree of genetic isolation between native and cultivated trees. Local meteorological data was utilized to understand the effect of wind on the pollen and seed dispersal patterns. Gender information of the adult trees was considered for calculating the dispersal distances. We found that the majority of the studied seeds (55-64%) and seedlings (75-98%) in the forest patch were fathered and mothered by the trees within the same patch. However, we determined a considerable amount of pollen flow (26-45%) from outside of the study site, representing background pollination at the landscape level. Limited pollen flow was observed from neighbouring cultivated trees (2%). Both pollen and seeds were dispersed in all directions in accordance with the local wind directions. Whereas there was no positive correlation between pollen dispersal distance and wind speed, the correlation between seed dispersal distance and wind speed was significant (0.71, p < 0.001), indicating that strong wind favours long-distance dispersal of ash seeds. Finally, we discussed the implications of establishing gene conservation stands and the use of enrichment planting in the face of ash dieback.

Exploring public perceptions of solutions to tree diseases in the UK: Implications for policy-makers

Tree diseases are on the increase in many countries and the implications of their appearance can be political, as well as ecological and economic. Preventative policy approaches to tree diseases are difficult to formulate because dispersal pathways for pest and pathogens are numerous, poorly known and likely to be beyond human management control. Genomic techniques could offer the quickest and most predictable approach to developing a disease tolerant native ash. The population of European Ash (Fraxinus Excelsior) has suffered major losses in the last decade, due to the onset of Hymenoscyphus fraxineus (previously called Chalara Fraxinea) commonly known in the UK as ash dieback. This study presents evidence on the public acceptability of tree-breed solutions to the spread of Chalara, with the main aim to provide science and policy with an up-stream 'steer' on the likely public acceptability of different tree breeding solutions. The findings showed that whilst there was a firm anti-GM and 'we shouldn't tamper with nature' attitude among UK publics, there was an equally firm and perhaps slightly larger pragmatic attitude that GM (science and technology) should be used if there is a good reason to do so, for example if it can help protect trees from disease and help feed the world. The latter view was significantly stronger among younger age groups (Millennials), those living in urban areas and when the (GM)modified trees were destined for urban and plantation, rather than countryside settings. Overall, our findings suggest that the UK government could consider genomic solutions to tree breeding with more confidence in the future, as large and influential publics appear to be relaxed about the use of genomic techniques to increase tolerance of trees to disease.

Developing publicly acceptable tree health policy: public perceptions of tree-breeding solutions to ash dieback among interested publics in the UK

The UK needs to develop effective policy responses to the spread of tree pathogens and pests. This has been given the political urgency following the media and other commentary associated with the arrival of a disease that causes 'dieback' of European Ash (Fraxinus excelsior) - a tree species with deep cultural associations. In 2014 the UK government published a plant biosecurity strategy and linked to this invested in research to inform policy. This paper reports the findings of a survey of informed UK publics on the acceptability of various potential strategies to deal with ash dieback, including "no action". During the summer of 2015, we conducted a face-to-face survey of 1152 respondents attending three major countryside events that attract distinct publics interested in the countryside: landowners & land managers; naturalists and gardeners. We found that UK publics who are likely to engage discursively and politically (through letter writing, petitions etc.) with the issue of ash dieback a) care about the issue, b) want an active response, c) do not really distinguish between ash trees in forestry or ecological settings, and d) prefer traditional breeding solutions. Further that e) younger people and gardeners are open to GM breeding techniques, but f) the more policy-empowered naturalists are more likely to be anti-GM. We suggest that these findings provide three 'steers' for science and policy: 1) policy needs to include an active intervention component involving the breeding of disease-tolerant trees, 2) that the development of disease tolerance using GM-technologies could be part of a tree-breeding policy, and 3) there is a need for an active dialogue with publics to manage expectations on the extent to which science and policy can control tree disease or, put another way, to build acceptability for the prospect that tree diseases may have to run their course.

Genome sequence and genetic diversity of European ash trees

Ash trees (genus Fraxinus, family Oleaceae) are widespread throughout the Northern Hemisphere, but are being devastated in Europe by the fungus Hymenoscyphus fraxineus, causing ash dieback, and in North America by the herbivorous beetle Agrilus planipennis. Here we sequence the genome of a low-heterozygosity Fraxinus excelsior tree from Gloucestershire, UK, annotating 38,852 protein-coding genes of which 25% appear ash specific when compared with the genomes of ten other plant species. Analyses of paralogous genes suggest a whole-genome duplication shared with olive (Olea europaea, Oleaceae). We also re-sequence 37 F. excelsior trees from Europe, finding evidence for apparent long-term decline in effective population size. Using our reference sequence, we re-analyse association transcriptomic data, yielding improved markers for reduced susceptibility to ash dieback. Surveys of these markers in British populations suggest that reduced susceptibility to ash dieback may be more widespread in Great Britain than in Denmark. We also present evidence that susceptibility of trees to H. fraxineus is associated with their iridoid glycoside levels. This rapid, integrated, multidisciplinary research response to an emerging health threat in a non-model organism opens the way for mitigation of the epidemic.