# Modification of Early Response of Vitis vinifera to Pathogens Relating to Esca Disease and Biocontrol Agent Vintec® Revealed By Untargeted Metabolomics on Woody Tissues

Justine Chervin, Ana Romeo-Oliván, Sylvie Fournier, Virginie Puech Pagès, Bernard Dumas, Alban Jacques, Guillaume Marti

► To cite this version:

Justine Chervin, Ana Romeo-Oliván, Sylvie Fournier, Virginie Puech Pagès, Bernard Dumas, et al. Modification of Early Response of Vitis vinifera to Pathogens Relating to Esca Disease and Biocontrol Agent Vintec® Revealed By Untargeted Metabolomics on Woody Tissues. Frontiers in Microbiology, 2022, 13, pp.835463. 10.3389/fmicb.2022.835463. hal-03652523

HAL Id: hal-03652523

https://ut3-toulouseinp.hal.science/hal-03652523

Submitted on 26 Apr 2022

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
<figure><img src="image_2.png" /></figure>

Check for
updates

# Modification of Early Response of Vitis vinifera to Pathogens Relating to Esca Disease and Biocontrol Agent Vintec® Revealed By Untargeted Metabolomics on Woody Tissues

Justine Chervin<fnref n="1" />,<fnref n="2" />,<fnref n="3" />, Ana Romeo-Oliván⁴, Sylvie Fournier<fnref n="1" />,<fnref n="2" />,<fnref n="3" />, Virginie Puech-Pages<fnref n="1" />,<fnref n="2" />,<fnref n="3" />, Bernard Dumas<fnref n="1" />, Alban Jacques⁴ and Guillaume Marti<fnref n="1" />,<fnref n="2" />,<fnref n="3" />*

OPEN ACCESS

**Edited by:**

Florence Fontaine,
Université de Reims
Champagne-Ardenne, France

**Reviewed by:**

Marta R. M. Lima,
Virginia Tech, United States
Julie Chong,
University of Upper Alsace, France

**\*Correspondence:**

Guillaume Marti
guillaume.marti@univ-tlse3.fr

**Specialty section:**

This article was submitted to
Microbe and Virus Interactions with
Plants,
a section of the journal
Frontiers in Microbiology

Received: 14 December 2021

Accepted: 19 January 2022

Published: 02 March 2022

**Citation:**

Chervin J, Romeo-Oliván A,
Fournier S, Puech-Pages V, Dumas B,
Jacques A and Marti G (2022)
Modification of Early Response of Vitis
vinifera to Pathogens Relating to Esca
Disease and Biocontrol Agent Vintec®
Revealed By Untargeted
Metabolomics on Woody Tissues.
Front. Microbiol. 13:835463.
doi: 10.3389/fmicb.2022.835463

Esca disease is one of the most destructive grapevine trunk diseases. *Phaeoacremonium minimum* and *Phaeomoniella chlamydospora* are two of the known fungal pathogens associated with this disease. Today, biocontrol agents against Esca are mainly based on the use of the strain of the mycoparasite fungal genus *Trichoderma* such as the Vintec® product. The aim of this study was to investigate early response of woody tissues to Esca pathogens and identify metabolites that could be correlated with a biocontrol activity within a complex woody matrix. An untargeted liquid chromatography–high-resolution mass spectrometry metabolomic approach coupled to a spectral similarity network was used to highlight clusters of compounds associated with the plant response to pathogens and biocontrol. Dereplication highlighted the possible role of glycerophospholipids and polyphenol compounds, the latest mainly belonging to stilbenoids. Antifungal activity of some relevant biomarkers, evaluated *in vitro* on *Phaeomoniella chlamydospora* and *Botrytis cinerea*, suggests that some of these compounds can play a role to limit the development of Esca pathogens in planta.

Keywords: Vitis vinifera, Esca disease, metabolomics, biomarkers, biocontrol

INTRODUCTION

Grapevine trunk diseases (GTDs) are of growing concern among worldwide viticulture. These diseases concern different organs of the plant, some causing the progressive death of vines in vineyards. These losses influence the productivity of the vineyard and, therefore, represent a considerable economic impact in the viticultural sector (Hofstetter et al., 2012; Fontaine et al., 2016). Grapevine trunk diseases are associated with the presence of different pathogenic fungal species, which can affect vines at different stages of their life cycle (Bertsch et al., 2013; Bruez et al., 2014).

<footnote n="1">¹ Laboratoire de Recherche en Sciences Végétales, Université de Toulouse, CNRS, UPS, Toulouse INP, Toulouse, France,</footnote>

<footnote n="2">² Metatoul-AgromiX Platform, LRSV, Université de Toulouse, CNRS, UPS, Toulouse INP, Toulouse, France,</footnote>

<footnote n="3">³ MetaboHUB-MetaToul, National Infrastructure of Metabolomics and Fluxomics, Toulouse, France, ⁴ Unité de Recherche Physiologie, Pathologie, et Génétique Végétales (PPGV), INP PURPAN, Université de Toulouse, Toulouse, France</footnote>
Esca is one of the most complex and destructive trunk diseases. The Esca disease comprises different syndromes, depending on the symptoms, the age of the plant, and the fungal species associated (Graniti et al., 2000; Surico, 2009; Bertsch et al., 2013). One of the syndromes, named Esca, is associated with white rot in the trunk, and it is caused by fungal species of the genus *Formitiporia* (Bruno and Sparapano, 2007). The three vascular syndromes—brown wood streaking (Mugnai et al., 1999), Petri disease, and grapevine leaf stripe disease (GLSD) (Bertsch et al., 2013)—are associated with *Phaeomoniella chlamydospora* (*P. chlamydospora*), *Phaeoacremonium minimum* (*P. minimum*), and other *Phaeoacremonium* ssp. (Luque et al., 2009; Úrbez-Torres et al., 2014). The presence of these pathogens in the wood affects not only the trunk but also leaves and berries. Symptoms on leaves start with chlorotic areas becoming necrotic and conferring a "tiger striped" appearance. On berries, dark purple spots on the epidermis can be observed (Mugnai et al., 1999). The last syndrome, Esca proper, corresponds to the co-occurrence of Esca and GLSD in the same plant (Surico, 2009).

In the absence of curative methods, since the banishment of sodium arsenite, there are numerous strategies to limit the occurrence of GTDs both in nurseries and in the field (Mondello et al., 2018a,b). One of these strategies includes the use of biocontrol agents (BCAs). Biocontrol agents represent promising alternatives to conventional phytosanitary products (Ram et al., 2018). *Trichoderma* spp. have been largely studied and are of great interest as biocontrol agents (Moya et al., 2020; Poveda, 2021). In viticulture, there exist two *Trichoderma*-based commercial solutions to control GTDs: Vintec® (Belchim Crop Protection) and Esquive® (Agrauxine by Lesaffre). Some microorganisms have also achieved good results in grapevine protection against GTDs, such as the oomycete *Pythium oligandrum* (Yacoub et al., 2016) or a number of bacteria (Trotel-Aziz et al., 2019; Haidar et al., 2020). *Pythium oligandrum* is known to have an antifungal effect against different plant fungal pathogens (Gerbore et al., 2014). Concerning GTDs, the presence of this oomycete reduces *P. chlamydospora* infection (Yacoub et al., 2016). The known modes of action of *P. oligandrum* are diverse and, as well as *Trichoderma* spp., it can enhance plant defense responses (Yacoub et al., 2020). In the case of bacteria, Trotel-Aziz et al. (2019) found that *Bacillus subtilis* protects grapevine against *Neofusicoccum parvum*, the causal agent of Botryosphaeria dieback, by the activation of detoxification process in the plant. In addition, Haidar et al. (2021) found that bacteria from the genus *Paenibacillus* produced antifungal metabolites and volatile organic compounds that reduced the infection of *N. parvum* in the grapevine's trunk.

*Trichoderma* species have been widely studied for this purpose. Their ability to protect plants resides in the combination of several factors: (i) the enhancement of nutrient uptake and mobilization, (ii) mycoparasitism and competition with the pathogens, and (iii) induction of plant systemic resistance (Vinale et al., 2008; Braun et al., 2018). Vintec®, a formulation of *Trichoderma atroviridae* SC1 (TASC1), is one of the main approved BCAs today to fight trunk diseases (Esca, black dead arm, eutypiose). Its efficiency has been successfully tested in both nursery and vineyard (Berbegal et al., 2020).

In order to better understand this disease, different studies were focused on the physiological modifications on Esca-affected vines (Fontaine et al., 2016). Early studies reported the accumulation of stress metabolites in the wood of Esca-diseased vines (Amalfitano et al., 2000, 2011). Modifications in phenolic compounds in roots showing an early decline were also observed (Del Rio et al., 2001). However, most available works are based on targeted approaches linked with polyphenol compound changes as described in Amalfitano et al. (2011); Lima et al. (2012), or Rusjan et al. (2017), for example.

