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

Caterina Villari, Arnaud Dowkiw, Rasmus Enderle, Marjan Ghasemkhani, Thomas Kirisits, Erik D. Kjaer, Diana Marciulyniene, Lea V. Mckinney, Berthold Metzler, Facundo Munoz, et al.

► To cite this version:

Caterina Villari, Arnaud Dowkiw, Rasmus Enderle, Marjan Ghasemkhani, Thomas Kirisits, et al. Advanced spectroscopy-based phenotyping offers a potential solution to the ash dieback epidemic. Scientific Reports, 2018, 8, 10.1038/s41598-018-35770-0. hal-02623741

HAL Id: hal-02623741

https://hal.inrae.fr/hal-02623741

Submitted on 26 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

<figure><img src="image_2.png" /></figure>

OPEN

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

Received: 19 January 2018

Accepted: 8 November 2018

Published online: 28 November 2018

Caterina Villari<fnref n="1" />,², Arnaud Dowkiw³, Rasmus Enderle⁴,⁵, Marjan Ghasemkhani⁶,
Thomas Kirisits⁷, Erik D. Kjær<fnref n="1" />,⁸, Diana Marčiulynienė⁹, Lea V. McKinney⁸,
Berthold Metzler⁴, Facundo Muñoz<fnref n="1" />,⁳, Lene R. Nielsen<fnref n="1" />,⁸, Alfas Pliūra⁹, Lars-Göran Stener¹⁰,
Vytautas Suchockas⁹, Luis Rodriguez-Saona¹¹, Pierluigi Bonello<fnref n="1" /> & Michelle Cleary⁶

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.

Invasive, alien tree pathogens threaten biodiversity, ecosystem services, and commercial forestry on a global scale<fnref n="1" />⁻³. Such bioinvasions are increasing at an unprecedented rate, in part because of higher global connectivity⁴, which facilitates unintended long-distance movement of pathogens into regions outside their historical distribution range⁵, and because of climate change, which causes host maladaptation⁶ and emergence of new pathogens

<footnote n="1">¹Department of Plant Pathology, The Ohio State University, 201 Kottman Hall, 2021 Coffey Road, 43210, Columbus, Ohio, United States. ²Warnell School of Forestry & Natural Resources, University of Georgia, 180 E Green Street, 30602, Athens, Georgia, United States. ³Institut National de la Recherche Agronomique (INRA), UMR 0588 BioForA, 2163 Avenue de la Pomme de Pin, 45075, Orléans Cedex 2, France. ⁴Forest Research Institute Baden-Wuerttemberg, Department Forest Protection, Wonnhaldestrasse 4, 79100, Freiburg, Germany. ⁵Institute for Plant Protection in Horticulture and Forests, Federal Research Centre for Cultivated Plants (Julius Kühn Institute), Messeweg 11/12, 38104, Braunschweig, Germany. ⁶Swedish University of Agricultural Sciences (SLU), Southern Swedish Forest Research Centre, Sundsvägen 3, 23053, Alnarp, Sweden. ⁷Institute of Forest Entomology, Forest Pathology and Forest Protection (IFFF), Department of Forest and Soil Sciences, University of Natural Resources and Life Sciences, Vienna (BOKU), Peter-Jordan-Straße 82, 1190, Vienna, Austria. ⁸Department of Geosciences and Natural Resource Management, University of Copenhagen, Rolighedsvej 23, 1958, Frederiksberg C., Denmark. ⁹Lithuanian Research Centre for Agriculture and Forestry, Institute of Forestry, Liepu 1, LT53101, Girionys, Kaunas district, Lithuania. ¹⁰SKOGFORSK - The Forest Research Institute, Ekebo 2250, 26890, Svalöv, Sweden. ¹¹Department of Food Science and Technology, The Ohio State University, Parker Food Science and Technology, 2015 Fyffe Road, 43210, Columbus, Ohio, United States. Correspondence and requests for materials should be addressed to C.V. (email: cvillari@uga.edu) or M.C. (email: Michelle.Cleary@slu.se)</footnote>
through alteration of host-microbe interactions⁷. These disturbances have the potential to cause massive and irreversible damage⁸ by eliminating keystone tree species in many areas of the world, and permanently altering trophic structures, nutrient dynamics¹ and primary productivity² of forest communities. Such wide ranging ecological impacts can compromise the maintenance of ecosystem services upon which humans rely, including those associated with reduced incidence of human morbidities⁹.

Genetically controlled resistance is crucial to successfully manage natural and urban forests for resilience against alien invasive pathogens¹⁰. Molecular markers are invaluable to enhance screening for resistance, and can be applied to practical breeding programs with high precision and reductions in cost and time. However, marker-assisted selection has not been used for genetic improvement of forest trees to the same extent as agricultural crops¹¹. Indeed, efficient marker development requires accurate field testing, which in trees can be very slow due to the long time lag intrinsic in symptom development¹². This lack of efficient tools for rapid resistance phenotyping significantly hinders our ability to screen natural populations for conventional breeding of trees. Our ability to readily and reliably detect superior genotypes could enhance the success of current restoration efforts, or protect trees from logging or other activities associated with forest and urban landscape management, in support of *in situ* conservation.

We hypothesized that disease-resistant genotypes can be identified through chemometrics using vibrational spectroscopy, which is based on the absorption of infrared radiation resulting from fundamental molecular (bond) vibrations. This technique can rapidly fingerprint a wide range of biological samples using multivariate statistical classification models that identify and delineate target classes¹³⁻¹⁵. Considering that plant resistance against pests and pathogens essentially relies on host chemistry¹⁶, vibrational spectroscopy-based techniques hold vast potential in distinguishing between plant chemical phenotypes (chemotypes) that are genetically and epigenetically driven, and vary in disease susceptibility. Among these techniques, Fourier-transform infrared (FT-IR) spectroscopy has so far shown promising results in applications involving forest trees¹³, including delin-eating resistance to invasive pathogens such as *Ophiostoma novo-ulmi*, causal agent of Dutch elm disease¹⁷,¹⁸, and *Phytophthora ramorum*, causal agent of sudden oak death¹⁹. In the latter case, the technique was able to success-fully model tree resistance prior to infection.

