# Unravelling the interplay of nitrogen nutrition and the Botrytis cinerea pectin lyase BcPNL1 in modulating Arabidopsis thaliana susceptibility

Antoine Davière, Aline Voxeur, Sylvie Jolivet, Luka Lelas, Samantha

Vernhettes, Marie-Christine Soulié, Mathilde Fagard

► To cite this version:

Antoine Davière, Aline Voxeur, Sylvie Jolivet, Luka Lelas, Samantha Vernhettes, et al.. Unravelling the interplay of nitrogen nutrition and the Botrytis cinerea pectin lyase BcPNL1 in modulating Arabidopsis thaliana susceptibility. Communications Biology, 2025, 8 (1), pp.204. 10.1038/s42003-025-07642-7. http://dx.doi.org/10.1038/s42003-025-07642-7

HAL Id: hal-04994861

https://hal.science/hal-04994861v1

Submitted on 22 Aug 2025

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>

# Unravelling the interplay of nitrogen nutrition and the *Botrytis cinerea* pectin lyase BcPNL1 in modulating *Arabidopsis thaliana* susceptibility

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

Check for updates

Antoine Davière<fnref n="1" />,²,⁴, Aline Voxeur<fnref n="1" />,⁴✉, Sylvie Jolivet<fnref n="1" />, Luka Lelas<fnref n="1" />,², Samantha Vernettes<fnref n="1" />, Marie-Christine Soulié<fnref n="1" />,³,⁴ & Mathilde Fagard<fnref n="1" />,⁴

In this study, we investigated the interplay between nitrogen nutrition and the pectin degradation dynamics during Arabidopsis and Botrytis interaction. Our findings revealed that infected detached leaves from nitrogen-sufficient plants released more pectin lyase (PNL)-derived oligogalacturonides compared to nitrogen-deficient ones. We then focused on BcPNL1, the most highly expressed Botrytis PNL upon infection. Using mutant strains lacking BcPNL1, we observed reduced pathogenicity, a delay in germination and a lag in triggering the plant defense response. Additionally, in nitrogen-sufficient detached leaves, the elevated expression of jasmonic acid repressor genes observed upon infection with the wild-type strain was abolished with the mutants. These results linked the increased production of BcPNL-derived products to the increased expression of jasmonic acid repressor genes, contributing partially to the higher susceptibility of nitrogen-sufficient detached leaves. These findings could lay the foundation for new strategies aimed at reconciling both crop resistance to pathogens and the improvement of nitrogen nutrition.

The ascomycete *Botrytis cinerea* is a widely spread necrotrophic fungus responsible for the grey mold disease in over 1000 plant species<fnref n="1" />. As a necrotrophic pathogen, its infectious process relies on several synergistic molecular mechanisms, involving the secretion of a vast array of molecules that include cell wall-degrading enzymes, toxins, reactive oxygen species, peptidases, and acidifying compounds, which have been reported to contribute to its pathogenicity².

The availability of the sequenced genome of this fungus³, has greatly facilitated genetic studies and the identification of putative pathogenicity factors. However, among the virulence factors, some are encoded by multi-gene families such as Carbohydrate-Active enZYmes (CaZY), which are subject to functional redundancy and compensation phenomena that make it challenging to pinpoint the roles of individual B. cinerea genes in pathogenicity. Consequently, only a few genes have been described as having a major role in pathogenicity². Among them, some transcription factors and kinases have been identified, such as sucrose non-fermenting protein kinase 1 (BcSNF1) and the pH regulator BcPACC⁴,⁵. Furthermore, chitin synthases like Bcchs3a are also important pathogenicity factors,

influencing fungal adhesion to the host and plant defense mechanisms⁶. In addition to these, several secreted factors, such as the botrydial and botcinic acid toxins whose production involve the *BcBOT* and *BcBOA* gene clusters, as well as cell death-inducing proteins like BcXIG1 and BcCRH1, play significant roles in pathogenicity⁷⁻⁹. However, recent papers have raised concerns about studying pathogenicity factors exclusively in optimum conditions since plant-pathogen interactions are strongly impacted by abiotic factors¹⁰,¹¹. Indeed for *B. cinerea*, we were able in a previous work to identify undescribed pathogenicity factors by varying plant nitrogen supply which significantly affect the interaction¹².

During infection, *B. cinerea* degrades the cell wall of the plants it infects. To achieve this, the pathogen releases numerous enzymes, with the most abundant belonging to the pectinase family¹³. Pectinases enable the pathogen to penetrate and spread within host tissues by degrading pectins, one of the major components of the plant cell wall barrier. The primary content of pectin, homogalacturonan (HGA), consists of α-1,4-linked residues of D-galacturonic acid (GalA), which can be methyl (Me)- and/or acetyl-esterified (Ac). Four types of HGA with different Degrees of

<footnote n="1">¹Université Paris-Saclay, INRAE, AgroParisTech, Institut Jean-Pierre Bourgin (IJPB), Versailles, France. ²Ecole Doctorale 567 Sciences du Végétal, Univ Paris-Sud, Univ Paris-Saclay, Orsay, France. ³Sorbonne Universités, UFR 927, Paris, France. ⁴These authors contributed equally: Antoine Davière, Aline Voxeur, Marie-Christine Soulié, Mathilde Fagard. ✉e-mail: aline.voxeur@inrae.fr</footnote>
Methylesterification (DM) can be distinguished: the fully methylesterified (DM close to 100%), the highly methylesterified (DM ~ 70%), the low methylesterified (DM ~ 30%) and the fully demethylesterified homogalacturonans also called polygalacturonic acids. The degradation of these different HGA by pectinases leads to the production of oligogalacturonides (OG) with varying degrees of polymerization (DP) and methyl-acetylesterification. Among pectinases, pectin lyases (PNLs, CAZy family PL1), of which five are predicted in *B. cinerea*, degrade fully methylesterified pectins in which leaf epidermal cells are enriched¹⁴. BcPNL1 is the most expressed pectinase during the prepenetration and penetration stages on hard, wax-coated surface¹⁵ and is likely the primary OG-producing pectinases of *B. cinerea* during early infection events fin tomato fruit¹⁶ and *Arabidopsis thaliana* leaves¹⁷ (Supplementary Fig. 1A). PNLs employ β-elimination resulting in highly methyl- and acetylesterified OGs of DP > 4 harboring an unsaturated bond at the non-reducing end¹⁸ (Supplementary Fig. 1C). Concomitantly, fungal esterases, including pectin methylesterases (PMEs, 3 predicted in *B. cinerea*, CAZy family CE8) deesterify both the large methylesterified oligosaccharides and highly methylesterified pectins. Some fungi can also deacetylate pectins using pectin acetylesterases (PAEs) but no PAE is predicted in the *B. cinerea* genome¹⁹,²⁰. Next, depolymerases degrade partially and fully methylesterified pectins and further break down the large PNL-derived oligosaccharides into smaller ones. These pectinases are categorized into two groups: hydrolases with polygalacturonases (Endo-PGs and Exo-PGs, 6 and 3 predicted respectively, CAZy family GH28) which can cleave demethylesterified sites of HGA backbone releasing both low or un-methylesterified OGs (Supplementary Fig. 1C) and lyases with pectate lyases which cleave fully demethylesterified pectins (Endo-PLs and Exo-PLs, 2 predicted, CAZy family PL1) releasing unsaturated non-methylesterified OGs. It is worth to note that BcPG1, which is the most expressed PG upon growth of the lesion beyond the inoculation spot²¹, where pectins are less methylesterified¹⁴, and BcPG2 are crucial pathogenicity factors²²,²³ (Supplementary Fig. 1B). PGs collectively represent 20% of the total protein content in the early secretome of *B. cinerea*¹³.

Plants perceive the oligogalacturonides (OGs) produced by these pectinases as signals of damaged or modified self, enabling them to act as Damage-Associated Molecular Patterns (DAMPs) in Pattern-Triggered Immunity (PTI)²⁴. This signaling pathway triggers several modifications in plant cells, leading to the accumulation of defense proteins and the production of hormones that enhance the defense reaction. Recently, Davidson et al. (2017)²⁵ showed that short OGs (DP3) can efficiently trigger hormonal immunity mediated by jasmonic acid (JA) and salicylic acid (SA). These short OGs resulting from PG activity effectively protected *A. thaliana* from *Pectobacterium carotovorum* infection and are detected in *A. thaliana* leaves infected by *B. cinerea*¹⁷. Interestingly, research has also shown that OG oxidases (OGOX) secreted by the plant itself can inactivate OG-induced defense, thus preventing excessive plant cell death and infection during *B. cinerea* infection²⁶. The hydrogen peroxide ($H_2O_2$) produced during OG oxidation mediates the activation of other defense processes, such as lignification or the inactivation of auxins, through the action of class III plant peroxidases²⁷. Finally, early reports on OG activity demonstrated that long-chain OGs as a result of PG activity on fully demethylesterified HGA (DP10-15) serve as defense elicitors and can induce phytoalexin production²⁸,²⁹. However, we failed to detect these long-chain OGs in *A. thaliana* leaves infected by *B. cinerea*¹⁷ as well as in infected tomato fruit¹⁶. Instead, we detected unsaturated methyl- and acetylesterified OGs with long DP capable of suppressing JA-mediated defense activation by up-regulating the pool of transcripts encoding negative regulators of JA signaling¹⁷. However, the subsequent degradation of these OGs by PMEs and PGs leads notably to the production of unsaturated methyl and acetylesterified tetragalacturonides ($Gal_4MeAc-H_2O$), which, in the end, strongly elicits JA-related plant defense. It is worth to note that, to date, neither the impact of OG methylesterification and DP on their activity is fully clear yet, nor is it known whether the presence of unsaturated bonds affects OG activity.