Untargeted approaches as transcriptomics and metabolomics are of prime interest to study the overall modifications induced in plants during disease outcomes. Recently, several omic studies were assessed in the context of Botryosphaeria dieback issue. A metabolomic study based on wood response with the pathosystem *Neofusicoccum parvum* revealed significant modifications of primary and secondary metabolites after inoculation, in particular a decrease in sugars and an increase in stilbenes and flavonoids (Labois et al., 2020). Another study also reported the importance of lipids in understanding this disease and showed that in diseased vines, phytochemicals, in particular stilbenes and flavonoids, mainly characterized the brown stripe area, whereas lipids, in particular phospholipids, mainly characterized the adjacent white area (Lemaitre-Guillier et al., 2020).

Regarding the Esca issue, a first report on the pathosystem *Vitis vinifera* / Esca upon untargeted metabolomics through nuclear magnetic resonance spectroscopy-based metabolomic showed that a double stress induced by fungi and drought led to the change of several amino acids (Lima et al., 2017b). Nowadays, the majority of studies focus on leaf metabolomes such as the investigation of phenolic modification in leaves of *Vitis vinifera* cv. Alvarinho revealing an apparent accumulation of phenolic compounds in diseased leaves (Lima et al., 2017a). Another work focusing on alterations in leaf metabolome reported that changes, especially concerning phenylpropanoids but also carbohydrates and amino acids, could already be observed prior to symptom appearance (Magnin-Robert et al., 2017). Lipidomic analysis of leaves of Esca-affected grapevine was also conducted (Goufo and Cortez, 2020) and notably showed a general increase in most lipids that correlated with symptom progression.

The present study investigates the response of grapevine woody tissue to fungal infection with the biocontrol agent (BCA) Vintec® (IV), a mixture of both fungi *P. chlamydospora* and *P. minimum* (IPP), as it is usually observed in a natural context, and their combination (IVPP), at an early stage of infection by an ultra-high-performance liquid chromatography-high-resolution mass spectrometry (UHPLC-HRMS)-based metabolomic approach.

The aim of this UHPLC-HRMS approach is to highlight metabolite changes that could be further considered as early biomarkers of infection. This is crucial since fungi associated with Esca disease are known to colonize grapevines during several years without any detectable symptoms in aerial organs like leaves (Bruez et al., 2016).

For this study, a pathosystem was set up to mimic the natural events by co-inoculation of pathogens with and without a
biocontrol agent in the trunk of canes. The aim is to be as close as
possible to the vineyard to propose biomarkers to monitor both
Esca pre-emergence and biocontrol effects.

# MATERIALS AND METHODS

Standards

Individual solutions of natural products (NPs) compounds from Metasci (Metasci, Toronto, Canada) and Chemfaces (Wuhan, Hubei) were prepared at 100 µg/ml according to suppliers' recommendations. Belchim Crop Protection (Lissieu, France) supplied the BCA Vintec®.

Fungal Material

*Phaeoacremonium minimum* (CBS 100398), *Phaeomoniella chlamydospora* (CBS 239.74), and *Botrytis cinerea* (CBS 131.28) strains came from the fungal biodiversity center (Centraalbureau voor Schimmelcultures, CBS. Utrecht, the Netherlands). Fungal isolates were grown in individual Petri plates containing the malt extract-agar (MEA) medium for five weeks at 26°C in obscurity before use.

Plant Culture

One-year-old canes of *Vitis vinifera* L. cv Cabernet-Sauvignon clone 15 purchased from a French nursery (Pépinières Daydé, Montans, Midi-Pyrénées, France) of the 2019 and 2021 seasons were divided in cuttings containing two dormant nodes. Cuttings were then kept at 5°C for one night. Before planting, their surface was disinfected. At first, cuttings were soaked for 30 s in a 10-L water bath with 5-ml bleach (2.6% active chlorine) and then rinsed twice with 10 L of tap water. Then, cuttings were stored overnight at 4°C in a 10-L water bath containing 5-ml 8-hydroxyquinoline sulfate (0.05% Cryptonol®, Chauvin, France), rinsed three times with 10 L of demineralized water. They were planted in moist autoclaved (121°C, 15 min) rock wool to allow disbudding and rooting. At this point, plants were watered twice a week. Four weeks later, cuttings were transferred from the rock wool into individual pots containing 100 g of autoclaved soil (PAM2, PROVEEN Substrates). Plants were finally grown in a plant growth chamber at 25°C, with 45% of humidity, a photoperiod of 12 h/12 h, and watered three times weekly.

Plant Inoculation

The following conditions were studied, and for each, 16 plants
were used: Injured/Not inoculated (INi); Injured/Vintec®
(IV); Injured/P. minimum + P. chlamydospora (IPP);
Injured/Vintec® + P. minimum + P. chlamydospora (IVPP).

The injured plants were drilled at the internode, as described
in Reis et al. (2019). Plants injured but not inoculated with
Vintec® received 20 µl in the hole of sterile distilled water; plants
injured and inoculated with the biocontrol received 20 µl in
the hole of a suspension of the product Vintec® at 2 g L⁻¹ of
formulated product 5 days before infection with Esca fungi (T-
5d). The injury was covered with a strip of Parafilm (American
National Can, Chicago, IL).

Inoculation with fungi *P. minimum* and *P. chlamydospora* at
T0 was performed by drilling again all the plants at the same

spot, and a plug of agar-agar colonized by fungi was directly
introduced. All the plants were grown in a plant growth chamber
at 25°C, with 45% of humidity and a photoperiod of 12 h/12 h
for 3 or 6 weeks.

Wood Sample Extraction for
Metabolomic Analysis and Fungal
Colonization Quantification

For metabolomic analysis, samples were harvested at 3 weeks (T3) and 6 weeks (T6) after treatment. We performed an analysis on samples collected above and behind the inoculation site, and these preliminary data did not show a systemic response. This observation agreed with Pierron et al. (2015), showing that wood responses to *P. minimum* (*P. aleophilum*) were mainly restricted to the inoculation point six weeks post-infection. Thus, around 1 cm of wood in the inoculation part (I) was collected and immediately dipped in liquid nitrogen. Likewise, wood samples of the same treatment were pooled in 6 tubes containing 3 pieces of wood from 3 different plants. Each sample was then grounded with a Mixer Mill MM 300 (Retsh, Eragny sur Oise, France) by applying 28 oscillations per second for 1 min. Around 100 mg of powder was then placed in a FastPrep tube (MP Biomedicals Lysing Matrix D, Illkirch, France) and kept at -80°C.

For sample extraction, 1 ml of 70% EtOH was added per
100 mg followed by three cycles of 20 s at 6 m/s in the
FastPrep while keeping samples in ice between each cycle. After a
centrifugation at 4°C and 12,000 rpm for 10 min, the supernatant
was collected and a second extraction was achieved on the
residue. The two supernatants were pooled, transferred to vials,
and kept at -20°C before injection. An extraction blank (without
plant material) and quality control (QC) samples were also
prepared for extraction and analytical validation.

A set of three plants of each treatment was reserved for
evaluating fungal colonization six weeks post-inoculation via
qPCR quantifications. DNA extraction was performed using
a CTAB/CIA-based protocol combined with a Qiagen DNA
extraction kit (DNeasy plant mini Kit, Qiagen, United States),
as described in *Romeo-Oliván et al.* (2021). qPCR reactions
were conducted with the GoTaq® qPCR System (Promega) and
ABI 7500 Real-Time PCR Cycler (Applied Biosystems, Foster
City, United States) following this cycling program: 15 min
at 95°C (denaturation), 40 cycles of 15 s at 95°C and 45 s
at 62°C (annealing and extension), and 40 min from 60 to
95°C (melting analysis). Specific primers targeting the β-tubulin
gene of *P. chlamydospora* and *P. minimum* were previously
described (Pouzoulet *et al.*, 2013) as well as those amplifying
the endochitinase 42 gene of *T. atroviridae* SC1 for Vintec®
(Savazzini *et al.*, 2008). They were used at 0.5 µM in a
final reaction volume of 10 µl. The number of gene copies
of *P. chlamydospora* and *P. minimum* was estimated using a
standard curve (Pouzoulet *et al.*, 2013). Two-way ANOVA and
Dunnett's *post hoc* test (p ≤ 0.05) were used to compare the gene
copy number of the various modalities using GraphPad Prism 8
version 8.3.0 (San Diego, California).

For Vintec®, we only considered the presence/absence
of amplification.
# Ultra-High-Performance Liquid Chromatography-High-Resolution Mass Spectrometry Profiling

Ultra-high-performance liquid chromatography-high-resolution MS (UHPLC-HRMS) analyses were performed on a Q Exactive Plus quadrupole (Orbitrap) mass spectrometer, equipped with a heated electrospray probe (HESI II) coupled to a U-HPLC Ultimate 3000 RSLC system (Thermo Fisher Scientific, Hemel Hempstead, United Kingdom). Separation was done on a Luna Omega Polar C18 column (150 mm × 2.1 mm i.d., 1.6 µm, Phenomenex, Sartrouville, France) equipped with a guard column. The mobile phase A (MPA) was water with 0.05% formic acid (FA), and the mobile phase B (MPB) was acetonitrile with 0.05% FA. The solvent gradient was 0 min, 100% MPA; 1 min, 100% MPA; 22 min, 100% MPB; 25 min, 100% MPB; 25.5 min, 100% MPA; and 28 min, 100% MPA. The flow rate was 0.3 ml/min, the column temperature was set to 40°C, the autosampler temperature was set to 5°C, and the injection volume was fixed to 5 µl. Mass detection was performed in positive ionization (PI) mode at resolution 35,000 power [full width at half-maximum (fwhm) at 400 m/z] for MS1 and 17,500 for MS2 with an automatic gain control (AGC) target of $1 \times 10^6$ for full scan MS1 and $1 \times 10^5$ for MS2. Ionization spray voltages were set to 3.5 kV, and the capillary temperature was kept at 256°C. The mass scanning range was m/z 100-1500. Each full MS scan was followed by data-dependent acquisition of MS/MS spectra for the six most intense ions using stepped normalized collision energy of 20, 40, and 60 eV.