At present, there is rising concern that European ash (Fraxinus excelsior), an important keystone species in natural plant communities protected under EU legislation²⁰, may be functionally extirpated from European for-ests by the alien invasive fungus, *Hymenoscyphus fraxineus*, causal agent of ash dieback (ADB)²¹,²². A low pro-portion of the natural population of European ash expresses heritable, quantitative resistance that is stable across environments²³⁻²⁵, and genomic solutions to uncover mechanisms associated with disease resistance have been explored as a means to accelerate breeding of trees with resistance against ADB²⁶,²⁷. In a recent paper, Sollars et al.²⁸ found an association between disease susceptibility and the levels of two putative iridoid glycosides in the leaves of *F. excelsior*, which suggests that chemotypes associated with this pathosystem can be targeted and tested using FT-IR spectroscopy and chemometric models. Consequently, our goal was to determine the feasibility and efficacy of FT-IR to phenotype European ash for resistance to ADB.

## Results and Discussion

We analysed the Fourier transform mid-infrared spectral region¹⁴,²⁹ of phenolic extracts¹⁶ of uninfected tissue samples from 76 different genotypes collected across six European countries (Austria, Denmark, France, Germany, Lithuania and Sweden) (Fig. 1) that were previously phenotyped as having either low, intermediate or high susceptibility to *H. fraxineus*²⁴,²⁵,³⁰⁻³⁶ (Table 1). We targeted phenolics due to the established role of this class of secondary metabolites in general plant defense¹⁶,³⁷, and the successful prior phenotyping efforts using this technique in other pathosystems¹⁹. Samples comprised both leaves and twig bark, which included a thin outer bark (<1 mm thick), cortex, secondary phloem, and cambial tissues. To test the robustness and reliability of our approach, we ensured broad heterogeneity of the samples within each tissue type by sampling trees at different developmental stages and environmental conditions across locations. Up to three ramets per genotype were sampled, for a total of 134 trees across the six countries (Table 1). To be practically useful, chemotypes should be associated with constitutive composition and levels of specialized phytochemicals. Therefore, it was critical to ensure that *H. fraxineus* was not present in plant tissues. Despite observations at the time of sample collection of some ripe apothecia on pseudosclerotial leaf rachises in the litter of moist microsites at the Austrian location, all samples from all locations were confirmed free of the pathogen via PCR³⁸. All trees may have harboured infections in other parts of the tree, prior to our observations in this study, and thus we cannot exclude possible systemic effects of these other infections on the chemotypes. However, previous work in the coast live oak – sudden oak death pathosystem has shown that healthy tissues sampled away from active cankers were chemically indistinguishable from healthy tissues taken from asymptomatic and presumably uninfected trees³⁹.

We processed FT-IR chemical fingerprints using a soft independent modelling of class analogy (SIMCA) chemometric approach to discriminate between different resistance phenotypes¹³. In order to resolve overlapping peaks, minimize background and improve the model predictions, spectral data were pre-processed via the standard normal variate function⁴⁰, and smoothed and transformed into their second derivative using a Savitzky-Golay polynomial filter⁴¹. We initially analysed both leaf and twig bark tissues, for a total of 131 and 112 samples, respectively. Preliminary observations of the SIMCA 3D class projection plots of resistant and susceptible samples showed that geographic location of the trees strongly affected the chemistry of the leaves (Fig. 2). This difference may be attributable to a higher sensitivity of foliage to solar irradiation⁴², temperature⁴³ and nutrient or water availability⁴⁴ associated with the different latitudes, geographic locations, and microclimates. However, this effect was not evident in twig bark tissues (Fig. 2), therefore, we built the chemometric model using only the FT-IR spectra of the twig bark samples. We do not exclude the possibility that resistance is expressed at the leaf level³¹; rather, we suggest that environmental variation masked any possible chemical signature in the leaf tissues. Indeed, in previous studies²⁸, leaf constitutive iridoid glycosides could discriminate among European
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Map of the sampling sites. Uninfected bark and leaf samples were collected from a total of 76 *Fraxinus excelsior* genotypes of known susceptibility to *Hymenophyllum fraxineus* in six European countries: Austria, Denmark, France, Germany, Lithuania, and Sweden. Sampling sites (red mapping pins) are overlaid on the natural distribution map of *F. excelsior* (sky-blue) (EUFORGEN 2009, www.euforgen.org.).</figcaption></figure>