To gain deeper insight into the strategies employed by *B. cinerea* to successfully infect plants, we investigated the role of BcPNLs in

pathogenicity within the context of environmental modulation of the plant-pathogen interaction, specifically through nitrogen fertilization. Our research provides compelling evidence that nitrogen fertilization influences pectin degradation in *A. thaliana* detached leaves infected by *B. cinerea*. Through mutant and complementation analyses, we demonstrate the crucial role of BcPNL1 in increasing the susceptibility of plants grown with high nitrogen nutrition by inducing the expression of repressors of JA signaling. However, BcPNL1 alone cannot fully explain the connection between nitrogen fertilization and disease susceptibility. Our data confirm the involvement of BcPNL1 in producing unsaturated methyl-acetylesterified OGs with long DP, which are capable of suppressing the JA-mediated defense activation. In contrast, the $\Delta Bcpgnl1$ mutant shows a delay in germination and in the establishment of plant defense response. In high nitrogen nutrition at a late phase of infection with the mutant, there is still a down-regulation of JA-repressor genes which are correlated with the absence of unsaturated methyl-acetylesterified OGs with long DPs. Altogether, our findings reinforce the fundamental role of BcPNL1 in modulating plant susceptibility under high N conditions.

## Results

### Host nitrogen nutrition impacts pectin degradability during *A. thaliana*/*B. cinerea* interaction

First, we assessed if the host nitrogen nutrition impacts the OG turnover (i.e OG production and subsequent degradation into monomers) during infection by analyzing the OG produced in presence of the WT strain B05.10 on detached leaves from plants grown either with 0.5 mM (low N) or 10 mM $NO_3^-$ (high N) which correspond physiologically to nitrogen-deficient or -sufficient regimes respectively¹². We decided to conduct all experiments on living detached leaves. Although this biological model does not reflect actual infection conditions as the disease development has been shown to be slower than on attached leaves³⁰, it offers the possibility to incubate plant leaves in a minimal liquid medium containing fungal spores that facilitates the collection of large amounts of OGs, allowing their structural characterization by a high-resolution mass spectrometry (HRMS)¹⁷. To do so, we incubated plant leaves in a liquid medium containing fungal spores of the WT strain B05.10, and subsequently analyzed the cell wall breakdown products released into the medium using HRMS¹⁷. It is worth noting that the minimal liquid medium used contains 25 mM $NO_3^-$ Consequently, the phenomena we observed correspond to the impact of N nutrition on the plant, notably its cell wall structure, that in turn affects the infection process rather than the impact of N availability on fungal virulence. Since only a few OGs have been shown to be detected at 12 h post-inoculation (hpi), we chose to perform this experiment a bit later, at 15 hpi. Interestingly, we found about 10% more PNL products at high N compared to low N with notably a higher quantity of PNL products of DP6, 7, and 9 at high N compared to low N by about 45%, 26% and 85%, respectively. This resulted in an increased proportion of PNL products and a decreased proportion of PG products, particularly for DP3, 5, and 6 (Fig. 1A). The discrimination between OGs resulting from PNL or PG activity is possible thanks to the specific β-elimination-mediated cleavage of HGA employed by PNLs (Supplementary Fig. 1C).

We next assessed if host nitrogen nutrition could impact fungal pectinase gene expression. We monitored fungal pectinase gene expression at 1-day post-inoculation (dpi), a time point when BcPG1-encoding gene expression has already been shown to be impacted by plant host nitrogen nutrition¹². As depicted in Supplementary Fig. 2, BcPG1,2 and BcPNL1 stand out as the most highly expressed genes within the BcPG and BcPNL gene families, respectively, at both high and low N. BcPNL2, 3, 4 were more than 10-fold less expressed compared to BcPNL1 and BcPNL5 was not expressed at all. Moreover, we observed that at 1 dpi on detached leaves and at 15 hpi in liquid culture with detached leaves, the expressions of BcPNL1 were respectively decreased by 1.5 at high N in infected detached leaves and not significantly different (Supplementary Fig. 2A, B) in liquid culture although we detected more PNL products (Fig. 1A). As a consequence, the variations in PG and PNL products cannot be fully explained by transcriptional regulation of the PG and PNL encoding genes. We then hypothesized that this
<figure><img src="image_4.png" /><figcaption>**Fig. 1 | Impact of host nitrate nutrition on Botrytis cinerea pectinolytic capacities and plant defense regulation.** A LC-MS quantitative analysis of OG production after incubation of spores from the WT strain of the fungus for 15 h with *A. thaliana* Col-0 leaves from plants cultivated at high or low N. The graph represents the peak area of each type of product over the total of all detected OGs expressed as %. B Transcript accumulation during infection in low N and high N conditions at 1-day post-inoculation on detached leaves are represented for plant defense genes. Plant</figcaption></figure>

gene targets were normalized with *AtUBI4* and similar results were obtained using *AtAPT1*. Three independent experiments were conducted with similar results. Data are expressed as mean normalized expression in arbitrary units (a.u.) and are the means of triplicates (±SE). Statistical differences represented between high and low N are the results of an ANOVA followed by a Fisher's LSD test in (A) and a two-sample t-test in (B): \*p < 0.05, \*\*p < 0.01, \*\*\*p < 0.001, \*\*\*\*p < 0.0001.

increased PNL activity in high N conditions could be caused by differences in pectin chemistry. One important feature for pectin degradation is its esterification status, so we tested whether PNL activity could be affected by methyl- and acylesterification levels. Interestingly, we could only detect PNL products when the fungus was grown *in vitro* with methyl-acylesterified pectins and not with methylesterified pectins suggesting that the presence of acylesterification is a key factor modulating PNL activity (Supplementary Fig. 3A, B). Then, in a time course experiment on

detached leaves from 12 hpi to 18 hpi, we observed that methyl-acylesterified PNL products and methylesterified PG products were respectively increased and reduced in high N conditions compared to low N whereas unesterified and acylesterified products were similar (Supplementary Fig. 3C). Altogether, these results suggest that the differences in PNLs mediated degradation upon fungal infection between high and low N detached leaves are due to differences in the organisation and/or abundance of methyl-acylesterified regions in HGA.
<figure><img src="image_5.png" /><figcaption>**Fig. 2 |** Targeted gene replacement of *BcPNL1* in *Botrytis cinerea*. Representation of the deletion strategy used to obtain *BcPNL1* mutants (A) and complemented strains (B). Arrows indicate primers P1–P6 used for diagnostic PCR and primers F1 and F2 used to obtain the *fenR* probe to confirm absence of ectopic integration by Southern blot (Supplementary Fig. 4). C Transcript accumulation of the *BcPNL1* gene at 1 dpi</figcaption></figure>

on *A. thaliana* Col-0 leaves (high N) with the different strains of *B. cinerea*. *BcActA* was used for normalization of *BcPNL1* expression. Data are expressed as mean normalized expression in arbitrary units (a.u.) and are the means of triplicates (±SE) Letters represent the results of a two-sample t-test, *p* < 0.05.

Since PNL products of DP > 4 have been shown to up-regulate *JOX* and *JAZ* genes¹⁷ respectively involved in repression of the JA pathway by inactivation of JA through hydroxylation³¹,³² and repression of JA-induced genes³³, we searched for differences in nitrogen-dependent *JAZ* and *JOX* regulation during infection. Since it has been shown that the majority of changes in plant gene expression occurs by 24 hpi when the pathogen has penetrated the leaf epidermis, and that, naturally, OG production precedes transcriptomic activation³⁴, we chose to delay the transcriptomic analysis slightly compared to the OG analysis and performed it at 1 dpi. We found a significant higher expression of *AtJAZ3,4* and *AtJOX3* in high N compared to low N at 1 dpi with increases of 1.7-, 2.0-, and 1.7-fold, respectively (Fig. 1B). This was consistent with the higher production of high DP PNL products upon high N detached leaves infection (Fig. 1B). Accordingly, the JA marker gene *AtPDF1.2* was less expressed in high N conditions which was also previously observed in Soulie et al.¹². Since SA and JA defense signalling pathways are mutually antagonistic, we also searched for SA marker genes³⁵. Upon infection, we observed a 6-fold increase in *AtPR1* expression at high N and a 2-fold increase in the transcription factor *AtWRKY70* which was shown to be activated by SA and repressed by JA although its regulation is more complex³⁶,³⁷.