# Data Processing

Ultra-high-performance liquid chromatography-high-resolution mass spectrometry (UHPLC-HRMS) raw data were processed with MS-DIAL version 4.70 (Tsugawa et al., 2015) for mass signal extraction between 100 and 1,500 Da from 0.5 to 18.5 min. Respectively, MS1 and MS2 tolerance were set to 0.01 and 0.05 Da in the centroid mode. The optimized detection threshold was set to $1.5 \times 10^6$ concerning MS1 and 10 for MS2. Peaks were aligned on a QC reference file with a retention time tolerance of 0.1 min and a mass tolerance of 0.015 Da. Peak annotation was performed with an in-house database built on an MS-FINDER model (Tsugawa et al., 2016).

MS-DIAL data were then cleaned with the MS-CleanR workflow (Fraisier-Vannier et al., 2020) by selecting all filters with a minimum blank ratio set to 0.8, a maximum relative standard deviation (RSD) set to 40, and a relative mass defect (RMD) ranging from 50 to 3,000. The maximum mass difference for feature relationships detection was set to 0.005 Da and the maximum RT difference to 0.025 min. The Pearson correlation links were considered with correlation $\ge$ 0.8 and statistically significant with $\alpha$ = 0.05. Two peaks were kept in each cluster with the most intense and the most connected. The kept features (m/z × RT pairs) were annotated with MS-FINDER version 3.52. The MS1 and MS2 tolerances were, respectively, set to 5 and 15 ppm. Formula finders were only processed with C, H, O, N, P, and S atoms. Databases (DBs) based on *Vitis* (genus), *Vitaceae* (family), *Trichoderma* (genus of the BCA strain), and *Togniniaceae* (the family of the two pathogenic fungi) were

constituted with the dictionary of natural product (DNP, CRC press, DNP on DVD v. 28.2). The internal generic DBs from MS-FINDER used were KNApSAcK, PlantCyc, HMDB, LipidMaps, NANPDB, and UNPD. Annotation prioritization was done by ranking *Vitis* DB, followed by *Vitaceae* DB, *Trichoderma* and *Togniniaceae* DBs, and finally generic DBs using the final MS-CleanR step.

# Statistical Analysis

Statistical analyses were made with Orange 3.30.1 (Demsar et al., 2013). All data were scaled by unit variance scaling (UV) before multivariate analysis. ANOVA models were proposed to rank features and highlight biomarkers of plant response to the various studied treatments. Models confronting only two treatments at once were built taking as reference the sample Injured/Not inoculated (INi). A random forest model was used to validate the comparisons based on the cross-validation algorithm. Two-way ANOVA and Dunnett's *post hoc* test ($p \le 0.05$) were used to check biomarker significance for both kinetics in comparison to references INi with GraphPad Prism 8 version 8.3.0 (San Diego, California).

The heatmap of revealed biomarkers based on average peak area was built with the web-interface MetaboAnalyst version 4.0 (Chong et al., 2019).

# Mass Spectral Similarity Network

The.msp PI and metadata files generated at the end of the MS-CleanR workflow were imported into MetGem version 1.2.2. A mass spectral similarity network was built with a cosine score cutoff fixed at 0.6, a maximum of ten connections between nodes, and at least four matched peaks. The resulting network was imported into Cytoscape12 version 3.7.2 to tune visualization. Node color was based on chemical classes, and the size of revealed biomarkers was increased. The edge width was deepened according to the cosine score.

# In vitro Bioassays

Antifungal activity of pure available standard metabolites was investigated *in vitro* in 96-well plate cultures by measuring fungal growth. Standard solutions were previously dissolved in 8% DMSO in water with concentration ranging from 0.1 to 1.6 mg/ml. Three final concentrations in each well were investigated for each standard, based on antifungal activity results found in the literature: 100, 10, and 1 µg/ml for $\alpha$-viniferin (Langcake and Pryce, 1977); 50, 5, and 0.5 µg/ml for viniferin (Langcake, 1981), and finally 25, 2.5, and 0.25 µg/ml for pterostilbene (Jeandet et al., 2002). As no information was found on oxyresveratrol, a higher range of concentration was studied: 400, 40, and 4 µg/ml.

To obtain *Botrytis cinerea* and *P. chlamydospora* spores, two plugs of agar colonized by the mycelium of each were cut, placed in 1-ml sterile water, and vortexed. The spore concentration was adjusted under a microscope using a Malassez counting chamber to $10^3$ spores/ml with water.

Triplicates using *Botrytis cinerea* spores as the positive control as well as *P. chlamydospora* spores as the target pathogen were prepared for each standard concentration. To the well, 130 µl of malt extract medium (0.01 g/ml) was added, followed by 20 µl
of the spore suspension and 50 µl of standard solutions. For the negative control, 20 µl of the spore suspension was replaced by 20 µl of sterile distilled water. Plates were then incubated at 26°C for three days. The absorbance was measured at 600 nm in a microplate reader (Infinite 200 Pro Nano, Tecan, Austria). IC₅₀ (concentrations causing an inhibition of 50% of fungal growth) were then calculated using GraphPad Prism 8 version 8.3.0 (San Diego, California). Control replicates were prepared to subtract the compound absorbance.

# RESULTS

## Fungal Inoculation on Vitis vinifera

The dataset was built on a reductionist model to study the local effect of the BCA Vintec® in the metabolic regulation of a plant challenged or not with the two fungal pathogens and highlight early biomarkers of pathogens and biocontrol presence (Figure 1).

*Phaeomoniella chlamydospora* and *P. minimum* colonization was estimated via quantifications by qPCR of their β-tubulin gene (Figure 2). The results indicated that the wood was still colonized six weeks post-infection by the two fungal pathogens in IPP and IVPP samples. Moreover, a significant reduction of colonization by *P. chlamydospora* and *P. minimum* due to the BCA Vintec® was observed. As expected, control samples (INi) and samples inoculated with Vintec® (IV) were free of pathogenic DNA. The presence of *T. atroviridae* SC1 for Vintec® pretreatment was also assessed via qPCR, confirming the colonization by *T. atroviridae* in the biocontrol-treated IV and IVPP samples (Supplementary Table 1). Altogether, these results corroborated that the data obtained further are related to the presence of either the pathogens, the BCA, or the combination of the three.

## Metabolomic Analysis of Wood Samples

Ultra-high-performance liquid chromatography–high-resolution mass spectrometry (UHPLC-HRMS) profiles of all the extracts afforded 239 features (m/z-RT pairs) in the PI mode after

<figure><img src="image_3.png" /><figcaption>**FIGURE 1** | Grapevine model used to study Esca disease in laboratory conditions: (A) drilling injury site, (B) inoculation of biocontrol agent (BCA) Vintec® or sterile water, 5 days before fungal PP inoculation (T-5d), (C) fungal PP inoculation (T0), and (D) sampling site for metabolomic study after 0, 3, and 6 culture weeks (T0, T3, and T6) or fungal colonization rate (T6).</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>**FIGURE 2** | Gene copy number per mg of wood of *Phaeoacremonium minimum* (black) and *Phaeomoniella chlamydospora* (gray) 6 weeks post-infection (T6w). INi = control, injured samples; IV = samples inoculated with Vintec®; IPP = samples infected with *P. chlamydospora* and *P. minimum*; IVPP = samples pre-inoculated with Vintec® and infected with *P. chlamydospora* and *P. minimum*. N = 3. * Significance: 0.1234 (ns), 0.0332(*), 0.0021(**), < 0.0001(****).</figcaption></figure>
application of the MS-CleanR workflow. Among these features, 85% of them were annotated using different DBs: 3% (*Trichoderma*), 1% (Togniniaceae), 4% (Vitaceae), 42% (*Vitis*), and 35% (generic internal MS-FINDER DBs). Fungal metabolites were found in samples inoculated, and their intensity was between 5 × 10⁵ and 3 × 10⁶. As the background threshold was estimated to 1.5 × 10⁶, they appeared in the minority not only in terms of number but also in terms of intensity. Thus, separation of samples in principal component analysis (PCA) was mainly due to plant metabolites, and we only focused on them in the rest of the paper.

PCA score plots of respective kinetics T3 (Figure 3A) and T6 (Figure 3B), providing an unsupervised overview of UHPLC-HRMS fingerprints, marked sample-to-sample variability particularly for PP-infected wood at T3 and T6. It, respectively, displayed 42% and 32.7% of total explained variance using the first two principal components. Quality control samples

of the whole dataset (Supplementary Figure 1) were centered on the PCA plot, demonstrating the high reproducibility of the analysis.