<table><caption>**Table 1.** Detailed information on the *Fraxinus excelsior* genotypes analysed in the study. *For all countries except Lithuania, genotypes in each trial originated in the same country where the trial is located. Genotypes in the Lithuanian trial originated in the Czech Republic, Germany, Ireland and Lithuania. **Literature reference number. Sus., susceptible; Int., intermediate; Res., resistant.</caption><thead><tr><th rowspan="2">Country*</th><th rowspan="2">Location</th><th colspan="4">Trial Details</th><th rowspan="2">Date of sample collection in 2015</th><th colspan="3">Number of genotypes sampled per susceptibility class (for clonal trials, the number of ramets per genotype is given in parentheses)</th></tr><tr><th>Type of genetic trial</th><th>Coordinates</th><th>Elevation (m asl)</th><th>Year established</th><th>Spacing of trees (m)</th><th>Sus.</th><th>Int.</th><th>Res.</th></tr></thead><tbody><tr><td>Austria<sup>31,36**</sup></td><td>Feldkirchen an der Donau</td><td>Clonal seed orchard</td><td>48°20'11.4" N 14°02'53.2" E</td><td>264</td><td>1993</td><td>7.5 x 8.6</td><td>9 June</td><td>7 (2)</td><td></td><td>7 (2)</td></tr><tr><td>Denmark<sup>32</sup></td><td>Tuse næs</td><td>Clonal seed orchard</td><td>55°45' 58.0" N 11°42' 47.4" E</td><td>22</td><td>1998</td><td>3.0 x 6.0</td><td>2, 4 June</td><td>3 (3)</td><td>2 (3)</td><td>3 (3)</td></tr><tr><td>France<sup>24</sup></td><td>Devecey</td><td>Provenance and progeny trial</td><td>47°19'31.5" N 06°01'54.1" E</td><td>250</td><td>1995</td><td>4.0 x 4.0</td><td>18 June</td><td>7</td><td>3</td><td>7</td></tr><tr><td>Germany<sup>35</sup></td><td>Weisweil</td><td>Provenance trial</td><td>48°11'29.7" N 07°42'02.5" E</td><td>173</td><td>2005</td><td>2.0 x 2.0</td><td>19 May</td><td>5</td><td></td><td>5</td></tr><tr><td>Lithuania<sup>34</sup></td><td>Sasnava</td><td>Clonal collection</td><td>54°37'32.1" N 23°33'55.5" E</td><td>100</td><td>2012</td><td>6.0 x 5.4</td><td>2 June</td><td>2 (3), 2 (2)</td><td>2 (2)</td><td>3 (3), 1 (2), 1 (1)</td></tr><tr><td>Sweden<sup>30</sup></td><td>Snogeholm</td><td>Clonal seed orchard</td><td>55°32'33.8" N 13°42'22.7" E</td><td>50</td><td>1992</td><td>3.5 x 3.5</td><td>28 May</td><td>3 (1), 1 (2)</td><td>5 (2)</td><td>5 (2), 2 (2)</td></tr><tr><td colspan="3">Total number of genotypes</td><td></td><td></td><td></td><td></td><td></td><td>30</td><td>12</td><td>34</td></tr><tr><td colspan="3">Total number of trees</td><td></td><td></td><td></td><td></td><td></td><td>50</td><td>23</td><td>61</td></tr></tbody></table>
<figure><img src="image_4.png" /><figcaption>**Figure 2.** SIMCA 3D class projections. SIMCA 3D class projection plots for spectral data of *Fraxinus excelsior* leaf and twig bark tissue phenolic extracts analysed as a function of the resistance phenotype, but visualized according to either the sample geographic location or its resistance phenotype. Spectral data were pre-processed using the standard normal variate function and then smoothed and transformed into their second derivative. Two technical replicates were analysed separately. Clouds of black points indicate the 95% confidence interval for each class (i.e., resistance phenotypes) in each principal component direction (i.e., PC1, PC2 and PC3) projected into the three-factor principal component hyper-plane.</figcaption></figure>

ash resistance phenotypes. Secoiridoids have also been shown to be upregulated in the leaves of *F. excelsior* in response to treatment with a *H. fraxineus* toxin⁴⁵.

To optimize the model, clustering patterns of the resistance phenotypes were sought in the SIMCA 3D class projection plot (Fig. 2), while spectral regions with the highest discriminating power were identified using SIMCA discriminating power plots and Coomans plots⁴⁶. We used Coomans plots and 3D class projection plots also to identify outliers, which were removed from the model. After trimming the data in this way, the complete data set comprised a minimum of two and up to seven genotypes per susceptibility class per country, and up to three ramets per clone, for a total of 71 samples (42 resistant, 7 intermediate, and 22 susceptible). Of these, 75% of samples within each geographic location and susceptibility class were selected for the training data set to build our calibration model, while the remaining 25% were used as the testing data set for model validation¹³. Samples were randomly assigned to the groups; in doing so, 10 clones were represented by ramets in both the training and testing data sets, while eight clones were unique to the testing population. Therefore, the testing dataset verified the ability of the model to classify the resistance phenotype for both clones that were independent of the training population and those that were biological replicates of clones included in the training population. Incremental analysis of different models showed that inclusion of intermediate phenotypes in the calibration compromised the formulation and strength of the chemometric classification. This was to be expected, as SIMCA models work best when built on well-defined groups¹³,¹⁹. We hence decided to exclude intermediates from the training set and move all of them to the testing set.

The SIMCA calibration model that best discriminated between resistant and susceptible ash trees was a 3-factor model (Fig. 3a) obtained by including spectral regions from ~748 to 798 cm⁻¹ and from ~879 to 947 cm⁻¹ wavenumber (Fig. 3b), which primarily correspond to the C–H wagging of substituted benzenes of aromatic compounds⁴⁷. The highest discriminating power peak (discriminating power of 47.7) corresponded to ~895 cm⁻¹ wavenumber, which may also correspond to the wagging of the hydrogen on the C-1 position of the cellulose glucose ring⁴⁸. However, since the extraction protocol we adopted is highly specific for phenolic compounds⁴⁹, the presence of polysaccharides in the analysed extracts is very unlikely, and the ~895 cm⁻¹ wavenumber almost certainly corresponds to the C–H wagging of the aromatic hydrocarbon groups of phenols⁴⁷. The second highest discriminating power peak (discriminating power of 25.4) corresponded to ~770 cm⁻¹ wavenumber. Interclass distance between the groups was 2.2074, indicating good separation between the phenotypes, and 100% of extracts from both resistant and susceptible ash trees were correctly classified, showing high specificity of the model. The low number of factors included in the model (i.e., three) argues against model overfitting⁵⁰,
<figure><img src="image_5.png" /><figcaption>b</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>c</figcaption></figure>