## Targeted gene deletion of the *BcPNL1* gene in *B. cinerea* leads to reduced pathogenicity at high N and low N

Next, to draw a clear connection between the higher accumulation of PNL products and the downregulation of the JA defense signaling pathway in high N compared to low N detached leaves, we investigated the function of *BcPNL1* which stands out as the most highly expressed gene within the *BcPNL* gene family as depicted in Supplementary Fig. 2. To investigate its function, we replaced the *BcPNL1* coding region with a deletion cassette containing the fenhexamid resistance gene, which was amplified

from pTEL-fenh³⁸. We employed PEG/CaCl₂-mediated transformation of *B. cinerea* protoplasts, in conjunction with the recently described CRISPR-Cas9 system, utilizing the ribonucleoprotein complex direct delivery strategy (Fig. 2A)³⁸. *In vitro*, we selected fenhexamid-resistant transformants.

The absence of the wild-type *BcPNL1* gene and the correct integration of the cassette in each isolate were determined through PCR and confirmed via Southern blot using a DIG-labeled probe. Supplementary Fig. 4A demonstrates that out of 20 resistant transformants, we successfully identified three homokaryotic mutant strains, designated *ΔBcPNL1.1*, *ΔBcPNL1.2*, and *ΔBcPNL1.3*. Others, such as the *ΔBcPNL1.4* transformant, were found to be heterokaryotic. Nevertheless, all tested strains were devoid of ectopic integration, confirming the precision of the CRISPR-Cas9 system (as shown in Supplementary Fig. 4B). Subsequently, we selected the mutant strains *ΔBcPNL1.1*, *ΔBcPNL1.2*, and *ΔBcPNL1.3* for further phenotypic characterization.

We next complemented the *ΔBcPNL1.1* mutant with the *BcPNL1* gene either under the control of the native promoter (1.5 kb upstream of the ORF) or the strong *olicC* promoter from *Aspergillus nidulans*, resulting in the *ΔBcPNL1.1c1* and the *ΔBcPNL1.1c2* strains respectively (Fig. 2B and Supplementary Fig. 4C). The complementation cassettes were integrated at the same intergenic region of chromosome 3 using CRISPR-Cas9³⁹. RT-qPCR analysis of the complemented lines revealed that the level of *BcPNL1* expression was only partially restored compared to the B05.10 strain (Fig. 2C).

First, we investigated whether the deletion of *BcPNL1* would affect the growth of the fungus *in vitro* when cultivated on solid medium in the presence of different types of substrates. Although we observed slight differences between the mutant strains, we observed no significant impact of *BcPNL1* deletion either on radial growth with Czapeck minimal medium supplemented with glucose or citrus highly methylesterified pectin as the sole carbon sources or on *in vitro* spore production on PDA (Supplementary Fig. 5).
<figure><img src="image_6.png" /><figcaption>**Fig. 3** | The Δ*Bcpn1l* mutant strains showed a strongly reduced pathogenicity on detached leaves that is partially restored in the complemented strains. A Col-0 *A. thaliana* detached leaves at 2dpi. Scale bar = 1 cm. B Plot of the lesions at 2 and 3 dpi on Col-0 *A. thaliana* detached leaves. Plants were either grown with high or low nitrate supply (low N, high N). C Plot of the lesions on *A. thaliana* detached leaves at 2 and 3dpi in high N with B05.10, the Δ*Bcpn1l* mutant and the two complemented strains Δ*Bcpn1l*1c1 and Δ*Bcpn1l*1c2. D Coton blue coloration of the fungus on onion epidermis. Onion epidermises were inoculated with spore droplets of either B05.10 or the Δ*Bcpn1l* mutant strains in ½ PDB and incubated for 8 h (top panel) or 24 h (lower panel) before coloration and observation under the microscope. Red arrows show penetrating hyphae. Since results were the same with Δ*Bcpn1l*1/2/3,</figcaption></figure>

only the observations for Δ*Bcpn1l*1 are presented. Scale bar = 200 µm. E *S. lyco-
persicum* detached leaves, 2 dpi. F Plot of the lesions on tomato detached leaves at 2
and 3 dpi. Tomato plants were grown in high N (10 mM NO₃⁻). Scale bar = 1 cm.
G Plot of the lesions at 2 and 3 dpi on Wassilewskija *A. thaliana* detached leaves.
H Plot of the lesions at 2 and 3 dpi on Col-0 *A. thaliana* detached leaves in high N
condition with B05.10, the Δ*Bcpn1l*1 mutant and two other pectinase mutants of the
fungus Δ*Bcpg1* and Δ*Bcpme1/2*.
Letters represent the results of an ANOVA test
followed by a Fisher's LSD for multiple comparison (*p* < 0.05). Lesion sizes are
expressed in cm² or normalized between different experiments by the largest lesion
and expressed as % of the largest lesion, each boxplot represents at least 15 different
leaves.

Next, we assessed whether the deletion of *BcPNL1* could affect pathogenicity on high and low N detached leaves. The *ΔBcpn1l* mutants exhibited a substantial reduction in pathogenicity on Columbia (Col-0) under both nitrogen nutrition conditions after 2 and 3 dpi compared to B05.10 infection (Fig. 3A, B). qPCR quantification of fungal DNA during the infection confirmed the reduced pathogenicity of these mutants (Supplementary Fig. 6). Furthermore,

we observed an increased susceptibility of high N detache,d leaves
infected either by B05.10 or Δ*Bcpn1l* mutants (Fig. 3B) suggesting
that BcPNL1 pectinolytic activity is not sufficient to explain the
impact of host nitrogen nutrition on disease severity. Moreover, we
observed a partial restoration of the pathogenicity with both com-
plemented strains on Col-0 *A. thaliana* detached leaves (Fig. 3C)
which is consistent with their *BcPNL1* expression level (Fig. 3C). To
characterize the origin of the pathogenicity defect, we analyzed the germination and penetration abilities of the Δ*Bcpn11* mutants on onion epidermis (Fig. 3D). Interestingly, the Δ*Bcpn11* mutants exhibit a delay in spore germination compared to the wild-type B05.10 strain, although penetration was still observed at the 24 hpi time point as for the wild-type.

We then confirmed that the role of *BcPNL1* in pathogenicity is not restricted to Col-0 *A. thaliana* infection since detached leaves of tomato plants and Wassilewskija (Ws) *A. thaliana* plants infected with the *ΔBcpn11* mutants also presented smaller lesions compared to B05.10 (Fig. 3E–G). Lastly, we compared the pathogenicity on Col-0 detached leaves of *ΔBcpn11*, *ΔBcpg1* and *ΔBcpme1ΔBcpme2* and observed that the impact of *BcPNL1* deletion was more dramatic than that of *BcPG1* (Fig. 3H), which is consistent with the fact that the digestion of highly methyl-acetylesterified pectins by *BcPNL1* occurs earlier (i.e pre- and penetration stage) to *BcPG1* activity, of which the corresponding gene reaches maximum expression upon lesion spreading²¹. In contrary to the increased pathogenicity observed previously with the *ΔBcpme1ΔBcpme2* mutant on Ws detached leaves¹⁷, we found here a WT-like pathogenicity level on Col-0 as observed previously on grapevine and tomato leaves⁴⁰ indicating some host specificity in the role of *B. cinerea* pectin methylesterases.

### BcPNL1 is involved in the reduced cell death and oxidative burst responses during A. thaliana/ B. cinerea interaction

To characterize the plant response to the hypovirulent Δ*Bcpn11* mutants, we analyzed the oxidative burst by measuring H₂O₂ production using a 3,3'-diaminobenzidine (DAB) staining at an early time point after inoculation. At 1 dpi, high and low N Col-0 leaves infected with B05.10 strain exhibited a strong DAB staining at the inoculation site and its vicinity, while leaves infected with the Δ*Bcpn11* mutants showed only slight staining at the inoculation site (Fig. 4). We also assessed plant cell death upon infection using trypan blue (TB) staining. At the same time point (1 dpi), TB staining revealed fungal growth as well as small, darker spots indicating plant cell death, but this was observed only for the B05.10 strain at high and low N (Fig. 4B, C). These results indicate that at this early time point, the plant had strongly responded to the pathogenic B05.10 strain and lesser to the mutants.