More generally, a separation of samples can be seen according to treatment type. Non-treated wood composed of INi class and wood infected with fungi composed of IPP and IVPP classes are separated by PC1 in both kinetics. Interestingly, wood samples only treated with Vintec® are separated by PC1 for kinetics T3 and by PC2 for kinetics T6.

Using the MS-CleanR workflow (Fraisier-Vannier et al., 2020), a cleaner mass spectral similarity network was built to highlight common chemical classes related to plant response to pathogens and Vintec® (Figure 4). Potential biomarkers were revealed by ANOVA models to obtain a classification of the features (m/z-RT pairs) for each treatment in comparison with the reference.

All the two-to-two comparisons models were validated by random forest models based on the ROC curve coupled with

<figure><img src="image_5.png" /><figcaption>**FIGURE 3** | PCA score plot of ESI-PI data from wood extract of Vitis cuttings of kinetics T3 (A) and kinetics T6 (B) (in color). Colored circles arbitrarily enclose treatment types: Injured/Not inoculated (INi) samples; Injured/Vintec® (IV); Injured/*P. chlamydospora* + *P. minimum* (IPP) and Injured/*P. chlamydospora* + *P. minimum* + Vintec® (IVPP).</figcaption></figure>
<figure><img src="image_6.png" /></figure>

the area under the curve (AUC). We observed a high AUC score around 0.9 (Supplementary Table 2) indicating a good prediction ability.

Sixteen significant biomarkers correlated to the highest ANOVA coefficient scores were tentatively annotated by interrogating simultaneously *Vitis*, *Vitaceae*, *Trichoderma*, *Togniniaceae*, and local databases of natural products in MS-FINDER (Table 1). Three main biomarker clusters were observed, and annotation results gathered from the aligned peak list of LCMS features allowed the highlight of three main chemical classes: stilbenoid, flavonoid, and glycerophospholipids. The molecular formulae of these compounds matched structures already identified in the *Vitis* genus, except the three glycerophospholipids found in generic databases. No biomarkers were annotated with fungal databases.

The majority of the biomarkers belong to the stilbenoid class. They are resveratrol monomers (compounds 1, 2, and 3), dimers (compounds 4, 5, and 6), and trimers (compounds 12, 13, 14, and 15). Four of them displayed the identification level 1 according to the Metabolomics Standards Initiative (MSI) (Sumner et al., 2007) as their retention time, and MS/MS fragments perfectly match with the standards injected in the same conditions as samples (RT similarity score > 850; similarity

spectrum score > 800): oxyresveratrol, pterostilbene, α-viniferin, and viniferin. For all other compounds, the identification level 2 was proposed according to the MSI.

The three glycerophospholipids were annotated as lysophosphatidylcholines lysoPC (18:2) (9) and lysoPC (18:3) (8), and as lysophosphatidylethanolamine LysoPE (18:3) (7).

Two flavonoids were annotated as procyanidin derivatives (10 and 16), while the last one was annotated as an apigenin derivative (11).

After tentative annotation, we compared their distribution among treated conditions. On the heatmap representing the average peak intensity by class of each biomarker (Figure 5), it can be seen that these significant features are overproduced in the presence of pathogenic fungi alone (IPP) and in combination with BCA Vintec® (IVPP), validating our approach. In addition, the hierarchical clustering analysis, based on the sample class, separates well the non-infected wood samples, including those with Vintec® alone, from the pathogen-inoculated infected wood samples.

Among the different metabolites, flavonoids are mainly present at kinetics T3 and T6 with IPP and IVPP (metabolites 10, 11, and 16) and glycerophospholipids at kinetics T3 with IPP and IVPP (metabolites 7, 8, and 9). Suffruticosol A (13
<table><caption>**TABLE 1** | Summary of annotated compounds most differentially present in the tested conditions.</caption><thead><tr><th>N°</th><th>m/z</th><th>RT (min)</th><th>Molecular Formula</th><th>Error (ppm)</th><th>Chemical class<sup>a</sup></th><th>Putative annotation</th></tr></thead><tbody><tr><td>1</td><td>245.08081</td><td>10.608</td><td>C<sub>14</sub>H<sub>12</sub>O<sub>4</sub></td><td>4.5781</td><td>Stilbenoid</td><td>Oxyresveratrol<sup>b</sup></td></tr><tr><td>2</td><td>257.11743</td><td>16.784</td><td>C<sub>16</sub>H<sub>16</sub>O<sub>3</sub></td><td>3.4537</td><td>Stilbenoid</td><td>Pterostilbene<sup>b</sup></td></tr><tr><td>3</td><td>393.15448</td><td>10.065</td><td>C<sub>20</sub>H<sub>24</sub>O<sub>8</sub></td><td>2.5715</td><td>Stilbenoid</td><td>1-(3,5-Dihydroxyphenyl)-2-(4-hydroxyphenyl)ethane;<br/>4'-O-β-D-glucopyranoside<sup>c</sup></td></tr><tr><td>4</td><td>455.14893</td><td>10.863</td><td>C<sub>28</sub>H<sub>20</sub>O<sub>5</sub></td><td>2.3772</td><td>Stilbenoid</td><td>Leachianol F<sup>c</sup></td></tr><tr><td>5</td><td>455.14883</td><td>11.139</td><td>C<sub>28</sub>H<sub>22</sub>O<sub>6</sub></td><td>2.5969</td><td>Stilbenoid</td><td>Viniferin<sup>b</sup></td></tr><tr><td>6</td><td>455.14868</td><td>13.502</td><td>C<sub>28</sub>H<sub>22</sub>O<sub>6</sub></td><td>2.9265</td><td>Stilbenoid</td><td>Ampelopsin D<sup>c</sup></td></tr><tr><td>7</td><td>476.27747</td><td>16.766</td><td>C<sub>30</sub>H<sub>37</sub>NO<sub>4</sub></td><td>4.2820</td><td>Glycerophospholipids</td><td>LysoPE[18:3(9Z,12Z,15Z)/0:0]<sup>d</sup></td></tr><tr><td>8</td><td>518.32422</td><td>16.922</td><td>C<sub>26</sub>H<sub>48</sub>NO<sub>7</sub>P</td><td>1.9158</td><td>Glycerophospholipids</td><td>LysoPC [18:3(9Z,12Z,15Z)]<sup>d</sup></td></tr><tr><td>9</td><td>520.33972</td><td>17.781</td><td>C<sub>26</sub>H<sub>50</sub>NO<sub>7</sub>P</td><td>2.1966</td><td>Glycerophospholipids</td><td>LysoPC [18:2(9Z,12Z)]<sup>d</sup></td></tr><tr><td>10</td><td>593.12939</td><td>9.237</td><td>C<sub>30</sub>H<sub>24</sub>O<sub>13</sub></td><td>1.1363</td><td>Flavonoid</td><td>Procyanidin A1<sup>c</sup></td></tr><tr><td>11</td><td>609.12372</td><td>8.31</td><td>C<sub>30</sub>H<sub>24</sub>O<sub>14</sub></td><td>2.0669</td><td>Flavonoid</td><td>Apigenin 7-glycosides; 7-O-[3,4-Dihydroxy-E-cinnamoyl-(→4)-β-D-glucuronopyranoside]<sup>c</sup></td></tr><tr><td>12</td><td>679.19611</td><td>13.571</td><td>C<sub>42</sub>H<sub>30</sub>O<sub>9</sub></td><td>1.8345</td><td>Stilbenoid</td><td>α-viniferin<sup>b</sup></td></tr><tr><td>13</td><td>679.19611</td><td>14.085</td><td>C<sub>42</sub>H<sub>30</sub>O<sub>9</sub></td><td>1.8345</td><td>Stilbenoid</td><td>Suffruticosol A<sup>c</sup></td></tr><tr><td>14</td><td>681.21185</td><td>14.168</td><td>C<sub>42</sub>H<sub>32</sub>O<sub>9</sub></td><td>0.0866</td><td>Stilbenoid</td><td>Diptoindonesin B<sup>c</sup></td></tr><tr><td>15</td><td>681.21155</td><td>11.753</td><td>C<sub>42</sub>H<sub>32</sub>O<sub>9</sub></td><td>0.5270</td><td>Stilbenoid</td><td>cis-Miyabenol C<sup>c</sup></td></tr><tr><td>16</td><td>745.13983</td><td>10.039</td><td>C<sub>37</sub>H<sub>28</sub>O<sub>17</sub></td><td>1.6010</td><td>Flavonoid</td><td>3,3',4',5,5',7-Hexahydroxyflavan-(2→7,4→8)-<br/>3,3',4',5,7-pentahydroxyflavan;<br/>3'-O-(3,4,5-Trihydroxybenzoyl)<sup>c</sup></td></tr></tbody></table>

<figure><img src="image_7.png" /><figcaption>**FIGURE 5** | Heatmap of the annotated biomarkers at kinetics 3 weeks (T3) and 6 weeks (T6). The color key is based on the average peak intensity of each feature by class: red color for higher peak intensity and blue color for lower peak intensity. Injured/Not inoculated (INi) samples; Injured/Vintec® (IV); Injured/P. chlamydospora + P. minimum (IPP) and Injured/P. chlamydospora + P. minimum + Vintec® (IVPP). Report to **Table 1** for metabolite annotation.</figcaption></figure>

and α-viniferin (12), two resveratrol trimers, are overall equally expressed with IPP and IVPP at both kinetics. Cis-miyabenol C (15) and ampelopsin D (6) are mainly overproduced with IVPP at kinetics T3, while pterostilbene (2) is mainly overproduced with

IVPP at kinetics T6. Oxyresveratrol (1), the bibenzyl derivative (3), and diptoindonesin B (14) are mainly overproduced with IPP at kinetics T6. Finally, all other annotated resveratrol dimers (metabolites 4 and 5) are overproduced with IPP and

<footnote><sup>a</sup>Determined with ClassyFire (Djoumbou Feunang et al., 2016).</footnote>

<footnote><sup>b</sup>Compounds confirmed by injection of standard in the same analytical conditions.</footnote>

<footnote><sup>c</sup>Putative annotation based on experimental HRMS, MS/MS, and in silico fragmentation matches restricting interrogation to Vitis genus and Vitaceae family databases.</footnote>

<footnote><sup>d</sup>Putative annotation based on experimental HRMS, MS/MS, and in silico fragmentation matches by interrogation of natural product databases within the MS-FINDER software.</footnote>

<footnote>* Significance in comparison to INi for each kinetics independently [0.0332(*), 0.0021(**), 0.0002(***), < 0.0001(****)].</footnote>
<figure><img src="image_8.png" /><figcaption>**FIGURE 6** | Two-way ANOVA and Dunnett's *post hoc* test (p ≤ 0.05) of metabolites overproduced with the combination Vintec®-pathogens. 0.0021(**), 0.0002(***), < 0.0001(****).</figcaption></figure>

IVPP at kinetics T6. ANOVA comparisons can be found in **Supplementary Figure 2**.

Interestingly, all biomarkers overproduced with IVPP were not significantly overproduced in the presence of BCA Vintec® alone, as well as for kinetics T3 and for kinetics T6. Thus, it resulted that Vintec® treatment alone induces a weak grapevine defense response, which seemed rather correlated with plant reaction to pathogenic fungal infection. Nevertheless, some metabolites, in particular oxyresveratrol (1), the bibenzyl derivative (3), and diptoindonesin B (14), were less produced with pathogens in combination with Vintec® than with pathogens alone.

Thus, Vintec® seems to reduce some metabolic response to pathogen attack.

From another side, Vintec® increases plant response to pathogen treatment (Figure 6), in particular regarding the expression of stilbenoids pterostilbene (2) at T6 and cis-miyabenol C (15) at T3, as well as flavonoids 10, 11, and 16 at T3. To a lesser extent, it also stimulates the production of lysoPCs 8 and 9 at T6 and ampelopsin D at T3.

The reproducibility of the experiment and the raw materials were confirmed for kinetics T3 with canes of *Vitis vinifera* L. cv Cabernet-Sauvignon clone 15 during the 2021 season (data not shown). UHPLC-HRMS analysis revealed a metabolic response
quite constant as 11 of the 16 highlighted biomarkers were also overproduced in IPP and IVPP samples comparatively to control or BCA Vintec® treatment of the 2021 season (Supplementary Figure 3).

## Validation of Stilbenoids Antifungal Activity

Antifungal activity of four stilbenoids standards was investigated on *Botrytis cinerea* and *P. chlamydospora* fungal growth. IC₅₀ was determined (Table 2), and the results pointed out that all stilbenes showed an antifungal activity on both fungi, *P. chlamydospora* being the more susceptible. Pterostilbene and viniferin had the highest antifungal activities for both fungi, whereas oxyresveratrol presented the lowest activity.

## DISCUSSION

The objective of this study was to characterize the grapevine wood metabolome changes, during a 6-week kinetics, after the inoculation of the Esca pathogens *Phaeomoniella chlamydospora* and *Phaeoacremonium minimum*, and of the commercialized BCA Vintec® formulated with *Trichoderma atroviridae* SC1, within the framework of Esca disease. This work was expected to identify some biomarkers of Esca attack that could be used to further evaluate BCA efficiency in vineyard conditions. Even if BCA can exhibit many different modes of action on targeted disease, the objective was to demonstrate the prime interest of metabolomics on the semi-ligneous and quite challenging Esca pathosystem.

From an analytical point of view, the PI mode was preferred as MS/MS spectra were of better quality than in the NI mode, thus improving feature annotation using MS-FINDER and network feature clusterization. In order to highlight significant features based on pathogens and BCA presence, a classification model between treatment and UHPLC-HRMS fingerprints was realized by random forest. In addition, a mass spectral similarity network was built and allowed the distinction of three main significant clusters from three chemical classes: flavonoids, glycerophospholipids, and mainly stilbenoids as resveratrol monomers, dimers, and trimers. The majority of highlighted biomarkers were annotated from a *Vitis* database. Interestingly, none of them were fungal metabolites even if nearly 5% of all features were annotated with fungal databases. Indeed, the peak intensity of these metabolites was lower than the one of plant metabolites, and they were often close to the background noise demonstrating the prevalence of wood metabolites using this untargeted metabolomic approach. However, their presence is a

very good indicator of the presence of pathogens as well as BCA. These results can stimulate further investigation to follow the presence/absence of these different partners in the time course of Esca disease, which show a slow decaying evolution. This could be particularly helpful to understand the first steps of this disease especially when there are no symptoms.

Stilbenoids belong to the phenylpropanoid pathway and include stilbenes, bibenzyls, and phenanthrenes derivatives. They are found as constitutive compounds in several plant families especially Vitaceae (Rivière et al., 2012). As mentioned above, they are synthesized through the phenylpropanoid pathway as well as the flavonoids. However, there is an early divergence between these classes at the p-coumaroylCoA level that can enter either the stilbenoid way through the key stilbene synthase or the flavonoids way thanks to chalcone synthase isomerase (Ali et al., 2011). They are widely considered as phytoalexins as they play a protective role in plant defense to fungal pathogen infection, such as the detected pterostilbene in the presence of *Plasmopara viticola* (Langcake et al., 1979; Pezet et al., 2004). These compounds are also of great interest since they exhibit many biological properties, in particular cardioprotection, neuroprotection, antidiabetic properties, depigmentation, anti-inflammation, cancer prevention, and treatment (Akinwumi et al., 2018). Several other particular stilbenoids were previously reported in literature and isolated from diseased grapevine wood: the resveratrol dimers leachianol F as well as the resveratrol trimers α-viniferin and cis-miyabenol C (Amalfitano et al., 2011). It was also previously shown that infection of *Vitis vinifera* with *Phaeomoniella chlamydospora* induced changes in phenolic compounds with a significant increase in their production after this fungal inoculation. The production of viniferin-type stilbenes was, in particular, observed, absent from the control condition (Lima et al., 2012).

The role of stilbenes as plant defense compounds was strengthened by antifungal bioassays on four available standards. Our results on positive control *Botrytis cinerea* support the data from literature. It was previously shown that ED₅₀ (effective dose to obtain 50% mortality) of pterostilbene on *B. cinerea* was around 18 to 20 µg/ml (Jeandet et al., 2002). Regarding α-viniferin, an IC₅₀ around 97 µg/ml was observed on *B. cinerea* (Langcake and Pryce, 1977). *P. chlamydospora* appeared to be more sensitive than *B. cinerea*. Pterostilbene and viniferin presented the highest antifungal activities for both fungi, whereas oxyresveratrol showed the lowest activity. A detailed analysis of stilbenoids antifungal activity on various fungi involved in wood disease showed a weak sensibility of *P. chlamydospora* and *P. aleophilum*. *Phaeomoniella chlamydospora* growth was inhibited by pterostilbene and two resveratrol tetramers, vitisins A and B, whereas no stilbenoids had an impact on *P. aleophilum* growth (Lambert et al., 2012). Nevertheless, they found that pterostilbene was the most active stilbenoids agreeing with our results. Still, the apparent divergence regarding sensitivity could be explained by the use of different fungal isolates and by the experimental conditions as our analysis was realized on fungal spores and not mycelium. These results support the key role of stilbenoids in the framework of wood disease.

<table><caption>**TABLE 2** | IC₅₀ of the standard solutions on *Botrytis cinerea* and *Phaeomoniella chlamydospora* spores 3 days after onset of treatment.</caption><thead><tr><th></th><th>IC<sub>50</sub> Botrytis (µg/ml)</th><th>IC<sub>50</sub> Pch (µg/ml)</th></tr></thead><tbody><tr><td>Oxyresveratrol (1)</td><td>162</td><td>40.0</td></tr><tr><td>Pterostilbene (2)</td><td>20.3</td><td>7.3</td></tr><tr><td>Viniferin (5)</td><td>37.5</td><td>8.6</td></tr><tr><td>α-viniferin (12)</td><td>92.8</td><td>13.5</td></tr></tbody></table>
Flavonoids are also a class of compounds known to play a role in plant defense against pathogens and other environmental stresses (Treutter, 2005). They are often constitutive compounds of the plant with a production enhanced by the stress. In *Vitis vinifera*, they are mainly distributed in leaves, stems, and canes (Goufo et al., 2020). In addition, accumulation of flavonoids, namely procyanidin derivatives, was also highlighted in wood parts of Esca-infected grapevines (Rusjan et al., 2017).

The third class of compounds identified was phospholipids. It was also previously reported that phospholipids lysophosphatidylcholines (lysoPCs), in particular lysoPCs C18:2 and C18:3 identified in this study, increase after pathogen inoculation (Cho et al., 2012). Regarding the leaves of Esca-affected grapevine, a lipidomic study also reported a progressive increase in lipids, including lysophospholipids as lysoPCs C18:2 and C18:3, with symptom progression (Goufo and Cortez, 2020). These compounds are enzymatically produced through phospholipase A activity linked to the signaling pathways of jasmonic acid leading to defense induction (Chandra et al., 1996). However, the precise role of lysoPCs in the interaction between wood tissues and wood pathogens needs further investigation. Nevertheless, lipids appeared to be good biomarkers of grapevine disease as a higher accumulation of them was also observed in asymptomatic diseased grapevine (Lemaitre-Guillier et al., 2020).

*Trichoderma* spp. strains are known to reduce infections caused by trunk disease pathogens, such as *Phaeomoniella chlamydospora* (Pertot et al., 2016) and *Phaeoacremonium minimum* (Carro-Huerga et al., 2020) in a preventive way, as we confirmed here at 6 weeks post-inoculation. Different modes of action have been described for *Trichoderma* spp. to explain their capacity to protect plants against different biotic and abiotic stresses. These explain (i) the competition for resources and space, (ii) the production of toxic compounds, and (iii) the stimulation of plant natural defenses. The ability of *Trichoderma* spp. to enhance plant defenses against different organisms has been already described several times in different pathosystems (Perazzolli et al., 2011; Mastouri et al., 2012; Mathys et al., 2012; Hermosa and Rubio, 2013; Banani et al., 2014; Singh et al., 2019). For instance, a recent study of the application on *T. harzianum* on tomato reported a protective effect against root-knot nematode by the overproduction of reactive oxygen species (ROS), secondary metabolites, and defense-related hormones (Yan et al., 2021). It has also been observed that the application of different *Trichoderma* spp. in eggplant resulted in a higher total phenolic content during the infection with *S. sclerotirum* (Pratap Singh et al., 2021).

Metabolomic results revealed that Vintec® composed of a *Trichoderma atroviride* (TASC1) strain induces a weak metabolomic response alone, sufficient enough to differentiate the treated and not treated plants by metabolomics comparison, but therefore does not have plant defense stimulating effect.

Nevertheless, BCA Vintec® seems to have not only an effect on pathogens but also a priming effect after infection with fungi. Indeed, some metabolites were more produced in IPP samples compared to IVPP samples, especially at T6. This observation could indicate a weaker virulence of pathogens in the presence of Vintec® corroborating BCA capacity to compete and

interact with these pathogens. This effect appears within a certain period probably due to the time necessary for the sufficient strain development. On the other hand, Vintec® increases plant response to challenge inoculation with a clear stimulation of the phenylpropanoid pathway with increasing amounts of stilbenoid pterostilbene as well as an increase in flavonoids. Thus, these compounds appear to be good BCA indicators as they are still observable after 3 weeks and even 6 weeks after infection.

## CONCLUSION

A UHPLC-HRMS based metabolomic study was performed and allowed the dereplication of plant defense metabolites in the framework of Esca disease. Comparison with injured reference (INi) was conducted and significant metabolic changes were observed. Sixteen biomarkers were highlighted and annotated based on MS/MS fragmentation patterns. Four of them were clearly identified thanks to the injection of corresponding standards. These compounds belong mainly to the stilbenoid chemical class, as well as flavonoids and glycerophospholipids known to be involved in plant defense. They were mainly accumulated after two wood treatments: inoculation of the pathogenic fungi, *Phaeoacremonium minimum* and *Phaeomoniella chlamydospora*, and in combination with BCA Vintec®. One stilbenoid and three flavonoids were highly associated with the synergy between the BCA and the pathogens. Indeed, they were mainly overproduced when Vintec® was in the presence of fungi. Thus, these four compounds appear to be good indicators of BCA effects in field condition.

## DATA AVAILABILITY STATEMENT

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: 10.5281/zenodo.579519.

## AUTHOR CONTRIBUTIONS

GM, VP-P, SF, BD, and AJ supervised the project and revised the article. JC was in charge of the plant extraction, the analytical workflow, and the data processing and interpretation. AR-O was in charge of the plant culture and treatment and realized the bioassays. JC and AR-O wrote the article. All authors have read and agreed the published version of the manuscript.

## FUNDING

Financial support was received from the French National Infrastructure for Metabolomics and Fluxomics, Grant MetaboHUB-ANR-11-INBS-0010, and the PSPC SOLSTICE Project (SOLutionS pour des Traitements Intégrés dans une Conduite Environnementale) managed by Belchim Crop
Protection and partly funded by the French state within the framework of the Programme d'Investissements d'Avenir.

## ACKNOWLEDGMENTS

We thank Wilfried Remus-Borel and Gabriel D' Enjoy from Belchim Crop Protection for fruitful discussions and for supplying the Vintec® product. MetaToul (Toulouse metabolomics & fluxomics facilities, www.metatoul.fr) is part

of the French National Infrastructure for Metabolomics and Fluxomics MetaboHUB-ANR-11-INBS-0010.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.835463/full#supplementary-material

## REFERENCES

Akinwumi, B., Bordun, K.-A., and Anderson, H. (2018). Biological activities of stilbenoids. *Int. J. Mol. Sci.* 19:792. doi: 10.3390/ijms19030792

Ali, M. B., Howard, S., Chen, S., Wang, Y., Yu, O., Kovacs, L. G., et al. (2011). Berry skin development in norton grape: distinct patterns of transcriptional regulation and flavonoid biosynthesis. *BMC Plant Biol.* 11:7. doi: 10.1186/1471-2229-11-7

Amalfitano, C., Agrelli, D., Arrigo, A., Mugnai, L., Surico, G., and Evidente, A. (2011). Stilbene Polyphenols in the brown red wood of *Vitis vinifera* Cv. sangiovese affected by "Esca Proper". *Phytopathol. Mediter.* 50, S224-S235.

Amalfitano, C., Evidente, A., Surico, G., Tegli, S., Bertelli, E., and Mugnai, L. (2000). Phenols and stilbene polyphenols in the wood of esca-diseased grapevines. *Phytopathol. Mediter.* 39:6.

Banani, H., Roatti, B., Ezzahi, B., Giovannini, O., Gessler, G., Pertot, I., et al. (2014). Characterization of resistance mechanisms activated by *Trichoderma harzianum* T39 and benzothiadiazole to downy mildew in different grapevine cultivars. *Plant Pathol.* 63, 334-343. doi: 10.1111/ppa.12089

Berbegal, M., Ramón-Albalat, A., León, M., and Armengol, J. (2020). Evaluation of long-term protection from nursery to vineyard provided by trichoderma atroviride SC1 against fungal grapevine trunk pathogens. *Pest Manag. Sci.* 76, 967-977. doi: 10.1002/ps.5605

Bertsch, C., Ramírez-Suero, M., Magnin-Robert, M., Larignon, P., Chong, J., Abou-Mansour, E., et al. (2013). Grapevine trunk diseases: complex and still poorly understood: grapevine trunk diseases. *Plant Pathol.* 62, 243-265. doi: 10.1111/j.1365-3059.2012.02674.x

Braun, H., Woitsch, L., Hetzer, B., Geisen, R., Zange, B., and Schmidt-Heydt, M. (2018). *Trichoderma harzianum*: inhibition of mycotoxin producing fungi and toxin biosynthesis. *Int. J. Food Microbiol.* 280, 10-16. doi: 10.1016/j.ijfoodmicro. 2018.04.021

Bruez, E., Baumgartner, K., Bastien, S., Travadon, R., Guérin-Dubrana, L., and Rey, P. (2016). Various fungal communities colonise the functional wood tissues of old grapevines externally free from grapevine trunk disease symptoms: fungal microflora of gtd-free old vines. *Austral. J. Grape Wine Res.* 22, 288-295. doi: 10.1111/ajgw.12209

Bruez, E., Vallance, J., Gerbore, J., Lecomte, P., Da Costa, J. P., Guerin-Dubrana, L., et al. (2014). Analyses of the temporal dynamics of fungal communities colonizing the healthy wood tissues of esca leaf-symptomatic and asymptomatic vines. *PLoS One* 9:e95928. doi: 10.1371/journal.pone.0095928

Bruno, G., and Sparapano, L. (2007). Effects of three esca-associated fungi on *Vitis vinifera* L.: V. changes in the chemical and biological profile of xylem sap from diseased Cv. sangiovese vines. *Physiol. Mol. Plant Pathol.* 71, 210-229. doi: 10.1016/j.pmpp.2008.02.005

Carro-Huerga, G., Compant, S., Gorfer, M., Cardoza, R. E., Schmoll, M., Gutierrez, S., et al. (2020). Colonization of *Vitis vinifera* L. by the *Endophyte trichoderma* Sp. Strain T154: biocontrol activity against phaeoacremonium minimum. *Front. Plant Sci.* 11:15. doi: 10.3389/fpls.2020.01170

Chandra, S., Heinstein, P. F., and Low, P. S. (1996). Activation of phospholipase a by plant defense elicitors. *Plant Physiol.* 110, 979-986. doi: 10.1104/pp.11 0.3.979

Cho, K., Kim, Y., Wi, S. J., Seo, J. B., Kwon, J., Chung, J. H., et al. (2012). Nontargeted metabolite profiling in compatible pathogen-inoculated tobacco

(Nicotiana tabacum L. Cv. Wisconsin 38) Using UPLC-Q-TOF/MS. *J. Agric. Food Chem.* 60, 11015-11028. doi: 10.1021/jf303702j

Chong, J., Wishart, D. S., and Xia, J. (2019). Using metaboanalyst 4.0 for comprehensive and integrative metabolomics data analysis. *Curr. Protoc. Bioinformatics* 68:e86. doi: 10.1002/cpbi.86

Del Rio, J. A., Gonzalez, A., Fuster, M. D., Botia, J. M., Gomez, P., Frias, V., et al. (2001). Tylose formation and changes in phenolic compounds of grape roots infected with *Phaeomoniella chlamydospora* and Phaeoacremonium Species. *Phytopathol. Mediter.* 40, S394-S399.

Demsar, J., Curk, T., Erjavec, A., Demsar, J., Curk, T., Erjave, A., et al. (2013). Orange: data mining toolbox in python. *J. Mach. Learn. Res.* 14, 2349-2353.

Djoumbou Feunang, Y., Eisner, R., Knox, C., Chepelev, L., Hastings, J., Owen, G., et al. (2016). ClassyFire: automated chemical classification with a comprehensive, computable taxonomy. *J. Cheminform.* 8:61. doi: 10.1186/ s13321-016-0174-y

Fontaine, F., Pinto, C., Vallet, J., Clément, C., Gomes, A. C., and Spagnolo, A. (2016). The effects of grapevine trunk diseases (GTDs) on vine physiology. *Eur. J. Plant Pathol.* 144, 707-721. doi: 10.1007/s10658-015-0770-0

Fraisier-Vannier, O., Chervin, J., Cabanac, G., Puech, V., Fournier, S., Durand, V., et al. (2020). MS-CleanR: a feature-filtering workflow for untargeted LC-MS based metabolomics. *Anal. Chem.* 92, 9971-9981. doi: 10.1021/acs.analchem. 0c01594

Gerbore, J., Benhamou, N., Vallance, J., Le Floch, G., Grizard, D., Regnault-Roger, C., et al. (2014). Biological control of plant pathogens: advantages and limitations seen through the case study of *Pythium oligandrum*. *Environ. Sci. Pollut. Res.* 21, 4847-4860. doi: 10.1007/s11356-013-1807-6

Goufo, P., and Cortez, I. (2020). A lipidomic analysis of leaves of esca-affected grapevine suggests a role for galactolipids in the defense response and appearance of foliar symptoms. *Biology* 9:268. doi: 10.3390/biology90 90268

Goufo, P., Singh, R. K., and Cortez, I. (2020). A reference list of phenolic compounds (including stilbenes) in grapevine (*Vitis vinifera* L.) roots, woods, canes, stems, and leaves. *Antioxidants* 9:398. doi: 10.3390/antiox9050398

Graniti, A., Surico, G., and Mugnai, L. (2000). Esca of grapevine: a disease complex or a complex of diseases?*. *Phytopathol. Mediter.* 39:5.

Haidar, R., Amira, Y., Roudet, J., Marc, F., and Patrice, R. (2021). Application Methods and Modes of Action of *Pantoea agglomerans* and *Paenibacillus* Sp. to control the grapevine trunk disease-pathogen, *Neofusicoccum parvum*. *OENO One* 55, 1-16. doi: 10.20870/oeno-one.2021.55.3.4530

Haidar, R., Yacoub, A., Pinard, A., Roudet, J., Fermaud, M., and Rey, P. (2020). Synergistic effects of water deficit and wood-inhabiting bacteria on pathogenicity of the grapevine trunk pathogen *Neofusicoccum parvum*. *Phytopathol. Mediter.* 59, 473-484. doi: 10.14601/Phyto-11990

Hermosa, R., and Rubio, M. B. (2013). The contribution of trichoderma to balancing the costs of plant growth and defense. *Int. Microbiol.* 16, 69-80. doi: 10.2436/20.1501.01.181

Hofstetter, V., Buyck, B., Croll, D., Viret, O., Couloux, A., and Gindro, K. (2012). What if esca disease of grapevine were not a fungal disease?* *Fungal Divers.* 54, 51-67. doi: 10.1007/s13225-012-0171-z

Jeandet, P., Douillet-Breuil, A.-C., Bessis, R., Debord, S., Sbaghi, M., and Adrian, M. (2002). Phytoalexins from the Vitaceae: biosynthesis, phytoalexin gene expression in transgenic plants, antifungal activity, and metabolism. *J. Agric. Food Chem.* 50, 2731-2741. doi: 10.1021/jf011429s
Labois, C., Wilhelm, K., Laloue, H., Tarnus, C., Bertsch, C., Goddard, M.-L., et al. (2020). Wood metabolomic responses of wild and cultivated grapevine to infection with *Neofusicoccum parvum*, a trunk disease pathogen. *Metabolites* 10:232. doi: 10.3390/metabo10060232

Lambert, C., Bisson, J., Waffo-Téguo, P., Papastamoulis, Y., Richard, T., Corio-Costet, M.-F., et al. (2012). Phenolics and their antifungal role in grapevine wood decay: focus on the Botryosphaeriaceae family. *J. Agric. Food Chem.* 60, 11859–11868. doi: 10.1021/jf303290g

Langcake, P. (1981). Disease resistance of *Vitis* Spp. and the production of the stress metabolites resveratrol, ε-Viniferin, α-Viniferin and pterostilbene. *Physiol. Plant Pathol.* 18, 213–226. doi: 10.1016/S0048-4059(81)8 0043-4

Langcake, P., Cornford, C. A., and Pryce, R. J. (1979). Identification of pterostilbene as a Phytoalexin from *Vitis vinifera* leaves. *Phytochemistry* 18, 1025–1027. doi: 10.1016/S0031-9422(00)91470-5

Langcake, P., and Pryce, R. J. (1977). A new class of phytoalexins from grapevines. *Experientia* 33, 151–152. doi: 10.1007/BF02124034

Lemaitre-Guillier, C., Fontaine, F., Roullier-Gall, C., Harir, M., Magnin-Robert, M., Clément, C., et al. (2020). Cultivar- and Wood area-dependent metabolic fingerprints of grapevine infected by Botryosphaeria Dieback. *Phytopathology* 110, 1821–1837. doi: 10.1094/PHYTO-02-20-0055-R

Lima, M. R. M., Felgueiras, M. L., Cunha, A., Chicau, G., Ferreres, F., and Dias, A. C. P. (2017a). Differential phenolic production in leaves of *Vitis vinifera* Cv. alvarinho affected with esca disease. *Plant Physiol. Biochem.* 112, 45–52. doi: 10.1016/j.plaphy.2016.12.020

Lima, M. R. M., Machado, A. F., and Gubler, W. D. (2017b). Metabolomic study of chardonnay grapevines double stressed with esca-associated fungi and drought. *Phytopathology* 107, 669–680. doi: 10.1094/PHYTO-11-16- 0410-R

Lima, M. R. M., Ferreres, F., and Dias, A. C. P. (2012). Response of *Vitis vinifera* Cell Cultures to *Phaeomoniella chlamydospora*: changes in phenolic production, oxidative state and expression of defence-related genes. *Eur. J. Plant Pathol.* 132, 133–146. doi: 10.1007/s10658-011-9857-4

Luque, J., Martos, S., Aroca, A., Raposo, R., and Garcia-Figueres, F. (2009). Symptoms and fungi associated with declining mature grapevine plants in Northeast Spain. *J. Plant Pathol.* 91, 381–390.

Magnin-Robert, M., Adrian, M., Trouvelot, S., Spagnolo, A., Jacquens, L., Letousey, P., et al. (2017). Alterations in grapevine leaf metabolism occur prior to esca apoplexy appearance. *Mol. Plant Microbe Interact.* 30, 946–959. doi: 10.1094/ MPMI-02-17-0036-R

Mastouri, F., Björkman, T., and Harman, G. E. (2012). *Trichoderma harzianum* enhances antioxidant defense of tomato seedlings and resistance to water deficit. *Mol. Plant Microbe Interact.* 25, 1264–1271. doi: 10.1094/MPMI-09-11- 0240

Mathys, J., De Cremer, K., Timmermans, P., Van Kerckhove, S., Lievens, B., Vanhaecke, M., et al. (2012). Genome-wide characterization of ISR induced in *Arabidopsis thaliana* by *Trichoderma hamatum* T382 against botrytis cinerea infection. *Front. Plant Sci.* 3:108. doi: 10.3389/fpls.2012.00108

Mondello, V., Larignon, P., Armengol, J., Kortekamp, A., Vaczy, K., Prezman, F., et al. (2018a). Management of grapevine trunk diseases: knowledge transfer, current strategies and innovative strategies adopted in Europe. *Phytopathol. Mediter.* 57, 369–383. doi: 10.14601/Phytopathol_Mediterr-23942

Mondello, V., Songy, A., Battiston, E., Pinto, C., Coppin, G., Trotel-Aziz, P., et al. (2018b). Grapevine trunk diseases: a review of fifteen years of trials control with chemicals and biocontrol agents. *Plant Dis.* 102, 1189–1217. doi: 10.1094/PDIS- 08-17-1181-FE

Moya, P., Barrera, V., Cipollone, J., Bedoya, C., Kohan, L., Toledo, A., et al. (2020). New isolates of *Trichoderma* Spp. as biocontrol and plant growth-promoting agents in the pathosystem pyrenophora teres-barley in Argentina. *Biol. Control* 141:104152. doi: 10.1016/j.biocontrol.2019.10 4152

Mugnai, L., Graniti, A., and Surico, G. (1999). Esca (Black Measles) and brown wood-streaking: two old and elusive diseases of grapevines. *Plant Dis.* 83, 404–418. doi: 10.1094/PDIS.1999.83.5.404

Perazzolli, M., Roatti, B., Bozza, E., and Pertot, I. (2011). *Trichoderma harzianum* T39 induces resistance against downy mildew by priming for defense without costs for grapevine. *Biol. Control* 58, 74–82. doi: 10.1016/j.biocontrol.2011.04. 006

Pertot, I., Prodrutti, D., Colombini, A., and Pasini, L. (2016). *Trichoderma atroviride* SC1 prevents *Phaeomoniella chlamydospora* and *Phaeoacremonium aleophilum* infection of grapevine plants during the grafting process in nurseries. *BioControl* 61, 257–267.

Pezet, R., Gindro, K., Viret, O., and Richter, H. (2004). Effects of resveratrol, viniferins and pterostilbene on *Plasmopara viticola* zoospore mobility and disease development. *Vitis. J. Grapevine Res.* 43, 145–148.

Pierron, R., Gorfer, M., Berger, H., Jacques, A., Sessitsch, A., Strauss, J., et al. (2015). Deciphering the niches of colonisation of *Vitis vinifera* L. by the esca-associated fungus *Phaeoacremonium aleophilum* using a Gfp marked strain and cutting systems. *PLoS One* 10:e0126851. doi: 10.1371/journal.pone.0126851

Pouzoulet, J., Mailhac, N., Couderc, C., Besson, X., Daydé, J., Lummerzheim, M., et al. (2013). A method to detect and quantify phaeomoniella chlamydospora and *Phaeoacremonium aleophilum* DNA in grapevine-wood samples. *Appl. Microbiol. Biotechnol.* 97, 10163–10175. doi: 10.1007/s00253-013-5299-6

Poveda, J. (2021). Trichoderma as biocontrol agent against pests: new uses for a mycoparasite. *Biol. Control* 159:104634. doi: 10.1016/j.biocontrol.2021.1 04634

Pratap Singh, S., Keswani, C., Pratap Singh, S., Sansinenea, E., and Xuan Hoat, T. (2021). Trichoderma Spp. Mediated induction of systemic defense response in brinjal against *Sclerotinia sclerotiorum*. *Curr. Res. Microb. Sci* 2:100051. doi: 10.1016/j.crmicr.2021.100051

Ram, R. M., Keswani, C., Bisen, K., Tripathi, R., Singh, S. P., and Singh, H. B. (2018). "Biocontrol technology: eco-friendly approaches for sustainable agriculture," in *Omics Technologies and Bio-Engineering*, Vol. 2, eds D. Barh and V. Azevedo (San Diego, CA: Academic Press), 177–190. doi: 10.1016/B978-0- 12-815870-8.00010-3

Reis, P., Pierron, R., Larignon, P., Lecomte, P., Abou-Mansour, E., Farine, S., et al. (2019). Vitis methods to understand and develop strategies for diagnosis and sustainable control of grapevine trunk diseases. *Phytopathology* 109, 916–931. doi: 10.1094/PHYTO-09-18-0349-RVW

Rivière, C., Pawlus, A. D., and Mérillon, J.-M. (2012). Natural stilbenoids: distribution in the plant kingdom and chemotaxonomic interest in vitaceae. *Nat. Prod. Rep.* 29:1317. doi: 10.1039/c2np20049J

Romeo-Oliván, A., Pagés, M., Breton, C., Lagarde, F., Cros, H., Yobrégat, O., et al. (2021). Ozone dissolved in water: an innovative tool for the production of young plants in grapevine nurseries? *Ozone Sci. Eng.* 1–15. doi: 10.1080/ 01919512.2021.1984203

Rusjan, D., Persic, M., Likar, M., Biniari, K., and Mikulic-Petkovsek, M. (2017). Phenolic responses to esca-associated fungi in differently decayed grapevine woods from different trunk parts of "Cabernet Sauvignon". *J. Agric. Food Chem.* 65, 6615–6624. doi: 10.1021/acs.jafc.7b02188

Savazzini, F., Longa, C. M. O., Pertot, I., and Gessler, C. (2008). Real-time PCR for detection and quantification of the biocontrol agent *Trichoderma atroviride* strain SC1 in soil. *J. Microbiol. Methods* 73, 185–194. doi: 10.1016/j.mimet.2008. 02.004

Singh, B. N., Padmanabh, D., Sarma, B. K., and Singh, H. B. (2019). *Trichoderma Asperellum* T42 induces local defense against *Xanthomonas oryzae* Pv. Oryzae under nitrate and ammonium nutrients in tobacco. *RSC Adv.* 9, 39793–39810. doi: 10.1039/C9RA06802C

Sumner, L. W., Amberg, A., Barrett, D., Beale, M. H., Beger, R., Daykin, C. A., et al. (2007). Proposed minimum reporting standards for chemical analysis: chemical analysis working group (CAWG) metabolomics standards initiative (MSI). *Metabolomics* 3, 211–221. doi: 10.1007/s11306-007-0082-2

Surico, G. (2009). Towards a redefinition of the diseases within the esca complex of grapevine. *Phytopathol. Mediter.* 48, 5–10.

Treutter, D. (2005). Significance of flavonoids in plant resistance and enhancement of their biosynthesis. *Plant Biol.* 7, 581–591. doi: 10.1055/s-2005-8 73009

Trotel-Aziz, P., Abou-Mansour, E., Courteaux, B., Rabenoelina, F., Clément, C., Fontaine, F., et al. (2019). Bacillus subtilis PTA-271 counteracts botryosphaeria dieback in grapevine, triggering immune responses and detoxification of fungal phytoxins. *Front. Plant Sci.* 10:25. doi: 10.3389/fpls.2019.00025

Tsugawa, H., Cajka, T., Kind, T., Ma, Y., Higgins, B., Ikeda, K., et al. (2015). MS-DIAL: data-independent MS/MS deconvolution for comprehensive metabolome analysis. *Nat. Methods* 12, 523–526. doi: 10.1038/nmeth.3393

Tsugawa, H., Kind, T., Nakabayashi, R., Yukihira, D., Tanaka, W., Cajka, T., et al. (2016). Hydrogen rearrangement rules: computational MS/MS fragmentation
and structure elucidation using MS-FINDER software. *Anal. Chem.* 88, 7946-7958. doi: 10.1021/acs.analchem.6b00770

Úrbez-Torres, J. R., Haag, P., Bowen, P., and O’Gorman, D. T. (2014). Grapevine trunk diseases in British Columbia: incidence and characterization of the fungal pathogens associated with Esca and Petri diseases of grapevine. *Plant Dis.* 98, 469–482. doi: 10.1094/PDIS-05-13-0523-RE

Vinale, F., Sivasithamparam, K., Ghisalberti, E. L., Marra, R., Woo, S. L., and Lorito, M. (2008). Trichoderma–plant–pathogen interactions. *Soil Biol. Biochem.* 40, 1–10. doi: 10.1016/j.soilbio.2007.07.002

Yacoub, A., Gerbore, J., Magnin, N., Chambon, P., Dufour, M.-C., Corio-Costet, M.-F., et al. (2016). Ability of *Pythium oligandrum* strains to protect *Vitis vinifera* L., by inducing plant resistance against *Phaeomoniella chlamydospora*, a pathogen involved in Esca, a grapevine trunk disease. *Biol. Control* 92, 7–16. doi: 10.1016/j.biocontrol.2015.08.005

Yacoub, A., Haidar, R., Gerbore, J., Masson, C., Dufour, M.-C., Guyoneaud, R., et al. (2020). *Pythium Oligandrum* induces grapevine defence mechanisms against the trunk pathogen *Neofusicoccum parvum*. *Phytopathol. Mediter.* 59, 565–580. doi: 10.14601/Phyto-11270

Yan, Y., Mao, Q., Wang, Y., Zhao, J., Fu, Y., Yang, Z., et al. (2021). *Trichoderma harzianum* induces resistance to root-knot nematodes by increasing secondary metabolite synthesis and defense-related enzyme activity

in *Solanum lycopersicum* L. *Biol. Control* 158:104609. doi: 10.1016/j.biocontrol.2021.104609

**Conflict of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

**Publisher’s Note:** All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Copyright © 2022 Chervin, Romeo-Oliván, Fournier, Puech-Pages, Dumas, Jacques and Marti. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.