<figure><img src="image_7.png" /></figure>

**Figure 3.** SIMCA Coomans plots and discriminating power plot. Panel a, SIMCA Coomans plot showing the relative, dimension-free distance between the samples of the training data set used to build the 3-factor (@3) calibration model designed to discriminate between ash dieback resistant (red diamonds) and susceptible (blue diamonds) *Fraxinus excelsior* trees. X-axis represents the distance from the resistant class, while y-axis represents the distance from the susceptible class. Two technical replicates were analysed separately, for a total of 92 spectra, corresponding to 48 biological replicates. Dashed lines indicate critical sample residual thresholds. Panel b, SIMCA discriminating power plot of the 3-factor calibration model. The discriminating power (black line) is overlaid on the second derivative, smoothed and standard normal variate transformed spectra. The SIMCA calibration model that best discriminated between resistant (red lines) and susceptible (blue lines) ash trees included spectral regions from ~748 to 798 cm⁻¹ and from ~879 to 947 cm⁻¹ wavenumber (highlighted in yellow). The black arrows point to regions of the spectra where the discrimination between the two resistance phenotypes is evident by visual inspection. Panel c, SIMCA Coomans plot showing the relative, dimension-free distance between the samples of the testing data set used to validate the 3-factor (@3) model. In addition to resistant (red diamonds) and susceptible (blue diamonds) trees, the testing data set included accessions of intermediate phenotype (green diamonds), based on field observations. Two technical replicates were analysed separately, for a total of 44 spectra, corresponding to 23 biological replicates randomly selected from each of the six European countries. Dashed lines indicate critical sample residual thresholds.
with factors 1–3 explaining 97.87%, 1.21% and 0.53%, respectively, of the variability in the susceptible class, and 97.10%, 1.42%, and 1.00%, respectively, of the variability in the resistant class.

Validation of the chemical marker-based statistical model is a crucial step to ensure accurate predictions and confirm the applicability of the model at a larger scale¹³. We validated our SIMCA model with the testing data set, and 100% of the ten resistant and six susceptible ash ramets randomly selected from each of the six European countries (including nine representing biological replicates of the clones included in the training data population) were correctly identified as belonging solely to their phenotype group (Fig. 3c). Four of the intermediate samples, corresponding to four different genotypes, were identified as potentially belonging to both classes, which correctly reflects their intermediate phenotype. However, the remaining three samples, corresponding to two different genotypes, were identified as belonging solely to the resistant group (Fig. 3c). One of the intermediate genotypes had different ramets that were classified in either groups (resistant and intermediate). Such overestimation of the resistant samples demonstrates that while the model may not be perfect with intermediates, it clearly separated the most susceptible samples from the rest. Furthermore, the overall accuracy of the model was 87%, which is rather impressive, given the vast heterogeneity of host genotype, developmental stage, and environmental conditions represented by the tested ash populations.

Our results show that European ash possesses readily exploitable levels of resistance that can be detected using FT-IR spectroscopy. This work represents a major advancement in the application of marker-assisted technology for tree breeding and offers a strong proof of concept of a novel solution in the fight against ADB. Compared to labor-intensive and time-consuming traditional tree phenotyping techniques based on artificial inoculations or natural infection assays, or even nascent nucleic acid-based phenotyping⁵¹, FT-IR spectroscopy can significantly accelerate the process of selecting resistant phenotypes and limit the need for growing out large segregating progenies as in conventional breeding programs. This technique can be applicable at a landscape level for screening large naïve populations for disease resistance, at dramatically reduced costs than traditional selection. With this rapid phenotyping approach, breeding for resistance, *ex situ* and *in situ* conservation, restoration, and long-term sustainable management of this threatened tree species becomes a feasible and realistic strategy to mitigate the ADB epidemic. This method can also pave the way for rapid progress in managing the global forest health crises resulting from many other chronic and emerging forest diseases where there is evidence of some level of resistance in the naïve host populations.

## Methods

**Plant material.** Uninfected bark (including a thin outer bark (<1 mm thick), cortex, secondary phloem, and cambial tissues) and leaf samples were collected from a total of 76 *F. excelsior* genotypes with known susceptibility to *H. fraxineus* in six European countries: Austria, Denmark, France, Germany, Lithuania, and Sweden. Source material originated from genetic trials established as either clonal seed orchards³⁰⁻³²,³⁶ or for testing ash provenance³⁵ or progeny²⁴,²⁵,³³,³⁴ (Table 1). Sample collection was performed between 19 May and 18 June 2015. At each site, a minimum of three and up to seven genotypes were selected from the extreme ends of the spectrum of host susceptibility¹⁹ (low vs high), while a minimum of two and up to five genotypes of intermediate resistance were selected when available (i.e., from Denmark, France, Lithuania, and Sweden). All selections were based on a relative measurement and ranking of disease severity at the crown level (i.e. dieback intensity) as determined in previous assessments²⁴,³⁰⁻³²,³⁴⁻³⁶. In the case of clonal trials, up to three ramets per clone were sampled, for a total of 134 trees (Table 1). High heterogeneity of tree developmental stage and environmental conditions was purposely included in order to test the robustness and broad applicability of the model. The current year's shoots were harvested from each individual tree. Due to the occurrence of a late frost event in Lithuania, epicormic shoots instead of crown shoots were harvested from both resistant and susceptible genotypes. Similarly, epicormic shoots were the only ones present and suitable for harvesting from some of the susceptible genotypes across the different countries, which provided a further source of variation.