### ΔBcpn11 mutants produce less PNL and PG products at high N and low N

Next, we investigated the impact of *BcPNL1* deletion on cell wall degradation. The analysis of pectin degradation products in the presence of the Δ*Bcpn11* mutants at 15 hpi revealed on one hand a drastically reduced production of PNL-derived OGs compared to the WT strain, with both high N and low N detached leaves (Fig. 5A), confirming the importance of BcPNL1 on the total fungal PNL activity. On the second hand, PG products were also strongly reduced upon infection with the mutants, indicating that the PNL activity is essential upon infection to regulate PG and PME activity. Next, we performed an OG analysis on high N or low N detached leaves infected by Δ*Bcpn11.1c1* and the Δ*Bcpn11.1c2* strains at 18 hpi. We chose 18 hpi in order to assess in the meantime if any compensation mechanism would have been triggered in Δ*Bcpn11.1* mutant at a later time point. At this later time point than our previous analysis (Fig. 1A), high DP OGs couldn't be detected but we could still observe that the production of OGs was partially restored with the complemented strains compared to the Δ*Bcpn11.1* mutant (Fig. 5B). No clear conclusion could be formulated regarding PNL products for which only small DP were detected suggesting that, at this time point, most of the PNL products were already degraded. However, we observed with the Δ*Bcpn11.1c* strains the restoration of PG products accumulation for DP3 and DP4 in both N conditions and DP5 in low N that confirms that PNL activity is essential prior to PG and PME activity.

In order to assess if BcPNL1 activity could regulate the other pectinases at the transcriptional level or if the decrease of PG products only results from

the inability of PMEs or PGs to have access to their substrate without prior PNL activity, we assessed the level of expression of pectinase genes in Δ*Bcpn11* mutants at 1 dpi (Fig. 5C). Interestingly, at 1dpi all the tested pectinase genes were downregulated in the Δ*Bcpn11* mutants compared to the WT strain in both N conditions. At 2 dpi, although less pronounced, we could still observe downregulation of the three *BcPNLs*, *BcPG4*, *BcPG5*, and *BcPG6* in the mutants whereas it was more variable for the *BcPMEs* and *BcPG1*. Interestingly, *BcPG2* and more clearly *BcPG3* were upregulated in the Δ*Bcpn11* mutants suggesting a delayed compensation mechanism to degrade pectin.

<figure><img src="image_7.png" /><figcaption>**Fig. 4 |** Cell death and oxidative burst reaction during infection with WT and mutant *B. cinerea* strains. Detached leaves of *A. thaliana* were infected with drops either only with ½ PDB (Mock) or drops containing spores of WT or Δ*Bcpn11* mutant strains of *B. cinerea*. A One day after inoculation, treated leaves were cut in half and each half was immediately colored with Trypan Bleu (TB) for cell death visualisation (on the left of each image) or with 3,3'-diaminobenzidine (DAB) for oxidative burst visualisation (on the right in each image). Scale bars = 500 µm. Zoom of a WT (B) or mutant (C) infected leave colored with TB. Black arrows show the fungus and the red arrow shows a spot of induced necrosis. Scale bar = 50 µm. Images were taken with a ZEISS Axio Zoom V16.</figcaption></figure>

Taken together, these results suggest a lower aggressiveness of the Δ*Bcpn11* mutants compared to B05.10 at high and low N, consistent with the findings from the pathogenicity tests, the cell death and oxidative burst responses and the OGs analysis.

### The impact of host nitrogen nutrition on plant JA-related gene expression is partially altered during infection with the Δ*Bcpn11* mutants

Next, we investigated the expression of plant defense-related genes during the infection (Fig. 6). At 1 dpi, the expression of plant defense-related genes, such as the JA marker gene *AtPDF1.2a*, the chitinase *AtPR4*, and *AtPAD3* coding for the final enzyme involved in camalexin biosynthesis⁴¹, was
<figure><img src="image_8.png" /><figcaption>**Fig. 5** | Impact of *BcPNL1* mutation on *B. cinerea* pectinolytic capacities. A LC-MS analysis of OG production from leaves of *A. thaliana* plants (cultivated at high or low N) incubated for 15 h in liquid medium supplemented with spores of the wild-type and mutant strains of *B. cinerea*. A two-way ANOVA test was performed followed by a Fisher's LSD test between wild type and mutant strains: \*p < 0.05, \*\*p < 0.01, \*\*\*p < 0.001. B LC-MS analysis of OG production from leaves of *A. thaliana* plants (cultivated at high or low N) incubated for 18 h in liquid medium supplemented with spores of the WT, mutant and complemented strains of *B. cinerea*. Letters represent the results of a two-way ANOVA test followed by a Fisher's</figcaption></figure>

LSD test between strain for each DP group (p > 0.05). DP: degree of polymerisation. C Heatmap representing the log2 fold change of fungal pectinase transcript accumulation between the B05.10 reference strain and the Δ*Bcpnl1* mutants during infection in high and low N conditions at 1 dpi and 2 dpi on detached leaves. Fungal gene targets were normalized with *BcActA* and similar results were obtained using *BcUBI* for normalization. Three independent experiments were conducted with similar results. Data are means of triplicates (±SE). Differences that were not significant according to an unpaired t-test are represented in gray: *p* < 0.05.
<figure><img src="image_9.png" /><figcaption>**Fig. 6** | Differences in plant defense gene expression during infection of detached leaves from high and low N *A. thaliana* plants with the B05.10 reference strain and the *ΔBcpn1* mutant strains. Transcript accumulation of plant defense related genes obtained by RT-qPCR from detached leaves after 1 and 2 dpi with the WT and ΔBcpn1 mutant strains in high N and low N conditions. Letters are the result of an</figcaption></figure>

ANOVA followed by a fisher LSD test (*p* < 0.05). Plant gene targets were normalized with *AtUBI4*. Data are expressed as mean normalized expression in arbitrary units (a.u.) and are the means of triplicates (±SE). Similar results were obtained using *AtAPT1* for normalization. Three independent experiments were conducted with similar results.
decreased in detached leaves infected with the ΔBcpn1l mutants compared to the WT strain. However, at 2 dpi, transcript levels in detached leaves infected with the ΔBcpn1l mutants increased in all cases compared to 1 dpi, particularly under high N conditions. In this scenario, *AtPDF1.2a* expression was notably more elevated upon infection with the ΔBcpn1l mutants than with the WT. These initial data suggest that the activation of defense mechanisms, particularly through JA production, is delayed in presence of ΔBcpn1l mutants and modulated by nitrate content.

We also observed that the expression of *AtJOX3* and the transcription factors *AtJAZ1*, 3, and 4 was down-regulated in detached leaves infected with the mutants compared to the WT at 1 dpi under all conditions (Fig. 6). Interestingly, this down-regulation persisted at 2 dpi, but only under high N conditions. On the other hand, the expression of the *AtPR1* gene was very low at 1 dpi across all conditions, but an increase was observed particularly with the WT strain and under high N concentration. With the mutants, however, this response was diminished under high N, suggesting a decrease in antimicrobial activity in these conditions⁴².

The expression profile of the *AtPR2* gene, coding for a β-1,3-glucanase with antifungal properties also involved in callose deposition and SA-dependent defense responses⁴³,⁴⁴, indicates that the presence of ΔBcpn1l mutants on detached leaves triggered a plant response linked to SA at 1 dpi in high N which is absent in mock or upon inoculation with the WT. At 2 dpi, *AtPR2* expression was induced in the high N mock compared to 1 dpi, which may be an effect of using detached leaves. By contrast, *AtPR2* expression was downregulated upon infection compared to the mock in a similar way with the WT and the mutants. (Fig. 6). Finally, the expression of the *AtWRY70* gene, which is involved in SA and JA/ethylene signaling pathways³⁶,³⁷, was repressed upon infection compared to the mock at 1 dpi suggesting the presence of JA. This gene repression was notably higher with the WT compared to the ΔBcpn1l mutants in both N regimes. At 2 dpi, *AtWRY70* expression was still repressed by the infection but interestingly the repression was higher with ΔBcpn1l mutants compared to the WT under high N conditions which might reflect a delay in *AtWRKY70* repression with ΔBcpn1l mutants. This indicates that SA and JA/ET signaling is modified during the kinetics of *B. cinerea* infection to tend to produce JA and this process is influenced by N levels and the absence of *BcPNL1*.

All these differences in gene expression profiles in presence of the WT strain or the mutants suggest a delay in the establishment of the plant's defense mechanisms in the presence of ΔBcpn1l mutants, which is more pronounced under high N conditions. This weak repression of the JA biosynthesis pathway is correlated to the absence of high DP unsaturated OGs in ΔBcpn1l mutants under high N conditions as proposed in the model presented in Fig. 7.