Leaves collected in the field were labelled according to country, trial, family, ramet (if pertinent) and susceptibility status, placed in a plastic bag and immediately stored on dry ice. Bark was either dissected from the internodes of the twig supporting the leaf rachis with a sterile razor blade or scalpel, and similarly labelled and stored on dry ice, or whole twigs were immediately frozen on dry ice in the field and tissues were dissected later in the lab. All samples were then transported or shipped frozen to the Swedish University of Agricultural Sciences in Alnarp, Sweden and stored at -20 °C until further processing.

**Sample processing and phenolic extraction and purification.** To avoid any potential degradation or oxidation of the tissues, samples were constantly kept frozen at -20 °C or in liquid nitrogen during all processing steps preceding phenolic extraction. Bark and leaf tissues were finely ground in liquid nitrogen and 200 ± 1 mg aliquots of either tissue type were placed in individual 2 ml microcentrifuge tubes.

Phenolic compounds were extracted and purified following the protocol described by Wrolstad *et al*.⁴⁹, with modifications. Ground tissue was submerged in 700 µl of 70% HPLC-grade acetone (Sigma-Aldrich Sweden AB, Stockholm, Sweden) in Milli-Q (Merck Chemicals and Life Science AB, Solna, Sweden) purified water (v/v). Samples were vortexed at maximum speed for 10 sec and then subjected to sonication for 30 min at room temperature, followed by centrifugation at 1,600 rcf for 8 min. The supernatant was transferred to a new 2 ml tube and twice the volume of HPLC-grade chloroform (Sigma-Aldrich Sweden AB) was added. Vials were inverted by hand two times and then vortexed at maximum speed for 15 sec prior to centrifugation at 10,000 rcf for 2 minutes at 10 °C. Working on ice, the aqueous phase, containing the phenolic compounds, was collected and transferred to a new 2 ml screw-cap tube with O-ring seal. Samples were then lyophilized and shipped at room temperature to the Department of Plant Pathology, The Ohio State University. Working on ice, lyophilized crude extracts were resuspended in 1 ml of Milli-Q (Millipore, Bedford, Massachusetts, USA) purified water and sonicated for 30 min in water and ice to allow complete solution of the pellet prior to the phenolic solid-phase purification step.
Phenolic compounds were purified on an Xpertek SPE Snap Cap 300 mg C₁₈ silica cartridge (Cobert Associates, Inc., Saint Louis, Missouri, USA). Cartridges were conditioned with HPLC-grade methanol (Fisher Scientific, Pittsburgh, Pennsylvania, USA) and equilibrated with Milli-Q water. After forcing extracts through the cartridge, cartridges were washed twice with Milli-Q water and then the adsorbed compounds were eluted in HPLC-grade methanol. Before use, extracts were evaporated to dryness in a Vacufuge vacuum concentrator (Eppendorf, Westbury, New York, USA) and re-dissolved in HPLC-grade methanol to a final concentration of 10 times the original one. Purified and concentrated phenolic extracts were stored at -20 °C until further analysis.

**Molecular detection of *Hymenoscyphus fraxineus* in leaves and bark.** The absence of *H. fraxineus* from all analysed leaves and twig bark tissues at the time of sampling was verified by PCR using species-specific ITS primers³⁸. DNA was extracted in 2 ml screw cap tubes containing 30 mg of homogenized tissue using the E.Z.N.A. SP Plant DNA Kit (Omega Bio-tek, Doraville, Georgia, USA), according to manufacturer's instructions. Amplifications were performed in 10 µl reaction volumes containing 0.05 µl of Taq Polymerase (5 u/µl); 1.0 µl of 10x Buffer; 1.0 µl of dNTPs (2 mM); 0.15 µl of MgCl₂ (25 mM); 0.2 µl of each primer *Chafra F/Chafra R*, 5.0 µl of sample DNA (0.5 ng/µl) and 2.4 µl of milli-Q water. Each PCR run included a no-template negative control of water, and DNA of *H. fraxineus* isolate nf4 collected in Sweden as positive control. The PCR cycling conditions included an initial denaturation step at 95 °C for 5 min followed by 35 amplification cycles of denaturation at 95 °C for 30 s, annealing at 56 °C for 30 s, and extension at 72 °C for 1 min. The reaction was terminated by an extension step at 72 °C for 7 min. PCR products were visualized by gel electrophoresis on a 1.5% agarose gel in TE buffer (Sigma-Aldrich Sweden AB) stained with GelRed™. The GeneRuler mix (Fermentas, Burlington, Ontario, Canada) was used as size standard.

**FT-IR analysis.** FT-IR spectroscopy of the phenolic extracts was carried out on an Excalibur 3500GX FT-IR spectrometer (Digilab, Randolph, Massachusetts, USA), equipped with a potassium bromide beamsplitter and a MIRacle triple-bounce zinc selenide crystal (Pike Technologies, Madison, Wisconsin, USA) attenuated total reflectance (ATR) accessory. Either 6 µl or 5 µl of concentrated extracts were placed on the ATR crystal and vacuum dried for approximately 1 min or 40 sec for leaves and twig bark samples, respectively. Spectra were collected in the mid-infrared region, over a wavenumber range of 4000 to 700 cm⁻¹. Resolution was set at 4 cm⁻¹ and an interferogram of 64 scans was co-added for each sample. Two technical replicates were collected for each sample, and spectra were displayed in terms of absorbance using Win-IR Pro Software (Agilent Technologies, Santa Clara, California, USA)¹⁹. Total number of spectral data collected, net of any tissue processing losses, was 262 for leaf samples and 224 for twig bark tissues, corresponding to 131 and 112 samples, respectively.

We carried out SIMCA multivariate analysis of the spectral data using the chemometric modelling software Pirouette (v. 4.5, Infometrix Inc., Bothell, Washington, USA). SIMCA classification technique develops 3D principal components models for each training group (in this case resistant, intermediate and susceptible ash trees) and identifies the most important variables for the discrimination of groups, while preserving relevant information and reducing noise⁵². Spectral data were pre-processed via the standard normal variate function to remove multiplicative scatter and particle size interference⁴⁰, and then smoothed and transformed into their second derivative using a 35-points Savitzky-Golay polynomial filter to increase the signal-to-noise ratio, minimize background and reduce overlapping bands⁴¹. Technical replicates were analysed separately¹⁹.

Preliminary analyses with the spectra of resistant and susceptible leaf samples clearly showed that clustering patterns of the SIMCA 3D class projection plot were mainly driven by the geographic location of the trees, while no pattern was associated with different resistance phenotypes (Fig. 2). Chemical composition of leaves is strongly affected by environmental factors such as nutrient and water availability⁴⁴, low temperatures⁴³ or solar irradiation⁴², and this high variability might have masked any potential chemical signal associated with resistance. We hence concluded that leaves were not ideal and were therefore excluded from further analysis. On the other hand, preliminary analyses with the spectra of resistant and susceptible twig bark samples showed that clustering patterns of the SIMCA 3D class projection plot were mainly driven by the resistance level of phenotypes, while geographic location did not show any strong patterns (Fig. 2).

Optimization of the SIMCA model was obtained by initially including the whole collected wavenumber range (4000 to 700 cm⁻¹) of all 224 spectral data for twig bark tissues, and then progressively refining the model both by reducing the spectral range to those regions with the highest discriminating power (Fig. 3c), and by removing outliers. Each incrementally refined model was evaluated by observing clustering patterns of the resistance phenotypes in the 3D class projection and Coomans plots, and by evaluating the discriminating power plot⁴⁶. Outliers were visually identified on the Coomans plot and 3D class projection plot. Based on their chemical signature, all nine Swedish and all six Lithuanian genotypes originally classified as intermediate and susceptible, based on field observations of the extent of dieback, were classified by the model as resistant and intermediate, respectively. This discrepancy in the classification, together with a general trend of overestimating resistance, may be attributable to some random variation among ramets, or most likely the use of slightly different parameters in the evaluation of susceptible phenotypes in Sweden³⁰ and Lithuania²⁵. The inclusion of uncertain phenotypes in both the training and testing sets would have compromised the strength of the chemometric model and its validation; thus intermediate and susceptible Swedish and Lithuanian genotypes were excluded from further analyses. Incremental refinement of the training models showed that, while overall clustering of the groups was not changing, the inclusion of intermediates significantly reduced the discriminating power. This is an expected outcome, as SIMCA models are known to work best when built on well-defined groups¹³,¹⁹ (in this case, the extreme ends of the spectrum of host susceptibility). Therefore, we decided to exclude intermediates from the training set and move them to the testing set instead. We opted for this approach, as we wanted to verify if the model was still able to correctly classify all intermediate samples included in the analyses. The total number of samples included in the training model, net of any trimming, was 48 (32 resistant and 16 susceptible) and 92 spectral data, while the
testing data set comprised 10 resistant, seven intermediate and six susceptible samples, for a total of 23 samples and 44 spectral data. Model accuracy was calculated as the percentage of correctly identified samples relative to the total number of samples included in the testing data set.

## Data Availability

The datasets generated and analysed during the current study are available from the corresponding authors upon reasonable request.

## References

1. Loo, J. A. Ecological impacts of non-indigenous invasive fungi as forest pathogens. *Biol. Invasions* **11**, 81–96 (2009).

2. Mack, R. N. et al. Biotic Invasions: Causes, epidemiology, global consequences, and control. *Ecol. Appl.* **10**, 689–710 (2000).

3. Boyd, I. L., Freer-Smith, P. H., Gilligan, C. A. & Godfray, H. C. J. The consequence of tree pests and diseases for ecosystem services. *Science* **342**, 1235773 (2013).

4. Santini, A. et al. Biogeographical patterns and determinants of invasion by forest pathogens in Europe. *New Phytol.* **197**, 238–250 (2013).

5. Prospero, S. & Cleary, M. Effects of host variability on the spread of invasive forest diseases. *Forests* **8**, 80 (2017).

6. Sturrock, R. N. et al. Climate change and forest diseases. *Plant Pathol.* **60**, 133–149 (2011).

7. Desprez-Louostau, M.-L., Marçais, B., Nageleisen, L.-M., Piou, D. & Vannini, A. Interactive effects of drought and pathogens in forest trees. *Ann. For. Sci.* **63**, 597–612 (2006).

8. Anagnostakis, S. L. The effect of multiple importations of pests and pathogens on a native tree. *Biol. Invasions* **3**, 245–254 (2001).

9. Donovan, G. H. et al. The relationship between trees and human health: evidence from the spread of the emerald ash borer. *Am. J. Prev. Med.* **44**, 139–145 (2013).

10. Pecori, F., Ghelardini, L., Luchi, N., Pepori, A. & Santini, A. Lights and shadows of a possible strategy to cope with alien and destructive forest pathogens: The example of breeding for resistance to Dutch elm disease in Italy. *Baltic For.* **23**, 255–263 (2017).

11. Kovalchuk, A. et al. Antimicrobial defenses and resistance in forest trees: challenges and perspectives in a genomic era. *Annu. Rev. Phytopathol.* **51**, 221–244 (2013).

12. Harfouche, A. et al. Accelerating the domestication of forest trees in a changing world. *Trends Plant Sci.* **17**, 64–72 (2012).

13. Conrad, A. O. & Bonello, P. Application of infrared and raman spectroscopy for the identification of disease resistant trees. *Front. Plant Sci.* **6**, 1152 (2015).