## Discussion

*Botrytis cinerea*, like all pectin-degrading fungal pathogens, relies on its pectinases to break the cell wall barrier and release nutrients that support fungal growth and development. However, the impact of abiotic factors such as nitrogen nutrition on this cell wall degradation process, as well as the role of PNLs as pathogenicity factors in *B. cinerea* pathogenesis, had never been reported. In this study, we were able to point towards the involvement of BcPNL1, encoded by the most expressed PNL genes in the early phase of infection (Fig. 1 and Supplementary Fig. 1) and from which the *in vitro* PNL activity was previously confirmed¹⁷ in orchestrating part of the interaction between host nitrate nutrition and disease severity.

Thanks to the recently described CRISPR-Cas9-mediated gene replacement in *B. cinerea*³⁸, we have obtained and then characterized ΔBcpn1l mutants (Fig. 2). *In vitro*, the mutants presented normal growth with glucose or highly methylesterified pectin as the sole carbon source (Supplementary Fig. 5). Similarly, Δpnl1 mutants of *Penicillium digitatum* displayed no growth defects *in vitro*¹⁵. However, it would be interesting to test the growth of the mutants with more or less acetylated pectins given the preference of BcPNL1 for methyl-acetylated pectins (Supplementary Fig. 3). Interestingly, on onion epidermis we observed a delayed spore germination

of ΔBcpn1l mutants (Figs. 3D, 7A). In the transcriptome analysis of germinating spores on a plant mimicking surface, Leroch et al. (2013)¹⁵ showed that BcPNL1 expression is strongly induced after only one hour of incubation in agreement with our finding that BcPNL1 is the most expressed PNL at 1 dpi (Supplementary Figs. 1 and 2). These observations suggest that BcPNL1 might be implicated in early events promoting germination upon perception of the plant surface.

In other phytopathogens such as *Pectobacterium carotovorum*, *P. digitatum* and *Verticillium dahlia*, PNL activities represent important determinants of pathogenicity⁴⁵⁻⁴⁷. Concerning *B. cinerea*, to our knowledge no study had clearly investigated the role of BcPNL genes. On detached leaves, we observed a reduced pathogenicity of ΔBcpn1l mutants restored in the complemented strains (Fig. 3) demonstrating that PNL activity is an important determinant of *B. cinerea* pathogenicity too. Furthermore, the reduced pathogenicity of the ΔBcpn1l mutants on different host plants also suggests that this pathogenicity factor is not linked to host specificity. Supporting this hypothesis, Cotoras & Silva (2005) observed in WT *B. cinerea* strains isolated from tomato and grape a similar level of secretion of PNLs. In other fungal pathogens such as *V. dahlia* and *P. digitatum*, mutants of PNLs were also hypovirulent on their respective hosts⁴⁵,⁴⁶,⁴⁸,⁴⁹. Authors not only mention their role in promoting penetration through cell wall degradation but also their necrotizing activity. In our study, according to trypan blue colorations, initial necrosis of the plant cells seemed to be reduced in presence of the ΔBcpn1l mutants (Fig. 4). Movahedi & Heale (1990) showed that a PNL purified from *B. cinerea* germinating spores had a cell death activity on carrot cell culture and was able to produce elicitors from cell wall components. Since BcPNL1 was the only PNL gene that was strongly expressed in the early interaction, it might be the same enzyme that the authors purified and characterized as a necrosis-inducing protein. Finally, we observed that the increased susceptibility of high N plants to *B. cinerea* was still observed with ΔBcpn1l mutants suggesting that BcPNL1 pectinolytic activity alone is not sufficient to explain the impact of host nitrogen nutrition on disease severity.

We also showed that more PNL products were produced in high N plants infected with the B05.10 wild-type strain of *B. cinerea*, and this higher production correlates with a higher up-regulation of *AtJOX3* and *AtJAZs* (Fig. 1). These genes are involved in repression of the JA pathway which is strongly implicated in efficient defense strategy against *B. cinerea*²⁰,⁵¹ and are induced by long methyl esterified PNL products, accumulated upon plant infection with the ΔBcpme1ΔBcpme2 mutant, hypervirulent on *A. thaliana*. With the ΔBcpme1ΔBcpme2 mutant, these high DP (>4) PNL products, which cannot be further degraded by PGs induce an up-regulation of genes involved in JA pathway repression such as *AtJOX3* and *AtJAZs*¹⁷. Therefore, to explain why high DP PNL products were accumulated in high N infected plants, we assessed if the host nitrogen nutrition could impact the transcriptional regulation of BcPNLs and/or BcPMEs. However, we did not observe any significant up-regulation of these genes in high N condition that would have explained the accumulation of high DP PNLs OGs (Supplementary Fig. 2). Interestingly we found that pectin acetylesterification is required for PNL activity and that more methyl-acetylesterified OGs were produced upon inoculation of high N detached leaves in liquid compared to low N (Supplementary Fig. 3). Consequently, although we cannot rule out the possibility of post-transcriptional or post-translational regulation of PNLs, our results suggest that the increased release of high DP PNLs OGs in the high N condition is due to a higher level of methyl-acetylesterified regions allowing PNL activity as proposed in the model (Fig. 7). It was shown in *A. thaliana* that an acetyl-transferase mutant and a transgenic line producing fungal PAEs with reduced pectin acetylation level are more resistant to *B. cinerea* infection²⁰,⁵². Thereby, a control of PNL activity by the pectin acetylation level could contribute to its impact on disease severity.

Consistently, in the ΔBcpn1l mutants, pectin degradation capacities were strongly affected both at high and low N. Far less PNL-derived cell wall
<figure><img src="image_10.png" /><figcaption>**Fig. 7** | **Proposed model placing BcPNL1 in the interplay between plant nitrogen nutrition and disease susceptibility.** A The pectin lyase BcPNL1 is involved in early germination events on the plant surface. B, C Plant nitrogen nutrition affects early degradation of the plant cell wall by BcPNL1 probably because of differences in pectin acetylation levels which is a key factor for BcPNL1 activity. In this model, the</figcaption></figure>

stronger acetylation level of pectins in high N leads to an increased BcPNL1 activity and the accumulation of high DP methyl-acetylated and unsaturated OGs. This correlates on the plant side with earlier and stronger activation of JA repressor genes in high N, absent upon infection with the *ΔBcpnl1* mutants confirming the role of BcPNL1-derived OGs in inhibiting plant defenses through unknown receptors.

products were detected upon infection with the *ΔBcpnl1* mutants (Fig. 5) which corroborates the absence of functional compensation observed with the other pectin lyases at the transcriptional level (Fig. 5C). Only *BcPG3* was strongly upregulated in the mutants at 2 dpi and *BcPG2* to a lesser extent although they were downregulated at 1 dpi as all the tested pectinase genes and probably because of a delayed development of the fungus. Interestingly, *BcPG3* was shown to have a better activity on low methylesterified pectin (DM 7%) compared to polygalacturonic acid²² which can partly explain why it is compensating the *BcPNL1* deletion at the transcriptional level and not the others. Indeed, it is known that pectinases are regulated at the transcriptional level by pectin and pectic components¹⁸, thus, *BcPNL1* products released upon degradation might induce genes coding for enzymes able to achieve the next degradation step. Overall, deletion of *BcPNL1* seemed to modify deeply the kinetics of pectin degradation and OG turnover during infection which can explain the delayed pathogenicity of the mutants at high and low N.

During host-pathogen interactions, OGs produced by pectinases are important DAMPs which modulate immunity in plants¹⁷,²⁴,²⁵. In this study, the pectin lyase mutant *ΔBcpnl1* impaired in pectin degradation presented a germination delay associated with a delay in plant defense induction

(Figs. 3D, 6, 7). Interestingly, at a later stage (2 dpi) the plant defense response, particularly the JA defense pathway, was stronger with the mutant pointing out its incapacity to repress early JA induction (Figs. 6, 7). This delayed hyper-induction of plant defenses caused by BcPNL1 deletion was exacerbated in high N conditions where the JA pathway repression is stronger in the WT and the BcPNL1 derived OGs more accumulated. These data underline on one hand, that high DP unsaturated OG produced by BcPNL1 are capable of suppressing plant JA-mediated defense activation and on the other hand that these OGs have a role in the nitrogen-dependent down-regulation of plant defense. In addition, using a large collection of fungal isolates, it was found that the level of *BcPNL1* expression is positively correlated to the disease symptoms on *A. thaliana* Col-0 but not with the *npr1* and *coi-1* mutant lines affected in the SA and JA pathway respectively⁵³ reinforcing our finding that BcPNL1 is involved in plant defense modulation. It would be of interest to investigate the exact mechanism by which those specific BcPNL1-derived OGs are perceived and trigger the downstream plant defense repression. A recent study showed that the previously designated OG receptors that belong to the wall-associated kinase family (WAK, 5 members) are in fact dispensable for the OG-mediated immune response⁵⁴ indicating that there is still much to discover on OG signalling.
To conclude, we were able in this work to describe the pectin lyase BcPNL1 which, to the best of our knowledge, is an unreported pathogenicity factor of *B. cinerea* crucial for pectin degradation. Its activity appears to be affected by pectin acetylation changes induced by host nitrogen nutrition resulting in a greater production of high DP PNL-derived OGs at high N which partially hijack the DAMP-mediated immunity.

## Methods

### Fungal and plant material

*B. cinerea* strain B05.10 collected from *Vitis* in Germany⁵⁵ was used throughout this work as the WT reference strain. Strains of *B. cinerea* were routinely cultured on Potato dextrose agar (PDA, Difco) or Malt Agar (1% malt extract, 1.2% agar) at 23 °C under constant light. For fungal growth rate measurements, we used Czapeck minimal media containing 2.55 g.l⁻¹ NaNO₃, 0.5 g.l⁻¹ KCl, 0.5 g.l⁻¹ MgSO₄·7H₂O, 10 mg.l⁻¹ FeSO₄·7H₂O, 1 g.l⁻¹ K₂HPO₄, 10 mg.l⁻¹ Na₂MoO₄·2H₂O and either 30 g.l⁻¹ glucose or 10 g.l⁻¹ of highly methylesterified pectin from citrus peel (DM > 70%, P9561, Sigma-Aldrich) as the sole carbon source. *A. thaliana* ecotypes Col-0 and Ws-0 were obtained from the INRAE-Versailles collection. For all experiments, plants were cultured in growth chambers (65% relative humidity, 8 h photoperiod, 21 °C) for 6 weeks on nonsterile sand either supplied with a nutrient solution containing 0.5 mM NO₃⁻ (Low N) or 10 mM (High N). *Solanum lycopersicum* plants (Monalbo) were cultured in the same conditions for 6 weeks but were only supplied with 10 mM NO₃⁻ nutritive solution.

### Construction of the ΔBcpln1 mutants and complemented strains

ΔBcpln1 mutants and complemented strains were obtained by using CRISPR-Cas9 mediated homologous recombination (HR) (Fig. 2)³⁸. Briefly, fresh spores of the B05.10 strain (1.10⁸ conidia) were cultured overnight in 100 ml of HA liquid medium (1% malt extract, 0.4% glucose, 0.4% yeast extract, pH 5.5). The obtained fungal biomass was washed twice in KCl/NaPi buffer (0.6 M KCl, 100 mM sodium phosphate pH 5.8) and then incubated in KCl/NaPi buffer with 1.35% (w/v) protoplasting enzymes (Glucanex) 1 h at 28°C on a rotary shaker at 60 rpm. Protoplast were then washed twice with TMS buffer (1 M sorbitol, 10 mM MOPS, pH 6.3) and concentrated in TMSC buffer (TMS + 50 mM CaCl₂) at the desired concentration (2.10⁷ conidia.100 µl⁻¹). PEG-mediated transformation was achieved by addition of PEG solution (0.6 g ml⁻¹ PEG 3350, 1 M sorbitol, 10 mM MOPS, pH 6.3) to the protoplasts previously incubated for 10 min on ice with the ribonucloprotein (RNP) complex and the DNA repair template (10 µg). The RNP complex was obtained by incubating for 30 min at 37 °C in 15 µl of EngesCas9 specific buffer (NEB), 6 µg of the commercial EngenCas9 (NEB) with 2 µg of the freshly synthesized SgRNA1 obtained with the HiScribe™ T7 High Yield RNA Synthesis Kit (NEB). DNA repair templates containing the fenhexamid resistance cassette surrounded by 60 bp of the target gene for HR were obtained by PCR. For the ΔBcpln1 mutants, primer X1 and X2 were used on the pTEL-fenh plasmid kindly provided by Pr. Matthias Hahn (Kaiserslautern, Germany). For the complemented strain c1, primers X3 and X4 were used to amplify the coding sequence of *PNL1* together with the native promoter (1.5 kb) and the terminator (0.5 kb). For the complemented strains c2 the oliC promoter and the gluc terminator were amplified from the pndh-ogg plasmid⁵⁶, kindly provided by Muriel Viaud (INRAE, Saclay), using the primer pairs X5-X6 and X7-X8 respectively whereas the *PNL1* coding sequence was amplified with X9 and X10. To amplify the hygr cassette, the X11 and X12 primers were used on the pndh-ogg. Finally, the two or four fragments respectively were then assembled by overlapping tail PCR using X13 and X14 to form the final fragment carrying the 60 bp for HR in an intergenic region of chromosome 3³⁹. The protoplasts were then plated on SH agar medium (0.6 M sucrose, 5 mM Tris-HCl pH 6.5, 1 mM (NH₄)₂H₂PO₄, 9 g.l⁻¹ bacto agar) supplemented with 1 µg.ml⁻¹ fenhexamid or 40 µg.ml⁻¹ hygromycin B for regeneration. Resistant transformants were then subcultured with the selection pressure on HA agar (HA, 1.5% agar) and once again on PDA and MA for stock solution of spores and DNA extraction respectively. Finally,

correct homokaryotic integration was verified by PCR (Supplementary Fig. 4A and C) on genomic DNA extracted⁵⁷. Unique integration of the cassette in the mutants was then controlled by Southern blot (Supplementary Fig. 4B) using DIG-High Prime DNA Labeling and Detection Starter Kit I (Sigma). Primers used for these experiments are listed in Supplementary Table 1. Selection of the SgRNA protospacer sequences was done with the online CRISPOR tool⁵⁸.

### Detached leaves infection experiments

For all infection experiments, spores of each strain were harvested from 2 weeks-old PDA cultures. For pathogenicity assays, RT-qPCR experiments and colorations, detached leaves were inoculated on their upper surface with 20 µl drops containing 5.10⁵ conidia. ml⁻¹ in ½ Potato dextrose broth (PDB, Difco). Leaves were then kept in a Petri dish under high humidity conditions and daylight at 23°C. Fungal DNA quantification in infected leaves was performed according to⁵⁹.

### DAB and BT coloration experiments

For coloration assays, leaves were infected as described above and cut in half 24 h after inoculation prior to coloration. Detection of H₂O₂ production was achieved by 3,3'-diaminobenzidine (DAB) staining. Briefly, samples were immersed in an aqueous solution containing 1 mg.ml⁻¹ DAB (pH 3.6, Sigma) and incubated overnight at RT on a rotary shaker. To discolor the leaves, chlorophyll was extracted by incubation in 75% ethanol. To detect cell death, Lactophenol Trypan Blue (TB, Sigma) was performed⁶⁰. Pictures were taken on a ZEISS Axio Zoom V16.

### Gene expression analysis

For RT-qPCR, in each condition total RNA was extracted from a pool of 8 leaves after 1 and 2 dpi and extracted using TRIzol reagent (Invitrogen). Three independent biological replicates were analyzed. Reverse transcription was performed using an oligo-dT20 for a primer and Superscript II RNaseH-reverse transcriptase (Invitrogen). PCR was performed on a Applied Biosystems™ QuantStudio™ 5 Real-Time PCR System with SYBR Green PCR MasterMix (Eurogentec). Each reaction was performed on a 1:20 dilution of the cDNA, synthesized as described above, in a total reaction of 5 µl. Gene expression values were normalized to expression of the *A. thaliana* *AtAPT1* and *AtUBI4* genes or *B. cinerea* *BcACTA* and *BcUBI*⁶¹. We obtained similar results with both genes, therefore results with only one reference gene are shown. Specific primer sets are given in Supplementary Table 2.

### Analysis of OG production