14. Ellis, D. I. & Goodacre, R. Metabolic fingerprinting in disease diagnosis: biomedical applications of infrared and Raman spectroscopy. *Analyst* **131**, 875–885 (2006).

15. Altangerel, N. et al. *In vivo* diagnostics of early abiotic plant stress response via Raman spectroscopy. *Proc. Natl. Acad. Sci.* **114**, 3393–3396 (2017).

16. Pichersky, E. & Lewinsohn, E. Convergent evolution in plant specialized metabolism. *Ann. Rev. Plant Biol.* **62**, 549–566 (2011).

17. Martín, J., Solla, A., Woodward, S. & Gil, L. Fourier transform-infrared spectroscopy as a new method for evaluating host resistance in the Dutch elm disease complex. *Tree Physiol.* **25**, 1331–1338 (2005).

18. Martin, J. A., Solla, A., Coimbra, M. A. & Gil, L. Metabolic fingerprinting allows discrimination between *Ulmus pumila* and *U. minor*, and between *U. minor* clones of different susceptibility to Dutch elm disease. *Forest Pathol.* **38**, 244–256 (2008).

19. Conrad, A. O., Rodriguez-Saona, L. E., McPherson, B. A., Wood, D. L. & Bonello, P. Identification of *Quercus agrifolia* (coast live oak) resistant to the invasive pathogen *Phytophthora ramorum* in native stands using Fourier-transform infrared (FT-IR) spectroscopy. *Front. Plant Sci.* **5**, 521 (2014).

20. European Commission. Interpretation Manual of European Union habitats EUR27 DG Environment, Nature and Biodiversity. (Bruxelles, Belgium, 2007).

21. Baral, H.-O. & Bemmann, M. *Hymenoscyphus fraxineus* vs. *Hymenoscyphus albidus* – A comparative light microscopic study on the causal agent of European ash dieback and related foliicolous, stroma-forming species. *Mycology* **5**, 228–290 (2014).

22. Gross, A., Holdenrieder, O., Pautasso, M., Queloz, V. & Sieber, T. N. *Hymenoscyphus pseudoalbidus*, the causal agent of European ash dieback. *Mol. Plant Pathol.* **15**, 5–21 (2014).

23. Kjær, E. D., McKinney, L. V., Nielsen, L. R., Hansen, L. N. & Hansen, J. K. Adaptive potential of ash (*Fraxinus excelsior*) populations against the novel emerging pathogen *Hymenoscyphus pseudoalbidus*. *Evol. Appl.* **5**, 219–228 (2012).

24. Muñoz, F., Marçais, B., Dufour, J. & Dowkiw, A. Rising out of the ashes: additive genetic variation for crown and collar resistance to *Hymenoscyphus fraxineus* in *Fraxinus excelsior*. *Phytopathology* **106**, 1535–1543 (2016).

25. Pliūra, A., Lygis, V. & Suchockas, V. Performance of twenty four European *Fraxinus excelsior* populations in three Lithuanian progeny trials with a special emphasis on resistance to *Chalara fraxinea*. *Baltic For.* **17**, 17–34 (2011).

26. Downie, J. A. Ash dieback epidemic in Europe: How can molecular technologies help? *PLOS Pathog.* **13**, e1006381 (2017).

27. McKinney, L. V. et al. The ash dieback crisis: genetic variation in resistance can prove a long-term solution. *Plant Pathol.* **63**, 485–499 (2014).

28. Sollars, E. S. A. et al. Genome sequence and genetic diversity of European ash trees. *Nature* **541**, 212–216 (2017).

29. Cozzolino, D. Use of infrared spectroscopy for in-field measurement and phenotyping of plant properties: instrumentation, data analysis, and examples. *App. Spectrosc. Rev.* **49**, 564–584 (2014).

30. Stener, L.-G. Clonal differences in susceptibility to the dieback of *Fraxinus excelsior* in southern Sweden. *Scand. J. Forest Res.* **28**, 205–216 (2013).

31. Kirisits, T. & Freinschlag, C. Ash dieback caused by *Hymenoscyphus pseudoalbidus* in a seed plantation of *Fraxinus excelsior* in Austria. *J. Agr. Extension Rural Dev.* **4**, 184–191 (2012).

32. McKinney, L. V., Nielsen, L. R., Hansen, J. K. & Kjær, E. D. Presence of natural genetic resistance in *Fraxinus excelsior* (Oleraceae) to *Chalara fraxinea* (Ascomycota): an emerging infectious disease. *Heredity* **106**, 788–797 (2011).

33. Pliūra, A. & Baliuckas, V. Genetic variation in adaptive traits of progenies of Lithuanian and western European populations of *Fraxinus excelsior* L. *Baltic For.* **13**, 28–38 (2007).

34. Pliūra, A., Marčiulynienė, D., Bakys, R. & Suchockas, V. Dynamics of genetic resistance to *Hymenoscyphus pseudoalbidus* in juvenile *Fraxinus excelsior* clones. *Baltic For.* **20**, 10–27 (2014).

35. Enderle, R., Peters, F., Nakou, A. & Metzler, B. Temporal development of ash dieback symptoms and spatial distribution of collar rots in a provenance trial of *Fraxinus excelsior*. *Eur. J. For. Res.* **132**, 865–876 (2013).

36. Heinze, B., Tiefenbacher, H., Litschauer, R. & Kirisits, T. Ash dieback in Austria – history, current situation and outlook in *Dieback of European Ash (Fraxinus spp.) – Consequences and Guidelines for Sustainable Management* (eds Vasaitis, R. & Enderle, R.) 33–52 (SLU Service/Repro, 2017).