For OG production experiments, triplicates of 10 detached leaves were directly immersed in 35 mL of a *B. cinerea* spores suspension containing 3.10⁵ conidia.ml⁻¹ in Gamborg's B5 basal medium, 2% fructose and 10 mM phosphate buffer pH 6.4¹⁷. Samples were incubated at 23 °C on a rotary shaker at 80 rpm. After 12, 15, or 18 h of incubation of the leaves in the suspension of spores an equal volume of ethanol 96% was added in order to precipitate the largest molecules. For assays performed on pectins, GAT103 and GAT300 were used as methyl- and methyl/acylesterified pectins, respectively, https://www.elicityl-oligotech.com/ and incubated overnight in Gamborg medium containing *B. cinerea* spores at 3.10⁵ spores.ml⁻¹. For the analysis of OG production, samples were diluted to 1 mg.ml⁻¹ in 50 mM ammonium formate containing 0.1% formic acid. Chromatographic separation utilized the ACQUITY UPLC Protein BEH SEC Column (125 Å, 1.7 µm, 4.6 mm × 300 mm; Waters Corporation, Milford, MA, USA) with elution in 50 mM ammonium formate and 0.1% formic acid at a flow rate of 400 l.min⁻¹ and a column oven temperature of 40°C and was made using an HPLC system (UltiMate 3000 RS HPLC system, Thermo Scientific, Waltham, MA, USA). The system was coupled to an Impact II Ultra-High Resolution Qq-Time-Of-Flight (UHR-QqTOF) spectrometer (Bruker Daltonics, Bremen, Germany) equipped with an electrospray ionisation (ESI) source in
negative mode with the end plate offset set voltage to 500 V, capillary voltage to 4000 V, nebulizer to 40 psi, dry gas to 8 l/min and dry temperature of 180 °C. The Compass 1.8 software (Bruker Daltonics) was used to acquire the data. Peak annotation was performed based on accurate mass annotation, isotopic pattern, and MS/MS analysis¹⁷.

## Statistics and Reproducibility

All statistical analyses were done using the RStudio software (2024.09.1, Build 394). All bar graphs represent the mean value and the standard error as errors bars, points represent individual sample values. In the box plots median and interquartile range are indicated, with median shown as black line inside the boxes and interquartile range shown as the length of the boxes. For oligoprofiling data, an ANOVA followed by a Fisher's LSD post hoc test were conducted for multiple comparisons. For data presented with the aim of comparing low N and high N conditions or WT and mutants, a two-sample unpaired T-test was conducted, if both comparisons were done at the same time a two-way ANOVA (or a three-way ANOVA if the time factor was added) followed by a Fisher's LSD post hoc test was used. Differences were considered significant when the p value was lower than 0.05. All gene expression analyses were done with three technical replicates and one of the three independent experiment showing similar data was presented. For oligoprofiling data, three or four technical replicates and one of at least three independent experiment showing similar data was presented. For pathogenicity tests, at least 15 technical replicates each representing one inoculated detached leaf were done for each individual experiment repeated at least three times.

## Reporting summary

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

## Data availability

All data supporting the findings of this study are provided within the paper and its Supplementary Information. Source data for the charts and graphs presented in the main figures are available as Supplementary Data. The raw LC-MS datasets can be accessed at https://doi.org/10.6084/m9.fishshare.28290053.v1. Materials and all other data are available from the corresponding authors upon reasonable request.

Received: 19 March 2024; Accepted: 31 January 2025; Published online: 08 February 2025

## References

1. Fillinger, S. & Elad, Y. Botrytis—The Fungus, the Pathogen and its Management in Agricultural Systems 1–486. https://doi.org/10.1007/978-3-319-23371-0 (Springer, 2016).

2. Cheung, N., Tian, L., Liu, X. & Li, X. The destructive fungal pathogen Botrytis cinerea—insights from genes studied with mutant analysis. Pathogens 9, 1–46 (2020).

3. Amselem, J. et al. Genomic analysis of the necrotrophic fungal pathogens sclerotinia sclerotiorum and botrytis cinerea. PLoS Genet. 7, e10022230 (2011).

4. Lengyel, S. et al. Snf1 Kinase Differentially Regulates Botrytis cinerea Pathogenicity according to the Plant Host. Microorganisms 10, 1–20 (2022).

5. Rascle, C. et al. The pH regulator PacC: a host-dependent virulence factor in Botrytis cinerea. Environ. Microbiol. Rep. 10, 555–568 (2018).

6. Arbelet, D. et al. Disruption of the Bcchs3a chitin synthase gene in botrytis cinerea is responsible for altered adhesion and overstimulation of host plant immunity. Mol. Plant Microbe Interact. 23, 1324–1334 (2020).

7. Bi, K. et al. The Botrytis cinerea Crh1 transglycosylase is a cytoplasmic effector triggering plant cell death and defense response. Nat. Commun. 12, 2166 (2021).

8. Dalmais, B. et al. The Botrytis cinerea phytotoxin botcinic acid requires two polyketide synthases for production and has a redundant role in virulence with botrydial. Mol. Plant Pathol. 12, 564–579 (2011).

9. Zhu, W. et al. BcXYG1, a secreted xyloglucanase from Botrytis cinerea, triggers both cell death and plant immune responses. Plant Physiol. 175, 438–456 (2017).

10. Farjad, M. et al. Plant nitrate supply regulates Erwinia amylovora virulence gene expression in Arabidopsis. Mol. Plant Pathol. 22, 1332–1346 (2021).

11. Zarattini, M. et al. Every cloud has a silver lining: How abiotic stresses affect gene expression in plant-pathogen interactions. J. Exp. Bot. 72, 1020–1033 (2021).

12. Soulie, M. C. et al. Plant nitrogen supply affects the Botrytis cinerea infection process and modulates known and novel virulence factors. Mol. Plant Pathol. 21, 1436–1450 (2020).

13. Espino, J. J. et al. The Botrytis cinerea early secretome. Proteomics 10, 3020–3034 (2010).

14. Reynoud, N. et al. The complex architecture of plant cuticles and its relation to multiple biological functions. Front. Plant Sci. 12, 782773 (2021).

15. Leroch, M. et al. Transcriptome profiling of Botrytis cinerea conidial germination reveals upregulation of infection-related genes during the prepenetration stage. Eukaryot. Cell 12, 614–626 (2013).

16. An, H. J. et al. Determination of pathogen-related enzyme action by mass spectrometry analysis of pectin breakdown products of plant cell walls. Anal. Biochem 338, 71–82 (2005).

17. Voxeur, A. et al. Oligogalacturonide production upon Arabidopsis thaliana-Botrytis cinerea interaction. Proc. Natl. Acad. Sci. USA 116, 19743–19752 (2019).

18. Yadav, S., Yadav, P. K., Yadav, D. & Yadav, K. D. S. Pectin lyase: a review. Process Biochem. 44, 1–10 (2009).

19. Kong, G. et al. Pectin acylesterase PAE5 is associated with the virulence of plant pathogenic oomycete Peronophythora litchii. Physiol. Mol. Plant Pathol. 106, 16–22 (2019).

20. Pogorelko, G. et al. Arabidopsis and Brachypodium distachyon transgenic plants expressing Aspergillus nidulans acylesterases have decreased degree of polysaccharide acetylation and increased resistance to pathogens. Plant Physiol. 162, 9–23 (2013).

21. Blanco-Ulate, B. et al. Genome-wide transcriptional profiling of Botrytis cinerea genes targeting plant cell walls during infections of different hosts. Front. Plant Sci. 5, 1–16 (2014).

22. Kars, I. et al. Necrotizing activity of five Botrytis cinerea endopolygalacturonases produced in Pichia pastoris. Plant J. 43, 213–225 (2005).

23. Ten Have, A., Mulder, W., Visser, J. & Van Kan, J. A. L. The endopolygalacturonase gene Bcpg1 is required to full virulence of Botrytis cinerea. Mol. Plant Microbe Interact. 11, 1009–1016 (1998).

24. Pontiggia, D., Benedetti, M., Costantini, S., De Lorenzo, G. & Cervone, F. Dampening the DAMPs: how plants maintain the homeostasis of cell wall molecular patterns and avoid hyper-immunity. Front. Plant Sci. 11, 613259 (2020).

25. Davidsson, P. et al. Short oligogalacturonides induce pathogen resistance-associated gene expression in Arabidopsis thaliana. BMC Plant Biol. 17, 19 (2017).

26. Benedetti, M. et al. Four Arabidopsis berberine-bridge enzyme-like proteins are specific oxidases that inactivate the elicitor-active oligogalacturonides. Plant J. 94, 260–273 (2018).

27. Scortica, A. et al. Berberine bridge enzyme-like oligosaccharide oxidases act as enzymatic transducers between microbial glycoside hydrolases and plant peroxidases. Mol. Plant Microbe Interact. 35, 881–886 (2022).

28. Ferrari, S. et al. Resistance to Botrytis cinerea induced in arabidopsis by elicitors is independent of salicylic acid, ethylene, or jasmonate signaling but requires PHYTOALEXIN DEFICIENT3. Plant Physiol. 144, 367–379 (2007).
29. Hahn, M. G., Darvill, A. G. & Albersheim, P. Host-Pathogen Interactions XIX. THE ENDOGENOUS ELICITOR, A FRAGMENT OF A PLANT CELL WALL POLYSACCHARIDE THAT ELICITS PHYTOALEXIN ACCUMULATION IN SOYBEANS. *Plant Physiol.* **68**, 1161–1169 (1981).

30. Eizner, E. et al. Characterization of Botrytis–plant interactions using PathTrack®—an automated system for dynamic analysis of disease development. *Mol. Plant Pathol.* **18**, 503–512 (2017).

31. Caarls, L. et al. Arabidopsis JASMONATE-INDUCED OXYGENASES down-regulate plant immunity by hydroxylation and inactivation of the hormone jasmonic acid. *Proc. Natl. Acad. Sci. USA* **114**, 6388–6393 (2017).

32. Smirnova, E. et al. Jasmonic acid oxidase 2 hydroxylates jasmonic acid and represses basal defense and resistance responses against Botrytis cinerea Infection. *Mol. Plant* **10**, 1159–1173 (2017).

33. Zhang, F. et al. Structural basis of JAZ repression of MYC transcription factors in jasmonate signalling. *Nature* **525**, 269–273 (2015).

34. Windram, O. et al. Arabidopsis defense against Botrytis cinerea: chronology and regulation deciphered by high-resolution temporal transcriptomic analysis. *Plant Cell* **24**, 3530–3557 (2012).

35. Pieterse, C. M. J., Van der Does, D., Zamioudis, C., Leon-Reyes, A. & Van Wees, S. C. M. Hormonal modulation of plant immunity. *Annu Rev. Cell Dev. Biol.* **28**, 489–521 (2012).

36. Li, J., Brader, G. & Palva, E. T. The WRKY70 transcription factor: a node of convergence for jasmonate-mediated and salicylate-mediated signals in plant defense. *Plant Cell* **16**, 319–331 (2004).

37. Ren, C. M. et al. Transcription factor WRKY70 displays important but no indispensable roles in jasmonate and salicylic acid signaling. *J. Integr. Plant Biol.* **50**, 630–637 (2008).

38. Leisen, T. et al. CRISPR/Cas with ribonucleoprotein complexes and transiently selected telomere vectors allows highly efficient marker-free and multiple genome editing in Botrytis cinerea. *PLoS Pathog.* **16**, e1008326 (2020).

39. Jeblick, T. et al. Botrytis hypersensitive response inducing protein 1 triggers noncanonical PTI to induce plant cell death. *Plant Physiol.* **191**, 125–141 (2023).

40. Kars, I., McCalman, M., Wagemakers, L. & Van Kan, J. A. L. Functional analysis of Botrytis cinerea pectin methylesterase genes by PCR-based targeted mutagenesis: Bcpme1 and Bcpme2 are dispensable for virulence of strain B05.10. *Mol. Plant Pathol.* **6**, 641–652 (2005).

41. Ferrari, S., Plotnikova, J. M., De Lorenzo, G. & Ausubel, F. M. Arabidopsis local resistance to Botrytis cinerea involves salicylic acid and camalexin and requires EDS4 and PAD2, but not SID2, EDS5 or PAD4. *Plant J.* **35**, 193–205 (2003).

42. Han, Z., Xiong, D., Schneiter, R. & Tian, C. The function of plant PR1 and other members of the CAP protein superfamily in plant-pathogen interactions. *Mol. Plant Pathol.* **24**, 651–668 (2023).

43. Oide, S. et al. A novel role of PR2 in abscisic acid (ABA) mediated, pathogen-induced callose deposition in Arabidopsis thaliana. *N. Phytol.* **200**, 1187–1199 (2013).

44. Zhang, J. et al. Expression of pathogenesis related genes in response to salicylic acid, methyl jasmonate and 1-aminocyclopropane-1-carboxylic acid in Malus hupehensis (Pamp.) Rehd. BMC Res. Notes. https://www.biomedcentral.com/1756-0500/3/208 (2010).

45. López-Pérez, M., Ballester, A. R. & González-Candelas, L. Identification and functional analysis of Penicillium digitatum genes putatively involved in virulence towards citrus fruit. *Mol. Plant Pathol.* **16**, 262–275 (2015).

46. Chen, J. Y. et al. Characterization of the Verticillium dahliae exoproteome involves in pathogenicity from cotton-containing medium. *Front. Microbiol* **7**, 1–15 (2016).

47. Pirhonen, M., Saarilahti, H., Karlsson, M.-B. & Tapio Palvo, E. Identification of pathogenicity determinants of Erwinia carotovora subsp. carotovora by transposon mutagenesis. *Mol. Plant Microbes Interact.* **4**, 276–283 (1991).

48. Cotoras, M. & Silva, E. Differences in the initial events of infection of Botrytis cinerea strains isolated from tomato and grape. *Mycologia* **97**, 485–492 (2005).

49. Movahedi, S. & Heale, J. B. The roles of aspartic proteinase and endo-pectin lyase enzymes in the primary stages of infection and pathogenesis of various host tissues by different isolates of Botrytis cinerea Pers ex. *Pers. Physiol. Mol. Plant Pathol.* **36**, 303–324 (1990).

50. AbuQamar, S., Moustafa, K. & Tran, L. S. P. Mechanisms and strategies of plant defense against Botrytis cinerea. *Crit. Rev. Biotechnol.* **37**, 262–274 (2017).

51. Thomma, B. P. H. J. et al. Separate jasmonate-dependent and salicylate-dependent defense-response pathways in arabidopsis are essential for resistance to distinct microbial pathogens. *Proc. Natl. Acad. Sci. USA* **95**, 15107–15111 (1998).

52. Manabe, Y. et al. Loss-of-function mutation of REDUCED WALL ACETYLATION2 in Arabidopsis leads to reduced cell wall acetylation and increased resistance to Botrytis cinerea. *Plant Physiol.* **155**, 1068–1078 (2011).

53. Zhang, W. et al. Plant-necrotroph co-transcriptome networks illuminate a metabolic battlefield. *Elife* **8**, 1–32 (2019).

54. Herold, L. et al. Arabidopsis WALL-ASSOCIATED KINASES are not required for oligogalacturonide-induced signaling and immunity. *Plant Cell* **37**, koae317 (2024).

55. Quidde, T., Osbourn, A. E. & Tudzynski, P. Detoxification of α-tomatine by Botrytis cinerea. *Physiol. Mol. Plant Pathol.* **52**, 151–165 (1998).

56. Schumacher, J. Tools for Botrytis cinerea: new expression vectors make the gray mold fungus more accessible to cell biology approaches. *Fungal Genet. Biol.* **49**, 483–497 (2012).

57. Feng, J. et al. An inexpensive method for extraction of genomic DNA from fungal mycelia. *Can. J. Plant Pathol.* **32**, 396–401 (2010).

58. Concordet, J. P. & Haeussler, M. CRISPOR: Intuitive guide selection for CRISPR/Cas9 genome editing experiments and screens. *Nucleic Acids Res.* **46**, W242–W245 (2018).

59. Gachon, C. & Saindrenan, P. Real-time PCR monitoring of fungal development in Arabidopsis thaliana infected by Alternaria brassicicola and Botrytis cinerea. *Plant Physiol. Biochem.* **42**, 367–371 (2004).

60. Chassot, C., Buchala, A., Schoonbeek, H. J., Métraux, J. P. & Lamotte, O. Wounding of Arabidopsis leaves causes a powerful but transient protection against Botrytis infection. *Plant J.* **55**, 555–567 (2008).

61. Ren, H. et al. Selection of reliable reference genes for gene expression studies in Botrytis cinerea. *J. Microbiol Methods* **142**, 71–75 (2017).

## Acknowledgements

This work has benefited from the support of IJPB's Plant Observatory technological platforms and from the following funding French National Research Agency ANR-22-CE43-0013-WALLDERIVE (A.V., S.V., and L.L.). M.F., A.D., and M.C.S. were awarded the Eqolysin Grant (Award Number: BAP2022_35). A.D. was awarded a grant from the doctoral school SEVE of U. P. Saclay. IJPB received support from Saclay Plant Sciences-SPS (ANR-17-EUR-0007).

## Author contributions

Conceptualization: A.D., A.V., S.V., M.C.S., and M.F. Methodology: A.D., A.V., M.C.S., M.F., L.L., and S.J. Investigation: A.D. and A.V. Visualization: A.D. and A.V. Supervision: M.F., M.C.S., and A.V. Writing—original draft: A.D. and A.V. Writing—review & editing: A.D., A.V., S.V., M.C.S., and M.F.

## Competing interests

The authors declare no competing interests.
**Additional information**

**Supplementary information** The online version contains supplementary material available at
https://doi.org/10.1038/s42003-025-07642-7.

**Correspondence** and requests for materials should be addressed to
Aline Voxeur.

**Peer review information** *Communications Biology* thanks the anonymous reviewers for their contribution to the peer review of this work. Primary Handling Editors: David Favero and Manuel Breuer. A peer review file is available.

**Reprints and permissions information** is available at
http://www.nature.com/reprints

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

**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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence 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 licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2025