37. Lattanzio, V., Lattanzio, V. M. T. & Cardinali, A. Role of phenolics in the resistance mechanisms of plants against fungal pathogens and insects in *Phytochemistry Advances in Research* (ed. Imperato, F.) 23–67 (Research Signpost, 2006).

38. Johansson, S. B. K., Vasaitis, R., Ihrmark, K., Barklund, P. & Stenlid, J. Detection of *Chalara fraxinea* from tissue of *Fraxinus excelsior* using species-specific ITS primers. *Forest Pathol.* **40**, 111–115 (2010).

39. Ockels, F. S., Eyles, A., McPherson, B. A., Wood, D. L. & Bonello, P. Phenolic chemistry of coast live oak response to *Phytophthora ramorum* infection. *J. Chem. Ecol.* **33**, 1721–1732 (2007).
40. Barnes, R. J., Dhanoa, M. S. & Lister, S. J. Standard normal variate transformation and de-trending of near-infrared diffuse reflectance spectra. *Appl. Spectrosc.* **43**, 772–777 (1989).

41. Savitzky, A. & Golay, M. J. Smoothing and differentiation of data by simplified least squares procedures. *Anal. Chem.* **36**, 1627–1639 (1964).

42. Sullivan, J. H., Gitz, D. C., Peek, M. S. & McElrone, A. J. Response of three eastern tree species to supplemental UV-B radiation: leaf chemistry and gas exchange. *Agr. For. Entomol.* **120**, 219–228 (2003).

43. St Clair, S. B., Monson, S. D., Smith, E. A., Cahill, D. G. & Calder, W. J. Altered leaf morphology, leaf resource dilution and defense chemistry induction in frost-defoliated aspen (*Populus tremuloides*). *Tree Physiol.* **29**, 1259–1268 (2009).

44. Cunningham, S. A., Summerhayes, B. & Westoby, M. Evolutionary divergences in leaf structure and chemistry, comparing rainfall and soil nutrient gradients. *Ecol. Monogr.* **69**, 569–588 (1999).

45. Cleary, M. R. et al. Genotypes of *Fraxinus excelsior* with different susceptibility to the ash dieback pathogen *Hymenoscyphus pseudoalbidus* and their response to the phytotoxin viridiol – A metabolomic and microscopic study. *Phytochemistry* **102**, 115–125 (2014).

46. Coomans, D. & Broeckaert, I. *Potential pattern recognition in chemical and medical decision making*. (John Wiley & Sons, Inc., 1986).

47. Larkin, P. *Infrared and Raman spectroscopy: principles and spectral interpretation*. (Elsevier, 2011).

48. Müller, U., Rätzsch, M., Schwanninger, M., Steiner, M. & Zöbl, H. Yellowing and IR-changes of spruce wood as result of UV-irradiation. *J. Photoch. Photobio. B* **69**, 97–105 (2003).

49. Wrolstad, R. E. et al. *Handbook of food analytical chemistry: Pigments, colorants, flavors, texture, and bioactive food components*. (Wiley Interscience, 2005).

50. Lupoi, J. S. et al. High-throughput prediction of eucalypt lignin syringyl/guaiacyl content using multivariate analysis: a comparison between mid-infrared, near-infrared, and Raman spectroscopies for model development. *Biotechnol. Biofuels* **7**, 93 (2014).

51. Harper, A. L. et al. Molecular markers for tolerance of European ash (*Fraxinus excelsior*) to dieback disease identified using Associative Transcriptomics. *Sci. Rep.* **6**, 19335 (2016).

52. Mark, H. Data analysis: Multilinear regression and principal component analysis in *Handbook of Near-Infrared Analysis* (eds Burns, D. A. & Ciurczak, E. W.) 129–184 (Taylor and Francis, 2001).

## Acknowledgements

We thank C. Boldrini, V. Guérin, F. Oskay and A. Stradner for assistance with field sampling and logistics, and H. Konrad and C. Jasser for permission to sample the Austrian seed orchard. Research was funded by Skogssällskapet (The Swedish Forest Society Foundation) (Grant no. 2016-015); Kungl. Skogs- och Lantbruksakademien (The Royal Swedish Academy of Agriculture and Forestry) (Grant no. GFS2015-0189); EU Cost Action FP1103 FRAXBACK; The Danish Council for Independent Research (Grant no. 6111-00254B); Villum Foundation (Grant no. VKR 023062); French Ministry of Agriculture (Programme 149 – Action 13 – Sous action 32); EU COST Action TD1209 ALIEN CHALLENGE; Lithuanian Science Council (National Research Program project no. SIT-4/2015 MIŠKOEKOKAITA); ForstBW; The Austrian Federal Ministry of Agriculture, Forestry, Environment and Water Management (BMLFUW research project no. 100343); The Provincial Governments of Lower Austria, Styria, Upper Austria, Carinthia, Salzburg and Burgenland; The Forest Office and Urban Agriculture (MA 49) of the Vienna City Administration; The Austrian Federal Forests (ÖBf AG); The Foundation “120 Jahre Universität für Bodenkultur”; and state and federal funds appropriated to the Ohio Agricultural Research and Development Center, The Ohio State University.

## Author Contributions

P.B. ideated the experiment. C.V., P.B. and M.C. designed the experiment. A.D., R.E., T.K., E.K., L.M., B.M., F.M., L.N., A.P., L.-G.S., V.S. and M.C. collected the samples on trees that they had categorized from previous studies. C.V., M.G. and M.C. processed the samples and extracted phenolics. D.M. performed molecular analyses. C.V. performed FT-IR analyses and chemometric modelling. L.R.-S. helped with the chemometric modelling. C.V., P.B. and M.C. discussed results and wrote the paper, and all other authors contributed equally to the final version of the manuscript.

## Additional Information

**Competing Interests:** The authors declare no competing interests.

**Publisher’s note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_8.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2018